Giant Island Mice Exhibit Widespread Gene Expression Changes in Key Metabolic Organs

Abstract Island populations repeatedly evolve extreme body sizes, but the genomic basis of this pattern remains largely unknown. To understand how organisms on islands evolve gigantism, we compared genome-wide patterns of gene expression in Gough Island mice, the largest wild house mice in the world, and mainland mice from the WSB/EiJ wild-derived inbred strain. We used RNA-seq to quantify differential gene expression in three key metabolic organs: gonadal adipose depot, hypothalamus, and liver. Between 4,000 and 8,800 genes were significantly differentially expressed across the evaluated organs, representing between 20% and 50% of detected transcripts, with 20% or more of differentially expressed transcripts in each organ exhibiting expression fold changes of at least 2×. A minimum of 73 candidate genes for extreme size evolution, including Irs1 and Lrp1, were identified by considering differential expression jointly with other data sets: 1) genomic positions of published quantitative trait loci for body weight and growth rate, 2) whole-genome sequencing of 16 wild-caught Gough Island mice that revealed fixed single-nucleotide differences between the strains, and 3) publicly available tissue-specific regulatory elements. Additionally, patterns of differential expression across three time points in the liver revealed that Arid5b potentially regulates hundreds of genes. Functional enrichment analyses pointed to cell cycling, mitochondrial function, signaling pathways, inflammatory response, and nutrient metabolism as potential causes of weight accumulation in Gough Island mice. Collectively, our results indicate that extensive gene regulatory evolution in metabolic organs accompanied the rapid evolution of gigantism during the short time house mice have inhabited Gough Island.

Download Full-text

Differential gene expression analysis of ankylosing spondylitis shows deregulation of the HLA-DRB, HLA-DQB, ITM2A, and CTLA4 genes

Egyptian Journal of Medical Human Genetics ◽

10.1186/s43042-021-00161-0 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Rowan AlEjielat ◽

Anas Khaleel ◽

Amneh H. Tarkhan

Keyword(s):

Gene Expression ◽

Ankylosing Spondylitis ◽

Differentially Expressed Genes ◽

Gene Network ◽

Disease Diagnosis ◽

Receptor Activation ◽

Functional Enrichment ◽

Differentially Expressed ◽

Inflammatory Disorder ◽

Data Set

Abstract Background Ankylosing spondylitis (AS) is a rare inflammatory disorder affecting the spinal joints. Although we know some of the genetic factors that are associated with the disease, the molecular basis of this illness has not yet been fully elucidated, and the genes involved in AS pathogenesis have not been entirely identified. The current study aimed at constructing a gene network that may serve as an AS gene signature and biomarker, both of which will help in disease diagnosis and the identification of therapeutic targets. Previously published gene expression profiles of 16 AS patients and 16 gender- and age-matched controls that were profiled on the Illumina HumanHT-12 V3.0 Expression BeadChip platform were mined. Patients were Portuguese, 21 to 64 years old, were diagnosed based on the modified New York criteria, and had Bath Ankylosing Spondylitis Disease Activity Index scores > 4 and Bath Ankylosing Spondylitis Functional Index scores > 4. All patients were receiving only NSAIDs and/or sulphasalazine. Functional enrichment and pathway analysis were performed to create an interaction network of differentially expressed genes. Results ITM2A, ICOS, VSIG10L, CD59, TRAC, and CTLA-4 were among the significantly differentially expressed genes in AS, but the most significantly downregulated genes were the HLA-DRB6, HLA-DRB5, HLA-DRB4, HLA-DRB3, HLA-DRB1, HLA-DQB1, ITM2A, and CTLA-4 genes. The genes in this study were mostly associated with the regulation of the immune system processes, parts of cell membrane, and signaling related to T cell receptor and antigen receptor, in addition to some overlaps related to the IL2 STAT signaling, as well as the androgen response. The most significantly over-represented pathways in the data set were associated with the “RUNX1 and FOXP3 which control the development of regulatory T lymphocytes (Tregs)” and the “GABA receptor activation” pathways. Conclusions Comprehensive gene analysis of differentially expressed genes in AS reveals a significant gene network that is involved in a multitude of important immune and inflammatory pathways. These pathways and networks might serve as biomarkers for AS and can potentially help in diagnosing the disease and identifying future targets for treatment.

Download Full-text

High heterogeneity undermines generalization of differential expression results in RNA-Seq analysis

Human Genomics ◽

10.1186/s40246-021-00308-5 ◽

2021 ◽

Vol 15 (1) ◽

Author(s):

Weitong Cui ◽

Huaru Xue ◽

Lei Wei ◽

Jinghua Jin ◽

Xuewen Tian ◽

...

Keyword(s):

Gene Expression ◽

Differential Expression ◽

Small Sample ◽

Differentially Expressed ◽

Cancer Type ◽

Rna Seq ◽

Sample Sizes ◽

Large Sample ◽

Expression Levels ◽

Gene Expression Levels

Abstract Background RNA sequencing (RNA-Seq) has been widely applied in oncology for monitoring transcriptome changes. However, the emerging problem that high variation of gene expression levels caused by tumor heterogeneity may affect the reproducibility of differential expression (DE) results has rarely been studied. Here, we investigated the reproducibility of DE results for any given number of biological replicates between 3 and 24 and explored why a great many differentially expressed genes (DEGs) were not reproducible. Results Our findings demonstrate that poor reproducibility of DE results exists not only for small sample sizes, but also for relatively large sample sizes. Quite a few of the DEGs detected are specific to the samples in use, rather than genuinely differentially expressed under different conditions. Poor reproducibility of DE results is mainly caused by high variation of gene expression levels for the same gene in different samples. Even though biological variation may account for much of the high variation of gene expression levels, the effect of outlier count data also needs to be treated seriously, as outlier data severely interfere with DE analysis. Conclusions High heterogeneity exists not only in tumor tissue samples of each cancer type studied, but also in normal samples. High heterogeneity leads to poor reproducibility of DEGs, undermining generalization of differential expression results. Therefore, it is necessary to use large sample sizes (at least 10 if possible) in RNA-Seq experimental designs to reduce the impact of biological variability and DE results should be interpreted cautiously unless soundly validated.

Download Full-text

Behavioral sensitization induced by methamphetamine causes differential alterations in gene expression and histone acetylation of the prefrontal cortex in rats

BMC Neuroscience ◽

10.1186/s12868-021-00616-5 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Hui Li ◽

Jing-An Chen ◽

Qian-Zhi Ding ◽

Guan-Yi Lu ◽

Ning Wu ◽

...

Keyword(s):

Gene Expression ◽

Prefrontal Cortex ◽

Histone Acetylation ◽

Differentially Expressed Genes ◽

Behavioral Sensitization ◽

Functional Enrichment ◽

Differentially Expressed ◽

Adult Rats ◽

Mrna Microarray ◽

H3 Acetylation

Abstract Background Methamphetamine (METH) is one of the most widely abused illicit substances worldwide; unfortunately, its addiction mechanism remains unclear. Based on accumulating evidence, changes in gene expression and chromatin modifications might be related to the persistent effects of METH on the brain. In the present study, we took advantage of METH-induced behavioral sensitization as an animal model that reflects some aspects of drug addiction and examined the changes in gene expression and histone acetylation in the prefrontal cortex (PFC) of adult rats. Methods We conducted mRNA microarray and chromatin immunoprecipitation (ChIP) coupled to DNA microarray (ChIP-chip) analyses to screen and identify changes in transcript levels and histone acetylation patterns. Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, were performed to analyze the differentially expressed genes. We then further identified alterations in ANP32A (acidic leucine-rich nuclear phosphoprotein-32A) and POU3F2 (POU domain, class 3, transcription factor 2) using qPCR and ChIP-PCR assays. Results In the rat model of METH-induced behavioral sensitization, METH challenge caused 275 differentially expressed genes and a number of hyperacetylated genes (821 genes with H3 acetylation and 10 genes with H4 acetylation). Based on mRNA microarray and GO and KEGG enrichment analyses, 24 genes may be involved in METH-induced behavioral sensitization, and 7 genes were confirmed using qPCR. We further examined the alterations in the levels of the ANP32A and POU3F2 transcripts and histone acetylation at different periods of METH-induced behavioral sensitization. H4 hyperacetylation contributed to the increased levels of ANP32A mRNA and H3/H4 hyperacetylation contributed to the increased levels of POU3F2 mRNA induced by METH challenge-induced behavioral sensitization, but not by acute METH exposure. Conclusions The present results revealed alterations in transcription and histone acetylation in the rat PFC by METH exposure and provided evidence that modifications of histone acetylation contributed to the alterations in gene expression caused by METH-induced behavioral sensitization.

Download Full-text

Integrative transcriptome data mining for identification of core lncRNAs in breast cancer

PeerJ ◽

10.7717/peerj.7821 ◽

2019 ◽

Vol 7 ◽

pp. e7821 ◽

Cited By ~ 7

Author(s):

Xiaoming Zhang ◽

Jing Zhuang ◽

Lijuan Liu ◽

Zhengguo He ◽

Cun Liu ◽

...

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Prognostic Value ◽

Early Stage ◽

Expression Profiles ◽

Diagnostic Value ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Differentially Expressed ◽

Data Set

Background Cumulative evidence suggests that long non-coding RNAs (lncRNAs) play an important role in tumorigenesis. This study aims to identify lncRNAs that can serve as new biomarkers for breast cancer diagnosis or screening. Methods First, the linear fitting method was used to identify differentially expressed genes from the breast cancer RNA expression profiles in The Cancer Genome Atlas (TCGA). Next, the diagnostic value of all differentially expressed lncRNAs was evaluated using a receiver operating characteristic (ROC) curve. Then, the top ten lncRNAs with the highest diagnostic value were selected as core genes for clinical characteristics and prognosis analysis. Furthermore, core lncRNA-mRNA co-expression networks based on weighted gene co-expression network analysis (WGCNA) were constructed, and functional enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The differential expression level and diagnostic value of core lncRNAs were further evaluated by using independent data set from Gene Expression Omnibus (GEO). Finally, the expression status and prognostic value of core lncRNAs in various tumors were analyzed based on Gene Expression Profiling Interactive Analysis (GEPIA). Results Seven core lncRNAs (LINC00478, PGM5-AS1, AL035610.1, MIR143HG, RP11-175K6.1, AC005550.4, and MIR497HG) have good single-factor diagnostic value for breast cancer. AC093850.2 has a prognostic value for breast cancer. AC005550.4 and MIR497HG can better distinguish breast cancer patients in early-stage from the advanced-stage. Low expression of MAGI2-AS3, LINC00478, AL035610.1, MIR143HG, and MIR145 may be associated with lymph node metastasis in breast cancer. Conclusion Our study provides candidate biomarkers for the diagnosis and prognosis of breast cancer, as well as a bioinformatics basis for the further elucidation of the molecular pathological mechanism of breast cancer.

Download Full-text

Diagnostic model of combined ceRNA and DNA methylation related genes in esophageal carcinoma

PeerJ ◽

10.7717/peerj.8831 ◽

2020 ◽

Vol 8 ◽

pp. e8831 ◽

Cited By ~ 1

Author(s):

Xiaojiao Guan ◽

Yao Yao ◽

Guangyao Bao ◽

Yue Wang ◽

Aimeng Zhang ◽

...

Keyword(s):

Gene Expression ◽

Esophageal Cancer ◽

Esophageal Carcinoma ◽

Differentially Expressed Genes ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Diagnostic Model ◽

Link Type ◽

Key Genes

Esophageal cancer is a common malignant tumor in the world, and the aim of this study was to screen key genes related to the development of esophageal cancer using a variety of bioinformatics analysis tools and analyze their biological functions. The data of esophageal squamous cell carcinoma from the Gene Expression Omnibus (GEO) were selected as the research object, processed and analyzed to screen differentially expressed microRNAs (miRNAs) and differential methylation genes. The competing endogenous RNAs (ceRNAs) interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, String, and Cytoscape. Biofunctional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of the screened genes and the survival of the patients were verified. By analyzing GSE59973 and GSE114110, we found three down-regulated and nine up-regulated miRNAs. The gene expression matrix of GSE120356 was calculated by Pearson correlation coefficient, and the 11696 pairs of ceRNA relation were determined. In the ceRNA network, 643 lncRNAs and 147 mRNAs showed methylation difference. Functional enrichment analysis showed that these differentially expressed genes were mainly concentrated in the FoxO signaling pathway and were involved in the corresponding cascade of calcineurin. By analyzing the clinical data in The Cancer Genome Atlas (TCGA) database, it was found that four lncRNAs had an important impact on the survival and prognosis of esophageal carcinoma patients. QRT-PCR was also conducted to identify the expression of the key lncRNAs (RNF217-AS1, HCP5, ZFPM2-AS1 and HCG22) in ESCC samples. The selected key genes can provide theoretical guidance for further research on the molecular mechanism of esophageal carcinoma and the screening of molecular markers.

Download Full-text

Comparative transcriptome and histomorphology analysis of testis tissues from mulard and Pekin ducks

Archives Animal Breeding ◽

10.5194/aab-63-303-2020 ◽

2020 ◽

Vol 63 (2) ◽

pp. 303-313

Author(s):

Li Li ◽

Linli Zhang ◽

Zhenghong Zhang ◽

Nemat O. Keyhani ◽

Qingwu Xin ◽

...

Keyword(s):

Gene Expression ◽

Differential Expression ◽

Candidate Genes ◽

Transcriptome Sequencing ◽

Differentially Expressed ◽

Pekin Duck ◽

Comparative Transcriptome ◽

Qrt Pcr ◽

Pekin Ducks ◽

Testis Tissue

Abstract. Testicular transcriptomes were analyzed to characterize the differentially expressed genes between mulard and Pekin ducks, which will help establish gene expression datasets to assist in further determination of the mechanisms of genetic sterility in mulard ducks. Paraffin sections were made to compare the developmental differences in testis tissue between mulard and Pekin ducks. Comparative transcriptome sequencing of testis tissues was performed, and the expression of candidate genes was verified by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). In mulard ducks, spermatogonia and spermatocytes were arranged in a disordered manner, and no mature sperm were observed in the testis tissue. However, different stages of development of sperm were observed in seminiferous tubules in the testis tissue of Pekin ducks. A total of 43.84 Gb of clean reads were assembled into 193 535 UniGenes. Of these, 2131 transcripts exhibited differential expression (false discover rate <0.001 and fold change ≥2), including 997 upregulated and 1134 downregulated transcripts in mulard ducks as compared to those in Pekin duck testis tissues. Several upregulated genes were related to reproductive functions, including ryanodine receptor 2 (RYR2), calmodulin (CALM), argininosuccinate synthase and delta-1-pyrroline-5-carboxylate synthetase ALDH18A1 (P5CS). Downregulated transcripts included the testis-specific serine/threonine-protein kinase 3, aquaporin-7 (AQP7) and glycerol kinase GlpK (GK). The 10 related transcripts involved in the developmental biological process were identified by GO (Gene Ontology) annotation. The KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways indicated that peroxisome proliferator-activated receptors (PPARs) and calcium signaling pathways were significantly (P<0.001) associated with normal testis physiology. The differential expression of select genes implicated in reproductive processes was verified by qRT-PCR, which was consistent with the expression trend of transcriptome sequencing (RNA-seq). Differentially expressed candidate genes RYR2, CALM, P5CS, AQP7 and GK were identified by transcriptional analysis in mulard and Pekin duck testes. These were important for the normal development of the male duck reproductive system. These data provide a framework for the further exploration of the molecular and genetic mechanisms of sterility in mulard ducks. Highlights. The mulard duck is an intergeneric sterile hybrid offspring resulting from mating between Muscovy and Pekin ducks. The transcriptomes of testis tissue from mulard and Pekin ducks were systematically characterized, and differentially expressed genes were screened, in order to gain insights into potential gonad gene expression mechanisms contributing to genetic sterility in mulard ducks.

Download Full-text

Association of differentially expressed genes and autoantibody type in patients with systemic sclerosis

Rheumatology ◽

10.1093/rheumatology/keaa447 ◽

2020 ◽

Author(s):

Jun Inamo

Keyword(s):

Gene Expression ◽

Differentially Expressed Genes ◽

Cellular Response ◽

Rna Polymerase Iii ◽

Skin Lesions ◽

Functional Enrichment ◽

Differentially Expressed ◽

Study Cohort ◽

Healthy Controls ◽

Specific Autoantibody

Abstract Objectives The aims of this study were to investigate the relationship between the type of autoantibody and gene expression profile in skin lesions from patients with SSc, and to identify specific dysregulated pathways in SSc patients compared with healthy controls. Methods Sixty-one patients with SSc from the Genetics vs Environment in Scleroderma Outcome Study cohort and 36 healthy controls were included in this study. Differentially expressed genes were extracted and functional enrichment and pathway analysis were conducted. Results Compared with healthy controls, lists containing 2, 71, 10, 144 and 78 differentially expressed genes were created for patients without specific autoantibody, ACA, anti-U1 RNP antibody (RNP), anti-RNA polymerase III antibody (RNAP) and anti-topoisomerase I antibody (ATA), respectively. While part of the enriched pathways overlapped, distinct pathways were identified except in those patients lacking specific autoantibody. The distinct enriched pathways included ‘keratinocyte differentiation’ for ACA, ‘nuclear factor κB signalling’ and ‘cellular response to TGF-β stimulus’ for RNAP, ‘interferon α/β signalling’ for RNP, and ‘cellular response to stress’ for ATA. Cell type signature score analysis revealed that macrophages/monocytes, endothelial cells and fibroblasts were associated with ACA, RNAP, ATA and the severity of the SSc skin lesions. Conclusion Pathogenic pathways were identified according to the type of autoantibody by leveraging gene expression data of patients and controls from a multicentre cohort. The current study may promote the search for new therapeutic targets for SSc.

Download Full-text

Influence of fatty acid diets on gene expression in rat mammary epithelial cells

Physiological Genomics ◽

10.1152/physiolgenomics.00007.2009 ◽

2009 ◽

Vol 38 (1) ◽

pp. 80-88 ◽

Cited By ~ 12

Author(s):

M. Medvedovic ◽

R. Gear ◽

J. M. Freudenberg ◽

J. Schneider ◽

R. Bornschein ◽

...

Keyword(s):

Gene Expression ◽

Cell Cycle ◽

Fatty Acids ◽

Fatty Acid ◽

Epithelial Cells ◽

Differential Expression ◽

Mammary Epithelial Cells ◽

Mammary Epithelial ◽

Differentially Expressed ◽

Cell Cycle Genes

Background: This study examines the impact of dietary fatty acids on regulation of gene expression in mammary epithelial cells before and during puberty. Methods: Diets primarily consisted of n-9 monounsaturated fatty acids (olive oil), n-6 polyunsaturated fatty acids (safflower), saturated acids (butter), and the reference AIN-93G diet (soy oil). The dietary regimen mimics the repetitive nature of fatty acid exposure in Western diets. Diet-induced changes in gene expression were examined in laser capture microdissected mammary ductal epithelial cells at day of weaning and end of puberty. PCNA immunohistochemistry analysis compared proliferation rates between diets. Results: Genes differentially expressed between each test diets and the reference diet were significantly enriched by cell cycle genes. Some of these genes were involved in activation of the cell cycle pathway or the G2/M check point pathway. Although there were some differences in the level of differential expression, all diets showed qualitatively the same pattern of differential expression compared to the reference diet. Cluster analysis identified an expanded set of cell cycle as well as immunity and sterol metabolism related clusters of differentially expressed genes. Conclusion: Fatty acid-enriched diets significantly upregulated proliferation above normal physiological levels during puberty. Higher cellular proliferation during puberty caused by enriched fatty acid diets poses a potential increase risk of mammary cancer in later life. The human homologs of 27 of 62 cell cycle rat genes are included in a human breast cancer cluster of 45 cell cycle genes, further emphasizing the importance of our findings in the rat model.

Download Full-text

Differential Expression of MicroRNAs between Eutopic and Ectopic Endometrium in Ovarian Endometriosis

Journal of Biomedicine and Biotechnology ◽

10.1155/2010/369549 ◽

2010 ◽

Vol 2010 ◽

pp. 1-29 ◽

Cited By ~ 79

Author(s):

Nicoletta Filigheddu ◽

Ilaria Gregnanin ◽

Paolo E. Porporato ◽

Daniela Surico ◽

Beatrice Perego ◽

...

Keyword(s):

Gene Expression ◽

Differential Expression ◽

Posttranscriptional Regulation ◽

Noncoding Rna ◽

Critical Role ◽

Regulation Of Gene Expression ◽

Differentially Expressed ◽

Ovarian Endometriosis ◽

Eutopic Endometrium ◽

Ectopic Endometrium

Endometriosis, defined as the presence of endometrial tissue outside the uterus, is a common gynecological disease with poorly understood pathogenesis. MicroRNAs are members of a class of small noncoding RNA molecules that have a critical role in posttranscriptional regulation of gene expression by repression of target mRNAs translation. We assessed differentially expressed microRNAs in ectopic endometrium compared with eutopic endometrium in 3 patients through microarray analysis. We identified 50 microRNAs differentially expressed and the differential expression of five microRNAs was validated by real-time RT-PCR in other 13 patients. We identifiedin silicotheir predicted targets, several of which match the genes that have been identified to be differentially expressed in ectopicversuseutopic endometrium in studies of gene expression. A functional analysis of the predicted targets indicates that several of these are involved in molecular pathways implicated in endometriosis, thus strengthening the hypothesis of the role of microRNAs in this pathology.

Download Full-text

Identification of Potentially Functional Circular RNAs hsa_circ_0070934 and hsa_circ_0004315 as Prognostic Factor of Hepatocellular Carcinoma by Integrated Bioinformatics Analysis

10.21203/rs.3.rs-948198/v1 ◽

2021 ◽

Author(s):

Pejman Morovat ◽

Saman Morovat ◽

Arash M. Ashrafi ◽

Shahram Teimourian

Keyword(s):

Gene Expression ◽

Hepatocellular Carcinoma ◽

Prognostic Factor ◽

Sequence Data ◽

Functional Enrichment ◽

Differentially Expressed ◽

Circular Rnas ◽

Ppi Network ◽

Hub Genes ◽

Cerna Network

Abstract Hepatocellular carcinoma (HCC) is one of the most prevalent cancers worldwide, which has a high mortality rate and poor treatment outcomes with yet unknown molecular basis. It seems that gene expression plays a pivotal role in the pathogenesis of the disease. Circular RNAs (circRNAs) can interact with microRNAs (miRNAs) to regulate gene expression in various malignancies by acting as competitive endogenous RNAs (ceRNAs). However, the potential pathogenesis roles of the ceRNA network among circRNA/miRNA/mRNA in HCC are unclear. In this study, first, the HCC circRNA expression data were obtained from three Gene Expression Omnibus microarray datasets (GSE164803, GSE94508, GSE97332), and the differentially expressed circRNAs (DECs) were identified using R limma package. Also, the liver hepatocellular carcinoma (LIHC) miRNA and mRNA sequence data were retrieved from TCGA, and differentially expressed miRNAs (DEMIs) and mRNAs (DEGs) were determined using the R DESeq2 package. Second, CSCD website was used to uncover the binding sites of miRNAs on DECs. The DECs' potential target miRNAs were revealed by conducting an intersection between predicted miRNAs from CSCD and downregulated DEMIs. Third, some related genes were uncovered by intersecting targeted genes predicted by miRWalk and targetscan online tools with upregulated DEGs. The ceRNA network was then built using the Cytoscape software. The functional enrichment and the overall survival time of these potential targeted genes were analyzed, and a PPI network was constructed in the STRING database. Network visualization was performed by Cytoscape, and ten hub genes were detected using the CytoHubba plugin tool. Four DECs (hsa_circ_0000520, hsa_circ_0008616, hsa_circ_0070934, hsa_circ_0004315) were obtained and six miRNAs (hsa-miR-542-5p, hsa-miR-326, hsa-miR-511-5p, hsa-miR-195-5p, hsa-miR-214-3p, and hsa-miR-424-5p) which are regulated by the above DECs were identified. Then 543 overlapped genes regulated by six miRNAs mentioned above were predicted. Functional enrichment analysis showed that these genes are mostly associated with cancer regulation functions. Ten hub genes (TTK،AURKB, KIF20A، KIF23، CEP55، CDC6، DTL، NCAPG، CENPF، PLK4) have been screened from the PPI network of the 204 survival-related genes. KIF20A, NCAPG, TTK, PLK4, and CDC6 were selected for the highest significant p-values. In the end, a circRNA-miRNA-mRNA regulatory axis was established for five final selected hub genes. This study implies the potential pathogenesis of the obtained network and proposes that the two DECs (has_circ_0070934 and has_circ_0004315) may be important prognostic factor for HCC.

Download Full-text