Human MicroRNA targets

A preliminary study of calcium channel-associated mRNA and miRNA networks in post-traumatic epileptic rats

The calcium channels are the main pathogenesis and therapeutic target for post-traumatic epilepsy (PTE). However, differentially expressed miRNAs (DEMs) and mRNAs associated with calcium channels in PTE and their interactions are poorly understood. We produced a PTE model in rats and conducted RNA-seq in PTE rats. Gene annotation was used to verify differentially expressed mRNAs related to calcium channels. RNAhybrid, PITA, and Miranda prediction were used to build the miRNA-mRNA pairs. Furthermore, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used for the functional enrichment analysis of DEMs. The quantification changes of mRNA and miRNA were verified by RT-qPCR. There were 431 identified differentially expressed genes (DEGs) in PTE rats compared with the sham group, of which five mRNAs and 7 miRNAs were related to calcium channels. The miRNA-mRNA network suggested a negative correlation between 11 pairs of miRNA-mRNA involved in the p53 signaling pathway, HIF-1 signaling pathway. RT-qPCR verified three upregulated mRNAs in PTE rats, associated with 7 DEMs negatively related to them, respectively. This study has revealed the changes in miRNA-mRNA pairs associated with calcium channels in PTE, which might contribute to the further interpretation of potential underlying molecular mechanisms of PTE and the discovery of promising diagnostics.

Construction of chronic glomerulonephritis‑related lncRNA‑mRNA regulatory network and lncRNA‑-miRNA‑mRNA ceRNA network by bioinformatics analysis

Long non-coding RNAs (lncRNAs) are ncRNA transcripts >200 nucleotides that are important genetic regulators. LncRNAs can directly regulate mRNA through a lncRNA-mRNA regulatory mode and can also regulate mRNA through competitive binding to micro (mi)RNA, which is generally known as the competitive endogenous RNA (ceRNA) network. The present study evaluated the functional roles and regulatory networks of lncRNAs in chronic glomerulonephritis (CGN). The proliferative ability of mouse glomerular mesangial cells (GMCs) induced by different concentrations of lipopolysaccharide (LPS) was assessed using the Cell Counting Kit-8 assay, and RNA sequencing (RNA-seq) was performed to identify differentially expressed lncRNAs in LPS-induced GMCs. Based on the sequencing results, six lncRNAs were selected for validation using reverse transcription-quantitative PCR (RT-qPCR). Furthermore, the lncRNA-mRNA regulatory network and the lncRNA-miRNA-mRNA ceRNA network were constructed to assess the role and mechanism of CGN-related lncRNAs. To elucidate the biological functions of lncRNAs, Gene Ontology (GO) biological process term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on all mRNAs involved in the lncRNA-mRNA regulatory network and in the ceRNA network. A total of 1,532 differentially expressed lncRNAs, including 594 upregulated lncRNAs and 938 downregulated lncRNAs, were identified using RNA-seq. The results of RT-qPCR validation were consistent with RNA-seq results. An lncRNA-mRNA regulatory network, including 236 lncRNAs and 556 mRNAs, and a ceRNA network, including 6 lncRNAs, 18 miRNAs and 419 mRNAs, were successfully constructed. The GO biological process term enrichment and KEGG pathway enrichment analyses demonstrated that those lncRNAs were often related to inflammatory response and substance metabolism. The present study identified key CGN-related lncRNAs in LPS-induced GMCs, and further demonstrated a global view of the lncRNA-mRNA regulatory network and ceRNA network involved in CGN. These results offered novel insights into the roles of lncRNAs in the pathogenesis of CGN and identified potential diagnostic biomarkers.

Identification of markers for predicting prognosis and endocrine metabolism in nasopharyngeal carcinoma by miRNA-mRNA network mining and machine learning

Background

Nasopharyngeal cancer (NPC) has a high incidence in Southern China and Asia, and its survival is extremely poor in advanced patients. MiRNAs play critical roles in regulating gene expression and serve as therapeutic targets in cancer. This study sought to disclose key miRNAs and target genes responsible for NPC prognosis and endocrine metabolism.

Materials and methods

Three datasets (GSE32960, GSE70970, and GSE102349) of NPC samples came from Gene Expression Omnibus (GEO). Limma and WGCNA were applied to identify key prognostic miRNAs. There were 12 types of miRNA tools implemented to study potential target genes (mRNAs) of miRNAs. Univariate Cox regression and stepAIC were introduced to construct risk models. Pearson analysis was conducted to analyze the correlation between endocrine metabolism and RiskScore. Single-sample gene set enrichment analysis (ssGSEA), MCP-counter, and ESTIMATE were performed for immune analysis. The response to immunotherapy was predicted by TIDE and SubMap analyses.

Results

Two key miRNAs (miR-142-3p and miR-93) were closely involved in NPC prognosis. The expression of the two miRNAs was dysregulated in NPC cell lines. A total of 125 potential target genes of the key miRNAs were screened, and they were enriched in autophagy and mitophagy pathways. Five target genes (E2F1, KCNJ8, SUCO, HECTD1, and KIF23) were identified to construct a prognostic model, which was used to divide patients into high group and low group. RiskScore was negatively correlated with most endocrine-related genes and pathways. The low-risk group manifested higher immune infiltration, anticancer response, more activated immune-related pathways, and higher response to immunotherapy than the high-risk group.

Conclusions

This study revealed two key miRNAs that were highly contributable to NPC prognosis. We delineated the specific links between key miRNAs and prognostic mRNAs with miRNA-mRNA networks. The effectiveness of the five-gene model in predicting NPC prognosis as well as endocrine metabolism provided a guidance for personalized immunotherapy in NPC patients.

miRNA-559 and MTDH as possible diagnostic markers of psoriasis: Role of PTEN/AKT/FOXO pathway in disease pathogenesis

Psoriasis is a persistent, inflammatory, autoimmune skin disorder which can be elicited by genetic and environmental factors. Several microRNAs (miRNAs) that are abnormally expressed in psoriasis have emerged as an interesting candidate in psoriasis pathogenesis. However, the expression profile and function of miRNA-559, and its direct target metadherin (MTDH), in psoriasis need to be further illuminated. This study intended to assess miRNA-559 and MTDH levels in skin and sera of psoriatic patients and to investigate their clinical significance in an attempt for developing novel distinct tools for early diagnosis of psoriasis. Moreover, this study aimed at exploring participation of miRNA-559 in regulating MTDH/PTEN/AKT pathway in psoriasis. Expression levels of miRNA-559, AKT, FOXO1 and PTEN were measured by real-time qRT-PCR, whereas MTDH and p27 levels were assessed by ELISA in lesional, non-lesional tissues and serum of 20 psoriatic patients and 20 matching controls. Correlation study was conducted between different parameters. The diagnostic performance of miRNA-559 and MTDH in psoriasis was estimated by receiver operating characteristic (ROC) curve analysis. Expression of miRNA-559 in psoriatic patients was significantly downregulated in both lesional tissues and serum as compared to controls. Conversely, MTDH protein level showed significant increase in both tissues and serum of psoriatic patients and was inversely correlated with miRNA-559 level. Meanwhile, levels of PTEN, AKT and FOXO1 were dramatically changed in psoriatic patients compared to controls. Furthermore, serum miRNA-559 and MTDH displayed comparable diagnostic accuracy in discriminating psoriatic patients from controls. Yet, miRNA-559 demonstrated superior diagnostic performance than MTDH in psoriasis diagnosis. Together, the current findings provide the first suggestion of a new mechanism by which downregulation of miRNA-559 might induce proliferation in psoriasis through modulating PTEN/AKT/FOXO1 pathway by positive regulation of MTDH. Thus, miRNA-559 and MTDH might be proposed as promising diagnostic biomarkers of psoriasis.

Characterizations of microRNAs involved in the molecular mechanisms underlying the therapeutic effects of noni (Morinda citrifolia L.) fruit juice on hyperuricemia in mice

Background

Hyperuricemia is generally defined as the high level of serum uric acid and is well known as an important risk factor for the development of various medical disorders. However, the medicinal treatment of hyperuricemia is frequently associated with multiple side-effects.

Methods

The therapeutic effect of noni (Morinda citrifolia L.) fruit juice on hyperuricemia and the underlying molecular mechanisms were investigated in mouse model of hyperuricemia induced by potassium oxonate using biochemical and high-throughput RNA sequencing analyses.

Results

The levels of serum uric acid (UA) and xanthine oxidase (XOD) in mice treated with noni fruit juice were significantly decreased, suggesting that the noni fruit juice could alleviate hyperuricemia by inhibiting the XOD activity and reducing the level of serum UA. The contents of both serum creatinine and blood urine nitrogen of the noni fruit juice group were significantly lower than those of the model group, suggesting that noni fruit juice promoted the excretion of UA without causing deleterious effect on the renal functions in mice. The differentially expressed microRNAs involved in the pathogenesis of hyperuricemia in mice were identified by RNA sequencing with their target genes further annotated based on both Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases to explore the metabolic pathways and molecular mechanisms underlying the therapeutic effect on hyperuricemia by noni fruit juice.

Conclusion

Our study provided strong experimental evidence to support the further investigations of the potential application of noni fruit juice in the treatment of hyperuricemia.

Molecular and biochemical investigations of the anti-fatigue effects of tea polyphenols and fruit extracts of Lycium ruthenicum Murr. on mice with exercise-induced fatigue

Background: The molecular mechanisms regulating the therapeutic effects of plant-based ingredients on the exercise-induced fatigue (EIF) remain unclear. The therapeutic effects of both tea polyphenols (TP) and fruit extracts of Lycium ruthenicum (LR) on mouse model of EIF were investigated. Methods: The variations in the fatigue-related biochemical factors, i.e., lactate dehydrogenase (LDH), superoxide dismutase (SOD), tumor necrosis factor-α (TNF-α), interleukin-1β (IL-1β), interleukin-2 (IL-2), and interleukin-6 (IL-6), in mouse models of EIF treated with TP and LR were determined. The microRNAs involved in the therapeutic effects of TP and LR on the treatment of mice with EIF were identified using the next-generation sequencing technology. Results: Our results revealed that both TP and LR showed evident anti-inflammatory effect and reduced oxidative stress. In comparison with the control groups, the contents of LDH, TNF-α, IL-6, IL-1β, and IL-2 were significantly decreased and the contents of SOD were significantly increased in the experimental groups treated with either TP or LR. A total of 23 microRNAs (21 upregulated and 2 downregulated) identified for the first time by the high-throughput RNA sequencing were involved in the molecular response to EIF in mice treated with TP and LR. The regulatory functions of these microRNAs in the pathogenesis of EIF in mice were further explored based on Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses with a total of over 20,000-30,000 target genes annotated and 44 metabolic pathways enriched in the experimental groups based on GO and KEGG databases, respectively. Conclusion: Our study revealed the therapeutic effects of TP and LR and identified the microRNAs involved in the molecular mechanisms regulating the EIF in mice, providing strong experimental evidence to support further agricultural development of LR as well as the investigations and applications of TP and LR in the treatment of EIF in humans, including the professional athletes.

The PIWI-specific insertion module helps load longer piRNAs for translational activation essential for male fertility

PIWI-clade proteins harness piRNAs of 24-33 nt in length. Of great puzzles are how PIWI-clade proteins incorporate piRNAs of different sizes and whether the size matters to PIWI/piRNA function. Here we report that a PIWI-Ins module unique in PIWI-clade proteins helps define the length of piRNAs. Deletion of PIWI-Ins in Miwi shifts MIWI to load with shorter piRNAs and causes spermiogenic failure in mice, demonstrating the functional importance of this regulatory module. Mechanistically, we show that longer piRNAs provide additional complementarity to target mRNAs, thereby enhancing the assembly of the MIWI/eIF3f/HuR super-complex for translational activation. Importantly, we identify a c.1108C>T (p.R370W) mutation of HIWI (human PIWIL1) in infertile men and demonstrate in Miwi knock-in mice that this genetic mutation impairs male fertility by altering the property of PIWI-Ins in selecting longer piRNAs. These findings reveal a critical role of PIWI-Ins-ensured longer piRNAs in fine-tuning MIWI/piRNA targeting capacity, proven essential for spermatid development and male fertility.

Cerebral Semaphorin3D is a novel risk factor for age-associated cognitive impairment

BACKGROUND

We previously reported that miR-195 exerts neuroprotection by inhibiting Sema3A and cerebral miR-195 levels decreased with age, both of which urged us to explore the role of miR-195 and miR-195-regulated Sema3 family members in age-associated dementia.

METHODS

miR-195a KO mice were used to assess the effect of miR-195 on aging and cognitive functions. Sema3D was predicted as a miR-195 target by TargetScan and then verified by luciferase reporter assay, while effects of Sema3D and miR-195 on neural senescence were assessed by beta-galactosidase and dendritic spine density. Cerebral Sema3D was over-expressed by lentivirus and suppressed by si-RNA, and effects of over-expression of Sema3D and knockdown of miR-195 on cognitive functions were assessed by Morris Water Maze, Y-maze, and open field test. The effect of Sema3D on lifespan was assessed in Drosophila. Sema3D inhibitor was developed using homology modeling and virtual screening. One-way and two-way repeated measures ANOVA were applied to assess longitudinal data on mouse cognitive tests.

RESULTS

Cognitive impairment and reduced density of dendritic spine were observed in miR-195a knockout mice. Sema3D was identified to be a direct target of miR-195 and a possible contributor to age-associated neurodegeneration as Sema3D levels showed age-dependent increase in rodent brains. Injection of Sema3D-expressing lentivirus caused significant memory deficits while silencing hippocampal Sema3D improved cognition. Repeated injections of Sema3D-expressing lentivirus to elevate cerebral Sema3D for 10 weeks revealed a time-dependent decline of working memory. More importantly, analysis of the data on the Gene Expression Omnibus database showed that Sema3D levels were significantly higher in dementia patients than normal controls (p < 0.001). Over-expression of homolog Sema3D gene in the nervous system of Drosophila reduced locomotor activity and lifespan by 25%. Mechanistically, Sema3D might reduce stemness and number of neural stem cells and potentially disrupt neuronal autophagy. Rapamycin restored density of dendritic spines in the hippocampus from mice injected with Sema3D lentivirus. Our novel small molecule increased viability of Sema3D-treated neurons and might improve autophagy efficiency, which suggested Sema3D could be a potential drug target. Video Abstract CONCLUSION: Our results highlight the importance of Sema3D in age-associated dementia. Sema3D could be a novel drug target for dementia treatment.

Computational Analysis of miR-140 and miR-135 as Potential Targets to Develop Combinatorial Therapeutics for Degenerative Tendinopathy

Background

Degenerative tendinopathy, a condition causing movement restriction due to high pain, highly impacts productivity and quality of life. The healing process is a complex phenomenon and involves a series of intra-cellular and inter-cellular processes. Proliferation and differentiation of the tenocyte is a major and essential process to heal degenerative tendinopathy. The recent development in microRNA (miRNA)-mediated reprogramming of the cellular function through specific pathways opened door for the development of new regenerative therapeutics. Based on information about gene expression and regulation of tendon injury and healing, we attempted to evaluate the combinatorial effect of selected miRNAs for better healing of degenerative tendinopathy.

Methods

The present study was designed to evaluate the combinatorial effect of two miRNAs (has-miR-140 and has-miR-135) in the healing process of the tendon. Publicly available information/data were retrieved from appropriate platforms such as PubMed. Only molecular data, directly associated with tendinopathies, including genes/proteins and miRNAs, were used in this study. The miRNAs involved in tendinopathy were analyzed by a Bioinformatics tools (e.g., TargetScan, miRDB, and the RNA22v2). Interactive involvement of the miRNAs with key proteins involved in tendinopathy was predicted by the Insilco approach.

Results

Based on information available in the public domain, tendon healing-associated miRNAs were predicted to explore their therapeutic potentials. Based on computation analysis, focusing on the potential regulatory effect on tendon healing, the miR-135 and miR-140 were selected for this study. These miRNAs were found as key players in tendon healing through Rho-associated coiled-coil containing protein kinase 1 (ROCK1), IGF-1/PI3K/Akt, PIN, and Wnt signaling pathways. It was also predicted that these miRNAs may reprogram the cells to induce proliferation and differentiation activity. Many miRNAs are likely to regulate genes important for the tendinopathy healing process, and the result of this study allows an approach for miRNA-mediated regeneration of the tenocyte for tendon healing. Based on computational analysis, the role of these miRNAs in different pathways was established, and the results provided insights into the combinatorial approach of miRNA-mediated cell reprogramming.

Conclusions

In this study, the association between miRNAs and the disease was evaluated to correlate the tendinopathy genes and the relevant role of different miRNAs in their regulation. Through this study, it was established that the synergistic effect of more than one miRNA on directed reprogramming of the cell could be helpful in the regeneration of damaged tissue. It is anticipated that this study will be helpful for the design of miRNA cocktails for the orchestration of cellular reprogramming events.

Eriocheir sinensis feminization-1c (Fem-1c) and Its Predicted miRNAs Involved in Sexual Development and Regulation

Feminization-1c (Fem-1c) is important for sex differentiation in the model organism Caenorhabditis elegans. In our previous study, the basic molecular characteristics of the Fem-1c gene (EsFem-1c) in Eriocheir sinensis (Henri Milne Edwards, 1854) were cloned to determine the relationship with sex differentiation. In this study, the genomic sequence of EsFem-1c contained five exons and four introns, with an exceptionally long 3'UTR sequence. The qRT-PCR results of EsFem-1c demonstrated lower tissue expression in the androgenic gland of the intersex crab than the normal male crab, implying that EsFem-1c plays a role in crab AG development. RNA interference experiments and morphological observations of juvenile and mature crabs indicated that EsFem-1c influences sexual development in E. sinensis. A dual-luciferase reporter assay disclosed that tcf-miR-315-5p effectively inhibits the translation of the EsFem-1c gene, influencing male development. An intriguing finding was that miRNA tcf-miR-307 could increase EsFem-1c expression by binding to the alternative splicing region with a length of 248 bp (ASR-248) in the 3'UTR sequence. The present research contributes to a better understanding of the molecular regulation mechanism of EsFem-1c and provides a resource for future studies of the miRNA-mediated regulation of sexual development and regulation in E. sinensis.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Effect of miR-493-5p on proliferation and differentiation of myoblast by targeting ANKRD17

The hypertrophy and conversion of postnatal muscle fibers largely determine the yield and quality of meat, which is closely related to the economic value of pigs. MicroRNA (miRNA), as a kind of endogenous noncoding RNA molecule, is widely involved in myogenesis of livestock and poultry. The longissimus dorsi tissues of Lantang pigs at 1 and 90 days (LT1D and LT90D) were collected and profiled by miRNA-seq. We found 1871 and 1729 miRNA candidates in LT1D and LT90D samples, and 794 miRNAs were shared. We identified 16 differentially expressed miRNAs between two tested groups and explored the function of miR-493-5p inmyogenesis. The miR-493-5p promoted the proliferation and inhibited the differentiation of myoblasts. Using GO and KEGG analyses of 164 target genes of miR-493-5p, we found that ATP2A2, PPP3CA, KLF15, MED28, and ANKRD17 genes were related to muscle development. RT-qPCR detection showed that the expression level of ANKRD17 was highly expressed in LT1D libraries, and the double luciferase report test preliminarily proved that miR-493-5p and ANKRD17 have a directly targeting relationship. We established miRNA profiles for the longissimus dorsi tissues of 1-day-old and 90-day-old Lantang pigs and found that miR-493-5p was differentially expressed and associated with myogenesis by targeting ANKRD17 gene. Our results should serve as a reference for future studies on pork quality.

Dysregulated MicroRNAs in the Pathogenesis of Systemic Lupus Erythematosus: A Comprehensive Review

Systemic lupus erythematosus is a chronic autoimmune disease of which clinical presentation is vastly heterogeneous, ranging from mild skin rashes to severe renal diseases. Treatment goal of this illness is to minimize disease activity and prevent further organ damage. In recent years, much research has been done on the epigenetic aspects of SLE pathogenesis, for among the various factors known to contribute to the pathogenic process, epigenetic factors, especially microRNAs, bear the most therapeutic potential that can be altered unlike congenital genetic factors. This article reviews and updates what has been discovered so far about the pathogenesis of lupus, while focusing on the dysregulation of microRNAs in lupus patients in comparison to healthy controls along with the potentially pathogenic roles of the microRNAs commonly reported to be either upregulated or downregulated. Furthermore, this review includes microRNAs of which results are controversial, suggesting possible explanations for such discrepancies and directions for future research. Moreover, we aimed to emphasize the point that had been overlooked so far in studies regarding microRNA expression levels; that is, which specimen was used to assess the dysregulation of microRNAs. To our surprise, a vast number of studies have not considered this factor and have analyzed the potential role of microRNAs in general. Despite extensive investigations done on microRNA levels, their significance and potential role remain a mystery, which calls for further studies on this particular subject in regard of which specimen is used for assessment.

Identification of CircRNA signature associated with tumor immune infiltration to predict therapeutic efficacy of immunotherapy

Circular RNAs (circRNAs) play important roles in the regulation of cancer. However, the clinical implications and regulatory networks of circRNAs in cancer patients receiving immune checkpoint blockades (ICB) have not been fully elucidated. Here, we characterize circRNA expression profiles in two independent cohorts of 157 ICB-treated advanced melanoma patients and reveal overall overexpression of circRNAs in ICB non-responders in both pre-treatment and early during therapy. Then, we construct circRNA-miRNA-mRNA regulatory networks to reveal circRNA-related signaling pathways in the context of ICB treatment. Further, we construct an ICB-related circRNA signature (ICBcircSig) score model based on progression-free survival-related circRNAs to predict immunotherapy efficacy. Mechanistically, the overexpression of ICBcircSig circTMTC3 and circFAM117B could increase PD-L1 expression via the miR-142-5p/PD-L1 axis, thus reducing T cell activity and leading to immune escape. Overall, our study characterizes circRNA profiles and regulatory networks in ICB-treated patients, and highlights the clinical utility of circRNAs as predictive biomarkers of immunotherapy.

CTLA4 mRNA is downregulated by miR-155 in regulatory T cells, and reduced blood CTLA4 levels are associated with poor prognosis in metastatic melanoma patients

Cytotoxic T lymphocyte-associated antigen-4 (CTLA-4) is an immune checkpoint expressed in regulatory T (Treg) cells and activated T lymphocytes. Despite its potential as a treatment strategy for melanoma, CTLA-4 inhibition has limited efficacy. Using data from The Cancer Genome Atlas (TCGA) melanoma database and another dataset, we found that decreased CTLA4 mRNA was associated with a poorer prognosis in metastatic melanoma. To investigate further, we measured blood CTLA4 mRNA in 273 whole-blood samples from an Australian cohort and found that it was lower in metastatic melanoma than in healthy controls and associated with worse patient survival. We confirmed these findings using Cox proportional hazards model analysis and another cohort from the US. Fractionated blood analysis revealed that Treg cells were responsible for the downregulated CTLA4 in metastatic melanoma patients, which was confirmed by further analysis of published data showing downregulated CTLA-4 surface protein expression in Treg cells of metastatic melanoma compared to healthy donors. Mechanistically, we found that secretomes from human metastatic melanoma cells downregulate CTLA4 mRNA at the post-transcriptional level through miR-155 while upregulating FOXP3 expression in human Treg cells. Functionally, we demonstrated that CTLA4 expression inhibits the proliferation and suppressive function of human Treg cells. Finally, miR-155 was found to be upregulated in Treg cells from metastatic melanoma patients compared to healthy donors. Our study provides new insights into the underlying mechanisms of reduced CTLA4 expression observed in melanoma patients, demonstrating that post-transcriptional silencing of CTLA4 by miRNA-155 in Treg cells may play a critical role. Since CTLA-4 expression is downregulated in non-responder melanoma patients to anti-PD-1 immunotherapy, targeting miRNA-155 or other factors involved in regulating CTLA4 expression in Treg cells without affecting T cells could be a potential strategy to improve the efficacy of immunotherapy in melanoma. Further research is needed to understand the molecular mechanisms regulating CTLA4 expression in Treg cells and identify potential therapeutic targets for enhancing immune-based therapies.

Microgravity alters the expressions of DNA repair genes and their regulatory miRNAs in space-flown Caenorhabditis elegans

During spaceflight, multiple unique hazardous factors, particularly microgravity and space radiation, can induce different types of DNA damage, which pose a constant threat to genomic integrity and stability of living organisms. Although organisms have evolved different kinds of conserved DNA repair pathways to eliminate this DNA damage on Earth, the impact of space microgravity on the expressions of these DNA repair genes and their regulatory miRNAs has not been fully explored. In this study, we integrated all existing datasets, including both transcriptional and miRNA microarrays in wild-type (WT) Caenorhabditis elegans that were exposed to the treatments of spaceflight (SF), spaceflight control with a 1g centrifugal device (SC), and ground control (GC) in three space experiments with the periods of 4, 8 and 16.5 days. The results of principal component analysis showed the gene expression patterns for five major DNA repair pathways (i.e., non-homologous end joining (NHEJ), homologous recombination (HR), mismatch repair (MMR), nucleotide excision repair (NER), and base excision repair (BER)) were well separated and clustered between SF/GC and SC/GC treatments after three spaceflights. In the 16.5-days space experiment, we also selected the datasets of dys-1 mutant and ced-1 mutant of C. elegans, which respectively presented microgravity-insensitivity and radiosensitivity. Compared to the WT C. elegans flown in the 16.5-days spaceflight, the separation distances between SF and SC samples were significantly reduced in the dys-1 mutant, while greatly enhanced in the ced-1 mutant for five DNA repair pathways. By comparing the results of differential expression analysis in SF/GC versus SC/GC samples, we found the DNA repair genes annotated in the pathways of BER and NER were prominently down-regulated under microgravity during both the 4- and 8-days spaceflights. While, under microgravity, the genes annotated in MMR were dominatingly up-regulated during the 4-days spaceflight, and those annotated in HR were mainly up-regulated during the 8-days spaceflight. And, most of the DNA repair genes annotated in the pathways of BER, NER, MMR, and HR were up-regulated under microgravity during the 16.5-days spaceflight. Using miRNA-mRNA integrated analysis, we determined the regulatory networks of differentially expressed DNA repair genes and their regulatory miRNAs in WT C. elegans after three spaceflights. Compared to GC conditions, the differentially expressed miRNAs were analyzed under SF and SC treatments of three spaceflights, and some altered miRNAs that responded to SF and SC could regulate the expressions of corresponding DNA repair genes annotated in different DNA repair pathways. In summary, these findings indicate that microgravity can significantly alter the expression patterns of DNA repair genes and their regulatory miRNAs in space-flown C. elegans. The alterations of the expressions of DNA repair genes and the dominating DNA repair pathways under microgravity are possibly related to the spaceflight period. In addition, the key miRNAs are identified as the post-transcriptional regulators to regulate the expressions of various DNA repair genes under microgravity. These altered miRNAs that responded to microgravity can be implicated in regulating diverse DNA repair processes in space-flown C. elegans.

Exploring microRNAs in craniofacial regenerative medicine

microRNAs (miRs) have been reported over the decades as important regulators in bone development and bone regeneration. They play important roles in maintaining the stem cell signature as well as regulating stem cell fate decisions. Thus, delivering miRs and miR inhibitors to the defect site is a potential treatment towards craniofacial bone defects. However, there are challenges in translation of basic research to clinics, including the efficiency, specificity, and efficacy of miR manipulation methods and the safety of miR delivery systems. In this review, we will compare miR oligonucleotides, mimics and antagomirs as therapeutic reagents to treat disease and regenerate tissues. Newer technology will be discussed as well as the efficiency and efficacy of using these technologies to express or inhibit miRs in treating and repairing oral tissues. Delivery of these molecules using extracellular vesicles and nanoparticles can achieve different results and depending on their composition will elicit specific effects. We will highlight the specificity, toxicity, stability, and effectiveness of several miR systems in regenerative medicine.

Identification and comparative analysis of miRNA transcriptomes after allograft and xenograft transplantation in Pinctada fucata martensii

Effective immune regulation after transplantation during pearl production is crucial for the cultivation of high-quality pearls. MicroRNAs (miRNAs) play an important role in a variety of physiological processes. To understand the regulatory rules of miRNAs after transplantation in Pinctada funcata martensii, we constructed 13 miRNA transcriptomes, including the control group (Con), allograft (Al), and xenograft (Xe) transplantation at six time points (6, 12, and 24 h and 3, 6, and 12 days), in which the xenografted mantle tissue was from Pinctada maxima. We identified 159 differentially expressed miRNAs (DEMs) and found that these DEMs showed high expression at 12 h, 24 h, and 3 days after transplantation. A total of 130 DEMs, such as Let-7, were present in the Al and Xe groups; miR-34 and 16 other DEMs were specifically present in the Al group; miR-216b and 13 other DEMs were specifically present in the Xe group. Compared with the Con group, the target genes of DEMs in the Al group were significantly enriched in protein complex, cytoskeleton, and macromolecular complex, and the Xe group was significantly enriched in ribonucleoside metabolic process, nucleoside binding, and cell division. Compared with the Al group, the target genes in the Xe group were significantly enriched in response to DNA damage stimulation. Overall, multiple pathways associated with cellular activity were enriched in higher numbers of genes in the Xe group than in the Al group. These findings enriched the information on immune regulatory mechanisms at the expression level of miRNAs in P. f. martensii after transplantation.

Global expression and functional analysis of human piRNAs during HSV-1 infection

Piwi-interacting RNAs (piRNAs) are a class of non-coding RNAs that play a key role in spermatogenesis. However, little is known about their expression characterization and role in somatic cells infected with herpes simplex virus type 1 (HSV-1). In this study, we systematically investigated the cellular piRNA expression profiles of HSV-1-infected human lung fibroblasts. Compared with the control group, 69 differentially expressed piRNAs were identified in the infection group, among which 52 were up-regulated and 17 were down-regulated. The changes in the expression of 8 piRNAs were further verified by RT-qPCR with a similar expression trend. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the target genes of piRNAs were mainly involved in antiviral immunity and various human disease-related signaling pathways. Furthermore, we tested the effects of four up-regulated piRNAs on viral replication by transfecting piRNA mimics. The results showed that the virus titers of the group transfected with piRNA-hsa-28,382 (alias piR-36,233) mimic decreased significantly, and that of the group transfected piRNA-hsa-28,190 (alias piR-36,041) mimic significantly increased. Overall, our results revealed the expression characteristics of piRNAs in HSV-1-infected cells. We also screened two piRNAs that potentially regulate HSV-1 replication. These results may promote a better understanding of the regulatory mechanism of pathophysiological changes induced by HSV-1 infection.

Integration of miRNA and mRNA expression profiles in Asian spongy moth Lymantria dispar in response to cyantraniliprole

The Asian spongy moth, Lymantria dispar, is a worldwide forest pest that damages >500 plant species. Nowadays, chemical control is the most widely used method because of its rapidity and effectiveness, but the insecticide resistance is a growing concern for spongy moth. As important post-transcriptional regulators of gene expression, whether microRNAs (miRNAs) are involved in insecticide tolerance is little known in spongy moth. Therefore, an integrated analysis of miRNA and mRNA was performed on L. dispar larvae treated with cyantraniliprole. Compared to the control group, a total of 432 differentially expressed genes (DEGs) and 23 differentially expressed miRNAs (DEMs) were identified in L. dispar larvae under cyantraniliprole exposure. Among them, twelve DEGs encoding detoxification enzymes/proteins were further analyzed. Twenty-one genes related to insecticide tolerance were predicted by 11 DEMs, of which 25 miRNA-mRNA interactions were identified. In the miRNA-mRNA network, novel-miR-4 and mmu-miR-3475-3p were involved in the response of L. dispar to cyantraniliprole stress by regulating five genes associated with detoxification, respectively. The P450 gene CYP4C1 (c34384.graph_c0) was the only DEG related to detoxification in the network, which was regulated by novel-miR-4. The expression levels of ten DEMs were confirmed by quantitative reverse transcription PCR (RT-qPCR) and the trends were consistent with miRNA-seq. This study identified some candidate miRNAs and mRNAs related to cyantraniliprole tolerance in L. dispar, which provides valuable transcriptomic information for revealing the molecular mechanisms of insect tolerance and developing novel insecticides.

Identification and characterization of immune-related microRNAs in hybrid snakehead(Channa maculata♀ × Channa argus♂)after treated by Echinacea purpurea (Linn.) Moench

Echinacea purpurea (Linn.) Moench (EP) is a globally popular herbal medicine, which showed effects on growth promotion, antioxidant and immunomodulatory activities in fish culture world widely. However, there are few studies about the effects on miRNAs by EP in fish. The hybrid snakehead fish (Channa maculate♀ × Channa argus ♂) was new important economic specie of freshwater aquaculture in China with high market value and demand while there were only a few reports about its miRNAs. To overview immune-related miRNAs of the hybrid snakehead fish and to further understand the immune regulating mechanism of EP, we herein constructed and analyzed three small RNA libraries of immune tissues including liver, spleen and head kidney of the fish with or without EP treatment via Illumina high-throughput sequencing technology. Results showed that EP can affect the immune activities of fish by the miRNA-regulated ways. Totally, 67 (47 up and 20 down) miRNAs in liver, 138 (55 up and 83 down) miRNAs in spleen, and 251 (15 up and 236 down) miRNAs in spleen were detected, as well as 30, 60, 139 kinds of immune-related miRNAs belonging to 22, 35 and 66 families of the three tissues respectively. The expressions of 8 immune-related miRNA family members were found in all the three tissues, including miR-10, miR-133, miR-22 and etc. Some miRNAs have been identified involved in the innate and adaptive immune responses, such as the miR-125, miR-138, and miR-181 family. Ten miRNA families with antioxidant target genes were also discovered, including miR-125, miR-1306, and miR-138, etc. Results from Gene Ontology (GO) and KEGG pathway analysis further confirmed there are a majority immune response targets of the miRNAs involved in the EP treatment process. Our study deepened understanding roles of miRNAs in fish immune system and provides new ideas for the study of immune mechanism of EP.

Cardioprotective effects of circ_0002612 in myocardial ischemia/reperfusion injury correlate with disruption of miR-30a-5p-dependent Ppargc1a inhibition

INTRODUCTION

Novel mechanistic insights into the effects of circular RNAs (circRNAs) on the physiology and pathology of cardiovascular diseases are under increasingly active investigation. This study defined the cardioprotective role and mechanistic actions of circ_0002612 in myocardial ischemia/reperfusion injury (MI/RI).

METHODS

MI/RI was induced in mice by ligation of the left anterior descending (LAD) artery followed by reperfusion, and the in vitro model was established in cultured cardiomyocytes under hypoxia/reoxygenation (H/R) conditions. Interaction among circ_0002612, miR-30a-5p, Ppargc1a, and NLRP3 was predicted by bioinformatics analysis and further experimentally identified. Gain- and loss-of-function experiments were performed to evaluate the effect of the circ_0002612/miR-30a-5p/Ppargc1a/NLRP3 axis on the cardiac function and myocardial infarction of I/R-injured mice, as well as viability and apoptosis of H/R-challenged cardiomyocytes.

RESULTS

In the myocardial tissues of MI/RI mice, miR-30a-5p was negatively correlated with circ_0002612 or Ppargc1a, but circ_0002612 was positively correlated with the expression of Ppargc1a. circ_0002612 competitively bound to miR-30a-5p to release expression of its target gene Ppargc1a. circ_0002612 promoted cardiomyocyte viability while suppressing the apoptosis by impairing the miR-30a-5p-mediated inhibition of Ppargc1a. Additionally, Ppargc1a inhibited the expression of NLRP3 and consequently facilitated cardiomyocyte proliferation while suppressing cell apoptosis. By inhibiting the expression of NLRP3, circ_0002612 protected mice from MI/RI.

CONCLUSION

Overall, this study demonstrates the cardioprotective role of circ_0002612 against MI/RI, which may be a viable target for MI/RI.

Competing endogenous RNA network construction based on long non-coding RNAs, microRNAs, and mRNAs related to fat deposition in Songliao black swine

China has a long history of pig breeding and a number of local breeds. The Songliao Black pig, bred in China in 2009, shows high variation in backfat thickness and therefore is well-suited to fat deposition research. Fat deposition is a complex trait, and the underlying regulatory factors are not fully characterized. In this study, the molecular basis of fat deposition traits was evaluated by comparisons between three individuals with extremely high-backfat thickness and three with extremely low-backfat thickness selected from 53 gilts. Subcutaneous adipose tissues of the back were collected for strand-specific library RNA sequencing (RNA-seq) and small RNA-seq. We identified 13 184 mRNAs, 2046 long non-coding (lnc)RNAs, and 494 micro (mi)RNAs by high-throughput sequencing. Furthermore, we detected 150 differentially expressed mRNAs, 66 differentially expressed lncRNAs, and eight differentially expressed miRNAs. A functional enrichment analysis indicated that these genes are involved in multiple fat metabolism-related pathways, including positive regulation of fat cell differentiation, and fat digestion and absorption. We used various algorithms (miRanda, TargetScan, and RNAhybrid) to predict targeting relationships and constructed a competing endogenous RNA network containing seven lncRNAs, three miRNAs, and six mRNAs. All these genes were differentially expressed between the extremely high and low backfat thickness groups or enriched in pathways related to fat metabolism. Our results provide insight into the regulatory mechanisms by which non-coding RNAs and their target genes influence backfat deposition in pigs. Furthermore, our newly constructed competing endogenous RNA (lncRNA-miRNA-mRNA) network provides a basis for further exploration of fat deposition traits and non-coding RNA functions.

Insight into the mechanism of DNA methylation and miRNA-mRNA regulatory network in ischemic stroke

BACKGROUND

Epigenetic changes, such as DNA methylation and miRNA-target gene mechanisms, have recently emerged as key provokers in Ischemic stroke (IS) onset. However, cellular and molecular events harboring these epigenetic alterations are poorly understood. Therefore, the present study aimed to explore the potential biomarkers and therapeutic targets for IS.

METHODS

miRNAs, mRNAs and DNA methylation datasets of IS were derived from the GEO database and normalized by PCA sample analysis. Differentially expressed genes (DEGs) were identified, and GO and KEGG enrichment analyses were performed. The overlapped genes were utilized to construct a protein-protein interaction network (PPI). Meanwhile, differentially expressed mRNAs and miRNAs interaction pairs were obtained from the miRDB, TargetScan, miRanda, miRMap and miTarBase databases. We constructed differential miRNA-target gene regulatory networks based on mRNA-miRNA interactions.

RESULTS

A total of 27 up-regulated and 15 down-regulated differential miRNAs were identified. Dataset analysis identified 1053 and 132 up-regulated and 1294 and 9068 down-regulated differentially expressed genes in the GSE16561 and GSE140275 datasets, respectively. Moreover, 9301 hypermethylated and 3356 hypomethylated differentially methylated sites were also identified. Moreover, DEGs were enriched in terms related to translation, peptide biosynthesis, gene expression, autophagy, Th1 and Th2 cell differentiation, primary immunodeficiency, oxidative phosphorylation and T cell receptor signaling pathway. MRPS9, MRPL22, MRPL32 and RPS15 were identified as hub genes. Finally, a differential miRNA-target gene regulatory network was constructed.

CONCLUSIONS

RPS15, along with hsa-miR-363-3p and hsa-miR-320e have been identified in the differential DNA methylation protein interaction network and miRNA-target gene regulatory network, respectively. These findings strongly posit the differentially expressed miRNAs as potential biomarkers to improve ischemic stroke diagnosis and prognosis.

Erythrocyte miRNA-92a-3p interactions with PfEMP1 as determinants of clinical malaria

Based on the recently added high throughput analysis data on small noncoding RNAs in modulating disease pathophysiology of malaria, we performed an integrative computational analysis for exploring the role of human-host erythrocytic microRNAs (miRNAs) and their influence on parasite survival and host homeostasis. An in silico analysis was performed on transcriptomic datasets accessed from PlasmoDB and Gene Expression Omnibus (GEO) repositories analyzed using miRanda, miRTarBase, mirDIP, and miRDB to identify the candidate miRNAs that were further subjected to network analysis using MCODE and DAVID. This was followed by immune infiltration analysis and screening for RNA degradation mechanisms. Seven erythrocytic miRNAs, miR-451a, miR-92a-3p, miR-16-5p, miR-142-3p, miR-15b-5p, miR-19b-3p, and miR-223-3p showed favourable interactions with parasite genes expressed during blood stage infection. The miR-92a-3p that targeted the virulence gene PfEMP1 showed drastic reduction during infection. Performing pathway analysis for the human-host gene targets for the miRNA identified TOB1, TOB2, CNOT4, and XRN1 genes that are associated to RNA degradation processes, with the exoribonuclease XRN1, highly enriched in the malarial samples. On evaluating the role of exoribonucleases in miRNA degradation further, the pattern of Plasmodium falciparum_XRN1 showed increased levels during infection thus suggesting a defensive role for parasite survival. This study identifies miR-92a-3p, a member of C13orf25/ miR-17-92 cluster, as a novel miRNA inhibitor of the crucial parasite genes responsible for symptomatic malaria. Evidence for a plausible link to chromosome 13q31.3 loci controlling the epigenetic disease regulation is also suggested.

MiRNA-Seq reveals key MicroRNAs involved in fat metabolism of sheep liver

There is a genetic difference between Hu sheep (short/fat-tailed sheep) and Tibetan sheep (short/thin-tailed sheep) in tail type, because of fat metabolism. Previous studies have mainly focused directly on sheep tail fat, which is not the main organ of fat metabolism. The function of miRNAs in sheep liver fat metabolism has not been thoroughly elucidated. In this study, miRNA-Seq was used to identify miRNAs in the liver tissue of three Hu sheep (short/fat-tailed sheep) and three Tibetan sheep (short/thin-tailed sheep) to characterize the differences in fat metabolism of sheep. In our study, Hu sheep was in a control group, we identified 11 differentially expressed miRNAs (DE miRNAs), including six up-regulated miRNAs and five down-regulated miRNAs. Miranda and RNAhybrid were used to predict the target genes of DE miRNAs, obtaining 3,404 target genes. A total of 115 and 67 GO terms as well as 54 and 5 KEGG pathways were significantly (padj &lt; 0.05) enriched for predicted 3,109 target genes of up-regulated and 295 target genes of down-regulated miRNAs, respectively. oar-miR-432 was one of the most up-regulated miRNAs between Hu sheep and Tibetan sheep. And SIRT1 is one of the potential target genes of oar-miR-432. Furthermore, functional validation using the dual-luciferase reporter assay indicated that the up-regulated miRNA; oar-miR-432 potentially targeted sirtuin 1 (SIRT1) expression. Then, the oar-miR-432 mimic transfected into preadipocytes resulted in inhibited expression of SIRT1. This is the first time reported that the expression of SIRT1 gene was regulated by oar-miR-432 in fat metabolism of sheep liver. These results could provide a meaningful theoretical basis for studying the fat metabolism of sheep.

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.

The effect of microRNA on protein variability and gene expression fidelity

Small regulatory RNA molecules such as microRNA modulate gene expression through inhibiting the translation of messenger RNA (mRNA). Such posttranscriptional regulation has been recently hypothesized to reduce the stochastic variability of gene expression around average levels. Here, we quantify noise in stochastic gene expression models with and without such regulation. Our results suggest that silencing mRNA posttranscriptionally will always increase, rather than decrease, gene expression noise when the silencing of mRNA also increases its degradation, as is expected for microRNA interactions with mRNA. In that regime, we also find that silencing mRNA generally reduces the fidelity of signal transmission from deterministically varying upstream factors to protein levels. These findings suggest that microRNA binding to mRNA does not generically confer precision to protein expression.

Insights into the mechanisms of triptolide nephrotoxicity through network pharmacology-based analysis and RNA-seq

INTRODUCTION

Triptolide (TPL) is a promising plant-derived compound for clinical therapy of multiple human diseases; however, its application was limited considering its toxicity.

METHODS

To explore the underlying molecular mechanism of TPL nephrotoxicity, a network pharmacology based approach was utilized to predict candidate targets related with TPL toxicity, followed by deep RNA-seq analysis to characterize the features of three transcriptional elements include protein coding genes (PCGs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) as well as their associations with nephrotoxicity in rats with TPL treatment.

RESULTS & DISCUSSION

Although the deeper mechanisms of TPL nephrotoxcity remain further exploration, our results suggested that c-Jun is a potential target of TPL and Per1 related circadian rhythm signaling is involved in TPL induced renal toxicity.

Trypanosoma cruzi dysregulates expression profile of piRNAs in primary human cardiac fibroblasts during early infection phase

Trypanosoma cruzi, the etiological agent of Chagas Disease, causes severe morbidity, mortality, and economic burden worldwide. Though originally endemic to Central and South America, globalization has led to increased parasite presence in most industrialized countries. About 40% of infected individuals will develop cardiovascular, neurological, and/or gastrointestinal pathologies. Accumulating evidence suggests that the parasite induces alterations in host gene expression profiles in order to facilitate infection and pathogenesis. The role of regulatory gene expression machinery during T. cruzi infection, particularly small noncoding RNAs, has yet to be elucidated. In this study, we aim to evaluate dysregulation of a class of sncRNAs called piRNAs during early phase of T. cruzi infection in primary human cardiac fibroblasts by RNA-Seq. We subsequently performed in silico analysis to predict piRNA-mRNA interactions. We validated the expression of these selected piRNAs and their targets during early parasite infection phase by stem loop qPCR and qPCR, respectively. We found about 26,496,863 clean reads (92.72%) which mapped to the human reference genome. During parasite challenge, 441 unique piRNAs were differentially expressed. Of these differentially expressed piRNAs, 29 were known and 412 were novel. In silico analysis showed several of these piRNAs were computationally predicted to target and potentially regulate expression of genes including SMAD2, EGR1, ICAM1, CX3CL1, and CXCR2, which have been implicated in parasite infection, pathogenesis, and various cardiomyopathies. Further evaluation of the function of these individual piRNAs in gene regulation and expression will enhance our understanding of early molecular mechanisms contributing to infection and pathogenesis. Our findings here suggest that piRNAs play important roles in infectious disease pathogenesis and can serve as potential biomarkers and therapeutic targets.

Bioactive glass-elicited stem cell-derived extracellular vesicles regulate M2 macrophage polarization and angiogenesis to improve tendon regeneration and functional recovery

Effective countermeasures for tendon injury remains unsatisfactory. Mesenchymal stem cell-derived extracellular vesicles (MSC-EVs)-based therapy via regulation of Mφ-mediated angiogenesis has emerged as a promising strategy for tissue regeneration. Still, approaches to tailor the functions of EVs to treat tendon injuries have been limited. We reported a novel strategy by applying MSC-EVs boosted with bioactive glasses (BG). BG-elicited EVs (EV_B) showed up-regulation of medicinal miRNAs, including miR-199b-3p and miR-125a-5p, which play a pivotal role in M2 Mφ-mediated angiogenesis. EV_B accelerated angiogenesis via the reprogrammed anti-inflammatory M2 Mφs compared with naïve MSC-EVs (EV_N). In rodent Achilles tendon rupture model, EV_B local administration activated anti-inflammatory responses via M2 polarization and led to a spatial correlation between M2 Mφs and newly formed blood vessels. Our results showed that EV_B outperformed EV_N in promoting tenogenesis and in reducing detrimental morphological changes without causing heterotopic ossification. Biomechanical test revealed that EV_B significantly improved ultimate load, stiffness, and tensile modulus of the repaired tendon, along with a positive correlation between M2/M1 ratio and biomechanical properties. On the basis of the boosted nature to reprogram regenerative microenvironment, EV_B holds considerable potential to be developed as a next-generation therapeutic modality for enhancing functional regeneration to achieve satisfying tendon regeneration.

Microarray profiling of circular RNAs in LPS-induced Beas-2B cells and identification of circ_2979 as a potential biomarker for ARDS

Acute respiratory distress syndrome (ARDS) is a major cause of high mortality and morbidity in critically ill patients. Circular RNAs (CircRNAs) are widely expressed in numerous tissues and are associated with various diseases. However, the role of circRNAs in ARDS remains unclear. In this study, we found that cell viability and proliferation were reduced in lipopolysaccharide (LPS)-induced Beas-2B cells. Microarray analysis identified 1131 differentially expressed circRNAs in LPS-treated Beas-2B cells, with 623 circRNAs significantly upregulated and 508 circRNAs strongly downregulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed significant enrichment and indicated potential functions and pathways of differentially expressed circRNAs. Reverse transcription-polymerase chain reaction (RT-PCR) analysis confirmed that expression of circ_2979, circ_5438, circ_4557 and circ_2066 in LPS-induced Beas-2B cells was consistent with the results obtained by RNA sequencing (RNA-seq). Additionally, we recruited 17 patients with ARDS and 13 healthy volunteers and found that expression of circ_2979 in serum was significantly increased in the patients with ARDS compared with healthy volunteers. Spearman's analyses indicated that circ_2979 was correlated with partial pressure of carbon dioxide in arterial blood (PaCO₂), the ratio of partial pressure of arterial oxygen to the fraction of inspired oxygen (PaO₂/FiO₂), interleukin 2 receptor (IL-2R) and fibrinogen (FIB). The results suggested that circRNAs may play an important role in the progression of ARDS, and that circ_2979 may serve as a diagnosis and prognosis biomarker for ARDS.

Elucidation of transfer RNAs as stress regulating agents and the experimental strategies to conceive the functional role of tRNA-derived fragments in plants

In plants, the transfer RNAs (tRNAs) exhibit their profound influence in orchestrating diverse physiological activities like cell growth, development, and response to several surrounding stimuli. The tRNAs, which were known to restrict their function solely in deciphering the codons, are now emerging as frontline defenders in stress biology. The plants that are constantly confronted with a huge panoply of stresses rely on tRNA-mediated stress regulation by altering the tRNA abundance, curbing the transport of tRNAs, fragmenting the mature tRNAs during stress. Among them, the studies on the generation of transfer RNA-derived fragments (tRFs) and their biological implication in stress response have attained huge interest. In plants, the tRFs hold stable expression patterns and regulate biological functions under diverse environmental conditions. In this review, we discuss the fate of plant tRNAs upon stress and thereafter how the tRFs are metamorphosed into sharp ammunition to wrestle with stress. We also address the various methods developed to date for uncovering the role of tRFs and their function in plants.

Green Tea Epigallocatechin Gallate Inhibits Preadipocyte Growth via the microRNA-let-7a/HMGA2 Signaling Pathway

SCOPE

This study investigates the effect of epigallocatechin gallate (EGCG) on white and beige preadipocyte growth and explores the involvement of the miR-let-7a/HMGA2 pathway.

METHODS AND RESULTS

3T3-L1 and D12 cells are treated with EGCG. The effect of EGCG on cell proliferation and viability is evaluated, as well as microRNA (miRNA)-related signaling pathways. EGCG inhibits 3T3-L1 and D12 preadipocyte growth, upregulates miR-let-7a expression, and downregulates high-mobility group AT-hook 2 (HMGA2) mRNA and protein levels in a time- and dose-dependent manner. In addition, overexpression of miR-let-7a significantly inhibits the growth of 3T3-L1 and D12 cells and decreases HMGA2 mRNA and protein levels. MiR-let-7a inhibitor antagonizes the inhibitory effects of EGCG on the number and viability of 3T3-L1 and D12 cells. Furthermore, miR-let-7a inhibitor reverses the EGCG-induced increase in miR-let-7a expression levels and decrease in HMGA2 mRNA and protein levels. HMGA2 overexpression induces an increase in cell number and viability and antagonizes EGCG-suppressed cell growth and HMGA2 expression in 3T3-L1 and D12 preadipocytes.

CONCLUSION

EGCG inhibits the growth of 3T3-L1 and D12 preadipocytes by modulating the miR-let-7a and HMGA2 pathways.

Systems Biology Approaches for the Improvement of Oncolytic Virus-Based Immunotherapies

Oncolytic virus (OV)-based immunotherapy is mainly dependent on establishing an efficient cell-mediated antitumor immunity. OV-mediated antitumor immunity elicits a renewed antitumor reactivity, stimulating a T-cell response against tumor-associated antigens (TAAs) and recruiting natural killer cells within the tumor microenvironment (TME). Despite the fact that OVs are unspecific cancer vaccine platforms, to further enhance antitumor immunity, it is crucial to identify the potentially immunogenic T-cell restricted TAAs, the main key orchestrators in evoking a specific and durable cytotoxic T-cell response. Today, innovative approaches derived from systems biology are exploited to improve target discovery in several types of cancer and to identify the MHC-I and II restricted peptide repertoire recognized by T-cells. Using specific computation pipelines, it is possible to select the best tumor peptide candidates that can be efficiently vectorized and delivered by numerous OV-based platforms, in order to reinforce anticancer immune responses. Beyond the identification of TAAs, system biology can also support the engineering of OVs with improved oncotropism to reduce toxicity and maintain a sufficient portion of the wild-type virus virulence. Finally, these technologies can also pave the way towards a more rational design of armed OVs where a transgene of interest can be delivered to TME to develop an intratumoral gene therapy to enhance specific immune stimuli.

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

In silico analysis reveals hypoxia-induced miR-210-3p specifically targets SARS-CoV-2 RNA

Human coronaviruses (HCoVs) until the emergence of SARS in 2003 were associated with mild cold and upper respiratory tract infections. The ongoing pandemic caused by SARS-CoV-2 has enhanced the potential for infection and transmission as compared to other known members of this family. MicroRNAs (miRNA) are 21-25 nucleotides long non-coding RNA that bind to 3' UTR of genes and regulate almost every aspect of cellular function. Several human miRNAs have been known to target viral genomes, mostly to downregulate their expression and sometimes to upregulate also. In some cases, host miRNAs could be sequestered by the viral genome to create a condition for favourable virus existence. The ongoing SARS CoV-2 pandemic is unique based on its transmissibility and severity and we hypothesised that there could be a unique mechanism for its pathogenesis. In this study, we exploited in silico approach to identify human respiratory system-specific miRNAs targeting the viral genome of three highly pathogenic HCoVs (SARS-CoV-2 Wuhan strain, SARS-CoV, and MERS-CoV) and three low pathogenic HCoVs (OC43, NL63, and HKU1). We identified ten common microRNAs that target all HCoVs studied here. In addition, we identified unique miRNAs which targeted specifically one particular HCoV. miR-210-3p was the single unique lung-specific miRNA, which was found to target the NSP3, NSP4, and NSP13 genes of SARS-CoV-2. Further miR-210-NSP3, miR-210-NSP4, and miR-210-NSP13 SARS-CoV-2 duplexes were docked with the hAGO2 protein (PDB ID 4F3T) which showed Z-score values of -1.9, -1.7, and -1.6, respectively. The role of miR-210-3p as master hypoxia regulator and inflammation regulation may be important for SARS-CoV-2 pathogenesis. Overall, this analysis advocates that miR-210-3p be investigated experimentally in SARS-CoV-2 infection.Communicated by Ramaswamy H. Sarma.

Long noncoding RNA TARL promotes antibacterial activity and prevents bacterial escape in Miichthys miiuy through suppression of TAK1 downregulation

Noncoding RNA (ncRNA) is an important regulatory factor that plays a major role in innate immunity. However, most studies on ncRNA have focused on mammals, resulting in a knowledge gap on ncRNA in lower vertebrates such as teleost fish. In this study, we identified a new long noncoding RNA (lncRNA), termed TAK1-related lncRNA (TARL), which can play a positive role in the antibacterial immunity of Miichthys miiuy to Vibrio anguillarum and V. harveyi. We also found a novel microRNA miR-2188-3p that could target TAK1 and inhibit the host antibacterial response and promote bacterial escape. We further found that the antibacterial effect inhibited by miR-2188-3p could be reversed with TARL. Moreover, V. anguillarum and V. harveyi are the two most susceptible Gram-negative pathogens of aquaculture fish, and the economic losses caused by these two bacteria are immeasurable every year. This study is the first to report on the ability of lncRNA to prevent the escape of V. anguillarum and V. harveyi in fish through the competing endogenous RNA (ceRNA) mechanism. Our results not only elucidate the ceRNA mechanism of the lncRNA in antibacterial immune responses but also provide new insights into the impact of lncRNA on host immunity and bacterial escape.

Uterine Flushing Fluid-Derived Let-7b Targets CXCL10 to Regulate Uterine Receptivity in Goats during Embryo Implantation

Exosomes have the ability to carry a wide range of chemicals, convey them to target cells or target regions, and act as "messengers." For the purpose of investigating embryo attachment, it is helpful to comprehend the range of exosomal mRNAs and miRNAs derived from the uterine flushing fluid before and after embryo attachment. In this study, we recovered exosomes from goat uterine rinsing fluid at 5, 15, and 18 days of gestation and used RNA-Seq to identify the mRNA and miRNA profiles of exosomes obtained from uterine rinsing fluid before and after embryo implantation. In total, 91 differently expressed miRNAs and 27,487 differentially expressed mRNAs were found. The target genes predicted by the differentially expressed miRNAs and the differentially expressed mRNAs were mainly membrane-related organelles with catalytic activity, binding activity, transcriptional regulation activity, and involved in metabolism, biological regulation, development, and other processes. This was revealed by GO analysis. Furthermore, KEGG analysis revealed that they were abundant in signaling pathways associated with embryo implantation, including the "PI3K-Akt signaling pathway," "Toll-like receptor signaling pathway," "TGF-beta signaling route," "Notch signaling pathway," and others. Moreover, our research has demonstrated, for the first time, that chi-let-7b-5p specifically targets the 3'UTR of CXCL10. Our research offers a fresh viewpoint on the mechanics of embryo attachment.

Integrative analysis of small non-coding RNAs predicts a piRNA/miRNA-CCND1/BRAF/HRH1/ATXN3 regulatory circuit that drives oncogenesis in glioblastoma

The high-grade astrocytoma, glioblastoma multiforme (GBM), is the most common primary tumour of the brain, known for being aggressive and developing drug resistance. The non-coding RNAs (ncRNAs), such as microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs), have critical functions in tumorigenesis and cancer drug resistance. Hence, we profiled miRNAs, piRNAs, and genes in U-87 MG GBM cells by next-generation sequencing and performed target prediction, pathway enrichment, protein-protein interaction, co-expression studies, and qRT-PCR validations to predict their possible roles in the malignancy. The study identified 335 miRNAs, 665 piRNAs, and 4286 genes differentially expressed (DE) in GBM. Among them 128 DE genes (DEGs) were targeted by both miRNAs and piRNAs, while 1817 and 192 were targeted solely by miRNAs or piRNAs, respectively. Interestingly, all the DEG targets enriched in cancer processes were overexpressed in GBM. Among these, BRAF was solely targeted by two piRNAs and this was found to be co-expressed with 19 sole targets of 5 miRNAs, including CCND1, and both were found to regulate cell proliferation in cancer. We conjectured that upregulated HRH1 and ATXN3 were targeted by both piRNAs and miRNAs, and along with BRAF and CCND1 might induce cell proliferation in GBM through G-protein-coupled receptor or Akt signalling pathways due to downregulation of the respective targeting small RNAs. These targets were also linked to the progression and overall survival of GBM patients, suggesting that they could be used as biomarkers. Overall, this study has identified a few novel ncRNA targets, which might aid in a better understanding of GBM pathogenesis.

Transcription analyses of differentially expressed mRNAs, lncRNAs, circRNAs, and miRNAs in the growth plate of rats with glucocorticoid-induced growth retardation

Background

Glucocorticoids (GCs) are commonly used to treat autoimmune diseases and malignancies in children and adolescents. Growth retardation is a common adverse effect of GC treatment in pediatric patients. Accumulating evidence indicates that non-coding RNAs (ncRNAs) are involved in the pathogenesis of glucocorticoid-induced growth retardation (GIGR), but the roles of specific ncRNAs in growth remain largely unknown.

Methods

In this study, 2-week-old male Sprague-Dawley rats had been treated with 2 mg/kg/d of dexamethasone for 7 or 14 days, after which the growth plate tissues were collected for high-throughput RNA sequencing to identify differentially expressed mRNAs, lncRNAs, circRNAs, and miRNAs in GIGR rats.

Results

Transcriptomic analysis identified 1,718 mRNAs, 896 lncRNAs, 60 circRNAs, and 72 miRNAs with different expression levels in the 7d group. In the 14d group, 1,515 mRNAs, 880 lncRNAs, 46 circRNAs, and 55 miRNAs with differential expression were identified. Four mRNAs and four miRNAs that may be closely associated with the development of GIGR were further validated by real-time quantitative fluorescence PCR. Function enrichment analysis indicated that the PI3K-Akt signaling pathway, NF-kappa B signaling pathway, and TGF-β signaling pathway participated in the development of the GIGR. Moreover, the constructed ceRNA networks suggested that several miRNAs (including miR-140-3p and miR-127-3p) might play an important role in the pathogenesis of GIGR.

Conclusions

These results provide new insights and important clues for exploring the molecular mechanisms underlying GIGR.

Altered miRNA Expression Profiles in the Serum of Beagle Dogs Experimentally Infected with Toxocara canis

Toxocara canis is a neglected roundworm, which can cause debilitating disease in dogs and humans worldwide. Serum is an excellent material for monitoring the occurrence of many diseases. However, no information is available on the expression of microRNAs (miRNAs) in the serum of dogs infected with T. canis. In this study, RNA-seq analysis was performed to identify the serum miRNA profiles in Beagle dogs infected with T. canis at different stages of infection. A total of 3, 25 and 25 differently expressed miRNAs (DEmiRNAs) were identified in dog serum at 24 h post-infection (hpi), 10 days post-infection (dpi) and 36 dpi, respectively, such as cfa-let-7g, cfa-miR-16, cfa-miR-92b, cfa-miR-93, cfa-miR-122, cfa-miR-485 and cfa-miR-451. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that these miRNAs could regulate the pathways related to parasitic infectious diseases and immune system, such as amoebiasis, toxoplasmosis, platelet activation, IL-17 signaling pathway and chemokine signaling pathway. These results provide a foundation to explore the underlying regulatory role of miRNAs in definitive hosts after T. canis infection.

Regulation of HIF-1 by MicroRNAs in Various Cardiovascular Diseases

Today, we see an increase in death due to cardiovascular diseases all over the world, which has a lot to do with the regulation of oxygen homeostasis. Also, hypoxia-inducing factor 1 (HIF-1) is considered a vital factor in hypoxia and its physiological and pathological changes. HIF- 1 is involved in cellular activities, including proliferation, differentiation, and cell death in endothelial cells (ECs) and cardiomyocytes. Similar to HIF-1α, which acts as a protective element against various diseases in the cardiovascular system, the protective role of microRNAs (miRNAs) has also been proved using animal models. The number of miRNAs identified in the regulation of gene expression responsive to hypoxia and the importance of investigating the involvement of the non-coding genome in cardiovascular diseases is increasing, which shows the issue's importance. In this study, the molecular regulation of HIF-1 by miRNAs is considered to improve therapeutic approaches in clinical diagnoses of cardiovascular diseases.

miR-2188-5p promotes GCRV replication by the targeted degradation of klf2a in Ctenopharyngodon idellus

Studies on host immunity evasion by aquatic viruses have largely focused on coding genes. There is accumulating evidence for the important biological functions of non-coding miRNAs in virus-host interactions. The regulatory functions of non-coding miRNAs in fish reovirus-host interactions remain unknown. Here, miR-2188-5p in grass carp (Ctenopharyngodon idellus), a miRNA specific to teleosts, was predicted to target the 3' UTR of the transcription factor klf2a. A correlation analysis and dual-luciferase reporter assay revealed that miR-2188-5p could induce the degradation of klf2a. The expression of miR-2188-5p induced the degradation of klf2a in a dose-dependent manner, suppressing the type I interferon response and promoting grass carp reovirus (GCRV) replication. As determined by a co-expression analysis, klf2a inhibited viral infection when miR-2188-5p was overexpressed. The targeted degradation of klf2a by miR-2188-5p could inhibit the type I interferon response and promote the replication of GCRV; however, this targeted degradation ability was insufficient to fully inhibit GCRV infection. These results provide novel insights into the regulatory effects and biological functions of non-coding miRNAs in fish-virus interactions.

miR-92b-5p negatively regulates IKK through targeting its ORF region in the innate immune responses of amphioxus (Branchiostoma belcheri)

Cephalochordate (Amphioxus), situated at a key phylogenetic position in the phylum Chordata, serves as a model organism for studying the origin and evolution of the vertebrate innate immune. In this study, five members of precursor miR-92 family (miR-92a-1, miR-92a-2, miR-92b, miR-92c and miR-92d) are identified in Branchiostoma belcheri, and their evolutionary conservation and potential molecular functions in innate immunity are analyzed. Among them, miR-92b-5p was validated in HEK293 cells to target the coding region but not classic 3'UTR of IKK (inhibitor of nuclear factor kappa-B kinase) mRNA, one integral component of MAPK and TLR4 immune signaling. Furthermore, the spatiotemporal expression patterns of miR-92b-5p and IKK were examined in different tissues or different time points (2 h, 4 h, 8 h, 12 h, 24 h and 48 h) post LPS stimulation at RNA and protein level in vivo. The seemingly inverse expression pattern between miR-92b-5p and IKK supports the involvement of miR-92b-5p in Branchiostoma belcheri innate immune response. In conclusion, our work not only illustrates the evolutionary pattern of Branchiostoma belcheri miR-92 family across chordates, but also reveals that miR-92b-5p could target IKK expression to regulate innate immune response.

DNA hypermethylation-induced miR-182 silence targets BCL2 and HOXA9 to facilitate the self-renewal of leukemia stem cell, accelerate acute myeloid leukemia progression, and determine the sensitivity of BCL2 inhibitor venetoclax

Rationale: microRNAs (miRNAs) are frequently deregulated and play important roles in the pathogenesis and progression of acute myeloid leukemia (AML). miR-182 functions as an onco-miRNA or tumor suppressor miRNA in the context of different cancers. However, whether miR-182 affects the self-renewal of leukemia stem cells (LSCs) and normal hematopoietic stem progenitor cells (HSPCs) is unknown. Methods: Bisulfite sequencing was used to analyze the methylation status at pri-miR-182 promoter. Lineage-negative HSPCs were isolated from miR-182 knockout (182KO) and wild-type (182WT) mice to construct MLL-AF9-transformed AML model. The effects of miR-182 depletion on the overall survival and function of LSC were analyzed in this mouse model in vivo. Results: miR-182-5p (miR-182) expression was lower in AML blasts than normal controls (NCs) with hypermethylation observed at putative pri-miR-182 promoter in AML blasts but unmethylation in NCs. Overexpression of miR-182 inhibited proliferation, reduced colony formation, and induced apoptosis in leukemic cells. In addition, depletion of miR-182 accelerated the development and shortened the overall survival (OS) in MLL-AF9-transformed murine AML through increasing LSC frequency and self-renewal ability. Consistently, overexpression of miR-182 attenuated AML development and extended the OS in the murine AML model. Most importantly, miR-182 was likely dispensable for normal hematopoiesis. Mechanistically, we identified BCL2 and HOXA9 as two key targets of miR-182 in this context. Most importantly, AML patients with miR-182 unmethylation had high expression of miR-182 followed by low protein expression of BCL2 and resistance to BCL2 inhibitor venetoclax (Ven) in vitro. Conclusions: Our results suggest that miR-182 is a potential therapeutic target for AML patients through attenuating the self-renewal of LSC but not HSPC. miR-182 promoter methylation could determine the sensitivity of Ven treatment and provide a potential biomarker for it.

Exosomal microRNAs regulate the heat stress response in sea cucumber Apostichopus japonicus

Exosomes are small extracellular vesicles that contain nucleic acids such as microRNAs and may participate in important biological processes. We made the initial report of exosomes from sea cucumber Apostichopus japonicus, that were classically cup-shaped and had an average size of 74.65 nm, and identified specific exosome biomarkers (HSP70, TSG101, and CD9). We explored changes in the global expression of microRNAs in exosomes from the commercially important A. japonicus under normal conditions and heat-stressed conditions for 3 and 7 d. We found that heat stress increased exosome production and modified the expression profiles of the microRNAs that they contained. Novel_mir31, novel_mir132, novel_mir26, miR-92_1, and novel_mir27 were commonly found to be differentially expressed in three comparison groups, indicating their importance in the heat stress response. The microRNA expression levels were validated by qPCR. Function analysis of the target genes of these microRNAs indicated they were involved mainly in replication and repair in the initial response of A. japonicus to heat stress exposure. Conversely, during acclimation to the high temperature conditions, the target genes of the differentially expressed microRNAs were primarily involved in metabolism adjustments. Our results will contribute to a better understanding of the regulatory roles of exosomes in sea cucumber, and provide insights into the functions of sea cucumber exosome-shuttled microRNAs against environmental stresses exacerbated by global warming.

A combined miRNA-piRNA signature in the serum and urine of rabbits infected with Toxoplasma gondii oocysts

BACKGROUND

Increasing evidence has shown that non-coding RNA (ncRNA) molecules play fundamental roles in cells, and many are stable in body fluids as circulating RNAs. Study on these ncRNAs will provide insights into toxoplasmosis pathophysiology and/or help reveal diagnostic biomarkers.

METHODS

We performed a high-throughput RNA-Seq study to comprehensively profile the microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) in rabbit serum and urine after infection with Toxoplasma gondii oocysts during the whole infection process.

RESULTS

Total RNA extracted from serum and urine samples of acutely infected [8 days post-infection (DPI)], chronically infected (70 DPI) and uninfected rabbits were subjected to genome-wide small RNA sequencing. We identified 2089 miRNAs and 2224 novel piRNAs from the rabbit sera associated with T. gondii infection. Meanwhile, a total of 518 miRNAs and 4182 novel piRNAs were identified in the rabbit urine associated with T. gondii infection. Of these identified small ncRNAs, 1178 and 1317 serum miRNAs and 311 and 294 urine miRNAs were identified as differentially expressed (DE) miRNAs in the acute and chronic stages of infections, respectively. A total of 1748 and 1814 serum piRNAs and 597 and 708 urine piRNAs were found in the acute and chronic infection stages, respectively. Of these dysregulated ncRNAs, a total of 88 common DE miRNAs and 120 DE novel piRNAs were found in both serum and urine samples of infected rabbits.

CONCLUSIONS

These findings provide valuable data for revealing the physiology of herbivore toxoplasmosis caused by oocyst infection. Circulating ncRNAs identified in this study are potential novel diagnostic biomarkers for the detection/diagnosis of toxoplasmosis in herbivorous animals.

Identification of Circular RNAs of Testis and Caput Epididymis and Prediction of Their Potential Functional Roles in Donkeys

Circular RNAs (circRNAs) are a class of noncoding RNAs with a covalently closed loop. Studies have demonstrated that circRNA can function as microRNA (miRNA) sponges or competing endogenous RNAs. Although circRNA has been explored in some species and tissues, the genetic basis of testis development and spermatogenesis in donkeys remain unknown. We performed RNA-seq to detect circRNA expression profiles of adult donkey testes. Length distribution and other characteristics were shown a total of 1971 circRNAs were differentially expressed and 12,648 and 6261 circRNAs were detected from the testis and caput epididymis, respectively. Among these circRNAs, 1472 circRNAs were downregulated and 499 circRNAs were upregulated in the testis. Moreover, KEGG pathway analyses and Gene Ontology were performed for host genes of circRNAs. A total of 39 upregulated circRNA host genes were annotated in spermatogenesis terms, including PIWIL2, CATSPERD, CATSPERB, SPATA6, and SYCP1. Other host genes were annotated in the focal adhesion, Rap1 signaling pathway. Downregulated expressed circRNA host genes participated in the TGF-β signaling pathway, GnRH signaling pathway, estrogen signaling pathway, and calcium signaling pathway. Our discoveries provide a solid foundation for identifying and characterizing critical circRNAs involved in testis development or spermatogenesis.

The Expression of miRNAs Involved in Long-Term Memory Formation in the CNS of the Mollusk Helix lucorum

Mollusks are unique animals with a relatively simple central nervous system (CNS) containing giant neurons with identified functions. With such simple CNS, mollusks yet display sufficiently complex behavior, thus ideal for various studies of behavioral processes, including long-term memory (LTM) formation. For our research, we use the formation of the fear avoidance reflex in the terrestrial mollusk Helix lucorum as a learning model. We have shown previously that LTM formation in Helix requires epigenetic modifications of histones leading to both activation and inactivation of the specific genes. It is known that microRNAs (miRNAs) negatively regulate the expression of genes; however, the role of miRNAs in behavioral regulation has been poorly investigated. Currently, there is no miRNAs sequencing data being published on Helix lucorum, which makes it impossible to investigate the role of miRNAs in the memory formation of this mollusk. In this study, we have performed sequencing and comparative bioinformatics analysis of the miRNAs from the CNS of Helix lucorum. We have identified 95 different microRNAs, including microRNAs belonging to the MIR-9, MIR-10, MIR-22, MIR-124, MIR-137, and MIR-153 families, known to be involved in various CNS processes of vertebrates and other species, particularly, in the fear behavior and LTM. We have shown that in the CNS of Helix lucorum MIR-10 family (26 miRNAs) is the most representative one, including Hlu-Mir-10-S5-5p and Hlu-Mir-10-S9-5p as top hits. Moreover, we have shown the involvement of the MIR-10 family in LTM formation in Helix. The expression of 17 representatives of MIR-10 differentially changes during different periods of LTM consolidation in the CNS of Helix. In addition, using comparative analysis of microRNA expression upon learning in normal snails and snails with deficient learning abilities with dysfunction of the serotonergic system, we identified a number of microRNAs from several families, including MIR-10, which expression changes only in normal animals. The obtained data can be used for further fundamental and applied behavioral research.