- Research article
- Open Access
- Open Peer Review
Integrative analyses of genes and microRNA expressions in human trisomy 21 placentas
BMC Medical Genomicsvolume 11, Article number: 46 (2018)
The most frequent chromosomal aneuploidy is trisomy 21 (T21) that is caused by an extra copy of chromosome 21. The imbalance of whole genome including genes and microRNAs contributes to the various phenotypes of T21. However, the integrative association between genes and microRNAs in the T21 placenta has yet to be determined.
We analyzed the expressions of genes and microRNAs in the whole genomes of chorionic villi cells from normal and T21 human fetal placentas based on our prior studies. The functional significances and interactions of the genes and microRNAs were predicted using bioinformatics tools.
Among 110 genes and 34 microRNAs showing significantly differential expression between the T21 and normal placentas, the expression levels of 17 genes were negatively correlated with those of eight microRNAs in the T21 group. Of these 17 genes, 10 with decreased expression were targeted by five up-regulated microRNAs, whereas seven genes with increased expression were targeted by three down-regulated microRNAs. These genes were significantly associated with hydrogen peroxide-mediated programmed cell death, cell chemotaxis, and protein self-association. They were also associated with T21 and its accompanying abnormalities. The constructed interactive signaling network showed that seven genes (three increased and four decreased expressions) were essential components of a dynamic signaling complex (P = 7.77e-16).
In this study, we have described the interplay of genes and microRNAs in the T21 placentas and their modulation in biological pathways related to T21 pathogenesis. These results may therefore contribute to further research about the interaction of genes and microRNAs in disease pathogenesis.
Trisomy 21 (T21) is the most frequent chromosome aneuploidy affecting 1 in 700 live births . Individuals with T21 have an increased risk of various congenital abnormalities, including eye, cardiac, gastrointestinal, renal and urinary tract defects . These defects are generally considered to originate from gene dosage imbalance between the trisomic genes on chromosome 21 and the disomic genes on other chromosomes. Therefore, studies of T21 have focused mainly on expression levels of the chromosome 21-derived genome in various tissues from subjects with the condition [3,4,5]. However, the severity and incidence of those phenotypic abnormalities are variable within the T21 population, possibly due to the genetic and epigenetic backgrounds of each individual.
MicroRNAs (miRNAs) are small (18~ 25-nucleotide-long) non-coding endogenous RNAs. They regulate expression of genes at the post-transcriptional level by regulating mRNA stability and translation [6,7,8,9]. More than 1000 miRNAs are expected to participate in regulating over 60% of all the genes [10, 11]. Hence, miRNAs seem to be involved in almost all cellular processes, such as cell apoptosis, differentiation, development, and proliferation [12, 13]. Moreover, changes in their expression levels are reported in various human diseases such as cancer, cardiovascular disease, mental retardation, fetal growth restriction, Alzheimer’s disease, and T21 [14,15,16,17,18,19]. For this reason, great attentions are currently devoted to miRNA research.
In our previous studies, we profiled expression levels of genome-wide miRNAs and genes in placental samples from normal and T21 fetuses using microarray analyses [20, 21]. Our results demonstrated that 34 miRNAs (16 up-regulated and 18 down-regulated) and 110 genes (77 up-regulated and 33 down-regulated) were significantly differentially expressed in the T21 placenta compared with that in normal placentas. Moreover, these miRNAs targeted 76 genes on chromosome 21, suggesting a relationship between genetic and epigenetic changes in the placentas of fetuses with T21. However, the association between genes and miRNA expressions in the whole genome has not yet been determined in the T21 placenta, and the functional significances of these genetic and epigenetic interactions are also unclear. Therefore, an integrative investigation of human genes and miRNAs in the whole genome might be important in understanding the complex genetic-epigenetic mechanisms involved in the pathogenesis of T21 associated abnormalities.
The miRNAs and genes differentially expressed between the placentas of normal and T21 fetuses were found in our previous data [20, 21], and in this study we could identify genes showing a negative correlation with miRNAs, and explored the biological function and molecular pathways of the identified genes using various bioinformatics tools.
The placenta cells were collected by the chorionic villi sampling (CVS) from first-trimester pregnant women. The written informed consents were obtained from participants in compliance with the Declaration of Helsinki. The institutional review board approval was received from the Ethics Committee at Cheil General Hospital (#CGH-IRB-2011-85). The fetal karyotype was analyzed by standard protocols using the Giemsa banding procedure. All trisomy samples used in this study were completely T21, and all normal samples were completely euploid.
Expression profiling of genome-wide whole genes and miRNAs
Expression of genes and miRNAs in the whole genomes was analyzed based on our previous studies [20, 21]. In brief, total RNA was extracted from normal and T21 fetal placentas. An RNA quantity, quality, and integrity number were measured by an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA) and a NanoVue Plus spectrophotometer (GE Healthcare, London, UK). Expression profiles of whole genes were determined using the Affymetrix GeneChip Human Genome U133 Plus 2.0 Array (Affymetrix Inc., Santa Clara, CA, USA) . The miRNAs expressions were profiled using Human miRNA Microarray kit, 8 × 60 K (based on miRBase release 16.0, Agilent Technologies) . Differences in expressions of genes and miRNAs between the T21 and normal groups were considered significant at a P-value of < 0.05. The Benjamini–Hochberg procedure was used to set the false discovery rate (FDR) at 0.05 .
Functional annotation of the candidate genes
All target genes of the 34 miRNAs differentially expressed in T21 were compared with 110 genes differentially expressed in T21 using the VENNY tool (http://bioinfogp.cnb.csic.es/tools/venny_old/index.html). The genes showing a negative correlation with miRNAs in terms of expression were selected as candidates for functional annotation. The web-based gene set analysis toolkit (http://www.webgestalt.org/webgestalt_2013) was used for gene ontology (GO) analysis, Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, and disease-associated analysis. The Search Tool for the Retrieval of Interacting Genes database (STRING; http://version10.string-db.org/) is a database to retrieve and display protein-protein interactions, including both physical and functional interactions. STRING was used to analyze an interaction of candidate genes.
The clinical characteristics were analyzed by the Mann-Whitney U-test and χ2-test. In all tests, a value of P < 0.05 was considered statistically significant (SPSS Inc., Chicago, IL, USA).
Study subjects were constructed with 10 women with euploid fetuses and seven women with T21 fetuses. The CVS from five euploid and three T21 placentas were used for expression profiling of genome-wide whole genes , and those from the other five euploid and four T21 placentas were used for miRNA expression analysis . Table 1 shows the clinical characteristics of the subjects. At the time of CVS, there were no significant differences between the two groups with regard to maternal and fetal characteristics (P > 0.05 for all).
Based on our previous studies, we had analyzed the expression levels of over 47,000 genes and 1349 miRNAs in the whole human genome. Thirty-four miRNAs were differentially expressed (16 up-regulated and 18 down-regulated) in the T21 placentas relative to that in the normal placentas. The number of predicted target genes of the up-regulated and down-regulated miRNAs was 7421 and 6058, respectively (Fig. 1). The 110 genes showing significant differential expression between the T21 and normal placentas were identified, of which 77 genes were up-regulated and 33 genes were down-regulated in the T21 group. The 110 genes and 34 miRNAs were selected for the analysis of genetic-epigenetic association in the T21 placenta. In the analysis of the miRNA-gene associations, 17 genes showed a negative correlation with eight miRNAs in terms of expression (Fig. 1 and Table 2). These 17 genes were selected as candidate genes for the functional annotation.
In the in-silico analysis using the 17 candidate genes, GO annotation and disease association analyses were performed by a statistical hypergeometric test (Table 3). In the “biological process” category of GO annotation, the candidate genes, HGF and MAP3K5, were significantly associated with hydrogen peroxide-mediated programmed cell death (adjP = 0.0008). The F2RL1, HGF, and JAM3 genes were associated with cell chemotaxis (adjP = 0.0435). Protein self-association (adjP = 0.0172) in the molecular function category of GO annotation was significantly associated with genes, AGA and DYRK1A. However, none of the candidate genes were associated with the “cellular component” category of GO annotation. The disease associations of the candidate genes are shown in Table 4. The most statistically significant association with candidate genes, AGA, DYRK1A, SETD4, and TTC3, was found in mental retardation (adjP = 0.0014). Besides this, the candidate genes were also significantly associated with T21, neurobehavioral manifestations, chromosome disorders, osteoarthritis, and fibrosis (adjP < 0.05 for all).
An interaction of the candidate genes was predicted by STRING tool (Fig. 2). The list of the identified candidate genes was used to reveal their functional interactions. Each node represents a protein, and each edge represents an interaction. Thicker lines represent stronger associations. On the basis of 17 genes showing a negative correlation with miRNAs, the part of the dynamic signaling complex in the constructed interaction network included three up-regulated (U2AF1, DYRK1A, and KSR1) and four down-regulated genes (MRPL43, F2RL1, TICAM2, and MAP3K5) (Fig. 2C).
To date, most studies have confirmed a gene dose effect of chromosome 21 in T21, because the main etiology of T21 has been known as an imbalance dosage of genes on chromosome 21 [23, 24]. However, the downstream consequences of T21 are complex. In other words, gene expression imbalances of chromosome 21 affect transcription factors, chromatin remodeling proteins, or related molecules on other chromosomes [25, 26]. These suggest that genomic dosage changes of chromosome 21 in T21 could be relatively subtle or massively disruptive to various genes on other chromosomes. Therefore, as well as primary gene dosage effects in pathogenesis T21, secondary (downstream) effects of disomic genes on other chromosomes are also likely to have a major role in T21. In gene expression process of the human genome, the miRNAs affect as key regulators. They induce mRNA degradation of target genes by perfect binding to target mRNAs or inhibit their translation by imperfect complementary binding to the 3′ untranslated region [27, 28]. Moreover, these miRNAs are involved in the occurrence and development of various diseases [14,15,16,17,18,19]. Therefore, understanding the whole genomic changes that contribute to the various phenotypes of T21 is becoming a major goal in T21 research. A comprehensive investigation of genes and miRNAs in the whole genome may improve our understanding of the genetic-epigenetic interactions of T21.
In this study, we investigated the genes and miRNAs with abnormal expression in placentas from T21 fetuses compared with that in euploid fetuses and found 17 genes that were negatively regulated by miRNAs in the T21 placentas. Among them, seven genes had increased and 10 had decreased expression. Of the seven genes with increased expression, four genes were on chromosome 21 and were target genes of three down-regulated miRNAs in the T21 placenta. The 10 genes with decreased expression in the T21 placenta were located on the various chromosomes and were target genes of five up-regulated miRNAs. Half of the eight miRNAs of interest are on the X-chromosome. Most of the tested samples were from males (Table 1). In male spermatogenesis, epigenetic changes play a crucial role in meiotic sex chromosome inactivation (MSCI) and escape. Escape from MSCI characterizes a set of miRNA genes such as mir-221, mir-374, mir-470 and mir-741 . Up to 86% of the X-linked miRNAs escape MSCI during male spermatogenesis . This is likely to have little impact in association between gene and miRNA expression of tested male samples.
Additionally, by in silico pathway-based exploratory analysis, we found an interaction network of three up-regulated genes (U2AF1, DYRK1A, and KSR1) and four down-regulated genes (MRPL43, F2RL1, TICAM2, and MAP3K5) in the T21 placenta. Our network showed the possibility of the processes involving DYRK1A and MAP3K5, being in the genetic-epigenetic mechanisms related to T21 pathophysiology. The DYRK1A protein is a member of the dual-specificity tyrosine-regulated kinases (DYRKs), and has the ability to phosphorylate serine/threonine and tyrosine residues. Its gene is located in the Down syndrome critical region of chromosome 21. DYRK1A overexpression alters both the phosphorylation of tau and alternative splicing factor, and causes an imbalance of 3R- and 4R-tau in the T21 brain . In neurons, the hyperphosphorylation and accumulation of tau into neurofibrillary tangles (NFTs) was found to characterize some neurodegenerative disorders, known as taupathies, Alzheimer’s disease being among them . Tau has an important impact on the organization of the cytoskeleton in neurons and, in particular, in the regulation of axonal transport. Therefore, it is considered a strong candidate gene for the neuronal degeneration associated with Down syndrome . Interestingly, researchers have identified molecules that can modulate splicing, selectively targeting DYRK1A and cyclin-dependent kinase-like 1 [34, 35], thus opening up new avenues for T21 therapy. Mitogen-activated protein kinase kinase kinase 5 (MAP3K5) acts as an essential molecule of the MAPK signal transduction. It plays an important role in the cascades of cellular responses evoked by changes in the environment, and mediates the signals that determine cell fate, such as differentiation and survival. In particular, by activating MAPKs, MAP3K5, mediates signaling pathways involved in both the differentiation and survival of neuronal cells. MAP3K5-null mice show impairment of long-term recognition memory, in addition to hyperactivity in a novel environment, and superior motor coordination . Therefore, MAP3K5 seems a good candidate for explaining the mechanisms underpinning intellectual disability and epilepsy . Our results showed that miRNAs negatively regulated the expression of their target genes DYRK1A and MAP3K5 in T21, likely through their transcriptional regulatory mechanisms of translational repression or mRNA degradation. Thus, our finding suggests that expression changes of miRNAs that target DYRK1A and MAP3K5 could lead to the changed levels of these two genes in T21, thereby playing a key part in the role of DYRK1A and MAP3K5 in T21 pathogenesis.
In this study, different CVS samples were used for gene expression profiling and for miRNA expression analysis. The amount of material from CVS available for analysis was limited. Because fetal placenta samples at the first-trimester pregnancy were very difficult to obtain, a small amount of chorionic villus were obtained per case. Therefore, as this study was limited by its small sample size, a larger-scale study is needed to clarify the findings.
To our best knowledge, this is the first study to survey whole genes and miRNAs in placentas of T21 fetuses. This study shows that 17 genes and eight miRNAs were differentially expressed between euploid and T21 fetuses and they were negatively regulated in T21. Furthermore, our results propose that many biological pathways that have been implicated in T21 and its complications are possibly regulated by these genes. Therefore, the present work provides a variety of information that may give to a better understanding of genetic-epigenetic modulations in T21.
Chorionic villus sampling
Mégarbané A, Ravel A, Mircher C, Sturtz F, Grattau Y, Rethoré MO, Delabar JM, Mobley WC. The 50th anniversary of the discovery of trisomy 21: the past, present, and future of research and treatment of trisomy 21. Genet Med. 2009;11:611–6.
Korenberg JR, Kawashima H, Pulst SM, Ikeuchi T, Ogasawara N, Yamamoto K, Schonberg SA, West R, Allen L, Magenis E, et al. Molecular definition of a region of chromosome 21 that causes features of the Down syndrome phenotype. Am J Hum Genet. 1990;47:236–46.
Li CM, Guo M, Salas M, Schupf N, Silverman W, Zigman WB, Husain S, Warburton D, Thaker H, Tycko B. Cell type-specific over-expression of chromosome 21 genes in fibroblasts and fetal hearts with trisomy 21. BMC Med Genet. 2006;7:24.
Giannone S, Strippoli P, Vitale L, Casadei R, Canaider S, Lenzi L, D'Addabbo P, Frabetti F, Facchin F, Farina A, Carinci P, Zannotti M. Gene expression profile analysis in human T lymphocytes from patients with Down syndrome. Ann Hum Genet. 2004;68:546–54.
Mao R, Wang X, Spitznagel EL Jr, Frelin LP, Ting JC, Ding H, Kim JW, Ruczinski I, Downey TJ, Pevsner J. Primary and secondary transcriptional effects in the developing human Down syndrome brain and heart. Genome Biol. 2005;6:R107.
Bushati N, Cohen SM. microRNA functions. Annu Rev Cell Dev Biol. 2007;23:175–205.
Kloosterman WP, Plasterk RH. The diverse functions of microRNAs in animal development and disease. Dev Cell. 2006;11:441–50.
MicroRNAs BDP. Genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R. Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev. 2006;20:515–24.
Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19:92–105.
Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–73.
Pillai RS, Bhattacharyya SN, Filipowicz W. Repression of protein synthesis by miRNAs: how many mechanisms? Trends Cell Biol. 2007;17:118–26.
Krützfeldt J, Stoffel M. MicroRNAs: a new class of regulatory genes affecting metabolism. Cell Metab. 2006;4:9–12.
Adams BD, Kasinski AL, Slack FJ. Aberrant regulation and function of microRNAs in cancer. Curr Biol. 2014;24:R762–76.
Condorelli G, Latronico MV, Cavarretta E. microRNAs in cardiovascular diseases: current knowledge and the road ahead. J Am Coll Cardiol. 2014;63:2177–87.
Szulwach KE, Jin P, Alisch RS. Noncoding RNAs in mental retardation. Clin Genet. 2009;75:209–19.
Huang L, Shen Z, Xu Q, Huang X, Chen Q, Li D. Increased levels of microRNA-424 are associated with the pathogenesis of fetal growth restriction. Placenta. 2013;34(7):624.
Provost P. MicroRNAs as a molecular basis for mental retardation, Alzheimer's and prion diseases. Brain Res. 2010;1338:58–66.
Siew WH, Tan KL, Babaei MA, Cheah PS, Ling KH. MicroRNAs and intellectual disability (ID) in Down syndrome, X-linked ID, and fragile X syndrome. Front Cell Neurosci. 2013;7:41.
Lim JH, Kim DJ, Lee DE, Han JY, Chung JH, Ahn HK, Lee SW, Lim DH, Lee YS, Park SY, Ryu HM. Genome-wide microRNA expression profiling in placentas of fetuses with Down syndrome. Placenta. 2015;36:322–8.
Lim JH, Han YJ, Kim HJ, Kwak DW, Park SY, Chun SH, Ryu HM. Genome-wide gene expression analysis in the placenta from fetus with trisomy 21. BMC Genomics. 2017;18:720.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57:289–300.
Korenberg JR, Chen XN, Schipper R, Sun Z, Gonsky R, Gerwehr S, Carpenter N, Daumer C, Dignan P, Disteche C, et al. Down syndrome phenotypes: the consequences of chromosomal imbalance. Proc Natl Acad Sci US A. 1994;91:4997–5001.
Pritchard MA, Kola I. The “gene dosage effect” hypothesis versus the “amplified developmental instability” hypothesis in Down syndrome. J Neural Transm Suppl. 1999;57:293–303.
FitzPatrick DR. Transcriptional consequences of autosomal trisomy: primary gene dosage with complex downstream effects. Trends Genet. 2005;21:249–53.
Lockstone HE, Harris LW, Swatton JE, Wayland MT, Holland AJ, Bahn S. Gene expression profiling in the adult Down syndrome brain. Genomics. 2007;90:647–60.
Filipowicz W, Bhattacharyya SN, Sonenberg N. Mechanisms of posttranscriptional regulation by microRNAs: are the answers in sight? Nat Rev Genet. 2008;9:102–14.
Du T, Zamore PD. microPrimer: the biogenesis and function of microRNA. Development. 2005;132:4645–52.
Berletch JB, Yang F, Disteche CM. Escape from X inactivation in mice and humans. Genome Biol. 2010;11:213.
Song R, Ro S, Michaels JD, Park C, McCarrey JR, Yan W. Many X-linked microRNAs escape meiotic sex chromosome inactivation. Nat Genet. 2009;41:488–93.
Bardoni B, Abekhoukh S, Zongaro S, Melko M. Intellectual disabilities, neuronal posttranscriptional RNA metabolism, and RNA-binding proteins: three actors for a complex scenario. Prog Brain Res. 2012;197:29–51.
Ittner LM, Götz J. Amyloid-β and tau--a toxic pas de deux in Alzheimer's disease. Nat Rev Neurosci. 2011;12:65–72.
Wegiel J, Kaczmarski W, Barua M, Kuchna I, Nowicki K, Wang KC, Wegiel J, Yang SM, Frackowiak J, Mazur-Kolecka B, Silverman WP, Reisberg B, Monteiro I, de Leon M, Wisniewski T, Dalton A, Lai F, Hwang YW, Adayev T, Liu F, Iqbal K, Iqbal IG, Gong CX. Link between DYRK1A overexpression and several-fold enhancement of neurofibrillary degeneration with 3-repeat tau protein Down syndrome. J Neuropathol Exp Neurol. 2011;70:36–50.
Giraud F, Alves G, Debiton E, Nauton L, Théry V, Durieu E, Ferandin Y, Lozach O, Meijer L, Anizon F, Pereira E, Moreau P. Synthesis, protein kinase inhibitory potencies, and in vitro antiproliferative activities of meridianin derivatives. J Med Chem. 2011;54:4474–89.
Rosenthal AS, Tanega C, Shen M, Mott BT, Bougie JM, Nguyen DT, Misteli T, Auld DS, Maloney DJ, Thomas CJ. An inhibitor of the Cdc2-like kinase 4 (Clk4). In: Probe reports from the NIH molecular libraries program [internet]. Bethesda (MD): National Center for biotechnology information (US); 2010–.2010 mar 29 [updated 2011 Mar 3].
Kumakura K, Nomura H, Toyoda T, Hashikawa K, Noguchi T, Takeda K, Ichijo H, Tsunoda M, Funatsu T, Ikegami D, Narita M, Suzuki T, Matsuki N. Hyperactivity in novelenvironment with increased dopamine and impaired noveltypreference in apoptosis signal-regulating kinase 1 (ASK1)-deficient mice. Neurosci Res. 2010;66:313–20.
Nguyen LS, Jolly L, Shoubridge C, Chan WK, Huang L, Laumonnier F, Raynaud M, Hackett A, Field M, Rodriguez J, Srivastava AK, Lee Y, Long R, Addington AM, Rapoport JL, Suren S, Hahn CN, Gamble J, Wilkinson MF, Corbett MA, Gecz J. Transcriptome profiling of UPF3B/NMD-deficient lymphoblastoid cells from patients with various forms of intellectual disability. Mol Psychiatry. 2012;17:1103–15.
The authors thank all staffs to participate in this study; Jin Woo Kim, Da Eun Lee, Bom Yi Lee, Ju Yeon Park, Do Jin Kim, Shin Young Kim, Yeon Woo Lee, Ah Rum Oh, Shin Young Lee, So Min Seo, Kyoung Mee Han, and Hwa Jin Choi.
The publication of this article was funded by grants (HI16C0628) from the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea.
Availability of data and materials
The datasets analyzed for the current study are not publicly available due to a patent application but are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved in accordance to the requirements of the Ethics Committee at Cheil General Hospital (#CGH-IRB-2011-85). Written informed consent was obtained from all participants in compliance with the Declaration of Helsinki.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.