mRNA Interactions Research Articles

An ensemble of stacking classifiers for improved prediction of miRNA–mRNA interactions

MicroRNAs (miRNAs) are small non-coding RNA molecules that play a crucial role in regulating gene expression at the post-transcriptional level by binding to potential target sites of messenger RNAs (mRNAs), facilitated by the Argonaute family of proteins. Selecting the conservative candidate target sites (CTS) is a challenging step, considering that most of the existing computational algorithms primarily focus on canonical site types, which is a time-consuming and inefficient utilization of miRNA target site interactions. We developed a stacking classifier algorithm that addresses the CTS selection criteria using feature-encoding techniques that generates feature vectors, including k-mer nucleotide composition, dinucleotide composition, pseudo-nucleotide composition, and sequence order coupling. This innovative stacking classifier algorithm surpassed previous state-of-the-art algorithms in predicting functional miRNA targets. We evaluated the performance of the proposed model on 10 independent test datasets and obtained an average accuracy of 79.77%, which is a significant improvement of 7.26 % over previous models. This improvement shows that the proposed method has great potential for distinguishing highly functional miRNA targets and can serve as a valuable tool in biomedical and drug development research.

SARS-CoV-2 Nsp1 mediated mRNA degradation requires mRNA interaction with the ribosome

ABSTRACT Nsp1 is a SARS-CoV-2 host shutoff factor that both represses cellular translation and promotes host RNA decay. However, it is unclear how these two activities are connected and interact with normal translation processes. Here, we performed mutational analyses of Nsp1, and these revealed that both the N and C terminal domains of Nsp1 are important for translational repression. Furthermore, we demonstrate that specific residues in the N terminal domain are required for cellular RNA degradation but not bulk translation shutoff of host mRNAs, thereby separating RNA degradation from translation repression. We also present evidence that Nsp1 mediated RNA degradation requires engagement of the ribosome with mRNA. First, we observe that cytosolic lncRNAs, which are not translated, escape Nsp1 mediated degradation. Second, inhibition of translation elongation with emetine does not prevent Nsp1 mediated degradation, while blocking translation initiation before 48S ribosome loading reduces mRNA degradation. Taken together, we suggest that Nsp1 represses translation and promotes mRNA degradation only after ribosome engagement with the mRNA. This raises the possibility that Nsp1 may trigger RNA degradation through pathways that recognize stalled ribosomes.

An miRNA with a major impact on stress

miRNAs regulate mRNAs, including those important for synapse function in the brain. Mucha and colleagues recently identified a novel miRNA–mRNA interaction in the basolateral amygdala that acts as a homeostatic counter to stress-induced anxiety and synaptic changes, suggesting miRNAs as potential avenues for therapeutic intervention in anxiety disorders.

Circular RNA expression profiles and circRNA-miRNA-mRNA crosstalk in pre-eclamptic placenta

Objective: To identify the expression profile of circular RNA (circRNA) in placenta of pre-eclampsia (PE) pregnant women by high-throughput sequencing, and to construct the circRNA-microRNA (miRNA)-messenger RNA (mRNA) interaction network, so as to reveal the related pathways and regulatory mechanisms of PE. Methods: The clinical data and placentas of 42 women with PE (PE group) and 30 normal pregnant women (control group) who delivered in West China Second University Hospital from November 2019 to June 2021 were collected. (1) High-throughput sequencing was used to establish the differentially expressed circRNA profiles in placental tissues of 5 pairs of PE group and the control group. (2) Real-time quantitative PCR (qRT-PCR) was used to verify the expression levels of 6 differentially expressed circRNAs in placental tissues of PE group and control group. (3) Bioinformatics analysis was used to predict the target miRNA and analyze the co-expressed mRNA to construct a competitive endogenous RNA (ceRNA) network. The differentially expressed circRNAs were analyzed by Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathways. (4) Logistic regression analysis, Pearson correlation and Kendall's tau-b correlation analysis were used to test the correlation between the three differentially expressed circRNAs and the risk of PE and clinical characteristics. (5) circRNA_05393 was selected for subsequent functional study. Small interfering RNA (siRNA) and overexpression plasmid were used to knock down or increase the expression level of circRNA_05393 in trophoblast cell line HTR-8/SVneo cells, respectively. Transwell assay was used to detect the migration and invasion ability of the trophoblasts in vitro. Cell counting kit-8 assay was used to detect the proliferation ability of the trophoblasts. Results: (1) Seventy-two differentially expressed circRNAs were identified by high-throughput sequencing, of which 35 were up-regulated and 37 were down-regulated. (2) qRT-PCR showed that compared with the control group, circRNA_00673 (1.306±0.168 vs 2.059±0.242; t=2.356, P=0.021) and circRNA_07796 (1.275±0.232 vs 1.954±0.230; t=2.018, P=0.047) were significantly increased, while circRNA_05393 (1.846±0.377 vs 0.790±0.094; t=3.138, P=0.002) was significantly decreased. (3) The circRNA-miRNA-mRNA interaction network contained 3 circRNAs, 8 miRNAs and 53 mRNAs. GO functional annotation analysis showed that the biological process was mainly enriched in iron ion homeostasis, membrane depolarization during action potential and neuronal action potential. In terms of cellular components, they were mainly enriched in cytoskeleton and membrane components. In terms of molecular function, they were mainly enriched in the activity of voltage-gated sodium channel and basic amino acid transmembrane transporter. KEGG pathway enrichment analysis showed that mRNAs in the interaction network were mainly enriched in complement and coagulation cascade, glycine, serine and threonine metabolism, p53 signaling pathway and peroxisome proliferators-activated receptors (PPAR) signaling pathway. (4) Logistic regression analysis showed that down-regulation of circRNA_05393 expression was a risk factor for PE (OR=0.044, 95%CI: 0.003-0.596; P=0.019). Correlation analysis showed that circRNA_05393 was significantly correlated with systolic blood pressure and diastolic blood pressure in PE pregnant women (both P<0.05). (5) Knock down or overexpression of circRNA_05393 significantly reduced or increased the migration and invasion abilities of HTR-8/SVneo cells (all P<0.05), but had no significant effect on the ability of tube formation and proliferation (all P>0.05). Conclusions: The construction of circRNA expression profile in placenta and the exploration of circRNA-miRNA-mRNA interaction network provide the possibility to reveal the regulatory mechanism of specific circRNA involved in PE. Inhibition of circRNA_05393 may induce the progression of PE by reducing the migration and invasion of trophoblasts.

Integration analysis of miRNA–mRNA expression exploring their potential roles in intrahepatic cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is the second common primary hepatic malignancy tumor. In this study, an integrative analysis of differentially expressed genes (DEGs) and miRNAs from the ICC onset and adjacent normal tissues were performed to explore the regulatory roles of miRNA–mRNA interaction. A total of 1018 DEGs and 39 miRNAs were likely involved in ICC pathogenesis, suggesting the changes in cell metabolism in ICC development. The built network indicated that 30 DEGs were regulated by 16 differentially expressed miRNA. The screened DEGs and miRNA together were probably considered the biomarkers of ICC, and their important roles in ICC pathogenesis remain to be elucidated. This study could provide a good basis to uncover the regulatory mechanism of miRNA and mRNAs in ICC pathogenesis.

Uterus proliferative period ceRNA network of Yunshang black goat reveals candidate genes on different kidding number trait.

Pregnancy loss that occurs in the uterus is an important and widespread problem in humans and farm animals and is also a key factor affecting the fecundity of livestock. Understanding the differences in the fecundity of goats may be helpful in guiding the breeding of goats with high fecundity. In this study, we performed RNA sequencing and bioinformatics analysis to study the uterus of Yunshang black goats with high and low fecundity in the proliferative period. We identified mRNAs, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) by analyzing the uterine transcriptomes. The target genes of the identified miRNAs and lncRNAs were predicted, and miRNA-mRNA interaction and competitive endogenous RNA (ceRNA) networks were constructed. By comparisons between low- and high-fecundity groups, we identified 1,674 differentially expressed mRNAs (914 were upregulated, and 760 were downregulated), 288 differentially expressed lncRNAs (149 were upregulated, and 139 were downregulated), and 17 differentially expressed miRNAs (4 were upregulated, and 13 were downregulated). In addition, 49 miRNA-mRNA pairs and 45 miRNA-lncRNA pairs were predicted in the interaction networks. We successfully constructed a ceRNA interaction network with 108 edges that contained 19 miRNAs, 11 mRNAs, and 73 lncRNAs. Five candidate genes (PLEKHA7, FAT2, FN1, SYK, and ITPR2) that were annotated as cell adhesion or calcium membrane channel protein were identified. Our results provide the overall expression profiles of mRNAs, lncRNAs, and miRNAs in the goat uterus during the proliferative period and are a valuable reference for studies into the mechanisms associated with the high fecundity, which may be helpful to guide goat to reduce pregnancy loss.

Value of negatively correlated miR-205-5p/HMGB3 and miR-96-5p/FOXO1 on the diagnosis of breast cancer and benign breast diseases

BackgroundMicroRNA (miRNA) and mRNA levels in matching specimens were used to identify miRNA–mRNA interactions. We aimed to integrate transcriptome, immunophenotype, methylation, mutation, and survival data analyses to examine the profiles of miRNAs and target mRNAs and their associations with breast cancer (BC) diagnosis. MethodsBased on the Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA), differentially expressed miRNAs and targeted mRNAs were screened from experimentally verified miRNA-target interaction databases using Pearson's correlation analysis. We used real-time quantitative reverse transcription polymerase chain reaction to verify BC and benign disease samples, and logistic regression analysis was used to establish a diagnostic model based on miRNAs and target mRNAs. Receiver operating characteristic curve analysis was performed to test the ability to recognize the miRNA-mRNA pairs. Next, we investigated the complex interactions between miRNA-mRNA regulatory pairs and phenotypic hallmarks. ResultsWe identified 27 and 359 dysregulated miRNAs and mRNAs, respectively, based on the GEO and TCGA databases. Using Pearson's correlation analysis, 10 negative miRNA-mRNA regulatory pairs were identified after screening both databases, and the related miRNA and target mRNA levels were assessed in 40 BC tissues and 40 benign breast disease tissues. Two key regulatory pairs (miR-205-5p/High mobility group box 3 (HMGB3) and miR-96-5p/Forkhead Box O1 (FOXO1)) were selected to establish the diagnostic model. They also had utility in survival and clinical analyses. ConclusionsA diagnostic model including two miRNAs and their respective target mRNAs was established to distinguish between BC and benign breast diseases. These markers play essential roles in BC pathogenesis.

Synchronization of fractional-order repressilatory genetic oscillators with time delay

Genetic oscillators have been widely used in modeling key processes of biological systems, especially cell cycles and circadian rhythms. In particular, repressilatory genetic oscillators have been employed in modeling the dynamics of mRNA and protein interactions with transcriptional and translational feedback loops at the molecular level. In addition, synchronization of these oscillators is crucial for understanding the underlying mechanisms of the associated biological processes. In this paper, models of fractional-order genetic oscillators and their coupling are established, where the aspects of time delay, coupling strength, noise, and stability are all taken into consideration. Communication in the proposed coupling model is based on quorum sensing. The synchronization of the fractional-order repressilator model has been examined through simulations which show three main findings. Firstly, the synchronization of the fractional-order repressilator model can be optimized through coupling weight selection. Secondly, the synchronization can be enhanced by increasing the fractional order and decreasing the time delay and the noise intensity. Finally, transitions between the states of the fractional-order repressilatory oscillator can be achieved through varying the fractional order. The simulation results verify the biological relevance of the genetic oscillator models, and their potential for explaining the underlying mechanisms of the associated biological processes.

Screening prognostic markers for hepatocellular carcinoma based on pyroptosis-related lncRNA pairs

BackgroundPyroptosis is closely related to cancer prognosis. In this study, we tried to construct an individualized prognostic risk model for hepatocellular carcinoma (HCC) based on within-sample relative expression orderings (REOs) of pyroptosis-related lncRNAs (PRlncRNAs).MethodsRNA-seq data of 343 HCC samples derived from The Cancer Genome Atlas (TCGA) database were analyzed. PRlncRNAs were detected based on differentially expressed lncRNAs between sample groups clustered by 40 reported pyroptosis-related genes (PRGs). Univariate Cox regression was used to screen out prognosis-related PRlncRNA pairs. Then, based on REOs of prognosis-related PRlncRNA pairs, a risk model for HCC was constructed by combining LASSO and stepwise multivariate Cox regression analysis. Finally, a prognosis-related competing endogenous RNA (ceRNA) network was built based on information about lncRNA–miRNA–mRNA interactions derived from the miRNet and TargetScan databases.ResultsHierarchical clustering of HCC patients according to the 40 PRGs identified two groups with a significant survival difference (Kaplan–Meier log-rank, p = 0.026). Between the two groups, 104 differentially expressed lncRNAs were identified (|log2(FC)|> 1 and FDR < 5%). Among them, 83 PRlncRNA pairs showed significant associations between their REOs within HCC samples and overall survival (Univariate Cox regression, p < 0.005). An optimal 11-PRlncRNA-pair prognostic risk model was constructed for HCC. The areas under the curves (AUCs) of time-dependent receiver operating characteristic (ROC) curves of the risk model for 1-, 3-, and 5-year survival were 0.737, 0.705, and 0.797 in the validation set, respectively. Gene Set Enrichment Analysis showed that inflammation-related interleukin signaling pathways were upregulated in the predicted high-risk group (p < 0.05). Tumor immune infiltration analysis revealed a higher abundance of regulatory T cells (Tregs) and M2 macrophages and a lower abundance of CD8 + T cells in the high-risk group, indicating that excessive pyroptosis might occur in high-risk patients. Finally, eleven lncRNA–miRNA–mRNA regulatory axes associated with pyroptosis were established.ConclusionOur risk model allowed us to determine the robustness of the REO-based PRlncRNA prognostic biomarkers in the stratification of HCC patients at high and low risk. The model is also helpful for understanding the molecular mechanisms between pyroptosis and HCC prognosis. High-risk patients may have excessive pyroptosis and thus be less sensitive to immune therapy.

Symbiosis preservation: Putative regulation of fatty acyl-CoA reductase by miR-31a within the symbiont harboring bacteriome through tsetse evolution.

Tsetse flies are the sole vectors of African trypanosomes. In addition to trypanosomes, tsetse harbor obligate Wigglesworthia glossinidia bacteria that are essential to tsetse biology. The absence of Wigglesworthia results in fly sterility, thus offering promise for population control strategies. Here, microRNA (miRNAs) and mRNA expression are characterized and compared between the exclusive Wigglesworthia-containing bacteriome and adjacent aposymbiotic tissue in females of two evolutionarily distant tsetse species (Glossina brevipalpis and G. morsitans). A total of 193 miRNAs were expressed in either species, with 188 of these expressed in both species, 166 of these were novel to Glossinidae, and 41 miRNAs exhibited comparable expression levels between species. Within bacteriomes, 83 homologous mRNAs demonstrated differential expression between G. morsitans aposymbiotic and bacteriome tissues, with 21 of these having conserved interspecific expression. A large proportion of these differentially expressed genes are involved in amino acid metabolism and transport, symbolizing the essential nutritional role of the symbiosis. Further bioinformatic analyses identified a sole conserved miRNA::mRNA interaction (miR-31a::fatty acyl-CoA reductase) within bacteriomes likely catalyzing the reduction of fatty acids to alcohols which comprise components of esters and lipids involved in structural maintenance. The Glossina fatty acyl-CoA reductase gene family is characterized here through phylogenetic analyses to further understand its evolutionary diversification and the functional roles of members. Further research to characterize the nature of the miR-31a::fatty acyl-CoA reductase interaction may find novel contributions to the symbiosis to be exploited for vector control.

Abstract 1450: Global changes in RBPs and mRNA interactions during prostate cancer progression

Abstract Prostate cancer (PC) is a leading cause of cancer in men in the United States and worldwide, disproportionately affecting men of African ancestry. The etiology and molecular mechanism contributing to lethal human PC remain poorly understood. RNA binding proteins (RBPs) controlled all aspects of RNA biogenesis, and dysregulated RBPs are linked to human diseases, including cancer. However, the global changes in RBP-mRNA interactions during PC progression are unknown. Here, we have investigated the patterns of RBPs and polyA-mRNA interaction in the castration-sensitive LNCaP and castration-resistant C4-2 PC cell lines. Our enhanced-RNA interactome capture assay combined with mass spectroscopy (MS) revealed that LNCaP and C4-2 cells showed distinct RBP-mRNA interaction patterns in steady-state and androgen-treated conditions. Also, we have studied the interaction between Yes-associated protein 1 (YAP1) and nucleophosmin (NPM1) RBP. Androgen increased NPM1 and mRNA interaction, and NPM1 is part of the YAP1 multiprotein complex, as revealed by MS-based proteomics approaches. Our proximity ligation assay showed that androgen hormone signaling regulated the interaction between YAP1 and NPM1. Besides, our GST-pulldown demonstrated that NPM1 binds the proline-rich region domain of YAP1, whereas the WW/SH3 domain of YAP1 mediates the interaction in C4-2 cells. Genetic silencing of NPM1 reduced PC cell growth in culture. Furthermore, the upregulation of NPM1 correlates with prostate adenocarcinoma and metastasis. These observations suggest that changes in RBP and RNA interaction patterns may play a critical role in PC progression and recurrence. Citation Format: Abdulrahman M. Dwead, Marwah M. Al-Mathkour, Bekir Cinar. Global changes in RBPs and mRNA interactions during prostate cancer progression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1450.

SInterBase: a comprehensive database of E.coli sRNA-mRNA interactions.

sInterBase is a comprehensive and easy-to-operate web-based platform for mining experimentally identified sRNA-mRNA interactions in E.coli. Interactions in the database are annotated with an interaction duplex and a set of descriptive features. sInterBase provides advanced functionality, such as flexible search based on various criteria, statistical analysis via charts, browsing, and downloading interactions for further use. sInterBase is available at https://sinterbase.cs.bgu.ac.il/. Supplementary data are available at Bioinformatics online.

Integrative pan-cancer analysis of cuproplasia-associated genes for the genomic and clinical characterization of 33 tumors.

The molecular mechanisms driving tumorigenesis have continually been the focus of researchers. Cuproplasia is defined as copper-dependent cell growth and proliferation, including its primary and secondary roles in tumor formation and proliferation through signaling pathways. In this study, we analyzed the differences in the expression of cuproplasia-associated genes (CAGs) in pan-cancerous tissues and investigated their role in immune-regulation and tumor prognostication. Raw data from 11,057 cancer samples were acquired from multiple databases. Pan-cancer analysis was conducted to analyze the CAG expression, single-nucleotide variants, copy number variants, methylation signatures, and genomic signatures of micro RNA (miRNA)-messenger RNA (mRNA) interactions. The Genomics of Drug Sensitivity in Cancer and the Cancer Therapeutics Response Portal databases were used to evaluate drug sensitivity and resistance against CAGs. Using single-sample Gene Set Enrichment Analysis (ssGSEA) and Immune Cell Abundance Identifier database, immune cell infiltration was analyzed with the ssGSEA score as the standard. Aberrantly expressed CAGs were found in multiple cancers. The frequency of single-nucleotide variations in CAGs ranged from 1% to 54% among different cancers. Furthermore, the correlation between CAG expression in the tumor microenvironment and immune cell infiltration varied among different cancers. ATP7A and ATP7B were negatively correlated with macrophages in 16 tumors including breast invasive carcinoma and esophageal carcinoma, while the converse was true for MT1A and MT2A. In addition, we established cuproplasia scores and demonstrated their strong correlation with patient prognosis, immunotherapy responsiveness, and disease progression (P < 0.05). Finally, we identified potential candidate drugs by matching gene targets with existing drugs. This study reports the genomic characterization and clinical features of CAGs in pan-cancers. It helps clarify the relationship between CAGs and tumorigenesis, and may be helpful in the development of biomarkers and new therapeutic agents.

CircPLXNA2 Affects the Proliferation and Apoptosis of Myoblast through circPLXNA2/gga-miR-12207-5P/MDM4 Axis

CircRNAs are newly identified special endogenous RNA molecules that covalently close a loop by back-splicing with pre-mRNA. In the cytoplasm, circRNAs would act as molecular sponges to bind with specific miRNA to promote the expression of target genes. However, knowledge of circRNA functional alternation in skeletal myogenesis is still in its infancy. In this study, we identified a circRNA-miRNA-mRNA interaction network in which the axis may be implicated in the progression of chicken primary myoblasts' (CPMs) myogenesis by multi-omics (i.e., circRNA-seq and ribo-seq). In total, 314 circRNA-miRNA-mRNA regulatory axes containing 66 circRNAs, 70 miRNAs, and 24 mRNAs that may be relevant to myogenesis were collected. With these, the circPLXNA2-gga-miR-12207-5P-MDM4 axis aroused our research interest. The circPLXNA2 is highly differentially expressed during differentiation versus proliferation. It was demonstrated that circPLXNA2 inhibited the process of apoptosis while at the same time stimulating cell proliferation. Furthermore, we demonstrated that circPLXNA2 could inhibit the repression of gga-miR-12207-5p to MDM4 by directing binding to gga-miR-12207-5p, thereby restoring MDM4 expression. In conclusion, circPLXNA2 could function as a competing endogenous RNA (ceRNA) to recover the function of MDM4 by directing binding to gga-miR-12207-5p, thereby regulating the myogenesis.

MicroRNA Target Identification: Revisiting Accessibility and Seed Anchoring.

By pairing to messenger RNAs (mRNAs for short), microRNAs (miRNAs) regulate gene expression in animals and plants. Accurately identifying which mRNAs interact with a given miRNA and the precise location of the interaction sites is crucial to reaching a more complete view of the regulatory network of an organism. Only a few experimental approaches, however, allow the identification of both within a single experiment. Computational predictions of miRNA-mRNA interactions thus remain generally the first step used, despite their drawback of a high rate of false-positive predictions. The major computational approaches available rely on a diversity of features, among which anchoring the miRNA seed and measuring mRNA accessibility are the key ones, with the first being universally used, while the use of the second remains controversial. Revisiting the importance of each is the aim of this paper, which uses Cross-Linking, Ligation, And Sequencing of Hybrids (CLASH) datasets to achieve this goal. Contrary to what might be expected, the results are more ambiguous regarding the use of the seed match as a feature, while accessibility appears to be a feature worth considering, indicating that, at least under some conditions, it may favour anchoring by miRNAs.

Diversity and Differential Expression of MicroRNAs in the Human Skeletal Muscle with Distinct Fiber Type Composition.

The ratio of fast- and slow-twitch fibers in human skeletal muscle is variable and largely determined by genetic factors. In this study, we investigated the contribution of microRNA (miRNA) in skeletal muscle fiber type composition. The study involved biopsy samples of the vastus lateralis muscle from 24 male participants with distinct fiber type ratios. The miRNA study included samples from five endurance athletes and five power athletes with the predominance of slow-twitch (61.6-72.8%) and fast-twitch (69.3-80.7%) fibers, respectively. Total and small RNA were extracted from tissue samples. Total RNA sequencing (N = 24) revealed 352 differentially expressed genes between the groups with the predominance of fast- and slow-twitch muscle fibers. Small RNA sequencing showed upregulation of miR-206, miR-501-3p and miR-185-5p, and downregulation of miR-499a-5p and miR-208-5p in the group of power athletes with fast-twitch fiber predominance. Two miRtronic miRNAs, miR-208b-3p and miR-499a-5p, had strong correlations in expression with their host genes (MYH7 and MYH7B, respectively). Correlations between the expression of miRNAs and their experimentally validated messenger RNA (mRNA) targets were calculated, and 11 miRNA-mRNA interactions with strong negative correlations were identified. Two of them belonged to miR-208b-3p and miR-499a-5p, indicating their regulatory links with the expression of CDKN1A and FOXO4, respectively.

MiRNA Expression Profiles of Mouse Round Spermatids in GRTH/DDX25-Mediated Spermiogenesis: mRNA-miRNA Network Analysis.

GRTH/DDX25 is a testis-specific DEAD-box family of RNA helicase, which plays an essential role in spermatogenesis and male fertility. There are two forms of GRTH, a 56 kDa non-phosphorylated form and a 61 kDa phosphorylated form (pGRTH). GRTH-KO and GRTH Knock-In (KI) mice with R242H mutation (lack pGRTH) are sterile with a spermatogenic arrest at step 8 of spermiogenesis due to failure of round spermatids (RS) to elongate. We performed mRNA-seq and miRNA-seq analysis on RS of WT, KI, and KO to identify crucial microRNAs (miRNAs) and mRNAs during RS development by establishing a miRNA-mRNA network. We identified increased levels of miRNAs such as miR146, miR122a, miR26a, miR27a, miR150, miR196a, and miR328 that are relevant to spermatogenesis. mRNA-miRNA target analysis on these DE-miRNAs and DE-mRNAs revealed miRNA target genes involved in ubiquitination process (Ube2k, Rnf138, Spata3), RS differentiation, and chromatin remodeling/compaction (Tnp1/2, Prm1/2/3, Tssk3/6), reversible protein phosphorylation (Pim1, Hipk1, Csnk1g2, Prkcq, Ppp2r5a), and acrosome stability (Pdzd8). Post-transcriptional and translational regulation of some of these germ-cell-specific mRNAs by miRNA-regulated translation arrest and/or decay may lead to spermatogenic arrest in KO and KI mice. Our studies demonstrate the importance of pGRTH in the chromatin compaction and remodeling process, which mediates the differentiation of RS into elongated spermatids through miRNA-mRNA interactions.

Plant non-coding RNAs function in pollen development and male sterility.

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Targeting HuR-Vav3 mRNA interaction prevents Pseudomonas aeruginosa adhesion to the cystic fibrosis airway epithelium

Cystic fibrosis (CF) is characterized by chronic bacterial infections leading to progressive bronchiectasis and respiratory failure. Pseudomonas aeruginosa (Pa) is the predominant opportunistic pathogen infecting the CF airways. The guanine nucleotide exchange factor Vav3 plays a critical role in Pa adhesion to the CF airways by inducing luminal fibronectin deposition that favors bacteria trapping. Here we report that Vav3 overexpression in CF is caused by upregulation of the mRNA-stabilizing protein HuR. We found that HuR accumulates in the cytoplasm of CF airway epithelial cells and that it binds to and stabilizes Vav3 mRNA. Interestingly, disruption of the HuR-Vav3 mRNA interaction improved the CF epithelial integrity, inhibited the formation of the fibronectin-made bacterial docking platforms, and prevented Pa adhesion to the CF airway epithelium. These findings indicate that targeting HuR represents a promising antiadhesive approach in CF that can prevent initial stages of Pa infection in a context of emergence of multidrug-resistant pathogens.

MiR-197-3p Promotes Osteosarcoma Stemness and Chemoresistance by Inhibiting SPOPL.

First-line treatment for osteosarcoma includes chemotherapy and surgery. However, the five-year survival rate of refractory osteosarcoma remains unsatisfactory. Osteosarcoma cancer stem cells, possessing stemness and chemoresistance, are one of the critical causes of poor response to chemotherapy. Elucidating regulatory signaling pathways of osteosarcoma cancer stem cells may provide a rationale for improving regimens against chemoresistant osteosarcoma. Methotrexate (MTX)-resistant osteosarcoma cells were established. microRNA expression profiles were used for detecting differentially expressed microRNA in resistant clones and the parental cells. microRNA target databases were employed to predict potential microRNA and mRNA interactions. Flow cytometry was performed to measure stem cell marker Prominin-1 (CD133)-positive cells. Immunofluorescence staining was applied to detect CD133 expression. miR-197-3p mimic or anti-miR-197-3p stably transfected cells were used to generate xenograft models. In the study, we found that miR-197-3p was increased in MTX-resistant cell lines. Overexpression of miR-197-3p enhanced the expression of cancer stem cell markers CD133, Octamer-binding protein 4 (OCT4), Transcription factor SOX-2 (SOX2), and Homeobox protein NANOG (NANOG), as well as chemoresistance-associated genes ATP-dependent translocase ABCB1 (ABCB1) and Broad substrate specificity ATP-binding cassette transporter ABCG2 (ABCG2), whereas miR-197-3p knockdown inhibited stemness and recovered sensitivity to MTX. We also classified the tumor suppressor Speckle-type POZ protein-like (SPOPL) as a target of miR-197-3p. The miR-197-3p mutation that could not combine SPOPL promoter regions was unable to sustain stemness or chemoresistance. Collectively, we discovered miR-197-3p conferred osteosarcoma stemness and chemotherapy resistance by targeting SPOPL, prompting promising therapeutic candidates for refractory osteosarcoma treatment.