Spatial Mapping of Hematopoietic Clones in Human Bone Marrow.

Clonal hematopoiesis (CH) is the expansion of somatically mutated cells in the hematopoietic compartment of individuals without hematopoietic dysfunction. Large CH clones (i.e., >2% variant allele fraction) predispose to hematologic malignancy, but CH is detected at lower levels in nearly all middle-aged individuals. Prior work has extensively characterized CH in peripheral blood, but the spatial distribution of hematopoietic clones in human bone marrow is largely undescribed. To understand CH at this level, we developed a method for spatially aware somatic mutation profiling and characterized the bone marrow of a patient with polycythemia vera. We identified the complex clonal distribution of somatic mutations in the hematopoietic compartment, the restriction of somatic mutations to specific subpopulations of hematopoietic cells, and spatial constraints of these clones in the bone marrow. This proof of principle paves the way to answering fundamental questions regarding CH spatial organization and factors driving CH expansion and malignant transformation in the bone marrow. SIGNIFICANCE: CH occurs commonly in humans and can predispose to hematologic malignancy. Although well characterized in blood, it is poorly understood how clones are spatially distributed in the bone marrow. To answer this, we developed methods for spatially aware somatic mutation profiling to describe clonal heterogeneity in human bone marrow. See related commentary by Austin and Aifantis, p. 139.

Droplet Digital PCR for Oncogenic KMT2A Fusion Detection.

Acute myeloid leukemia (AML) is an aggressive blood cancer diagnosed in approximately 120,000 individuals worldwide each year. During treatment for AML, detecting residual disease is essential for prognostication and treatment decision-making. Currently, methods for detecting residual AML are limited to identifying approximately 1:100 to 1:1000 leukemic cells (morphology and DNA sequencing) or are difficult to implement (flow cytometry). AML arising after chemotherapy or radiation exposure is termed therapy-related AML (t-AML) and is exceptionally aggressive and treatment resistant. t-AML is often driven by oncogenic fusions that result from prior treatments that introduce double-strand DNA breaks. The most common t-AML-associated translocations affect KMT2A. There are at least 80 known KMT2A fusion partners, but approximately 80% of fusions involve only five partners-AF9, AF6, AF4, ELL, and ENL. We present a novel droplet digital PCR assay targeting the most common KMT2A-rearrangements to enable detection of rare AML cells harboring these fusions. This assay was benchmarked in cell lines and patient samples harboring oncogenic KMT2A fusions and demonstrated a limit of detection of approximately 1:1,000,000 cells. Future application of this assay could improve disease detection and treatment decision-making for patients with t-AML with KMT2A fusions and premalignant oncogenic fusion detection in at-risk individuals after chemotherapy exposure.

KDM6B protects T-ALL cells from NOTCH1-induced oncogenic stress.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematopoietic neoplasm resulting from the malignant transformation of T-cell progenitors. While activating NOTCH1 mutations are the dominant genetic drivers of T-ALL, epigenetic dysfunction plays a central role in the pathology of T-ALL and can provide alternative mechanisms to oncogenesis in lieu of or in combination with genetic mutations. The histone demethylase enzyme KDM6A (UTX) is also recurrently mutated in T-ALL patients and functions as a tumor suppressor. However, its gene paralog, KDM6B (JMJD3), is never mutated and can be significantly overexpressed, suggesting it may be necessary for sustaining the disease. Here, we used mouse and human T-ALL models to show that KDM6B is required for T-ALL development and maintenance. Using NOTCH1 gain-of-function retroviral models, mouse cells genetically deficient for Kdm6b were unable to propagate T-ALL. Inactivating KDM6B in human T-ALL patient cells by CRISPR/Cas9 showed KDM6B-targeted cells were significantly outcompeted over time. The dependence of T-ALL cells on KDM6B was proportional to the oncogenic strength of NOTCH1 mutation, with KDM6B required to prevent stress-induced apoptosis from strong NOTCH1 signaling. These studies identify a crucial role for KDM6B in sustaining NOTCH1-driven T-ALL and implicate KDM6B as a novel therapeutic target in these patients.

Effect of Clonal Hematopoiesis on Cardiovascular Disease in People Living with HIV.

Hematopoietic stem cells (HSCs) with age-associated somatic mutations that disproportionally contribute to hematopoiesis generate the condition known as clonal hematopoiesis (CH). While CH conveys increased risk of hematologic cancer, there is also strong association between CH and cardiovascular disease (CVD). Accumulating evidence suggests that inflammation mechanistically links CH to CVD, and we hypothesized that CH may be a predictive biomarker of CVD in conditions of chronic inflammation. One such patient population comprises people living with HIV (PLWH) who also have substantially increased incidences of CVD and CH . We studied the association between CH and CVD in PLWH using samples from ACTG Study A5001 (or ALLRT), a prospective clinical trial of HIV-infected persons with long-term follow-up. We observed a positive association between CH and CVD in PLWH independent of traditional CVD risk factors. Moreover, in CVD cases, the CH clone was identifiable in the blood years before CVD diagnosis, unlike in PLWH with CH who did not have CVD. With the life span of PLWH increasing because of advances in treatment, our results indicate that the presence of CH and its clonal dynamics could be used as a prognostic biomarker of the risk for CVD in PLWH.

Cancer therapy shapes the fitness landscape of clonal hematopoiesis.

Acquired mutations are pervasive across normal tissues. However, understanding of the processes that drive transformation of certain clones to cancer is limited. Here we study this phenomenon in the context of clonal hematopoiesis (CH) and the development of therapy-related myeloid neoplasms (tMNs). We find that mutations are selected differentially based on exposures. Mutations in ASXL1 are enriched in current or former smokers, whereas cancer therapy with radiation, platinum and topoisomerase II inhibitors preferentially selects for mutations in DNA damage response genes (TP53, PPM1D, CHEK2). Sequential sampling provides definitive evidence that DNA damage response clones outcompete other clones when exposed to certain therapies. Among cases in which CH was previously detected, the CH mutation was present at tMN diagnosis. We identify the molecular characteristics of CH that increase risk of tMN. The increasing implementation of clinical sequencing at diagnosis provides an opportunity to identify patients at risk of tMN for prevention strategies.

The evolutionary dynamics and fitness landscape of clonal hematopoiesis.

Somatic mutations acquired in healthy tissues as we age are major determinants of cancer risk. Whether variants confer a fitness advantage or rise to detectable frequencies by chance remains largely unknown. Blood sequencing data from ~50,000 individuals reveal how mutation, genetic drift, and fitness shape the genetic diversity of healthy blood (clonal hematopoiesis). We show that positive selection, not drift, is the major force shaping clonal hematopoiesis, provide bounds on the number of hematopoietic stem cells, and quantify the fitness advantages of key pathogenic variants, at single-nucleotide resolution, as well as the distribution of fitness effects (fitness landscape) within commonly mutated driver genes. These data are consistent with clonal hematopoiesis being driven by a continuing risk of mutations and clonal expansions that become increasingly detectable with age.

The Zebrafish Xenograft Platform-A Novel Tool for Modeling KSHV-Associated Diseases.

Kaposi's sarcoma associated-herpesvirus (KSHV, also known as human herpesvirus-8) is a gammaherpesvirus that establishes life-long infection in human B lymphocytes. KSHV infection is typically asymptomatic, but immunosuppression can predispose KSHV-infected individuals to primary effusion lymphoma (PEL); a malignancy driven by aberrant proliferation of latently infected B lymphocytes, and supported by pro-inflammatory cytokines and angiogenic factors produced by cells that succumb to lytic viral replication. Here, we report the development of the first in vivo model for a virally induced lymphoma in zebrafish, whereby KSHV-infected PEL tumor cells engraft and proliferate in the yolk sac of zebrafish larvae. Using a PEL cell line engineered to produce the viral lytic switch protein RTA in the presence of doxycycline, we demonstrate drug-inducible reactivation from KSHV latency in vivo, which enabled real-time observation and evaluation of latent and lytic phases of KSHV infection. In addition, we developed a sensitive droplet digital PCR method to monitor latent and lytic viral gene expression and host cell gene expression in xenografts. The zebrafish yolk sac is not well vascularized, and by using fluorogenic assays, we confirmed that this site provides a hypoxic environment that may mimic the microenvironment of some human tumors. We found that PEL cell proliferation in xenografts was dependent on the host hypoxia-dependent translation initiation factor, eukaryotic initiation factor 4E2 (eIF4E2). This demonstrates that the zebrafish yolk sac is a functionally hypoxic environment, and xenografted cells must switch to dedicated hypoxic gene expression machinery to survive and proliferate. The establishment of the PEL xenograft model enables future studies that exploit the innate advantages of the zebrafish as a model for genetic and pharmacologic screens.

Clonal hematopoiesis and risk of acute myeloid leukemia.

Nearly all adults harbor acute myeloid leukemia (AML)-related clonal hematopoietic mutations at a variant allele fraction (VAF) of ≥0.0001, yet relatively few develop hematologic malignancies. We conducted a nested analysis in the Nurses' Health Study and Health Professionals Follow-Up Study blood subcohorts, with up to 22 years of follow up to investigate associations of clonal mutations of ≥0.0001 allele frequency with future risk of AML. We identified 35 cases with AML that had pre-diagnosis peripheral blood samples and matched two controls without history of cancer per case by sex, age, and ethnicity. We conducted blinded error-corrected sequencing on all study samples and assessed variant-associated risk using conditional logistic regression. We detected AML-associated mutations in 97% of all participants (598 mutations, 5.8/person). Individuals with mutations ≥0.01 variant allele fraction had a significantly increased AML risk (OR 5.4, 95%CI: 1.8-16.6), as did individuals with higher-frequency clones and those with DNMT3A R882H/C mutations. The risk of lower-frequency clones was less clear. In the 11 case-control sets with samples banked ten years apart, clonal mutations rarely expanded over time. Our findings are consistent with published evidence that detection of clonal mutations ≥0.01 VAF identifies individuals at increased risk for AML. Further study of larger populations, mutations co-occurring within the same pre-leukemic clone and other risk factors (lifestyle, epigenetics, etc.), are still needed to fully elucidate the risk conferred by low-frequency clonal hematopoiesis in asymptomatic adults.

Rare Event Detection Using Error-corrected DNA and RNA Sequencing.

Conventional next-generation sequencing techniques (NGS) have allowed for immense genomic characterization for over a decade. Specifically, NGS has been used to analyze the spectrum of clonal mutations in malignancy. Though far more efficient than traditional Sanger methods, NGS struggles with identifying rare clonal and subclonal mutations due to its high error rate of ~0.5-2.0%. Thus, standard NGS has a limit of detection for mutations that are >0.02 variant allele fraction (VAF). While the clinical significance for mutations this rare in patients without known disease remains unclear, patients treated for leukemia have significantly improved outcomes when residual disease is <0.0001 by flow cytometry. In order to mitigate this artefactual background of NGS, numerous methods have been developed. Here we describe a method for Error-corrected DNA and RNA Sequencing (ECS), which involves tagging individual molecules with both a 16 bp random index for error-correction and an 8 bp patient-specific index for multiplexing. Our method can detect and track clonal mutations at variant allele fractions (VAFs) two orders of magnitude lower than the detection limit of NGS and as rare as 0.0001 VAF.

Clonal haematopoiesis harbouring AML-associated mutations is ubiquitous in healthy adults.

Clonal haematopoiesis is thought to be a rare condition that increases in frequency with age and predisposes individuals to haematological malignancy. Recent studies, utilizing next-generation sequencing (NGS), observed haematopoietic clones in 10% of 70-year olds and rarely in younger individuals. However, these studies could only detect common haematopoietic clones->0.02 variant allele fraction (VAF)-due to the error rate of NGS. To identify and characterize clonal mutations below this threshold, here we develop methods for targeted error-corrected sequencing, which enable the accurate detection of clonal mutations as rare as 0.0003 VAF. We apply these methods to study serially banked peripheral blood samples from healthy 50-60-year-old participants in the Nurses' Health Study. We observe clonal haematopoiesis, frequently harbouring mutations in DNMT3A and TET2, in 95% of individuals studied. These clonal mutations are often stable longitudinally and present in multiple haematopoietic compartments, suggesting a long-lived haematopoietic stem and progenitor cell of origin.

Pooled Sequencing of Candidate Genes Implicates Rare Variants in the Development of Asthma Following Severe RSV Bronchiolitis in Infancy.

Severe infection with respiratory syncytial virus (RSV) during infancy is strongly associated with the development of asthma. To identify genetic variation that contributes to asthma following severe RSV bronchiolitis during infancy, we sequenced the coding exons of 131 asthma candidate genes in 182 European and African American children with severe RSV bronchiolitis in infancy using anonymous pools for variant discovery, and then directly genotyped a set of 190 nonsynonymous variants. Association testing was performed for physician-diagnosed asthma before the 7th birthday (asthma) using genotypes from 6,500 individuals from the Exome Sequencing Project (ESP) as controls to gain statistical power. In addition, among patients with severe RSV bronchiolitis during infancy, we examined genetic associations with asthma, active asthma, persistent wheeze, and bronchial hyperreactivity (methacholine PC20) at age 6 years. We identified four rare nonsynonymous variants that were significantly associated with asthma following severe RSV bronchiolitis, including single variants in ADRB2, FLG and NCAM1 in European Americans (p = 4.6x10-4, 1.9x10-13 and 5.0x10-5, respectively), and NOS1 in African Americans (p = 2.3x10-11). One of the variants was a highly functional nonsynonymous variant in ADRB2 (rs1800888), which was also nominally associated with asthma (p = 0.027) and active asthma (p = 0.013) among European Americans with severe RSV bronchiolitis without including the ESP. Our results suggest that rare nonsynonymous variants contribute to the development of asthma following severe RSV bronchiolitis in infancy, notably in ADRB2. Additional studies are required to explore the role of rare variants in the etiology of asthma and asthma-related traits following severe RSV bronchiolitis.

Internal derangements of the temporomandibular joint: A review of the anatomy, diagnosis, and management.

Internal derangements of the temporomandibular joint are conditions in which the articular disc has become displaced from its original position the condylar head. Relevant anatomic structures and their functional relationships are briefly discussed. The displacement of the disc can result in numerous presentations, with the most common being disc displacement with reduction (with or without intermittent locking), and disc displacement without reduction (with or without limited opening). These are described in this article according to the standardized Diagnostic Criteria for Temporomandibular Disorders, as well as the less common posterior disc displacement. Appropriate management usually ranges from patient education and monitoring to splints, physical therapy, and medications. In rare and select cases, surgery may be necessary. However, in for the majority of internal derangements, the prognosis is good, particularly with conservative care.

Role of TP53 mutations in the origin and evolution of therapy-related acute myeloid leukaemia.

Therapy-related acute myeloid leukaemia (t-AML) and therapy-related myelodysplastic syndrome (t-MDS) are well-recognized complications of cytotoxic chemotherapy and/or radiotherapy. There are several features that distinguish t-AML from de novo AML, including a higher incidence of TP53 mutations, abnormalities of chromosomes 5 or 7, complex cytogenetics and a reduced response to chemotherapy. However, it is not clear how prior exposure to cytotoxic therapy influences leukaemogenesis. In particular, the mechanism by which TP53 mutations are selectively enriched in t-AML/t-MDS is unknown. Here, by sequencing the genomes of 22 patients with t-AML, we show that the total number of somatic single-nucleotide variants and the percentage of chemotherapy-related transversions are similar in t-AML and de novo AML, indicating that previous chemotherapy does not induce genome-wide DNA damage. We identified four cases of t-AML/t-MDS in which the exact TP53 mutation found at diagnosis was also present at low frequencies (0.003-0.7%) in mobilized blood leukocytes or bone marrow 3-6 years before the development of t-AML/t-MDS, including two cases in which the relevant TP53 mutation was detected before any chemotherapy. Moreover, functional TP53 mutations were identified in small populations of peripheral blood cells of healthy chemotherapy-naive elderly individuals. Finally, in mouse bone marrow chimaeras containing both wild-type and Tp53(+/-) haematopoietic stem/progenitor cells (HSPCs), the Tp53(+/-) HSPCs preferentially expanded after exposure to chemotherapy. These data suggest that cytotoxic therapy does not directly induce TP53 mutations. Rather, they support a model in which rare HSPCs carrying age-related TP53 mutations are resistant to chemotherapy and expand preferentially after treatment. The early acquisition of TP53 mutations in the founding HSPC clone probably contributes to the frequent cytogenetic abnormalities and poor responses to chemotherapy that are typical of patients with t-AML/t-MDS.

Population-based rare variant detection via pooled exome or custom hybridization capture with or without individual indexing.

BACKGROUND: Rare genetic variation in the human population is a major source of pathophysiological variability and has been implicated in a host of complex phenotypes and diseases. Finding disease-related genes harboring disparate functional rare variants requires sequencing of many individuals across many genomic regions and comparing against unaffected cohorts. However, despite persistent declines in sequencing costs, population-based rare variant detection across large genomic target regions remains cost prohibitive for most investigators. In addition, DNA samples are often precious and hybridization methods typically require large amounts of input DNA. Pooled sample DNA sequencing is a cost and time-efficient strategy for surveying populations of individuals for rare variants. We set out to 1) create a scalable, multiplexing method for custom capture with or without individual DNA indexing that was amenable to low amounts of input DNA and 2) expand the functionality of the SPLINTER algorithm for calling substitutions, insertions and deletions across either candidate genes or the entire exome by integrating the variant calling algorithm with the dynamic programming aligner, Novoalign. RESULTS: We report methodology for pooled hybridization capture with pre-enrichment, indexed multiplexing of up to 48 individuals or non-indexed pooled sequencing of up to 92 individuals with as little as 70 ng of DNA per person. Modified solid phase reversible immobilization bead purification strategies enable no sample transfers from sonication in 96-well plates through adapter ligation, resulting in 50% less library preparation reagent consumption. Custom Y-shaped adapters containing novel 7 base pair index sequences with a Hamming distance of ≥2 were directly ligated onto fragmented source DNA eliminating the need for PCR to incorporate indexes, and was followed by a custom blocking strategy using a single oligonucleotide regardless of index sequence. These results were obtained aligning raw reads against the entire genome using Novoalign followed by variant calling of non-indexed pools using SPLINTER or SAMtools for indexed samples. With these pipelines, we find sensitivity and specificity of 99.4% and 99.7% for pooled exome sequencing. Sensitivity, and to a lesser degree specificity, proved to be a function of coverage. For rare variants (≤2% minor allele frequency), we achieved sensitivity and specificity of ≥94.9% and ≥99.99% for custom capture of 2.5 Mb in multiplexed libraries of 22-48 individuals with only ≥5-fold coverage/chromosome, but these parameters improved to ≥98.7 and 100% with 20-fold coverage/chromosome. CONCLUSIONS: This highly scalable methodology enables accurate rare variant detection, with or without individual DNA sample indexing, while reducing the amount of required source DNA and total costs through less hybridization reagent consumption, multi-sample sonication in a standard PCR plate, multiplexed pre-enrichment pooling with a single hybridization and lesser sequencing coverage required to obtain high sensitivity.

Multiple oral Candida infections in patients with Sjogren's syndrome -- prevalence and clinical and drug susceptibility profiles.

OBJECTIVE: To determine the prevalence of oral candidiasis and multiple oral Candida infections in patients with primary Sjögren's syndrome (SS), and the clinical and drug susceptibility profile. METHODS: Thirty patients with primary SS were enrolled in our study. The diagnosis of oral candidiasis was based on the clinical manifestation, and confirmed by a concentrated rinse culture. Candida spp. assessment was accomplished using standard methods: Sabouraud dextrose agar with 50 mg/l chloramphenicol and CHROMagar were used for the rapid screening of clinical species, followed by the API 20C system for further species identification. In vitro antifungal drug susceptibility of Candida isolates was determined by the minimal inhibitory concentrations. RESULTS: In our study, 87% (26/30) of subjects had oral candidiasis, in which 42% (11/26) had multiple Candida spp. infection. Although C. albicans remains the predominant isolate, other rare species such as C. tropicalis, C. glabrata, C. parapsilosis, and C. krusei were present, alone or in combination. Chronic atrophic candidiasis is the most common clinical type of oral candidiasis in patients with SS. The susceptibilities of the 44 Candida isolates to 7 antifungal agents varied dramatically. The resistance to azoles was remarkable, and the phenomenon of cross-resistance between itraconazole and fluconazole was observed. CONCLUSION: Patients with primary SS carry a high risk of oral candidiasis and a high frequency of multiple Candida infections. The azole resistance patterns of Candida spp. support the necessity for drug susceptibility testing as a routine procedure for patients with oral Candida infections.

A new strategy for genome assembly using short sequence reads and reduced representation libraries.

We have developed a novel approach for using massively parallel short-read sequencing to generate fast and inexpensive de novo genomic assemblies comparable to those generated by capillary-based methods. The ultrashort (<100 base) sequences generated by this technology pose specific biological and computational challenges for de novo assembly of large genomes. To account for this, we devised a method for experimentally partitioning the genome using reduced representation (RR) libraries prior to assembly. We use two restriction enzymes independently to create a series of overlapping fragment libraries, each containing a tractable subset of the genome. Together, these libraries allow us to reassemble the entire genome without the need of a reference sequence. As proof of concept, we applied this approach to sequence and assembled the majority of the 125-Mb Drosophila melanogaster genome. We subsequently demonstrate the accuracy of our assembly method with meaningful comparisons against the current available D. melanogaster reference genome (dm3). The ease of assembly and accuracy for comparative genomics suggest that our approach will scale to future mammalian genome-sequencing efforts, saving both time and money without sacrificing quality.

Use of masseteric and deep temporal nerve blocks for reduction of mandibular dislocation.

A patient presented with a unilateral dislocated condyle that was resistant to reduction by simple manual manipulation because of elevator muscle spasm and severe muscle and temporomandibular joint pain. A technique involving a masseteric nerve block and a temporal nerve block was used, allowing a quick, safe, and minimally painful reduction. The method used for delivering these nerve blocks is described here.