Silverman, Edwin K.

Abstract Rationale Constantly exposed to the external environment and mutagens such as tobacco smoke, human lungs have one of the highest somatic mutation rates among all human organs. However, the relationship of these mutations to lung disease and function is not known. Objectives To identify the prevalence and significance of clonal somatic mutations in chronic lung diseases. Methods We analyzed the clonal somatic mutations from 1,251 samples of normal and diseased noncancerous lung tissue RNA sequencing with paired whole-genome sequencing from the Lung Tissue Research Consortium. We examined the associations of somatic mutations with lung function, disease status, and computationally deconvoluted cell types in two of the most common diseases represented in our dataset, chronic obstructive pulmonary disease (COPD; 29%) and idiopathic pulmonary fibrosis (IPF; 13%). Measurements and Main Results Clonal somatic mutational burden was associated with reduced lung function in both COPD and IPF. We identified an increased prevalence of clonal somatic mutations in individuals with IPF compared with normal control subjects and individuals with COPD independent of age and smoking status. IPF clonal somatic mutations were enriched in disease-related and airway epithelial–expressed genes such as MUC5B in IPF. Patients who were MUC5B risk variant carriers had increased odds of developing somatic mutations of MUC5B that were explained by increased expression of MUC5B. Conclusions Our identification of an increased prevalence of clonal somatic mutation in diseased lung that correlates with airway epithelial gene expression and disease severity highlights for the first time the role of somatic mutational processes in lung disease genetics.

Abstract Rationale In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case–control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10−95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10−87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Abstract Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Sonic hedgehog signaling regulates processes of embryonic development across multiple tissues, yet factors regulating context-specific Shh signaling remain poorly understood. Exome sequencing of families with polymicrogyria (disordered cortical folding) revealed multiple individuals with biallelic deleterious variants in TMEM161B , which encodes a multi-pass transmembrane protein of unknown function. Tmem161b null mice demonstrated holoprosencephaly, craniofacial midline defects, eye defects, and spinal cord patterning changes consistent with impaired Shh signaling, but were without limb defects, suggesting a CNS-specific role of Tmem161b. Tmem161b depletion impaired the response to Smoothened activation in vitro and disrupted cortical histogenesis in vivo in both mouse and ferret models, including leading to abnormal gyration in the ferret model. Tmem161b localizes non-exclusively to the primary cilium, and scanning electron microscopy revealed shortened, dysmorphic, and ballooned ventricular zone cilia in the Tmem161b null mouse, suggesting that the Shh-related phenotypes may reflect ciliary dysfunction. Our data identify TMEM161B as a regulator of cerebral cortical gyration, as involved in primary ciliary structure, as a regulator of Shh signaling, and further implicate Shh signaling in human gyral development.

Abstract With aging, somatic mutations in hematopoietic stem and progenitor cells (HSPC) can give rise to clonal hematopoiesis of indeterminate potential (CHIP), a premalignant state associated with diverse age-related diseases. Here we report the largest multi-ancestry genome-wide analysis of CHIP to date (N = 1,018,305), including individuals of African (N = 85,978), Admixed American (N = 191,371), East Asian (N = 13,532), European (N = 694,015), and South Asian (N = 13,193) ancestry. Multi-ancestry meta-analyses identified 72 genome-wide significant loci, including 44 novel associations implicating genes such as AFF1 , ATF7IP , ATP8B4 , BCL2 , CEBPA , CYRIA , DNM2 , ELF1 , NKX2-3 , PIK3CB , PRDM16 , RPN1 , TERC , and TRIM4 . Notably, variants at MECOM and PHF20L1 showed opposite allelic effects between DNMT3A - and non- DNMT3A -driven CHIP, highlighting driver-specific germline influences. Implicated loci converge on pathways regulating telomere maintenance, cell-cycle control, hematopoietic transcription, DNA damage response and immune signaling. These findings support a model in which germline variation both expands the HSPC pool and biases clonal selection in a driver-dependent manner, providing a mechanistic basis for inter-individual heterogeneity in CHIP. Phenome-wide association analyses further linked CHIP to hematologic, neoplastic, and circulatory traits, with enrichment across hematopoietic and non-hematopoietic cell types. Together, this work expands the genomic and phenomic landscape of CHIP and reveals germline–somatic interactions that shape clonal evolution during aging.

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program ( n = 63,302) and UK Biobank ( n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

ObjectiveTo investigate the independent relationship of rheumatoid arthritis (RA) to the type and severity of pulmonary patterns on spirometry compared to the pulmonary patterns in general population controls.MethodsIn this cross‐sectional study, we investigated the association of RA with pulmonary function measures on spirometry among subjects in the UK Biobank who underwent spirometry for research purposes. RA cases were identified based on self‐report and current disease‐modifying antirheumatic drug/glucocorticoid use. Controls were subjects without RA from the general population. Outcome measures included continuous forced expiratory volume in 1 second percent predicted (FEV1%) and forced vital capacity percent predicted (FVC%), type of spirometric pattern (restrictive or obstructive), and severity of the restrictive or obstructive pattern. We used multivariable regression to estimate the effects in RA cases compared to the effects in controls, adjusting for age, sex, body mass index, and smoking status/pack‐years.ResultsAmong 350,776 analyzed subjects who underwent spirometry (mean age 56.3 years; 55.8% female; 45.5% ever smokers), we identified 2,008 cases of treated RA. In multivariable analyses, RA was associated with lower FEV1% (β = –2.93 [95% confidence interval (95% CI) –3.63, –2.24]), FVC% (β = –2.08 [95% CI –2.72, –1.45]), and FEV1/FVC (β = –0.008 [95% CI –0.010, –0.005]) compared to controls. RA was additionally associated with restrictive patterns (odds ratio [OR] 1.36 [95% CI 1.21, 1.52]) and obstructive patterns (OR 1.21 [95% CI 1.07, 1.37]) independent of confounders, and was most strongly associated with severe restrictive and obstructive patterns.ConclusionRA is associated with increased odds of restrictive and obstructive patterns, and this relationship is not explained by confounders, including smoking status. In addition to restrictive lung disease, clinicians should also be aware that airway obstruction may be a pulmonary manifestation of RA.