Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency.

[Figure: see text].

Robust identification of deletions in exome and genome sequence data based on clustering of Mendelian errors.

Multiple tools have been developed to identify copy number variants (CNVs) from whole exome (WES) and whole genome sequencing (WGS) data. Current tools such as XHMM for WES and CNVnator for WGS identify CNVs based on changes in read depth. For WGS, other methods to identify CNVs include utilizing discordant read pairs and split reads and genome-wide local assembly with tools such as Lumpy and SvABA, respectively. Here, we introduce a new method to identify deletion CNVs from WES and WGS trio data based on the clustering of Mendelian errors (MEs). Using our Mendelian Error Method (MEM), we identified 127 deletions (inherited and de novo) in 2,601 WES trios from the Pediatric Cardiac Genomics Consortium, with a validation rate of 88% by digital droplet PCR. MEM identified additional de novo deletions compared with XHMM, and a significant enrichment of 15q11.2 deletions compared with controls. In addition, MEM identified eight cases of uniparental disomy, sample switches, and DNA contamination. We applied MEM to WGS data from the Genome In A Bottle Ashkenazi trio and identified deletions with 97% specificity. MEM provides a robust, computationally inexpensive method for identifying deletions, and an orthogonal approach for verifying deletions called by other tools.

Robust identification of mosaic variants in congenital heart disease.

Mosaicism due to somatic mutations can cause multiple diseases including cancer, developmental and overgrowth syndromes, neurodevelopmental disorders, autoinflammatory diseases, and atrial fibrillation. With the increased use of next generation sequencing technology, multiple tools have been developed to identify low-frequency variants, specifically from matched tumor-normal tissues in cancer studies. To investigate whether mosaic variants are implicated in congenital heart disease (CHD), we developed a pipeline using the cancer somatic variant caller MuTect to identify mosaic variants in whole-exome sequencing (WES) data from a cohort of parent/affected child trios (n = 715) and a cohort of healthy individuals (n = 416). This is a novel application of the somatic variant caller designed for cancer to WES trio data. We identified two cases with mosaic KMT2D mutations that are likely pathogenic for CHD, but conclude that, overall, mosaicism detectable in peripheral blood or saliva does not account for a significant portion of CHD etiology.

Increased frequency of de novo copy number variants in congenital heart disease by integrative analysis of single nucleotide polymorphism array and exome sequence data.

RATIONALE: Congenital heart disease (CHD) is among the most common birth defects. Most cases are of unknown pathogenesis. OBJECTIVE: To determine the contribution of de novo copy number variants (CNVs) in the pathogenesis of sporadic CHD. METHODS AND RESULTS: We studied 538 CHD trios using genome-wide dense single nucleotide polymorphism arrays and whole exome sequencing. Results were experimentally validated using digital droplet polymerase chain reaction. We compared validated CNVs in CHD cases with CNVs in 1301 healthy control trios. The 2 complementary high-resolution technologies identified 63 validated de novo CNVs in 51 CHD cases. A significant increase in CNV burden was observed when comparing CHD trios with healthy trios, using either single nucleotide polymorphism array (P=7×10(-5); odds ratio, 4.6) or whole exome sequencing data (P=6×10(-4); odds ratio, 3.5) and remained after removing 16% of de novo CNV loci previously reported as pathogenic (P=0.02; odds ratio, 2.7). We observed recurrent de novo CNVs on 15q11.2 encompassing CYFIP1, NIPA1, and NIPA2 and single de novo CNVs encompassing DUSP1, JUN, JUP, MED15, MED9, PTPRE SREBF1, TOP2A, and ZEB2, genes that interact with established CHD proteins NKX2-5 and GATA4. Integrating de novo variants in whole exome sequencing and CNV data suggests that ETS1 is the pathogenic gene altered by 11q24.2-q25 deletions in Jacobsen syndrome and that CTBP2 is the pathogenic gene in 10q subtelomeric deletions. CONCLUSIONS: We demonstrate a significantly increased frequency of rare de novo CNVs in CHD patients compared with healthy controls and suggest several novel genetic loci for CHD.

CANOES: detecting rare copy number variants from whole exome sequencing data.

We present CANOES, an algorithm for the detection of rare copy number variants from exome sequencing data. CANOES models read counts using a negative binomial distribution and estimates variance of the read counts using a regression-based approach based on selected reference samples in a given dataset. We test CANOES on a family-based exome sequencing dataset, and show that its sensitivity and specificity is comparable to that of XHMM. Moreover, the method is complementary to Gaussian approximation-based methods (e.g. XHMM or CoNIFER). When CANOES is used in combination with these methods, it will be possible to produce high accuracy calls, as demonstrated by a much reduced and more realistic de novo rate in results from trio data.

Remote Cardiac Safety Monitoring through the Lens of the FDA Biomarker Qualification Evidentiary Criteria Framework: A Case Study Analysis.

Clinical safety findings remain one of the reasons for attrition of drug candidates during clinical development. Cardiovascular liabilities are not consistently detected in early-stage clinical trials and often become apparent when drugs are administered chronically for extended periods of time. Vital sign data collection outside of the clinic offers an opportunity for deeper physiological characterization of drug candidates and earlier safety signal detection. A working group representing expertise from biopharmaceutical and technology sectors, US Food and Drug Administration (FDA) public-private partnerships, academia, and regulators discussed and presented a remote cardiac monitoring case study at the FNIH Biomarkers Consortium Remote Digital Monitoring for Medical Product Development workshop to examine applicability of the biomarker qualification evidentiary framework by the FDA. This use case examined the components of the framework, including the statement of need, the context of use, the state of the evidence, and the benefit/risk profile. Examination of results from 2 clinical trials deploying 510(k)-cleared devices for remote cardiac data collection demonstrated the need for analytical and clinical validity irrespectively of the regulatory status of a device of interest, emphasizing the importance of data collection method assessment in the context of intended use. Additionally, collection of large amounts of ambulatory data also highlighted the need for new statistical methods and contextual information to enable data interpretation. A wider adoption of this approach for drug development purposes will require collaborations across industry, academia, and regulatory agencies to establish methodologies and supportive data sets to enable data interpretation and decision-making.

De Novo Damaging Variants, Clinical Phenotypes, and Post-Operative Outcomes in Congenital Heart Disease.

BACKGROUND: De novo genic and copy number variants are enriched in patients with congenital heart disease, particularly those with extra-cardiac anomalies. The impact of de novo damaging variants on outcomes following cardiac repair is unknown. METHODS: We studied 2517 patients with congenital heart disease who had undergone whole-exome sequencing as part of the CHD GENES study (Congenital Heart Disease Genetic Network). RESULTS: Two hundred ninety-four patients (11.7%) had clinically significant de novo variants. Patients with de novo damaging variants were 2.4 times more likely to have extra-cardiac anomalies (P=5.63×10-12). In 1268 patients (50.4%) who had surgical data available and underwent open-heart surgery exclusive of heart transplantation as their first operation, we analyzed transplant-free survival following the first operation. Median follow-up was 2.65 years. De novo variants were associated with worse transplant-free survival (hazard ratio, 3.51; P=5.33×10-04) and longer times to final extubation (hazard ratio, 0.74; P=0.005). As de novo variants had a significant interaction with extra-cardiac anomalies for transplant-free survival (P=0.003), de novo variants conveyed no additional risk for transplant-free survival for patients with these anomalies (adjusted hazard ratio, 1.96; P=0.06). By contrast, de novo variants in patients without extra-cardiac anomalies were associated with worse transplant-free survival during follow-up (hazard ratio, 11.21; P=1.61×10-05) than that of patients with no de novo variants. Using agnostic machine-learning algorithms, we identified de novo copy number variants at 15q25.2 and 15q11.2 as being associated with worse transplant-free survival and 15q25.2, 22q11.21, and 3p25.2 as being associated with prolonged time to final extubation. CONCLUSIONS: In patients with congenital heart disease undergoing open-heart surgery, de novo variants were associated with worse transplant-free survival and longer times on the ventilator. De novo variants were most strongly associated with adverse outcomes among patients without extra-cardiac anomalies, suggesting a benefit for preoperative genetic testing even when genetic abnormalities are not suspected during routine clinical practice. Registration: URL: https://www.clinicaltrials.gov. Unique identifier: NCT01196182.

Genomic analyses implicate noncoding de novo variants in congenital heart disease.

A genetic etiology is identified for one-third of patients with congenital heart disease (CHD), with 8% of cases attributable to coding de novo variants (DNVs). To assess the contribution of noncoding DNVs to CHD, we compared genome sequences from 749 CHD probands and their parents with those from 1,611 unaffected trios. Neural network prediction of noncoding DNV transcriptional impact identified a burden of DNVs in individuals with CHD (n = 2,238 DNVs) compared to controls (n = 4,177; P = 8.7 × 10-4). Independent analyses of enhancers showed an excess of DNVs in associated genes (27 genes versus 3.7 expected, P = 1 × 10-5). We observed significant overlap between these transcription-based approaches (odds ratio (OR) = 2.5, 95% confidence interval (CI) 1.1-5.0, P = 5.4 × 10-3). CHD DNVs altered transcription levels in 5 of 31 enhancers assayed. Finally, we observed a DNV burden in RNA-binding-protein regulatory sites (OR = 1.13, 95% CI 1.1-1.2, P = 8.8 × 10-5). Our findings demonstrate an enrichment of potentially disruptive regulatory noncoding DNVs in a fraction of CHD at least as high as that observed for damaging coding DNVs.

Continuous Monitoring Using a Wearable Device Detects Activity-Induced Heart Rate Changes After Administration of Amphetamine.

Wearable digital devices offer potential advantages over traditional methods for the collection of health-related information, including continuous collection of dense data while study subjects are ambulatory or in remote settings. We assessed the utility of collecting continuous actigraphy and cardiac monitoring by deploying two US Food and Drug Administration (FDA) 510(k)-cleared devices in a phase I clinical trial of a novel compound, which included the use of an amphetamine challenge. The Phillips Actiwatch Spectrum Pro (Actiwatch) was used to assess mobility and sleep. The Preventice BodyGuardian (BodyGuardian) was used for monitoring heart rate (HR) and respiratory rate (RR), via single-lead electrocardiogram (ECG) recordings, together with physical activity. We measured data collection rates, compared device readouts with conventional measures, and monitored changes in HR measures during the amphetamine challenge. Completeness of data collection was good for the Actiwatch (96%) and lower for the BodyGuardian (80%). A good correlation was observed between device and in-clinic measures for HR (r = 0.99; P < 0.001), but was poor for RR (r = 0.39; P = 0.004). Manual reviews of selected ECG strips corresponding to HR measures below, within, and above the normal range were consistent with BodyGuardian measurements. The BodyGuardian device detected clear HR responses after amphetamine administration while subjects were physically active, whereas conventional measures collected at predefined timepoints while subjects were resting and supine did not. Wearable digital technology shows promise for monitoring human subjects for physiologic changes and pharmacologic responses, although fit-for-purpose evaluation and validation continues to be important prior to the wider deployment of these devices.

Paternal-age-related de novo mutations and risk for five disorders.

There are established associations between advanced paternal age and offspring risk for psychiatric and developmental disorders. These are commonly attributed to genetic mutations, especially de novo single nucleotide variants (dnSNVs), that accumulate with increasing paternal age. However, the actual magnitude of risk from such mutations in the male germline is unknown. Quantifying this risk would clarify the clinical significance of delayed paternity. Using parent-child trio whole-exome-sequencing data, we estimate the relationship between paternal-age-related dnSNVs and risk for five disorders: autism spectrum disorder (ASD), congenital heart disease, neurodevelopmental disorders with epilepsy, intellectual disability and schizophrenia (SCZ). Using Danish registry data, we investigate whether epidemiologic associations between each disorder and older fatherhood are consistent with the estimated role of dnSNVs. We find that paternal-age-related dnSNVs confer a small amount of risk for these disorders. For ASD and SCZ, epidemiologic associations with delayed paternity reflect factors that may not increase with age.

Contribution of rare inherited and de novo variants in 2,871 congenital heart disease probands.

Congenital heart disease (CHD) is the leading cause of mortality from birth defects. Here, exome sequencing of a single cohort of 2,871 CHD probands, including 2,645 parent-offspring trios, implicated rare inherited mutations in 1.8%, including a recessive founder mutation in GDF1 accounting for ∼5% of severe CHD in Ashkenazim, recessive genotypes in MYH6 accounting for ∼11% of Shone complex, and dominant FLT4 mutations accounting for 2.3% of Tetralogy of Fallot. De novo mutations (DNMs) accounted for 8% of cases, including ∼3% of isolated CHD patients and ∼28% with both neurodevelopmental and extra-cardiac congenital anomalies. Seven genes surpassed thresholds for genome-wide significance, and 12 genes not previously implicated in CHD had >70% probability of being disease related. DNMs in ∼440 genes were inferred to contribute to CHD. Striking overlap between genes with damaging DNMs in probands with CHD and autism was also found.

Single-Cell Resolution of Temporal Gene Expression during Heart Development.

Activation of complex molecular programs in specific cell lineages governs mammalian heart development, from a primordial linear tube to a four-chamber organ. To characterize lineage-specific, spatiotemporal developmental programs, we performed single-cell RNA sequencing of >1,200 murine cells isolated at seven time points spanning embryonic day 9.5 (primordial heart tube) to postnatal day 21 (mature heart). Using unbiased transcriptional data, we classified cardiomyocytes, endothelial cells, and fibroblast-enriched cells, thus identifying markers for temporal and chamber-specific developmental programs. By harnessing these datasets, we defined developmental ages of human and mouse pluripotent stem-cell-derived cardiomyocytes and characterized lineage-specific maturation defects in hearts of mice with heterozygous mutations in Nkx2.5 that cause human heart malformations. This spatiotemporal transcriptome analysis of heart development reveals lineage-specific gene programs underlying normal cardiac development and congenital heart disease.

Loss of RNA expression and allele-specific expression associated with congenital heart disease.

Congenital heart disease (CHD), a prevalent birth defect occurring in 1% of newborns, likely results from aberrant expression of cardiac developmental genes. Mutations in a variety of cardiac transcription factors, developmental signalling molecules and molecules that modify chromatin cause at least 20% of disease, but most CHD remains unexplained. We employ RNAseq analyses to assess allele-specific expression (ASE) and biallelic loss-of-expression (LOE) in 172 tissue samples from 144 surgically repaired CHD subjects. Here we show that only 5% of known imprinted genes with paternal allele silencing are monoallelic versus 56% with paternal allele expression-this cardiac-specific phenomenon seems unrelated to CHD. Further, compared with control subjects, CHD subjects have a significant burden of both LOE genes and ASE events associated with altered gene expression. These studies identify FGFBP2, LBH, RBFOX2, SGSM1 and ZBTB16 as candidate CHD genes because of significantly altered transcriptional expression.

Interpreting de novo Variation in Human Disease Using denovolyzeR.

Spontaneously arising (de novo) genetic variants are important in human disease, yet every individual carries many such variants, with a median of 1 de novo variant affecting the protein-coding portion of the genome. A recently described mutational model provides a powerful framework for the robust statistical evaluation of such coding variants, enabling the interpretation of de novo variation in human disease. Here we describe a new open-source software package, denovolyzeR, that implements this model and provides tools for the analysis of de novo coding sequence variants.

HEART DISEASE. Titin mutations in iPS cells define sarcomere insufficiency as a cause of dilated cardiomyopathy.

Human mutations that truncate the massive sarcomere protein titin [TTN-truncating variants (TTNtvs)] are the most common genetic cause for dilated cardiomyopathy (DCM), a major cause of heart failure and premature death. Here we show that cardiac microtissues engineered from human induced pluripotent stem (iPS) cells are a powerful system for evaluating the pathogenicity of titin gene variants. We found that certain missense mutations, like TTNtvs, diminish contractile performance and are pathogenic. By combining functional analyses with RNA sequencing, we explain why truncations in the A-band domain of TTN cause DCM, whereas truncations in the I band are better tolerated. Finally, we demonstrate that mutant titin protein in iPS cell-derived cardiomyocytes results in sarcomere insufficiency, impaired responses to mechanical and ß-adrenergic stress, and attenuated growth factor and cell signaling activation. Our findings indicate that titin mutations cause DCM by disrupting critical linkages between sarcomerogenesis and adaptive remodeling.

De novo mutations in congenital heart disease with neurodevelopmental and other congenital anomalies.

Congenital heart disease (CHD) patients have an increased prevalence of extracardiac congenital anomalies (CAs) and risk of neurodevelopmental disabilities (NDDs). Exome sequencing of 1213 CHD parent-offspring trios identified an excess of protein-damaging de novo mutations, especially in genes highly expressed in the developing heart and brain. These mutations accounted for 20% of patients with CHD, NDD, and CA but only 2% of patients with isolated CHD. Mutations altered genes involved in morphogenesis, chromatin modification, and transcriptional regulation, including multiple mutations in RBFOX2, a regulator of mRNA splicing. Genes mutated in other cohorts examined for NDD were enriched in CHD cases, particularly those with coexisting NDD. These findings reveal shared genetic contributions to CHD, NDD, and CA and provide opportunities for improved prognostic assessment and early therapeutic intervention in CHD patients.

Digital Droplet PCR: CNV Analysis and Other Applications.

Digital droplet PCR (ddPCR) is an assay that combines state-of-the-art microfluidics technology with TaqMan-based PCR to achieve precise target DNA quantification at high levels of sensitivity and specificity. Because quantification is achieved without the need for standard assays in an easy to interpret, unambiguous digital readout, ddPCR is far simpler, faster, and less error prone than real-time qPCR. The basic protocol can be modified with minor adjustments to suit a wide range of applications, such as CNV analysis, rare variant detection, SNP genotyping, and transcript quantification. This unit describes the ddPCR workflow in detail for the Bio-Rad QX100 system, but the theory and data interpretation are generalizable to any ddPCR system.

Targeted and genome-wide sequencing reveal single nucleotide variations impacting specificity of Cas9 in human stem cells.

CRISPR/Cas9 has demonstrated a high-efficiency in site-specific gene targeting. However, potential off-target effects of the Cas9 nuclease represent a major safety concern for any therapeutic application. Here, we knock out the Tafazzin gene by CRISPR/Cas9 in human-induced pluripotent stem cells with 54% efficiency. We combine whole-genome sequencing and deep-targeted sequencing to characterise the off-target effects of Cas9 editing. Whole-genome sequencing of Cas9-modified hiPSC clones detects neither gross genomic alterations nor elevated mutation rates. Deep sequencing of in silico predicted off-target sites in a population of Cas9-treated cells further confirms high specificity of Cas9. However, we identify a single high-efficiency off-target site that is generated by a common germline single-nucleotide variant (SNV) in our experiment. Based on in silico analysis, we estimate a likelihood of SNVs creating off-target sites in a human genome to be ~1.5-8.5%, depending on the genome and site-selection method, but also note that mutations might be generated at these sites only at low rates and may not have functional consequences. Our study demonstrates the feasibility of highly specific clonal ex vivo gene editing using CRISPR/Cas9 and highlights the value of whole-genome sequencing before personalised CRISPR design.

JNK signaling coordinates integrin and actin functions during Drosophila embryogenesis.

Epithelial movements are key morphogenetic events in animal development. They are driven by multiple mechanisms, including signal-dependent changes in cytoskeletal organization and in cell adhesion. Such processes must be controlled precisely and coordinated to accurately sculpt the three-dimensional form of the developing organism. By observing the Drosophila epidermis during embryonic development using confocal time-lapse microscopy, we have investigated how signaling through the Jun-N-terminal kinase (JNK) pathway governs the tissue sheet movements that result in dorsal closure (DC). We find that JNK controls the polymerization of actin into a cable at the epidermal leading edge as previously suggested, as well as the joining (zipping) of the contralateral epithelial cell sheets. Here, we show that zipping is mediated by regulation of the integrins myospheroid and scab. Our data demonstrate that JNK signaling regulates a set of target genes that cooperate to facilitate epithelial movement and closure.