Whole genome sequencing delineates regulatory, copy number, and cryptic splice variants in early onset cardiomyopathy

Cardiomyopathy (CMP) is a heritable disorder. Over 50% of cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole-genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine percent of cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen percent harbored high-risk regulatory variants in promoters and enhancers of CMP genes (odds ratio 2.25, p = 6.70 × 10−7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (odds ratio 6.7–58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early onset CMP.This project was supported by the Ted Rogers Centre for Heart Research (SM, JE), the Canadian Institutes of Health Research (PJT 175034) (SM, JE) and by the Canadian Institutes of Health Research (ENP 161429), under the frame of ERA PerMed (SM). SM holds the Heart and Stroke Foundation of Canada & Robert M Freedom Chair in Cardiovascular Science. SWS holds the GlaxoSmithKline Endowed Chair in Genome Sciences at the Hospital for Sick Children and the University of Toronto. PGM holds a Canada Research Chair Tier 2 in Non-coding Disease Mechanisms. PGM acknowledges the support of the Government of Canada’s New Frontiers in Research Fund (NFRF), [NFRFE-2018-01305]. EO holds the Bitove Family Professorship of Adult Congenital Heart Disease. MM holds a Ramon y Cajal grant from the Spanish Ministry of Science and Innovation (RYC-2017-22249). WO is supported by funding from Fundació La Marató (321/C/2019). JB is funded by a Frans Van de Werf fellowship for clinical cardiovascular research, and by a senior clinical investigator fellowship of the FWO Flanders. KM was a National Science Foundation Graduate Research Fellow under grant no. DGE1144152 during the majority of the project. CS is the recipient of a National Health and Medical Research Council (NHMRC) Practitioner Fellowship (1154992). JI is the recipient of an NHMRC Career Development Fellowship (1162929). RDB is the recipient of a New South Wales Health Cardiovascular Disease Senior Scientist Grant. PSD is supported by the DBT/Wellcome Trust- Indian Alliance. We acknowledge the Labatt Family Heart Centre Biobank at the Hospital for Sick Children for access to DNA samples, and The Centre for Applied Genomics at the Hospital for Sick Children for performing WGS. We thank Xiucheng Cui and Emanuela Pannia for performing the zebrafish experiments at the SickKids Zebrafish Genetics and Disease Models Core (CRISPR-Cas9 and gRNA syntheses, zebrafish embryo microinjections, gRNA PCR validation, qRT-PCR, cardiac imaging). This research was made possible through access to the data and findings generated by the 100,000 Genomes Project. The 100,000 Genomes Project is managed by Genomics England Limited (a wholly owned company of the Department of Health and Social Care). The 100,000 Genomes Project is funded by the National Institute for Health Research and NHS England. The Wellcome Trust, Cancer Research UK and the Medical Research Council have also funded research infrastructure. The 100,000 Genomes Project uses data provided by patients and collected by the National Health Service as part of their care and support. We thank members of the ICGC/PCAWG working groups for generating the variant calls used in our case-control burden analyses.Peer Reviewed"Article signat per 38 autors/es: Robert Lesurf, Abdelrahman Said, Oyediran Akinrinade, Jeroen Breckpot, Kathleen Delfosse, Ting Liu, Roderick Yao, Gabrielle Persad, Fintan McKenna, Ramil R. Noche, Winona Oliveros, Kaia Mattioli, Shreya Shah, Anastasia Miron, Qian Yang, Guoliang Meng, Michelle Chan Seng Yue, Wilson W. L. Sung, Bhooma Thiruvahindrapuram, Jane Lougheed, Erwin Oechslin, Tapas Mondal, Lynn Bergin, John Smythe, Shashank Jayappa, Vinay J. Rao, Jayaprakash Shenthar, Perundurai S. Dhandapany, Christopher Semsarian, Robert G. Weintraub, Richard D. Bagnall, Jodie Ingles, Genomics England Research Consortium, Marta Melé, Philipp G. Maass, James Ellis, Stephen W. Scherer & Seema Mital"Postprint (published version

Collaboratively charting the gene-to-phenotype network of human congenital heart defects

Background How to efficiently integrate the daily practice of molecular biologists, geneticists, and clinicians with the emerging computational strategies from systems biology is still much of an open question. Description We built on the recent advances in Wiki-based technologies to develop a collaborative knowledge base and gene prioritization portal aimed at mapping genes and genomic regions, and untangling their relations with corresponding human phenotypes, congenital heart defects (CHDs). This portal is not only an evolving community repository of current knowledge on the genetic basis of CHDs, but also a collaborative environment for the study of candidate genes potentially implicated in CHDs - in particular by integrating recent strategies for the statistical prioritization of candidate genes. It thus serves and connects the broad community that is facing CHDs, ranging from the pediatric cardiologist and clinical geneticist to the basic investigator of cardiogenesis. Conclusions This study describes the first specialized portal to collaboratively annotate and analyze gene-phenotype networks. Of broad interest to the biological community, we argue that such portals will play a significant role in systems biology studies of numerous complex biological processes. CHDWiki is accessible at http://www.esat.kuleuven.be/~bioiuser/chdwikistatus: publishe

Phenotype-specific effect of chromosome 1q21.1 rearrangements and GJA5 duplications in 2436 congenital heart disease patients and 6760 controls

Recurrent rearrangements of chromosome 1q21.1 that occur via non-allelic homologous recombination have been associated with variable phenotypes exhibiting incomplete penetrance, including congenital heart disease (CHD). However, the gene or genes within the ∼1 Mb critical region responsible for each of the associated phenotypes remains unknown. We examined the 1q21.1 locus in 948 patients with tetralogy of Fallot (TOF), 1488 patients with other forms of CHD and 6760 ethnically matched controls using single nucleotide polymorphism genotyping arrays (Illumina 660W and Affymetrix 6.0) and multiplex ligation-dependent probe amplification. We found that duplication of 1q21.1 was more common in cases of TOF than in controls [odds ratio (OR) 30.9, 95% confidence interval (CI) 8.9-107.6); P = 2.2 × 10−7], but deletion was not. In contrast, deletion of 1q21.1 was more common in cases of non-TOF CHD than in controls [OR 5.5 (95% CI 1.4-22.0); P = 0.04] while duplication was not. We also detected rare (n = 3) 100-200 kb duplications within the critical region of 1q21.1 in cases of TOF. These small duplications encompassed a single gene in common, GJA5, and were enriched in cases of TOF in comparison to controls [OR = 10.7 (95% CI 1.8-64.3), P = 0.01]. These findings show that duplication and deletion at chromosome 1q21.1 exhibit a degree of phenotypic specificity in CHD, and implicate GJA5 as the gene responsible for the CHD phenotypes observed with copy number imbalances at this locu

Phenotype-specific effect of chromosome 1q21.1 rearrangements and GJA5 duplications in 2436 congenital heart disease patients and 6760 controls

Recurrent rearrangements of chromosome 1q21.1 that occur via non-allelic homologous recombination have been associated with variable phenotypes exhibiting incomplete penetrance, including congenital heart disease (CHD). However, the gene or genes within the ∼1 Mb critical region responsible for each of the associated phenotypes remains unknown. We examined the 1q21.1 locus in 948 patients with tetralogy of Fallot (TOF), 1488 patients with other forms of CHD and 6760 ethnically matched controls using single nucleotide polymorphism genotyping arrays (Illumina 660W and Affymetrix 6.0) and multiplex ligation-dependent probe amplification. We found that duplication of 1q21.1 was more common in cases of TOF than in controls [odds ratio (OR) 30.9, 95% confidence interval (CI) 8.9–107.6); P = 2.2 × 10−7], but deletion was not. In contrast, deletion of 1q21.1 was more common in cases of non-TOF CHD than in controls [OR 5.5 (95% CI 1.4–22.0); P = 0.04] while duplication was not. We also detected rare (n = 3) 100–200 kb duplications within the critical region of 1q21.1 in cases of TOF. These small duplications encompassed a single gene in common, GJA5, and were enriched in cases of TOF in comparison to controls [OR = 10.7 (95% CI 1.8–64.3), P = 0.01]. These findings show that duplication and deletion at chromosome 1q21.1 exhibit a degree of phenotypic specificity in CHD, and implicate GJA5 as the gene responsible for the CHD phenotypes observed with copy number imbalances at this locus

Overt Cleft Palate Phenotype and TBX1 Genotype Correlations in Velo-cardio-facial/DiGeorge/22q11.2 Deletion Syndrome Patients

Velo-cardio-facial syndrome/DiGeorge syndrome, also known as 22q11.2 deletion syndrome (22q11DS) is the most common microdeletion syndrome, with an estimated incidence of 1/2,000 – 1/4,000 live births. Approximately 9–11% of patients with this disorder have an overt cleft palate (CP), but the genetic factors responsible for CP in the 22q11DS subset are unknown. The TBX1 gene, a member of the T-box transcription factor gene family, lies within the 22q11.2 region that is hemizygous in patients with 22q11DS. Inactivation of one allele of Tbx1 in the mouse does not result in CP, but inactivation of both alleles does. Based on these data, we hypothesized that DNA variants in the remaining allele of TBX1 may confer risk to CP in patients with 22q11DS. To test the hypothesis, we evaluated TBX1 exon sequencing (n = 360) and genotyping data (n = 737) with respect to presence (n = 54) or absence (n = 683) of CP in patients with 22q11DS. Two upstream SNPs (rs4819835 and rs5748410) showed individual evidence for association but they were not significant after correction for multiple testing. Associations were not identified between DNA variants and haplotypes in 22q11DS patients with CP. Overall, this study indicates that common DNA variants in TBX1 may be nominally causative for CP in patients with 22q11DS. This raises the possibility that genes elsewhere on the remaining allele of 22q11.2 or in the genome could be relevant

Genotype and Cardiovascular Phenotype Correlations With TBX1 in 1,022 Velo-Cardio-Facial/Digeorge/22q11.2 Deletion Syndrome Patients

Haploinsufficiency of TBX1, encoding a T-box transcription factor, is largely responsible for the physical malformations in velo-cardio-facial /DiGeorge/22q11.2 deletion syndrome (22q11DS) patients. Cardiovascular malformations in these patients are highly variable, raising the question as to whether DNA variations in the TBX1 locus on the remaining allele of 22q11.2 could be responsible. To test this, a large sample size is needed. The TBX1 gene was sequenced in 360 consecutive 22q11DS patients. Rare and common variations were identified. We did not detect enrichment in rare SNP (single nucleotide polymorphism) number in those with or without a congenital heart defect. One exception was that there was increased number of very rare SNPs between those with normal heart anatomy compared to those with right-sided aortic arch or persistent truncus arteriosus, suggesting potentially protective roles in the SNPs for these phenotype-enrichment groups. Nine common SNPs (minor allele frequency, MAF \u3e 0.05) were chosen and used to genotype the entire cohort of 1,022 22q11DS subjects. We did not find a correlation between common SNPs or haplotypes and cardiovascular phenotype. This work demonstrates that common DNA variations in TBX1 do not explain variable cardiovascular expression in 22q11DS patients, implicating existence of modifiers in other genes on 22q11.2 or elsewhere in the genome

Genome-wide association study identifies loci on 12q24 and 13q32 associated with Tetralogy of Fallot

We conducted a genome-wide association study to search for risk alleles associated with Tetralogy of Fallot (TOF), using a northern European discovery set of 835 cases and 5159 controls. A region on chromosome 12q24 was associated (P = 1.4 × 10−7) and replicated convincingly (P = 3.9 × 10−5) in 798 cases and 2931 controls [per allele odds ratio (OR) = 1.27 in replication cohort, P = 7.7 × 10−11 in combined populations]. Single nucleotide polymorphisms in the glypican 5 gene on chromosome 13q32 were also associated (P = 1.7 × 10−7) and replicated convincingly (P = 1.2 × 10−5) in 789 cases and 2927 controls (per allele OR = 1.31 in replication cohort, P = 3.03 × 10−11 in combined populations). Four additional regions on chromosomes 10, 15 and 16 showed suggestive association accompanied by nominal replication. This study, the first genome-wide association study of a congenital heart malformation phenotype, provides evidence that common genetic variation influences the risk of TO

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

RATIONALE: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. OBJECTIVE: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. METHODS AND RESULTS: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10 CONCLUSIONS: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 nea