- Research
- Open access
- Published:
Transcriptomic analysis of stem cells from chorionic villi uncovers the impact of chromosomes 2, 6 and 22 in the clinical manifestations of Down syndrome
Stem Cell Research & Therapy volume 14, Article number: 265 (2023)
Abstract
Background
Down syndrome (DS) clinical multisystem condition is generally considered the result of a genetic imbalance generated by the extra copy of chromosome 21. Recent discoveries, however, demonstrate that the molecular mechanisms activated in DS compared to euploid individuals are more complex than previously thought. Here, we utilize mesenchymal stem cells from chorionic villi (CV) to uncover the role of comprehensive functional genomics-based understanding of DS complexity.
Methods
Next-generation sequencing coupled with bioinformatic analysis was performed on CV obtained from women carrying fetuses with DS (DS-CV) to reveal specific genome-wide transcriptional changes compared to their euploid counterparts. Functional assays were carried out to confirm the biological processes identified as enriched in DS-CV compared to CV (i.e., cell cycle, proliferation features, immunosuppression and ROS production).
Results
Genes located on chromosomes other than the canonical 21 (Ch. 2, 6 and 22) are responsible for the impairment of life-essential pathways, including cell cycle regulation, innate immune response and reaction to external stimuli were found to be differentially expressed in DS-CV. Experimental validation confirmed the key role of the biological pathways regulated by those genes in the etiology of such a multisystem condition.
Conclusions
NGS dataset generated in this study highlights the compromised functionality in the proliferative rate and in the innate response of DS-associated clinical conditions and identifies DS-CV as suitable tools for the development of specifically tailored, personalized intervention modalities.
Background
Down syndrome (DS) is the most common chromosome-related disorder, with an incidence of 1 in 1000 live births worldwide according to the World Health Organization [1]. As a multisystem disease, DS is caused, in most cases, by meiotic non-disjunction of the maternal chromosome 21 and results in a third copy of the entire chromosome 21 in somatic cells [2,3,4]. The extra copy of the chromosome 21 gene set and the consequent gene dosage imbalance results in altered morphological development, leading to dysmorphic features (e.g., craniofacial abnormalities, hypotonia and cognitive impairment) [5] and often, associated pathological conditions. Common conditions include respiratory tract infections, autoimmune disorders, diabetes mellitus, celiac disease, thyroid impairment, gastrointestinal cardiac defects and, in some cases, leukemia [6].
While DS has been recognized as a clinical entity for about 150 years, effective therapies to improve speech and articulation, prevent cognitive decline, recover immune function and to prevent leukemia initiation remain to be developed for DS patients [7]. Only recently the in-depth molecular analyses have begun to reveal the utility of using molecular genetics and genomics as a promising exploratory approach to better understand DS pathogenesis [6, 8, 9]. Specifically, the current thrust of research on DS is based on the premise that understanding mechanisms at the genetic and molecular level will provide a rational basis for the development of effective therapeutic and preventive interventions [10]. The transcriptome of DS-derived cells offers a more comprehensive overview of the molecular traits associated to DS and could assist in the identification of candidate genes (within or outside of chromosome 21) to be potentially considered as therapeutic targets [11].
To date, molecular studies have been primarily focused on the gain of function genes secondary to the extra chromosome 21 and their roles in determining the specific DS phenotype, compared to an euploid individuals [12, 13]. Several studies identified differentially expressed genes (DEGs) spread over all of the chromosomes, not limited to chromosome 21, which contribute to the DS phenotype [8, 14, 15]. Among them, a study identified complex intrachromosomal functional gene re-adjustments with 65 up-regulated and 111 down-regulated genes in chromosome 21 [16], which include those that are related to central nervous system development (8 genes) and extracellular matrix organization (11 genes), influencing the normal development and communication among neuronal cells. Another study has shown that the trisomy not only determines an imbalance in gene expression associated to chromosome 21, but it also induced differences in gene transcripts situated on other chromosomes [13]. The transcriptomic analysis in monozygotic twins discordant for trisomy 21 revealed significant gene expression changes in other chromosome (namely 1, 6, 11, 16, 17 and 19), indicating that the molecular mechanisms activated in DS compared to euploid individuals is more complex than previously thought [16]. DS studies in human patients are limited for obvious ethical reasons, resulting in the use of differential models that have been proposed to study this syndrome. Mouse models or induced pluripotent stem cells (iPSCs) cells are a few examples of the strategies developed so far to provide new insights into the DS phenotype and represent suitable tools for the identification of potential targets for therapeutic options. On the other hand, modified iPSCs present therapeutic limitations related to the viral vector integration into the host genome and causing genomic and epigenetic changes [17]. Mesenchymal stem cells (MSCs) from extraembryonic tissues are easily collected before, during and after the birth, without requiring invasive procedures or experiencing ethical limitation in comparison with those associated with their adult and embryonic counterparts [18]. In this study, we propose MSCs isolated from a gestational tissue, the chorionic villi (CV), sampled during prenatal diagnosis, for a personalized and comprehensive functional genomics-based understanding of this complex clinical condition which will provide new knowledge-based insights for future clinical treatment of DS. We used next-generation sequencing (NGS) analysis to compare the whole transcriptomic profiles of MSC populations isolated from CV derived from women carrying fetuses diagnosed with DS (DS-CV) with whole transcriptomes of their euploid counterparts (CV). Experimental validation of the biological processes identified as enriched in DS-CV compared to CV included the assessment of crucial pathways related to the typical MSC phenotype (i.e., cell cycle and proliferation features) as well as those responsible for the marked reactivity to external stimuli (e.g., production of reactive oxygen species, ROS) and enhanced immunosuppressive capability attributed to DS-CV.
Methods
Establishment of cell cultures
Chorionic villi (CV) samples were obtained from pregnant women undergoing prenatal diagnosis between the 10th and the 13th week of gestation at the Cytogenetic Laboratory Children's Hospital Salesi (Ancona, Italy) upon informed written consent for the use of tissue for research purposes. The study was approved by the Regional Institutional Review Board (Comitato Etico Regione Marche) and was conducted in accordance with the Declaration of Helsinki. CV were carefully separated from maternal decidua using sterile fine forceps to avoid the contamination and washed with 1 × phosphate-buffered saline (PBS) to remove any blood clots and subjected to mechanical treatment to obtain small fragments (3 mm2). Finally, the sample was transferred into flasks where the cells were allowed to adhere and migrate out of the tissue. Primary MSCs isolated from healthy women (CV) and women carrying fetuses with diagnosed DS (DS-CV) were expanded in standard culture medium consisting of high-glucose Dulbecco's modified eagle medium (HG-DMEM, Corning) supplemented with 10% Fetal Bovine Serum (FBS), 1% L-glutamine and 2% antimitotic-antibiotics (A/A, EuroClone S.p.A). Human monocytic cells (THP-1 cells) used in this study were purchased from ATCC. THP-1 cell culture was established in Roswell Park Memorial Institute-1640 (RPMI-1640) media (Gibco) supplemented with 10% FBS (Corning), 1% L-glutamine and 2% A/A. Cell cultures were maintained at 5% CO2 atmosphere and 37◦C temperature.
RNA sequence analysis performed using next-generation sequencing (NGS)
RNA extraction and quality check
Total RNA was extracted from CV- and DS-CV at passage 2 (P2) using total RNA Purification Plus Kit (Norgen) and processed by the Functional Genomic Center (University of Verona). RNA concentration and integrity were assessed using the RNA 6000 Nano Kit (Agilent Technologies). RNA samples that showed an integrity number (RIN) > 9 were used for the NGS.
RNA-Seq library preparation
RNA-Seq library preparation was performed using the TruSeq stranded mRNA kit (Illumina) from 1 μg of RNA per sample. Library size was assessed by capillary electrophoretic analysis with the Agilent 4200 Tape station. RNA libraries were analyzed on an Illumina NextSeq 500 sequencer using 75nt single reads. The quality of the reads was checked using the software FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), discarding those reporting more than 50 bp with low scores. Subsequently, through the software Scythe (https://github.com/vsbuffalo/scythe), adaptor contamination was removed, and the reads showing low quality were modified at their endings using the software Sickle (https://github.com/vsbuffalo/sickle). Filtered reads were aligned to the Human reference genome GRCh38 using HISAT2 [19]. Differential analysis was performed using the R bioconductor package DESeq2 [20], contrasting the 3 CV samples with the 3 ds-CV samples. Raw and processed data are deposited in the Gene Expression Omnibus (GEO) database with accession number GSE184450.
Heatmap and volcano plot
To visualize the differentially expressed genes (DEGs) between CV and DS-CV, heatmap and volcano plot were created using ggplot2 package [21]. Chord diagrams were generated using the GOPlot R package [22]. The criteria of DEGs were adjusted p-value (Padj) < 0.05. In the heatmap and in the volcano plot, the DEGs are displayed as color-coded; red represents over-expression while blue under-expression. In the heatmap the Euclidean distance of samples was used, while in the volcano plot mean of normalized counts was used.
Pathway and gene ontology
Pathway analysis on all significantly dysregulated genes was performed using the WebGestaltR [23] package, with gene set enrichment analysis (GSEA) [24] enrichment method. The whole genome was used as the background, with gene sets REACTOME pathway [25], KEGG pathway [26], GO biological processes [27] and CytogenicBand analyzed for enrichment.
Protein–protein interactions
Protein–protein interactions were obtained using the Cytoscape StringApp [28], which uses STRING-DB [29] for protein–protein interaction data. For each pair of chromosomes considered, the significantly dysregulated genes were imported into the StringApp, which automatically generated protein networks based on interactions between proteins annotated in STRING-DB. The largest cluster of interacting proteins was selected, and enrichment analysis was performed within the StringApp, for Gene Ontology Biological Process [27], Gene Ontology Cellular Component [27] and REACTOME pathway sets [30], with the whole genome used as the background.
Characterization of MSCs from CV and DS-CV
Morphology
Morphology was evaluated through fluorescence microscopy (OLYMPUS BX51, equipped with the Spot Advanced software). ActinGreen (Life Technologies) and Hoechst (Sigma-Aldrich) were used to highlight the cytoskeleton and the nucleus, respectively. Briefly, cells were seeded at 1.4 × 104/well in 8 chamber-slides (Corning® BioCoat™ CultureSlides) and let adhere. Cells were then washed twice in PBS and immediately fixed with 4% paraformaldehyde (PFA) in PBS for 15 min. Subsequently, cells were washed for three times in PBS and permeabilized with 0.1% Triton X-100 (Sigma-Aldrich) (in PBS) for 10–15 min. After a second round of washing, MSCs were blocked in 1% bovine serum albumin (BSA)/PBS for 30 min at room temperature. Cells were then incubated with ActinGreen for 20 min at room temperature in the dark. Five minutes before visualization, cells were incubated with Hoechst (1 μg/ml).
Expression of MSC-associated markers by flow cytometry
Cell phenotype was analyzed according to the published criteria [31]. A total of 5 × 105 cells (at P2) were collected, fixed with 75% ethanol and incubated for 20 min with 0.5% BSA/PBS to block non-specific sites. Samples were then centrifuged at 500 × g for 5 min at 20 °C and stained with fluorescently labeled antibodies. Antibodies include phycoerythrin (PE)-conjugated anti-ecto-5′-nucleotidase (PE-CD73; BioLegend), fluorescein isothiocyanate (FITC)-conjugated anti-thymocyte differentiation antigen 1 (FITC-CD90; BioLegend) and fluorescein isothiocyanate (FITC)-conjugated anti-glycoprotein CD44 (FITC-BioLegend), and allophycocyanin (APC)-conjugated anti-integrin β1 (APC-CD29; BioLegend). Cells were incubated for 45 min at room temperature in the dark, then washed twice with filtered PBS to remove the excess of antibody and analyzed using Guava Easycyte Millipore flow cytometer with GUAVASOFT 2.2.3 software.
Differentiative potential
To test their multipotent differentiative potential, CV and DS-CV at P3 were seeded at a density of 1 × 103/cm2 in six-well tissue culture dishes. For osteogenesis, cells were cultured in HG-DMEM, supplemented with 10% FBS, 100 U/ml penicillin, 100 μg/ml streptomycin, 0.25 μg/ml amphotericin B, 2 mM/l L-glutamine, 10 mM β-glycerophosphate (Sigma, 50020), 0.1 μM dexamethasone (Sigma, D2915) and 250 μM ascorbic acid (Sigma, A8960). For adipogenic differentiation, cells were cultured in HG-DMEM, supplemented with 10% FBS, 100 U/ml penicillin, 100 μg/ml streptomycin, 0.25 μg/ml amphotericin B, 2 mM/l L-glutamine, 10 μg/ml insulin (Sigma I-6634), 150 μM indomethacin (Sigma I-7378), 1 μM dexamethasone and 500 μM IBMX (3-isobutyl-methyl-xanthine, Sigma I-7018). Induced cells were incubated for two weeks at 38.5 °C with 5% CO2. Non-induced control cells were cultured for the same time with standard medium (HG-DMEM supplemented with 10% FBS, 100 U/ml penicillin, 100 μg/ml streptomycin, 0.25 μg/ml amphotericin B, 2 mM/l L-glutamine). The presence of calcium deposits in differentiated cells was verified after 7 and 14 days of induction by von Kossa staining, whereas Oil Red-O staining (Sigma, O0625) was used to identify lipid droplets in the cytoplasm. Expression of specific genes was performed as reported in section “Gene expression analysis” to confirm occurred differentiation.
Clonogenicity
CV and DS-CV-derived MSCs (at P0) were seeded at different densities (namely, 3.5 × 102, 1 × 103, 3.5 × 103 cells/cm2) in six-well plates for 2 weeks in stem cell standard culture medium. At the end of the 2 weeks culture period, colonies were fixed with 1% paraformaldehyde (PFA), stained with Giemsa at room temperature (RT) and washed twice. Colonies formed by 15–20 nucleated were counted with an inverted microscope (Meiji Techno).
Proliferation of MSCs from CV and DS-CV
Growth curve were produced by seeding CV- and DS-CV (at P2) at the density 9.5 × 103 cells/well into six-well tissue culture polystyrene dishes (EuroClone). Every 2 days, over a 15-day culture period, cells were trypsinized and counted as previously reported [32]. Analysis was performed in triplicate.
Cell cycle analysis
To determine possible differences in cell cycle between MSCs from CV- and DS-CV, the distribution of cells in the major phases of the cycle (subG0, G1, S, G2/M) was examined using flow cytometry. One million cells (per group) were harvested (at P2) and were fixed by overnight incubation with cold 70% ethanol. Cells were subsequently resuspended with a staining solution containing propidium iodide (PI, 40 μg/ml) and RNase (100 μg/ml) in PBS. At least 20,000 events/sample were acquired through a Guava Millipore cytometer. The percentage of cells in each phase was established using FlowJo software, through a red channel.
Assessment of CV and DS-CV immunosuppressive potential
To assess the immunosuppressive potential of CV- and DS-CV, cells (at P2) were seeded at the density of 2 × 104 in 24-well plates and cultured for 24 h at 37 °C before stimulation with pro-inflammatory cytokines. Soluble recombinant human tumor necrosis factor-alpha (TNF-α) and interferon gamma (IFN-γ) (Peprotech) were used to simulate an inflammatory environment at the concentration of 20 ng/ml. After 24- and 48-h cells were harvested and processed for molecular analysis as reported in Gene expression analysis paragraph.
Reactive oxygen species (ROS) detection on MSCs from CV and DS-CV
Cellular oxidative stress was evaluated through the reactive 2′,7′-dichlorodihydrofluorecein (H2-DCFDA, Sigma-Aldrich) assay. CV- and DS-CV were trypsinized, washed two times with cold PBS, suspended with PBS containing the probe at the working concentration of 10 μM and incubated for 30 min at 37°. Cells were subsequently washed twice in 1 × PBS and stained with 10 μg/ml PI to exclude dead cells from the analysis. The mean fluorescence intensity of cells emitting in the green (H2-DCFDA) and red (PI) channel was analyzed using Guava Easycyte Millipore flow cytometer with GUAVASOFT 2.2.3 software.
Assessment of CV and DS-CV immunomodulatory potential
CV and DS-CV MSCs (at P2) were seeded at the density of 3 × 103 cells/cm2 in T75 flasks and cultured to 70–80% confluence. Cells were then washed thoroughly 3 times with Hank's balanced salt solution (HBSS) (without Ca2+ and Mg2+) (Gibco) and replenished with standard culture medium supplemented with 10% exosome depleted FBS (Gibco). Cells were maintained in standard conditions (5% CO2; 37 °C) for 24 and 48 h. The media conditioned by CV and DS-CV MSCs was then collected, centrifuged at 500 × g for 10 min at 4 °C to remove cell debris and frozen in aliquots at − 80 °C until further analyses. The role CV and DS-CV MSCs play in modulating the behavior/phenotype of immune cells was assessed on monocytic cells THP-1. THP-1 cells were plated at density 5 × 105 in 24-well plates and exposed for 24 and 48 h to the conditioned media produced by CV and DS-CV. At each time point cells were harvested and processed for molecular analysis, as reported in the gene expression analysis paragraph.
Gene expression analysis
Gene expression evaluation was performed using the human specific oligonucleotide primers described in Table 1. Primers were designed using open-source Primer-BLAST (National Center for Biotechnology Information, NCBI) across an exon–exon junction to avoid genomic DNA amplification. Manual corrections were made to improve specificity. Total RNA was isolated using TRI-reagent (Sigma) and purified from genomic DNA through digestion with DNAse (Sigma). RNA concentration and purity were measured using a NanoDrop ND1000 spectrophotometer (NanoDrop Technologies). The cDNA was synthesized from 500 ng of total RNA using the PrimeScript RT Master Mix (TAKARA) and qPCR was run in the Step One Plus Real-time PCR thermocycler (Applied Biosystems) using the commercially available PowerUp SYBER Green Master Mix (Applied Biosystems).
Statistical analysis
Statistical analysis was performed using GraphPad Instat 3.00 for Windows (GraphPad Software). Three replicates for each experiment (growth curve, colony forming unit, cell cycle, differentiative potential and flow cytometric analysis) were performed. Results are reported as mean ± standard deviation (SD). One-way analysis of variance for multiple comparisons by the Student–Newman–Keuls multiple comparison test was used to assess differences between groups. Differences were considered statistically significant for p values < 0.05. For quantitative PCR data, nonparametric tests were used.
Results
Transcriptome analysis reveals DEGs between CV and CV-DS
To investigate transcriptomic differences between CV and DS-CV, we conducted a whole transcriptome NGS analysis (RNA-Seq) on three biological replicates for each MSC population. Quality control analyses were performed on raw Illumina reads using FastQC to confirm the quality of the sequencing data (Table 2). Large percentages of uniquely mapped reads that align to exactly one location within the reference genome (> 98.9%) confirmed the high quality of the library preparation [33]. Following current RNA-Seq guidelines, these parameters were deemed acceptable for pursuing further analysis of RNA-Seq data obtained from all samples (Additional file 1).
Plotting the Euclidean distances revealed fundamental differences and similarities between CV and DS-CV biological samples (Fig. 1a). All three CV biological replicates were found within optimal proximity to one another and significantly separated from their DS-CV counterparts. On the contrary, only two of the DS-CV samples grouped together, with the third replicate isolated from CV and DS-CV samples. To further investigate this phenomenon, we performed a hierarchical clustering heatmap including only genes expressed on chromosome 21. In this case, all three DS-CV samples clustered together (Fig. 1b). DS-CV samples on chromosome 2, 6 and 22 did not clustered, as shown in the correlation plots reported in Additional file 2: Fig. S1.
Differential analysis between CV and DS-CV samples identified a total of 1238 significantly altered transcripts (adjusted p-value [padj] < 0.05; |log2foldchange|> 1) (Fig. 1c, Additional file 3: Table S1). Specifically, 529 DEGs were up-regulated in DS-CV samples compared with CV, whereas 705 DEGs were down-regulated (Fig. 1c and Additional file 3: Table S1). Interestingly, only 30 up-regulated and 6 down-regulated transcripts originated from genes located in chromosome 21 (5.7% and 0.9% of the total DEGs, respectively). To gain further insight into which pathways were significantly enriched in DS-CV samples, we performed extensive in silico GSEA [24]. GSEA highlighted a significant down-regulation of processes related to cell cycle when testing the enrichment of GO biological processes [29, 34] and reactome pathways [35] (FDR < 0.05) (Fig. 1d and Additional file 4: Fig. S2). GSEA also pinpointed the presence of enriched cytobands (FDR < 0.1), including a down-regulated chromosome region 2q31.1 (Additional file 5: Table S2) and another such chromosome region 6p21.32 32 (Additional file 6: Table S3), along with an up-regulated chromosome region 22q11 (Additional file 7: Table S4) (Fig. 1d). The chromosome 2q31.1 band harbored a significant number of cell cycle-related genes including Cdca7 (Fold change [FC] = -3.9), Hoxd8 (FC = -5.4) and Klhl23 (FC = -4.2). In addition, the metabolism-related glutathione catabolic process (FDR = 3.8e−4) and cofactor metabolic process (FDR = 0.003) genes were enriched. On cytoband chromosome 6p21.32, we identified three significantly down-regulated members of the major histocompatibility complex (MHC-II) genes, Hla-dra (FC = -4.8), Hla-dpa1 (FC = 6.0) and Hla-dpb1 (FC = 6.3) in DS-CV samples. In addition, we also explored protein–protein interactions between significantly dysregulated genes on chromosome 21 and the chromosomes containing significantly enriched cytobands (2, 6 and 22) in our DS-CV samples, generated networks of genes that were dysregulated within our DS-CV samples. A significant portion of the detected DEGs were primarily linked to chromosomes 2, 22 and 6 (Fig. 1d and Additional file 8: Fig. S3). Chromosomes 2 and 22 including all their gene clusters involved in cell proliferation and cell cycle regulation showed a significantly altered activity that merited further investigation. Similar analysis has also been conducted for the gene clusters related to the immune response on chromosome 6 (Additional file 5: Table S2, Additional file 6: Table S3, Additional file 7: Table S4).
DEGs located on chromosome 2 and chromosome 22 influence the DS-CV proliferative properties
Figure 2a depicts a network diagram of the largest networks of interacting proteins dysregulated on chromosome 2 and chromosome 21 in DS-CV samples. These networks contain a total of 24 nodes and 33 edges, with 4 nodes up-regulated and 1 down-regulated on chromosome 21 and 2 up-regulated and 17 down-regulated on chromosome 2. The annotations of protein–protein interactions performed by the STRING cytoscape platform identified pathways involving these genes [35]. The five enriched pathways with the lowest FDR and corresponding genes in the network are shown in Fig. 2b. Genes identified in these pathways are located only on chromosome 22. Highlighted pathways were predominantly related to cell cycle, cell cycle checkpoints and key processes involved in the cell cycle, indicating that its function is altered in DS-CV cells when compared with CV cells. The largest network of interacting dysregulated genes on chromosome 21 and 22 is shown in Fig. 2c. The network contains a total of 6 nodes and 7 edges, with all nodes down-regulated on both chromosomes 21 and 22. The enriched pathways and corresponding genes are shown in Fig. 2d.
Functional studies show striking proliferative differences between CV and DS-CV
As shown in Fig. 3, significant differences (p < 0.05 and p < 0.01) were found between CV and DS-CV in their proliferative properties. The growth curve shows a similar lag phase (until day 3) between the two lines analyzed. Following the lag phase (days 4–13), we noticed a significant discrepancy in the growth rate of CV compared to DS-CV. At day 6 the number of cells produced by CV was almost twice that observed in DS-CV (4.8 × 104 ± 7.5 × 103 vs 2.8 × 104 ± 2.0 × 103, respectively). This trend was maintained until the end of the log phase when the cell growth plateaued. At the latest time point (day 15), a significant difference (p < 0.05) was still evident between the two groups, with 5.6 × 104 ± 1.2 × 104 and 2.8 × 104 ± 1.1 × 104 for CV and DS-CV, respectively (Fig. 3a). Cell cycle analysis reveals important growth characteristics in eukaryotic cells. Our data indicated a statistically significant increase in the percentage of cells in G0/G1 (+ 5%) and G2/M (+ 6%) phase in DS-CV compared to CV (Fig. 3b). No differences were observed for the S phase S (+ 1.6% in DS-CV compared to CV). Furthermore, while CV and DS-CV display a fibroblast-like morphology, slight shape differences became noticeable: DS-CV MSCs were more elongated compared to their euploid counterparts (Fig. 3c). Flow cytometry revealed that in both cell line preparations, over 90% of cells were positive for MSC-associated markers (CD90, CD29, CD44 and CD73) at protein level, a significant increase (p < 0.01) in the percentage of CD73 + DS-CV compared to CV (98.8 ± 0.5 vs 91 ± 1.1, respectively) (Fig. 3d).
CV and CV-DS show clonogenic and differentiative potential
Before proceeding with subsequent analyses, CV and CV-DS were also tested for clonogenicity and plasticity. The number of cell colonies formed was counted at P0 after seeding cells at different density/cm2. DS-CV display a clonogenic capacity similar to their euploid counterparts, with a direct correlation between the increase of CFU frequency and the increase of cell seeding density (Table 3).
Multipotent differentiation was evaluated toward the osteogenic and adipogenic lineages (Additional file 9: Fig. S4). Over two weeks of culture in osteogenic induction medium, CV and DS-CV distinctly changed their morphology and were surrounded by calcium deposits positive to von Kossa staining (Additional file 9: Fig. S4a). Expression of osteogenesis associated genes confirmed the occurred differentiation in both cell types with euploid MSCs showing a more marked specification process. Expression of the transcription factor Runx-2 was assessed around 7.08 ± 1.23 and 4.22 ± 0.32 for CV and DS-CV, while the expression of the osteogenic marker Bglap showed values that assessed around 9.93 ± 1.32 and 3.13 ± 0.98, respectively. In controls, the same changes were not observed. A similar trend was observed with regard to adipogenesis. Both cell lines were found to differentiate to adipocytes (Additional file 9: Fig. S4b), as shown by the presence of lipid vacuoles. Gene expression confirmed occurred differentiation and significant extents between the two cell lines. Expression of the transcription factor PPAR-γ was found to be significantly increased compared to control, with values assessed at 257.53 ± 8.76 and 79.32 ± 9.21 for CV and DS-CV, respectively, while the expression of the adipogenesis associated marker adiponectin was found to be 11.54 ± 2.47 and 7.40 ± 2.12 for the same samples.
DEGs located on chromosome 6 confer a marked potential to react to external stimuli
Figure 4a shows the largest network of DEGs located on chromosome 21 and chromosome 6 containing a total of 60 nodes and 110 edges, with 20 nodes up-regulated and 1 down-regulated on chromosome 21 and 20 up-regulated and 6 down-regulated on chromosome 6. Immune-related genes located on chromosomes 6 and 21 showed the most significant transcription dysregulations (Fig. 4b). Specifically, Ifnar1 (FC = 0.815, P adj = 0.005) and Ifnar2 (FC = 1.013, Padj = 0.004) (chromosome 21) were both up-regulated. An up-regulation of all components of the MHC (class I for Hla-b and class II for the other four), including Hla-b (FC = 2.836, Padj = 1.24e−6) and Hla-dma (FC = 1.980, Padj = 0.012) and down-regulation of Hla-dra (FC = -4.757, Padj = 0.028), Hla-dpa1 (FC = -6.012, Padj = 0.0104) and Hla-dpb1 (FC = -6.278, Padj = 0.0029) were all from chromosome 6. Significantly enriched genes and processes corresponding to chromosomes 21/6 network are also shown. The enrichment of genes located in the MHC protein complex (FDR = 3.11e−6) was further confirmed by enrichment analysis. The innate immune response (FDR = 3.4e−2) and cellular response to cytokine stimulus (FDR = 7.7e-3) pathways were also enriched, indicating that immune response function is altered in DS-CV cells compared with CV cells.
To confirm the dysregulation for the inflammatory genes found in the NGS data, we analyzed the cellular environment through the measure of intracellular ROS levels and expression of inflammatory genes. The measurement of intracellular ROS levels production serves as cell signaling molecules for normal biologic processes. DS-CV displayed a prompt and marked ROS production compared to their CV counterparts (Fig. 4c). Specifically, DS-CV showed significantly increased levels of ROS (670.4 ± 52.1 p < 0.01) compared not only to CV (502.5 ± 49.3) but also to CV treated with H2O2 (565 ± 17 p < 0.5). The exposure to an inflammatory environment, induced by the presence of tumor necrosis factor-α and interferon-γ (TNF-α and IFN-γ), resulted in a statistically significant increase in the expression of MSC-associated markers compared to their baseline levels. These changes were remarkable in DS-CV. Specifically, compared to CV after 48 h of treatment, the expression levels found in DS-CV were as follows: 1.73 ± 0.05-fold (p < 0.05) for Cd105, 3.61 ± 0.32-fold (p < 0.05) for Cd90, 5.07 ± 0.18-fold (p < 0.01) for Cd73, 5.85 ± 0.37-fold (p < 0.01) for Cd44 and 8.90 ± 0.18-fold (p < 0.01) for Cd29, respectively (Fig. 4d). MSCs from DS-CV were also more prone than CV to react to pro-inflammatory cytokines (TNF-α and IFN-γ) as demonstrated by the evaluation of the immunosuppressive genes tested. Twenty-four hours after the treatment started, the levels of mRNA of Pge2 (Prostaglandin E2) and Tgf-β (Transforming growth factor beta 1) were up-regulated when compared to CV: 3.21 ± 0.25-fold (p < 0.05) and 2.92 ± 0.009-fold (p < 0.05), respectively (Fig. 4e). Moreover, after 48 h of treatment, the expression of immunosuppressive genes was increased progressively from a median fold of 3.83 (± 0.31), 2.38(± 0.28), 3.21(± 0.31) increase for Pge2, Cox2 (cyclooxygenase 2), Tgf-β (p < 0.05), respectively. Expression levels of the pro-inflammatory molecule Tnf-α were also found up-regulated in DS-CV compared to CV with values assessed around 14 ± 0.5-fold (p < 0.01) and 103(± 0.06)-fold increase at 24 and 48 h, respectively.
Inflammatory genes increase on monocytic cells after DS-CV conditioned media
To confirm the immunosuppressive capability of CV and DS-CV and understand whether this process could be mediated by the paracrine signals they release, media conditioned by the presence of CV and DS-CV (CM CV and CM DS-CV, respectively) was added to the THP1 monocytes. The effect of CM CV and CM DS-CV was assessed by analyzing the expression of anti- and pro-inflammatory genes by THP-1 at 24 and 48 h. As can be observed in Fig. 5a, after 24 h of exposure to CM DS-CV, a statistically significant increase in the expression of anti- (Pge2, Tgf-α and Cox2) inflammatory genes was observed in THP-1 compared to the treatment with CM CV. A 16 ± 0.3-fold increase was observed for Pge2 (p < 0.01), a 7.4 ± 0.1-fold increase (p < 0.01) was found for Tgf-β and 2.2 ± 0.3-fold increase (p < 0.05) was assessed for Cox2. No statistical differences between the two groups were observed in the expression of pro-inflammatory genes, Interleukin -6 (Il-6) and Interleukin-1β (Il-1β). These values were found to drastically drop after 48 h of exposure in CM DS-CV although significant differences were still maintained compared to their CM CV counterparts for the expression of Il-1β (1.62 ± 0.1-fold p < 0.05), Pge2 (1.99 ± 0.4-fold p < 0.01) and Tgf-β (2.4 ± 0.3-fold p < 0.01) (Fig. 5b). On the other hand, an increase in the expression levels of Il-6 (2.22 ± 0.47-fold, p < 0.05) was noted in THP-1 treated with CM DS-CV, compared to those exposed to CM CV.
Discussion
In this study, we proposed that MSCs derived from a gestational source, the chorionic villi, isolated from fetuses diagnosed with DS, could be used as a means of performing deep, functional genomic-informed studies of the underlying clinical complexity of Down syndrome phenotypes. In particular, the availability of CV samples during prenatal diagnosis allows for the identification of patient-specific genetic alterations and opens the opportunity to develop therapeutic approaches that might optimize health outcomes in these individuals [12].
NGS analysis unraveled molecular differences between CV-DS and their euploid counterparts that may contribute to the pathogenesis of DS and influence the patient’s susceptibility to one of the many associated disorders (e.g., heart diseases, leukemia and neurological disorder such as epilepsy, and Alzheimer’s disease) that compromise the quality of life of DS patients [33]. These differences include the participation and the correlation of several chromosomes in addition to the chromosome 21, which play a role in the genetic imbalance associated with this condition. While early transcriptomic studies performed on primary culture of amniocytes [6, 36] have already revealed significant transcriptional alterations in the regions of other chromosomes, including 1, 6, 11, 16, 17 and 19 [16], the present study proposes additional chromosomes as responsible for the specific signatures associated with DS, namely chromosome 2 and 22, and highlights the role of new DEGs located on chromosome 6 as responsible for proliferation, inflammation and immune response in addition to chromosome 21 which to characterize to develop of additional comorbidities in Down syndrome. Our data showed the presence of DEGs in different chromosomes, such as the chromosome 2, 6 and 22, in line with observations from Letourneau et al. who suggested that the global gene expression in DS is based on differences in chromatin topology generating gene expression dysregulation domains (GEDDs) [37]. Other groups consider DEG heterogeneity as a stochastic mechanism occurring among individuals and within the same individual, which tightly depends on the cell type and tissue analyzed [12, 38]. Several factors contribute to such heterogeneity, including the trisomy-induced changes which alter cellular function and its interaction with surrounding cells, ultimately leading to a secondary distortion of function or gene expression in neighboring cells, and epigenetic factors [12]. The latter have been proposed to be involved in correct chromosome organization and loss of chromatin-accessibility that ultimately play a role in the neurodevelopmental pathogenesis of DS [39]. All these factors highlight the differences in DEGs obtained in different studies and they could also emphasize the important role of other chromosomes in the DS-associated genetic imbalance, besides the well-known extra copy of chromosome 21 and the recently suggested chromosome 6 [37].
Our findings show that most DEGs identified in CV-DS compared to CV are mainly involved in the regulation of the innate immune response and cell response to oxidative stress, which are located on chromosomes 6 and 21. In particular, NGS analysis showed an up-regulation of interferon receptors (IFNRs)-associated genes (Ifnr1 and Infr2) in chromosome 21 and Ifngr1 in chromosome 6 [37]. To date, IFN ligands are produced by cells in response to a variety of insults (such as IFN-α, β or γ) activating the signaling pathway for the expression of pro-inflammatory cytokines (e.g., IL-1β, TNF-α and IL-6) [38]. In vitro experiments confirmed these observations, as demonstrated by the marked increase of the levels of inflammatory genes in CV-DS compared to CV following the treatment with pro-inflammatory cytokines. This aspect is of a particular interest, as hyperactive IFN signaling has been shown to have profound negative impacts on human development, making it vulnerable at mostly viral or bacterial infections and thus determining DS vulnerability [11]. In addition, a study conducted by Sullivan et al. revealed a condition of chronic autoinflammation in DS with elevated levels of important inflammatory cytokines linked to IFN signaling, such as IL-6, TNF-α and monocyte chemoattractant protein-1 (MCP-1) [38], which might explain why DS individuals are more prone to develop additional comorbidities [40].
Further investigations have shed light on the involvement of transcriptional variations of App, Ifnr1 and Infr2 genes in brain neuropathologies and mental retardation in DS [41]. RNA-Seq performed by us revealed the up-regulation of another set of genes mapped on chromosome 6 that are involved in mitochondrial dysfunctions [42], the superoxide dismutase 2 genes (Sod2). It is known that Cu/ZnSOD activity is increased in DS patients, with the consequent production of more hydrogen peroxide than catalase and glutathione peroxidase can catabolize. This, in turn, gives rise to an oxidative stress positive feedback that induces mitochondrial dysfunction and impairs the respiratory complex enzymes, ultimately leading to an increased ROS production [43,44,45]. As reported in the literature, pathways involved with an excessive production of ROS by cells during aerobic respiration play important roles in the pathogenesis of a plethora of conditions, including cardiovascular diseases, hypertension and atherosclerosis [42]. The alteration of cellular redox status and the increase of the ROS generation has been associated to the up-regulation of β-amyloid precursor protein (APP) and to enhance DS individual’s vulnerability to various stress factors [46]. The App gene is located on chromosome 21 and its protein is expressed in glial cells and thus, in neuronal axons. APP protein is known to be localized in mitochondria and it is involved in regulation of secretory pathways. The overexpression of APP in neuronal cells is a common factor in DS patients and it is related to extra copy of chromosome 21 [47]. According to literature, we found an up-regulation in the APP-encoding gene present on chromosome 21. The disruption of normal signaling function of APP is reported to cause cell cycle abnormalities in neurons and to be associated with neurodegeneration and consequent dementia in DS. Indeed, approximately 50 to 70% of individuals with DS develop dementia by age 60 to 70 years [46]. Different studies report the neurological disease [46] as Alzheimer [48] is due to early amyloid plaques formation as consequent overexpression related to DS clinical conditions.
DEGs located on chromosome 6 also include Tnf and Hla, besides Sod2 and Ifngr1. The up-regulation of these genes is correlated with the induction of signal transducer and activator of transcription (STATs) and mitogen-activated protein kinase (MAPKs) pathways [49], which are associated with inflammatory conditions. The overexpression of inflammatory genes creates the chronic inflammation status within cells [50]. Interestingly, our transcriptomic profile showed a down-regulation of genes associated with MHC-II class genes and the Kinesin family member C1 (Kifc1) which are again located on chromosome 6. We found three Hla-genes (Hla-dra, Hla-dpa1 and Hla-dpb1) related with MHC which are involved in the immune suppression [49,50,51] and play a crucial role in tumor immune escape [52]. A low expression of HLA-class II molecules, as in our study with the DS sample, has been correlated with a higher risk of metastasis and associated to several tumors, such as sarcomas, hepatocellular carcinomas [53], B cell lymphomas, leukemia and autoimmune disorders [54]. Furthermore, the decreased susceptibility of the HLA-class II molecules has been linked to predisposition to develop type 1 diabetes a common pathology in DS [55]. Through GSEA analysis we found an up-regulation of several genes coding the histone cluster which make the nucleosome structure (Hist1-h1c, Hist1-h2ac, Hist1-h2bo, Hist1-h2bn, Hist1-h2bc and Hist1-h4h) located on chromosome 6. According to the literature, a correct nucleosome structure formation is required for the accurate assembly of chromatin as well as for the proper cell division and growth. Their overexpression has been shown to block transcription in vitro by triggering chromatin aggregation and to increase the incidence of mitotic chromosome loss, DNA damage, cell growth inhibition, inflammatory status promotion and cell toxicity [56]. The overexpression of histone gene Hist1-h1c along with the Sirt1 gene has also been correlated to the development of diabetes, which is a common condition secondary to DS [57]. As mentioned above, the overexpression of inflammatory genes can define a framework of increased production of cytokines and signaling activity for the immune system and increase the risk for DS individuals to develop an autoimmune disease [58]. Concomitantly, the GSEA analysis (Go process, KEGG and Reactome) showed the presence of DEGs associated with the cell cycle, on chromosomes 2, 6 and 22. Interestingly, an altered expression of centromere protein (Cenp) genes and cell cycle checkpoints (i.e., Cenpm, Mcm5, Gtse1 and Cdc45) was found. CENPs are proteins expressed in the formation of centromeric structure and key elements for the formation of a correct kinetochore, essential for mitosis and whose dysregulation induces the mitosis alteration causing the alteration of the whole cell cycle [59, 60].
Two more DEGs were identified in DS-CV, Bub1 and Cdk1. During mitosis, Bub1 is required to regulate chromosome segregation and maintain telomeric genome integrity [61]. It also plays an important role in the TGF-β signaling [62]. Cdk1 gene is involved in cell cycle progression and mitochondrial metabolism [63]. Our functional studies have confirmed these observations, showing a significant reduction of the proliferative potential and an impairment in cell cycle in DS-CV samples compared to their euploid counterparts. Experimental data are further corroborated by other gene sets we identified on chromosome 22. They include the Cdc45 gene, a component of the CMG (CDC45-MCMs-GINS) complex, whose expression is related to cycle progression. In case of loss of the expression levels, Cdc45 determines the inhibition of DNA replication and G1-phase arrest with consequently inhibition of the cell proliferation [64]. In addition to the molecular justification provided above, impairment in cell cycle could also be associated with the increasing expression of the CD73 we found in DS-CV. To date, CD73 is an ectoenzyme involved in several cell activities, including tissue homeostasis and support of epithelial cell transport [65]. Cd73 and Cenp genes could interfere with the normal cell cycle, confirming the alteration in the cell cycle observed in DS-CV [66]. Taken together, our data reveal potential gene patterns which could influence the occurrence of diseases related to DS. The pathways we found dysregulated in this study are involved in cell cycle and in proliferation activity, which have been linked with neuropathologies, including Alzheimer's disease, which represents one of the typical DS comorbidities [67]. Furthermore, the expression of CD73 strongly correlates with the responses to immunity and cancer, since this protein is involved in the reinforcement of lymphocyte–endothelium interactions, inhibition of macrophage and mesenchymal cell-mediated inflammation [68]. CD73 induction is associated with the increase of the expression of pro-inflammatory cytokines (i.e., IL-6 and IL-27), high levels of oxidative stress and hypoxia, among others [69]. With several molecular moieties being involved in the greater reactivity of DS-CV, we hypothesized that the DS-CV phenotype could be reflected in their capability to release bioactive factors acting as mediators of cell-to-cell communication. To test this hypothesis, we investigated the immune-suppressive and immune-modulatory properties of DS-CV in comparison with CV by exposing monocytic cells (THP-1 cells) to media conditioned by each cell type. Gene expression performed on THP-1 cells grown in the presence of DS-CV CM, showed a significant increase in the levels of pro (Il-6 and Il1β)- and anti (Pge2 Tgf-β and Cox2)-inflammatory genes at 24 and 48 h, suggesting a potentially different molecular composition in the paracrine factors they release [28]. Among them, extracellular vesicles (EVs) play a crucial role in cellular signaling as they are taken up by recipient cells, determining the activation or inhibition of several biological processes including cell proliferation, differentiation, apoptosis and the immune response [70]. MSC-derived EVs are extremely relevant as diagnostic tools and for the development of therapeutic cell-free drugs [71, 72]. The reaction to external stimuli leads to a feedback response via cellular or exosomal signaling consequently to damages or altered pathways. Gauthier et al. detected an enhanced endosomal release into extracellular space via exosomes as a protective effect under pathological conditions [73]. According to this study, our GSEA analysis found DEGs on chromosome 2 and 6 that are involved in the regulation of endoplasmic reticulum (ER) with a key role in protein folding. The inflammatory state observed during Down syndrome state could therefore influence the ER creating an imbalance between the protein-folding load and the capacity of the ER, causing unfolded or misfolded proteins [74]. Based on the tight correlation between the inflammatory response and the ER stress, we observed in DS-CV an excess of metabolic factors being released such as lipids, glucose and cytokines and neurotransmitters. In addition, ER controls cellular metabolism by regulating protein synthesis and secretion, as well as triglyceride and cholesterol biosynthesis. The ER stress could then trigger an excessive lipid accumulation as well as an increase of glucose levels associated with both obesity and diabetes [75]. Furthermore, protein aggregation resulted from ER stress has been associated with the Alzheimer and Parkinson’s disease and confirmed by the GSEA analysis we performed, identifying the App gene among those involved in the endoplasmic reticulum which is reported to play a role in the development of the Alzheimer’s disease, another typical secondary condition in DS [76].
Conclusions
This study provides a deeper investigation on the role of chromosomes 21, 6, 2 and 22 play in determining the complex pathological phenotype of DS individuals. The transcriptomic analysis performed on MSCs isolated from CV obtained for prenatal diagnosis of genetic disorders allowed us to expand our knowledge of this complex clinical condition. Our data demonstrate that while trisomy 21 does not interfere with standard stem cell features, it does promote an imbalance in the expression of genes that are associated with essential cellular processes, impairing cell cycle at a structural and functional standpoint, and dysregulating the immune response under specific stimuli, ultimately leading to chronic disorders commonly associated with DS. The newly identified DEGs and their locations represent potential new targets for the development of therapeutic strategies against DS morbidities. Furthermore, the insights obtained from the in vitro evaluation of the paracrine signals released by DS-CV compared to CV provide a solid basis for their use as diagnostic biomarkers as their production and content is expected to reflect divergences across genetically induced pathological conditions. While datasets generated here are promising they call for a subsequent investigation based on a larger number of patient-derived samples to provide stronger evidence and draw a definitive conclusion.
Availability of data and materials
Raw and processed datasets have been deposited in NCBI’s Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo) with accession reference GSE184450. Reviewer token: gtmjicaaxfetzup).
Abbreviations
- A/A:
-
Antimycotic–antibiotics
- ADIPQ:
-
Adiponectin
- APC:
-
Allophycocyanin
- APP:
-
β-Amyloid precursor protein
- BGLAP:
-
Bone gamma-carboxyglutamate protein
- BSA:
-
Bovine serum albumin
- CENP:
-
Centromere protein
- CFU-F:
-
Colony forming unit—fibroblastic like
- COX2:
-
Cyclooxygenase 2
- CV:
-
Chorionic villi
- DEGs:
-
Differentially expressed genes
- DS-CV:
-
Down syndrome chorionic villi
- DS:
-
Down syndrome
- ER:
-
Endoplasmic reticulum
- EV:
-
Extracellular vesicles
- FBS:
-
Fetal bovine serine
- FDR:
-
False discovery rates
- FITC:
-
Fluorescein isothiocyanate
- GEDD:
-
Gene expression dysregulation domains
- GEO:
-
Gene expression omnibus
- GSEA:
-
Gene set enrichment analysis
- H2-DCFDA:
-
2′,7′-Dichlorodihydrofluorecein
- HBSS:
-
Hank’s balanced salt solution
- HG-DMEM:
-
High-glucose Dulbecco’s modified eagle medium
- IL-1β:
-
Interleukin 1 beta
- IL-6:
-
Interleukin 6
- INF:
-
Interferon gamma
- iPSC:
-
Induced pluripotent stem cell
- KIFC1:
-
Kinesin family member C1
- MAPK:
-
Mitogen-activated protein kinase
- MCP-1:
-
Monocyte chemoattractant protein-1
- MHC-II:
-
Major histocompatibility complex
- MSC:
-
Mesenchymal stem cell
- NCBI:
-
National Center for Biotechnology Information
- NGS:
-
Next-generation sequencing
- Padj:
-
Adjusted p-value
- PBS:
-
Phosphate-buffered saline
- PE:
-
Phycoerythrin
- PFA:
-
Paraformaldehyde
- PGE2:
-
Prostaglandin E2
- PI:
-
Propidium iodide
- PPAR-γ:
-
Peroxisome proliferator-activated receptor gamma
- RIN:
-
Integrity number
- ROS:
-
Reactive of reactive oxygen species
- RPMI:
-
Roswell Park Memorial Institute-1640
- RUNX-2:
-
RUNX family transcription factor 2
- SD:
-
Standard deviation
- SOD2:
-
Superoxide dismutase 2
- STAT:
-
Signal transducer and activator of transcription
- TFG-β:
-
Transforming growth factor beta 1
- THP-1:
-
Acute monocytic leukemia
- TNF-α:
-
Tumor necrosis factor-alpha
References
Antonarakis SE, Lyle R, Dermitzakis ET, Reymond A, Deutsch S. Chromosome 21 and Down syndrome: from genomics to pathophysiology. Nat Rev Genet. 2004;5(10):725–38.
Antonarakis SE. 10 years of genomics, chromosome 21, and Down syndrome. Genomics. 1998;51(1):1–16.
Lejeune J, Gauthier M, Turpin R. Human chromosomes in tissue cultures. C R Hebd Seances Acad Sci. 1959;248(4):602–3.
Ross JA, Spector LG, Robison LL, Olshan AF. Epidemiology of leukemia in children with Down syndrome. Pediatr Blood Cancer. 2005;44(1):8–12.
Liu B, Filippi S, Roy A, Roberts I. Stem and progenitor cell dysfunction in human trisomies. EMBO Rep. 2015;16(1):44–62.
Silva CR, Biselli-Perico JM, Zampieri BL, Silva WA Jr, de Souza JE, Burger MC, et al. Differential expression of inflammation-related genes in children with Down syndrome. Mediators Inflamm. 2016;2016:6985903.
Antonarakis SE, Epstein CJ. The challenge of Down syndrome. Trends Mol Med. 2006;12(10):473–9.
Hartley D, Blumenthal T, Carrillo M, DiPaolo G, Esralew L, Gardiner K, et al. Down syndrome and Alzheimer’s disease: common pathways, common goals. Alzheimers Dement. 2015;11(6):700–9.
Granese B, Scala I, Spatuzza C, Valentino A, Coletta M, Vacca RA, et al. Validation of microarray data in human lymphoblasts shows a role of the ubiquitin-proteasome system and NF-kB in the pathogenesis of Down syndrome. BMC Med Genomics. 2013;5(6):24.
Potier MC, Rivals I, Mercier G, Ettwiller L, Moldrich RX, Laffaire J, et al. Transcriptional disruptions in Down syndrome: a case study in the Ts1Cje mouse cerebellum during post-natal development. J Neurochem. 2006;97(Suppl 1):104–9.
Sullivan KD, Lewis HC, Hill AA, Pandey A, Jackson LP, Cabral JM, et al. Trisomy 21 consistently activates the interferon response. Elife. 2016;29(5): e16220.
Pelleri MC, Cattani C, Vitale L, Antonaros F, Strippoli P, Locatelli C, et al. Integrated quantitative transcriptome maps of human trisomy 21 tissues and cells. Front Genet. 2018;9:125.
Roper RJ, Reeves RH. Understanding the basis for Down syndrome phenotypes. PLoS Genet. 2006;2(3): e50.
Antonaros F, Zenatelli R, Guerri G, Bertelli M, Locatelli C, Vione B, et al. The transcriptome profile of human trisomy 21 blood cells. Hum Genomics. 2021;15(1):25.
Kong XD, Liu N, Xu XJ. Bioinformatics analysis of biomarkers and transcriptional factor motifs in Down syndrome. Braz J Med Biol Res. 2014;47(10):834–41.
You SH, Lee YS, Chang YJ, Lin CY, Wang TH, Chang YL. Gene expression profiling of amniotic fluid mesenchymal stem cells of monozygotic twins discordant for trisomy 21. Gene. 2020;15(738): 144461.
Brigida AL, Siniscalco D. Induced pluripotent stem cells as a cellular model for studying Down syndrome. J Stem Cells Regen Med. 2016;12(2):54–60.
Parolini O, Alviano F, Bagnara GP, Bilic G, Buhring HJ, Evangelista M, et al. Concise review: isolation and characterization of cells from human term placenta: outcome of the first international Workshop on Placenta Derived Stem Cells. Stem Cells. 2008;26(2):300–11.
Lindholm D, Wootz H, Korhonen L. ER stress and neurodegenerative diseases. Cell Death Differ. 2006;13(3):385–92.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37(8):907–15.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Abraham A, Pedregosa F, Eickenberg M, Gervais P, Mueller A, Kossaifi J, et al. Machine learning for neuroimaging with scikit-learn. Front Neuroinform. 2014;8:14.
Walter W, Sanchez-Cabo F, Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics. 2015;31(17):2912–4.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019;47(W1):W199–205.
Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48(D1):D498–503.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44(D1):D457–62.
Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: network analysis and visualization of proteomics data. J Proteome Res. 2019;18(2):623–32.
The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. 2021;49(D1):D325–34.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999;27(1):29–34.
Dominici M, Le Blanc K, Mueller I, Slaper-Cortenbach I, Marini F, Krause D, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy. 2006;8(4):315–7.
Bull MJ. Down syndrome. N Engl J Med. 2020;382(24):2344–52.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
Gillespie M, Jassal B, Stephan R, Milacic M, Rothfels K, Senff-Ribeiro A, et al. The reactome pathway knowledgebase 2022. Nucleic Acids Res. 2022;50(D1):D687–92.
FitzPatrick DR, Ramsay J, McGill NI, Shade M, Carothers AD, Hastie ND. Transcriptome analysis of human autosomal trisomy. Hum Mol Genet. 2002;11(26):3249–56.
Letourneau A, Santoni FA, Bonilla X, Sailani MR, Gonzalez D, Kind J, et al. Domains of genome-wide gene expression dysregulation in Down’s syndrome. Nature. 2014;508(7496):345–50.
Ahlfors H, Anyanwu N, Pakanavicius E, Dinischiotu N, Lana-Elola E, Watson-Scales S, et al. Gene expression dysregulation domains are not a specific feature of Down syndrome. Nat Commun. 2019;10(1):2489.
Meharena HS, Marco A, Dileep V, Lockshin ER, Akatsu GY, Mullahoo J, et al. Down-syndrome-induced senescence disrupts the nuclear architecture of neural progenitors. Cell Stem Cell. 2022;29(1):116–30.
Sullivan KD, Evans D, Pandey A, Hraha TH, Smith KP, Markham N, et al. Trisomy 21 causes changes in the circulating proteome indicative of chronic autoinflammation. Sci Rep. 2017;7(1):14818.
Wilcock DM, Griffin WS. Down’s syndrome, neuroinflammation, and Alzheimer neuropathogenesis. J Neuroinflammation. 2013;16(10):84.
Desai SS. Down syndrome: a review of the literature. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 1997;84(3):279–85.
Rachidi M, Lopes C. Mental retardation in Down syndrome: from gene dosage imbalance to molecular and cellular mechanisms. Neurosci Res. 2007;59(4):349–69.
Fukai T, Ushio-Fukai M. Superoxide dismutases: role in redox signaling, vascular function, and diseases. Antioxid Redox Signal. 2011;15(6):1583–606.
Lee M, Hyun D, Jenner P, Halliwell B. Effect of overexpression of wild-type and mutant Cu/Zn-superoxide dismutases on oxidative damage and antioxidant defences: relevance to Down’s syndrome and familial amyotrophic lateral sclerosis. J Neurochem. 2001;76(4):957–65.
Roat E, Prada N, Ferraresi R, Giovenzana C, Nasi M, Troiano L, et al. Mitochondrial alterations and tendency to apoptosis in peripheral blood cells from children with Down syndrome. FEBS Lett. 2007;581(3):521–5.
Valenti D, Manente GA, Moro L, Marra E, Vacca RA. Deficit of complex I activity in human skin fibroblasts with chromosome 21 trisomy and overproduction of reactive oxygen species by mitochondria: involvement of the cAMP/PKA signalling pathway. Biochem J. 2011;435(3):679–88.
Cheon MS, Dierssen M, Kim SH, Lubec G. Protein expression of BACE1, BACE2 and APP in Down syndrome brains. Amino Acids. 2008;35(2):339–43.
Wilkins HM, Swerdlow RH. Amyloid precursor protein processing and bioenergetics. Brain Res Bull. 2017;133:71–9.
Jiang Y, Mullaney KA, Peterhoff CM, Che S, Schmidt SD, Boyer-Boiteau A, et al. Alzheimer’s-related endosome dysfunction in Down syndrome is Abeta-independent but requires APP and is reversed by BACE-1 inhibition. Proc Natl Acad Sci U S A. 2010;107(4):1630–5.
Kessal K, Liang H, Rabut G, Daull P, Garrigue JS, Docquier M, et al. Conjunctival inflammatory gene expression profiling in dry eye disease: correlations With HLA-DRA and HLA-DRB1. Front Immunol. 2018;9:2271.
Park SK, Dahmer MK, Quasney MW. MAPK and JAK-STAT signaling pathways are involved in the oxidative stress-induced decrease in expression of surfactant protein genes. Cell Physiol Biochem. 2012;30(2):334–46.
Campoli M, Ferrone S. HLA antigen changes in malignant cells: epigenetic mechanisms and biologic significance. Oncogene. 2008;27(45):5869–85.
Leite FA, Lira RC, Fedatto PF, Antonini SR, Martinelli CE Jr, de Castro M, et al. Low expression of HLA-DRA, HLA-DPA1, and HLA-DPB1 is associated with poor prognosis in pediatric adrenocortical tumors (ACT). Pediatr Blood Cancer. 2014;61(11):1940–8.
Meissner M, Konig V, Hrgovic I, Valesky E, Kaufmann R. Human leucocyte antigen class I and class II antigen expression in malignant fibrous histiocytoma, fibrosarcoma and dermatofibrosarcoma protuberans is significantly downregulated. J Eur Acad Dermatol Venereol. 2010;24(11):1326–32.
Rimsza LM, Roberts RA, Miller TP, Unger JM, LeBlanc M, Braziel RM, et al. Loss of MHC class II gene and protein expression in diffuse large B-cell lymphoma is related to decreased tumor immunosurveillance and poor patient survival regardless of other prognostic factors: a follow-up study from the Leukemia and Lymphoma Molecular Profiling Project. Blood. 2004;103(11):4251–8.
Aitken RJ, Mehers KL, Williams AJ, Brown J, Bingley PJ, Holl RW, et al. Early-onset, coexisting autoimmunity and decreased HLA-mediated susceptibility are the characteristics of diabetes in Down syndrome. Diabetes Care. 2013;36(5):1181–5.
Parssinen J, Alarmo EL, Khan S, Karhu R, Vihinen M, Kallioniemi A. Identification of differentially expressed genes after PPM1D silencing in breast cancer. Cancer Lett. 2008;259(1):61–70.
Wang W, Wang Q, Wan D, Sun Y, Wang L, Chen H, et al. Histone HIST1H1C/H1.2 regulates autophagy in the development of diabetic retinopathy. Autophagy. 2017;13(5):941–54.
Wang L, Yang B, Jiang H, Yu G, Feng M, Lu X, et al. The molecular mechanism study of insulin in promoting wound healing under high-glucose conditions. J Cell Biochem. 2019;120(9):16244–53.
Carroll CW, Silva MC, Godek KM, Jansen LE, Straight AF. Centromere assembly requires the direct recognition of CENP-A nucleosomes by CENP-N. Nat Cell Biol. 2009;11(7):896–902.
Hori T, Okada M, Maenaka K, Fukagawa T. CENP-O class proteins form a stable complex and are required for proper kinetochore function. Mol Biol Cell. 2008;19(3):843–54.
Logarinho E, Bousbaa H. Kinetochore-microtubule interactions “in check” by Bub1, Bub3 and BubR1: the dual task of attaching and signalling. Cell Cycle. 2008;7(12):1763–8.
Kim T, Gartner A. Bub1 kinase in the regulation of mitosis. Anim Cells Syst (Seoul). 2021;25(1):1–10.
Xie B, Wang S, Jiang N, Li JJ. Cyclin B1/CDK1-regulated mitochondrial bioenergetics in cell cycle progression and tumor resistance. Cancer Lett. 2019;28(443):56–66.
Fu Y, Lv Z, Kong D, Fan Y, Dong B. High abundance of CDC45 inhibits cell proliferation through elevation of HSPA6. Cell Prolif. 2022;55(7): e13257.
Rosner M, Kowalska A, Freilinger A, Prusa AR, Marton E, Hengstschlager M. Cell cycle and cell size regulation in Down syndrome cells. J Neural Transm Suppl. 2003;67:51–8.
Gao ZW, Dong K, Zhang HZ. The roles of CD73 in cancer. Biomed Res Int. 2014;2014: 460654.
Minor M, Alcedo KP, Battaglia RA, Snider NT. Cell type- and tissue-specific functions of ecto-5’-nucleotidase (CD73). Am J Physiol Cell Physiol. 2019;317(6):C1079–92.
Allard B, Longhi MS, Robson SC, Stagg J. The ectonucleotidases CD39 and CD73: novel checkpoint inhibitor targets. Immunol Rev. 2017;276(1):121–44.
Barile L, Vassalli G. Exosomes: therapy delivery tools and biomarkers of diseases. Pharmacol Ther. 2017;174:63–78.
Abusamra AJ, Zhong Z, Zheng X, Li M, Ichim TE, Chin JL, et al. Tumor exosomes expressing Fas ligand mediate CD8+ T-cell apoptosis. Blood Cells Mol Dis. 2005;35(2):169–73.
Yu S, Liu C, Su K, Wang J, Liu Y, Zhang L, et al. Tumor exosomes inhibit differentiation of bone marrow dendritic cells. J Immunol. 2007;178(11):6867–75.
Gauthier SA, Perez-Gonzalez R, Sharma A, Huang FK, Alldred MJ, Pawlik M, et al. Enhanced exosome secretion in Down syndrome brain—a protective mechanism to alleviate neuronal endosomal abnormalities. Acta Neuropathol Commun. 2017;5(1):65.
Zhang K, Kaufman RJ. From endoplasmic-reticulum stress to the inflammatory response. Nature. 2008;454(7203):455–62.
Hotamisligil GS. Inflammation and metabolic disorders. Nature. 2006;444(7121):860–7.
Funding
B.C. received salary support from the Sêr Cymru II program, funded by the European Commission through the Horizon 2020 Marie Skłodowska-Curie Actions (MSCA) COFUND scheme and the Welsh European Funding Office (WEFO) under the European Regional Development Fund (ERDF) and Houston Methodist Research Institute (USA). This funding body played no role in the design of the study and collection, analysis and interpretation of data and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
SV conducted, designed the experiments, interpreted data and wrote the manuscript; DWJ and MQ conducted the bioinformatic analysis; GM contributed to data interpretation and manuscript drafting; AP conducted experimental analyses; RHF assisted paper drafting and provided advice; RSC and LWF supervised the bioinformatic analysis; BC conceived and designed the project, supervised research activities and wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was conducted on cells isolated at the Cytogenetic Laboratory Children's Hospital Salesi (Ancona, Italy) in accordance with the Declaration of Helsinki. Samples were obtained for prenatal diagnosis upon informed consent signature from patients for the use of tissue for research purposes (based on d.l.gs. 196/2003—101/2018 and privacy protection 8/2016) and ethical approval by the Regional Institutional Review Board (Comitato Etico Regione Marche, Ref. UE 679/2016) on 07/29/2019).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Quantitative PCR between CV and DS-CV performed to evaluate the expression gene levels of collagen I at 7 and 14 days of osteogenic induction. Results show an increased expression in both CV and DS-CV. Data are showed as the mean (±SD) of three biological replicates.
Additional file 2: Figure S1.
Heatmaps of Euclidean distance between CV and DS-CV. a Heatmap of Euclidean distance scores of normalized genes counts on chromosome 2 for DS-CV versus CV. Dendrograms show hierarchical clustering results. b Heatmap of Euclidean distance scores of normalized genes counts on chromosome 6 for DS-CV versus CV. c Heatmap of Euclidean distance scores of normalized genes counts on chromosome 22 for DS-CV versus CV. Dendrograms show hierarchical clustering results.
Additional file 3: Table S1.
The list of differential expressed genes (DEGs) in CV and DS-CV samples.
Additional file 4: Figure S2.
Gene Set Enrichment Analysis (GSEA) for all genes (GO process, KEGG pathway and Reactome pathway). GSEA highlighted ranks genes by log fold change to identify up- and downregulated pathways and processes. In GO process (a), KEGG (b) and Reactome pathway (c) cell cycle processes or cell cycle pathways related are significantly down-regulated.
Additional file 5: Table S2.
Gene set enrichment analysis (GSEA) for CHr21_2_ (GO process, KEGG pathway, Reactome pathway).
Additional file 6: Table S3.
Gene set enrichment analysis (GSEA) for CHr21_6_ (GO process, KEGG pathway, Reactome pathway).
Additional file 7: Table S4.
Gene set enrichment analysis (GSEA) for CHr21_22_ (GO process, KEGG pathway, Reactome pathway).
Additional file 8: Figure S3.
Chromosome cytoband enrichment. The GSEA applied to all differentially expressed genes identify the enriched genomic locations (chromosome cytobands) from the msigdb hallmarks collections. a Three cytobands outside of chromosome 21 were found to be enriched with FDR < 0.1. The band chr22q11.25 was up-regulated and the bands, chr2q31.1 and chr6p21.32 were down-regulated. b Part of the detected DEGs mainly linked to chromosomes 2, 22 and 6.
Additional file 9: Figure S4.
CV and DS-CV differentiative potential. Representative images showing CV and DS-CV undergoing osteogenic (a) and adipogenic (b) differentiation after 7 and 14 days of induction. Osteogenesis was confirmed by von Kossa staining to highlight mineral deposition and adipogenesis by Oil-red-O positive cytoplasmic neutral lipids. Expression levels of osteogenic markers, (the transcription factor RUNX-2 and osteocalcin, Bglap) and (the nuclear receptor PPAR-g and adiponectin, Adipq) were determined by quantitative RT-PCR after 14 days of induction. Data were normalized to the reference gene (Gapdh) and represented as fold change compared with the expression untreated CV or DS-CV. Values are mean ± SD (n = 3). Asterisks depict highly significant (*p < 0.01 and **p < 0.001).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Vaiasicca, S., Melone, G., James, D.W. et al. Transcriptomic analysis of stem cells from chorionic villi uncovers the impact of chromosomes 2, 6 and 22 in the clinical manifestations of Down syndrome. Stem Cell Res Ther 14, 265 (2023). https://doi.org/10.1186/s13287-023-03503-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13287-023-03503-4