High-fidelity, large-scale targeted profiling of microsatellites.

Microsatellites are highly mutable sequences that can serve as markers for relationships among individuals or cells within a population. The accuracy and resolution of reconstructing these relationships depends on the fidelity of microsatellite profiling and the number of microsatellites profiled. However, current methods for targeted profiling of microsatellites incur significant "stutter" artifacts that interfere with accurate genotyping, and sequencing costs preclude whole-genome microsatellite profiling of a large number of samples. We developed a novel method for accurate and cost-effective targeted profiling of a panel of > 150,000 microsatellites per sample, along with a computational tool for designing large-scale microsatellite panels. Our method addresses the greatest challenge for microsatellite profiling - "stutter" artifacts - with a low-temperature hybridization capture that significantly reduces these artifacts. We also developed a computational tool for accurate genotyping of the resulting microsatellite sequencing data that uses an ensemble approach integrating three microsatellite genotyping tools, which we optimize by analysis of de novo microsatellite mutations in human trios. Altogether, our suite of experimental and computational tools enables high-fidelity, large-scale profiling of microsatellites, which may find utility in diverse applications such as lineage tracing, population genetics, ecology, and forensics.

Ultra-Rapid Droplet Digital PCR Enables Intraoperative Tumor Quantification.

The diagnosis and treatment of tumors often depends on molecular-genetic data. However, rapid and iterative access to molecular data is not currently feasible during surgery, complicating intraoperative diagnosis and precluding measurement of tumor cell burdens at surgical margins to guide resections. To address this gap, we developed Ultra-Rapid droplet digital PCR (UR-ddPCR), which can be completed in 15 minutes from tissue to result with an accuracy comparable to standard ddPCR. We demonstrate UR-ddPCR assays for the IDH1 R132H and BRAF V600E clonal mutations that are present in many low-grade gliomas and melanomas, respectively. We illustrate the clinical feasibility of UR-ddPCR by performing it intraoperatively for 13 glioma cases. We further combine UR-ddPCR measurements with UR-stimulated Raman histology intraoperatively to estimate tumor cell densities in addition to tumor cell percentages. We anticipate that UR-ddPCR, along with future refinements in assay instrumentation, will enable novel point-of-care diagnostics and the development of molecularly-guided surgeries that improve clinical outcomes.

DNA mismatch and damage patterns revealed by single-molecule sequencing.

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

Single duplex DNA sequencing with CODEC detects mutations with high sensitivity.

Detecting mutations from single DNA molecules is crucial in many fields but challenging. Next-generation sequencing (NGS) affords tremendous throughput but cannot directly sequence double-stranded DNA molecules ('single duplexes') to discern the true mutations on both strands. Here we present Concatenating Original Duplex for Error Correction (CODEC), which confers single duplex resolution to NGS. CODEC affords 1,000-fold higher accuracy than NGS, using up to 100-fold fewer reads than duplex sequencing. CODEC revealed mutation frequencies of 2.72 × 10-8 in sperm of a 39-year-old individual, and somatic mutations acquired with age in blood cells. CODEC detected genome-wide, clonal hematopoiesis mutations from single DNA molecules, single mutated duplexes from tumor genomes and liquid biopsies, microsatellite instability with 10-fold greater sensitivity and mutational signatures, and specific tumor mutations with up to 100-fold fewer reads. CODEC enables more precise genetic testing and reveals biologically significant mutations, which are commonly obscured by NGS errors.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing.

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Biochemical characterization of two novel mutations in the human high-affinity choline transporter 1 identified in a patient with congenital myasthenic syndrome.

Congenital myasthenic syndrome (CMS) is a heterogeneous condition associated with 34 different genes, including SLC5A7, which encodes the high-affinity choline transporter 1 (CHT1). CHT1 is expressed in presynaptic neurons of the neuromuscular junction where it uses the inward sodium gradient to reuptake choline. Biallelic CHT1 mutations often lead to neonatal lethality, and less commonly to non-lethal motor weakness and developmental delays. Here, we report detailed biochemical characterization of two novel mutations in CHT1, p.I294T and p.D349N, which we identified in an 11-year-old patient with a history of neonatal respiratory distress, and subsequent hypotonia and global developmental delay. Heterologous expression of each CHT1 mutant in human embryonic kidney cells showed two different mechanisms of reduced protein function. The p.I294T CHT1 mutant transporter function was detectable, but its abundance and half-life were significantly reduced. In contrast, the p.D349N CHT1 mutant was abundantly expressed at the cell membrane, but transporter function was absent. The residual function of the p.I294T CHT1 mutant may explain the non-lethal form of CMS in this patient, and the divergent mechanisms of reduced CHT1 function that we identified may guide future functional studies of the CHT1 myasthenic syndrome. Based on these in vitro studies that provided a diagnosis, treatment with cholinesterase inhibitor together with physical and occupational therapy significantly improved the patient's strength and quality of life.

Serial enrichment of heteroduplex DNA using a MutS-magnetic bead system.

Numerous applications in molecular biology and genomics require characterization of mutant DNA molecules present at low levels within a larger sample of non-mutant DNA. This is often achieved either by selectively amplifying mutant DNA, or by sequencing all the DNA followed by computational identification of the mutant DNA. However, selective amplification is challenging for insertions and deletions (indels). Additionally, sequencing all the DNA in a sample may not be cost effective when only the presence of a mutation needs to be ascertained rather than its allelic fraction. The MutS protein evolved to detect DNA heteroduplexes in which the two DNA strands are mismatched. Prior methods have utilized MutS to enrich mutant DNA by hybridizing mutant to non-mutant DNA to create heteroduplexes. However, the purity of heteroduplex DNA these methods achieve is limited because they can only feasibly perform one or two enrichment cycles. We developed a MutS-magnetic bead system that enables rapid serial enrichment cycles. With six cycles, we achieve complete purification of heteroduplex indel DNA originally present at a 5% fraction and over 40-fold enrichment of heteroduplex DNA originally present at a 1% fraction. This system may enable novel approaches for enriching mutant DNA for targeted sequencing.

An HNRNPK-specific DNA methylation signature makes sense of missense variants and expands the phenotypic spectrum of Au-Kline syndrome.

Au-Kline syndrome (AKS) is a neurodevelopmental disorder associated with multiple malformations and a characteristic facial gestalt. The first individuals ascertained carried de novo loss-of-function (LoF) variants in HNRNPK. Here, we report 32 individuals with AKS (26 previously unpublished), including 13 with de novo missense variants. We propose new clinical diagnostic criteria for AKS that differentiate it from the clinically overlapping Kabuki syndrome and describe a significant phenotypic expansion to include individuals with missense variants who present with subtle facial features and few or no malformations. Many gene-specific DNA methylation (DNAm) signatures have been identified for neurodevelopmental syndromes. Because HNRNPK has roles in chromatin and epigenetic regulation, we hypothesized that pathogenic variants in HNRNPK may be associated with a specific DNAm signature. Here, we report a unique DNAm signature for AKS due to LoF HNRNPK variants, distinct from controls and Kabuki syndrome. This DNAm signature is also identified in some individuals with de novo HNRNPK missense variants, confirming their pathogenicity and the phenotypic expansion of AKS to include more subtle phenotypes. Furthermore, we report that some individuals with missense variants have an "intermediate" DNAm signature that parallels their milder clinical presentation, suggesting the presence of an epi-genotype phenotype correlation. In summary, the AKS DNAm signature may help elucidate the underlying pathophysiology of AKS. This DNAm signature also effectively supported clinical syndrome delineation and is a valuable aid for variant interpretation in individuals where a clinical diagnosis of AKS is unclear, particularly for mild presentations.

Applications of Single-Cell DNA Sequencing.

Over the past decade, genomic analyses of single cells-the fundamental units of life-have become possible. Single-cell DNA sequencing has shed light on biological questions that were previously inaccessible across diverse fields of research, including somatic mutagenesis, organismal development, genome function, and microbiology. Single-cell DNA sequencing also promises significant future biomedical and clinical impact, spanning oncology, fertility, and beyond. While single-cell approaches that profile RNA and protein have greatly expanded our understanding of cellular diversity, many fundamental questions in biology and important biomedical applications require analysis of the DNA of single cells. Here, we review the applications and biological questions for which single-cell DNA sequencing is uniquely suited or required. We include a discussion of the fields that will be impacted by single-cell DNA sequencing as the technology continues to advance.

Transcriptome sequencing identifies a noncoding, deep intronic variant in CLCN7 causing autosomal recessive osteopetrosis.

BACKGROUND: Over half of children with rare genetic diseases remain undiagnosed despite maximal clinical evaluation and DNA-based genetic testing. As part of an Undiagnosed Diseases Program applying transcriptome (RNA) sequencing to identify the causes of these unsolved cases, we studied a child with severe infantile osteopetrosis leading to cranial nerve palsies, bone deformities, and bone marrow failure, for whom whole-genome sequencing was nondiagnostic. METHODS: We performed transcriptome (RNA) sequencing of whole blood followed by analysis of aberrant transcript isoforms and osteoclast functional studies. RESULTS: We identified a pathogenic deep intronic variant in CLCN7 creating an unexpected, frameshifting pseudoexon causing complete loss of function. Functional studies, including osteoclastogenesis and bone resorption assays, confirmed normal osteoclast differentiation but loss of osteoclast function. CONCLUSION: This is the first report of a pathogenic deep intronic variant in CLCN7, and our approach provides a model for systematic identification of noncoding variants causing osteopetrosis-a disease for which molecular-genetic diagnosis can be pivotal for potentially curative hematopoietic stem cell transplantation. Our work illustrates that cryptic splice variants may elude DNA-only sequencing and supports broad first-line use of transcriptome sequencing for children with undiagnosed diseases.

A Modern History RAAS Inhibition and Beta Blockade for Heart Failure to Underscore the Non-equivalency of ACEIs and ARBs.

Beta blockers and renin-angiotensin-aldosterone-inhibitors (RAAS-i) including angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin II receptor blockers (ARBs) have been a mainstay of guideline-based medical therapy for heart failure with reduced ejection fraction (HFrEF) for decades. However, initial evidence supporting each of the aforenoted class of drug for heart failure indications was largely found independently of the other two classes with the exception of the addition of BBs to ACEIs. In the initial ACEI trials for HFrEF, few participants were on BBs as BBs were seen as contraindicated in HFrEF at the time. The seminal BB in HFrEF trials had high prevalence of ACEIs use as ACEIs for HF were standard of care by then, but ARBs as a class were still in their infancy. We closely examine the evidence for combinations of BB and ACEIs versus ARBs in HFrEF. In doing so, we demonstrate the lack of evidence for consideration of ARBs to be interchangeable with ACEIs when used in combination with BB and provide evidence that calls in to question the validity of assuming benefits from each drug class are independently cumulative, widening the gap between ACEIs and ARBs when used with BBs. Modern guidelines should emphasize this lack of evidence for the combination use of ARB and BB in HFrEF, except for candesartan. Even as practice moves towards the widespread uptake of angiotensin receptor-neprilysin inhibitors (which contain the ARB valsartan) in heart failure, the distinction has important implications for the ongoing role of combination therapy with BB, which thus far has been assumed, but not proven.

DPH1 syndrome: two novel variants and structural and functional analyses of seven missense variants identified in syndromic patients.

DPH1 variants have been associated with an ultra-rare and severe neurodevelopmental disorder, mainly characterized by variable developmental delay, short stature, dysmorphic features, and sparse hair. We have identified four new patients (from two different families) carrying novel variants in DPH1, enriching the clinical delineation of the DPH1 syndrome. Using a diphtheria toxin ADP-ribosylation assay, we have analyzed the activity of seven identified variants and demonstrated compromised function for five of them [p.(Leu234Pro); p.(Ala411Argfs*91); p.(Leu164Pro); p.(Leu125Pro); and p.(Tyr112Cys)]. We have built a homology model of the human DPH1-DPH2 heterodimer and have performed molecular dynamics simulations to study the effect of these variants on the catalytic sites as well as on the interactions between subunits of the heterodimer. The results show correlation between loss of activity, reduced size of the opening to the catalytic site, and changes in the size of the catalytic site with clinical severity. This is the first report of functional tests of DPH1 variants associated with the DPH1 syndrome. We demonstrate that the in vitro assay for DPH1 protein activity, together with structural modeling, are useful tools for assessing the effect of the variants on DPH1 function and may be used for predicting patient outcomes and prognoses.

Correction: DPH1 syndrome: two novel variants and structural and functional analyses of seven missense variants identified in syndromic patients.

Following the publication of the article, it was noted that the last column in Table 1, the total % should have read 5/8 (62.5) for the 'Epilepsy' row, and not 5.7 (71.4). This has now been amended in the HTML and PDF of the original article.

Integrated genome and transcriptome sequencing identifies a noncoding mutation in the genome replication factor DONSON as the cause of microcephaly-micromelia syndrome.

While next-generation sequencing has accelerated the discovery of human disease genes, progress has been largely limited to the "low hanging fruit" of mutations with obvious exonic coding or canonical splice site impact. In contrast, the lack of high-throughput, unbiased approaches for functional assessment of most noncoding variants has bottlenecked gene discovery. We report the integration of transcriptome sequencing (RNA-seq), which surveys all mRNAs to reveal functional impacts of variants at the transcription level, into the gene discovery framework for a unique human disease, microcephaly-micromelia syndrome (MMS). MMS is an autosomal recessive condition described thus far in only a single First Nations population and causes intrauterine growth restriction, severe microcephaly, craniofacial anomalies, skeletal dysplasia, and neonatal lethality. Linkage analysis of affected families, including a very large pedigree, identified a single locus on Chromosome 21 linked to the disease (LOD > 9). Comprehensive genome sequencing did not reveal any pathogenic coding or canonical splicing mutations within the linkage region but identified several nonconserved noncoding variants. RNA-seq analysis detected aberrant splicing in DONSON due to one of these noncoding variants, showing a causative role for DONSON disruption in MMS. We show that DONSON is expressed in progenitor cells of embryonic human brain and other proliferating tissues, is co-expressed with components of the DNA replication machinery, and that Donson is essential for early embryonic development in mice as well, suggesting an essential conserved role for DONSON in the cell cycle. Our results demonstrate the utility of integrating transcriptomics into the study of human genetic disease when DNA sequencing alone is not sufficient to reveal the underlying pathogenic mutation.

Ticagrelor and Rosuvastatin Have Additive Cardioprotective Effects via Adenosine.

BACKGROUND: Ticagrelor inhibits the equilibrative-nucleoside-transporter-1 and thereby, adenosine cell re-uptake. Ticagrelor limits infarct size (IS) in non-diabetic rats and the effect is adenosine-dependent. Statins, via ecto-5'-nucleotidase activation, also increase adenosine levels and limit IS. HYPOTHESIS: Ticagrelor and rosuvastatin have additive effects on myocardial adenosine levels, and therefore, on IS and post-reperfusion activation of the NLRP3-inflammasome. METHODS: Diabetic ZDF rats received via oral gavage; water (control), ticagrelor (150 mg/kg/d), prasugrel (7.5 mg/kg/d), rosuvastatin (5 mg/kg/d), ticagrelor + rosuvastatin and prasugrel + rosuvastatin for 3d. On day 4, rats underwent 30 min coronary artery occlusion and 24 h of reperfusion. Two additional groups received, ticagrelor + rosuvastatin or water in combination with CGS15943 (CGS, an adenosine receptor antagonist, 10 mg/kg i.p. 1 h before ischemia). RESULTS: Both ticagrelor and rosuvastatin increased myocardial adenosine levels with an additive effect of the combination whereas prasugrel had no effect. Similarly, both ticagrelor and rosuvastatin significantly reduced IS with an additive effect of the combination whereas prasugrel had no effect. The effect on IS was adenosine dependent as CGS15943 reversed the effect of ticagrelor + rosuvastatin. The ischemia-reperfusion injury increased myocardial mRNA levels of NLRP3, ASC, IL-1ß and IL-6. Ticagrelor and rosuvastatin, but not prasugrel, significantly decreased these pro-inflammatory mediators with a trend to an additive effect of the combination. The combination also increased the levels of anti-inflammatory 15-epilipoxin A4. CONCLUSIONS: Ticagrelor and rosuvastatin when given in combination have an additive effect on local myocardial adenosine levels in the setting of ischemia reperfusion. This translates into an additive cardioprotective effect mediated by adenosine-induced effects including downregulation of pro- but upregulation of anti-inflammatory mediators.

Resolving rates of mutation in the brain using single-neuron genomics.

Whether somatic mutations contribute functional diversity to brain cells is a long-standing question. Single-neuron genomics enables direct measurement of somatic mutation rates in human brain and promises to answer this question. A recent study (Upton et al., 2015) reported high rates of somatic LINE-1 element (L1) retrotransposition in the hippocampus and cerebral cortex that would have major implications for normal brain function, and suggested that these events preferentially impact genes important for neuronal function. We identify aspects of the single-cell sequencing approach, bioinformatic analysis, and validation methods that led to thousands of artifacts being interpreted as somatic mutation events. Our reanalysis supports a mutation frequency of approximately 0.2 events per cell, which is about fifty-fold lower than reported, confirming that L1 elements mobilize in some human neurons but indicating that L1 mosaicism is not ubiquitous. Through consideration of the challenges identified, we provide a foundation and framework for designing single-cell genomics studies.

Somatic mutation in single human neurons tracks developmental and transcriptional history.

Neurons live for decades in a postmitotic state, their genomes susceptible to DNA damage. Here we survey the landscape of somatic single-nucleotide variants (SNVs) in the human brain. We identified thousands of somatic SNVs by single-cell sequencing of 36 neurons from the cerebral cortex of three normal individuals. Unlike germline and cancer SNVs, which are often caused by errors in DNA replication, neuronal mutations appear to reflect damage during active transcription. Somatic mutations create nested lineage trees, allowing them to be dated relative to developmental landmarks and revealing a polyclonal architecture of the human cerebral cortex. Thus, somatic mutations in the brain represent a durable and ongoing record of neuronal life history, from development through postmitotic function.