Integrated Analysis Of miRNAs Research Articles

A statistical integrative analysis method for microRNA and gene expression

AbstractMicroRNA (miRNA), a type of non‐coding RNA, is known to suppress the expression of target genes and there is an increasing number of studies analyzing various miRNAs in cancer. In an integrated analysis of miRNA and gene expressions, it is often difficult to obtain a significant result with a small amount of data. In this study, we proposed a method that allows integrated analysis even with a small amount of miRNA and gene expression data, and examined its validity. The method is based on a group comparison between target genes having specific functions controlled by a specific miRNA and a background group of randomly selected genes. We conducted a simulation to verify the validity of the proposed method. The simulation results indicated that the method can detect significant differences in expression of target genes controlled by a miRNA showing significant expression change. We then applied the proposed method to real data derived from multiple myeloma cells. As a result, some functions were detected as significantly changed due to a significant change of miRNA.

Integrated analysis of miRNA, transcriptome, and degradome sequencing provides new insights into lipid metabolism in perilla seed

MicroRNAs (miRNA) are small noncoding RNAs that play a crucial as molecular regulators in lipid metabolism in various oil crops. Perilla (Perilla frutescens) is a specific oil crop known for its high alpha-linolenic acid (C18:3n3, ALA) content (>65 %) in their seed oils. In view of the regulatory mechanism of miRNAs in perilla remains unclear, we conducted miRNAs and transcriptome sequencing in two cultivars with distinct lipid compositions. A total of 525 unique miRNAs, including 142 differentially expressed miRNAs was identified in perilla seeds. The 318 miRNAs targeted 7,761 genes. Furthermore, we identified 112 regulated miRNAs and their 610 target genes involved in lipid metabolism. MiR159b and miR167a as the core nodes to regulate the expression of genes in oil biosynthesis (e.g., KAS, FATB, GPAT, FAD, DGK, LPAAT) and key regulatory TFs (e.g., MYB, ARF, DOF, SPL, NAC, TCP, and bHLH). The 1,219 miRNA-mRNA regulation modules were confirmed through degradome sequencing. Notably, pf-miR159b-MYBs and pf-miR167a-ARFs regulation modules were confirmed. They exhibited significantly different expression levels in two cultivars and believed to play important roles in oil biosynthesis in perilla seeds. This provides valuable insights into the functional analysis of miRNA-regulated lipid metabolism in perilla seeds.

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.

Integrated mRNA and miRNA transcriptomic analysis reveals the response of Rapana venosa to the metamorphic inducer (juvenile oysters)

Metamorphosis, as a critical developmental event, controls the population dynamics of most marine invertebrates, especially some carnivorous gastropods that feed on bivalves, whose population dynamics not only affect the maintenance of the ecological balance but also impact the protection of bivalve resources; therefore, the metamorphosis of carnivorous gastropods deserve attention. Here, we investigated the mechanism underlying the response of the carnivorous gastropod Rapana venosa to its metamorphic inducer juvenile oysters through integrated analysis of miRNA and mRNA profiles. According to the results, we speculated that the AMPK signaling pathway may be the critical regulator in the response to juvenile oysters in R. venosa competent larvae. The NF-kB and JAK-STAT signaling pathways that regulated apoptosis were also activated by the metamorphic inducer, which may result in the degeneration of the velum. Additionally, the significant changes in the expression of the SARP-19 precursor gene and protein cibby homolog 1-like gene may indicate that these signaling pathways also regulate growth and development during metamorphosis. This study provides further evidence that juvenile oysters can induce metamorphosis of R. venosa at the transcriptional level, which expands our understanding of the metamorphosis mechanism in carnivorous gastropods.

Integrated analysis of miRNAs and mRNA profiling reveals the potential roles of miRNAs in sheep hair follicle development

BackgroundMerino sheep exhibit high wool production and excellent wool quality. The fleece of Merino sheep is predominantly composed of wool fibers grown from hair follicles (HFs). The HF is a complex biological system involved in a dynamic process governed by gene regulation, and gene expression is regulated by microRNAs (miRNAs). miRNA inhibits posttranscriptional gene expression by specifically binding to target messenger RNA (mRNA) and plays an important role in regulating gene expression, the cell cycle and biological development sequences. The purpose of this study was to examine mRNA and miRNA binding to identify key miRNAs and target genes related to HF development. This will provide new and important insights into fundamental mechanisms that regulate cellular activity and cell fate decisions within and outside of the skin.ResultsWe analyzed miRNA data in skin tissues collected from 18 Merino sheep on four embryonic days (E65, E85, E105 and E135) and two postnatal days (D7 and D30) and identified 87 differentially expressed miRNAs (DE-miRNAs). These six stages were further divided into two longer developmental stages based on heatmap cluster analysis, and the results showed that DE-mRNAs in Stage A were closely related to HF morphogenesis. A coanalysis of Stage A DE-mRNAs and DE-miRNAs revealed that 9 DE-miRNAs and 17 DE-mRNAs presented targeting relationships in Stage A. We found that miR-23b and miR-133 could target and regulate ACVR1B and WNT10A. In dermal fibroblasts, the overexpression of miR-133 significantly reduced the mRNA and protein expression levels of ACVR1B. The overexpression of miR-23b significantly reduced the mRNA and protein expression levels of WNT10A.ConclusionThis study provides a new reference for understanding the molecular basis of HF development and lays a foundation for further improving sheep HF breeding. miRNAs and target genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine-wool sheep.

Integrated analysis of miRNA and mRNA expression profiles in testes of Landrace and Hezuo boars.

Precocious puberty is closely related to testicular development and spermatogenesis, and there is increasing evidence that miRNAs are involved in regulation of testicular development and spermatogenesis. However, little is known about the regulation of microRNAs (miRNAs) during precocious maturation in Hezuo (HZ) boars. In this study, serum Testosterone (T), Estradiol (E2), Follicle-stimulating hormone (FSH) and Luteinizing hormone (LH) levels were detected in HZ and Landrace (LC) boars in the postnatal period at 30, 90, 120, 180, and 240 days, and the testes of HZ and LC boars at 30 and 120 days were used for histological observation. In addition, we performed small RNA-Seq to identify miRNA at sexual immaturity (30-days-old) and maturity (120-days-old) of HZ boar testis (using LC boar as control) to reveal the key miRNA in regulation of precocious puberty. Hormone assay results showed that high levels of T, E2, FSH, and LH may be related to precocious sexual maturity of HZ boars, and that FSH may play an important function before sexual maturity. Histological observation showed that HZ boars developed earlier than LC boars and had reached sexual maturity at 120 days. Small RNA-Seq yielded a total of 359 exist miRNAs, 767 known miRNAs and 322 novel miRNAs in 12 samples; 549, 468, 133, and 247 differentially expressed (DE) miRNAs were identified between Ha vs. Hb, La vs. Lb, Ha vs. La, and Hb vs. Lb (log2 fold change >1 and p < 0.05). Enrichment analysis showed that target genes of these DE miRNAs were enriched in many gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways (such as PI3K-Akt, Hippo and Rap1 signaling pathways) were related to testicular development and spermatogenesis. Further screening, some miRNAs (such as ssc-miR-29b, ssc-miR-199b, ssc-miR-383, ssc-miR-149, ssc-miR-615, and ssc-miR-370) were possibly associated with precocious puberty. These results provide new light on miRNA regulatory mechanisms involved in precocious puberty.

Integrated analysis identified novel miRNAs and mRNA in endometriosis.

Endometriosis is a common gynecological disease that seriously affects women's health and quality of life. However, the pathogenesis of endometriosis remains uncertain. This study aims to find the key microRNAs (miRNAs) and mRNAs and further to elucidate the pathogenesis of endometriosis. Differentially expressed mRNAs (DEmRNAs) and the differentially expressed miRNAs (DEmiRNAs) were obtained by Gene Expression Omnibus (GEO) datasets integration analysis. Functional enrichment analysis of DEmRNAs and DEmRNAs targeted by DEmiRNAs was enforced using GeneCodis3. The DEmiRNA-DEmRNA interaction network was built using Cytoscape. The expression of candidate DEmRNA and DEmiRNA was verified using quantitative real time-polymerase chain reaction (QRT-PCR) and online datasets followed by diagnostic and immune cell infiltration analysis. A total of 835 (327 down-regulated and 508 up-regulated) DEmRNAs and 39 (24 down-regulated and 15 up-regulated) DEmiRNAs were identified between ectopic endometria (EC) group and eutopic endometria (EU) group. DEmRNAs targeted by DEmiRNAs were markedly enriched in cell adhesion molecules, pathways in cancer, leukocyte transendothelial migration, cytokine-cytokine receptor interaction and MAPK signaling pathway. The DEmiRNA-DEmRNA interaction network of up-regulated miRNAs was consisted of 15 miRNAs and 188 corresponding mRNAs. For down-regulated miRNAs, the DEmiRNA-DEmRNA interaction network was consisted of 24 miRNAs and 305 corresponding mRNAs. QRT-PCR validation results of IRF6, PTGER3, NTRK2, hsa-miR-449a and hsa-miR-873-5p were in line with the GEO analysis result. RF6, PTGER3 and NTRK2 had a potential diagnostic value for endometriosis. In addition, the infiltration of macrophages M2 and NK cells activated was the most significantly increased and reduced in ectopic endometrial, respectively. These identified DEmRNAs and DEmiRNAs may be may be associated with the pathogenesis of endometriosis. The integrated analysis of miRNA and mRNA expression profiles may provide a new perspective for understanding the mechanisms of endometriosis and developing new treatments.

Integrated analysis of miRNA and mRNA expression profiles in diabetic mouse kidney treated to Korean Red Ginseng

Integrated analysis of miRNA and mRNA expression profiles in diabetic mouse kidney treated to Korean Red Ginseng

Integrated analysis of miRNA and mRNA expression profiles in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis.

Insect gut microbiota has been reported to participate in regulating host multiple biological processes including metabolism and reproduction. However, the corresponding molecular mechanisms remain largely unknown. Recent studies suggest that microRNAs (miRNAs) are involved in complex interactions between the gut microbiota and the host. Here, we used next-generation sequencing technology to characterize miRNA and mRNA expression profiles and construct the miRNA-gene regulatory network in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis. A total of 3016 differentially expressed genes (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Based on the integrated analysis of miRNA and mRNA sequencing data, 229 negatively correlated miRNA-gene pairs were identified from the miRNA-mRNA network. Gene ontology enrichment analysis indicated that DEMs could target several genes involved in the metabolic process, oxidation-reduction process, oogenesis, and insulin signaling pathway. Finally, real-time quantitative polymerase chain reaction further verified the accuracy of RNA sequencing results. In conclusion, our study provides the profiles of miRNA and mRNA expressions under antibiotics treatment and provides an insight into the roles of miRNAs and their target genes in the interaction between the gut microbiota and its host.

Integrated Analysis of Microarray, Small RNA, and Degradome Datasets Uncovers the Role of MicroRNAs in Temperature-Sensitive Genic Male Sterility in Wheat.

Temperature-sensitive genic male sterile (TGMS) line Beijing Sterility 366 (BS366) has been utilized in hybrid breeding for a long time, but the molecular mechanism underlying male sterility remains unclear. Expression arrays, small RNA, and degradome sequencing were used in this study to explore the potential role of miRNA in the cold-induced male sterility of BS366. Microspore observation showed defective cell plates in dyads and tetrads and shrunken microspores at the vacuolated stage. Differential regulation of Golgi vesicle transport, phragmoplast formation, sporopollenin biosynthesis, pollen exine formation, and lipid metabolism were observed between cold and control conditions. Pollen development was significantly represented in the 352 antagonistic miRNA-target pairs in the integrated analysis of miRNA and mRNA profiles. The specific cleavage of ARF17 and TIR1 by miR160 and miR393 were found in the cold-treated BS366 degradome, respectively. Thus, the cold-mediated miRNAs impaired cell plate formation through repression of Golgi vesicle transport and phragmoplast formation. The repressed expression of ARF17 and TIR1 impaired pollen exine formation. The results of this study will contribute to our understanding of the roles of miRNAs in male sterility in wheat.

Integrative Analysis of miRNAs and Their Targets Involved in Ray Floret Growth in Gerbera hybrida.

MicroRNAs (miRNAs) are involved in regulating many aspects of plant growth and development at the post-transcriptional level. Gerbera (Gerbera hybrida) is an important ornamental crop. However, the role of miRNAs in the growth and development of gerbera is still unclear. In this study, we used high-throughput sequencing to analyze the expression profiles of miRNAs in ray floret during inflorescence opening. A total of 164 miRNAs were obtained, comprising 24 conserved miRNAs and 140 novel miRNAs. Ten conserved and 15 novel miRNAs were differentially expressed during ray floret growth, and 607 differentially expressed target genes of these differentially expressed miRNAs were identified using psRNATarget. We performed a comprehensive analysis of the expression profiles of the miRNAs and their targets. The changes in expression of five miRNAs (ghy-miR156, ghy-miR164, ghy-miRn24, ghy-miRn75 and ghy-miRn133) were inversely correlated with the changes in expression of their eight target genes. The miRNA cleavage sites in candidate target gene mRNAs were determined using 5′-RLM-RACE. Several miRNA-mRNA pairs were predicted to regulate ray floret growth and anthocyanin biosynthesis. In conclusion, the results of small RNA sequencing provide valuable information to reveal the mechanisms of miRNA-mediated ray floret growth and anthocyanin accumulation in gerbera.

MicroRNA expression profile of chicken cecum in different stages during Histomonas meleagridis infection

BackgroundHistomonas meleagridis is an anaerobic, intercellular parasite, which infects gallinaceous birds such as turkeys and chickens. In recent years, the reemergence of Histomoniasis has caused serious economic losses as drugs to treat the disease have been banned. At present, H. meleagridis research focuses on virulence, gene expression analysis, and the innate immunity of the host. However, there are no studies on the differentially expressed miRNAs (DEMs) associated with the host inflammatory and immune responses induced by H. meleagridis. In this research, high-throughput sequencing was used to analyze the expression profile of cecum miRNA at 10 and 15 days post-infection (DPI) in chickens infected with Chinese JSYZ-F strain H. meleagridis.ResultsCompared with the controls, 94 and 127 DEMs were found in cecum of infected chickens at 10 DPI (CE vs CC) and 15 DPI (CEH vs CCH), respectively, of which 60 DEMs were shared at two-time points. Gene Ontology (GO) functional enrichment analysis of the target genes of DEMs indicated that 881 and 1027 GO terms were significantly enriched at 10 and 15 DPI, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG, www.kegg.jp/kegg/kegg1.html) pathway enrichment analysis of the target genes of DEMs demonstrated that 5 and 3 KEGG pathways were significantly enriched at 10 and 15 DPI, respectively. For previous uses, the Kanehisa laboratory have happily provided permission. The integrated analysis of miRNA–gene network revealed that the DEMs played important roles in the host inflammatory and immune responses to H. meleagridis infection by dynamically regulating expression levels of inflammation and immune-related cytokines.ConclusionThis article not only suggested that host miRNA expression was dynamically altered by H. meleagridis and host but also revealed differences in the regulation of T cell involved in host responses to different times H. meleagridis infection.

Integrated analysis of expression profiles with meat quality traits in cattle

MicroRNAs (miRNAs) play a vital role in improving meat quality by binding to messenger RNAs (mRNAs). We performed an integrated analysis of miRNA and mRNA expression profiling between bulls and steers based on the differences in meat quality traits. Fat and fatty acids are the major phenotypic indices of meat quality traits to estimate between-group variance. In the present study, 90 differentially expressed mRNAs (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Eighty-three potential DEG targets and 18 DEMs were used to structure a negative interaction network, and 75 matching target genes were shown in this network. Twenty-six target genes were designated as intersection genes, screened from 18 DEMs, and overlapped with the DEGs. Seventeen of these genes enriched to 19 terms involved in lipid metabolism. Subsequently, 13 DEGs and nine DEMs were validated using quantitative real-time PCR, and seven critical genes were selected to explore the influence of fat and fatty acids through hub genes and predict functional association. A dual-luciferase reporter and Western blot assays confirmed a predicted miRNA target (bta-miR-409a and PLIN5). These findings provide substantial evidence for molecular genetic controls and interaction among genes in cattle.

Small RNA profiling for identification of microRNAs involved in regulation of seed development and lipid biosynthesis in yellowhorn

BackgroundYellowhorn (Xanthoceras sorbifolium), an endemic woody oil-bearing tree, has become economically important and is widely cultivated in northern China for bioactive oil production. However, the regulatory mechanisms of seed development and lipid biosynthesis affecting oil production in yellowhorn are still elusive. MicroRNAs (miRNAs) play crucial roles in diverse aspects of biological and metabolic processes in seeds, especially in seed development and lipid metabolism. It is still unknown how the miRNAs regulate the seed development and lipid biosynthesis in yellowhorn.ResultsHere, based on investigations of differences in the seed growth tendency and embryo oil content between high-oil-content and low-oil-content lines, we constructed small RNA libraries from yellowhorn embryos at four seed development stages of the two lines and then profiled small RNA expression using high-throughput sequencing. A total of 249 known miRNAs from 46 families and 88 novel miRNAs were identified. Furthermore, by pairwise comparisons among the four seed development stages in each line, we found that 64 miRNAs (53 known and 11 novel miRNAs) were differentially expressed in the two lines. Across the two lines, 15, 11, 10, and 7 differentially expressed miRNAs were detected at 40, 54, 68, and 81 days after anthesis, respectively. Bioinformatic analysis was used to predict a total of 2654 target genes for 141 differentially expressed miRNAs (120 known and 21 novel miRNAs). Most of these genes were involved in the fatty acid biosynthetic process, regulation of transcription, nucleus, and response to auxin. Using quantitative real-time PCR and an integrated analysis of miRNA and mRNA expression, miRNA-target regulatory modules that may be involved in yellowhorn seed size, weight, and lipid biosynthesis were identified, such as miR172b-ARF2 (auxin response factor 2), miR7760-p3_1-AGL61 (AGAMOUS-LIKE 61), miR319p_1-FAD2–2 (omega-6 fatty acid desaturase 2–2), miR5647-p3_1-DGAT1 (diacylglycerol acyltransferase 1), and miR7760-p5_1-MED15A (Mediator subunit 15a).ConclusionsThis study provides new insights into the important regulatory roles of miRNAs in the seed development and lipid biosynthesis in yellowhorn. Our results will be valuable for dissecting the post-transcriptional and transcriptional regulation of seed development and lipid biosynthesis, as well as improving yellowhorn in northern China.

Integrated analysis of dysregulated microRNA and mRNA expression in intestinal epithelial cells following ethanol intoxication and burn injury

Gut barrier dysfunction is often implicated in pathology following alcohol intoxication and burn injury. MicroRNAs (miRNAs) are negative regulators of gene expression that play a central role in gut homeostasis, although their role after alcohol and burn injury is poorly understood. We performed an integrated analysis of miRNA and RNA sequencing data to identify a network of interactions within small intestinal epithelial cells (IECs) which could promote gut barrier disruption. Mice were gavaged with ~ 2.9 g/kg ethanol and four hours later given a ~ 12.5% TBSA full thickness scald injury. One day later, IECs were harvested and total RNA extracted for RNA-seq and miRNA-seq. RNA sequencing showed 712 differentially expressed genes (DEGs) (padj < 0.05) in IECs following alcohol and burn injury. Furthermore, miRNA sequencing revealed 17 differentially expressed miRNAs (DEMs) (padj < 0.1). Utilizing the miRNet, miRDB and TargetScan databases, we identified both validated and predicted miRNA gene targets. Integration of small RNA sequencing data with mRNA sequencing results identified correlated changes in miRNA and target expression. Upregulated miRNAs were associated with decreased proliferation (miR-98-3p and miR-381-3p) and cellular adhesion (miR-29a-3p, miR-429-3p and miR3535), while downregulated miRNAs were connected to upregulation of apoptosis (Let-7d-5p and miR-130b-5p) and metabolism (miR-674-3p and miR-185-5p). Overall, these findings suggest that alcohol and burn injury significantly alters the mRNA and miRNA expression profile of IECs and reveals numerous miRNA–mRNA interactions that regulate critical pathways for gut barrier function after alcohol and burn injury.

Integrated Bioinformatics Analysis to Identify Alternative Therapeutic Targets for Alzheimer's Disease: Insights from a Synaptic Machinery Perspective.

Alzheimer's disease (AD), the most common type of dementia, is a serious neurodegenerative disease that has no cure yet, but whose symptoms can be alleviated with available medications. Therefore, early and accurate diagnosis of the disease and elucidation of the molecular mechanisms involved in the progression of pathogenesis are critically important. This study aimed to identify dysregulated miRNAs and their target mRNAs through the integrated analysis of miRNA and mRNA expression profiling in AD patients versus unaffected controls. Expression profiles in postmortem brain samples from AD patients and healthy individuals were extracted from the Gene Expression Omnibus database and were analyzed using bioinformatics approaches to identify gene ontologies, pathways, and networks. Finally, the module analysis of the PPI network and hub gene selection was carried out. A total of five differentially expressed miRNAs were extracted from the miRNA dataset, and 4312 differentially expressed mRNAs were obtained from the mRNA dataset. By comparing the DEGs and the putative targets of the altered miRNAs, 116 (3 upregulated and 113 downregulated) coordinated genes were determined. Also, six hub genes (SNAP25, GRIN2A, GRIN2B, DLG2, ATP2B2, and SCN2A) were identified by constructing a PPI network. The results of the present study provide insight into mechanisms such as synaptic machinery and neuronal communication underlying AD pathogenesis, specifically concerning miRNAs.

Integrated analysis of miRNA and mRNA expression profiles identifies potential regulatory interactions during sexual development of Pacific oyster Crassostrea gigas

One of the most important physiological process in the life cycle of most bivalves is the sexual reproduction required for species propagation. Understanding the molecular regulatory mechanisms of gender determination and gonad development have potential application in bivalve aquaculture, but there is still limited knowledge about this issue. The Pacific oyster Crassostrea gigas is a commercially important bivalve in aquaculture worldwide, as well as a representative specie to study the sexual development in bivalves for its complicated sexuality, including dioecy, sex change and rare hermaphroditism. In this study, an integrated analysis of miRNA and mRNA expression profiles of oyster gonads at different gametogenesis stages was performed to explore their possibly regulatory patterns in sexual development of C. gigas. Based on small RNA sequencing, a total of 122 miRNAs were identified, including 53 conserved and 69 novel miRNAs. Six miRNAs exhibited sex-specific pattern that expressed exclusively from either male or female gonads with increasing abundance throughout the gonadal development. Moreover, 287 putative targets of the six sex-biased miRNAs were predicted by combining previous transcriptome data. GO and KEGG enrichment analysis of target genes revealed that these sex-biased miRNAs might participate in the regulation of gametogenesis. Further, potential miRNA-mRNA interactions based on six sex-biased miRNAs were collected. The homolog of a key participator in testis development, meiotic recombination protein SPO11 in C. gigas (CgSPO11), was predicted as a target of female-biased miRNA lgi-miR-96a. This negative miRNA-mRNA expression relationship was validated in additional samples by quantitative real-time PCR (qRT-PCR) and the regulation of CgSPO11 by lgi-miR-96a was preliminarily confirmed by a dual-luciferase reporter assay. This study provides information regarding the gonad miRNAs of C. gigas and leads to the discovery of miRNA-mRNA interactions in gonad development, which will facilitate further investigations in molecular mechanisms of bivalve sex differentiation.

Integrative analysis of miRNAs and mRNAs revealed regulation of lipid metabolism in dairy cattle.

Lipid metabolism in bovine mammary epithelial cells has been the primary focus of the research of milk fat percentage of dairy cattle. Functional microRNAs can affect lipid metabolism by regulating the expression of candidate genes. The purpose of the study was to screen and identify differentially expressed miRNAs, candidate genes, and co-regulatory pathways related to the metabolism of milk fat. To achieve this aim, we used miRNA and transcriptome data from the mammary epithelial cells of dairy cattle with high (H, 4.85%) and low milk fat percentages (L, 3.41%) during mid-lactation. One hundred ninety differentially expressed genes and 33 differentially expressed miRNAs were significantly enriched in related regulatory networks, of which 27 candidate genes regulated by 18 differentially expressed miRNAs significantly enriched in pathways related to lipid metabolism (p < 0.05). Target relationships between PDE4D and bta-miR-148a, PEG10 and bta-miR-877, SOD3 and bta-miR-2382-5p, and ADAMTS1 and bta-miR-2425-5p were verified using luciferase reporter assays and quantitative RT-PCR. The detection of triglyceride production in BMECs showed that bta-miR-21-3p and bta-miR-148a promote triglyceride synthesis, whereas bta-miR-124a, bta-miR-877, bta-miR-2382-5p, and bta-miR-2425-5p inhibit triglyceride synthesis. The conjoint analysis could identify functional miRNAs and regulatory candidate genes involved in lipid metabolism within the co-expression networks of the dairy cattle mammary system, which contributes to the understanding of potential regulatory mechanisms of genetic element and gene signaling networks involved in milk fat metabolism.

Integrated analysis of miRNA and mRNA transcriptomic reveals antler growth regulatory network.

The growth of antler is driven by endochondral ossification in the growth center of the apical region. Antler grows faster than cancer tissues, but it can be stably regulated and regenerated periodically. To elucidate the molecular mechanisms of how antler grows rapidly without carcinogenesis, in this study, we used RNA-seq technology to evaluate the changes of miRNA and mRNA profiles in antler at four different developmental stages, including 15, 60, 90, and 110days. We identified a total of 55004 unigenes and 246 miRNAs of which, 10182, 13258, 10740 differentially expressed (DE) unigenes and 35, 53, 27 DE miRNAs were identified in 60-day vs. 15-day, 90-day vs. 60-day, and 110-day vs. 90-day. GO and KEGG pathway analysis indicated that DE unigenes and DE miRNA were mainly associated with chondrogenesis, osteogenesis and inhibition of oncogenesis, that were closely related to antler growth. The interaction networks of mRNA-mRNA and miRNA-mRNA related to chondrogenesis, osteogenesis and inhibition of oncogenesis of antler were constructed. The results indicated that mRNAs (COL2A1, SOX9, WWP2, FGFR1, SPARC, LOX, etc.) and miRNAs (miR-145, miR-199a-3p, miR-140, miR-199a-5p, etc.) might have key roles in chondrogenesis and osteogenesis of antler. As well as mRNA (TP53, Tpm3 and ATP1A1, etc.) and miRNA (miR-106a, miR-145, miR-1260b and miR-2898, etc.) might play important roles in inhibiting the carcinogenesis of antler. In summary, we constructed the mRNA-mRNA and miRNA-mRNA regulatory networks related to chondrogenesis, osteogenesis and inhibition of oncogenesis of antler, and identified key candidate mRNAs and miRNAs among them. Further developments and validations may provide a reference for in-depth analysis of the molecular mechanism of antler growth without carcinogenesis.

Integrated miRNA and mRNA expression profiling reveals dysregulated miRNA-mRNA regulatory networks in eosinophilic and non-eosinophilic chronic rhinosinusitis with nasal polyps.

The precise mechanisms underlying pathogenesis of different subtypes of chronic rhinosinusitis with nasal polyps (CRSwNP) are still unclear. Emerging evidence indicates that microRNAs may play a role in the pathogenesis of CRSwNP. This study aimed to identify the dysregulated microRNA-messenger RNA (miRNA-mRNA) regulatory networks in eosinophilic (E) and non-eosinophilic (non-E) CRSwNP. Whole-transcriptome sequencing was performed on nasal tissues of patients with ECRSwNP and non-ECRSwNP, and control subjects. An integrated analysis of miRNA and mRNA expression was conducted to identify key mRNAs and miRNAs involved in the pathogenesis of ECRSwNP and non-ECRSwNP. The miRNAs of interest and their target genes were validated using quantitative real-time polymerase chain reaction (PCR). A group of differentially expressed mRNAs (DE-mRNAs) and miRNAs (DE-miRs) were identified in ECRSwNP patients vs control subjects, non-ECRSwNP patients vs control subjects, and non-ECRSwNP vs ECRSwNP patients, respectively. Pathway enrichment analysis showed distinct immune and inflammatory functions associated with DE-mRNAs and target genes of DE-miRs in ECRSwNP vs control and non-ECRSwNP vs control groups. The miRNA-mRNA regulatory networks constructed with Cytoscape highlighted the roles of miR-154, miR-221, and miR-223 family miRNAs relating to both ECRSwNP and non-ECRSwNP, and the roles of the let-7 and miR-34/449 families in the development of non-ECRSwNP. Assessment using real-time PCR for the expression of miRNAs and target genes demonstrated highly consistent data with the RNA sequencing data. ECRSwNP and non-ECRSwNP patients express distinct miRNA-mRNA regulatory networks compared with control subjects, thus providing potential targets for future development of novel therapeutic approaches for the management of CRSwNP.