Effect of deletion of the protein kinase PRKD1 on development of the mouse embryonic heart.

Congenital heart disease (CHD) is the most common congenital anomaly, with an overall incidence of approximately 1% in the United Kingdom. Exome sequencing in large CHD cohorts has been performed to provide insights into the genetic aetiology of CHD. This includes a study of 1891 probands by our group in collaboration with others, which identified three novel genes-CDK13, PRKD1, and CHD4, in patients with syndromic CHD. PRKD1 encodes a serine/threonine protein kinase, which is important in a variety of fundamental cellular functions. Individuals with a heterozygous mutation in PRKD1 may have facial dysmorphism, ectodermal dysplasia and may have CHDs such as pulmonary stenosis, atrioventricular septal defects, coarctation of the aorta and bicuspid aortic valve. To obtain a greater appreciation for the role that this essential protein kinase plays in cardiogenesis and CHD, we have analysed a Prkd1 transgenic mouse model (Prkd1em1) carrying deletion of exon 2, causing loss of function. High-resolution episcopic microscopy affords detailed morphological 3D analysis of the developing heart and provides evidence for an essential role of Prkd1 in both normal cardiac development and CHD. We show that homozygous deletion of Prkd1 is associated with complex forms of CHD such as atrioventricular septal defects, and bicuspid aortic and pulmonary valves, and is lethal. Even in heterozygotes, cardiac differences occur. However, given that 97% of Prkd1 heterozygous mice display normal heart development, it is likely that one normal allele is sufficient, with the defects seen most likely to represent sporadic events. Moreover, mRNA and protein expression levels were investigated by RT-qPCR and western immunoblotting, respectively. A significant reduction in Prkd1 mRNA levels was seen in homozygotes, but not heterozygotes, compared to WT littermates. While a trend towards lower PRKD1 protein expression was seen in the heterozygotes, the difference was only significant in the homozygotes. There was no compensation by the related Prkd2 and Prkd3 at transcript level, as evidenced by RT-qPCR. Overall, we demonstrate a vital role of Prkd1 in heart development and the aetiology of CHD.

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.

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

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Whole Exome Sequencing Reveals the Major Genetic Contributors to Nonsyndromic Tetralogy of Fallot.

RATIONALE: Familial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, nonsyndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease phenotype. Rare genetic variants have been identified as important contributors to the risk of congenital heart disease, but relatively small numbers of TOF cases have been studied to date. OBJECTIVE: We used whole exome sequencing to assess the prevalence of unique, deleterious variants in the largest cohort of nonsyndromic TOF patients reported to date. METHODS AND RESULTS: Eight hundred twenty-nine TOF patients underwent whole exome sequencing. The presence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database and a scaled combined annotation-dependent depletion score of ≥20. The clustering of variants in 2 genes, NOTCH1 and FLT4, surpassed thresholds for genome-wide significance (assigned as P<5×10-8) after correction for multiple comparisons. NOTCH1 was most frequently found to harbor unique, deleterious variants. Thirty-one changes were observed in 37 probands (4.5%; 95% CI, 3.2%-6.1%) and included 7 loss-of-function variants 22 missense variants and 2 in-frame indels. Sanger sequencing of the unaffected parents of 7 cases identified 5 de novo variants. Three NOTCH1 variants (p.G200R, p.C607Y, and p.N1875S) were subjected to functional evaluation, and 2 showed a reduction in Jagged1-induced NOTCH signaling. FLT4 variants were found in 2.4% (95% CI, 1.6%-3.8%) of TOF patients, with 21 patients harboring 22 unique, deleterious variants. The variants identified were distinct to those that cause the congenital lymphoedema syndrome Milroy disease. In addition to NOTCH1, FLT4 and the well-established TOF gene, TBX1, we identified potential association with variants in several other candidates, including RYR1, ZFPM1, CAMTA2, DLX6, and PCM1. CONCLUSIONS: The NOTCH1 locus is the most frequent site of genetic variants predisposing to nonsyndromic TOF, followed by FLT4. Together, variants in these genes are found in almost 7% of TOF patients.

Acetylation of TBX5 by KAT2B and KAT2A regulates heart and limb development.

TBX5 plays a critical role in heart and forelimb development. Mutations in TBX5 cause Holt-Oram syndrome, an autosomal dominant condition that affects the formation of the heart and upper-limb. Several studies have provided significant insight into the role of TBX5 in cardiogenesis; however, how TBX5 activity is regulated by other factors is still unknown. Here we report that histone acetyltransferases KAT2A and KAT2B associate with TBX5 and acetylate it at Lys339. Acetylation potentiates its transcriptional activity and is required for nuclear retention. Morpholino-mediated knockdown of kat2a and kat2b transcripts in zebrafish severely perturb heart and limb development, mirroring the tbx5a knockdown phenotype. The phenotypes found in MO-injected embryos were also observed when we introduced mutations in the kat2a or kat2b genes using the CRISPR-Cas system. These studies highlight the importance of KAT2A and KAT2B modulation of TBX5 and their impact on heart and limb development.

A genome-wide association study of congenital cardiovascular left-sided lesions shows association with a locus on chromosome 20.

Congenital heart defects involving left-sided lesions (LSLs) are relatively common birth defects with substantial morbidity and mortality. Previous studies have suggested a high heritability with a complex genetic architecture, such that only a few LSL loci have been identified. We performed a genome-wide case-control association study to address the role of common variants using a discovery cohort of 778 cases and 2756 controls. We identified a genome-wide significant association mapping to a 200 kb region on chromosome 20q11 [P= 1.72 × 10-8 for rs3746446; imputed Single Nucleotide Polymorphism (SNP) rs6088703 P= 3.01 × 10-9, odds ratio (OR)= 1.6 for both]. This result was supported by transmission disequilibrium analyses using a subset of 541 case families (lowest P in region= 4.51 × 10-5, OR= 1.5). Replication in a cohort of 367 LSL cases and 5159 controls showed nominal association (P= 0.03 for rs3746446) resulting in P= 9.49 × 10-9 for rs3746446 upon meta-analysis of the combined cohorts. In addition, a group of seven SNPs on chromosome 1q21.3 met threshold for suggestive association (lowest P= 9.35 × 10-7 for rs12045807). Both regions include genes involved in cardiac development-MYH7B/miR499A on chromosome 20 and CTSK, CTSS and ARNT on chromosome 1. Genome-wide heritability analysis using case-control genotyped SNPs suggested that the mean heritability of LSLs attributable to common variants is moderately high ([Formula: see text] range= 0.26-0.34) and consistent with previous assertions. These results provide evidence for the role of common variation in LSLs, proffer new genes as potential biological candidates, and give further insight to the complex genetic architecture of congenital heart disease.

Tropomyosin 1: Multiple roles in the developing heart and in the formation of congenital heart defects.

Tropomyosin 1 (TPM1) is an essential sarcomeric component, stabilising the thin filament and facilitating actin's interaction with myosin. A number of sarcomeric proteins, such as alpha myosin heavy chain, play crucial roles in cardiac development. Mutations in these genes have been linked to congenital heart defects (CHDs), occurring in approximately 1 in 145 live births. To date, TPM1 has not been associated with isolated CHDs. Analysis of 380 CHD cases revealed three novel mutations in the TPM1 gene; IVS1+2T>C, I130V, S229F and a polyadenylation signal site variant GATAAA/AATAAA. Analysis of IVS1+2T>C revealed aberrant pre-mRNA splicing. In addition, abnormal structural properties were found in hearts transfected with TPM1 carrying I130V and S229F mutations. Phenotypic analysis of TPM1 morpholino-treated embryos revealed roles for TPM1 in cardiac looping, atrial septation and ventricular trabeculae formation and increased apoptosis was seen within the heart. In addition, sarcomere assembly was affected and altered action potentials were exhibited. This study demonstrated that sarcomeric TPM1 plays vital roles in cardiogenesis and is a suitable candidate gene for screening individuals with isolated CHDs.

Rare variants in NR2F2 cause congenital heart defects in humans.

Congenital heart defects (CHDs) are the most common birth defect worldwide and are a leading cause of neonatal mortality. Nonsyndromic atrioventricular septal defects (AVSDs) are an important subtype of CHDs for which the genetic architecture is poorly understood. We performed exome sequencing in 13 parent-offspring trios and 112 unrelated individuals with nonsyndromic AVSDs and identified five rare missense variants (two of which arose de novo) in the highly conserved gene NR2F2, a very significant enrichment (p = 7.7 × 10(-7)) compared to 5,194 control subjects. We identified three additional CHD-affected families with other variants in NR2F2 including a de novo balanced chromosomal translocation, a de novo substitution disrupting a splice donor site, and a 3 bp duplication that cosegregated in a multiplex family. NR2F2 encodes a pleiotropic developmental transcription factor, and decreased dosage of NR2F2 in mice has been shown to result in abnormal development of atrioventricular septa. Via luciferase assays, we showed that all six coding sequence variants observed in individuals significantly alter the activity of NR2F2 on target promoters.

High-content screening identifies small molecules that remove nuclear foci, affect MBNL distribution and CELF1 protein levels via a PKC-independent pathway in myotonic dystrophy cell lines.

Myotonic dystrophy (DM) is a multi-system neuromuscular disorder for which there is no treatment. We have developed a medium throughput phenotypic assay, based on the identification of nuclear foci in DM patient cell lines using in situ hybridization and high-content imaging to screen for potentially useful therapeutic compounds. A series of further assays based on molecular features of DM have also been employed. Two compounds that reduce and/or remove nuclear foci have been identified, Ro 31-8220 and chromomycin A3. Ro 31-8220 is a PKC inhibitor, previously shown to affect the hyperphosphorylation of CELF1 and ameliorate the cardiac phenotype in a DM1 mouse model. We show that the same compound eliminates nuclear foci, reduces MBNL1 protein in the nucleus, affects ATP2A1 alternative splicing and reduces steady-state levels of CELF1 protein. We demonstrate that this effect is independent of PKC activity and conclude that this compound may be acting on alternative kinase targets within DM pathophysiology. Understanding the activity profile for this compound is key for the development of targeted therapeutics in the treatment of DM.

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 TOF.

Contribution of global rare copy-number variants to the risk of sporadic congenital heart disease.

Previous studies have shown that copy-number variants (CNVs) contribute to the risk of complex developmental phenotypes. However, the contribution of global CNV burden to the risk of sporadic congenital heart disease (CHD) remains incompletely defined. We generated genome-wide CNV data by using Illumina 660W-Quad SNP arrays in 2,256 individuals with CHD, 283 trio CHD-affected families, and 1,538 controls. We found association of rare genic deletions with CHD risk (odds ratio [OR] = 1.8, p = 0.0008). Rare deletions in study participants with CHD had higher gene content (p = 0.001) with higher haploinsufficiency scores (p = 0.03) than they did in controls, and they were enriched with Wnt-signaling genes (p = 1 × 10(-5)). Recurrent 15q11.2 deletions were associated with CHD risk (OR = 8.2, p = 0.02). Rare de novo CNVs were observed in ~5% of CHD trios; 10 out of 11 occurred on the paternally transmitted chromosome (p = 0.01). Some of the rare de novo CNVs spanned genes known to be involved in heart development (e.g., HAND2 and GJA5). Rare genic deletions contribute ~4% of the population-attributable risk of sporadic CHD. Second to previously described CNVs at 1q21.1, deletions at 15q11.2 and those implicating Wnt signaling are the most significant contributors to the risk of sporadic CHD. Rare de novo CNVs identified in CHD trios exhibit paternal origin bias.

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.

Mutation in myosin heavy chain 6 causes atrial septal defect.

Atrial septal defect is one of the most common forms of congenital heart malformation. We identified a new locus linked with atrial septal defect on chromosome 14q12 in a large family with dominantly inherited atrial septal defect. The underlying mutation is a missense substitution, I820N, in alpha-myosin heavy chain (MYH6), a structural protein expressed at high levels in the developing atria, which affects the binding of the heavy chain to its regulatory light chain. The cardiac transcription factor TBX5 strongly regulates expression of MYH6, but mutant forms of TBX5, which cause Holt-Oram syndrome, do not. Morpholino knock-down of expression of the chick MYH6 homolog eliminates the formation of the atrial septum without overtly affecting atrial chamber formation. These data provide evidence for a link between a transcription factor, a structural protein and congenital heart disease.

Low-frequency intermediate penetrance variants in the ROCK1 gene predispose to Tetralogy of Fallot.

BACKGROUND: Epidemiological studies indicate a substantial excess familial recurrence of non-syndromic Tetralogy of Fallot (TOF), implicating genetic factors that remain largely unknown. The Rho induced kinase 1 gene (ROCK1) is a key component of the planar cell polarity signalling pathway, which plays an important role in normal cardiac development. The aim of this study was to investigate the role of genetic variation in ROCK1 on the risk of TOF. RESULTS: ROCK1 was sequenced in a discovery cohort of 93 non-syndromic TOF probands to identify rare variants. TagSNPs were selected to capture commoner variation in ROCK1. Novel variants and TagSNPs were genotyped in a discovery cohort of 458 TOF cases and 1331 healthy controls, and positive findings were replicated in a further 209 TOF cases and 1290 healthy controls. Association between genotypes and TOF was assessed using LAMP.A rare SNP (c.807C > T; rs56085230) discovered by sequencing was associated with TOF risk (p = 0.006) in the discovery cohort. The variant was also significantly associated with the risk of TOF in the replication cohort (p = 0.018). In the combined cohorts the odds ratio for TOF was 2.61 (95% CI 1.58-4.30); p < 0.0001. The minor allele frequency of rs56085230 in the cases was 0.02, and in the controls it was 0.007. The variant accounted for 1% of the population attributable risk (PAR) of TOF. We also found significant association with TOF for an uncommon TagSNP in ROCK1, rs288979 (OR 1.64 [95% CI 1.15-2.30]; p = 1.5x10⁻5). The minor allele frequency of rs288979 in the controls was 0.043, and the variant accounted for 11% of the PAR of TOF. These association signals were independent of each other, providing additional internal validation of our result. CONCLUSIONS: Low frequency intermediate penetrance (LFIP) variants in the ROCK1 gene predispose to the risk of TOF.

Expanded CUG repeats Dysregulate RNA splicing by altering the stoichiometry of the muscleblind 1 complex.

To understand the role of the splice regulator muscleblind 1 (MBNL1) in the development of RNA splice defects in myotonic dystrophy I (DM1), we purified RNA-independent MBNL1 complexes from normal human myoblasts and examined the behavior of these complexes in DM1 myoblasts. Antibodies recognizing MBNL1 variants (MBNL1(CUG)), which can sequester in the toxic CUG RNA foci that develop in DM1 nuclei, were used to purify MBNL1(CUG) complexes from normal myoblasts. In normal myoblasts, MBNL1(CUG) bind 10 proteins involved in remodeling ribonucleoprotein complexes including hnRNP H, H2, H3, F, A2/B1, K, L, DDX5, DDX17, and DHX9. Of these proteins, only MBNL1(CUG) colocalizes extensively with DM1 CUG foci (>80% of foci) with its partners being present in <10% of foci. Importantly, the stoichiometry of MBNL1(CUG) complexes is altered in DM1 myoblasts, demonstrating an increase in the steady state levels of nine of its partner proteins. These changes are recapitulated by the expression of expanded CUG repeat RNA in Cos7 cells. Altered stoichiometry of MBNL1(CUG) complexes results from aberrant protein synthesis or stability and is unlinked to PKCα function. Modeling these changes in normal myoblasts demonstrates that increased levels of hnRNP H, H2, H3, F, and DDX5 independently dysregulate splicing in overlapping RNA subsets. Thus expression of expanded CUG repeats alters the stoichiometry of MBNL1(CUG) complexes to allow both the reinforcement and expansion of RNA processing defects.

Muscleblind-like proteins: similarities and differences in normal and myotonic dystrophy muscle.

In myotonic dystrophy, muscleblind-like protein 1 (MBNL1) protein binds specifically to expanded CUG or CCUG repeats, which accumulate as discrete nuclear foci, and this is thought to prevent its function in the regulation of alternative splicing of pre-mRNAs. There is strong evidence for the role of the MBNL1 gene in disease pathology, but the roles of two related genes, MBNL2 and MBNL3, are less clear. Using new monoclonal antibodies specific for each of the three gene products, we found that MBNL2 decreased during human fetal development and myoblast culture, while MBNL1 was unchanged. In Duchenne muscular dystrophy muscle, MBNL2 was elevated in immature, regenerating fibres compared with mature fibres, supporting some developmental role for MBNL2. MBNL3 was found only in C2C12 mouse myoblasts. Both MBNL1 and MBNL2 were partially sequestered by nuclear foci of expanded repeats in adult muscle and cultured cells from myotonic dystrophy patients. In adult muscle nucleoplasm, both proteins were reduced in myotonic dystrophy type 1 compared with an age-matched control. In normal human myoblast cultures, MBNL1 and MBNL2 always co-distributed but their distribution could change rapidly from nucleoplasmic to cytoplasmic. Functional differences between MBNL1 and MBNL2 have not yet been found and may prove quite subtle. The dominance of MBNL1 in mature, striated muscle would explain why ablation of the mouse mbnl1 gene alone is sufficient to cause a myotonic dystrophy.

Systems genetics analysis identifies calcium-signaling defects as novel cause of congenital heart disease.

BACKGROUND: Congenital heart disease (CHD) occurs in almost 1% of newborn children and is considered a multifactorial disorder. CHD may segregate in families due to significant contribution of genetic factors in the disease etiology. The aim of the study was to identify pathophysiological mechanisms in families segregating CHD. METHODS: We used whole exome sequencing to identify rare genetic variants in ninety consenting participants from 32 Danish families with recurrent CHD. We applied a systems biology approach to identify developmental mechanisms influenced by accumulation of rare variants. We used an independent cohort of 714 CHD cases and 4922 controls for replication and performed functional investigations using zebrafish as in vivo model. RESULTS: We identified 1785 genes, in which rare alleles were shared between affected individuals within a family. These genes were enriched for known cardiac developmental genes, and 218 of these genes were mutated in more than one family. Our analysis revealed a functional cluster, enriched for proteins with a known participation in calcium signaling. Replication in an independent cohort confirmed increased mutation burden of calcium-signaling genes in CHD patients. Functional investigation of zebrafish orthologues of ITPR1, PLCB2, and ADCY2 verified a role in cardiac development and suggests a combinatorial effect of inactivation of these genes. CONCLUSIONS: The study identifies abnormal calcium signaling as a novel pathophysiological mechanism in human CHD and confirms the complex genetic architecture underlying CHD.

CDK12 inhibition reduces abnormalities in cells from patients with myotonic dystrophy and in a mouse model.

Myotonic dystrophy type 1 (DM1) is an RNA-based disease with no current treatment. It is caused by a transcribed CTG repeat expansion within the 3' untranslated region of the dystrophia myotonica protein kinase (DMPK) gene. Mutant repeat expansion transcripts remain in the nuclei of patients' cells, forming distinct microscopically detectable foci that contribute substantially to the pathophysiology of the condition. Here, we report small-molecule inhibitors that remove nuclear foci and have beneficial effects in the HSALR mouse model, reducing transgene expression, leading to improvements in myotonia, splicing, and centralized nuclei. Using chemoproteomics in combination with cell-based assays, we identify cyclin-dependent kinase 12 (CDK12) as a druggable target for this condition. CDK12 is a protein elevated in DM1 cell lines and patient muscle biopsies, and our results showed that its inhibition led to reduced expression of repeat expansion RNA. Some of the inhibitors identified in this study are currently the subject of clinical trials for other indications and provide valuable starting points for a drug development program in DM1.

Knockdown of alpha myosin heavy chain disrupts the cytoskeleton and leads to multiple defects during chick cardiogenesis.

Atrial septal defects are a common congenital heart defect in humans. Although mutations in different genes are now frequently being described, little is known about the processes and mechanisms behind the early stages of atrial septal development. By utilizing morpholino-induced knockdown in the chick we have analysed the role of alpha myosin heavy chain during early cardiogenesis in a temporal manner. Upon knockdown of alpha myosin heavy chain, three different phenotypes of the atrial septum were observed: (1) the atrial septum failed to initiate, (2) the septum was initiated but was growth restricted, or (3) incorrect specification occurred resulting in multiple septa forming. In addition, at a lower frequency, decreased alpha myosin heavy chain was found to give rise to an abnormally looped heart or an enlarged heart. Staining of the actin cytoskeleton indicated that many of the myofibrils in the knockdown hearts were not as mature as those observed in the controls, suggesting a mechanism for the defects seen. Therefore, these data suggest a role for alpha myosin heavy chain in modelling of the early heart and the range of defects to the atrial septum suggest roles in its initiation, specification and growth during development.