Development of a comprehensive cardiovascular disease genetic risk assessment test.

Background: Despite monogenic and polygenic contributions to cardiovascular disease (CVD), genetic testing is not widely adopted, and current tests are limited by the breadth of surveyed conditions and interpretation burden. Methods: We developed a comprehensive clinical genome CVD test with semi-automated interpretation. Monogenic conditions and risk alleles were selected based on the strength of disease association and evidence for increased disease risk, respectively. Non-CVD secondary findings genes, pharmacogenomic (PGx) variants and CVD polygenic risk scores (PRS) were assessed for inclusion. Test performance was modeled using 2,594 genomes from the 1000 Genomes Project, and further investigated in 20 previously tested individuals. Results: The CVD genome test is composed of a panel of 215 CVD gene-disease pairs, 35 non-CVD secondary findings genes, 4 risk alleles or genotypes, 10 PGx genes and a PRS for coronary artery disease. Modeling of test performance using samples from the 1000 Genomes Project revealed ~6% of individuals with a monogenic finding in a CVD-associated gene, 6% with a risk allele finding, ~1% with a non-CVD secondary finding, and 93% with CVD-associated PGx variants. Assessment of blinded clinical samples showed complete concordance with prior testing. An average of 4 variants were reviewed per case, with interpretation and reporting time ranging from 9-96 min. Conclusions: A genome sequencing based CVD genetic risk assessment can provide comprehensive genetic disease and genetic risk information to patients with CVD. The semi-automated and limited interpretation burden suggest that this testing approach could be scaled to support population-level initiatives.

Spinal muscular atrophy diagnosis and carrier screening from genome sequencing data.

PURPOSE: Spinal muscular atrophy (SMA), caused by loss of the SMN1 gene, is a leading cause of early childhood death. Due to the near identical sequences of SMN1 and SMN2, analysis of this region is challenging. Population-wide SMA screening to quantify the SMN1 copy number (CN) is recommended by the American College of Medical Genetics and Genomics. METHODS: We developed a method that accurately identifies the CN of SMN1 and SMN2 using genome sequencing (GS) data by analyzing read depth and eight informative reference genome differences between SMN1/2. RESULTS: We characterized SMN1/2 in 12,747 genomes, identified 1568 samples with SMN1 gains or losses and 6615 samples with SMN2 gains or losses, and calculated a pan-ethnic carrier frequency of 2%, consistent with previous studies. Additionally, 99.8% of our SMN1 and 99.7% of SMN2 CN calls agreed with orthogonal methods, with a recall of 100% for SMA and 97.8% for carriers, and a precision of 100% for both SMA and carriers. CONCLUSION: This SMN copy-number caller can be used to identify both carrier and affected status of SMA, enabling SMA testing to be offered as a comprehensive test in neonatal care and an accurate carrier screening tool in GS sequencing projects.

Pathogenic variants in USP7 cause a neurodevelopmental disorder with speech delays, altered behavior, and neurologic anomalies.

PURPOSE: Haploinsufficiency of USP7, located at chromosome 16p13.2, has recently been reported in seven individuals with neurodevelopmental phenotypes, including developmental delay/intellectual disability (DD/ID), autism spectrum disorder (ASD), seizures, and hypogonadism. Further, USP7 was identified to critically incorporate into the MAGEL2-USP7-TRIM27 (MUST), such that pathogenic variants in USP7 lead to altered endosomal F-actin polymerization and dysregulated protein recycling. METHODS: We report 16 newly identified individuals with heterozygous USP7 variants, identified by genome or exome sequencing or by chromosome microarray analysis. Clinical features were evaluated by review of medical records. Additional clinical information was obtained on the seven previously reported individuals to fully elucidate the phenotypic expression associated with USP7 haploinsufficiency. RESULTS: The clinical manifestations of these 23 individuals suggest a syndrome characterized by DD/ID, hypotonia, eye anomalies,feeding difficulties, GERD, behavioral anomalies, and ASD, and more specific phenotypes of speech delays including a nonverbal phenotype and abnormal brain magnetic resonance image findings including white matter changes based on neuroradiologic examination. CONCLUSION: The consistency of clinical features among all individuals presented regardless of de novo USP7 variant type supports haploinsufficiency as a mechanism for pathogenesis and refines the clinical impact faced by affected individuals and caregivers.

Copy-number variants in clinical genome sequencing: deployment and interpretation for rare and undiagnosed disease.

PURPOSE: Current diagnostic testing for genetic disorders involves serial use of specialized assays spanning multiple technologies. In principle, genome sequencing (GS) can detect all genomic pathogenic variant types on a single platform. Here we evaluate copy-number variant (CNV) calling as part of a clinically accredited GS test. METHODS: We performed analytical validation of CNV calling on 17 reference samples, compared the sensitivity of GS-based variants with those from a clinical microarray, and set a bound on precision using orthogonal technologies. We developed a protocol for family-based analysis of GS-based CNV calls, and deployed this across a clinical cohort of 79 rare and undiagnosed cases. RESULTS: We found that CNV calls from GS are at least as sensitive as those from microarrays, while only creating a modest increase in the number of variants interpreted (~10 CNVs per case). We identified clinically significant CNVs in 15% of the first 79 cases analyzed, all of which were confirmed by an orthogonal approach. The pipeline also enabled discovery of a uniparental disomy (UPD) and a 50% mosaic trisomy 14. Directed analysis of select CNVs enabled breakpoint level resolution of genomic rearrangements and phasing of de novo CNVs. CONCLUSION: Robust identification of CNVs by GS is possible within a clinical testing environment.

Congenital sodium diarrhea and chorioretinal coloboma with optic disc coloboma in a patient with biallelic SPINT2 mutations, including p.(Tyr163Cys).

Congenital sodium diarrhea is a rare and life-threatening disorder characterized by a severe, secretory diarrhea containing high concentrations of sodium, leading to hyponatremia and metabolic acidosis. It may occur in isolation or in association with systemic features such as facial dysmorphism, choanal atresia, imperforate anus, and corneal erosions. Mutations in the serine protease inhibitor, Kunitz-Type 2 (SPINT2) gene have been associated with congenital sodium diarrhea and additional syndromic features. We present a child with congenital sodium diarrhea, cleft lip and palate, corneal erosions, optic nerve coloboma, and intermittent exotropia who was found to have biallelic mutations in SPINT2. One mutation, c.488A > G, predicting p.(Tyr163Cys), has been previously associated with a syndromic form of congenital sodium diarrhea. The other mutation, c.166_167dupTA, predicting p.(Asn57Thrfs*24) has not previously been reported and is likely a novel pathogenic variant for this disorder. We found only one other report of an optic nerve coloboma associated with SPINT2 mutations and this occurred in a patient with congenital tufting enteropathy. Our patient confirms an association of ocular coloboma with presumed loss of SPINT2 function.

Attenuation of yeast UPR is essential for survival and is mediated by IRE1 kinase.

The unfolded protein response (UPR) activates Ire1, an endoplasmic reticulum (ER) resident transmembrane kinase and ribonuclease (RNase), in response to ER stress. We used an in vivo assay, in which disappearance of the UPR-induced spliced HAC1 messenger ribonucleic acid (mRNA) correlates with the recovery of the ER protein-folding capacity, to investigate the attenuation of the UPR in yeast. We find that, once activated, spliced HAC1 mRNA is sustained in cells expressing Ire1 carrying phosphomimetic mutations within the kinase activation loop, suggesting that dephosphorylation of Ire1 is an important step in RNase deactivation. Additionally, spliced HAC1 mRNA is also sustained after UPR induction in cells expressing Ire1 with mutations in the conserved DFG kinase motif (D828A) or a conserved residue (F842) within the activation loop. The importance of proper Ire1 RNase attenuation is demonstrated by the inability of cells expressing Ire1-D828A to grow under ER stress. We propose that the activity of the Ire1 kinase domain plays a role in attenuating its RNase activity when ER function is recovered.