ArCH: improving the performance of clonal hematopoiesis variant calling and interpretation.

MOTIVATION: The acquisition of somatic mutations in hematopoietic stem and progenitor stem cells with resultant clonal expansion, termed clonal hematopoiesis (CH), is associated with increased risk of hematologic malignancies and other adverse outcomes. CH is generally present at low allelic fractions, but clonal expansion and acquisition of additional mutations leads to hematologic cancers in a small proportion of individuals. With high depth and high sensitivity sequencing, CH can be detected in most adults and its clonal trajectory mapped over time. However, accurate CH variant calling is challenging due to the difficulty in distinguishing low frequency CH mutations from sequencing artifacts. The lack of well-validated bioinformatic pipelines for CH calling may contribute to lack of reproducibility in studies of CH. RESULTS: Here, we developed ArCH, an Artifact filtering Clonal Hematopoiesis variant calling pipeline for detecting single nucleotide variants and short insertions/deletions by combining the output of four variant calling tools and filtering based on variant characteristics and sequencing error rate estimation. ArCH is an end-to-end cloud-based pipeline optimized to accept a variety of inputs with customizable parameters adaptable to multiple sequencing technologies, research questions, and datasets. Using deep targeted sequencing data generated from six acute myeloid leukemia patient tumor: normal dilutions, 31 blood samples with orthogonal validation, and 26 blood samples with technical replicates, we show that ArCH improves the sensitivity and positive predictive value of CH variant detection at low allele frequencies compared to standard application of commonly used variant calling approaches. AVAILABILITY AND IMPLEMENTATION: The code for this workflow is available at: https://github.com/kbolton-lab/ArCH.

A study to assess the efficacy of enasidenib and risk-adapted addition of azacitidine in newly diagnosed IDH2-mutant AML.

ABSTRACT: Enasidenib (ENA) is an inhibitor of isocitrate dehydrogenase 2 (IDH2) approved for the treatment of patients with IDH2-mutant relapsed/refractory acute myeloid leukemia (AML). In this phase 2/1b Beat AML substudy, we applied a risk-adapted approach to assess the efficacy of ENA monotherapy for patients aged ≥60 years with newly diagnosed IDH2-mutant AML in whom genomic profiling demonstrated that mutant IDH2 was in the dominant leukemic clone. Patients for whom ENA monotherapy did not induce a complete remission (CR) or CR with incomplete blood count recovery (CRi) enrolled in a phase 1b cohort with the addition of azacitidine. The phase 2 portion assessing the overall response to ENA alone demonstrated efficacy, with a composite complete response (cCR) rate (CR/CRi) of 46% in 60 evaluable patients. Seventeen patients subsequently transitioned to phase 1b combination therapy, with a cCR rate of 41% and 1 dose-limiting toxicity. Correlative studies highlight mechanisms of clonal elimination with differentiation therapy as well as therapeutic resistance. This study demonstrates both efficacy of ENA monotherapy in the upfront setting and feasibility and applicability of a risk-adapted approach to the upfront treatment of IDH2-mutant AML. This trial is registered at www.clinicaltrials.gov as #NCT03013998.

Clonal Hematopoiesis in Clinical and Experimental Heart Failure With Preserved Ejection Fraction.

BACKGROUND: Clonal hematopoiesis (CH), which results from an array of nonmalignant driver gene mutations, can lead to altered immune cell function and chronic disease, and has been associated with worse outcomes in patients with heart failure (HF) with reduced ejection fraction. However, the role of CH in the prognosis of HF with preserved ejection fraction (HFpEF) has been understudied. This study aimed to characterize CH in patients with HFpEF and elucidate its causal role in a murine model. METHODS: Using a panel of 20 candidate CH driver genes and a variant allele fraction cutoff of 0.5%, ultradeep error-corrected sequencing identified CH in a cohort of 81 patients with HFpEF (mean age, 71±6 years; ejection fraction, 63±5%) and 36 controls without a diagnosis of HFpEF (mean age, 74±7 years; ejection fraction, 61.5±8%). CH was also evaluated in a replication cohort of 59 individuals with HFpEF. RESULTS: Compared with controls, there was an enrichment of TET2-mediated CH in the HFpEF patient cohort (12% versus 0%, respectively; P=0.02). In the HFpEF cohort, patients with CH exhibited exacerbated diastolic dysfunction in terms of E/e' (14.9 versus 11.7, respectively; P=0.0096) and E/A (1.69 versus 0.89, respectively; P=0.0206) compared with those without CH. The association of CH with exacerbated diastolic dysfunction was corroborated in a validation cohort of individuals with HFpEF. In accordance, patients with HFpEF, an age ≥70 years, and CH exhibited worse prognosis in terms of 5-year cardiovascular-related hospitalization rate (hazard ratio, 5.06; P=0.042) compared with patients with HFpEF and an age ≥70 years without CH. To investigate the causal role of CH in HFpEF, nonconditioned mice underwent adoptive transfer with Tet2-wild-type or Tet2-deficient bone marrow and were subsequently subjected to a high-fat diet/L-NAME (Nω-nitro-l-arginine methyl ester) combination treatment to induce features of HFpEF. This model of Tet2-CH exacerbated cardiac hypertrophy by heart weight/tibia length and cardiomyocyte size, diastolic dysfunction by E/e' and left ventricular end-diastolic pressure, and cardiac fibrosis compared with the Tet2-wild-type condition. CONCLUSIONS: CH is associated with worse heart function and prognosis in patients with HFpEF, and a murine experimental model of Tet2-mediated CH displays greater features of HFpEF.

Pilot Study of Donor-Engrafted Clonal Hematopoiesis Evolution and Clinical Outcomes in Allogeneic Hematopoietic Cell Transplantation Recipients Using a National Registry.

Improved treatment options, such as reduced-intensity conditioning (RIC), enable older patients to receive potentially curative allogeneic hematopoietic cell transplantation (HCT). This progress has led to increased use of older HLA-matched sibling donors. An unintended potential risk associated with older donors is transplantation of donor cells with clonal hematopoiesis (CH) into patients. We aimed to determine the prevalence of CH in older HLA-matched sibling donors pretransplantation and to assess the clinical impact of donor-engrafted CH on HCT outcomes. This was an observational study using donor peripheral blood samples from the Center for International Blood and Marrow Transplant Research repository, linked with corresponding recipient outcomes. To explore engraftment efficiency and evolution of CH mutations following HCT, recipient follow-up samples available through the Bone Marrow Transplant Clinical Trials Network (Protocol 1202) were included. Older donors and patients (both ≥55 years) receiving first RIC HCT for myeloid malignancies were eligible. DNA from archived donor blood samples was used for targeted deep sequencing to identify CH. The associations between donor CH status and recipient outcomes, including acute graft-versus-host disease (aGVHD), chronic GVHD (cGVHD), overall survival, relapse, nonrelapse mortality, disease-free survival, composite GVHD-free and relapse-free survival, and cGVHD-free and relapse-free survival, were analyzed. A total of 299 donors were successfully sequenced to detect CH. At a variant allele frequency (VAF) ≥2%, there were 44 CH mutations in 13.7% (41 of 299) of HLA-matched sibling donors. CH mostly involved DNMT3A (n = 27; 61.4%) and TET2 (n= 9; 20.5%). Post-HCT samples from 13 recipients were also sequenced, of whom 7 had CH+ donors. All of the donor CH mutations (n = 7/7; 100%) were detected in recipients at day 56 or day 90 post-HCT. Overall, mutation VAFs remained relatively constant up to day 90 post-HCT (median change, .005; range, -.008 to .024). Doubling time analysis of recipient day 56 and day 90 data showed that donor-engrafted CH mutations initially expand then decrease to a stable VAF; germline mutations had longer doubling times than CH mutations. The cumulative incidence of grade II-IV aGVHD at day 100 was higher in HCT recipients with CH+ donors (37.5% versus 25.1%); however, the risk for aGVHD by donor CH status did not reach statistical significance (hazard ratio, 1.35; 95% confidence interval, .61 to 3.01; P = .47). There were no statistically significant differences in the cumulative incidence of cGVHD or any secondary outcomes by donor CH status. In subset analysis, the incidence of cGVHD was lower in recipients of grafts from DNMT3A CH+ donors versus donors without DNMT3A CH (34.4% versus 57%; P = .035). Donor cell leukemia was not reported in any donor-recipient pairs. CH in older HLA-matched sibling donors is relatively common and successfully engrafts and persists in recipients. In a homogenous population (myeloid malignancies, older donors and recipients, RICr, non-cyclophosphamide-containing GVHD prophylaxis), we did not detect a difference in cGVHD risk or other secondary outcomes by donor CH status. Subgroup analyses suggest potential differential effects by clinical characteristics and CH mutations. Larger prospective studies are needed to robustly determine which subsets of patients and CH mutations elicit meaningful impacts on clinical outcomes.

KMT2A-D pathogenicity, prevalence, and variation according to a population database.

INTRODUCTION: The KMT2 family of genes is essential epigenetic regulators promoting gene expression. The gene family contains three subgroups, each with two paralogues: KMT2A and KMT2B; KMT2C and KMT2D; KMT2F and KMT2G. KMT2A-D are among the most frequent somatically altered genes in several different cancer types. Somatic KMT2A rearrangements are well-characterized in infant leukemia (IL), and growing evidence supports the role of additional family members (KMT2B, KMT2C, and KMT2D) in leukemogenesis. Enrichment of rare heterozygous frameshift variants in KMT2A and C has been reported in acute myeloid leukemia (AML), IL, and solid tumors. Currently, the non-synonymous variation, prevalence, and penetrance of these four genes are unknown. METHODS: This study determined the prevalence of pathogenic/likely pathogenic (P/LP) germline KMT2A-D variants in a cancer-free adult population from the Genome Aggregation Database (gnomAD). Two methods of variant interpretation were utilized: a manual genomic variant interpretation and an automated ACMG pipeline. RESULTS: The ACMG pipeline identified considerably fewer P/LP variants (n = 89) compared to the manual method (n = 660) in all 4 genes. Consequently, the total P/LP prevalence and allele frequency (AF) were higher in the manual method (1:112, AF = 4.46E-03) than in ACMG (1:832, AF = 6.01E-04). Multiple ancestry-exclusive P/LP variants were identified along with an increased frequency in males compared to females. Many of these variants identified in this population database are also associated with severe juvenile conditions. CONCLUSION: These data demonstrate that putatively functional germline variation in these developmentally important genes is more common than previously appreciated and identification in cancer-free adults may indicate incomplete penetrance for many of these variants. Future research should examine a genetic predisposing role in IL and other pediatric cancers.

Clonal Hematopoiesis Analyses in Clinical, Epidemiologic, and Genetic Aging Studies to Unravel Underlying Mechanisms of Age-Related Dysfunction in Humans.

Aging is characterized by increased mortality, functional decline, and exponential increases in the incidence of diseases such as cancer, stroke, cardiovascular disease, neurological disease, respiratory disease, etc. Though the role of aging in these diseases is widely accepted and considered to be a common denominator, the underlying mechanisms are largely unknown. A significant age-related feature observed in many population cohorts is somatic mosaicism, the detectable accumulation of somatic mutations in multiple cell types and tissues, particularly those with high rates of cell turnover (e.g., skin, liver, and hematopoietic cells). Somatic mosaicism can lead to the development of cellular clones that expand with age in otherwise normal tissues. In the hematopoietic system, this phenomenon has generally been referred to as "clonal hematopoiesis of indeterminate potential" (CHIP) when it applies to a subset of clones in which mutations in driver genes of hematologic malignancies are found. Other mechanisms of clonal hematopoiesis, including large chromosomal alterations, can also give rise to clonal expansion in the absence of conventional CHIP driver gene mutations. Both types of clonal hematopoiesis (CH) have been observed in studies of animal models and humans in association with altered immune responses, increased mortality, and disease risk. Studies in murine models have found that some of these clonal events are involved in abnormal inflammatory and metabolic changes, altered DNA damage repair and epigenetic changes. Studies in long-lived individuals also show the accumulation of somatic mutations, yet at this advanced age, carriership of somatic mutations is no longer associated with an increased risk of mortality. While it remains to be elucidated what factors modify this genotype-phenotype association, i.e., compensatory germline genetics, cellular context of the mutations, protective effects to diseases at exceptional age, it points out that the exceptionally long-lived are key to understand the phenotypic consequences of CHIP mutations. Assessment of the clinical significance of somatic mutations occurring in blood cell types for age-related outcomes in human populations of varied life and health span, environmental exposures, and germline genetic risk factors will be valuable in the development of personalized strategies tailored to specific somatic mutations for healthy aging.

Loss of KMT2C reprograms the epigenomic landscape in hPSCs resulting in NODAL overexpression and a failure of hemogenic endothelium specification.

Germline or somatic variation in the family of KMT2 lysine methyltransferases have been associated with a variety of congenital disorders and cancers. Notably, KMT2A-fusions are prevalent in 70% of infant leukaemias but fail to phenocopy short latency leukaemogenesis in mammalian models, suggesting additional factors are necessary for transformation. Given the lack of additional somatic mutation, the role of epigenetic regulation in cell specification, and our prior results of germline KMT2C variation in infant leukaemia patients, we hypothesized that germline dysfunction of KMT2C altered haematopoietic specification. In isogenic KMT2C KO hPSCs, we found genome-wide differences in histone modifications at active and poised enhancers, leading to gene expression profiles akin to mesendoderm rather than mesoderm highlighted by a significant increase in NODAL expression and WNT inhibition, ultimately resulting in a lack of in vitro hemogenic endothelium specification. These unbiased multi-omic results provide new evidence for germline mechanisms increasing risk of early leukaemogenesis.

High-dose AraC is essential for the treatment of ML-DS independent of postinduction MRD: results of the COG AAML1531 trial.

Myeloid leukemia in children with Down syndrome (ML-DS) is associated with young age and somatic GATA1 mutations. Because of high event-free survival (EFS) and hypersensitivity of the leukemic blasts to chemotherapy, the prior Children's Oncology Group protocol ML-DS protocol (AAML0431) reduced overall treatment intensity but lacking risk stratification, retained the high-dose cytarabine course (HD-AraC), which was highly associated with infectious morbidity. Despite high EFS of ML-DS, survival for those who relapse is rare. AAML1531 introduced therapeutic risk stratification based on the previously identified prognostic factor, measurable residual disease (MRD) at the end of the first induction course. Standard risk (SR) patients were identified by negative MRD using flow cytometry (<0.05%) and did not receive the historically administered HD-AraC course. Interim analysis of 114 SR patients revealed a 2-year EFS of 85.6% (95% confidence interval [CI], 75.7-95.5), which was significantly lower than for MRD- patients treated with HD-AraC on AAML0431 (P = .0002). Overall survival at 2 years was 91.0% (95% CI, 83.8-95.0). Twelve SR patients relapsed, mostly within 1 year from study entry and had a 1-year OS of 16.7% (95% CI, 2.7-41.3). Complex karyotypes were more frequent in SR patients who relapsed compared with those who did not (36% vs 9%; P = .0248). MRD by error-corrected sequencing of GATA1 mutations was piloted in 18 SR patients and detectable in 60% who relapsed vs 23% who did not (P = .2682). Patients with SR ML-DS had worse outcomes without HD-AraC after risk classification based on flow cytometric MRD.

A Non-stop identity complex (NIC) supervises enterocyte identity and protects from premature aging.

A hallmark of aging is loss of differentiated cell identity. Aged Drosophila midgut differentiated enterocytes (ECs) lose their identity, impairing tissue homeostasis. To discover identity regulators, we performed an RNAi screen targeting ubiquitin-related genes in ECs. Seventeen genes were identified, including the deubiquitinase Non-stop (CG4166). Lineage tracing established that acute loss of Non-stop in young ECs phenocopies aged ECs at cellular and tissue levels. Proteomic analysis unveiled that Non-stop maintains identity as part of a Non-stop identity complex (NIC) containing E(y)2, Sgf11, Cp190, (Mod) mdg4, and Nup98. Non-stop ensured chromatin accessibility, maintaining the EC-gene signature, and protected NIC subunit stability. Upon aging, the levels of Non-stop and NIC subunits declined, distorting the unique organization of the EC nucleus. Maintaining youthful levels of Non-stop in wildtype aged ECs safeguards NIC subunits, nuclear organization, and suppressed aging phenotypes. Thus, Non-stop and NIC, supervise EC identity and protects from premature aging.

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.

Humanized zebrafish enhance human hematopoietic stem cell survival and promote acute myeloid leukemia clonal diversity.

Xenograft models are invaluable tools in establishing the current paradigms of hematopoiesis and leukemogenesis. The zebrafish has emerged as a robust alternative xenograft model but, like mice, lack specific cytokines that mimic the microenvironment found in human patients. To address this critical gap, we generated the first humanized zebrafish that express human hematopoietic-specific cytokines (GM-CSF, SCF, and SDF1α). Termed GSS fish, these zebrafish promote survival, self-renewal and multilineage differentiation of human hematopoietic stem and progenitor cells and result in enhanced proliferation and hematopoietic niche-specific homing of primary human leukemia cells. Using error-corrected RNA sequencing, we determined that patient-derived leukemias transplanted into GSS zebrafish exhibit broader clonal representation compared to transplants into control hosts. GSS zebrafish incorporating error-corrected RNA sequencing establish a new standard for zebrafish xenotransplantation that more accurately recapitulates the human context, providing a more representative cost-effective preclinical model system for evaluating personalized response-based treatment in leukemia and therapies to expand human hematopoietic stem and progenitor cells in the transplant setting.

Germline Sequencing Identifies Rare Variants in Finnish Subjects with Familial Germ Cell Tumors.

PURPOSE: Pediatric germ cell tumors are rare, representing about 3% of childhood malignancies in children less than 15 years of age, presenting in neonates or adolescents with a greater incidence noted in older adolescents. Aberrations in primordial germ cell proliferation/differentiation can lead to a variety of neoplasms, including teratomas, embryonal carcinoma, choriocarcinoma, and yolk sac tumors. PATIENTS AND METHODS: Three Finnish families with varying familial germ cell tumors were identified, and whole-genome sequencing was performed using an Illumina sequencing platform. In total, 22 unique subjects across the three families were sequenced. Family 1 proband (female) was affected by malignant ovarian teratoma, Family 2 proband (female) was affected by sacrococcygeal teratoma with yolk sac tumor in the setting of Cornelia de Lange syndrome, and Family 3 proband (male) was affected by malignant testicular teratoma. Rare variants were identified using an autosomal recessive or de novo model of inheritance. RESULTS: For family 1 proband (female), an autosomal recessive or de novo model of inheritance identified variants of interest in the following genes: CD109, IKBKB, and CTNNA3, SUPT6H, MUC5AC, and FRG1. Family 2 proband (female) analysis identified gene variants of interest in the following genes: LONRF2, ANO7, HS6ST1, PRB2, and DNM2. Family 3 proband (male) analysis identified the following potential genes: CRIPAK, KRTAP5-7, and CACNA1B. CONCLUSION: Leveraging deep pedigrees and next-generation sequencing, rare germline variants were identified that were enriched in three families from Finland with a history of familial germ cell tumors. The data presented support the importance of germline mutations when analyzing complex cancers with a low somatic mutation landscape.

Strong concordance between RNA structural and single nucleotide variants identified via next generation sequencing techniques in primary pediatric leukemia and patient-derived xenograft samples.

Acute leukemia represents the most common pediatric malignancy comprising diverse subtypes with varying prognosis and treatment outcomes. New and targeted treatment options are warranted for this disease. Patient-derived xenograft (PDX) models are increasingly being used for preclinical testing of novel treatment modalities. A novel approach involving targeted error-corrected RNA sequencing using ArcherDX HemeV2 kit was employed to compare 25 primary pediatric acute leukemia samples and their corresponding PDX samples. A comparison of the primary samples and PDX samples revealed a high concordance between single nucleotide variants and gene fusions whereas other complex structural variants were not as consistent. The presence of gene fusions representing the major driver mutations at similar allelic frequencies in PDX samples compared to primary samples and over multiple passages confirms the utility of PDX models for preclinical drug testing. Characterization and tracking of these novel cryptic fusions and exonal variants in PDX models is critical in assessing response to potential new therapies.

Prevalence of clinically actionable disease variants in exceptionally long-lived families.

BACKGROUND: Phenotypic expression of pathogenic variants in individuals with no family history of inherited disorders remains unclear. METHODS: We evaluated the prevalence of pathogenic variants in 25 genes associated with Mendelian-inherited disorders in 3015 participants from 485 families in the Long Life Family Study (LLFS). Boot-strapping and Fisher's exact test were used to determine whether allele frequencies in LLFS were significantly different from the allele frequencies reported in publicly available genomic databases. RESULTS: The proportions of pathogenic autosomal dominant mutation carriers in BRCA1 and SDHC in LLFS study participants were similar to those reported in publicly available genomic databases (0.03% vs. 0.0008%, p = 1 for BRCA1, and 0.08% vs. 0.003%, p = 0.05 for SDHC). The frequency of carriers of pathogenic autosomal recessive variants in CPT2, ACADM, SUMF1, WRN, ATM, and ACADVL were also similar in LLFS as compared to those reported in genomic databases. The lack of clinical disease among LLFS participants with well-established pathogenic variants in BRCA1 and SDHC suggests that penetrance of pathogenic variants may be different in long lived families. CONCLUSION: Further research is needed to better understand the penetrance of pathogenic variants before expanding large scale genomic testing to asymptomatic individuals.

Error-corrected sequencing strategies enable comprehensive detection of leukemic mutations relevant for diagnosis and minimal residual disease monitoring.

BACKGROUND: Pediatric leukemias have a diverse genomic landscape associated with complex structural variants, including gene fusions, insertions and deletions, and single nucleotide variants. Routine karyotype and fluorescence in situ hybridization (FISH) techniques lack sensitivity for smaller genomic alternations. Next-generation sequencing (NGS) assays are being increasingly utilized for assessment of these various lesions. However, standard NGS lacks quantitative sensitivity for minimal residual disease (MRD) surveillance due to an inherently high error rate. METHODS: Primary bone marrow samples from pediatric leukemia (n = 32) and adult leukemia subjects (n = 5), cell line MV4-11, and an umbilical cord sample were utilized for this study. Samples were sequenced using molecular barcoding with targeted DNA and RNA library enrichment techniques based on anchored multiplexed PCR (AMP®) technology, amplicon based error-corrected sequencing (ECS) or a human cancer transcriptome assay. Computational analyses were performed to quantitatively assess limit of detection (LOD) for various DNA and RNA lesions, which could be systematically used for MRD assays. RESULTS: Matched leukemia patient samples were analyzed at three time points; diagnosis, end of induction (EOI), and relapse. Similar to flow cytometry for ALL MRD, the LOD for point mutations by these sequencing strategies was ≥0.001. For DNA structural variants, FLT3 internal tandem duplication (ITD) positive cell line and patient samples showed a LOD of ≥0.001 in addition to previously unknown copy number losses in leukemia genes. ECS in RNA identified multiple novel gene fusions, including a SPANT-ABL gene fusion in an ALL patient, which could have been used to alter therapy. Collectively, ECS for RNA demonstrated a quantitative and complex landscape of RNA molecules with 12% of the molecules representing gene fusions, 12% exon duplications, 8% exon deletions, and 68% with retained introns. Droplet digital PCR validation of ECS-RNA confirmed results to single mRNA molecule quantities. CONCLUSIONS: Collectively, these assays enable a highly sensitive, comprehensive, and simultaneous analysis of various clonal leukemic mutations, which can be tracked across disease states (diagnosis, EOI, and relapse) with a high degree of sensitivity. The approaches and results presented here highlight the ability to use NGS for MRD tracking.

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.

Engraftment of rare, pathogenic donor hematopoietic mutations in unrelated hematopoietic stem cell transplantation.

Clonal hematopoiesis is associated with various age-related morbidities. Error-corrected sequencing (ECS) of human blood samples, with a limit of detection of ≥0.0001, has demonstrated that nearly every healthy individual >50 years old harbors rare hematopoietic clones below the detection limit of standard high-throughput sequencing. If these rare mutations confer survival or proliferation advantages, then the clone(s) could expand after a selective pressure such as chemotherapy, radiotherapy, or chronic immunosuppression. Given these observations and the lack of quantitative data regarding clonal hematopoiesis in adolescents and young adults, who are more likely to serve as unrelated hematopoietic stem cell donors, we completed this pilot study to determine whether younger adults harbored hematopoietic clones with pathogenic mutations, how often those clones were transferred to recipients, and what happened to these clones over time after transplantation. We performed ECS on 125 blood and marrow samples from 25 matched unrelated donors and recipients. Clonal mutations, with a median variant allele frequency of 0.00247, were found in 11 donors (44%; median, 36 years old). Of the mutated clones, 84.2% of mutations were predicted to be molecularly pathogenic and 100% engrafted in recipients. Recipients also demonstrated de novo clonal expansion within the first 100 days after hematopoietic stem cell transplant (HSCT). Given this pilot demonstration that rare, pathogenic clonal mutations are far more prevalent in younger adults than previously appreciated, and they engraft in recipients and persist over time, larger studies with longer follow-up are necessary to correlate clonal engraftment with post-HSCT morbidity.

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.