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.

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

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.

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.

HDAC4 and 5 repression of TBX5 is relieved by protein kinase D1.

TBX5 is a T-box family transcription factor that regulates heart and forelimb development in vertebrates and functional deficiencies in this protein result in Holt-Oram syndrome. Recently, we have shown that acetylation of TBX5 potentiates its activity and is important for heart and limb development. Here we report that class II histone deacetylases HDAC4 and HDAC5 associate with TBX5 and repress its role in cardiac gene transcription. Both HDAC4 and HDAC5 deacetylate TBX5, which promotes its relocation to the cytoplasm and HDAC4 antagonizes the physical association and functional cooperation between TBX5 and MEF2C. We also show that protein kinase D1 (PRKD1) relieves the HDAC4/5-mediated repression of TBX5. Thus, this study reveals a novel interaction of HDAC4/5 and PRKD1 in the regulation of TBX5 transcriptional activity.

Global Increase in Circular RNA Levels in Myotonic Dystrophy.

Splicing aberrations induced as a consequence of the sequestration of muscleblind-like splicing factors on the dystrophia myotonica protein kinase transcript, which contains expanded CUG repeats, present a major pathomechanism of myotonic dystrophy type 1 (DM1). As muscleblind-like factors may also be important factors involved in the biogenesis of circular RNAs (circRNAs), we hypothesized that the level of circRNAs would be decreased in DM1. To test this hypothesis, we selected 20 well-validated circRNAs and analyzed their levels in several experimental systems (e.g., cell lines, DM muscle tissues, and a mouse model of DM1) using droplet digital PCR assays. We also explored the global level of circRNAs using two RNA-Seq datasets of DM1 muscle samples. Contrary to our original hypothesis, our results consistently showed a global increase in circRNA levels in DM1, and we identified numerous circRNAs that were increased in DM1. We also identified many genes (including muscle-specific genes) giving rise to numerous (>10) circRNAs. Thus, this study is the first to show an increase in global circRNA levels in DM1. We also provided preliminary results showing the association of circRNA level with muscle weakness and alternative splicing changes that are biomarkers of DM1 severity.

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.

Quantitative Methods to Monitor RNA Biomarkers in Myotonic Dystrophy.

Myotonic dystrophy type 1 (DM1) and type 2 (DM2) are human neuromuscular disorders associated with mutations of simple repetitive sequences in affected genes. The abnormal expansion of CTG repeats in the 3'-UTR of the DMPK gene elicits DM1, whereas elongated CCTG repeats in intron 1 of ZNF9/CNBP triggers DM2. Pathogenesis of both disorders is manifested by nuclear retention of expanded repeat-containing RNAs and aberrant alternative splicing. The precise determination of absolute numbers of mutant RNA molecules is important for a better understanding of disease complexity and for accurate evaluation of the efficacy of therapeutic drugs. We present two quantitative methods, Multiplex Ligation-Dependent Probe Amplification and droplet digital PCR, for studying the mutant DMPK transcript (DMPKexpRNA) and the aberrant alternative splicing in DM1 and DM2 human tissues and cells. We demonstrate that in DM1, the DMPKexpRNA is detected in higher copy number than its normal counterpart. Moreover, the absolute number of the mutant transcript indicates its low abundance with only a few copies per cell in DM1 fibroblasts. Most importantly, in conjunction with fluorescence in-situ hybridization experiments, our results suggest that in DM1 fibroblasts, the vast majority of nuclear RNA foci consist of a few molecules of DMPKexpRNA.

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.

Recent advances in congenital heart disease genomics.

Congenital heart disease is the most common congenital abnormality, and advances in medical care mean that this population of individuals is surviving for longer than ever before. It represents a significant healthcare challenge, as many patients require life-long care and individuals may ask about the likelihood of their children being affected. Whilst a number of genes have been identified previously from investigation of families with Mendelian inheritance patterns, sequencing the DNA from large cohorts of individuals with congenital heart disease is now providing fresh insights into the genetics of these conditions. This research has enabled novel gene discovery and uncovered the different genetic mechanisms underlying both isolated congenital heart disease and that which occurs in association with other medical problems. This article discusses the most recent advances in this field and the implications for patient care. In addition, we consider the challenges facing researchers in this field and emphasise the need for close working relationships between clinicians and researchers.

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.

Body composition and clinical outcome measures in patients with myotonic dystrophy type 1.

Myotonic Dystrophy type 1 (DM1) is the most common form of adult onset muscular dystrophy. In this study we compared body composition in DM1 and matched controls and evaluated the relationship between these parameters and clinical outcome measures in DM1 patients. In addition we established the sensitivity to change of these measures in a prospective 18 month longitudinal study of the DM1 patient cohort. Clinical data, manual muscle testing (MMT), quantitative muscle testing (QMT) of ankle dorsiflexion, bilateral grip dynamometry, 6 minute walk test and a DM1 functional rating scale (DM1-Activ) were collected at baseline (n = 38) and 18 month follow-up (n = 36). The case-control analysis was performed comparing baseline data with 31 anthropometrically matched controls. Dual-energy X-ray absorptiometry (DEXA) was used to obtain regional measurements of fat-free mass index (FFMI) and fat mass index (FMI) and demonstrated significant reduction of FFMI in the legs (left p = 0.004; right p = 0.017) and trunk (p < 0.0001) and increased FMI localised to the trunk (p < 0.0001) in DM1 patients compared to controls. Regional left and right arm FFMI and FMI significantly positively and negatively correlated with grip strength and both total FFMI (p = 0.0009) and FMI (p = 0.02) decreased and increased by 0.38 kg/m2 and 0.31 kg/m2 respectively after 18 month follow-up. DEXA is likely to provide a useful secondary outcome measurement of disease progression in addition to muscle strength and timed functional tasks in clinical trials.

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.

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.

Association between C677T polymorphism of methylene tetrahydrofolate reductase and congenital heart disease: meta-analysis of 7697 cases and 13,125 controls.

BACKGROUND: Association between the C677T polymorphism of the methylene tetrahydrofolate reductase (MTHFR) gene and congenital heart disease (CHD) is contentious. METHODS AND RESULTS: We compared genotypes between CHD cases and controls and between mothers of CHD cases and controls. We placed our results in context by conducting meta-analyses of previously published studies. Among 5814 cases with primary genotype data and 10 056 controls, there was no evidence of association between MTHFR C677T genotype and CHD risk (odds ratio [OR], 0.96 [95% confidence interval, 0.87-1.07]). A random-effects meta-analysis of all studies (involving 7697 cases and 13 125 controls) suggested the presence of association (OR, 1.25 [95% confidence interval, 1.03-1.51]; P=0.022) but with substantial heterogeneity among contributing studies (I(2)=64.4%) and evidence of publication bias. Meta-analysis of large studies only (defined by a variance of the log OR <0.05), which together contributed 83% of all cases, yielded no evidence of association (OR, 0.97 [95% confidence interval, 0.91-1.03]) without significant heterogeneity (I(2)=0). Moreover, meta-analysis of 1781 mothers of CHD cases (829 of whom were genotyped in this study) and 19 861 controls revealed no evidence of association between maternal C677T genotype and risk of CHD in offspring (OR, 1.13 [95% confidence interval, 0.87-1.47]). There was no significant association between MTHFR genotype and CHD risk in large studies from regions with different levels of dietary folate. CONCLUSIONS: The MTHFR C677T polymorphism, which directly influences plasma folate levels, is not associated with CHD risk. Publication biases appear to substantially contaminate the literature with regard to this genetic association.

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.