Modular transcriptional repertoire and MicroRNA target analyses characterize genomic dysregulation in the thymus of Down syndrome infants

Differential microRNA expression profile in blood of children with Down syndrome suggests a role in immunological dysfunction

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.

Immune Dysregulation and the Increased Risk of Complications and Mortality Following Respiratory Tract Infections in Adults With Down Syndrome

The risk of severe outcomes following respiratory tract infections is significantly increased in individuals over 60 years, especially in those with chronic medical conditions, i.e., hypertension, diabetes, cardiovascular disease, dementia, chronic respiratory disease, and cancer. Down Syndrome (DS), the most prevalent intellectual disability, is caused by trisomy-21 in ~1:750 live births worldwide. Over the past few decades, a substantial body of evidence has accumulated, pointing at the occurrence of alterations, impairments, and subsequently dysfunction of the various components of the immune system in individuals with DS. This associates with increased vulnerability to respiratory tract infections in this population, such as the influenza virus, respiratory syncytial virus, SARS-CoV-2 (COVID-19), and bacterial pneumonias. To emphasize this link, here we comprehensively review the immunobiology of DS and its contribution to higher susceptibility to severe illness and mortality from respiratory tract infections.

Premature Senescence and Increased Oxidative Stress in the Thymus of Down Syndrome Patients

Down syndrome (DS) patients prematurely show clinical manifestations usually associated with aging. Their immune system declines earlier than healthy individuals, leading to increased susceptibility to infections and higher incidence of autoimmune phenomena. Clinical features of accelerated aging indicate that trisomy 21 increases the biological age of tissues. Based on previous studies suggesting immune senescence in DS, we hypothesized that induction of cellular senescence may contribute to early thymic involution and immune dysregulation. Immunohistochemical analysis of thymic tissue showed signs of accelerated thymic aging in DS patients, normally seen in older healthy subjects. Moreover, our whole transcriptomic analysis on human Epcam-enriched thymic epithelial cells (hTEC), isolated from three DS children, which revealed disease-specific transcriptomic alterations. Gene set enrichment analysis (GSEA) of DS TEC revealed an enrichment in genes involved in cellular response to stress, epigenetic histone DNA modifications and senescence. Analysis of senescent markers and oxidative stress in hTEC and thymocytes confirmed these findings. We detected senescence features in DS TEC, thymocytes and peripheral T cells, such as increased β-galactosidase activity, increased levels of the cell cycle inhibitor p16, telomere length and integrity markers and increased levels of reactive oxygen species (ROS), all factors contributing to cellular damage. In conclusion, our findings support the key role of cellular senescence in the pathogenesis of immune defect in DS while adding new players, such as epigenetic regulation and increased oxidative stress, to the pathogenesis of immune dysregulation.

Network analysis of Down syndrome and SARS-CoV-2 identifies risk and protective factors for COVID-19

SARS-CoV-2 infection has spread uncontrollably worldwide while it remains unknown how vulnerable populations, such as Down syndrome (DS) individuals are affected by the COVID-19 pandemic. Individuals with DS have more risk of infections with respiratory complications and present signs of auto-inflammation. They also present with multiple comorbidities that are associated with poorer COVID-19 prognosis in the general population. All this might place DS individuals at higher risk of SARS-CoV-2 infection or poorer clinical outcomes. In order to get insight into the interplay between DS genes and SARS-cov2 infection and pathogenesis we identified the genes associated with the molecular pathways involved in COVID-19 and the host proteins interacting with viral proteins from SARS-CoV-2. We then analyzed the overlaps of these genes with HSA21 genes, HSA21 interactors and other genes consistently differentially expressed in DS (using public transcriptomic datasets) and created a DS-SARS-CoV-2 network. We detected COVID-19 protective and risk factors among HSA21 genes and interactors and/or DS deregulated genes that might affect the susceptibility of individuals with DS both at the infection stage and in the progression to acute respiratory distress syndrome. Our analysis suggests that at the infection stage DS individuals might be more susceptible to infection due to triplication of TMPRSS2, that primes the viral S protein for entry in the host cells. However, as the anti-viral interferon I signaling is also upregulated in DS, this might increase the initial anti-viral response, inhibiting viral genome release, viral replication and viral assembly. In the second pro-inflammatory immunopathogenic phase of the infection, the prognosis for DS patients might worsen due to upregulation of inflammatory genes that might favor the typical cytokine storm of COVID-19. We also detected strong downregulation of the NLRP3 gene, critical for maintenance of homeostasis against pathogenic infections, possibly leading to bacterial infection complications.

Inborn Errors of Adaptive Immunity in Down Syndrome

Down syndrome fits an immunophenotype of combined immunodeficiency with immunodysregulation, manifesting with increased susceptibility to infections, autoimmunity, autoinflammatory diseases, and hematologic malignancies. Qualitative and quantitative alterations in innate and adaptive immunity are found in most individuals with Down syndrome. However, there is substantial heterogeneity and no correlation between immunophenotype and clinical presentation. Previously, it was thought that the immunological changes in Down syndrome were caused by precocious aging. We emphasize in this review that the immune system in Down syndrome is intrinsically different from the very beginning. The overexpression of specific genes located on chromosome 21 contributes to immunodeficiency and immunodysregulation, but gene expression differs between genes located on chromosome 21 and depends on tissue and cell type. In addition, trisomy 21 results in gene dysregulation of the whole genome, reflecting the complex nature of this syndrome in comparison to well-known inborn errors of immunity that result from monogenic germline mutations. In this review, we provide an updated overview focusing on inborn errors of adaptive immunity in Down syndrome.

Cystathionine-β-Synthase: Molecular Regulation and Pharmacological Inhibition

Cystathionine-β-synthase (CBS), the first (and rate-limiting) enzyme in the transsulfuration pathway, is an important mammalian enzyme in health and disease. Its biochemical functions under physiological conditions include the metabolism of homocysteine (a cytotoxic molecule and cardiovascular risk factor) and the generation of hydrogen sulfide (H2S), a gaseous biological mediator with multiple regulatory roles in the vascular, nervous, and immune system. CBS is up-regulated in several diseases, including Down syndrome and many forms of cancer; in these conditions, the preclinical data indicate that inhibition or inactivation of CBS exerts beneficial effects. This article overviews the current information on the expression, tissue distribution, physiological roles, and biochemistry of CBS, followed by a comprehensive overview of direct and indirect approaches to inhibit the enzyme. Among the small-molecule CBS inhibitors, the review highlights the specificity and selectivity problems related to many of the commonly used "CBS inhibitors" (e.g., aminooxyacetic acid) and provides a comprehensive review of their pharmacological actions under physiological conditions and in various disease models.

Role of 3-Mercaptopyruvate Sulfurtransferase in the Regulation of Proliferation and Cellular Bioenergetics in Human Down Syndrome Fibroblasts

Down syndrome (trisomy of human chromosome 21) is a common genetic disorder. Overproduction of the gaseous mediator hydrogen sulfide (H₂S) has been implicated in the pathogenesis of neurological and metabolic deficits associated with Down syndrome. Several lines of data indicate that an important enzyme responsible for H₂S overproduction in Down syndrome is cystathionine-β-synthase (CBS), an enzyme localized on chromosome 21. The current study explored the possibility that a second H₂S-producing enzyme, 3-mercaptopyruvate sulfurtransferase (3-MST), may also contribute to the development of functional deficits of Down syndrome cells. Western blotting analysis demonstrated a significantly higher level of 3-MST protein expression in human Down syndrome fibroblasts compared to cells from healthy control individuals; the excess 3-MST was mainly localized to the mitochondrial compartment. Pharmacological inhibition of 3-MST activity improved mitochondrial electron transport and oxidative phosphorylation parameters (but did not affect the suppressed glycolytic parameters) and enhanced cell proliferation in Down syndrome cells (but not in healthy control cells). The findings presented in the current report suggest that in addition to the indisputable role of CBS, H₂S produced from 3-MST may also contribute to the development of mitochondrial metabolic and functional impairments in Down syndrome cells.

Age-related transcriptional modules and TF-miRNA-mRNA interactions in neonatal and infant human thymus

The human thymus suffers a transient neonatal involution, recovers and then starts a process of decline between the 1^st and 2^nd years of life. Age-related morphological changes in thymus were extensively investigated, but the genomic mechanisms underlying this process remain largely unknown. Through Weighted Gene Co-expression Network Analysis (WGCNA) and TF-miRNA-mRNA integrative analysis we studied the transcriptome of neonate and infant thymic tissues grouped by age: 0–30 days (A); 31days-6 months (B); 7–12 months (C); 13–18 months (D); 19-31months (E). Age-related transcriptional modules, hubs and high gene significance (HGS) genes were identified, as well as TF-miRNA-hub/HGS co-expression correlations. Three transcriptional modules were correlated with A and/or E groups. Hubs were mostly related to cellular/metabolic processes; few were differentially expressed (DE) or related to T-cell development. Inversely, HGS genes in groups A and E were mostly DE. In A (neonate) one third of the hyper-expressed HGS genes were related to T-cell development, against one-twentieth in E, what may correlate with the early neonatal depletion and recovery of thymic T-cell populations. This genomic mechanism is tightly regulated by TF-miRNA-hub/HGS interactions that differentially govern cellular and molecular processes involved in the functioning of the neonate thymus and in the beginning of thymic decline.

Identification of a thymus microRNA‑mRNA regulatory network in Down syndrome

The present bioinformatics analysis was performed using a multi‑step approach to identify a microRNA (miR)‑mRNA regulatory network in Down syndrome. miR (GSE69210) and mRNA (GSE70573) data was downloaded and collected from the thymic tissues of both Down syndrome and karyotypically normal subjects and placed in a public repository. Then, weighted gene co‑expression network analysis (WGCNA) was performed to screen for miRs and mRNAs associated with Down syndrome. Subsequently, differentially expressed miRs (DEmiRs) and mRNAs/differentially expressed genes (DEGs) were identified following screening and mapping to RNA data. Bidirectional hierarchical clustering analysis was then performed to distinguish DEmiRs and DEGs between Down syndrome samples and normal control samples. DEmiR targets were retrieved using the miRanda database and mapped to the mRNA module screen by WGCNA. A gene co‑expression network was constructed and subjected to functional enrichment analysis. During WGCNA, a total of 6 miR modules and 20 mRNA modules associated with Down syndrome were identified. Following mapping of these miRs and mRNAs to the miR and mRNA modules screened using WGNCA, a total of 12 DEmiRs and 237 DEGs were collected. Following comparison with DEmiR targets retrieved from the miRanda database, a total of 255 DEmiR‑DEG pairs, including 6 DEmiRs and 106 DEGs were obtained. At expression correlation coefficient >0.9, a total of 231 gene pairs were selected. These gene pairs were enriched in response to stress and response to stimuli following functional annotation and module division. An integrated analysis of miR and mRNA expression in the thymus in Down syndrome is reported in the present study. miR‑30c, miR‑145, miR‑183 and their targets may serve important roles in the pathogenesis and development of complications in Down syndrome. However, further experimental studies are required to verify these results.

iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data

BACKGROUND

RNA-seq is widely used for transcriptomic profiling, but the bioinformatics analysis of resultant data can be time-consuming and challenging, especially for biologists. We aim to streamline the bioinformatic analyses of gene-level data by developing a user-friendly, interactive web application for exploratory data analysis, differential expression, and pathway analysis.

RESULTS

iDEP (integrated Differential Expression and Pathway analysis) seamlessly connects 63 R/Bioconductor packages, 2 web services, and comprehensive annotation and pathway databases for 220 plant and animal species. The workflow can be reproduced by downloading customized R code and related pathway files. As an example, we analyzed an RNA-Seq dataset of lung fibroblasts with Hoxa1 knockdown and revealed the possible roles of SP1 and E2F1 and their target genes, including microRNAs, in blocking G1/S transition. In another example, our analysis shows that in mouse B cells without functional p53, ionizing radiation activates the MYC pathway and its downstream genes involved in cell proliferation, ribosome biogenesis, and non-coding RNA metabolism. In wildtype B cells, radiation induces p53-mediated apoptosis and DNA repair while suppressing the target genes of MYC and E2F1, and leads to growth and cell cycle arrest. iDEP helps unveil the multifaceted functions of p53 and the possible involvement of several microRNAs such as miR-92a, miR-504, and miR-30a. In both examples, we validated known molecular pathways and generated novel, testable hypotheses.

CONCLUSIONS

Combining comprehensive analytic functionalities with massive annotation databases, iDEP ( http://ge-lab.org/idep/ ) enables biologists to easily translate transcriptomic and proteomic data into actionable insights.

Minipuberty and Sexual Dimorphism in the Infant Human Thymus

AIRE expression in thymus is downregulated by estrogen after puberty, what probably renders women more susceptible to autoimmune disorders. Here we investigated the effects of minipuberty on male and female infant human thymic tissue in order to verify if this initial transient increase in sex hormones - along the first six months of life - could affect thymic transcriptional network regulation and AIRE expression. Gene co-expression network analysis for differentially expressed genes and miRNA-target analysis revealed sex differences in thymic tissue during minipuberty, but such differences were not detected in the thymic tissue of infants aged 7-18 months, i.e. the non-puberty group. AIRE expression was essentially the same in both sexes in minipuberty and in non-puberty groups, as assessed by genomic and immunohistochemical assays. However, AIRE-interactors networks showed several differences in all groups regarding gene-gene expression correlation. Therefore, minipuberty and genomic mechanisms interact in shaping thymic sexual dimorphism along the first six months of life.