Molecular signatures of cardiac defects in down syndrome lymphoblastoid cell lines suggest altered ciliome and hedgehog pathways

Forty percent of people with Down syndrome exhibit heart defects, most often an atrioventricular septal defect (AVSD) and less frequently a ventricular septal defect (VSD) or atrial septal defect (ASD). Lymphoblastoid cell lines (LCLs) were established from lymphocytes of individuals with trisomy 21, the chromosomal abnormality causing Down syndrome. Gene expression profiles generated from DNA microarrays of LCLs from individuals without heart defects (CHD(-); n?=?22) were compared with those of LCLs from patients with cardiac malformations (CHD(+); n?=?21). After quantile normalization, principal component analysis revealed that AVSD carriers could be distinguished from a combined group of ASD or VSD (ASD+VSD) carriers. From 9,758 expressed genes, we identified 889 and 1,016 genes differentially expressed between CHD(-) and AVSD and CHD(-) and ASD+VSD, respectively, with only 119 genes in common. A specific chromosomal enrichment was found in each group of affected genes. Among the differentially expressed genes, more than 65% are expressed in human or mouse fetal heart tissues (GEO dataset). Additional LCLs from new groups of AVSD and ASD+VSD patients were analyzed by quantitative PCR; observed expression ratios were similar to microarray results. Analysis of GO categories revealed enrichment of genes from pathways regulating clathrin-mediated endocytosis in patients with AVSD and of genes involved in semaphorin-plexin-driven cardiogenesis and the formation of cytoplasmic microtubules in patients with ASD-VSD. A pathway-oriented search revealed enrichment in the ciliome for both groups and a specific enrichment in Hedgehog and Jak-stat pathways among ASD+VSD patients. These genes or related pathways are therefore potentially involved in normal cardiogenesis as well as in cardiac malformations observed in individuals with trisomy 21.     

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.     

Large-scale gene-centric analysis identifies novel variants for coronary artery disease

Coronary artery disease (CAD) has a significant genetic contribution that is incompletely characterized. To complement genome-wide association (GWA) studies, we conducted a large and systematic candidate gene study of CAD susceptibility, including analysis of many uncommon and functional variants. We examined 49,094 genetic variants in ∼2,100 genes of cardiovascular relevance, using a customised gene array in 15,596 CAD cases and 34,992 controls (11,202 cases and 30,733 controls of European descent; 4,394 cases and 4,259 controls of South Asian origin). We attempted to replicate putative novel associations in an additional 17,121 CAD cases and 40,473 controls. Potential mechanisms through which the novel variants could affect CAD risk were explored through association tests with vascular risk factors and gene expression. We confirmed associations of several previously known CAD susceptibility loci (eg, 9p21.3:p<10⁻³³; LPA:p<10⁻¹⁹; 1p13.3:p<10⁻¹⁷) as well as three recently discovered loci (COL4A1/COL4A2, ZC3HC1, CYP17A1:p<5×10⁻⁷). However, we found essentially null results for most previously suggested CAD candidate genes. In our replication study of 24 promising common variants, we identified novel associations of variants in or near LIPA, IL5, TRIB1, and ABCG5/ABCG8, with per-allele odds ratios for CAD risk with each of the novel variants ranging from 1.06–1.09. Associations with variants at LIPA, TRIB1, and ABCG5/ABCG8 were supported by gene expression data or effects on lipid levels. Apart from the previously reported variants in LPA, none of the other ∼4,500 low frequency and functional variants showed a strong effect. Associations in South Asians did not differ appreciably from those in Europeans, except for 9p21.3 (per-allele odds ratio: 1.14 versus 1.27 respectively; P for heterogeneity = 0.003). This large-scale gene-centric analysis has identified several novel genes for CAD that relate to diverse biochemical and cellular functions and clarified the literature with regard to many previously suggested genes.     

Transchromosomic cell model of Down syndrome shows aberrant migration, adhesion and proteome response to extracellular matrix

Background  

Down syndrome (DS), caused by trisomy of human chromosome 21 (HSA21), is the most common genetic birth defect. Congenital heart defects (CHD) are seen in 40% of DS children, and >50% of all atrioventricular canal defects in infancy are caused by trisomy 21, but the causative genes remain unknown.     

Identification of novel candidate genes in heterotaxy syndrome patients with congenital heart diseases by whole exome sequencing

Heterotaxy syndrome (HS) involves dysfunction of multiple systems resulting from abnormal left-right (LR) body patterning. Most HS patients present with complex congenital heart diseases (CHD), the disability and mortality of HS patients are extremely high. HS has great heterogeneity in phenotypes and genotypes, which have rendered gene discovery challenging. The aim of this study was to identify novel genes that underlie pathogenesis of HS patients with CHD. Whole exome sequencing was performed in 25 unrelated HS cases and 100 healthy controls; 19 nonsynonymous variants in 6 novel candidate genes (FLNA, ITGA1, PCNT, KIF7, GLI1, KMT2D) were identified. The functions of candidate genes were further analyzed in zebrafish model by CRISPR/Cas9 technique. Genome-editing was successfully introduced into the gene loci of flna, kmt2d and kif7, but the phenotypes were heterogenous. Disruption of each gene disturbed normal cardiac looping while kif7 knockout had a more prominent effect on liver budding and pitx2 expression. Our results revealed three potential HS pathogenic genes with probably different molecular mechanisms.     

Whole Exome Sequencing in Patients with the Cuticular Drusen Subtype of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is the leading cause of irreversible blindness in elderly people worldwide. Cuticular drusen (CD) is a clinical subtype of AMD, which typically displays an earlier age at onset, and has a strong genetic component. Genetic studies support a role for rare sequence variants in CD susceptibility, and rare sequence variants in the CFH gene have been identified in 8.8% of CD cases. To further explore the role of rare variants in CD, we performed whole exome sequencing (WES) in 14 affected members of six families and 12 sporadic cases with CD. We detected rare sequence variants in CFH and FBLN5, which previously were shown to harbor rare variants in patients with CD. In addition, we detected heterozygous rare sequence variants in several genes encoding components of the extracellular matrix (ECM), including FBLN1, FBLN3/EFEMP1, FBLN5, FBLN6/HMCN1, FBN2, and COL15A1. Two rare pathogenic variants were identified in the COL15A1 gene: one in a sporadic case and another was found to segregate in a family with six affected individuals with CD. In addition, two rare pathogenic variants were identified in the FGL1 gene in three unrelated CD cases. These findings suggest that alterations in the ECM and in the coagulation pathway may play a role in the pathogenesis of CD. The identified candidate genes require further analyses in larger cohorts to confirm their role in the CD subtype of AMD. No evidence was found of rare sequence variants in a single gene that segregate with CD in the six families, suggesting that the disease is genetically heterogeneous.     

Unusual genotypes in the COL6A1 gene in parents of children with trisomy 21 and major congenital heart defects

Collagen type VI is a candidate for a role in the pathogenesis of congenital heart defects (CHD) in Down's syndrome. Three restriction fragment length polymorphisms of the COL6A1 gene were used to determine COL6A1 genotypes in 50 families of affected children with trisomy 21 (29 with congenital heart defects and 21 without) and 37 unrelated volunteers. We found seven unusual genotypes in the parents of affected children with Down's syndrome, five being unique to the parents of children with trisomy 21 and CHD. There were no unusual genotypes associated with other chromosome 21 loci. No single COL6A1 genotype was associated with CHD. Thus, the unusual genotypes unique to parents of affected children suggest that genetic variation in the COL6A1 gene region contributes to the pathogenesis of CHD in Down's syndrome.     

Heterotrisomy, a significant contributing factor to ventricular septal defect associated with Down syndrome?

Down syndrome (DS; trisomy 21) is associated with a wide range of variable clinical features, one of the most common being congenital heart defects (CHD). We used molecular genetic techniques to study the inheritance of genes on chromosome 21 in children with DS and CHD. Polymorphic markers on the long arm of chromosome 21 were analysed in 99 families who had a child with DS. Of these, 60 children had a CHD and 39 children had no CHD. Heterotrisomy describes the inheritance of an allele from each of three different grandparents. In some cases heterotrisomy will involve the inheritance of three different alleles. Heterotrisomic regions were defined as those showing retention of non-disjoining parental heterozygosity at polymorphic loci in the non-disjoined chromosomes of children with DS. Using polymorphic non-coding markers, we identified a consistent 9.6-cM minimum region (D21S167-HMG14) of heterotrisomy in children with DS and ventricular septal defect (VSD). Comparing individuals with DS and VSD to all others with DS (those either with no CHD or with any other CHD combined) shows the individuals with DS and VSD to have significantly more non-reduction or heterotrisomy in this region (P=0.006, Fisher's exact test, two-tailed). We postulate that heterotrisomy for a gene or genes in this region is a contributing factor to the pathogenesis of VSD in trisomy 21 either through the presence of three different specific alleles or through the presence of specific combinations of alleles.     

Decreased DGCR8 Expression and miRNA Dysregulation in Individuals with 22q11.2 Deletion Syndrome

Deletion of the 1.5–3 Mb region of chromosome 22 at locus 11.2 gives rise to the chromosome 22q11.2 deletion syndrome (22q11DS), also known as DiGeorge and Velocardiofacial Syndromes. It is the most common micro-deletion disorder in humans and one of the most common multiple malformation syndromes. The syndrome is characterized by a broad phenotype, whose characterization has expanded considerably within the last decade and includes many associated findings such as craniofacial anomalies (40%), conotruncal defects of the heart (CHD; 70–80%), hypocalcemia (20–60%), and a range of neurocognitive anomalies with high risk of schizophrenia, all with a broad phenotypic variability. These phenotypic features are believed to be the result of a change in the copy number or dosage of the genes located in the deleted region. Despite this relatively clear genetic etiology, very little is known about which genes modulate phenotypic variations in humans or if they are due to combinatorial effects of reduced dosage of multiple genes acting in concert. Here, we report on decreased expression levels of genes within the deletion region of chromosome 22, including DGCR8, in peripheral leukocytes derived from individuals with 22q11DS compared to healthy controls. Furthermore, we found dysregulated miRNA expression in individuals with 22q11DS, including miR-150, miR-194 and miR-185. We postulate this to be related to DGCR8 haploinsufficiency as DGCR8 regulates miRNA biogenesis. Importantly we demonstrate that the level of some miRNAs correlates with brain measures, CHD and thyroid abnormalities, suggesting that the dysregulated miRNAs may contribute to these phenotypes and/or represent relevant blood biomarkers of the disease in individuals with 22q11DS.     

Rare Variants in MYD88, IRAK4 and IKBKG and Susceptibility to Invasive Pneumococcal Disease: A Population-Based Case-Control Study

Although rare variants within the Toll-like receptor signalling pathway genes have been found to underlie human primary immunodeficiencies associated with selective predisposition to invasive pneumococcal disease (IPD), the contribution of variants in these genes to IPD susceptibility at the population level remains unknown. Complete re-sequencing of IRAK4, MYD88 and IKBKG genes was undertaken in 164 IPD cases from the UK and 164 geographically-matched population-based controls. 233 single-nucleotide variants (SNVs) were identified, of which ten were in coding regions. Four rare coding variants were predicted to be deleterious, two variants in MYD88 and two in IRAK4. The predicted deleterious variants in MYD88 were observed as two heterozygote cases but not seen in controls. Frequencies of predicted deleterious IRAK4 SNVs were the same in cases and controls. Our findings suggest that rare, functional variants in MYD88, IRAK4 or IKBKG do not significantly contribute to IPD susceptibility in adults at the population level.     

Novel Mutations in the WNT1, TMEM38B,P4HB,and PLS3 Genes in Four Unrelated Chinese Families With Osteogenesis Imperfecta

Objective: Osteogenesis imperfecta (OI) is a group of heritable fragile bone diseases, and the majority are caused by pathogenic variants in the COL1A1 and COL1A2 genes. We sought to identify the genetic causes and phenotypes of OI in Chinese patients without COL1A1 or COL1A2 mutations.Methods: Twenty-three patients who were diagnosed with sporadic OI but did not carry COL1A1/2 mutations were recruited, and their genomic DNA was analyzed using targeted next-generation sequencing of rare OI-related genes. The resulting damaging mutations in the probands and their parents were verified using Sanger sequencing. Moreover, the efficacy of long-term bisphosphonate treatment was evaluated in proband 1.Results: Compound heterozygous variants in the WNT1 and TMEM38B genes were identified in proband 1 and proband 2, respectively. A heterozygous mutation in the P4HB gene was identified in proband 3, and a hemizygous mutation in PLS3 was identified in proband 4. The unaffected parents of the probands (except the father of proband 4) with mutations in the WNT1, TMEM38B, and PLS3 genes were heterozygous carriers of each of the variants, respectively. Notably, proband 3 had the characteristic exophthalmos, flat nasal bridge and flat, wide forehead. None of the patients presented with dentinogenesis imperfecta or hearing loss. Furthermore, bisphosphonates exerted beneficial effects on proband 1, who carried the WNT1 mutations, by increasing bone mineral density Z-score, reshaping the compressed vertebrae and decreasing the fracture risk.Conclusion: We identified novel mutations and expanded the spectrum of phenotypes and genotypes of the extremely rare disorder OI.Abbreviations: BMD = bone mineral density; MIM = Mendelian Inheritance in Man; OI = osteogenesis imperfecta; PDI = protein disulfide isomerase     

Impact of variants on type-2 diabetes risk genes identified through genomewide association studies in polycystic ovary syndrome: a case–control study

Polycystic ovary syndrome (PCOS) is a common endocrine disorder in females, and is associated with altered metabolic processes in particular insulin resistance and diabetes mellitus. PCOS shares with type-2 diabetes (T2D) a number of features, including beta cell dysfunction, impaired glucose tolerance and dyslipidaemia. Recently, genomewide association studies (GWAS) have reported a number of genes with reproducible associations and susceptibilities to T2D. To address this, we examined the association between the T2D GWAS candidate genes (CDKAL1, CDKN2B, COL8A1, HHEX, IGF2BP2, KCNJ1, KCNQ1 and SLC30A8) and PCOS in Saudi women. A case–control study, includes 162 cases and 162 controls was enrolled. Genotyping was carried out by the allelic discrimination method. Our results showed that the variants including rs792837 of COL8A1, rs61873498 of KCNQ1 and rs13266634 of SLC30A8 genes to be significantly more frequent in PCOS patients than in controls. Our results suggest that COL8A1, KCNQ1 and SLC30A8, which are previously identified through GWAS as T2D-associated genes, are associated with PCOS.     

Altered Expression of Immune-Related Genes in Children with Down Syndrome

Individuals with Down syndrome (DS) have a high incidence of immunological alterations with increased susceptibility to bacterial and viral infections and high frequency of different types of hematologic malignancies and autoimmune disorders. In the current study, we profiled the expression pattern of 92 immune-related genes in peripheral blood mononuclear cells (PBMCs) of two different groups, children with DS and control children, to identify differentially expressed genes that might be of pathogenetic importance for the development and phenotype of the immunological alterations observed in individuals with DS. PBMCs samples were obtained from six DS individuals with karyotypically confirmed full trisomy 21 and six healthy control individuals (ages 2–6 years). Gene expression was profiled in duplicate according to the manufacturer''s instructions provided by commercially available TaqMan Human Immune Array representing 92 immune function genes and four reference genes on a 96-plex gene card. A set of 17 differentially expressed genes, not located on chromosome 21 (HSA21), involved in immune and inflammatory pathways was identified including 13 genes (BCL2, CCL3, CCR7, CD19, CD28, CD40, CD40LG, CD80, EDN1, IKBKB, IL6, NOS2 and SKI) significantly down-regulated and four genes (BCL2L1, CCR2, CCR5 and IL10) significantly up-regulated in children with DS. These findings highlight a list of candidate genes for further investigation into the molecular mechanism underlying DS pathology and reinforce the secondary effects of the presence of a third copy of HSA21.     

Whole-exome sequencing identifies novel pathogenic mutations and putative phenotype-influencing variants in Polish limb-girdle muscular dystrophy patients

Background

Limb girdle muscular dystrophies (LGMD) are a group of heterogeneous hereditary myopathies with similar clinical symptoms. Disease onset and progression are highly variable, with an elusive genetic background, and around 50% cases lacking molecular diagnosis.

Methods

Whole exome sequencing (WES) was performed in 73 patients with clinically diagnosed LGMD. A filtering strategy aimed at identification of variants related to the disease included integrative analysis of WES data and human phenotype ontology (HPO) terms, analysis of genes expressed in muscle, analysis of the disease-associated interactome and copy number variants analysis.

Results

Genetic diagnosis was possible in 68.5% of cases. On average, 36.3 rare variants in genes associated with various muscle diseases per patient were found that could relate to the clinical phenotype. The putative causative mutations were mostly in LGMD-associated genes, but also in genes not included in the current LGMD classification (DMD, COL6A2, and COL6A3). In three patients, mutations in two genes were suggested as the joint cause of the disease (CAPN3+MYH7, COL6A3+CACNA1S, DYSF+MYH7). Moreover, a variety of phenotype-influencing variants were postulated, including in patients with an identified already known primary pathogenic mutation.

Conclusions

We hypothesize that LGMD could be better described as oligogenic disorders in which dominant clinical presentation can result from the combined effect of mutations in a set of genes. In this view, the inter- and intrafamilial variability could reflect a specific genetic background and the presence of sets of phenotype-influencing or co-causative mutations in genes that either interact with the known LGMD-associated genes or are a part of the same pathways or structures.

     

Identification of Novel Candidate Genes for Early-Onset Colorectal Cancer Susceptibility

Approximately 25–30% of colorectal cancer (CRC) cases are expected to result from a genetic predisposition, but in only 5–10% of these cases highly penetrant germline mutations are found. The remaining CRC heritability is still unexplained, and may be caused by a hitherto-undefined set of rare variants with a moderately penetrant risk. Here we aimed to identify novel risk factors for early-onset CRC using whole-exome sequencing, which was performed on a cohort of CRC individuals (n = 55) with a disease onset before 45 years of age. We searched for genes that were recurrently affected by rare variants (minor allele frequency ≤0.001) with potentially damaging effects and, subsequently, re-sequenced the candidate genes in a replication cohort of 174 early-onset or familial CRC individuals. Two functionally relevant genes with low frequency variants with potentially damaging effects, PTPN12 and LRP6, were found in at least three individuals. The protein tyrosine phosphatase PTP-PEST, encoded by PTPN12, is a regulator of cell motility and LRP6 is a component of the WNT-FZD-LRP5-LRP6 complex that triggers WNT signaling. All variants in LRP6 were identified in individuals with an extremely early-onset of the disease (≤30 years of age), and two of the three variants showed increased WNT signaling activity in vitro. In conclusion, we present PTPN12 and LRP6 as novel candidates contributing to the heterogeneous susceptibility to CRC.     

Analysis of pathogenic variants in 605 Chinese children with non-syndromic cardiac conotruncal defects based on targeted sequencing

ObjectiveDeleterious genetic variants comprise one cause of cardiac conotruncal defects (CTDs). Genes associated with CTDs are gradually being identified. In the present study, we aimed to explore the profile of genetic variants of CTD-associated genes in Chinese patients with non-syndromic CTDs.MethodsThirty-nine CTD-related genes were selected after reviewing published articles in NCBI, HGMD, OMIM, and HPO. In total, 605 patients with non-syndromic CTDs and 300 healthy controls, all of Han ethnicity, were recruited. High-throughput targeted sequencing was used to detect genetic variants in the protein-coding regions of genes. We performed rigorous variant-level filtrations to identify potentially damaging variants (Dvars) using prediction programs including CADD, SIFT, PolyPhen-2, and MutationTaster.ResultDvars were detected in 66.7% (26/39) of the targeted CTD-associated genes. In total, 11.07% (67/605) of patients with non-syndromic CTDs were found to carry one or more Dvars in targeted CTD-associated genes. Dvars in FOXH1, TBX2, NFATC1, FOXC2, and FOXC1 were common in the CTD cohort (1.5% [9/605], 1.2% [7/605], 1.2% [7/605], 1% [6/605], and 0.5% [3/605], respectively).ConclusionTargeted exon sequencing is a cost-effective approach for the genetic diagnosis of CTDs. Our findings contribute to an understanding of the genetic architecture of non-syndromic CTDs.     

Genetic associations with coronary heart disease: Meta-analyses of 12 candidate genetic variants

Aims

The aim of this study was to evaluate the combined contribution of 12 genetic variants to the risk of coronary heart disease (CHD).

Methods

Through a comprehensive literature search for genetic variants involved in the CHD association study, we harvested a total of 10 genes (12 variants) for the current meta-analyses. These genes consisted of GPX1 (rs1050450), PPARD (rs2016520), ALOX15 (rs34210653), SELPLG (rs2228315), FCGR2A (rs1801274), CCL5 (rs2107538), CYP1A1 (rs4646903), TP53 (rs1042522), CX37 (rs1764391), and PECAM1 (rs668, rs12953, and rs1131012).

Results

A total of 45 studies among 23,314 cases and 28,430 controls were retrieved for the meta-analyses of 12 genetic variants. The results showed a significant association between the GPX1 rs1050450 polymorphism and CHD (odd ratio (OR) = 1.61, 95% confidence interval (CI) = 1.25–2.07, P = 0.0002). Other meta-analyses of the rest 11 variants suggested a lack of association with the risk of CHD.

Conclusion

Our results confirmed that GPX1 rs1050450 was associated with susceptibility to CHD in Chinese and Indian populations.