Multiethnic meta-analysis identifies ancestry-specific and cross-ancestry loci for pulmonary function

Latin American Trans-ancestry INitiative for OCD genomics (LATINO): Study protocol

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder. Worldwide, its prevalence is ~2% and its etiology is mostly unknown. Identifying biological factors contributing to OCD will elucidate underlying mechanisms and might contribute to improved treatment outcomes. Genomic studies of OCD are beginning to reveal long-sought risk loci, but >95% of the cases currently in analysis are of homogenous European ancestry. If not addressed, this Eurocentric bias will result in OCD genomic findings being more accurate for individuals of European ancestry than other ancestries, thereby contributing to health disparities in potential future applications of genomics. In this study protocol paper, we describe the Latin American Trans-ancestry INitiative for OCD genomics (LATINO, https://www.latinostudy.org). LATINO is a new network of investigators from across Latin America, the United States, and Canada who have begun to collect DNA and clinical data from 5000 richly phenotyped OCD cases of Latin American ancestry in a culturally sensitive and ethical manner. In this project, we will utilize trans-ancestry genomic analyses to accelerate the identification of OCD risk loci, fine-map putative causal variants, and improve the performance of polygenic risk scores in diverse populations. We will also capitalize on rich clinical data to examine the genetics of treatment response, biologically plausible OCD subtypes, and symptom dimensions. Additionally, LATINO will help elucidate the diversity of the clinical presentations of OCD across cultures through various trainings developed and offered in collaboration with Latin American investigators. We believe this study will advance the important goal of global mental health discovery and equity.

GUIDE deconstructs genetic architectures using association studies

Genome-wide association studies have revealed that the genetic architecture of most complex traits is characterized by a large number of distinct effects scattered across the genome. Functional enrichment analyses of these results suggest that the associations for any given complex trait are not purely random. Thus, we set out to leverage the genetic association results from many traits with a view to identifying the set of modules, or latent factors, that mediate these associations. The identification of such modules may aid in disease classification as well as the elucidation of complex disease mechanisms. We propose a method, Genetic Unmixing by Independent Decomposition (GUIDE), to estimate a set of statistically independent latent factors that best express the patterns of association across many traits. The resulting latent factors not only have desirable mathematical properties, such as sparsity and a higher variance explained (for both traits and variants), but are also able to single out and prioritize key biological features or pathophysiological mechanisms underlying a given trait or disease. Moreover, we show that these latent factors can index biological pathways as well as epidemiological and environmental influences that compose the genetic architecture of complex traits.

Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function

BACKGROUND

Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function.

OBJECTIVES

We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions.

METHODS

We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings.

RESULTS

Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers.

CONCLUSION

Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Human genetic associations of the airway microbiome in chronic obstructive pulmonary disease

Little is known about the relationships between human genetics and the airway microbiome. Deeply sequenced airway metagenomics, by simultaneously characterizing the microbiome and host genetics, provide a unique opportunity to assess the microbiome-host genetic associations. Here we performed a co-profiling of microbiome and host genetics with the identification of over 5 million single nucleotide polymorphisms (SNPs) through deep metagenomic sequencing in sputum of 99 chronic obstructive pulmonary disease (COPD) and 36 healthy individuals. Host genetic variation was the most significant factor associated with the microbiome except for geography and disease status, with its top 5 principal components accounting for 12.11% of the microbiome variability. Within COPD individuals, 113 SNPs mapped to candidate genes reported as genetically associated with COPD exhibited associations with 29 microbial species and 48 functional modules (P < 1 × 10-5), where Streptococcus salivarius exhibits the strongest association to SNP rs6917641 in TBC1D32 (P = 9.54 × 10-8). Integration of concurrent host transcriptomic data identified correlations between the expression of host genes and their genetically-linked microbiome features, including NUDT1, MAD1L1 and Veillonella parvula, TTLL9 and Stenotrophomonas maltophilia, and LTA4H and Haemophilus influenzae. Mendelian randomization analyses revealed a potential causal link between PARK7 expression and microbial type III secretion system, and a genetically-mediated association between COPD and increased relative abundance of airway Streptococcus intermedius. These results suggest a previously underappreciated role of host genetics in shaping the airway microbiome and provide fresh hypotheses for genetic-based host-microbiome interactions in COPD.

Adam19 Deficiency Impacts Pulmonary Function: Human GWAS Follow-up in Mouse

Purpose

Over 550 loci have been associated with human pulmonary function in genome-wide association studies (GWAS); however, the causal role of most remains uncertain. Single nucleotide polymorphisms in a disintegrin and metalloprotease domain 19 (ADAM19) are consistently related to pulmonary function in GWAS. Thus, we used a mouse model to investigate the causal link between Adam19 and pulmonary function.

Methods

We created an Adam19 knockout (KO) mouse model and validated the gene targeting using RNA-Seq and RT-qPCR. Contrary to prior publications, the KO was not neonatal lethal. Thus, we phenotyped the Adam19 KO.

Results

KO mice had lower body weight and shorter tibial length than wild type (WT). Dual-energy X-ray Absorptiometry indicated lower soft weight, fat weight, and bone mineral content in KO mice. In lung function analyses using flexiVent, compared to WT, Adam19 KO had decreased baseline respiratory system elastance, minute work of breathing, tissue damping, tissue elastance, and forced expiratory flow at 50% forced vital capacity but higher FEV0.1 and FVC. Adam19 KO had attenuated tissue damping and tissue elastance in response to methacholine following LPS exposure. Adam19 KO also exhibited attenuated neutrophil extravasation into the airway after LPS administration compared to WT. RNA-Seq analysis of KO and WT lungs identified several differentially expressed genes (Cd300lg, Kpna2, and Pttg1) implicated in lung biology and pathogenesis. Gene set enrichment analysis identified negative enrichment for TNF pathways.

Conclusion

Our murine findings support a causal role of ADAM19, implicated in human GWAS, in regulating pulmonary function.

Multi-omics Integration Identifies Genes Influencing Traits Associated with Cardiovascular Risks: The Long Life Family Study

The Long Life Family Study (LLFS) enrolled 4,953 participants in 539 pedigrees displaying exceptional longevity. To identify genetic mechanisms that affect cardiovascular risks in the LLFS population, we developed a multi-omics integration pipeline and applied it to 11 traits associated with cardiovascular risks. Using our pipeline, we aggregated gene-level statistics from rare-variant analysis, GWAS, and gene expression-trait association by Correlated Meta-Analysis (CMA). Across all traits, CMA identified 64 significant genes after Bonferroni correction (p ≤ 2.8×10-7), 29 of which replicated in the Framingham Heart Study (FHS) cohort. Notably, 20 of the 29 replicated genes do not have a previously known trait-associated variant in the GWAS Catalog within 50 kb. Thirteen modules in Protein-Protein Interaction (PPI) networks are significantly enriched in genes with low meta-analysis p-values for at least one trait, three of which are replicated in the FHS cohort. The functional annotation of genes in these modules showed a significant over-representation of trait-related biological processes including sterol transport, protein-lipid complex remodeling, and immune response regulation. Among major findings, our results suggest a role of triglyceride-associated and mast-cell functional genes FCER1A, MS4A2, GATA2, HDC, and HRH4 in atherosclerosis risks. Our findings also suggest that lower expression of ATG2A, a gene we found to be associated with BMI, may be both a cause and consequence of obesity. Finally, our results suggest that ENPP3 may play an intermediary role in triglyceride-induced inflammation. Our pipeline is freely available and implemented in the Nextflow workflow language, making it easily runnable on any compute platform (https://nf-co.re/omicsgenetraitassociation).

Genetic and observational associations of lung function with gastrointestinal tract diseases: pleiotropic and mendelian randomization analysis

BACKGROUND

The two-way communications along the gut-lung axis influence the immune function in both gut and lung. However, the shared genetic characteristics of lung function with gastrointestinal tract (GIT) diseases remain to be investigated.

METHODS

We first investigated the genetic correlations between three lung function traits and four GIT diseases. Second, we illustrated the genetic overlap by genome-wide pleiotropic analysis (PLACO) and further pinpointed the relevant tissue and cell types by partitioning heritability. Furthermore, we proposed pleiotropic genes as potential drug targets by drug database mining. Finally, we evaluated the causal relationships by epidemiologic observational study and Mendelian randomization (MR) analysis.

RESULTS

We found lung function and GIT diseases were genetically correlated. We identified 258 pleiotropic loci, which were enriched in gut- and lung-specific regions marked by H3K4me1. Among these, 16 pleiotropic genes were targets of drugs, such as tofacitinib and baricitinib targeting TYK2 for the treatment of ulcer colitis and COVID-19, respectively. We identified a missense variant in TYK2, exhibiting a shared causal effect on FEV1/FVC and inflammatory bowel disease (rs12720356, PPLACO=1.38 × 10- 8). These findings suggested TYK2 as a promising drug target. Although the epidemiologic observational study suggested the protective role of lung function in the development of GIT diseases, no causalities were found by MR analysis.

CONCLUSIONS

Our study suggested the shared genetic characteristics between lung function and GIT diseases. The pleiotropic variants could exert their effects by modulating gene expression marked by histone modifications. Finally, we highlighted the potential of pleiotropic analyses in drug repurposing.

Pulmonary function and structure abnormalities in children and young adults with osteogenesis imperfecta point to intrinsic and extrinsic lung abnormalities

PURPOSE

Pulmonary disease is the major cause of morbidity and mortality in osteogenesis imperfecta (OI). We investigated the contribution of intrinsic lung factors to impaired pulmonary function in children and young adults with OI types III, IV, VI.

METHODS

Patients with type III (n=8), IV (n=21), VI (n=5), VII (n=2) or XIV (n=1) OI (mean age 23.6 years) prospectively underwent pulmonary function tests (PFTs) and thoracic CT and radiographs.

RESULTS

PFT results were similar using arm span or ulnar length as height surrogates. PFTs were significantly lower in type III than type IV or VI OI. All patients with type III and half of type IV OI had lung restriction; 90% of patients with OI had reduced gas exchange. Patients with COL1A1 variants had significantly lower forced expiratory flow (FEF)25%-75% compared with those with COL1A2 variants. PFTs correlated negatively with Cobb angle or age. CT scans revealed small airways bronchial thickening (100%, 86%, 100%), atelectasis (88%, 43%, 40%), reticulations (50%, 29%, 20%), ground glass opacities (75%, 5%, 0%), pleural thickening (63%, 48%, 20%) or emphysema (13%, 19%, 20%) in type III, IV or VI OI, respectively.

CONCLUSION

Both lung intrinsic and extrinsic skeletal abnormalities contribute to OI pulmonary dysfunction. Most young adult patients have restrictive disease and abnormal gas exchange; impairment is greater in type III than type IV OI. Decreased FEF25%-75% and thickening of small bronchi walls indicate a critical role for small airways. Lung parenchymal abnormalities (atelectasis, reticulations) and pleural thickening were also detected. Clinical interventions to mitigate these impairments are warranted.

TRIAL REGISTRATION NUMBER

NCT03575221.

Identifying QTLs involved in hybrid performance and heterotic group complementarity: new GWAS models applied to factorial and admixed diallel maize hybrid panels

KEY MESSAGE

An original GWAS model integrating the ancestry of alleles was proposed and allowed the detection of background specific additive and dominance QTLs involved in heterotic group complementarity and hybrid performance. Maize genetic diversity is structured into genetic groups selected and improved relative to each other. This process increases group complementarity and differentiation over time and ensures that the hybrids produced from inter-group crosses exhibit high performances and heterosis. To identify loci involved in hybrid performance and heterotic group complementarity, we introduced an original association study model that disentangles allelic effects from the heterotic group origin of the alleles and compared it with a conventional additive/dominance model. This new model was applied on a factorial between Dent and Flint lines and a diallel between Dent-Flint admixed lines with two different layers of analysis: within each environment and in a multiple-environment context. We identified several strong additive QTLs for all traits, including some well-known additive QTLs for flowering time (in the region of Vgt1/2 on chromosome 8). Yield trait displayed significant non-additive effects in the diallel panel. Most of the detected Yield QTLs exhibited overdominance or, more likely, pseudo-overdominance effects. Apparent overdominance at these QTLs contributed to a part of the genetic group complementarity. The comparison between environments revealed a higher stability of additive QTL effects than non-additive ones. Several QTLs showed variations of effects according to the local heterotic group origin. We also revealed large chromosomic regions that display genetic group origin effects. Altogether, our results illustrate how admixed panels combined with dedicated GWAS modeling allow the identification of new QTLs that could not be revealed by a classical hybrid panel analyzed with traditional modeling.

The Physiological and Pathological Role of Acyl-CoA Oxidation

Fatty acid metabolism, including β-oxidation (βOX), plays an important role in human physiology and pathology. βOX is an essential process in the energy metabolism of most human cells. Moreover, βOX is also the source of acetyl-CoA, the substrate for (a) ketone bodies synthesis, (b) cholesterol synthesis, (c) phase II detoxication, (d) protein acetylation, and (d) the synthesis of many other compounds, including N-acetylglutamate-an important regulator of urea synthesis. This review describes the current knowledge on the importance of the mitochondrial and peroxisomal βOX in various organs, including the liver, heart, kidney, lung, gastrointestinal tract, peripheral white blood cells, and other cells. In addition, the diseases associated with a disturbance of fatty acid oxidation (FAO) in the liver, heart, kidney, lung, alimentary tract, and other organs or cells are presented. Special attention was paid to abnormalities of FAO in cancer cells and the diseases caused by mutations in gene-encoding enzymes involved in FAO. Finally, issues related to α- and ω- fatty acid oxidation are discussed.

CSE reduces OTUD4 triggering lung epithelial cell apoptosis via PAI-1 degradation

Ovarian tumor family deubiquitinase 4 (OTUD4), a member of the OTU deubiquitinating enzyme, is implicated to decrease in cancer to regulate cell apoptosis. However, the role of OTUD4 in cigarette smoke induced epithelial cell apoptosis and its mechanism have not been elucidated. In this study, we showed that OTUD4 protein reduced in CSE treated mice and airway epithelial cells. OTUD4 silence aggravated cell apoptosis and emphysematous change in the lung tissue of cigarette smoke extract (CSE) treated mice. Additionally, restoration of OTUD4 in the lung of mice alleviated CSE induced apoptosis and emphysematous morphology change. The effect of OTUD4 on cell apoptosis was also confirmed in vitro. Through protein profile screening, we identified that OTUD4 may interact with plasminogen activator inhibitor 1(PAI-1). We further confirmed that OTUD4 interacted with PAI-1 for de-ubiquitination and inhibiting CSE induced PAI-1 degradation. Furthermore, the protective role of OTUD4 in airway epithelial cells apoptosis was blocked by PAI-1 deactivation. Taken together, our data suggest that OTUD4 regulates cigarette smoke (CS)-triggered airway epithelial cell apoptosis via modulating PAI-1 degradation. Targeting OUTD4/PAI-1 signaling might potentially provide a therapeutic target against the lung cell apoptosis in cigarette smoke (CS)-induced emphysema.

Systemic Markers of Lung Function and Forced Expiratory Volume in 1 Second Decline across Diverse Cohorts

Rationale: Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by airway obstruction and accelerated lung function decline. Our understanding of systemic protein biomarkers associated with COPD remains incomplete. Objectives: To determine what proteins and pathways are associated with impaired pulmonary function in a diverse population. Methods: We studied 6,722 participants across six cohort studies with both aptamer-based proteomic and spirometry data (4,566 predominantly White participants in a discovery analysis and 2,156 African American cohort participants in a validation). In linear regression models, we examined protein associations with baseline forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC). In linear mixed effects models, we investigated the associations of baseline protein levels with rate of FEV1 decline (ml/yr) in 2,777 participants with up to 7 years of follow-up spirometry. Results: We identified 254 proteins associated with FEV1 in our discovery analyses, with 80 proteins validated in the Jackson Heart Study. Novel validated protein associations include kallistatin serine protease inhibitor, growth differentiation factor 2, and tumor necrosis factor-like weak inducer of apoptosis (discovery β = 0.0561, Q = 4.05 × 10-10; β  = 0.0421, Q = 1.12 × 10-3; and β = 0.0358, Q = 1.67 × 10-3, respectively). In longitudinal analyses within cohorts with follow-up spirometry, we identified 15 proteins associated with FEV1 decline (Q < 0.05), including elafin leukocyte elastase inhibitor and mucin-associated TFF2 (trefoil factor 2; β = -4.3 ml/yr, Q = 0.049; β = -6.1 ml/yr, Q = 0.032, respectively). Pathways and processes highlighted by our study include aberrant extracellular matrix remodeling, enhanced innate immune response, dysregulation of angiogenesis, and coagulation. Conclusions: In this study, we identify and validate novel biomarkers and pathways associated with lung function traits in a racially diverse population. In addition, we identify novel protein markers associated with FEV1 decline. Several protein findings are supported by previously reported genetic signals, highlighting the plausibility of certain biologic pathways. These novel proteins might represent markers for risk stratification, as well as novel molecular targets for treatment of COPD.

Fine-mapping across diverse ancestries drives the discovery of putative causal variants underlying human complex traits and diseases

Genome-wide association studies (GWAS) of human complex traits or diseases often implicate genetic loci that span hundreds or thousands of genetic variants, many of which have similar statistical significance. While statistical fine-mapping in individuals of European ancestries has made important discoveries, cross-population fine-mapping has the potential to improve power and resolution by capitalizing on the genomic diversity across ancestries. Here we present SuSiEx, an accurate and computationally efficient method for cross-population fine-mapping, which builds on the single-population fine-mapping framework, Sum of Single Effects (SuSiE). SuSiEx integrates data from an arbitrary number of ancestries, explicitly models population-specific allele frequencies and LD patterns, accounts for multiple causal variants in a genomic region, and can be applied to GWAS summary statistics. We comprehensively evaluated SuSiEx using simulations, a range of quantitative traits measured in both UK Biobank and Taiwan Biobank, and schizophrenia GWAS across East Asian and European ancestries. In all evaluations, SuSiEx fine-mapped more association signals, produced smaller credible sets and higher posterior inclusion probability (PIP) for putative causal variants, and captured population-specific causal variants.

The complexity in DNA methylation analysis of allergic diseases

PURPOSE OF REVIEW

This review aims to report all the recent studies that are implicated in DNA methylation analysis in the field of allergy and to underline the complexity of the study methodologies and results.

RECENT FINDINGS

Although the growing number of DNA methylation studies have yet to point to a specific mechanism, herein we provide an overview of the majority of pathways considered to be implicated and highlight particular genes, like KNH2 , ATPAF2 and ZNF385A , for their potential as biomarkers.

SUMMARY

The epigenetic profile of respiratory allergic diseases, and particularly DNA methylation, has been investigated in various populations, so as to gain a better understanding of its role in pathogenesis. Through our analysis, multiple links are presented between differential DNA methylation loci and IgE sensitization, lung functionality and severity of the disease. Additionally, associations of this epigenetic change with maternal asthma, age, sex and environmental factors are described, thus uncovering specific gene families that, after further examination could be used as methylation biomarkers in cases of allergic disease.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Latin American Trans-ancestry INitiative for OCD genomics (LATINO): Study Protocol

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder. Worldwide, its prevalence is ~2% and its etiology is mostly unknown. Identifying biological factors contributing to OCD will elucidate underlying mechanisms and might contribute to improved treatment outcomes. Genomic studies of OCD are beginning to reveal long-sought risk loci, but >95% of the cases currently in analysis are of homogenous European ancestry. If not addressed, this Eurocentric bias will result in OCD genomic findings being more accurate for individuals of European ancestry than other ancestries, thereby contributing to health disparities in potential future applications of genomics. In this study protocol paper, we describe the Latin American Trans-ancestry INitiative for OCD genomics (LATINO, www.latinostudy.org). LATINO is a new network of investigators from across Latin America, the United States, and Canada who have begun to collect DNA and clinical data from 5,000 richly-phenotyped OCD cases of Latin American ancestry in a culturally sensitive and ethical manner. In this project, we will utilize trans-ancestry genomic analyses to accelerate the identification of OCD risk loci, fine-map putative causal variants, and improve the performance of polygenic risk scores in diverse populations. We will also capitalize on rich clinical data to examine the genetics of treatment response, biologically plausible OCD subtypes, and symptom dimensions. Additionally, LATINO will help elucidate the diversity of the clinical presentations of OCD across cultures through various trainings developed and offered in collaboration with Latin American investigators. We believe this study will advance the important goal of global mental health discovery and equity.

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 FEV₁/FVC ([Formula: see text] 0.020, SE 0.004, p 4.97 × 10^-08), with suggestive evidence of association with FEV₁ ([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.

Multi-omic association study identifies DNA methylation-mediated genotype and smoking exposure effects on lung function in children living in urban settings

Impaired lung function in early life is associated with the subsequent development of chronic respiratory disease. Most genetic associations with lung function have been identified in adults of European descent and therefore may not represent those most relevant to pediatric populations and populations of different ancestries. In this study, we performed genome-wide association analyses of lung function in a multiethnic cohort of children (n = 1,035) living in low-income urban neighborhoods. We identified one novel locus at the TDRD9 gene in chromosome 14q32.33 associated with percent predicted forced expiratory volume in one second (FEV1) (p = 2.4x10-9; βz = -0.31, 95% CI = -0.41- -0.21). Mendelian randomization and mediation analyses revealed that this genetic effect on FEV1 was partially mediated by DNA methylation levels at this locus in airway epithelial cells, which were also associated with environmental tobacco smoke exposure (p = 0.015). Promoter-enhancer interactions in airway epithelial cells revealed chromatin interaction loops between FEV1-associated variants in TDRD9 and the promoter region of the PPP1R13B gene, a stimulator of p53-mediated apoptosis. Expression of PPP1R13B in airway epithelial cells was significantly associated the FEV1 risk alleles (p = 1.3x10-5; β = 0.12, 95% CI = 0.06-0.17). These combined results highlight a potential novel mechanism for reduced lung function in urban youth resulting from both genetics and smoking exposure.

Early-life farm exposures and eczema among adults in the Agricultural Lung Health Study

Background

Several studies conducted in Europe have suggested a protective association between early-life farming exposures and childhood eczema or atopic dermatitis; few studies have examined associations in adults.

Objectives

To investigate associations between early-life exposures and eczema among 3217 adult farmers and farm spouses (mean age 62.8 years) in a case-control study nested within an US agricultural cohort.

Methods

We used sampling-weighted logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (95%CIs) for associations between early-life exposures and self-reported doctor-diagnosed eczema (273 cases) and polytomous logistic regression to estimate ORs (95%CIs) for a 4-level outcome combining information on eczema and atopy (specific IgE≥0.35). Additionally, we explored genetic and gene-environment associations with eczema.

Results

Although early-life farming exposures were not associated with eczema overall, several early-life exposures were associated with a reduced risk of having both eczema and atopy. Notably, results suggest stronger protective associations among individuals with both eczema and atopy than among those with either atopy alone or eczema alone. For example, ORs (95%CIs) for having a mother who did farm work while pregnant were 1.01 (0.60-1.69) for eczema alone and 0.80 (0.65-0.99) for atopy alone, but 0.54 (0.33-0.80) for having both eczema and atopy. A genetic risk score based on previously identified atopic dermatitis variants was strongly positively associated with eczema, and interaction testing suggested protective effects of several early-life farming exposures only in individuals at lower genetic risk.

Conclusions

In utero and childhood farming exposures are associated with decreased odds of having eczema with atopy in adults.

Global Biobank Meta-analysis Initiative: Powering genetic discovery across human disease

Biobanks facilitate genome-wide association studies (GWASs), which have mapped genomic loci across a range of human diseases and traits. However, most biobanks are primarily composed of individuals of European ancestry. We introduce the Global Biobank Meta-analysis Initiative (GBMI)-a collaborative network of 23 biobanks from 4 continents representing more than 2.2 million consented individuals with genetic data linked to electronic health records. GBMI meta-analyzes summary statistics from GWASs generated using harmonized genotypes and phenotypes from member biobanks for 14 exemplar diseases and endpoints. This strategy validates that GWASs conducted in diverse biobanks can be integrated despite heterogeneity in case definitions, recruitment strategies, and baseline characteristics. This collaborative effort improves GWAS power for diseases, benefits understudied diseases, and improves risk prediction while also enabling the nomination of disease genes and drug candidates by incorporating gene and protein expression data and providing insight into the underlying biology of human diseases and traits.

Best practices for multi-ancestry, meta-analytic transcriptome-wide association studies: Lessons from the Global Biobank Meta-analysis Initiative

The Global Biobank Meta-analysis Initiative (GBMI), through its diversity, provides a valuable opportunity to study population-wide and ancestry-specific genetic associations. However, with multiple ascertainment strategies and multi-ancestry study populations across biobanks, GBMI presents unique challenges in implementing statistical genetics methods. Transcriptome-wide association studies (TWASs) boost detection power for and provide biological context to genetic associations by integrating genetic variant-to-trait associations from genome-wide association studies (GWASs) with predictive models of gene expression. TWASs present unique challenges beyond GWASs, especially in a multi-biobank, meta-analytic setting. Here, we present the GBMI TWAS pipeline, outlining practical considerations for ancestry and tissue specificity, meta-analytic strategies, and open challenges at every step of the framework. We advise conducting ancestry-stratified TWASs using ancestry-specific expression models and meta-analyzing results using inverse-variance weighting, showing the least test statistic inflation. Our work provides a foundation for adding transcriptomic context to biobank-linked GWASs, allowing for ancestry-aware discovery to accelerate genomic medicine.

Whole Genome Association Study of the Plasma Metabolome Identifies Metabolites Linked to Cardiometabolic Disease in Black Individuals

Integrating genetic information with metabolomics has provided new insights into genes affecting human metabolism. However, gene-metabolite integration has been primarily studied in individuals of European Ancestry, limiting the opportunity to leverage genomic diversity for discovery. In addition, these analyses have principally involved known metabolites, with the majority of the profiled peaks left unannotated. Here, we perform a whole genome association study of 2,291 metabolite peaks (known and unknown features) in 2,466 Black individuals from the Jackson Heart Study. We identify 519 locus-metabolite associations for 427 metabolite peaks and validate our findings in two multi-ethnic cohorts. A significant proportion of these associations are in ancestry specific alleles including findings in APOE, TTR and CD36. We leverage tandem mass spectrometry to annotate unknown metabolites, providing new insight into hereditary diseases including transthyretin amyloidosis and sickle cell disease. Our integrative omics approach leverages genomic diversity to provide novel insights into diverse cardiometabolic diseases.

The Role of Polygenic Susceptibility on Air Pollution-Associated Asthma between German and Japanese Elderly Women

Polygenic susceptibility likely influences individual responses to air pollutants and the risk of asthma. We compared the role of polygenic susceptibility on air pollution-associated asthma between German and Japanese women. We investigated women that were enrolled in the German SALIA cohort (n = 771, mean age = 73 years) and the Japanese Shika cohort (n = 847, mean age = 67 years) with known asthma status. Adjusted logistic regression models were used to assess the associations between (1) particulate matter with a median aerodynamic diameter ≤ 2.5μm (PM2.5) and nitrogen dioxide (NO2), (2) polygenic risk scores (PRS), and (3) gene-environment interactions (G × E) with asthma. We found an increased risk of asthma in Japanese women after exposure to low pollutant levels (PM2.5: median = 12.7µg/m3, p-value < 0.001, NO2: median = 8.5µg/m3, p-value < 0.001) and in German women protective polygenic effects (p-value = 0.008). While we found no significant G × E effects, the direction in both groups was that the PRS increased the effect of PM2.5 and decreased the effect of NO2 on asthma. Our study confirms that exposure to low air pollution levels increases the risk of asthma in Japanese women and indicates polygenic effects in German women; however, there was no evidence of G × E effects. Future genome-wide G × E studies should further explore the role of ethnic-specific polygenic susceptibility to asthma.

Phenotypic Causal Inference Using Genome-Wide Association Study Data: Mendelian Randomization and Beyond

statistics for genome-wide association studies (GWAS) are increasingly available for downstream analyses. Meanwhile, the popularity of causal inference methods has grown as we look to gather robust evidence for novel medical and public health interventions. This has led to the development of methods that use GWAS summary statistics for causal inference. Here, we describe these methods in order of their escalating complexity, from genetic associations to extensions of Mendelian randomization that consider thousands of phenotypes simultaneously. We also cover the assumptions and limitations of these approaches before considering the challenges faced by researchers performing causal inference using GWAS data. GWAS summary statistics constitute an important data source for causal inference research that offers a counterpoint to nongenetic methods when triangulating evidence. Continued efforts to address the challenges in using GWAS data for causal inference will allow the full impact of these approaches to be realized.

Multi-ancestry fine-mapping improves precision to identify causal genes in transcriptome-wide association studies

Transcriptome-wide association studies (TWASs) are a powerful approach to identify genes whose expression is associated with complex disease risk. However, non-causal genes can exhibit association signals due to confounding by linkage disequilibrium (LD) patterns and eQTL pleiotropy at genomic risk regions, which necessitates fine-mapping of TWAS signals. Here, we present MA-FOCUS, a multi-ancestry framework for the improved identification of genes underlying traits of interest. We demonstrate that by leveraging differences in ancestry-specific patterns of LD and eQTL signals, MA-FOCUS consistently outperforms single-ancestry fine-mapping approaches with equivalent total sample sizes across multiple metrics. We perform TWASs for 15 blood traits using genome-wide summary statistics (average n_EA = 511 k, n_AA = 13 k) and lymphoblastoid cell line eQTL data from cohorts of primarily European and African continental ancestries. We recapitulate evidence demonstrating shared genetic architectures for eQTL and blood traits between the two ancestry groups and observe that gene-level effects correlate 20% more strongly across ancestries than SNP-level effects. Lastly, we perform fine-mapping using MA-FOCUS and find evidence that genes at TWAS risk regions are more likely to be shared across ancestries than they are to be ancestry specific. Using multiple lines of evidence to validate our findings, we find that gene sets produced by MA-FOCUS are more enriched in hematopoietic categories than alternative approaches (p = 2.36 × 10^-15). Our work demonstrates that including and appropriately accounting for genetic diversity can drive more profound insights into the genetic architecture of complex traits.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV₁, FVC, and FEV₁/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV₁, FVC, or FEV₁/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Heritability Analyses Uncover Shared Genetic Effects of Lung Function and Change over Time

Genetic influence on lung functions has been identified in previous studies; however, the relative longitudinal effects of genetic factors and their interactions with smoking on lung function remain unclear. Here, we identified the longitudinal effects of genetic variants on lung function by determining single nucleotide polymorphism (SNP) heritability and genetic correlations, and by analyzing interactions with smoking. Subject-specific means and annual change rates were calculated for eight spirometric measures obtained from 6622 Korean adults aged 40−69 years every two years for 14 years, and their heritabilities were estimated separately. Statistically significant (p < 0.05) heritability for the subject-specific means of all spirometric measures (8~32%) and change rates of forced expiratory volume in 1 s to forced vital capacity ratio (FEV1/FVC; 16%) and post-bronchodilator FEV1/FVC (17%) were detected. Significant genetic correlations of the change rate with the subject-specific mean were observed for FEV1/FVC (ρg = 0.64) and post-bronchodilator FEV1/FVC (ρg = 0.47). Furthermore, post-bronchodilator FEV1/FVC showed significant heritability of SNP-by-smoking interaction (hGXS2 = 0.4) for the annual change rate. The GWAS also detected genome-wide significant SNPs for FEV1 (rs4793538), FEV1/FVC (rs2704589, rs62201158, and rs9391733), and post-bronchodilator FEV1/FVC (rs2445936). We found statistically significant evidence of heritability role on the change in lung function, and this was shared with the effects on cross-sectional measurements. We also found some evidence of interaction with smoking for the change of lung function.

Bronchial epithelium epithelial-mesenchymal plasticity forms aberrant basaloid-like cells in vitro

Although epithelial-mesenchymal transition (EMT) is a common feature of fibrotic lung disease, its role in fibrogenesis is controversial. Recently, aberrant basaloid cells were identified in fibrotic lung tissue as a novel epithelial cell type displaying a partial EMT phenotype. The developmental origin of these cells remains unknown. To elucidate the role of EMT in the development of aberrant basaloid cells from the bronchial epithelium, we mapped EMT-induced transcriptional changes at the population and single-cell levels. Human bronchial epithelial cells grown as submerged or air-liquid interface (ALI) cultures with or without EMT induction were analyzed by bulk and single-cell RNA-Sequencing. Comparison of submerged and ALI cultures revealed differential expression of 8,247 protein coding (PC) and 1,621 long noncoding RNA (lncRNA) genes and revealed epithelial cell-type-specific lncRNAs. Similarly, EMT induction in ALI cultures resulted in robust transcriptional reprogramming of 6,020 PC and 907 lncRNA genes. Although there was no evidence for fibroblast/myofibroblast conversion following EMT induction, cells displayed a partial EMT gene signature and an aberrant basaloid-like cell phenotype. The substantial transcriptional differences between submerged and ALI cultures highlight that care must be taken when interpreting data from submerged cultures. This work supports that lung epithelial EMT does not generate fibroblasts/myofibroblasts and confirms ALI cultures provide a physiologically relevant system to study aberrant basaloid-like cells and mechanisms of EMT. We provide a catalog of PC and lncRNA genes and an interactive browser (https://bronc-epi-in-vitro.cells.ucsc.edu/) of single-cell RNA-Seq data for further exploration of potential roles in the lung epithelium in health and lung disease.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV₁] and its ratio to forced vital capacity [FEV₁/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV₁ and FEV₁/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10^-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Genetics of chronic obstructive pulmonary disease: understanding the pathobiology and heterogeneity of a complex disorder

Chronic obstructive pulmonary disease (COPD) is a deadly and highly morbid disease. Susceptibility to and heterogeneity of COPD are incompletely explained by environmental factors such as cigarette smoking. Family-based and population-based studies have shown that a substantial proportion of COPD risk is related to genetic variation. Genetic association studies have identified hundreds of genetic variants that affect risk for COPD, decreased lung function, and other COPD-related traits. These genetic variants are associated with other pulmonary and non-pulmonary traits, demonstrate a genetic basis for at least part of COPD heterogeneity, have a substantial effect on COPD risk in aggregate, implicate early-life events in COPD pathogenesis, and often involve genes not previously suspected to have a role in COPD. Additional progress will require larger genetic studies with more ancestral diversity, improved profiling of rare variants, and better statistical methods. Through integration of genetic data with other omics data and comprehensive COPD phenotypes, as well as functional description of causal mechanisms for genetic risk variants, COPD genetics will continue to inform novel approaches to understanding the pathobiology of COPD and developing new strategies for management and treatment.

Genetic Associations and Architecture of Asthma-COPD Overlap

BACKGROUND

Some people have characteristics of both asthma and COPD (asthma-COPD overlap), and evidence suggests they experience worse outcomes than those with either condition alone.

RESEARCH QUESTION

What is the genetic architecture of asthma-COPD overlap, and do the determinants of risk for asthma-COPD overlap differ from those for COPD or asthma?

STUDY DESIGN AND METHODS

We conducted a genome-wide association study in 8,068 asthma-COPD overlap case subjects and 40,360 control subjects without asthma or COPD of European ancestry in UK Biobank (stage 1). We followed up promising signals (P < 5 × 10^-6) that remained associated in analyses comparing (1) asthma-COPD overlap vs asthma-only control subjects, and (2) asthma-COPD overlap vs COPD-only control subjects. These variants were analyzed in 12 independent cohorts (stage 2).

RESULTS

We selected 31 independent variants for further investigation in stage 2, and discovered eight novel signals (P < 5 × 10^-8) for asthma-COPD overlap (meta-analysis of stage 1 and 2 studies). These signals suggest a spectrum of shared genetic influences, some predominantly influencing asthma (FAM105A, GLB1, PHB, TSLP), others predominantly influencing fixed airflow obstruction (IL17RD, C5orf56, HLA-DQB1). One intergenic signal on chromosome 5 had not been previously associated with asthma, COPD, or lung function. Subgroup analyses suggested that associations at these eight signals were not driven by smoking or age at asthma diagnosis, and in phenome-wide scans, eosinophil counts, atopy, and asthma traits were prominent.

INTERPRETATION

We identified eight signals for asthma-COPD overlap, which may represent loci that predispose to type 2 inflammation, and serious long-term consequences of asthma.

Predictive Equations for Adult Pulmonary Function in Zhejiang Province, China

Background

Accurate interpretation of lung function tests requires appropriate spirometry reference values derived from large-scale population-specific epidemiological surveys. The aim of this cross-sectional study was to establish normal spirometric values for the population of healthy, nonsmoking Han Chinese adults residing in Zhejiang province, China.

Methods

We measured lung function parameters such as forced expiratory volume in 1 s, forced vital capacity, peak expiratory flow, maximal midexpiratory flow, and diffusion capacity for carbon monoxide and considered age, height, and weight as independent factors that may modify these parameters. The clinical data were divided into the study arm and validation group. The study arms were used to construct predictive equations using stepwise multiple linear regression, and data from the validation group were used to assess the robustness of the equations.

Results

The 3866 participants were randomized into a study arm (n = 1,949) and a validation arm (n = 1,917). Lung function parameters had a negative association with age and a positive association with height. Data from the two groups were similar. Predictive equations were constructed from the study arm, and the validation group was used to test the feasibility of the reference equations.

Conclusions

The reference values we derived can be used to evaluate lung function in this cohort in both epidemiological studies and clinical practice.

CC16 polymorphisms in asthma, asthma subtypes, and asthma control in adults from the Agricultural Lung Health Study

BACKGROUND

The club cell secretory protein (CC16) has anti-inflammatory and antioxidant effects and is a potential early biomarker of lung damage. The CC16 single nucleotide polymorphism (SNP) rs3741240 risk allele (A) has been inconsistently linked to asthma; other tagging SNPs in the gene have not been explored. The aim was to determine whether CC16 tagging polymorphisms are associated with adult asthma, asthma subtypes or asthma control in the Agricultural Lung Health Study (ALHS).

METHODS

The ALHS is an asthma case-control study nested in the Agricultural Health Study cohort. Asthma cases were individuals with current doctor diagnosed asthma, likely undiagnosed asthma, or asthma-COPD overlap defined by questionnaire. We also examined asthma subtypes and asthma control. Five CC16 tagging SNPs were imputed to 1000 Genomes Integrated phase 1 reference panel. Logistic regression was used to estimate associations between CC16 SNPs and asthma outcomes adjusted for covariates.

RESULTS

The sample included 1120 asthma cases and 1926 controls of European ancestry, with a mean age of 63 years. The frequency of the risk genotype (AA) for rs3741240 was 12.5% (n = 382). CC16 rs3741240 was not associated with adult asthma outcomes. A tagging SNP in the CC16 gene, rs12270961 was associated with uncontrolled asthma (n = 208, OR_adj= 1.4, 95% CI 1.0, 1.9; p = 0.03).

CONCLUSION

This study, the largest study to investigate associations between CC16 tagging SNPs and asthma phenotypes in adults, did not confirm an association of rs3741240 with adult asthma. A tagging SNP in CC16 suggests a potential relationship with asthma control.

Interaction of Cigarette Smoking and Polygenic Risk Score on Reduced Lung Function

IMPORTANCE

The risk of airflow limitation and chronic obstructive pulmonary disease (COPD) is influenced by combinations of cigarette smoking and genetic susceptibility, yet it remains unclear whether gene-by-smoking interactions are associated with quantitative measures of lung function.

OBJECTIVE

To assess the interaction of cigarette smoking and polygenic risk score in association with reduced lung function.

DESIGN, SETTING, AND PARTICIPANTS

This UK Biobank cohort study included UK citizens of European ancestry aged 40 to 69 years with genetic and spirometry data passing quality control metrics. Data was analyzed from July 2020 to March 2021.

EXPOSURES

PRS of combined forced expiratory volume in 1 second (FEV1) and percent of forced vital capacity exhaled in the first second (FEV1/FVC), self-reported pack-years of smoking, ever- vs never-smoking status, and current- vs former- or never-smoking status.

MAIN OUTCOMES AND MEASURES

FEV1/FVC was the primary outcome. Models were used to test for interactions with models, including the main effects of PRS, different smoking variables, and their cross-product terms. The association between pack-years of smoking and FEV1/FVC were compared for those in the highest vs lowest decile of estimated genetic risk for low lung function.

RESULTS

We included 319 730 individuals, of whom 24 915 (8%) had moderate-to-severe COPD cases, and 44.4% were men. Participants had a mean (SD) age 56.5 of (8.02) years. The PRS and pack-years were significantly associated with lower FEV1/FVC (PRS: β, -0.03; 95% CI, -0.031 to -0.03; pack-years: β, -0.0064; 95% CI, -0.0064 to -0.0063) and the interaction term (β, -0.0028; 95% CI, -0.0029 to -0.0026). A stepwise increment in estimated effect sizes for these interaction terms was observed per 10 pack-years of smoking exposure. The interaction of PRS with 11 to 20, 31 to 40, and more than 50 pack-years categories were β (interaction) -0.0038 (95% CI, -0.0046 to -0.0031); -0.013 (95% CI, -0.014 to -0.012); and -0.017 (95% CI, -0.019 to -0.016), respectively. There was evidence of significant interaction between PRS with ever- or never- smoking status (β, interaction; -0.0064; 95% CI, -0.0068 to -0.0060) and current or not-current smoking (β, interaction; -0.0091; 95% CI, -0.0097 to -0.0084). For any given level of pack-years of smoking exposure, FEV1/FVC was significantly lower for individuals in the tenth decile (ie, highest risk) than the first decile (ie, lowest risk) of genetic risk. For every 20 pack-years of smoking, those in the tenth decile compared with the first decile of genetic risk showed nearly a 2-fold reduction in FEV1/FVC.

CONCLUSIONS AND RELEVANCE

COPD is characterized by diminished lung function, and our analyses suggest there is substantial interaction between genome-wide PRS and smoking exposures. While smoking was associated with decreased lung function across all genetic risk categories, the associations were strongest in individuals with higher estimated genetic risk.

Genetic loci for lung function in Japanese adults with adjustment for exhaled nitric oxide levels as airway inflammation indicator

Lung function reflects the ability of the respiratory system and is utilized for the assessment of respiratory diseases. Because type 2 airway inflammation influences lung function, genome wide association studies (GWAS) for lung function would be improved by adjustment with an indicator of the inflammation. Here, we performed a GWAS for lung function with adjustment for exhaled nitric oxide (FeNO) levels in two independent Japanese populations. Our GWAS with genotype imputations revealed that the RNF5/AGER locus including AGER rs2070600 SNP, which introduces a G82S substitution of AGER, was the most significantly associated with FEV1/FVC. Three other rare missense variants of AGER were further identified. We also found genetic loci with three candidate genes (NOS2, SPSB2 and RIPOR2) associated with FeNO levels. Analyses with the BioBank-Japan GWAS resource revealed genetic links of FeNO and asthma-related traits, and existence of common genetic background for allergic diseases and their biomarkers. Our study identified the genetic locus most strongly associated with airway obstruction in the Japanese population and three genetic loci associated with FeNO, an indicator of type 2 airway inflammation in adults.

The genetic case for cardiorespiratory fitness as a clinical vital sign and the routine prescription of physical activity in healthcare

BACKGROUND

Cardiorespiratory fitness (CRF) and physical activity (PA) are well-established predictors of morbidity and all-cause mortality. However, CRF is not routinely measured and PA not routinely prescribed as part of standard healthcare. The American Heart Association (AHA) recently presented a scientific case for the inclusion of CRF as a clinical vital sign based on epidemiological and clinical observation. Here, we leverage genetic data in the UK Biobank (UKB) to strengthen the case for CRF as a vital sign and make a case for the prescription of PA.

METHODS

We derived two CRF measures from the heart rate data collected during a submaximal cycle ramp test: CRF-vo2max, an estimate of the participants' maximum volume of oxygen uptake, per kilogram of body weight, per minute; and CRF-slope, an estimate of the rate of increase of heart rate during exercise. Average PA over a 7-day period was derived from a wrist-worn activity tracker. After quality control, 70,783 participants had data on the two derived CRF measures, and 89,683 had PA data. We performed genome-wide association study (GWAS) analyses by sex, and post-GWAS techniques to understand genetic architecture of the traits and prioritise functional genes for follow-up.

RESULTS

We found strong evidence that genetic variants associated with CRF and PA influenced genetic expression in a relatively small set of genes in the heart, artery, lung, skeletal muscle and adipose tissue. These functionally relevant genes were enriched among genes known to be associated with coronary artery disease (CAD), type 2 diabetes (T2D) and Alzheimer's disease (three of the top 10 causes of death in high-income countries) as well as Parkinson's disease, pulmonary fibrosis, and blood pressure, heart rate, and respiratory phenotypes. Genetic variation associated with lower CRF and PA was also correlated with several disease risk factors (including greater body mass index, body fat and multiple obesity phenotypes); a typical T2D profile (including higher insulin resistance, higher fasting glucose, impaired beta-cell function, hyperglycaemia, hypertriglyceridemia); increased risk for CAD and T2D; and a shorter lifespan.

CONCLUSIONS

Genetics supports three decades of evidence for the inclusion of CRF as a clinical vital sign. Given the genetic, clinical and epidemiological evidence linking CRF and PA to increased morbidity and mortality, regular measurement of CRF as a marker of health and routine prescription of PA could be a prudent strategy to support public health.

A large-scale genome-wide association analysis of lung function in the Chinese population identifies novel loci and highlights shared genetic aetiology with obesity

BACKGROUND

Lung function is a heritable complex phenotype with obesity being one of its important risk factors. However, knowledge of their shared genetic basis is limited. Most genome-wide association studies (GWASs) for lung function have been based on European populations, limiting the generalisability across populations. Large-scale lung function GWASs in other populations are lacking.

METHODS

We included 100 285 subjects from the China Kadoorie Biobank (CKB). To identify novel loci for lung function, single-trait GWAS analyses were performed on forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC in the CKB. We then performed genome-wide cross-trait analysis between lung function and obesity traits (body mass index (BMI), BMI-adjusted waist-to-hip ratio and BMI-adjusted waist circumference) to investigate the shared genetic effects in the CKB. Finally, polygenic risk scores (PRSs) of lung function were developed in the CKB and their interaction with BMI's association on lung function were examined. We also conducted cross-trait analysis in parallel with the CKB using up to 457 756 subjects from the UK Biobank (UKB) for replication and investigation of ancestry-specific effects.

RESULTS

We identified nine genome-wide significant novel loci for FEV1, six for FVC and three for FEV1/FVC in the CKB. FEV1 and FVC showed significant negative genetic correlation with obesity traits in both the CKB and UKB. Genetic loci shared between lung function and obesity traits highlighted important biological pathways, including cell proliferation, embryo, skeletal and tissue development, and regulation of gene expression. Mendelian randomisation analysis suggested significant negative causal effects of BMI on FEV1 and on FVC in both the CKB and UKB. Lung function PRSs significantly modified the effect of change in BMI on change in lung function during an average follow-up of 8 years.

CONCLUSION

This large-scale GWAS of lung function identified novel loci and shared genetic aetiology between lung function and obesity. Change in BMI might affect change in lung function differently according to a subject's polygenic background. These findings may open new avenues for the development of molecular-targeted therapies for obesity and lung function improvement.

Rare and low-frequency exonic variants and gene-by-smoking interactions in pulmonary function

Genome-wide association studies have identified numerous common genetic variants associated with spirometric measures of pulmonary function, including forced expiratory volume in one second (FEV₁), forced vital capacity, and their ratio. However, variants with lower minor allele frequencies are less explored. We conducted a large-scale gene-smoking interaction meta-analysis on exonic rare and low-frequency variants involving 44,429 individuals of European ancestry in the discovery stage and sought replication in the UK BiLEVE study with 45,133 European ancestry samples and UK Biobank study with 59,478 samples. We leveraged data on cigarette smoking, the major environmental risk factor for reduced lung function, by testing gene-by-smoking interaction effects only and simultaneously testing the genetic main effects and interaction effects. The most statistically significant signal that replicated was a previously reported low-frequency signal in GPR126, distinct from common variant associations in this gene. Although only nominal replication was obtained for a top rare variant signal rs142935352 in one of the two studies, interaction and joint tests for current smoking and PDE3B were significantly associated with FEV₁. This study investigates the utility of assessing gene-by-smoking interactions and underscores their effects on potential pulmonary function.

Genome-Wide Association Study of Pelvic Organ Prolapse Using the Michigan Genomics Initiative

OBJECTIVES

The aim of this study was to (1) replicate previously identified genetic variants significantly associated with pelvic organ prolapse and (2) identify new genetic variants associated with pelvic organ prolapse using a genome-wide association study.

METHODS

Using our institution's database linking genetic and clinical data, we identified 1,329 women of European ancestry with an International Classification of Diseases, Ninth Revision (ICD-9)/ICD-10 code for prolapse, 767 of whom also had Current Procedural Terminology (CPT)/ICD-9/ICD-10 procedure codes for prolapse surgery, and 16,383 women of European ancestry older than 40 years without a prolapse diagnosis code as controls. Patients were genotyped using the Illumina HumanCoreExome chip and imputed to the Haplotype Reference Consortium. We tested 20 million single nucleotide polymorphisms (SNPs) for association with pelvic organ prolapse adjusting for relatedness, age, chip version, and 4 principal components. We compared our results with 18 previously identified genome-wide significant SNPs from the UK Biobank, Commun Biol (2020;3:129), and Obstet Gynecol (2011;118:1345-1353).

RESULTS

No variants achieved genome-wide significance (P = 5 × 10-8). However, we replicated 4 SNPs with biologic plausibility at nominal significance (P ≤ 0.05): rs12325192 (P = 0.002), rs9306894 (P = 0.05), rs1920568 (P = 0.034), and rs1247943 (P = 0.041), which were all intergenic and nearest the genes SALL1, GDF7, TBX5, and TBX5, respectively.

CONCLUSIONS

Our replication of 4 biologically plausible previously reported SNPs provides further evidence for a genetic contribution to prolapse, specifically that rs12325192, rs9306894, rs1920568, and rs1247943 may contribute to susceptibility for prolapse. These and previously reported associations that have not yet been replicated should be further explored in larger, more diverse cohorts, perhaps through meta-analysis.

A polygenic score for acute vaso-occlusive pain in pediatric sickle cell disease

Individuals with monogenic disorders can experience variable phenotypes that are influenced by genetic variation. To investigate this in sickle cell disease (SCD), we performed whole-genome sequencing (WGS) of 722 individuals with hemoglobin HbSS or HbSβ0-thalassemia from Baylor College of Medicine and from the St. Jude Children's Research Hospital Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort study. We developed pipelines to identify genetic variants that modulate sickle hemoglobin polymerization in red blood cells and combined these with pain-associated variants to build a polygenic score (PGS) for acute vaso-occlusive pain (VOP). Overall, we interrogated the α-thalassemia deletion -α3.7 and 133 candidate single-nucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the β-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10-14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10-13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.

Interaction between Genetic Risk Scores for reduced pulmonary function and smoking, asthma and endotoxin

Rationale

Genome-wide association studies (GWASs) have identified numerous loci associated with lower pulmonary function. Pulmonary function is strongly related to smoking and has also been associated with asthma and dust endotoxin. At the individual SNP level, genome-wide analyses of pulmonary function have not identified appreciable evidence for gene by environment interactions. Genetic Risk Scores (GRSs) may enhance power to identify gene–environment interactions, but studies are few.

Methods

We analysed 2844 individuals of European ancestry with 1000 Genomes imputed GWAS data from a case–control study of adult asthma nested within a US agricultural cohort. Pulmonary function traits were FEV₁, FVC and FEV₁/FVC. Using data from a recent large meta-analysis of GWAS, we constructed a weighted GRS for each trait by combining the top (p value<5×10⁻⁹) genetic variants, after clumping based on distance (±250 kb) and linkage disequilibrium (r²=0.5). We used linear regression, adjusting for relevant covariates, to estimate associations of each trait with its GRS and to assess interactions.

Results

Each trait was highly significantly associated with its GRS (all three p values<8.9×10⁻⁸). The inverse association of the GRS with FEV₁/FVC was stronger for current smokers (p_interaction=0.017) or former smokers (p_interaction=0.064) when compared with never smokers and among asthmatics compared with non-asthmatics (p_interaction=0.053). No significant interactions were observed between any GRS and house dust endotoxin.

Conclusions

Evaluation of interactions using GRSs supports a greater impact of increased genetic susceptibility on reduced pulmonary function in the presence of smoking or asthma.

Lrp1 Regulation of Pulmonary Function. Follow-Up of Human GWAS in Mice

Human genome-wide association studies (GWASs) have identified more than 270 loci associated with pulmonary function; however, follow-up studies to determine causal genes at these loci are few. SNPs in low-density lipoprotein receptor-related protein 1 (LRP1) are associated with human pulmonary function in GWASs. Using murine models, we investigated the effect of genetic disruption of the Lrp1 gene in smooth muscle cells on pulmonary function in naive animals and after exposure to bacterial LPS or house dust mite extract. Disruption of Lrp1 in smooth muscle cells leads to an increase in tissue resistance, elastance, and tissue elastance at baseline. Furthermore, disruption of Lrp1 in smooth muscle increases airway responsiveness as measured by increased total lung resistance and airway resistance after methacholine. Immune cell counts in BAL fluid were increased in animals with Lrp1 disruption. The difference in airway responsiveness by genotype observed in naive animals was not observed after LPS or house dust mite extract exposure. To further explore the mechanisms contributing to changes in pulmonary function, we identified several ligands dysregulated with Lrp1 disruption in smooth muscle cells. These data suggest that dysregulation of LRP1 in smooth muscle cells affects baseline pulmonary function and airway responsiveness and helps establish LRP1 as the causal gene at this GWAS locus.

Fate or coincidence: do COPD and major depression share genetic risk factors?

Major depressive disorder (MDD) is a common comorbidity in chronic obstructive pulmonary disease (COPD), affecting up to 57% of patients with COPD. Although the comorbidity of COPD and MDD is well established, the causal relationship between these two diseases is unclear. A large-scale electronic health record clinical biobank and genome-wide association study summary statistics for MDD and lung function traits were used to investigate potential shared underlying genetic susceptibility between COPD and MDD. Linkage disequilibrium score regression was used to estimate genetic correlation between phenotypes. Polygenic risk scores (PRS) for MDD and lung function traits were developed and used to perform a phenome-wide association study (PheWAS). Multi-trait-based conditional and joint analysis identified single-nucleotide polymorphisms (SNPs) influencing both lung function and MDD. We found genetic correlations between MDD and all lung function traits were small and not statistically significant. A PRS-MDD was significantly associated with an increased risk of COPD in a PheWAS [odds ratio (OR) = 1.12, 95% confidence interval (CI): 1.09-1.16] when adjusting for age, sex and genetic ancestry, but this relationship became attenuated when controlling for smoking history (OR = 1.08, 95% CI: 1.04-1.13). No significant associations were found between the lung function PRS and MDD. Multi-trait-based conditional and joint analysis identified three SNPs that may contribute to both traits, two of which were previously associated with mood disorders and COPD. Our findings suggest that the observed relationship between COPD and MDD may not be driven by a strong shared genetic architecture.

A genome wide association study for lung function in the Korean population using an exome array

BACKGROUND/AIMS

Lung function is an objective indicator of diagnosis and prognosis of respiratory diseases. Many common genetic variants have been associated with lung function in multiple ethnic populations. We looked for coding variants associated with forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC) in the Korean general population.

METHODS

We carried out exome array analysis and lung function measurements of the FEV1 and FEV1/FVC in 7,524 individuals of the Korean population. We evaluated single variants with minor allele frequency greater than 0.5%. We performed look-ups for candidate coding variants associations in the UK Biobank, SpiroMeta, and CHARGE consortia.

RESULTS

We identified coding variants in the SMIM29 (C6orf1) (p = 1.2 × 10-5) and HMGA1 locus on chromosome 6p21, the GIT2 (p = 6.5 × 10-5) locus on chromosome 12q24, and the ARHGEF40 (p = 9.9 × 10-5) locus on chromosome 14q11 as having a significant association with lung function (FEV1). We also confirmed a previously reported association with lung function and chronic obstructive pulmonary disease in the FAM13A (p = 4.54 × 10-6) locus on chromosome 4q22, in TNXB (p = 1.30 × 10-6) and in AGER (p = 1.09 × 10-8) locus on chromosome 6p21.

CONCLUSION

Our exome array analysis identified that several protein coding variants were associated with lung function in the Korean population. Common coding variants in SMIM29 (C6orf1), HMGA1, GIT2, FAM13A, TNXB, AGER and low-frequency variant in ARHGEF40 potentially affect lung function, which warrant further study.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

Tractor uses local ancestry to enable the inclusion of admixed individuals in GWAS and to boost power

Admixed populations are routinely excluded from genomic studies due to concerns over population structure. Here, we present a statistical framework and software package, Tractor, to facilitate the inclusion of admixed individuals in association studies by leveraging local ancestry. We test Tractor with simulated and empirical two-way admixed African-European cohorts. Tractor generates accurate ancestry-specific effect-size estimates and P values, can boost genome-wide association study (GWAS) power and improves the resolution of association signals. Using a local ancestry-aware regression model, we replicate known hits for blood lipids, discover novel hits missed by standard GWAS and localize signals closer to putative causal variants.

The Emerging Physiological Role of AGMO 10 Years after Its Gene Identification

The gene encoding alkylglycerol monooxygenase (AGMO) was assigned 10 years ago. So far, AGMO is the only known enzyme capable of catalysing the breakdown of alkylglycerols and lyso-alkylglycerophospholipids. With the knowledge of the genetic information, it was possible to relate a potential contribution for mutations in the AGMO locus to human diseases by genome-wide association studies. A possible role for AGMO was implicated by genetic analyses in a variety of human pathologies such as type 2 diabetes, neurodevelopmental disorders, cancer, and immune defence. Deficient catabolism of stored lipids carrying an alkyl bond by an absence of AGMO was shown to impact on the overall lipid composition also outside the ether lipid pool. This review focuses on the current evidence of AGMO in human diseases and summarises experimental evidence for its role in immunity, energy homeostasis, and development in humans and several model organisms. With the progress in lipidomics platform and genetic identification of enzymes involved in ether lipid metabolism such as AGMO, it is now possible to study the consequence of gene ablation on the global lipid pool and further on certain signalling cascades in a variety of model organisms in more detail.

Lung Development Genes and Adult Lung Function

Rationale: Poor lung health in adult life may occur partly through suboptimal growth and development, as suggested by epidemiological evidence pointing to early life risk factors.

Objectives: To systematically investigate the effects of lung development genes on adult lung function.

Methods: Using UK Biobank data, we tested the association of 391 genes known to influence lung development with FVC and FEV₁/FVC. We split the dataset into two random subsets of 207,616 and 138,411 individuals, using the larger subset to select the most promising signals and the smaller subset for replication.

Measurements and Main Results: We identified 55 genes, of which 36 (16 for FVC, 19 for FEV₁/FVC, and one for both) had not been identified in the largest, most recent genome-wide study of lung function. Most of these 36 signals were intronic variants; expression data from blood and lung tissue showed that the majority affect the expression of the genes they lie within. Further testing of 34 of these 36 signals in the CHARGE and SpiroMeta consortia showed that 16 replicated after Bonferroni correction and another 12 replicated at nominal significance level. Of the 55 genes, 53 fell into four biological categories whose function is to regulate organ size and cell integrity (growth factors; transcriptional regulators; cell-to-cell adhesion; extracellular matrix), suggesting that these specific processes are important for adult lung health.

Conclusions: Our study demonstrates the importance of lung development genes in regulating adult lung function and influencing both restrictive and obstructive patterns. Further investigation of these developmental pathways could lead to druggable targets.