Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Cost-Effective and Scalable Clonal Hematopoiesis Assay Provides Insight into Clonal Dynamics.

Clonal hematopoiesis of indeterminate potential (CHIP) is a common age-related phenomenon in which hematopoietic stem cells acquire mutations in a select set of genes commonly mutated in myeloid neoplasia which then expand clonally. Current sequencing assays to detect CHIP mutations are not optimized for the detection of these variants and can be cost-prohibitive when applied to large cohorts or to serial sequencing. In this study, an affordable (approximately US $8 per sample), accurate, and scalable sequencing assay for CHIP is introduced and validated. The efficacy of the assay was demonstrated by identifying CHIP mutations in a cohort of 456 individuals with DNA collected at multiple time points in Vanderbilt University's biobank and quantifying clonal expansion rates over time. A total of 101 individuals with CHIP/clonal cytopenia of undetermined significance were identified, and individual-level clonal expansion rate was calculated using the variant allele fraction at both time points. Differences in clonal expansion rate by driver gene were observed, but there was also significant individual-level heterogeneity, emphasizing the multifactorial nature of clonal expansion. Additionally, mutation co-occurrence and clonal competition between multiple driver mutations were explored.

Solid Organ Transplant Recipients Exhibit More TET2-Mutant Clonal Hematopoiesis of Indeterminate Potential Not Driven by Increased Transplantation Risk.

PURPOSE: Solid organ transplant recipients comprise a unique population of immunosuppressed patients with increased risk of malignancy, including hematologic neoplasms. Clonal hematopoiesis of indeterminate potential (CHIP) represents a known risk factor for hematologic malignancy and this study describes the prevalence and patterns of CHIP mutations across several types of solid organ transplants. EXPERIMENTAL DESIGN: We use two national biobank cohorts comprised of >650,000 participants with linked genomic and longitudinal phenotypic data to describe the features of CHIP across 2,610 individuals who received kidney, liver, heart, or lung allografts. RESULTS: We find individuals with an allograft before their biobank enrollment had an increased prevalence of TET2 mutations (OR, 1.90; P = 4.0e-4), but individuals who received transplants post-enrollment had a CHIP mutation spectrum similar to that of the general population, without enrichment of TET2. In addition, we do not observe an association between CHIP and risk of incident transplantation among the overall population (HR, 1.02; P = 0.91). And in an exploratory analysis, we do not find evidence for a strong association between CHIP and rates of transplant complications such as rejection or graft failure. CONCLUSIONS: These results demonstrate that recipients of solid organ transplants display a unique pattern of clonal hematopoiesis with enrichment of TET2 driver mutations, the causes of which remain unclear and are deserving of further study.

Driver mutation zygosity is a critical factor in predicting clonal hematopoiesis transformation risk.

Clonal hematopoiesis (CH) can be caused by either single gene mutations (eg point mutations in JAK2 causing CHIP) or mosaic chromosomal alterations (e.g., loss of heterozygosity at chromosome 9p). CH is associated with a significantly increased risk of hematologic malignancies. However, the absolute rate of transformation on an annualized basis is low. Improved prognostication of transformation risk is urgently needed for routine clinical practice. We hypothesized that the co-occurrence of CHIP and mCAs at the same locus (e.g., transforming a heterozygous JAK2 CHIP mutation into a homozygous mutation through concomitant loss of heterozygosity at chromosome 9) might have important prognostic implications for malignancy transformation risk. We tested this hypothesis using our discovery cohort, the UK Biobank (n = 451,180), and subsequently validated it in the BioVU cohort (n = 91,335). We find that individuals with a concurrent somatic mutation and mCA were at significantly increased risk of hematologic malignancy (for example, In BioVU cohort incidence of hematologic malignancies is higher in individuals with co-occurring JAK2 V617F and 9p CN-LOH; HR = 54.76, 95% CI = 33.92-88.41, P < 0.001 vs. JAK2 V617F alone; HR = 44.05, 95% CI = 35.06-55.35, P < 0.001). Currently, the 'zygosity' of the CHIP mutation is not routinely reported in clinical assays or considered in prognosticating CHIP transformation risk. Based on these observations, we propose that clinical reports should include 'zygosity' status of CHIP mutations and that future prognostication systems should take mutation 'zygosity' into account.

Cost-effective and scalable clonal hematopoiesis assay provides insight into clonal dynamics.

Clonal hematopoiesis of indeterminate potential (CHIP) is a common age-related phenomenon that occurs when hematopoietic stem cells acquire mutations in a select set of genes commonly mutated in myeloid neoplasia which then expand clonally. Current sequencing assays to detect CHIP are not optimized for the detection of these variants and can be cost-prohibitive when applied to large cohorts or serial sequencing. Here, we present and validate a CHIP targeted sequencing assay that is affordable (∼$8/sample), accurate and highly scalable. To demonstrate the utility of this assay, we detected CHIP in a cohort of 456 individuals with DNA collected at multiple timepoints in the Vanderbilt BioVU biobank and quantified clonal expansion rates over time. A total of 101 individuals with CHIP were identified, and individual-level clonal expansion rate was calculated using the variant allele fraction (VAF) at both timepoints. Differences in clonal expansion rate by driver gene were observed, but there was also significant individual-level heterogeneity, emphasizing the multifactorial nature of clonal expansion. We further describe the mutation co-occurrence and clonal competition between multiple driver mutations.

Genetic Inhibition of APOL1 Pore-Forming Function Prevents APOL1-Mediated Kidney Disease.

SIGNIFICANCE STATEMENT: African Americans are at increased risk of CKD in part due to high-risk (HR) variants in the apolipoprotein L1 ( APOL1 ) gene, termed G1/G2. A different APOL1 variant, p.N264K , reduced the risk of CKD and ESKD among carriers of APOL1 HR variants to levels comparable with individuals with APOL1 low-risk variants in an analysis of 121,492 participants of African ancestry from the Million Veteran Program (MVP). Functional genetic studies in cell models showed that APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR mutations. Pharmacologic inhibitors that mimic this mutation blocking APOL1 -mediated pore formation may be able to prevent and/or treat APOL1 -associated kidney disease. BACKGROUND: African Americans are at increased risk for nondiabetic CKD in part due to HR variants in the APOL1 gene. METHODS: We tested whether a different APOL1 variant, p.N264K , modified the association between APOL1 HR genotypes (two copies of G1/G2) and CKD in a cross-sectional analysis of 121,492 participants of African ancestry from the MVP. We replicated our findings in the Vanderbilt University Biobank ( n =14,386) and National Institutes of Health All of Us ( n =14,704). Primary outcome was CKD and secondary outcome was ESKD among nondiabetic patients. Primary analysis compared APOL1 HR genotypes with and without p.N264K . Secondary analyses included APOL1 low-risk genotypes and tested for interaction. In MVP, we performed sequential logistic regression models adjusting for demographics, comorbidities, medications, and ten principal components of ancestry. Functional genomic studies expressed APOL1 HR variants with and without APOL1 p.N264K in cell models. RESULTS: In the MVP cohort, 15,604 (12.8%) had two APOL1 HR variants, of which 582 (0.5%) also had APOL1 p.N264K . In MVP, 18,831 (15%) had CKD, 4177 (3%) had ESKD, and 34% had diabetes. MVP APOL1 HR, without p.N264K , was associated with increased odds of CKD (odds ratio [OR], 1.72; 95% confidence interval [CI], 1.60 to 1.85) and ESKD (OR, 3.94; 95% CI, 3.52 to 4.41). In MVP, APOL1 p.N264K mitigated the renal risk of APOL1 HR, in CKD (OR, 0.43; 95% CI, 0.28 to 0.65) and ESKD (OR, 0.19; CI 0.07 to 0.51). In the replication cohorts meta-analysis, APOL1 p.N264K mitigated the renal risk of APOL1 HR in CKD (OR, 0.40; 95% CI, 0.18 to 0.92) and ESKD (OR, 0.19; 95% CI, 0.05 to 0.79). In the mechanistic studies, APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR variants. CONCLUSIONS: APOL1 p.N264K is associated with reduced risk of CKD and ESKD among carriers of APOL1 HR to levels comparable with individuals with APOL1 low-risk genotypes.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R 2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified TCL1A , NRIP1 , and TERT locus variants as modulators of mCA clonal expansion rate.

Genetic modification of inflammation- and clonal hematopoiesis-associated cardiovascular risk.

Clonal hematopoiesis of indeterminate potential (CHIP) is associated with an increased risk of cardiovascular diseases (CVDs), putatively via inflammasome activation. We pursued an inflammatory gene modifier scan for CHIP-associated CVD risk among 424,651 UK Biobank participants. We identified CHIP using whole-exome sequencing data of blood DNA and modeled as a composite, considering all driver genes together, as well as separately for common drivers (DNMT3A, TET2, ASXL1, and JAK2). We developed predicted gene expression scores for 26 inflammasome-related genes and assessed how they modify CHIP-associated CVD risk. We identified IL1RAP as a potential key molecule for CHIP-associated CVD risk across genes and increased AIM2 gene expression leading to heightened JAK2- and ASXL1-associated CVD risk. We show that CRISPR-induced Asxl1-mutated murine macrophages had a particularly heightened inflammatory response to AIM2 agonism, associated with an increased DNA damage response, as well as increased IL-10 secretion, mirroring a CVD-protective effect of IL10 expression in ASXL1 CHIP. Our study supports the role of inflammasomes in CHIP-associated CVD and provides evidence to support gene-specific strategies to address CHIP-associated CVD risk.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

A practical approach to curate clonal hematopoiesis of indeterminate potential in human genetic data sets.

Clonal hematopoiesis of indeterminate potential (CHIP) is a common form of age-related somatic mosaicism that is associated with significant morbidity and mortality. CHIP mutations can be identified in peripheral blood samples that are sequenced using approaches that cover the whole genome, the whole exome, or targeted genetic regions; however, differentiating true CHIP mutations from sequencing artifacts and germ line variants is a considerable bioinformatic challenge. We present a stepwise method that combines filtering based on sequencing metrics, variant annotation, and population-based associations to increase the accuracy of CHIP calls. We apply this approach to ascertain CHIP in ∼550 000 individuals in the UK Biobank complete whole exome cohort and the All of Us Research Program initial whole genome release cohort. CHIP ascertainment on this scale unmasks recurrent artifactual variants and highlights the importance of specialized filtering approaches for several genes, including TET2 and ASXL1. We show how small changes in filtering parameters can considerably increase CHIP misclassification and reduce the effect size of epidemiological associations. Our high-fidelity call set refines previous population-based associations of CHIP with incident outcomes. For example, the annualized incidence of myeloid malignancy in individuals with small CHIP clones is 0.03% per year, which increases to 0.5% per year among individuals with very large CHIP clones. We also find a significantly lower prevalence of CHIP in individuals of self-reported Latino or Hispanic ethnicity in All of Us, highlighting the importance of including diverse populations. The standardization of CHIP calling will increase the fidelity of CHIP epidemiological work and is required for clinical CHIP diagnostic assays.

Genetic investigation of the contribution of body composition to anorexia nervosa in an electronic health record setting.

Anorexia nervosa (AN) is a psychiatric disorder defined by anthropometric symptoms, such as low body weight, and cognitive-behavioral symptoms, such as restricted eating, fear of weight gain, and distorted body image. Recent studies have identified a genetic association between AN and metabolic/anthropometric factors, including body mass index (BMI). Although the reported associations may be under pleiotropic genetic influences, they may represent independent risk factors for AN. Here we examined the independent contributions of genetic predisposition to low body weight and polygenic risk (PRS) for AN in a clinical population (Vanderbilt University Medical Center biobank, BioVU). We fitted logistic and linear regression models in a retrospective case-control design (123 AN patients, 615 age-matched controls). We replicated the genetic correlations between PRSBMI and AN (p = 1.12 × 10-3, OR = 0.96), but this correlation disappeared when controlling for lowest BMI (p = 0.84, OR = 1.00). Additionally, we performed a phenome-wide association analysis of the PRSAN and found that the associations with metabolic phenotypes were attenuated when controlling for PRSBMI. These findings suggest that the genetic association between BMI and AN may be a consequence of the weight-related diagnostic criteria for AN and that genetically regulated anthropometric traits (like BMI) may be independent of AN psychopathology. If so, individuals with cognitive-behavioral symptomatology suggestive of AN, but with a higher PRSBMI, may be under-diagnosed given current diagnostic criteria. Furthermore, PRSBMI may serve as an independent risk factor for weight loss and weight gain during recovery.

Astrocytic GABA transporter 1 deficit in novel SLC6A1 variants mediated epilepsy: Connected from protein destabilization to seizures in mice and humans.

OBJECTIVE: Mutations in γ-aminobutyric acid (GABA) transporter 1 (GAT-1)-encoding SLC6A1 have been associated with myoclonic atonic epilepsy and other phenotypes. We determined the patho-mechanisms of the mutant GAT-1, in order to identify treatment targets. METHODS: We conducted whole-exome sequencing of patients with myoclonic atonic epilepsy (MAE) and characterized the seizure phenotypes and EEG patterns. We studied the protein stability and structural changes with homology modeling and machine learning tools. We characterized the function and trafficking of the mutant GAT-1 with 3H radioactive GABA uptake assay and confocal microscopy. We utilized different models including a knockin mouse and human astrocytes derived from induced pluripotent stem cells (iPSCs). We focused on astrocytes because of their direct impact of astrocytic GAT-1 in seizures. RESULTS: We identified four novel SLC6A1 variants associated with MAE and 2 to 4 Hz spike-wave discharges as a common EEG feature. Machine learning tools predicted that the variant proteins are destabilized. The variant protein had reduced expression and reduced GABA uptake due to endoplasmic reticular retention. The consistent observation was made in cortical and thalamic astrocytes from variant-knockin mice and human iPSC-derived astrocytes. The Slc6a+/A288V mouse, representative of MAE, had increased 5-7 Hz spike-wave discharges and absence seizures. INTERPRETATION: SLC6A1 variants in various locations of the protein peptides can cause MAE with similar seizure phenotypes and EEG features. Reduced GABA uptake is due to decreased functional GAT-1, which, in thalamic astrocytes, could result in increased extracellular GABA accumulation and enhanced tonic inhibition, leading to seizures and abnormal EEGs.

Genetic mosaicism, intrafamilial phenotypic heterogeneity, and molecular defects of a novel missense SLC6A1 mutation associated with epilepsy and ADHD.

BACKGROUND: Mutations in SLC6A1, encoding γ-aminobutyric acid (GABA) transporter 1 (GAT-1), have been recently associated with a spectrum of neurodevelopmental disorders ranging from variable epilepsy syndromes, intellectual disability (ID), autism and others. To date, most identified mutations are de novo. We here report a pedigree of two siblings associated with myoclonic astatic epilepsy, attention deficit hyperactivity disorder (ADHD), and ID. METHODS: Next-generation sequencing identified a missense mutation in the SLC6A1 gene (c.373G > A(p.Val125Met)) in the sisters but not in their shared mother who is also asymptomatic, suggesting gonadal mosaicism. We have thoroughly characterized the clinical phenotypes: EEG recordings identified features for absence seizures and prominent bursts of occipital intermittent rhythmic delta activity (OIRDA). The molecular pathophysiology underlying the clinical phenotypes was assessed using a multidisciplinary approach including machine learning, confocal microscopy, and high-throughput 3H radio-labeled GABA uptake assays in mouse astrocytes and neurons. RESULTS: The GAT-1(Val125Met) mutation destabilizes the global protein conformation and reduces transporter protein expression at total and cell surface. The mutant transporter protein was localized intracellularly inside the endoplasmic reticulum (ER) in both HEK293T cells and astrocytes which may directly contribute to seizures in patients. Radioactive 3H-labeled GABA uptake assay indicated the mutation reduced the function of the mutant GAT-1(Val125Met) to ~30% of the wildtype. CONCLUSIONS: The seizure phenotypes, ADHD, and impaired cognition are likely caused by a partial loss-of-function of GAT-1 due to protein destabilization resulting from the mutation. Reduced GAT-1 function in astrocytes and neurons may consequently alter brain network activities such as increased seizures and reduced attention.