scholarly journals Integrated Analysis of LncRNA-mRNA Coexpression in the Extracellular Matrix of Developing Deciduous Teeth in Miniature Pigs

2019 ◽  
Vol 2019 ◽  
pp. 1-9
Author(s):  
Yang Li ◽  
Guoqing Li ◽  
Fu Wang ◽  
Xiaoshan Wu ◽  
Zhifang Wu ◽  
...  
Keyword(s):  
Extracellular Matrix ◽  
Noncoding Rna ◽  
Tooth Development ◽  
Expression Profiles ◽  
Differentially Expressed ◽  
Deciduous Teeth ◽  
Human Tooth ◽  
Integrated Analysis ◽  
Miniature Pigs ◽  
Tooth Germs

Miniature pigs, a valuable alternative model for understanding human tooth development, have deciduous teeth from all four tooth families that are replaced once by permanent molars. The extracellular matrix (ECM) supports cells and maintains the integrity of tooth germs during tooth development. However, details on the role of the ECM in tooth development are poorly understood. Here, we performed long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles in the ECM components of deciduous tooth germs by RNA sequencing in miniature pigs. From the early cap to the late bell stages, we identified 4,562 and 3,238 differentially expressed genes (DEGs) from E40 to E50 and E50 to E60, respectively. In addition, a total of 1,464 differentially expressed lncRNAs from E40 to E50 and 969 differentially expressed lncRNAs from E50 to E60 were obtained. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs were enriched significantly for multiple signaling pathways, especially for the ECM pathway. We then outlined the detailed dynamic gene expression profiling of ECM components during deciduous molar development. Comparison of the cap and bell stages revealed that the structure and functions of the ECM dynamically changed. The ECM-related genes, including THBS1, COL4A5, COL4A6, COL1A1, CHAD, TNR, GP1BA, and ITGA3, were significantly changed, and some were shown to enrich during the bell stage development. Finally, we outlined the coexpression of lncRNAs and ECM properties during tooth development. We showed that the interplay of key lncRNAs could change ECM processes and influence the ECM establishment of tooth patterns to accomplish full tooth formation. These results might provide information to elucidate the regulation network of the lncRNA and ECM in tooth development.

scholarly journals An integrated analysis of the competing endogenous RNA network and co-expression network revealed seven hub long non-coding RNAs in osteoarthritis

2020 ◽  
Vol 9 (3) ◽  
pp. 90-98 ◽  
Author(s):  
Haitao Chen ◽  
Liaobin Chen
Keyword(s):  
Extracellular Matrix ◽  
Expression Profiles ◽  
Osteoclast Differentiation ◽  
Signalling Pathway ◽  
Differentially Expressed ◽  
Integrated Analysis ◽  
Knee Cartilage ◽  
Competing Endogenous Rna ◽  
Cerna Network ◽  
Non Coding Rnas

Aims This study aimed to uncover the hub long non-coding RNAs (lncRNAs) differentially expressed in osteoarthritis (OA) cartilage using an integrated analysis of the competing endogenous RNA (ceRNA) network and co-expression network. Methods Expression profiles data of ten OA and ten normal tissues of human knee cartilage were obtained from the Gene Expression Omnibus (GEO) database (GSE114007). The differentially expressed messenger RNAs (DEmRNAs) and lncRNAs (DElncRNAs) were identified using the edgeR package. We integrated human microRNA (miRNA)-lncRNA/mRNA interactions with DElncRNA/DEmRNA expression profiles to construct a ceRNA network. Likewise, lncRNA and mRNA expression profiles were used to build a co-expression network with the WGCNA package. Potential hub lncRNAs were identified based on an integrated analysis of the ceRNA network and co-expression network. StarBase and Multi Experiment Matrix databases were used to verify the lncRNAs. Results We detected 1,212 DEmRNAs and 49 DElncRNAs in OA and normal knee cartilage. A total of 75 dysregulated lncRNA-miRNA interactions and 711 dysregulated miRNA-mRNA interactions were obtained in the ceRNA network, including ten DElncRNAs, 69 miRNAs, and 72 DEmRNAs. Similarly, 1,330 dysregulated lncRNA-mRNA interactions were used to construct the co-expression network, which included ten lncRNAs and 407 mRNAs. We finally identified seven hub lncRNAs, named MIR210HG, HCP5, LINC00313, LINC00654, LINC00839, TBC1D3P1-DHX40P1, and ISM1-AS1. Subsequent enrichment analysis elucidated that these lncRNAs regulated extracellular matrix organization and enriched in osteoclast differentiation, the FoxO signalling pathway, and the tumour necrosis factor (TNF) signalling pathway in the development of OA. Conclusion The integrated analysis of the ceRNA network and co-expression network identified seven hub lncRNAs associated with OA. These lncRNAs may regulate extracellular matrix changes and chondrocyte homeostasis in OA progress. Cite this article: Bone Joint Res. 2020;9(3):90–98.

scholarly journals Identification of Circular RNA-MicroRNA-Messenger RNA Regulatory Network in Atrial Fibrillation by Integrated Analysis

2020 ◽  
Vol 2020 ◽  
pp. 1-13
Author(s):  
Tao Liu ◽  
Guoru Zhang ◽  
Yaling Wang ◽  
Mingyue Rao ◽  
Yang Zhang ◽  
...  
Keyword(s):  
Atrial Fibrillation ◽  
Regulatory Network ◽  
Messenger Rna ◽  
Noncoding Rna ◽  
Expression Profiles ◽  
Arrhythmogenic Right Ventricular Cardiomyopathy ◽  
Circular Rna ◽  
Host Gene ◽  
Differentially Expressed ◽  
Integrated Analysis

Background. Circular RNA (circRNA) is a noncoding RNA that forms a closed-loop structure, and its abnormal expression may cause disease. We aimed to find potential network for circRNA-related competitive endogenous RNA (ceRNA) in atrial fibrillation (AF). Methods. The circRNA, miRNA, and mRNA expression profiles in the heart tissue from AF patients were retrieved from the Gene Expression Omnibus database and analyzed comprehensively. Differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were identified, followed by the establishment of DEcircRNA-DEmiRNA-DEmRNA regulatory network. Functional annotation analysis of host gene of DEcircRNAs and DEmRNAs in ceRNA regulatory network was performed. In vitro experiment and electronic validation were used to validate the expression of DEcircRNAs, DEmiRNAs, and DEmRNAs. Results. A total of 1611 DEcircRNAs, 51 DEmiRNAs, and 1250 DEmRNAs were identified in AF. The DEcircRNA-DEmiRNA-DEmRNA network contained 62 circRNAs, 14 miRNAs, and 728 mRNAs. Among which, two ceRNA regulatory pairs of hsa-circRNA-100053-hsa-miR-455-5p-TRPV1 and hsa-circRNA-005843-hsa-miR-188-5p-SPON1 were identified. In addition, six miRNA-mRNA regulatory pairs including hsa-miR-34c-5p-INMT, hsa-miR-1253-DDIT4L, hsa-miR-508-5p-SMOC2, hsa-miR-943-ACTA1, hsa-miR-338-3p-WIPI1, and hsa-miR-199a-3p-RAP1GAP2 were also obtained. MTOR was a significantly enriched signaling pathway of host gene of DEcircRNAs. In addition, arrhythmogenic right ventricular cardiomyopathy, dilated cardiomyopathy, and hypertrophic cardiomyopathy were remarkably enriched signaling pathways of DEmRNAs in DEcircRNA-DEmiRNA-DEmRNA regulatory network. The expression validation of hsa-circRNA-402565, hsa-miR-34c-5p, hsa-miR-188-5p, SPON1, DDIT4L, SMOC2, and WIPI1 was consistent with the integrated analysis. Conclusion. We speculated that hsa-circRNA-100053-hsa-miR-455-5p-TRPV1 and hsa-circRNA-005843-hsa-miR-188-5p-SPON1 interaction pairs may be involved in AF.

scholarly journals Integrated analysis of lncRNA, miRNA and mRNA profiles reveals potential lncRNA functions during early HIV infection

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Lianwei Ma ◽  
Hui Zhang ◽  
Yue Zhang ◽  
Hailong Li ◽  
Minghui An ◽  
...  
Keyword(s):  
Hiv Infection ◽  
Immune Activation ◽  
Expression Profiles ◽  
Critical Role ◽  
Differentially Expressed ◽  
Integrated Analysis ◽  
Target Mrnas ◽  
Immunodeficiency Virus ◽  
Mirna Sequencing ◽  
Hiv 1

Abstract Background Long noncoding RNAs (lncRNAs) can regulate gene expression in a cis-regulatory fashion or as “microRNA sponges”. However, the expression and functions of lncRNAs during early human immunodeficiency virus (HIV) infection (EHI) remain unclear. Methods 3 HAART-naive EHI patients and 3 healthy controls (HCs) were recruited in this study to perform RNA sequencing and microRNA (miRNA) sequencing. The expression profiles of lncRNAs, mRNAs and miRNAs were obtained, and the potential roles of lncRNAs were analysed based on discovering lncRNA cis-regulatory target mRNAs and constructing lncRNA–miRNA–mRNA competing endogenous RNA (ceRNA) networks. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on 175 lncRNA-associated differentially expressed (DE) mRNAs to investigate the potential functions of DE lncRNAs in ceRNA networks. Results A total of 242 lncRNAs, 1240 mRNAs and 21 mature known miRNAs were determined as differentially expressed genes in HAART-naive EHI patients compared to HCs. Among DE lncRNAs, 44 lncRNAs were predicted to overlap with 41 target mRNAs, and 107 lncRNAs might regulate their nearby DE mRNAs. Two DE lncRNAs might regulate their cis-regulatory target mRNAs BTLA and ZAP70, respectively, which were associated with immune activation. In addition, the ceRNA networks comprised 160 DE lncRNAs, 21 DE miRNAs and 175 DE mRNAs. Seventeen DE lncRNAs were predicted to regulate HIF1A and TCF7L2, which are involved in the process of HIV-1 replication. Twenty DE lncRNAs might share miRNA response elements (MREs) with FOS, FOSB and JUN, which are associated with both immune activation and HIV-1 replication. Conclusions This study revealed that lncRNAs might play a critical role in HIV-1 replication and immune activation during EHI. These novel findings are helpful for understanding of the pathogenesis of HIV infection and provide new insights into antiviral therapy.

Integrated analysis profiles of long non–coding RNAs reveal potential biomarkers across brain regions in post–traumatic stress disorder

2020 ◽  
Author(s):  
Yaoyao Bian ◽  
Lili Yang ◽  
Zhongli Wang ◽  
Wen Li ◽  
Qing Wang ◽  
...  
Keyword(s):  
Post Traumatic Stress Disorder ◽  
Traumatic Stress ◽  
Stress Disorder ◽  
Expression Profiles ◽  
Differentially Expressed ◽  
Integrated Analysis ◽  
Post Traumatic Stress ◽  
Post Traumatic ◽  
Non Coding Rnas ◽  
Potential Biomarkers

Abstract Background Post–traumatic stress disorder (PTSD) is characterized by impaired fear extinction, excessive anxiety and depression. However, underlying mechanisms, especially the function roles of long non–coding RNAs (lncRNAs) involved in PTSD is still unclear. We argued that the lncRNAs, co–expressed mRNAs, as well as the associated pathways, are altered and may thus serve as potential biomarkers and key pathways related to PTSD.Methods The gene expression profiles of GSE68077 was downloaded from the GEO database, and the differentially expressed lncRNAs and mRNAs were identified. Gene ontology (GO) and Kyto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analysis were performed. Subsequently, protein–protein interaction (PPI) network was analyzed, and module analysis of the differentially expressed mRNAs was performed with Cytoscape software. Finally, lncRNAs–mRNAs co–expression network was constructed and core pair lncRNAs involved in PTSD were mapped.Results A total of 45 differentially expressed lncRNAs and 726 differentially expressed mRNAs were obtained. Among of which, 17 lncRNAs and 86 mRNAs were inter–regulated, and most of the lncRNAs–mRNAs co–expression showed positive correlations. The lncRNAs–mRNAs co–expressed network suggested the potentially functional roles of lncRNAs, regulated mRNAs and related pathways in PTSD. By implication of the core pair network, lncRNA–NONMMUT010120.2 synergistically up–regulated Ppargc1a and down–regulated Cir1, Slc38a9, Atp6v0a2. Moreover, lncRNA–NONMMUT023440.2, NONMMUT034155.2, NONMMUT105407.1 and NONMMUT149972.1 were co–expressed with 10 co–expressed mRNAs, which indicated that lncRNAs involved in PTSD might work by regulating the co–expressed mRNAs.

scholarly journals Construction and Comprehensive Analysis of Dysregulated Long Noncoding RNA-Associated Competing Endogenous RNA Network in Moyamoya Disease

2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Xuefeng Gu ◽  
Dongyang Jiang ◽  
Yue Yang ◽  
Peng Zhang ◽  
Guoqing Wan ◽  
...  
Keyword(s):  
Moyamoya Disease ◽  
Cerebral Artery ◽  
Regulatory Network ◽  
Noncoding Rna ◽  
Expression Profiles ◽  
Differentially Expressed ◽  
Competing Endogenous Rna ◽  
Cerna Network ◽  
And Control ◽  
Control Samples

Background. Moyamoya disease (MMD) is a rare cerebrovascular disease characterized by chronic progressive stenosis or occlusion of the bilateral internal carotid artery (ICA), the anterior cerebral artery (ACA), and the middle cerebral artery (MCA). MMD is secondary to the formation of an abnormal vascular network at the base of the skull. However, the etiology and pathogenesis of MMD remain poorly understood. Methods. A competing endogenous RNA (ceRNA) network was constructed by analyzing sample-matched messenger RNA (mRNA), long non-coding RNA (lncRNA), and microRNA (miRNA) expression profiles from MMD patients and control samples. Then, a protein-protein interaction (PPI) network was constructed to identify crucial genes associated with MMD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses were employed with the DAVID database to investigate the underlying functions of differentially expressed mRNAs (DEmRNAs) involved in the ceRNA network. CMap was used to identify potential small drug molecules. Results. A total of 94 miRNAs, 3649 lncRNAs, and 2294 mRNAs were differentially expressed between MMD patients and control samples. A synergistic ceRNA lncRNA-miRNA-mRNA regulatory network was constructed. Core regulatory miRNAs (miR-107 and miR-423-5p) and key mRNAs (STAT5B, FOSL2, CEBPB, and CXCL16) involved in the ceRNA network were identified. GO and KEGG analyses indicated that the DEmRNAs were involved in the regulation of the immune system and inflammation in MMD. Finally, two potential small molecule drugs, CAY-10415 and indirubin, were identified by CMap as candidate drugs for treating MMD. Conclusions. The present study used bioinformatics analysis of candidate RNAs to identify a series of clearly altered miRNAs, lncRNAs, and mRNAs involved in MMD. Furthermore, a ceRNA lncRNA-miRNA-mRNA regulatory network was constructed, which provides insights into the novel molecular pathogenesis of MMD, thus giving promising clues for clinical therapy.

scholarly journals Integrated analysis of the aging brain transcriptome and proteome in tauopathy

2020 ◽  
Vol 15 (1) ◽  
Author(s):  
Carl Grant Mangleburg ◽  
Timothy Wu ◽  
Hari K. Yalamanchili ◽  
Caiwei Guo ◽  
Yi-Chen Hsieh ◽  
...  
Keyword(s):  
Gene Expression ◽  
Human Brain ◽  
Expression Profiles ◽  
Gene Expression Profiles ◽  
Differentially Expressed ◽  
Integrated Analysis ◽  
Mutant Form ◽  
Brain Gene Expression ◽  
Brain Transcriptome ◽  
Age Dependent

Abstract Background Tau neurofibrillary tangle pathology characterizes Alzheimer’s disease and other neurodegenerative tauopathies. Brain gene expression profiles can reveal mechanisms; however, few studies have systematically examined both the transcriptome and proteome or differentiated Tau- versus age-dependent changes. Methods Paired, longitudinal RNA-sequencing and mass-spectrometry were performed in a Drosophila model of tauopathy, based on pan-neuronal expression of human wildtype Tau (TauWT) or a mutant form causing frontotemporal dementia (TauR406W). Tau-induced, differentially expressed transcripts and proteins were examined cross-sectionally or using linear regression and adjusting for age. Hierarchical clustering was performed to highlight network perturbations, and we examined overlaps with human brain gene expression profiles in tauopathy. Results TauWT induced 1514 and 213 differentially expressed transcripts and proteins, respectively. TauR406W had a substantially greater impact, causing changes in 5494 transcripts and 697 proteins. There was a ~ 70% overlap between age- and Tau-induced changes and our analyses reveal pervasive bi-directional interactions. Strikingly, 42% of Tau-induced transcripts were discordant in the proteome, showing opposite direction of change. Tau-responsive gene expression networks strongly implicate innate immune activation. Cross-species analyses pinpoint human brain gene perturbations specifically triggered by Tau pathology and/or aging, and further differentiate between disease amplifying and protective changes. Conclusions Our results comprise a powerful, cross-species functional genomics resource for tauopathy, revealing Tau-mediated disruption of gene expression, including dynamic, age-dependent interactions between the brain transcriptome and proteome.

scholarly journals Integrated Analysis of Long Noncoding RNA and Coding RNA Expression in Esophageal Squamous Cell Carcinoma

2013 ◽  
Vol 2013 ◽  
pp. 1-10 ◽  
Author(s):  
Wei Cao ◽  
Wei Wu ◽  
Fachun Shi ◽  
Xiaobing Chen ◽  
Lihua Wu ◽  
...  
Keyword(s):  
Squamous Cell Carcinoma ◽  
Esophageal Squamous Cell Carcinoma ◽  
Cell Carcinoma ◽  
Squamous Cell ◽  
Noncoding Rna ◽  
Cell Cancer ◽  
Differentially Expressed ◽  
Integrated Analysis ◽  
Aberrant Expression ◽  
A Genome

Tumorigenesis is a complex dynamic biological process that includes multiple steps of genetic and epigenetic alterations, aberrant expression of noncoding RNA, and changes in the expression profiles of coding genes. We call the collection of those perturbations in genome space the “cancer initiatome.” Long noncoding RNAs (lncRNAs) are pervasively transcribed in the genome and they have key regulatory functions in chromatin remodeling and gene expression. Spatiotemporal variation in the expression of lncRNAs has been observed in development and disease states, including cancer. A few dysregulated lncRNAs have been studied in cancers, but the role of lncRNAs in the cancer initiatome remains largely unknown, especially in esophageal squamous cell carcinoma (ESCC). We conducted a genome-wide screen of the expression of lncRNAs and coding RNAs from ESCC and matched adjacent nonneoplastic normal tissues. We identified differentially expressed lncRNAs and coding RNAs in ESCC relative to their matched normal tissue counterparts and validated the result using polymerase chain reaction analysis. Furthermore, we identified differentially expressed lncRNAs that are co-located and co-expressed with differentially expressed coding RNAs in ESCC and the results point to a potential interaction between lncRNAs and neighboring coding genes that affect ether lipid metabolism, and the interaction may contribute to the development of ESCC. These data provide compelling evidence for a potential novel genomic biomarker of esophageal squamous cell cancer.

scholarly journals Integrated analysis of miRNA and mRNA expression profiles in testes of Duroc and Meishan boars

2019 ◽  
Author(s):  
Haisheng Ding ◽  
Min Liu ◽  
Changfan Zhou ◽  
Xiangbin You ◽  
Tao Su ◽  
...  
Keyword(s):  
High Throughput Sequencing ◽  
Target Genes ◽  
Expression Profiles ◽  
Male Gamete ◽  
Mirna Expression Profile ◽  
Mapk Signaling ◽  
Differentially Expressed ◽  
Integrated Analysis ◽  
Testis Development ◽  
Differentially Expressed Mirnas

Abstract Background: MicroRNAs (miRNAs) are small non-coding RNAs playing vital roles in regulating posttranscriptional gene expression. Elucidating the expression regulation of miRNAs underlying pig testis development will contribute to a better understanding of boar fertility and spermatogenesis. Results: In this study, miRNA expression profile was investigated in testes of Duroc and Meishan boars at 20, 75, and 270 days of age by high-throughput sequencing. Forty-five differentially expressed miRNAs were identified from testes of Duroc and Meishan boars before and after puberty. Integrated analysis of miRNA and mRNA profiles predicted many miRNA-mRNA pairs. Gene ontology and biological pathway analyses revealed that predicted target genes of ssc-mir-423-5p, ssc-mir-34c, ssc-mir-107, ssc-165 mir-196b-5p, ssc-mir-92a, ssc-mir-320, ssc-mir-10a-5p, and ssc-mir-181b were involved in sexual reproduction, male gamete generation, and spermatogenesis, and GnRH, Wnt, and MAPK signaling pathway. Four significantly differentially expressed miRNAs and their predicted target genes were validated by quantitative real-time polymerase chain reaction, and phospholipase C beta 1 ( PLCβ1) gene was verified to be a target of ssc-mir-423-5p . Conclusions: This study provides an insight into the functional roles of miRNAs in testis development and spermatogenesis and offers useful resources for understanding differences in sexual function development caused by the change in miRNAs expression between Duroc and Meishan boars.

scholarly journals Circular RNA Expression Profiles and Features in NAFLD Mice: A Study Using RNA-Seq Data

2020 ◽  
Author(s):  
Xinlu Yuan ◽  
Jianjun Diao ◽  
Anqing Du ◽  
Song Wen ◽  
Ligang Zhou ◽  
...  
Keyword(s):  
Expression Profile ◽  
Noncoding Rna ◽  
Expression Profiles ◽  
Circular Rna ◽  
Differentially Expressed ◽  
Rna Seq ◽  
Mice Model ◽  
Specific Expression ◽  
Nonalcoholic Fatty Liver ◽  
Sequencing Method

Abstract Background: Nonalcoholic fatty liver disease (NAFLD) is primarily characterized by the hepatic cholesterol accumulation. Circular RNA (circRNA), one of noncoding RNA, involves in many liver diseases progression. However, no recent studies on circRNA expression profiles in NAFLD have been reported previously.Methods: A NAFLD mouse model was constructed by providing high-fat diet (HFD) for 32 weeks. The circRNAs expression profile in normal mice and NAFLD mice were determined using high-output RNA sequencing method and bioinformatics methods, while the differentially expressed circRNAs were confirmed using Sanger sequencing and qRT-PCR. The circRNA-miRNA network was also predicted. The biological functions of circRNAs were annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG).Results: The results demonstrated the successful construction of NAFLD mice model by immunohistology and serology assay. In total, 93 dysregulated circRNAs were observed, including 57 upregulated circRNAs and 36 downregulated circRNAs, in the NAFLD group. The circRNA-miRNA network revealed the complex interaction between circRNAs and its potential miRNA targets in NAFLD. The characteristic of tissue-specific expression in circRNA was demonstrated. The differentially expressed circRNAs with important biological function were also annotated using GO and KEGG. Both DDAH1 and VAV3 genes were found to be associated with the NAFLD development.Conclusions: Taken together, this study demonstrated the circRNAs expression profile and features in NAFLD, which may provide potential biological markers for the pathogenesis of NAFLD.

scholarly journals Integrated analysis of sex-biased mRNA and miRNA expression profiles in the gonad of the discus fish (Symphysodon aequifasciatus)

2018 ◽  
Author(s):  
yuanshuai Fu ◽  
Zhe Xu ◽  
Zaizhong Chen ◽  
Bin Wen ◽  
Jianzhong Gao
Keyword(s):  
Molecular Mechanisms ◽  
Target Genes ◽  
High Efficiency ◽  
Expression Profiles ◽  
Gonad Development ◽  
Differentially Expressed ◽  
Cdna Libraries ◽  
Integrated Analysis ◽  
Illumina Hiseq ◽  
Differentially Expressed Mirnas

The discus fish (Symphysodon aequifasciatus) is an ornamental fish that is well-known around the world. Phenotype investigation indicated that there are no significant differences in appearance between males and females of the discus fish. To better understand the sexual development mechanisms and obtain a high efficiency sex identification method in the artificial reproduction process of the discus fish, we constructed six cDNA libraries from three adult testes and three adult ovaries, and perform RNA-sequencing for identifying sex-biased candidate genes, microRNA (miRNA), and metabolic pathway using the Illumina Hiseq 4000. A total of 50,082 non-redundant genes (unigenes) were identified, of which 18,570 unigenes were significantly overexpressed in testes, and 11,182 unigenes were significantly overexpressed in ovaries, and 8 differentially expressed unigenes were validated by quantitative Real-Time PCR (qPCR). A total of 551 miRNAs were identified, of which 47 miRNAs were differentially expressed between testes and ovaries, and 7 differentially expressed miRNAs and one non-differential miRNA were validated by qPCR. Twenty-four of these differentially expressed miRNAs and their 15 predicted target genes constituted 41 important miRNA-mRNA interaction pairs, which may be important candidates for sex-related miRNAs and sex-related genes in the discus fish. Some of vital sex-related metabolic pathways were also identified that may play key roles in regulating gonad development of the discus fish. These results can provide important insights to better understand molecular mechanisms for sexual dimorphism in gonads development.

Sign in / Sign up

Export Citation Format

Share Document