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.

FGF20 and PGM2 variants are associated with childhood asthma in family-based whole-genome sequencing studies.

BACKGROUND: Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations. METHODS: We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants. RESULTS: Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22. CONCLUSION: Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

A statistical physics approach for disease module detection.

Extensive evidence indicates that the pathobiological processes of a complex disease are associated with perturbation in specific neighborhoods of the human protein-protein interaction (PPI) network (also known as the interactome), often referred to as the disease module. Many computational methods have been developed to integrate the interactome and omics profiles to extract context-dependent disease modules. Yet, existing methods all have fundamental limitations in terms of rigor and/or efficiency. Here, we developed a statistical physics approach based on the random-field Ising model (RFIM) for disease module detection, which is both mathematically rigorous and computationally efficient. We applied our RFIM approach to genome-wide association studies (GWAS) of ten complex diseases to examine its performance for disease module detection. We found that our RFIM approach outperforms existing methods in terms of computational efficiency, connectivity of disease modules, and robustness to the interactome incompleteness.

A comparison between similarity matrices for principal component analysis to assess population stratification in sequenced genetic data sets.

Genetic similarity matrices are commonly used to assess population substructure (PS) in genetic studies. Through simulation studies and by the application to whole-genome sequencing (WGS) data, we evaluate the performance of three genetic similarity matrices: the unweighted and weighted Jaccard similarity matrices and the genetic relationship matrix. We describe different scenarios that can create numerical pitfalls and lead to incorrect conclusions in some instances. We consider scenarios in which PS is assessed based on loci that are located across the genome ('globally') and based on loci from a specific genomic region ('locally'). We also compare scenarios in which PS is evaluated based on loci from different minor allele frequency bins: common (>5%), low-frequency (5-0.5%) and rare (<0.5%) single-nucleotide variations (SNVs). Overall, we observe that all approaches provide the best clustering performance when computed based on rare SNVs. The performance of the similarity matrices is very similar for common and low-frequency variants, but for rare variants, the unweighted Jaccard matrix provides preferable clustering features. Based on visual inspection and in terms of standard clustering metrics, its clusters are the densest and the best separated in the principal component analysis of variants with rare SNVs compared with the other methods and different allele frequency cutoffs. In an application, we assessed the role of rare variants on local and global PS, using WGS data from multiethnic Alzheimer's disease data sets and European or East Asian populations from the 1000 Genome Project.

Early Evidence of Chronic Obstructive Pulmonary Disease Obscured by Race-Specific Prediction Equations.

Rationale: The identification of early chronic obstructive pulmonary disease (COPD) is essential to appropriately counsel patients regarding smoking cessation, provide symptomatic treatment, and eventually develop disease-modifying treatments. Disease severity in COPD is defined using race-specific spirometry equations. These may disadvantage non-White individuals in diagnosis and care. Objectives: Determine the impact of race-specific equations on African American (AA) versus non-Hispanic White individuals. Methods: Cross-sectional analyses of the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) cohort were conducted, comparing non-Hispanic White (n = 6,766) and AA (n = 3,366) participants for COPD manifestations. Measurements and Main Results: Spirometric classifications using race-specific, multiethnic, and "race-reversed" prediction equations (NHANES [National Health and Nutrition Examination Survey] and Global Lung Function Initiative "Other" and "Global") were compared, as were respiratory symptoms, 6-minute-walk distance, computed tomography imaging, respiratory exacerbations, and St. George's Respiratory Questionnaire. Application of different prediction equations to the cohort resulted in different classifications by stage, with NHANES and Global Lung Function Initiative race-specific equations being minimally different, but race-reversed equations moving AA participants to more severe stages and especially between the Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage 0 and preserved ratio impaired spirometry groups. Classification using the established NHANES race-specific equations demonstrated that for each of GOLD stages 1-4, AA participants were younger, had fewer pack-years and more current smoking, but had more exacerbations, shorter 6-minute-walk distance, greater dyspnea, and worse BODE (body mass index, airway obstruction, dyspnea, and exercise capacity) scores and St. George's Respiratory Questionnaire scores. Differences were greatest in GOLD stages 1 and 2. Race-reversed equations reclassified 774 AA participants (43%) from GOLD stage 0 to preserved ratio impaired spirometry. Conclusions: Race-specific equations underestimated disease severity among AA participants. These effects were particularly evident in early disease and may result in late detection of COPD.

DRAGON: Determining Regulatory Associations using Graphical models on multi-Omic Networks.

The increasing quantity of multi-omic data, such as methylomic and transcriptomic profiles collected on the same specimen or even on the same cell, provides a unique opportunity to explore the complex interactions that define cell phenotype and govern cellular responses to perturbations. We propose a network approach based on Gaussian Graphical Models (GGMs) that facilitates the joint analysis of paired omics data. This method, called DRAGON (Determining Regulatory Associations using Graphical models on multi-Omic Networks), calibrates its parameters to achieve an optimal trade-off between the network's complexity and estimation accuracy, while explicitly accounting for the characteristics of each of the assessed omics 'layers.' In simulation studies, we show that DRAGON adapts to edge density and feature size differences between omics layers, improving model inference and edge recovery compared to state-of-the-art methods. We further demonstrate in an analysis of joint transcriptome - methylome data from TCGA breast cancer specimens that DRAGON can identify key molecular mechanisms such as gene regulation via promoter methylation. In particular, we identify Transcription Factor AP-2 Beta (TFAP2B) as a potential multi-omic biomarker for basal-type breast cancer. DRAGON is available as open-source code in Python through the Network Zoo package (netZooPy v0.8; netzoo.github.io).

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables.

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

Clonal Somatic Mutations in Chronic Lung Diseases Are Associated with Reduced Lung Function.

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.

Transgenerational Epigenetic Inheritance of Cardiovascular Diseases: A Network Medicine Perspective.

INTRODUCTION: The ability to identify early epigenetic signatures underlying the inheritance of cardiovascular risk, including trans- and intergenerational effects, may help to stratify people before cardiac symptoms occur. METHODS: Prospective and retrospective cohorts and case-control studies focusing on DNA methylation and maternal/paternal effects were searched in Pubmed from 1997 to 2023 by using the following keywords: DNA methylation, genomic imprinting, and network analysis in combination with transgenerational/intergenerational effects. RESULTS: Maternal and paternal exposures to traditional cardiovascular risk factors during critical temporal windows, including the preconceptional period or early pregnancy, may perturb the plasticity of the epigenome (mainly DNA methylation) of the developing fetus especially at imprinted loci, such as the insulin-like growth factor type 2 (IGF2) gene. Thus, the epigenome is akin to a "molecular archive" able to memorize parental environmental insults and predispose an individual to cardiovascular diseases onset in later life. Direct evidence for human transgenerational epigenetic inheritance (at least three generations) of cardiovascular risk is lacking but it is supported by epidemiological studies. Several blood-based association studies showed potential intergenerational epigenetic effects (single-generation studies) which may mediate the transmittance of cardiovascular risk from parents to offspring. DISCUSSION: In this narrative review, we discuss some relevant examples of trans- and intergenerational epigenetic associations with cardiovascular risk. In our perspective, we propose three network-oriented approaches which may help to clarify the unsolved issues regarding transgenerational epigenetic inheritance of cardiovascular risk and provide potential early biomarkers for primary prevention.

Zinc finger protein 33B demonstrates sex interaction with atopy-related markers in childhood asthma.

BACKGROUND: Sex differences related to immune responses can influence atopic manifestations in childhood asthma. While genome-wide association studies have investigated a sex-specific genetic architecture of the immune response, gene-by-sex interactions have not been extensively analysed for atopy-related markers including allergy skin tests, IgE and eosinophils in asthmatic children. METHODS: We performed a genome-wide gene-by-sex interaction analysis for atopy-related markers using whole-genome sequencing data based on 889 trios from the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS) and 284 trios from the Childhood Asthma Management Program (CAMP). We also tested the findings in UK Biobank participants with self-reported childhood asthma. Furthermore, downstream analyses in GACRS integrated gene expression to disentangle observed associations. RESULTS: Single nucleotide polymorphism (SNP) rs1255383 at 10q11.21 demonstrated a genome-wide significant gene-by-sex interaction (pinteraction=9.08×10-10) for atopy (positive skin test) with opposite direction of effects between females and males. In the UK Biobank participants with a history of childhood asthma, the signal was consistently observed with the same sex-specific effect directions for high eosinophil count (pinteraction=0.0058). Gene expression of ZNF33B (zinc finger protein 33B), located at 10q11.21, was moderately associated with atopy in girls, but not in boys. CONCLUSIONS: We report SNPs in/near a zinc finger gene as novel sex-differential loci for atopy-related markers with opposite effect directions in females and males. A potential role for ZNF33B should be studied further as an important driver of sex-divergent features of atopy in childhood asthma.

Epigenome-Wide Analysis of DNA Methylation and Optimism in Women and Men.

OBJECTIVE: Higher optimism is associated with reduced mortality and a lower risk of age-related chronic diseases. DNA methylation (DNAm) may provide insight into mechanisms underlying these relationships. We hypothesized that DNAm would differ among older individuals who are more versus less optimistic. METHODS: Using cross-sectional data from two population-based cohorts of women with diverse races/ethnicities ( n = 3816) and men (only White, n = 667), we investigated the associations of optimism with epigenome-wide leukocyte DNAm. Random-effects meta-analyses were subsequently used to pool the individual results. Significantly differentially methylated cytosine-phosphate-guanines (CpGs) were identified by the "number of independent degrees of freedom" approach: effective degrees of freedom correction using the number of principal components (PCs), explaining >95% of the variation of the DNAm data (PC-correction). We performed regional analyses using comb-p and pathway analyses using the Ingenuity Pathway Analysis software. RESULTS: We found that essentially all CpGs (total probe N = 359,862) were homogeneous across sex and race/ethnicity in the DNAm-optimism association. In the single CpG site analyses based on homogeneous CpGs, we identified 13 significantly differentially methylated probes using PC-correction. We found four significantly differentially methylated regions and two significantly differentially methylated pathways. The annotated genes from the single CpG site and regional analyses are involved in psychiatric disorders, cardiovascular disease, cognitive impairment, and cancer. Identified pathways were related to cancer, and neurodevelopmental and neurodegenerative disorders. CONCLUSION: Our findings provide new insights into possible mechanisms underlying optimism and health.

Benchmarking omics-based prediction of asthma development in children.

BACKGROUND: Asthma is a heterogeneous disease with high morbidity. Advancement in high-throughput multi-omics approaches has enabled the collection of molecular assessments at different layers, providing a complementary perspective of complex diseases. Numerous computational methods have been developed for the omics-based patient classification or disease outcome prediction. Yet, a systematic benchmarking of those methods using various combinations of omics data for the prediction of asthma development is still lacking. OBJECTIVE: We aimed to investigate the computational methods in disease status prediction using multi-omics data. METHOD: We systematically benchmarked 18 computational methods using all the 63 combinations of six omics data (GWAS, miRNA, mRNA, microbiome, metabolome, DNA methylation) collected in The Vitamin D Antenatal Asthma Reduction Trial (VDAART) cohort. We evaluated each method using standard performance metrics for each of the 63 omics combinations. RESULTS: Our results indicate that overall Logistic Regression, Multi-Layer Perceptron, and MOGONET display superior performance, and the combination of transcriptional, genomic and microbiome data achieves the best prediction. Moreover, we find that including the clinical data can further improve the prediction performance for some but not all the omics combinations. CONCLUSIONS: Specific omics combinations can reach the optimal prediction of asthma development in children. And certain computational methods showed superior performance than other methods.

X chromosome associations with chronic obstructive pulmonary disease and related phenotypes: an X chromosome-wide association study.

BACKGROUND: The association between genetic variants on the X chromosome to risk of COPD has not been fully explored. We hypothesize that the X chromosome harbors variants important in determining risk of COPD related phenotypes and may drive sex differences in COPD manifestations. METHODS: Using X chromosome data from three COPD-enriched cohorts of adult smokers, we performed X chromosome specific quality control, imputation, and testing for association with COPD case-control status, lung function, and quantitative emphysema. Analyses were performed among all subjects, then stratified by sex, and subsequently combined in meta-analyses. RESULTS: Among 10,193 subjects of non-Hispanic white or European ancestry, a variant near TMSB4X, rs5979771, reached genome-wide significance for association with lung function measured by FEV1/FVC ([Formula: see text] 0.020, SE 0.004, p 4.97 × 10-08), with suggestive evidence of association with FEV1 ([Formula: see text] 0.092, SE 0.018, p 3.40 × 10-07). Sex-stratified analyses revealed X chromosome variants that were differentially trending in one sex, with significantly different effect sizes or directions. CONCLUSIONS: This investigation identified loci influencing lung function, COPD, and emphysema in a comprehensive genetic association meta-analysis of X chromosome genetic markers from multiple COPD-related datasets. Sex differences play an important role in the pathobiology of complex lung disease, including X chromosome variants that demonstrate differential effects by sex and variants that may be relevant through escape from X chromosome inactivation. Comprehensive interrogation of the X chromosome to better understand genetic control of COPD and lung function is important to further understanding of disease pathology. Trial registration Genetic Epidemiology of COPD Study (COPDGene) is registered at ClinicalTrials.gov, NCT00608764 (Active since January 28, 2008). Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints Study (ECLIPSE), GlaxoSmithKline study code SCO104960, is registered at ClinicalTrials.gov, NCT00292552 (Active since February 16, 2006). Genetics of COPD in Norway Study (GenKOLS) holds GlaxoSmithKline study code RES11080, Genetics of Chronic Obstructive Lung Disease.

Use of the Spirometric "Fixed-Ratio" Underdiagnoses COPD in African-Americans in a Longitudinal Cohort Study.

BACKGROUND: COPD diagnosis is tightly linked to the fixed-ratio spirometry criteria of FEV1/FVC < 0.7. African-Americans are less often diagnosed with COPD. OBJECTIVE: Compare COPD diagnosis by fixed-ratio with findings and outcomes by race. DESIGN: Genetic Epidemiology of COPD (COPDGene) (2007-present), cross-sectional comparing non-Hispanic white (NHW) and African-American (AA) participants for COPD diagnosis, manifestations, and outcomes. SETTING: Multicenter, longitudinal US cohort study. PARTICIPANTS: Current or former smokers with ≥ 10-pack-year smoking history enrolled at 21 clinical centers including over-sampling of participants with known COPD and AA. Exclusions were pre-existing non-COPD lung disease, except for a history of asthma. MEASUREMENTS: Subject diagnosis by conventional criteria. Mortality, imaging, respiratory symptoms, function, and socioeconomic characteristics, including area deprivation index (ADI). Matched analysis (age, sex, and smoking status) of AA vs. NHW within participants without diagnosed COPD (GOLD 0; FEV1 ≥ 80% predicted and FEV1/FVC ≥ 0.7). RESULTS: Using the fixed ratio, 70% of AA (n = 3366) were classified as non-COPD, versus 49% of NHW (n = 6766). AA smokers were younger (55 vs. 62 years), more often current smoking (80% vs. 39%), with fewer pack-years but similar 12-year mortality. Density distribution plots for FEV1 and FVC raw spirometry values showed disproportionate reductions in FVC relative to FEV1 in AA that systematically led to higher ratios. The matched analysis demonstrated GOLD 0 AA had greater symptoms, worse DLCO, spirometry, BODE scores (1.03 vs 0.54, p < 0.0001), and greater deprivation than NHW. LIMITATIONS: Lack of an alternative diagnostic metric for comparison. CONCLUSIONS: The fixed-ratio spirometric criteria for COPD underdiagnosed potential COPD in AA participants when compared to broader diagnostic criteria. Disproportionate reductions in FVC relative to FEV1 leading to higher FEV1/FVC were identified in these participants and associated with deprivation. Broader diagnostic criteria for COPD are needed to identify the disease across all populations.

Alpha-1 Antitrypsin MZ Heterozygosity Is an Endotype of Chronic Obstructive Pulmonary Disease.

Rationale: Multiple studies have demonstrated an increased risk of chronic obstructive pulmonary disease (COPD) in heterozygous carriers of the AAT (alpha-1 antitrypsin) Z allele. However, it is not known if MZ subjects with COPD are phenotypically different from noncarriers (MM genotype) with COPD. Objectives: To assess if MZ subjects with COPD have different clinical features compared with MM subjects with COPD. Methods: Genotypes of SERPINA1 were ascertained by using whole-genome sequencing data in three independent studies. We compared outcomes between MM subjects with COPD and MZ subjects with COPD in each study and combined the results in a meta-analysis. We performed longitudinal and survival analyses to compare outcomes in MM and MZ subjects with COPD over time. Measurements and Main Results: We included 290 MZ subjects with COPD and 6,184 MM subjects with COPD across the three studies. MZ subjects had a lower FEV1% predicted and greater quantitative emphysema on chest computed tomography scans compared with MM subjects. In a meta-analysis, the FEV1 was 3.9% lower (95% confidence interval [CI], -6.55% to -1.26%) and emphysema (the percentage of lung attenuation areas <-950 HU) was 4.14% greater (95% CI, 1.44% to 6.84%) in MZ subjects. We found one gene, PGF (placental growth factor), to be differentially expressed in lung tissue from one study between MZ subjects and MM subjects. Conclusions: Carriers of the AAT Z allele (those who were MZ heterozygous) with COPD had lower lung function and more emphysema than MM subjects with COPD. Taken with the subtle differences in gene expression between the two groups, our findings suggest that MZ subjects represent an endotype of COPD.

Identifying chronic obstructive pulmonary disease from integrative omics and clustering in lung tissue.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a highly morbid and heterogenous disease. While COPD is defined by spirometry, many COPD characteristics are seen in cigarette smokers with normal spirometry. The extent to which COPD and COPD heterogeneity is captured in omics of lung tissue is not known. METHODS: We clustered gene expression and methylation data in 78 lung tissue samples from former smokers with normal lung function or severe COPD. We applied two integrative omics clustering methods: (1) Similarity Network Fusion (SNF) and (2) Entropy-Based Consensus Clustering (ECC). RESULTS: SNF clusters were not significantly different by the percentage of COPD cases (48.8% vs. 68.6%, p = 0.13), though were different according to median forced expiratory volume in one second (FEV1) % predicted (82 vs. 31, p = 0.017). In contrast, the ECC clusters showed stronger evidence of separation by COPD case status (48.2% vs. 81.8%, p = 0.013) and similar stratification by median FEV1% predicted (82 vs. 30.5, p = 0.0059). ECC clusters using both gene expression and methylation were identical to the ECC clustering solution generated using methylation data alone. Both methods selected clusters with differentially expressed transcripts enriched for interleukin signaling and immunoregulatory interactions between lymphoid and non-lymphoid cells. CONCLUSIONS: Unsupervised clustering analysis from integrated gene expression and methylation data in lung tissue resulted in clusters with modest concordance with COPD, though were enriched in pathways potentially contributing to COPD-related pathology and heterogeneity.

Genome-wide association analysis of COVID-19 mortality risk in SARS-CoV-2 genomes identifies mutation in the SARS-CoV-2 spike protein that colocalizes with P.1 of the Brazilian strain.

SARS-CoV-2 mortality has been extensively studied in relation to host susceptibility. How sequence variations in the SARS-CoV-2 genome affect pathogenicity is poorly understood. Starting in October 2020, using the methodology of genome-wide association studies (GWAS), we looked at the association between whole-genome sequencing (WGS) data of the virus and COVID-19 mortality as a potential method of early identification of highly pathogenic strains to target for containment. Although continuously updating our analysis, in December 2020, we analyzed 7548 single-stranded SARS-CoV-2 genomes of COVID-19 patients in the GISAID database and associated variants with mortality using a logistic regression. In total, evaluating 29,891 sequenced loci of the viral genome for association with patient/host mortality, two loci, at 12,053 and 25,088 bp, achieved genome-wide significance (p values of 4.09e-09 and 4.41e-23, respectively), though only 25,088 bp remained significant in follow-up analyses. Our association findings were exclusively driven by the samples that were submitted from Brazil (p value of 4.90e-13 for 25,088 bp). The mutation frequency of 25,088 bp in the Brazilian samples on GISAID has rapidly increased from about 0.4 in October/December 2020 to 0.77 in March 2021. Although GWAS methodology is suitable for samples in which mutation frequencies varies between geographical regions, it cannot account for mutation frequencies that change rapidly overtime, rendering a GWAS follow-up analysis of the GISAID samples that have been submitted after December 2020 as invalid. The locus at 25,088 bp is located in the P.1 strain, which later (April 2021) became one of the distinguishing loci (precisely, substitution V1176F) of the Brazilian strain as defined by the Centers for Disease Control. Specifically, the mutations at 25,088 bp occur in the S2 subunit of the SARS-CoV-2 spike protein, which plays a key role in viral entry of target host cells. Since the mutations alter amino acid coding sequences, they potentially imposing structural changes that could enhance viral infectivity and symptom severity. Our analysis suggests that GWAS methodology can provide suitable analysis tools for the real-time detection of new more transmissible and pathogenic viral strains in databases such as GISAID, though new approaches are needed to accommodate rapidly changing mutation frequencies over time, in the presence of simultaneously changing case/control ratios. Improvements of the associated metadata/patient information in terms of quality and availability will also be important to fully utilize the potential of GWAS methodology in this field.

Sex Differences in Airways at Chest CT: Results from the COPDGene Cohort.

Background The prevalence of chronic obstructive pulmonary disease (COPD) in women is fast approaching that in men, and women experience greater symptom burden. Although sex differences in emphysema have been reported, differences in airways have not been systematically characterized. Purpose To evaluate whether structural differences in airways may underlie some of the sex differences in COPD prevalence and clinical outcomes. Materials and Methods In a secondary analyses of a multicenter study of never-, current-, and former-smokers enrolled from January 2008 to June 2011 and followed up longitudinally until November 2020, airway disease on CT images was quantified using seven metrics: airway wall thickness, wall area percent, and square root of the wall thickness of a hypothetical airway with internal perimeter of 10 mm (referred to as Pi10) for airway wall; and lumen diameter, airway volume, total airway count, and airway fractal dimension for airway lumen. Least-squares mean values for each airway metric were calculated and adjusted for age, height, ethnicity, body mass index, pack-years of smoking, current smoking status, total lung capacity, display field of view, and scanner type. In ever-smokers, associations were tested between each airway metric and postbronchodilator forced expiratory volume in 1 second (FEV1)-to-forced vital capacity (FVC) ratio, modified Medical Research Council dyspnea scale, St George's Respiratory Questionnaire score, and 6-minute walk distance. Multivariable Cox proportional hazards models were created to evaluate the sex-specific association between each airway metric and mortality. Results In never-smokers (n = 420), men had thicker airway walls than women as quantified on CT images for segmental airway wall area percentage (least-squares mean, 47.68 ± 0.61 [standard error] vs 45.78 ± 0.55; difference, -1.90; P = .02), whereas airway lumen dimensions were lower in women than men after accounting for height and total lung capacity (segmental lumen diameter, 8.05 mm ± 0.14 vs 9.05 mm ± 0.16; difference, -1.00 mm; P < .001). In ever-smokers (n = 9363), men had greater segmental airway wall area percentage (least-squares mean, 52.19 ± 0.16 vs 48.89 ± 0.18; difference, -3.30; P < .001), whereas women had narrower segmental lumen diameter (7.80 mm ± 0.05 vs 8.69 mm ± 0.04; difference, -0.89; P < .001). A unit change in each of the airway metrics (higher wall or lower lumen measure) resulted in lower FEV1-to-FVC ratio, more dyspnea, poorer respiratory quality of life, lower 6-minute walk distance, and worse survival in women compared with men (all P < .01). Conclusion Airway lumen sizes quantified at chest CT were smaller in women than in men after accounting for height and lung size, and these lower baseline values in women conferred lower reserves against respiratory morbidity and mortality for equivalent changes compared with men. © RSNA, 2022 Online supplemental material is available for this article.

A polygenic risk score and age of diagnosis of COPD.

BACKGROUND: Genetic susceptibility may be associated with earlier onset of chronic obstructive pulmonary disease (COPD). We hypothesised that a polygenic risk score (PRS) for COPD would be associated with earlier age of diagnosis of COPD. METHODS: In 6647 non-Hispanic White (NHW) and 2464 African American (AA) participants from COPDGene, and 6812 participants from the Framingham Heart Study (FHS), we tested the relationship of the PRS and age of COPD diagnosis. Age at diagnosis was determined by: 1) self-reported age at COPD diagnosis or 2) age at visits when moderate-to-severe airflow limitation (Global Initiative for Chronic Obstructive Lung Disease (GOLD) grade 2-4) was observed on spirometry. We used Cox regression to examine the overall and time-dependent effects of the PRS on incident COPD. In the COPDGene study, we also examined the PRS's predictive value for COPD at age <50 years (COPD50) using logistic regression and area under the curve (AUC) analyses, with and without the addition of other risk factors present at early life (e.g. childhood asthma). RESULTS: In Cox models, the PRS demonstrated age-dependent associations with incident COPD, with larger effects at younger ages in both cohorts. The PRS was associated with COPD50 (OR 1.55 (95% CI 1.41-1.71) for NHW, OR 1.23 (95% CI 1.05-1.43) for AA and OR 2.47 (95% CI 2.12-2.88) for FHS participants). In COPDGene, adding the PRS to known early-life risk factors improved prediction of COPD50 in NHW (AUC 0.69 versus 0.74; p<0.0001) and AA (AUC 0.61 versus 0.64; p=0.04) participants. CONCLUSIONS: A COPD PRS is associated with earlier age of diagnosis of COPD and retains predictive value when added to known early-life risk factors.

Interstitial lung abnormalities are associated with decreased mean telomere length.

BACKGROUND: Interstitial lung abnormalities (ILA) share many features with idiopathic pulmonary fibrosis; however, it is not known if ILA are associated with decreased mean telomere length (MTL). METHODS: Telomere length was measured with quantitative PCR in the Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) and Age Gene/Environment Susceptibility Reykjavik (AGES-Reykjavik) cohorts and Southern blot analysis was used in the Framingham Heart Study (FHS). Logistic and linear regression were used to assess the association between ILA and MTL; Cox proportional hazards models were used to assess the association between MTL and mortality. RESULTS: In all three cohorts, ILA were associated with decreased MTL. In the COPDGene and AGES-Reykjavik cohorts, after adjustment there was greater than twofold increase in the odds of ILA when comparing the shortest quartile of telomere length to the longest quartile (OR 2.2, 95% CI 1.5-3.4, p=0.0001, and OR 2.6, 95% CI 1.4-4.9, p=0.003, respectively). In the FHS, those with ILA had shorter telomeres than those without ILA (-767 bp, 95% CI 76-1584 bp, p=0.03). Although decreased MTL was associated with chronic obstructive pulmonary disease (OR 1.3, 95% CI 1.1-1.6, p=0.01) in COPDGene, the effect estimate was less than that noted with ILA. There was no consistent association between MTL and risk of death when comparing the shortest quartile of telomere length in COPDGene and AGES-Reykjavik (HR 0.82, 95% CI 0.4-1.7, p=0.6, and HR 1.2, 95% CI 0.6-2.2, p=0.5, respectively). CONCLUSION: ILA are associated with decreased MTL.