Abstract
Down syndrome (DS), caused by trisomy of chromosome 21 (HSA21), results in a broad range of phenotypes. However, the determinants contributing to the complex and variable phenotypic expression of DS are still not fully known. Changes in microRNAs (miRNAs), short non-coding RNA molecules that regulate gene expression post-transcriptionally, have been associated with some DS phenotypes. Here, we investigated the genome-wide mature miRNA expression profile in peripheral blood mononuclear cells (PBMCs) of children with DS and controls and identified biological processes and pathways relevant to the DS pathogenesis. The expression of 754 mature miRNAs was profiled in PBMCs from six children with DS and six controls by RT-qPCR using TaqMan® Array Human MicroRNA Cards. Functions and signaling pathways analyses were performed using DIANA-miRPath v.3 and DIANA-microT-CDS software. Children with DS presented six differentially expressed miRNAs (DEmiRs): four overexpressed (miR-378a-3p, miR-130b-5p, miR-942-5p, and miR-424-3p) and two downregulated (miR-452-5p and miR-668-3p). HSA21-derived miRNAs investigated were not found to be differentially expressed between the groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed potential target genes involved in biological processes and pathways pertinent to immune response, e.g., toll-like receptors (TLRs) signaling, Hippo, and transforming growth factor β (TGF-β) signaling pathways. These results suggest that altered miRNA expression could be contributing to the well-known immunological dysfunction observed in individuals with DS.
Similar content being viewed by others
Introduction
Down syndrome (DS) or trisomy 21 (T21) is the most common live-birth human chromosomal disorder [1]. Individuals with DS present a broad range of phenotypes including dysmorphic features [2], intellectual disability [3], congenital heart diseases [4], neurological abnormalities [5], immunodeficiency [6], and increased risk for leukemia [7], among other comorbidities.
Different hypotheses have been proposed to explain how the presence of a third copy of human chromosome 21 (HSA21) results in the DS phenotype. According to the “gene dosage effect” hypothesis, the overexpression of HSA21 genes and the downstream consequences would be directly responsible for DS characteristics [8, 9]. Besides, studies have shown the occurrence of secondary effects as a consequence of the increased gene activity on HSA21 with dysregulation of several disomic genes located throughout the genome that would affect specific cellular processes [10, 11], suggesting the existence of complex molecular mechanisms regulating RNA and protein expression. However, the determinants contributing to the complex and variable phenotypic expression of DS are still not fully known. It is likely that the phenotypic manifestations of DS result from several sources including additional content, functional elements and variability of HSA21, in addition to variability of other chromosomes throughout the genome, the chromatin structure and epigenetic modifications [12].
In light of our current understanding of the importance of genome-wide abnormalities to the phenotype of DS, we profiled the genome-wide mature microRNA (miRNA) expression in peripheral blood mononuclear cells (PBMCs) of children with DS and controls. To assess the potential biological effects of the differentially expressed miRNAs (DEmiRs) observed in our study, the potential targets of these miRNAs were predicted. In addition, following the findings of a previous study [13] by our team that compared the expression of immune-related genes between the same case and control individuals studied here, we investigated the possible interactions between those immune-related differentially expressed genes (DEGs) and the DEmiRs observed in the present study.
Materials and methods
The study protocol was approved by the Research Ethics Committee of São José do Rio Preto Medical School (CEP-FAMERP, protocol 3340/2004 and CAAE number: 88278218.4.0000.5415), in São Paulo State, Brazil. After written informed consent was obtained from parents/legal guardians, six children with DS (three males and three females, median age 4.05; range 2.30–6.10 years) and six children without DS (one male and five females, median age 3.75; range 3.60–5.40 years) were included in the study. Data on age and sex of each study participant are available in Online Resource 1.
The inclusion criteria for case and control groups were ages between 2 and 6 years, no personal history of chronic infection (such as bronchitis, asthma, or recurrent pneumonia), or actual clinical manifestations suggestive of acute infection (including the flu, cough, fever, and/or antibiotic treatment within 10 days preceding blood draw). Specific inclusion criteria were the presence of karyotypically confirmed free trisomy 21 for the case group and the absence of diseases associated with clinical manifestation of DS, such as autoimmune diseases, for the control group.
Children with DS enrolled were attending routine follow-up at the Genetics and Down syndrome Outpatient Pediatric Clinic at Hospital de Base de São José do Rio Preto (HB), in São Paulo State, Brazil. The control group consisted of children attending follow-up at the Outpatient Pediatric Clinic at HB and children attending the daycare center at the Institute of Biosciences, Letters and Exact Sciences (IBILCE), São José do Rio Preto campus—São Paulo State University (UNESP). The personal data of all participants were kept confidential.
Up to 5 mL of whole blood samples were collected by venous puncture into EDTA-containing tubes. PBMCs were isolated using Ficoll® Paque Plus (Sigma-Aldrich) density gradient centrifugation. Total RNA samples were obtained using TRIzol® Reagent (Invitrogen) following the manufacturer’s instructions. RNA quality of the twelve samples used for genome-wide mature miRNA expression analysis was assessed by the Agilent 2100 BioAnalyzer using the Eukaryote Total RNA Pico assay (Agilent Technologies) and all samples presented preserved 18S and 28S ribosomal peaks and RNA integrity number (RIN) equal to or greater than 5 [14]. RNA concentration and purity of all samples were measured by spectrophotometry using NanoDrop ND-1000 (NanoDrop Technologies).
Mature miRNA expression profiling
Mature miRNAs expression profiling was performed by Real-Time qPCR using TaqMan® Array Human MicroRNA A + B Cards Set v3.0 (Applied Biosystems), which enables accurate quantitation of 754 human miRNAs along with three reference genes (MammU6 / U6 snRNA, RNU44, and RNU48) and one negative control (ath-miR159a), according to the manufacturer’s protocol (Online Resource 2).
Briefly, total RNA (50 ng) was reverse-transcribed using the TaqMan® microRNA reverse transcription kit with Megaplex™ RT Primers Human Pools A and B (Applied Biosystems). The cDNA samples were pre-amplified using Megaplex™ PreAmp Primers, Human Pools A and B (Applied Biosystems), on a Veriti Thermal Cycler (Applied Biosystems), following the manufacturer’s instructions.
Array reactions (A and B cards) were performed in duplicate on 7900HT Fast Real-Time PCR System (Applied Biosystems) according to the manufacturer’s recommended conditions. Raw cycle quantification (Cq) data were generated by SDS RQ Manager v.1.2.1 software (Applied Biosystems) and adjusted for basal fluorescence signal and the threshold for each miRNA.
Data were filtered and Cq values above 32 were considered as undetermined, according to the TaqMan® arrays manufacturer’s instructions (Applied Biosystems). Normal distribution of data was obtained from the mean value of case and control groups separately, and Cq values outside the 75% of the central mass of the data were considered as outliers and removed from further analysis. Subsequently, when the two groups were considered together, Cq values outside the 90% of the central mass were considered outliers and excluded. Only data from genes expressed in at least four of the six biological samples in each group met our criteria to be statistically analyzed.
All expression analyses were performed in R statistical environment [15] using packages freely available under the Bioconductor project [16]. The package HTqPCR [17] was used to preprocess the data and as an interface to other packages using the qPCR data as input. miRNA expression data were normalized by geometric averaging [18] of the reference genes mmU6/U6 snRNA, RNU44, and RNU48, and relative expression (fold change, FC) was calculated by the ΔΔCq method [19] using the control group as calibrator. The identification of DEmiRs was addressed using linear models for microarray data (LIMMA) using ΔCq data as input [20]. P values were adjusted for multiple testing issues using Benjamini and Hochberg’s method [21] to control the false discovery rate (FDR) and reduce the risk of false positive. Adjusted p values less than 0.05 were considered statistically significant. Graphical data were presented as FC calculated from the mean of ΔCq of each group.
In silico functional analysis of DEmiRs
Functional characterization analysis of DEmiRs was performed using the miRNA pathway analysis web-server DIANA-miRPath v3.0 [22] (microT threshold 0.8; P value < 0.05; FDR correction; Conservative Stats), an online software suite that uses predicted miRNA targets provided by the DIANA-microT-CDS algorithm and identifies Gene Ontology (GO) biological process and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.
The possible interactions between 17 DEGs observed in a previous study by our group [13] performed on the same individuals and the DEmiRs observed in the present study were investigated using DIANA-microT web server v5.0 [23, 24] (score threshold 0.5). Considering that miRNAs are classically known for the negative regulation of protein translation when there is a direct interaction with the target gene, only the target genes that displayed an inverse expression correlation with the miRNA were selected.
Results
Mature miRNA differential expression profiling
Of the 754 human mature miRNAs investigated using TaqMan® Arrays, 374 were detected in our experimental conditions and considered for further analysis. Initially, we observed 50 miRNAs differentially expressed between PBMCs of individuals with DS and controls (12 upregulated and 38 downregulated) (Online Resource 3). After FDR adjustment, four upregulated and two downregulated miRNAs maintained statistical significance: miR-378a-3p (FC mean = 2.53; P = 0.02), miR-130b-5p (FC mean = 2.71; P = 0.03), miR-942-5p (FC mean = 2.14; P = 0.03), miR-424-3p (FC mean = 2.36; P = 0.03), miR-452-5p (FC mean = 0.15; P = 0.04) and miR-668-3p (FC mean = 0.14; P = 0.01) (Fig. 1).
Eight mature forms of HSA21-derived miRNAs were investigated in the present study: miR-99a-5p and -3p, miR-155-5p and -3p, miR-125b-5p and miR-125b-2-3p, let-7c-5p and miR-802. However, only the four mature forms miR-99a-5p, miR-155-5p, miR-125b-5p, and let-7c-5p were detected in our experimental conditions. miR-155-5p was significantly overexpressed in DS in the initial analysis (Online Resource 3); however, after FDR correction, the statistical significance was lost. The other four HSA21-miRNAs did not present differential expression between DS and control groups (Fig. 2).
Target genes and in silico functional analysis of DEmiRs
Table 1 presents the number of potential target genes of the six DEmiRs observed in our group of children with DS. A total of 2,401 non-redundant targets were predicted (Online Resource 4). Predicted mRNA targets are widely distributed throughout the genome and 17 of them are located on HSA21 (Table 2).
The in silico functional analysis including the 2401 non-redundant predicted target genes showed enrichment of genes involved in biological processes relevant to immune response (e.g., toll-like receptors—TLRs signaling pathways), coagulation, mitotic cell cycle, DNA transcription, gene expression, mRNA processing, among others (P < 0.05) (Fig. 3; Online Resource 5). KEGG analysis also showed enrichment of target genes participating in immune response pathways such a Hippo and TGF-β signaling pathways, among other pathways (P < 0.05) (Fig. 4; Online Resource 6).
The analysis of interactions between the six DEmiRs observed in the present study and the 17 immune-related DEGs previously identified [13] in the same DS and control individuals showed possible interaction between five DEmiRs and 10 DEGs (Table 3).
Discussion
Currently, miRBase [25, 26] version 22.1 (accessed on 08/02/2021) reports 1917 entries for stem-loop miRNA on the human genome resulting in 2654 mature miRNAs. On HSA21, there are 30 precursors reported, resulting in 38 mature miRNAs. Of these precursors, only five HSA21 sequences are classified as “high confidence miRNAs” (i.e. present high levels of certainty that they do exist) or present relevant deep sequencing data (i.e. mean of counts reads per million) (Online Resource 7). The miRNAs miR-155, miR-802, miR-125b-2, let-7c, and miR-99a, located on HSA21, have been the most extensively studied and implicated in the DS pathogenesis [27,28,29,30,31,32]. These miRNAs present two mature forms (-5p and -3p), except miR-802 [25]. In the present study, we investigated the expression pattern of all -5p and -3p mature forms of miR-155, miR-125b-5p, and miR-99a, besides miR-802; however, only -5p mature forms were detected in our samples.
The four HSA21 mature miRNAs detected in the present study did not present significant differential expression patterns between children with DS and controls after FDR adjustment. Although the overexpression of HSA21-located miRNAs has been associated with DS phenotypes [27], including low blood pressure [28], dementia [29], and altered immune processes [30, 31], there are few studies concerning miRNA expression in blood cells from individuals with DS. Xu et al. [33] identified overexpression of miR-99a, let-7c, miR-125b-2, and miR-155 in lymphocytes from children with DS by high-throughput sequencing technology and validated their data using quantitative RT-PCR. Another similar study by the same group [34] showed the same miRNAs downregulated in fetal cord blood mononuclear cells from fetuses with DS. On the other hand, comparable expression of the miRNAs miR-99a, let-7c, and miR-125b-2 were observed between cultured PBMCs from children with DS and age-matched healthy donors while overexpression of the miR-155 was observed in DS cells compared with control ones, investigated by quantitative RT-PCR [31]. Considering the different study designs, in which different methodologies for quantifying miRNAs and analyzing results were used, additional studies are necessary to better elucidate the expression pattern of HSA21 miRNAs and its consequence for the DS phenotype.
The six DEmiRs identified in the present study between children with DS and controls are not located on HSA21; four of them were observed to be upregulated (miR-378a-3p, miR-130b-5p, miR-942-5p, and miR-424-3p) and two downregulated (miR-452-5p and miR-668-3p). Although some DEmiRs identified in the present study have been previously found to be dysregulated in other DS tissues [35,36,37], to our knowledge, this is the first study to report altered expression of these six DEmiRs in PBMCs from individuals with DS. Altered expression of miRNAs not located on HSA21 with possible involvement in DS phenotypes has been previously reported [33,34,35,36,37,38,39,40], supporting the existing hypothesis of secondary transcriptional changes as a consequence of the trisomy 21 [41].
According to the literature, sex can influence gene expression levels [42, 43], as well as the expression of miRNAs [44]. Considering that in the present study the sex was not matched between individuals with DS and controls, we investigated if there was difference in the expression of the six DEmiRs and the four HSA21 mature miRNAs between male and female participants regardless the study group (case or control). miR-452-5p and miR-668 were observed to be overexpressed in females (data not shown), indicating a possible sex bias in miRNA expression. Therefore, we compared the expression of these two miRNAs between case and control groups considering only female participants (data not shown) and obtained the same results from the total sample (both downregulated in the DS group), indicating that there is no sex bias in our results. Data on sex and age for each sample include in the study are available in the Online Resource 1.
Enrichment analyses of GO biological processes and KEGG pathways including the 2401 non-redundant predicted targets of the six DEmiRs point to the contribution of these genes to pathways related to the immune system, including the TLRs, TGF-β, and Hippo signaling pathways. TLRs are key elements of innate immunity as they participate in the recognition of pathogens, antigen presentation, apoptosis, and production of interferons by the virus-infected cells [45]. TLR signaling is dysregulated in children with T21 and this impairment is believed to contribute to the chronic inflammation and sepsis observed in DS [46]. The TGF-β signaling pathway is a pivotal regulator of immune responses, participating in the proliferation, survival, activation, differentiation, and repertoire diversity of B- and T- cells under normal conditions and immune challenges [47]. Moreover, TGF-β regulates the development and functionality of innate cells (dendritic cells, macrophages, granulocytes, and natural killer cells) and the peripheral tolerance against self-antigens [47]. In DS, increased TGF-β levels have been previously associated with transient abnormal myelopoiesis and complications of this condition [48] and pointed as a candidate biomarker for prenatal diagnosis [49].
Studies have revealed extensive roles for the Hippo signaling pathway in adaptive and innate immune systems regulation. Components of the Hippo pathway (MST1/2) are involved in T-cell development and differentiation, B cell homeostasis, and they play key roles in antigen-presenting cells, such as macrophages and dendritic cells [50]. Also, there is evidence of crosstalk between Hippo and TLR signaling pathways [51], along with other signaling networks involved in immune regulation [52]. To our knowledge, there are no studies concerning the role of the Hippo pathway in the DS phenotype and our study is the first to point out a possible association between Hippo pathway dysregulation and DS.
Several of the predicted HSA21-located target genes of our DEmiRs play a role in the immune system. BACH1 is involved in the regulation of the expression of macrophage-associated genes, B-cell development, and antigen presentation [53]; BRWD1 cooperates with networks to drive B-cell development [54]; ERG participates in the control of B lymphopoiesis [55]; GABPA is a regulator of cytokine secretion [56]; RUNX1 regulates TLR signaling pathways and inflammatory cytokine production [57]. Besides the fact that they are present in triplicate in individuals with DS, the altered expression of miRNAs that target these genes could be modulating their expression pattern giving further support to the role of these miRNAs in DS pathogenesis.
Previous studies from our research group showed dysregulation of immune and inflammation-related genes non-located on HSA21 in children with DS [13, 58]. The study of immune-related genes [13] was carried out in the same cohort of individuals with DS and controls included in the present study; thus, we deemed it interesting to investigate a possible interaction between the immune-related DEGs previously identified and the six DEmiRs identified in the present study. Fourteen possible miRNA–target interactions were found (Table 3). These immune-related genes are involved in key immune response processes, such as regulation of T and B cells, production of cytokines, antibodies and autoantibodies, activation of signaling pathways, and their dysregulation could negatively influence the immune response in DS, as previously discussed [13].
Immune dysregulation has been widely described in DS [6, 59,60,61]. In line with the contribution of immunological alteration and infectious diseases to the pathophysiology of DS, individuals with T21 present increased morbidity and mortality due to respiratory tract infections [62, 63]. Since the start of the COVID-19 pandemic caused by SARS-CoV-2, there is concern about the individuals with DS being a risk group for the disease [64]. In fact, studies have shown that they present an increased risk for severe disease course and COVID-19–related hospitalization and death [65,66,67]. Recent network analysis of DS and SARS-CoV-2 revealed that the susceptibility to infection, worse prognosis, and complications of COVID-19 could be related to genetic factors associated with DS that predispose to unfavorable conditions, such as the dysregulation of genes involved in the immune response [68].
Our findings suggest that altered miRNAs expression patterns could be contributing to the immunological dysfunction observed in children with DS by targeting genes involved in several key processes crucial for the immune system to function properly. We emphasize the importance of investigating and understanding the immune system of DS individuals who constitute a clinically vulnerable population so that they can have better support and clinical management, and consequently be better in a future pandemic scenario. Considering that one coding gene can be regulated by multiple miRNAs [69] and one miRNA can regulate multiple targets [70], studies performing high throughput screening of miRNA and mRNA expression patterns could contribute to the understanding of the gene expression regulation in DS. Possible miRNA–target interactions identified in the present study are candidates for further investigation and experimental validation to understand the impact of DEmiRs on the mRNA expression of specific genes and consequently on the pathophysiology of DS.
Data availability
The data that support the findings of this study are openly available in Open Science Framework at http://doi.org/10.17605/OSF.IO/P2JBT, reference number P2JB
References
Antonarakis SE, Skotko BG, Rafii MS, et al. Down syndrome. Nat Rev Dis Primers. 2020. https://doi.org/10.1038/s41572-019-0143-7.
Pavarino Bertelli EC, Biselli JM, Bonfim D, Goloni-Bertollo EM. Clinical profile of children with Down syndrome treated in a genetics outpatient service in the southeast of Brazil. Rev Assoc Med Bras. 2009. https://doi.org/10.1590/S0104-42302009000500017.
Grieco J, Pulsifer M, Seligsohn K, Skotko B, Schwartz A. Down syndrome: cognitive and behavioral functioning across the lifespan. Am J Med Genet C Semin Med Genet. 2015. https://doi.org/10.1002/ajmg.c.31439.
Stoll C, Dott B, Alembik Y, Roth MP. Associated congenital anomalies among cases with Down syndrome. Eur J Med Genet. 2015. https://doi.org/10.1016/j.ejmg.2015.11.003.
Vacca RA, Bawari S, Valenti D, et al. Down syndrome: neurobiological alterations and therapeutic targets. Neurosci Biobehav Rev. 2019. https://doi.org/10.1016/j.neubiorev.2019.01.001.
Huggard D, Doherty DG, Molloy EJ. Immune dysregulation in children with down syndrome. Front Pediatr. 2020. https://doi.org/10.3389/fped.2020.00073.
Saida S. Predispositions to leukemia in down syndrome and other hereditary disorders. Curr Treat Options Oncol. 2017. https://doi.org/10.1007/s11864-017-0485-x.
Li CM, Guo M, Salas M, et al. Cell type-specific over-expression of chromosome 21 genes in fibroblasts and fetal hearts with trisomy 21. BMC Med Genet. 2006. https://doi.org/10.1186/1471-2350-7-24.
Gardiner K. Gene-dosage effects in Down syndrome and trisomic mouse models. Genome Biol. 2004. https://doi.org/10.1186/gb-2004-5-10-244.
FitzPatrick DR. Transcriptional consequences of autosomal trisomy: primary gene dosage with complex downstream effects. Trends Genet. 2005. https://doi.org/10.1016/j.tig.2005.02.012.
Letourneau A, Santoni FA, Bonilla X, et al. Domains of genome-wide gene expression dysregulation in Down’s syndrome. Nature. 2014. https://doi.org/10.1038/nature13200.
Antonarakis SE. Down syndrome and the complexity of genome dosage imbalance. Nat Rev Genet. 2017. https://doi.org/10.1038/nrg.2016.154.
Zampieri BL, Biselli-Périco JM, de Souza JE, et al. Altered expression of immune-related genes in children with Down syndrome. PLoS ONE. 2014. https://doi.org/10.1371/journal.pone.0107218.
Becker C, Hammerle-Fickinger A, Riedmaier I, Pfaffl MW. mRNA and microRNA quality control for RT-qPCR analysis. Methods. 2010. https://doi.org/10.1016/j.ymeth.2010.01.010.
The R Project for Statistical Computing. https://www.r-project.org/. Accessed 02 Ago 2021.
Gentleman RC, Carey VJ, Bates DM, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004. https://doi.org/10.1186/gb-2004-5-10-r80.
Dvinge H, Bertone P. HTqPCR: high-throughput analysis and visualization of quantitative real-time PCR data in R. Bioinformatics. 2009. https://doi.org/10.1093/bioinformatics/btp578.
Vandesompele J, De Preter K, Pattyn F, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002. https://doi.org/10.1186/gb-2002-3-7-research0034.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001. https://doi.org/10.1006/meth.2001.1262.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004. https://doi.org/10.2202/1544-6115.1027.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995; https://www.jstor.org/stable/2346101
Vlachos IS, Zagganas K, Paraskevopoulou MD, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015. https://doi.org/10.1093/nar/gkv403.
Paraskevopoulou MD, Georgakilas G, Kostoulas N, et al. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013. https://doi.org/10.1093/nar/gkt393.
Reczko M, Maragkakis M, Alexiou P, Grosse I, Hatzigeorgiou AG. Functional microRNA targets in protein coding sequences. Bioinformatics. 2012. https://doi.org/10.1093/bioinformatics/bts043.
miRBase: the microRNA database. Release 22.1. http://www.mirbase.org/. Accessed 02 Ago 2021.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006. https://doi.org/10.1093/nar/gkj112.
Alexandrov PN, Percy ME, Lukiw WJ. Chromosome 21-encoded microRNAs (mRNAs): impact on Down’s syndrome and Trisomy-21 linked disease. Cell Mol Neurobiol. 2018. https://doi.org/10.1007/s10571-017-0514-0.
Sethupathy P, Borel C, Gagnebin M, et al. Human microRNA-155 on chromosome 21 differentially interacts with its polymorphic target in the AGTR1 3’ untranslated region: a mechanism for functional single-nucleotide polymorphisms related to phenotypes. Am J Hum Genet. 2007. https://doi.org/10.1086/519979.
Tili E, Mezache L, Michaille JJ, et al. microRNA 155 up regulation in the CNS is strongly correlated to Down’s syndrome dementia. Ann Diagn Pathol. 2018. https://doi.org/10.1016/j.anndiagpath.2018.03.006.
Li YY, Alexandrov PN, Pogue AI, et al. miRNA-155 upregulation and complement factor H deficits in Down’s syndrome. NeuroReport. 2012. https://doi.org/10.1097/WNR.0b013e32834f4eb4.
Farroni C, Marasco E, Marcellini V, et al. Dysregulated miR-155 and miR-125b are related to impaired B-cell responses in Down syndrome. Front Immunol. 2018. https://doi.org/10.3389/fimmu.2018.02683.
Izzo A, Manco R, de Cristofaro T, et al. Overexpression of chromosome 21 miRNAs may affect mitochondrial function in the hearts of Down syndrome fetuses. Int J Genomics. 2017. https://doi.org/10.1155/2017/8737649.
Xu Y, Li W, Liu X, et al. Identification of dysregulated microRNAs in lymphocytes from children with Down syndrome. Gene. 2013. https://doi.org/10.1016/j.gene.2013.07.055.
Xu Y, Li W, Liu X, Ma H, Tu Z, Dai Y. Analysis of microRNA expression profile by small RNA sequencing in Down syndrome fetuses. Int J Mol Med. 2013. https://doi.org/10.3892/ijmm.2013.1499.
Moreira-Filho CA, Bando SY, Bertonha FB, et al. Modular transcriptional repertoire and MicroRNA target analyses characterize genomic dysregulation in the thymus of Down syndrome infants. Oncotarget. 2016. https://doi.org/10.18632/oncotarget.7120.
Shi WL, Liu ZZ, Wang HD, et al. Integrated miRNA and mRNA expression profiling in fetal hippocampus with Down syndrome. J Biomed Sci. 2016. https://doi.org/10.1186/s12929-016-0265-0.
Lin H, Sui W, Li W, et al. Integrated microRNA and protein expression analysis reveals novel microRNA regulation of targets in fetal down syndrome. Mol Med Rep. 2016. https://doi.org/10.3892/mmr.2016.5775.
Salvi A, Vezzoli M, Busatto S, et al. Analysis of a nanoparticle-enriched fraction of plasma reveals miRNA candidates for Down syndrome pathogenesis. Int J Mol Med. 2019. https://doi.org/10.3892/ijmm.2019.4158.
Shaham L, Vendramini E, Ge Y, et al. MicroRNA-486-5p is an erythroid oncomiR of the myeloid leukemias of Down syndrome. Blood. 2015. https://doi.org/10.1182/blood-2014-06-581892.
Arena A, Iyer AM, Milenkovic I, et al. Developmental expression and dysregulation of miR-146a and miR-155 in Down’s syndrome and mouse models of Down’s syndrome and Alzheimer’s disease. Curr Alzheimer Res. 2017. https://doi.org/10.2174/1567205014666170706112701.
Salemi M, Cannarella R, Marchese G, et al. Role of long non-coding RNAs in Down syndrome patients: a transcriptome analysis study. Hum Cell. 2021. https://doi.org/10.1007/s13577-021-00602-3.
Lopes-Ramos CM, Chen CY, Kuijjer ML, et al. Sex differences in gene expression and regulatory networks across 29 human tissues. Cell Rep. 2020. https://doi.org/10.1016/j.celrep.2020.107795.
Oliva M, Muñoz-Aguirre M, Kim-Hellmuth S, et al. The impact of sex on gene expression across human tissues. Science. 2020. https://doi.org/10.1126/science.aba3066.
Cui C, Yang W, Shi J, et al. Identification and analysis of human sex-biased MicroRNAs. Genomics Proteomics Bioinform. 2018. https://doi.org/10.1016/j.gpb.2018.03.004.
Vidya MK, Kumar VG, Sejian V, Bagath M, Krishnan G, Bhatta R. Toll-like receptors: significance, ligands, signaling pathways, and functions in mammals. Int Rev Immunol. 2018. https://doi.org/10.1080/08830185.2017.1380200.
Huggard D, Koay WJ, Kelly L, et al. Altered toll-like receptor signalling in children with Down syndrome. Mediators Inflamm. 2019. https://doi.org/10.1155/2019/4068734.
Sanjabi S, Oh SA, Li MO. Regulation of the immune response by TGF-β: from conception to autoimmunity and infection. Cold Spring Harb Perspect Biol. 2017. https://doi.org/10.1101/cshperspect.a022236.
Maeda H, Go H, Imamura T, Sato M, Momoi N, Hosoya M. Plasma TGF-β1 levels are elevated in Down syndrome infants with transient abnormal myelopoiesis. Tohoku J Exp Med. 2016. https://doi.org/10.1620/tjem.240.1.
Bromage SJ, Lang AK, Atkinson I, Searle RF. Abnormal TGFbeta levels in the amniotic fluid of Down syndrome pregnancies. Am J Reprod Immunol. 2000. https://doi.org/10.1111/j.8755-8920.2000.440403.x.
Yamauchi T, Moroishi T. Hippo pathway in mammalian adaptive immune system. Cells. 2019. https://doi.org/10.3390/cells8050398.
Boro M, Singh V, Balaji KN. Mycobacterium tuberculosis-triggered Hippo pathway orchestrates CXCL1/2 expression to modulate host immune responses. Sci Rep. 2016. https://doi.org/10.1038/srep37695.
Wang S, Zhou L, Ling L, et al. The crosstalk between hippo-YAP pathway and innate immunity. Front Immunol. 2020. https://doi.org/10.3389/fimmu.2020.00323.
Zhang X, Guo J, Wei X, et al. Bach1: function, regulation, and involvement in disease. Oxid Med Cell Longev. 2018. https://doi.org/10.1155/2018/1347969.
Mandal M, Maienschein-Cline M, Maffucci P, et al. BRWD1 orchestrates epigenetic landscape of late B lymphopoiesis. Nat Commun. 2018. https://doi.org/10.1038/s41467-018-06165-6.
Ng AP, Coughlan HD, Hediyeh-Zadeh S, et al. An Erg-driven transcriptional program controls B cell lymphopoiesis. Nat Commun. 2020. https://doi.org/10.1038/s41467-020-16828-y.
Bhandage AK, Jin Z, Korol SV, et al. GABPA regulates release of inflammatory cytokines from peripheral blood mononuclear cells and CD4+ T cells and is immunosuppressive in type 1 diabetes. EBioMedicine. 2018. https://doi.org/10.1016/j.ebiom.2018.03.019.
Bellissimo DC, Chen CH, Zhu Q, et al. Runx1 negatively regulates inflammatory cytokine production by neutrophils in response to Toll-like receptor signaling. Blood Adv. 2020. https://doi.org/10.1182/bloodadvances.2019000785.
Silva CR, Biselli-Périco JM, Zampieri BL, et al. Differential expression of inflammation-related genes in children with down syndrome. Mediators Inflamm. 2016. https://doi.org/10.1155/2016/6985903.
Ram G, Chinen J. Infections and immunodeficiency in Down syndrome. Clin Exp Immunol. 2011. https://doi.org/10.1111/j.1365-2249.2011.04335.x.
Waugh KA, Araya P, Pandey A, et al. Mass cytometry reveals global immune remodeling with multi-lineage hypersensitivity to type i interferon in down syndrome. Cell Rep. 2019. https://doi.org/10.1016/j.celrep.2019.10.038.
Illouz T, Biragyn A, Iulita MF, et al. Immune dysregulation and the increased risk of complications and mortality following respiratory tract infections in adults with down syndrome. Front Immunol. 2021. https://doi.org/10.3389/fimmu.2021.621440.
O’Leary L, Hughes-McCormack L, Dunn K, Cooper SA. Early death and causes of death of people with Down syndrome: a systematic review. J Appl Res Intellect Disabil. 2018. https://doi.org/10.1111/jar.12446.
Landes SD, Stevens JD, Turk MA. Cause of death in adults with Down syndrome in the United States. Disabil Health J. 2020. https://doi.org/10.1016/j.dhjo.2020.100947.
Espinosa JM. Down syndrome and COVID-19: a perfect storm? Cell Rep Med. 2020. https://doi.org/10.1016/j.xcrm.2020.100019.
Malle L, Gao C, Hur C, et al. Individuals with Down syndrome hospitalized with COVID-19 have more severe disease. Genet Med. 2021. https://doi.org/10.1038/s41436-020-01004-w.
Clift AK, Coupland CAC, Keogh RH, Hemingway H, Hippisley-Cox J. COVID-19 mortality risk in Down syndrome: results from a cohort study of 8 million adults. Ann Intern Med. 2021. https://doi.org/10.7326/M20-4986.
Hüls A, Costa ACS, Dierssen M, et al. T21RS COVID-19 Initiative. Medical vulnerability of individuals with Down syndrome to severe COVID-19-data from the Trisomy 21 Research Society and the UK ISARIC4C survey. EClinicalMedicine. 2021. https://doi.org/10.1016/j.eclinm.2021.100769.
De Toma I, Dierssen M. Network analysis of Down syndrome and SARS-CoV-2 identifies risk and protective factors for COVID-19. Sci Rep. 2021. https://doi.org/10.1038/s41598-021-81451-w.
Wu S, Huang S, Ding J, et al. Multiple microRNAs modulate p21Cip1/Waf1 expression by directly targeting its 3’ untranslated region. Oncogene. 2010. https://doi.org/10.1038/onc.2010.34.
Zheng F, Liao YJ, Cai MY, et al. Systemic delivery of microRNA-101 potently inhibits hepatocellular carcinoma in vivo by repressing multiple targets. PLoS Genet. 2015. https://doi.org/10.1371/journal.pgen.1004873.
Acknowledgements
The authors thank Profa. Dra. Dirce Maria Carraro and Profa. Dra. Elisa Napolitano e Ferreira of Laboratory of Genomics and Molecular Biology, A.C. Camargo Hospital, São Paulo, Brazil, for providing the equipment 7900HT Fast Real-Time PCR System (Applied Biosystems) and technical assistance, and to the Hospital de Base of Sao Jose do Rio Preto and to IBILCE-UNESP for supporting the recruitment of participants. This work was funded by the São Paulo Research Foundation—FAPESP (2010/ 00153-3), Coordination for the Improvement of Higher Education Personnel (CAPES), and The Brazilian National Council for Scientific and Technological Development—CNPq (305641/2011-5).
Funding
This work was supported by the São Paulo Research Foundation—FAPESP (2010/00153-3). JMB and BLZ received fellowship from Coordination for the Improvement of Higher Education Personnel (CAPES). ECP is a fellow researcher of The Brazilian National Council for Scientific and Technological Development—CNPq (305641/2011-5).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Ethical approval
The study protocol was approved by the Research Ethics Committee of Sao José do Rio Preto Medical School (CEP-FAMERP, protocol 3340/2004 and CAAE number: 88278218.4.0000.5415), in Sao Paulo State, Brazil, and was in accordance with the Declaration of Helsinki ethical principles for medical research involving human subjects, as well as the Brazilian ethical guidelines for research involving human participants (National Health Council, Resolution 466/2012 and Resolution 510/2016).
Informed consent
Written informed consent was obtained from parents/legal guardians of all children included in the study.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Biselli, J.M., Zampieri, B.L., Biselli-Chicote, P.M. et al. Differential microRNA expression profile in blood of children with Down syndrome suggests a role in immunological dysfunction. Human Cell 35, 639–648 (2022). https://doi.org/10.1007/s13577-022-00672-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13577-022-00672-x