1. Introduction
Down syndrome (DS) is characterized by a high variability of clinical pictures and symptomatology that can involve many organs and tissues, in particular, the circulatory system, the central nervous system and the immune system. Such comorbidity in one or more organs and/or systems has been demonstrated to reduce the life expectancy of patients, with the development of phenotypic traits similar to those observed during aging in chromosomally normal subjects [
1].
Consistent with these observations, the incidence of congenital heart defects (CHDs) is 0.8% in the general population and rises to approximately 40–60% in DS [
2,
3]. One of the largest studies in patients with DS confirmed this percentage, detecting cardiovascular abnormalities in 342 out of 821 (42%) of children with DS born between 1985 and 2006 in the north-eastern region of England [
4]. A similar incidence rate (43.9%) has been recently reported by Santoro and coworkers in 230 Italian patients with DS [
5]. Atrioventricular septal defects (AVSDs) are the most common forms of CHDs in DS, followed by ventricular septal defects (VSDs) and tetralogy of Fallot (TF) [
6,
7]. Moreover, some asymptomatic adolescents and adults without CHDs at birth may develop acquired valve anomalies and electric defects [
8]. Despite the high incidence of CHDs in DS, the pathophysiology and the pathways behind them remain unknown, as well as the close relationship with DS.Gene mapping analysis studies on human chromosome 21 (HSA21) have revealed that the supernumerary genes related to DS onset (such as superoxide dismutase Cu/Zn, transcription factors, stress inducible factors, cytokine receptors and amyloid precursor protein) that codify molecules witha key role in different metabolic pathways involved in diverse processes, such as angiogenesis, as well as maintaining cardiovascular homeostasis [
9]. Even though the genetic imbalance caused by the triplication of HSA21 can surely contribute to the impact on thepathways involved in heart development, it cannot be sufficient to determine CHD onset [
10]. Thishas led to advancing various hypotheses on the genotype–phenotype relationship of DS with CHDs. The gene dosage imbalance hypothesis states that the increased copy number of genes that are on HSA21 can lead to an overall increase in gene and protein expression, resulting in the phenotypic manifestations associated with DS [
11].
The amplified developmental instability hypothesis suggests that gene over-dosages on HSA21 result in a general disruption of genomic regulation and function of genes, inducing the expressionof the traits associated with DS [
12]. The most accredited hypothesis, known as the critical region hypothesis, affirms that a small set of genes within the Down syndrome critical or chromosomal region (DSCR) is responsible for the development of common DS phenotypes [
13,
14].
However, data from animal models, such as mouse models, have contributed to evidence that CHDs in DS are likely associated with a cluster of genes. Variants of these genes (i.e., mutations and polymorphisms) may influence onset, phenotypic and severity heterogeneity of CHDs in DS. Liu et al., using the Ts1Cje mouse model for DS, identified a 3.7 Mb critical genomic region associated with DS heart defects [
15,
16]. This region is located between the
Tiam1 and
IL10RB loci and contains a cluster of four interferon receptor (IFNR) genes:
IFNAR1,
IFNAR2,
IFNGR2 and
IL10RB.
IFNAR1 and
IFNAR2 are implied in the expression of type I interferon (IFN-α/β) receptors. The IFN-γR2 chain is a component of the type II interferon receptor complex (IFN-γ) and IL-10Rβ participates in type III interferon(IFN-λ) receptor structure [
17]. These data might lead to the hypothesis that mutation and common polymorphisms (single nucleotide polymorphisms, SNPs) of gene, such as IFN receptor gene SNP, mapping at HSA21 might be responsible for the heterogeneity of expression of the diverse organ/system anomalies characterizing the Down syndrome phenotypes, such as CHDs. Increased expression of these receptors in patients with DS has been demonstrated to contribute to autoimmune diseases [
18]. Impaired activation of IFN signaling has also been reported in individuals with DS and periodontitis [
19]. To our knowledge, at this moment, no data are available on the role of IFNreceptor cluster genetic assetsin heart defect pathogenesis. On the other hand, in animal models, the increased dosage of genes on the homolog of HSA21 explains only inpart the increased risk to CHDs and suggests the involvement of other genes located on other chromosomes [
20,
21,
22]. A study, conducted on strains of mice that were trisomic for segments of MMU16, homologous to HSA21, demonstratedan overexpression of vascular endothelial growth factor (VEGF) in severe DS heart abnormalities, including an attenuated compact layer of myocardium, the overproduction of trabeculae, defective ventricular septation and remodeling of the outflow tract [
20]. Data obtained in knockout VEGFreceptor mice have therefore showna functional and structural defect in cardiac valve development, thereby suggesting an involvement of these mediators in heart valve formation [
21]. Furthermore, a recent meta-analysis reported a significant association of the rs1570360 andrs3025039
VEGF-A SNPs with the risk of CHD [
23]. In addition, an excess of deleterious variants in VEGF-A pathway genes in DS has been associated with Atrial septum defects ASDs [
24].
Based on these observations, we analyzed the frequencies of IFNAR2 rs2229207, IFNGR2 rs2834213, IL10RB rs2834167 and VEGFA rs3025039 SNPs in a group of young adults withDS compared to a group of healthy subjects, matched for age and gender, in order to evaluate whether these polymorphisms can influence the susceptibility of patients with DS to either CHDs or acquired defects.
3. Results
As reported in
Table 3,
IFNAR2 rs2229207 genotypes showed similar frequencies in DS and controls and were not associated with susceptibility to cardiac pathologies in DS (
Table 4).
Conversely, the
IFNGR2 rs2834213 GGG genotype, tagged as GG genotype, was more represented inpatients with DS than in controls (
Table 5). In addition, the rs2834213 GG genotype was significantly more represented in patients with DS affected by congenital cardiac pathologiesthan in controls (
Table 6). However, the distribution of
IFNGR2 rs2834213 genotype frequencies of the controls but not patientswith DS respected the Hardy–Weinberg equilibrium. On the other hand, statistical analyses reported in
Table 6 demonstrate that the DS intra-group analyses showed that this genotype was more represented in patients with DS compared to patients without cardiac pathologies and appeared significantly over-represented in patients with DS with CHDs compared to those that developed cardiac damage after birth. In contrast, the
IFNGR2 rs2834213-positive genotypes (AAA, AAG and AGG, indicated as A/*) and in particular the heterozygous genotypes, tagged as AG, appear to be significantly under-represented in the group of patients with congenital lesions. These results suggested that this SNP might be directly or through interaction with other genes implied in the susceptibility to congenital heartdamage in patients with DS.
The
IL10RB rs2834167 genotype distribution of controlsbut not patientswith DS respected the Hardy–Weinberg equilibrium. The homozygous allele A genotype frequency of rs2834167 was lower in patients with DS than in the control group (
Table 7). In addition, as reported in
Table 8, the frequency of
IL10RB rs2834167 G-positive genotypes (GGG, GGA and GAA, tagged as G/*), in particular heterozygous genotypes, was significantly increased in patients with DS. Moreover, the intragroup analyses evidenced a significant reduction of homozygousgenotypesand a significant increase in heterozygous genotypes in patients with cardiac pathologies compared both to the controls and CHD-unaffected individuals with DS (
Table 8). No significant differences were observed by stratifying results according to congenital or acquired cardiac pathologies.
Table 9 reports
VEGFA rs3025039 genotype frequencies of subjects affected by DS, controls and DS subgroups classified according the presence of congenitaland/or acquired cardiac pathologies. As reported in
Table 10, analyses of
VEGFA rs3025039 genotype distributions demonstrated a significantly increased frequency of the CC genotype in the patient groupwith respect to control subjects. No significant differences were found comparing the two subgroups of patients with or without cardiac pathologies. On the other hand, statistical analyses demonstrated that the CT genotype frequency among Down syndrome patients with congenital cardiac pathologies is significantly lower when compared to the CT frequency of patients with heart defects developed after birth (AwoCC).
Comprehensively, our results indicate that susceptibility to CHDs in patients with DS is influenced by common genetic variation (SNPs) of genes coding for IFN receptor mapping on critical regions of HSA21 and by genes located on other chromosomes, such as VEGFA.
4. Discussion
As is well known, the detection of heart defects, in particular congenital defects (CHDs), is a common finding in patients with DS. [
26]. Data from animal models indicate that a cluster of IFN receptor genes are closely related to a critical region located on chromosome 21 [
9,
15,
16].
We selected
INFAR2 rs2229207T>C,
INFGR2 rs2834213>G and
IL10RB rs2834167A>G SNPs to evaluate their relationship with the presence of congenital or acquired heart defects. The
INFAR2 rs2229207T>C SNP changes phenylalanine to serine at amino acid position 8 (F8S), inducing an increased response to IFNα and β [
27]. The intronic (intron 2)
INFGR2 rs2834213>G SNP is not located near a splice site (5582 nucleotides downstream of a splice donor site and 877 nucleotides upstream of a splice acceptor site) but, on the other hand, is in linkage disequilibrium withthree SNPs contained in a 300 bp promoter segment, forming a haplotype associated with high transcriptional activity in vitro [
28].
Our results evidenced that
IFNAR2 rs2229207 was not associated with significant susceptibility to cardiac pathologies in DS. Analyses of
IFNGR2 rs2834213 genotype distribution in the DS group showed a significant reduction of heterozygous genotypes compared to the expected frequency according to the Hardy–Weinberg (HW) equilibrium. Moreover, the G homozygous genotype was present in at least 30% of patients with DS and particularly in those with cardiac defects (
Table 4). Stratifying genotype frequencies of patients with cardiac defects according to the presence of congenital or of acquired without previous congenital cardiac pathologies, the G homozygous genotype was significantly increased in the first group. Our data suggested that the presence of an over-dosage of the G allele may be implied in the susceptibility to heart congenital damage in patients with DS.
Sustained type I and II IFN responses were documented in different cell lines derived from DS tissues, as fibroblasts, Epstain-Barr Virus-transformed B cells, monocytes and T cells [
29,
30,
31,
32].However, as was recently reported, the IFN activation pathway appears not to be essential for heart development [
33,
34,
35]. On the other hand, the activation of theses pathways, as a consequence ofinflammatory responses to intrauterine infections, induces interferon type I and type II production [
36,
37,
38] and activation of the JAK/STAT pathway, which, via its target gene p21, might directly inhibit the proliferation and selection of Sca1-positive cardiac stem cells with a lack of differentiation of cardiomyocyte progenitors. In addition, human fetal cardiac stem cells exposed to high levels of IFN-γ acquire immune competences that are conserved in adults and in particular conditions might predispose to acquired heart damage [
38,
39]. For example, an adult fibroblast cell line from patients with DS exposed to IFN-γ has been demonstrated to slowin proliferation when compared to normal fibroblasts [
29]. In addition, as recently stated by Qin et al. [
40], interferon-γ inhibits the differentiation of stem cells from different types of tissue by inhibiting the activation of Notch signaling. Moreover, type I inflammation-mediated dysregulation of Notch4 signaling is involved in endothelial dysfunction [
41]. Intriguingly, NOTCH4 is expressed in the developing heart and has previously been identified as a major player in early artery and endothelial-to-mesenchymal transformation, which is critical for endocardial cushiondifferentiation [
42,
43]. So, the
IFNGR2 rs2834213 G homozygous genotype might tag a genetically determined hyper-responsiveness to type II IFN, which in patients with DS, might be implied in susceptibility to CHDs by influencing crucial regulatory pathways in heart organogenesis.
Analyses of
IL10RB rs2834167 genotype distribution in the Down syndrome patient group showed a significant reduction of homozygous A genotype compared to the expected HW equilibrium frequency. The G-positive (homozygous and heterozygous) genotypes were present in more than 60% of patients with DS and, particularly, in patients with cardiac defects (
Table 7 and
Table 8). Intragroup analyses demonstrated a significant reduction of the AA genotype and a significant increase in AG genotypes in patients with DS with cardiac pathologies with respect to patients without cardiac defects. In this subgroup of patients, the odds ratio of association with the G-positive genotypes was very high (OR: 9.714 C.I.: 3.605–26.17). The rs2834167A>G SNP that inducesa lysine to glutamic acid changein the
IL10RB amino acidic sequence (K47E) seems to be associated with higher mRNA and cell-surface expression of IL-10Rβ [
44]. Increased IL10RB expression have been found to be associated with cardiovascular diseases in older heart patients withstroke and in patients with DS affected by periodontitis [
45,
46,
47].
IL10 RB engagement is the crucial event in IL-10/IL-10 receptor interaction to obtain the active receptor complex [
48]. The ligand–receptor interaction activates via IL-10Rα Janus kinase-1 and tyrosine kinase-2, which induces the activation of tyrosine phosphorylation and latent transcription factor-3 (STAT3). Alternatively, IL-10 receptor β engagement, that down regulates STAT3 pathways, is involved in the activation pathway involving p-AKT and p-mTOR [
47]. Cross talk between STAT3, AKT and mTOR modulates IL10-mediated inhibition of autophagy, one of the crucial mechanisms involved in cell development control and the clearance of damaged mitochondria by mitophagy [
49]. mTOR is hyper-activated in DS, leading to mitochondrial dysfunction in intra- and extra-uterine life [
50]. The increased levels of mTOR are actually implied in inhibition of autophagy induction, autophagosome formation and mitophagy [
49]. Mitochondrial dysfunction hasrecently been strictly associatedwith aberrant ciliogenesis [
51,
52]. Cilia are emerging as an important hub for signaling processes crucial during development that are also needed for heart formation [
52,
53,
54,
55]. In this view, a genetically determined IL10RB hyper-dosage might play a central role in the mTOR-mediated dysregulation of mitochondrial homeostasis and organogenesis in DS.
On the other hand, it is reasonably likely that both IFNGR2 rs2834213 and IL10RB rs2834167 risk genotypes represent the tags of other genetic variants on HSA21 relevant for CHDs. Actually, it is extremely improbable that a single polymorphism might be the principal agent responsible for susceptibility to CHDs in DS. Although the three copies ofthe genes of chromosome 21 certainly contribute to this risk, the increased gene dosage itself is not sufficient to cause pathologies, such as atrioventricular septal defects, that are extremely frequent in these patients. Indeed, half or more of those with DS have a structurally normal heart. Clearly, other factors in addition to trisomy 21 are required to cause CHDs in children with DS. However, the functional effects of IFNGR2 rs2834213 and IL10RB rs2834167 SNPs might be directly involved in alteration of the signaling cascade that in the presence of other critical genetic variants might lead to CHD development.
A recent paper [
22] suggests that rare variants in the Notch pathway and ciliome genes are associated with AVSDs in DS. These results provide further support for the involvement of the Notch pathway. In this view, our results allow us to hypothesize that minor genetic variants of the interferon receptor gene cluster (i.e.,
IFNGR2 rs2834213 G homozygous genotype and
IL10RB rs2834167 G positive genotypes) on HSA21 might be considered as tags ofa complex genetic constellation involving over-dosage of genetic susceptibility variantsonHSA21 and the presence of other susceptibility alleles on other chromosomes, resultingin CHDsfrequently, but not exclusively, associatedwith DS.
In particular,
VEGF gene variability may be important for several angiogenesis-associated diseases, such as tumors or coronary disease [
56]. Hyper-activation of the VEGFA pathway and an excess of deleterious variants in VEGFA and VEGFA pathway genes have been found to be associated with ASDs and TF [
23,
24,
57]. Our major findings on
VEGFA rs3025039 indicate that thehomozygous C genotype frequency is not significantly different among patients with DS affected or not by cardiac pathologies. In addition, T allele frequency is significantly lower in patients with CHDs and DS. The minor allele (T) of the rs3025039 polymorphism located inthe 3′UTR down-strand the coding region seems to be associated with lower plasma VEGF concentrations both in normal and in diseased subjects [
58,
59,
60]. However, regardless of the effects on cytokine production, gene expression studies confirmed that homeostasis genes involved in VEGFA signaling are important for normal heart septation and they are altered in either DS or euploid patients with CHDs [
61]. Our findings might therefore suggest that a genetically determined reduction of
VEGFA production might be protective against CHDs in DS individuals.