microRNA-mRNA Profile of Skeletal Muscle Differentiation and Relevance to Congenital Myotonic Dystrophy

microRNAs (miRNAs) regulate messenger RNA (mRNA) abundance and translation during key developmental processes including muscle differentiation. Assessment of miRNA targets can provide insight into muscle biology and gene expression profiles altered by disease. mRNA and miRNA libraries were generated from C2C12 myoblasts during differentiation, and predicted miRNA targets were identified based on presence of miRNA binding sites and reciprocal expression. Seventeen miRNAs were differentially expressed at all time intervals (comparing days 0, 2, and 5) of differentiation. mRNA targets of differentially expressed miRNAs were enriched for functions related to calcium signaling and sarcomere formation. To evaluate this relationship in a disease state, we evaluated the miRNAs differentially expressed in human congenital myotonic dystrophy (CMD) myoblasts and compared with normal control. Seventy-four miRNAs were differentially expressed during healthy human myocyte maturation, of which only 12 were also up- or downregulated in CMD patient cells. The 62 miRNAs that were only differentially expressed in healthy cells were compared with differentiating C2C12 cells. Eighteen of the 62 were conserved in mouse and up- or down-regulated during mouse myoblast differentiation, and their C2C12 targets were enriched for functions related to muscle differentiation and contraction.     

Characterization of MicroRNAs Associated with Reproduction in the Brown Planthopper,Nilaparvata lugens

Insects have a robust capacity to produce offspring for propagation, and the reproductive events of female insects have been achieved at the molecular and physiological levels via regulatory gene pathways. However, the roles of MicroRNAs (miRNAs) in the reproductive development of the brown planthopper (BPH), Nilaparvata lugens, remain largely unexplored. To understand the roles of miRNAs in reproductive development, miRNAs were identified by Solexa sequencing in short-winged (SW) female adults of BPH. Small RNA libraries derived from three developmental phases (1 day, 3 days, and 5 days after emergence) were constructed and sequenced. We identified 905 miRNAs, including 263 known and 642 novel miRNAs. Among them, a total of 43 miRNAs were differentially expressed in the three developmental phases, and 14,568 putative targets for 43 differentially expressed miRNAs (DEMs) were predicted by TargetScan and miRanda. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the predicted miRNA targets illustrated the putative roles for these DEMs in reproduction. The progress events were annotated, including oogenesis, lipid biosynthetic process, and related pathways such as apoptosis, ABC transporters, and amino acid metabolism. Four highly abundant DEMs (miR-9a-5p, miR-34-5p, miR-275-3p, and miR-317-3p) were further screened, and miR-34-5p was confirmed to be involved in the regulation of reproduction. Overexpression of miR-34-5p via injecting its mimics reduced fecundity and decreased Vg expression. Moreover, target genes prediction for miR-34-5p showed they might be involved in 20E signaling cascades, apoptosis, and gonadal development, including hormone receptor 4 (HR4), caspase-1 (Cp-1), and spermatogenesis-associated protein 20 (SPATA20). These findings provide a valuable resource for future studies on the role of miRNAs in BPH reproductive development.     

MicroRNA Expression Profiling of Lactating Mammary Gland in Divergent Phenotype Swine Breeds

MicroRNA (miRNA) plays a key role in development and specific biological processes, such as cell proliferation, differentiation, and apoptosis. Extensive studies of mammary miRNAs have been performed in different species and tissues. However, little is known about porcine mammary gland miRNAs. In this study, we report the identification and characterization of miRNAs in the lactating mammary gland in two distinct pig breeds, Jinhua and Yorkshire. Many miRNAs were detected as significantly differentially expressed between the two libraries. Among the differentially expressed miRNAs, many are known to be related to mammary gland development and lactation by interacting with putative target genes in previous studies. These findings suggest that miRNA expression patterns may contribute significantly to target mRNA regulation and influence mammary gland development and peak lactation performance. The data we obtained provide useful information about the roles of miRNAs in the biological processes of lactation and the mechanisms of target gene expression and regulation.     

Serum-Derived Neuronal Exosomal miRNAs as Biomarkers of Acute Severe Stress

Stress is the physical and psychological tension felt by an individual while adapting to difficult situations. Stress is known to alter the expression of stress hormones and cause neuroinflammation in the brain. In this study, miRNAs in serum-derived neuronal exosomes (nEVs) were analyzed to determine whether differentially expressed miRNAs could be used as biomarkers of acute stress. Specifically, acute severe stress was induced in Sprague-Dawley rats via electric foot-shock treatment. In this acute severe-stress model, time-dependent changes in the expression levels of stress hormones and neuroinflammation-related markers were analyzed. In addition, nEVs were isolated from the serum of control mice and stressed mice at various time points to determine when brain damage was most prominent; this was found to be 7 days after foot shock. Next-generation sequencing was performed to compare neuronal exosomal miRNA at day 7 with the neuronal exosomal miRNA of the control group. From this analysis, 13 upregulated and 11 downregulated miRNAs were detected. These results show that specific miRNAs are differentially expressed in nEVs from an acute severe-stress animal model. Thus, this study provides novel insights into potential stress-related biomarkers.     

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5′ RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.     

Identification of microRNAs in the Lyme Disease Vector Ixodes scapularis

MicroRNAs (miRNAs) are a class of small non-coding RNAs involved in many biological processes, including the immune pathways that control bacterial, parasitic, and viral infections. Pathogens probably modify host miRNAs to facilitate successful infection, so they might be useful targets for vaccination strategies. There are few data on differentially expressed miRNAs in the black-legged tick Ixodes scapularis after infection with Borrelia burgdorferi, the causative agent of Lyme disease in the United States. Small RNA sequencing and qRT-PCR analysis were used to identify and validate differentially expressed I. scapularis salivary miRNAs. Small RNA-seq yielded 133,465,828 (≥18 nucleotides) and 163,852,135 (≥18 nucleotides) small RNA reads from Borrelia-infected and uninfected salivary glands for downstream analysis using the miRDeep2 algorithm. As such, 254 miRNAs were identified across all datasets, 25 of which were high confidence and 51 low confidence known miRNAs. Further, 23 miRNAs were differentially expressed in uninfected and infected salivary glands: 11 were upregulated and 12 were downregulated upon pathogen infection. Gene ontology and network analysis of target genes of differentially expressed miRNAs predicted roles in metabolic, cellular, development, cellular component biogenesis, and biological regulation processes. Several Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including sphingolipid metabolism; valine, leucine and isoleucine degradation; lipid transport and metabolism; exosome biogenesis and secretion; and phosphate-containing compound metabolic processes, were predicted as targets of differentially expressed miRNAs. A qRT-PCR assay was utilized to validate the differential expression of miRNAs. This study provides new insights into the miRNAs expressed in I. scapularis salivary glands and paves the way for their functional manipulation to prevent or treat B. burgdorferi infection.     

Differential MicroRNA Expression in Porcine Endometrium Related to Spontaneous Embryo Loss during Early Pregnancy

Litter size is an important indicator to measure the production capacity of commercial pigs. Spontaneous embryo loss is an essential factor in determining sow litter size. In early pregnancy, spontaneous embryo loss in porcine is as high as 20–30% during embryo implantation. However, the specific molecular mechanism underlying spontaneous embryo loss at the end of embryo implantation remains unknown. Therefore, we comprehensively used small RNA sequencing technology, bioinformatics analysis, and molecular experiments to determine the microRNA (miRNA) expression profile in the healthy and arresting embryo implantation site of porcine endometrium on day of gestation (DG) 28. A total of 464 miRNAs were identified in arresting endometrium (AE) and healthy endometrium (HE), and 139 differentially expressed miRNAs (DEMs) were screened. We combined the mRNA sequencing dataset from the SRA database to predict the target genes of these miRNAs. A quantitative real-time PCR assay identified the expression levels of miRNAs and mRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed on differentially expressed target genes of DEMs, mainly enriched in epithelial development and amino acids metabolism-related pathways. We performed fluorescence in situ hybridization (FISH) and the dual-luciferase report gene assay to confirm miRNA and predicted target gene binding. miR-205 may inhibit its expression by combining 3′-untranslated regions (3′ UTR) of tubulointerstitial nephritis antigen-like 1 (TINAGL1). The resulting inhibition of angiogenesis in the maternal endometrium ultimately leads to the formation of arresting embryos during the implantation period. This study provides a reference for the effect of miRNA on the successful implantation of pig embryos in early gestation.     

Identification and Dynamic Regulation of microRNAs Involved in Salt Stress Responses in Functional Soybean Nodules by High-Throughput Sequencing

Both symbiosis between legumes and rhizobia and nitrogen fixation in functional nodules are dramatically affected by salt stress. Better understanding of the molecular mechanisms that regulate the salt tolerance of functional nodules is essential for genetic improvement of nitrogen fixation efficiency. microRNAs (miRNAs) have been implicated in stress responses in many plants and in symbiotic nitrogen fixation (SNF) in soybean. However, the dynamic regulation of miRNAs in functioning nodules during salt stress response remains unknown. We performed deep sequencing of miRNAs to understand the miRNA expression profile in normal or salt stressed-soybean mature nodules. We identified 110 known miRNAs belonging to 61 miRNA families and 128 novel miRNAs belonging to 64 miRNA families. Among them, 104 miRNAs were dramatically differentially expressed (>2-fold or detected only in one library) during salt stress. qRT-PCR analysis of eight miRNAs confirmed that these miRNAs were dynamically regulated in response to salt stress in functional soybean nodules. These data significantly increase the number of miRNAs known to be expressed in soybean nodules, and revealed for the first time a dynamic regulation of miRNAs during salt stress in functional nodules. The findings suggest great potential for miRNAs in functional soybean nodules during salt stress.     

miRNA Profiling and Its Role in Multi-Omics Regulatory Networks Connected with Somaclonal Variation in Cucumber (Cucumis sativus L.)

The role of miRNAs in connection with the phenomenon of somaclonal variation, which occurs during plant in vitro culture, remains uncertain. This study aims to investigate the possible role of miRNAs in multi-omics regulatory pathways in cucumber somaclonal lines. For this purpose, we performed sRNA sequencing (sRNA-seq) from cucumber fruit samples identified 8, 10 and 44 miRNAs that are differentially expressed between somaclones (S1, S2, S3 lines) and the reference B10 line of Cucumis sativus. For miRNA identification, we use ShortStack software designed to filter miRNAs from sRNAs according to specific program criteria. The identification of predicted in-silico targets revealed 2,886 mRNAs encoded by 644 genes. The functional annotation of miRNA’s target genes and gene ontology classification revealed their association with metabolic processes, response to stress, multicellular organism development, biosynthetic process and catalytic activity. We checked with bioinformatic analyses for possible interactions at the level of target proteins, differentially expressed genes (DEGs) and genes affected by genomic polymorphisms. We assume that miRNAs can indirectly influence molecular networks and play a role in many different regulatory pathways, leading to somaclonal variation. This regulation is supposed to occur through the process of the target gene cleavage or translation inhibition, which in turn affects the proteome, as we have shown in the example of molecular networks. This is a new approach combining levels from DNA-seq through mRNA-seq, sRNA-seq and in silico PPI in the area of plants’ somaclonal variation.     

An Integrative Study of Aortic mRNA/miRNA Longitudinal Changes in Long-Term LVAD Support

Studying the long-term impact of continuous-flow left ventricular assist device (CF-LVAD) offers an opportunity for a complex understanding of the pathophysiology of vascular changes in aortic tissue in response to a nonphysiological blood flow pattern. Our study aimed to analyze aortic mRNA/miRNA expression changes in response to long-term LVAD support. Paired aortic samples obtained at the time of LVAD implantation and at the time of heart transplantation were examined for mRNA/miRNA profiling. The number of differentially expressed genes (Pcorr < 0.05) shared between samples before and after LVAD support was 277. The whole miRNome profile revealed 69 differentially expressed miRNAs (Pcorr < 0.05). Gene ontology (GO) analysis identified that LVAD predominantly influenced genes involved in the extracellular matrix and collagen fibril organization. Integrated mRNA/miRNA analysis revealed that potential targets of miRNAs dysregulated in explanted samples are mainly involved in GO biological process terms related to dendritic spine organization, neuron projection organization, and cell junction assembly and organization. We found differentially expressed genes participating in vascular tissue engineering as a consequence of LVAD duration. Changes in aortic miRNA levels demonstrated an effect on molecular processes involved in angiogenesis.     

Characterization of Paraquat-Induced miRNA Profiling Response in hNPCs Undergoing Proliferation

Aberration during the development of the central nervous system (CNS) due to environmental factors underlies a variety of adverse developmental outcomes. Paraquat (PQ) is a widely studied neurotoxicant that perturbs the normal structure/function of adult CNS. Yet, the impacts of PQ exposure on the developing CNS remain unclear. miRNAs represent a class of small non-coding RNA molecules involved in the regulation of neural development. Thus in the present study, we analyzed the impacts of PQ on the miRNome of human neural progenitor cells (hNPCs) during proliferation by using the Exiqon miRCURY™ LNA Array. A total of 66 miRNAs were identified as differentially expressed in proliferating hNPCs upon PQ treatment. miRTarBase prediction identified 1465 mRNAs, including several genes (e.g., nestin, sox1, ngn1) previously proved to be associated with the neural proliferation and differentiation, as target genes of PQ-induced differentially expressed miRNAs. The database for annotation, visualization and integrated discovery (DAVID) bioinformatics analysis showed that target genes were enriched in regulation of cell proliferation and differentiation, cell cycle and apoptosis as well as tumor protein 53 (p53), Wnt, Notch and mitogen-activated protein kinases (MAPK) signaling pathways (p < 0.001). These findings were confirmed by real-time RT-PCR. Based on our results we conclude that PQ-induced impacts on the miRNA profiling of hNPCs undergoing proliferation may underlie the developmental neurotoxicity of PQ.     

Comprehensive Transcriptome Analysis of Hair Follicle Morphogenesis Reveals That lncRNA-H19 Promotes Dermal Papilla Cell Proliferation through the Chi-miR-214-3p/β-Catenin Axis in Cashmere Goats

Cashmere is initiated and develops in the fetal stages and the number and density of secondary hair follicles (SHFs) determine cashmere production and quality. Growing evidence indicates that both microRNA (miRNA) and long non-coding RNA (lncRNA) play an indispensable role in hair follicle (HF) growth and development. However, little is known about miRNAs, lncRNAs, and their functions as well as their interactions during cashmere initiation and development. Here, based on lncRNA and miRNA high-throughput sequencing and bioinformatics analysis, we identified 10,485 lncRNAs, 40,639 mRNAs, and 605 miRNAs in cashmere goat skin during HF induction, organogenesis, and cytodifferentiation stages. Among them, 521 lncRNAs, 5976 genes, and 204 miRNAs were differentially expressed (DE). KEGG analysis of DE genes indicated that ECM–receptor interaction and biosynthesis of amino acids were crucial for HF development. Notch, TGF-beta, and Wnt signaling pathways were also identified, which are conventional pathways associated with HF growth and development. Then, the ceRNA regulatory network was constructed, and the impact of lncRNA H19 was investigated in dermal papilla (DP) cells. The MTT, CCK-8, and EdU assays showed that the viability and proliferation of DP cells were promoted by H19, and mechanistic studies suggested that H19 performed its function through the chi-miR-214-3p/β-catenin axis. The present study created a resource for lncRNA, miRNA, and mRNA studies in cashmere morphogenesis. It could contribute to a better understanding of the molecular mechanism of ncRNAs involved in the regulation of HF growth and development.     

miRNA Expression Profiling in Subcutaneous Adipose Tissue of Monozygotic Twins Discordant for HIV Infection: Validation of Differentially Expressed miRNA and Bioinformatic Analysis

Combined AntiRetroviral Treatments (cARTs) used for HIV infection may result in varied metabolic complications, which in some cases, may be related to patient genetic factors, particularly microRNAs. The use of monozygotic twins, differing only for HIV infection, presents a unique and powerful model for the controlled analysis of potential alterations of miRNAs regulation consequent to cART treatment. Profiling of 2578 mature miRNA in the subcutaneous (SC) adipose tissue and plasma of monozygotic twins was investigated by the GeneChip^® miRNA 4.1 array. Real-time PCR and ddPCR experiments were performed in order to validate differentially expressed miRNAs. Target genes of deregulated miRNAs were predicted by the miRDB database (prediction score > 70) and enrichment analysis was carried out with g:Profiler. Processes in SC adipose tissue most greatly affected by miRNA up-regulation included (i) macromolecular metabolic processes, (ii) regulation of neurogenesis, and (iii) protein phosphorylation. Furthermore, KEGG analysis revealed miRNA up-regulation involvement in (i) insulin signaling pathways, (ii) neurotrophin signaling pathways, and (iii) pancreatic cancer. By contrast, miRNA up-regulation in plasma was involved in (i) melanoma, (ii) p53 signaling pathways, and (iii) focal adhesion. Our findings suggest a mechanism that may increase the predisposition of HIV⁺ patients to insulin resistance and cancer.     

Expression Profile of microRNAs during Development of the Hypopharyngeal Gland in Honey Bee,Apis mellifera

The hypopharyngeal gland is an important organ for honey bees to secrete royal jelly, and its secretory activity varies with the age of workers. However, by now, the regulation mechanism of hypopharyngeal gland development is still unclear. Here, the expression profiles of miRNAs in the hypopharyngeal gland of newly emerged workers, nurses, and foragers were investigated via small RNA sequencing. From these three stages, 81 known miRNAs and 135 novel miRNAs have been identified. A total of 85 miRNAs showed expression differences between different development stages, and their target genes were predicted to range from 1 to more than 10. Many of the differentially expressed miRNAs and target genes are related to growth and development or apoptosis. Moreover, dual-luciferase-reporter assays verified that novel-miR-11 directly targets the 3′-untranslated regions of LOC410685 (inactive tyrosine-protein kinase 7) and LOC725318 (uncharacterized protein). These results suggested that miRNAs were widely involved in the developmental regulation of the hypopharyngeal gland in honey bees.