Genome-wide joint meta-analysis of SNP and SNP-by-smoking interaction identifies novel loci for pulmonary function

Genotype × environment interactions in gene regulation and complex traits

Genotype × environment interactions (GxE) have long been recognized as a key mechanism underlying human phenotypic variation. Technological developments over the past 15 years have dramatically expanded our appreciation of the role of GxE in both gene regulation and complex traits. The richness and complexity of these datasets also required parallel efforts to develop robust and sensitive statistical and computational approaches. Although our understanding of the genetic architecture of molecular and complex traits has been maturing, a large proportion of complex trait heritability remains unexplained. Furthermore, there are increasing efforts to characterize the effect of environmental exposure on human health. We therefore review GxE in human gene regulation and complex traits, advocating for a comprehensive approach that jointly considers genetic and environmental factors in human health and disease.

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.

TNS1 and NRXN1 Genes Interacting With Early-Life Smoking Exposure in Asthma-Plus-Eczema Susceptibility

PURPOSE

Numerous genes have been associated with allergic diseases (asthma, allergic rhinitis, and eczema), but they explain only part of their heritability. This is partly because most previous studies ignored complex mechanisms such as gene-environment (G-E) interactions and complex phenotypes such as co-morbidity. However, it was recently evidenced that the co-morbidity of asthma-plus-eczema appears as a sub-entity depending on specific genetic factors. Besides, evidence also suggest that gene-by-early life environmental tobacco smoke (ETS) exposure interactions play a role in asthma, but were never investigated for asthma-plus-eczema. To identify genetic variants interacting with ETS exposure that influence asthma-plus-eczema susceptibility.

METHODS

To conduct a genome-wide interaction study (GWIS) of asthma-plus-eczema according to ETS exposure, we applied a 2-stage strategy with a first selection of single nucleotide polymorphisms (SNPs) from genome-wide association meta-analysis to be tested at a second stage by interaction meta-analysis. All meta-analyses were conducted across 4 studies including a total of 5,516 European-ancestry individuals, of whom 1,164 had both asthma and eczema.

RESULTS

Two SNPs showed significant interactions with ETS exposure. They were located in 2 genes, NRXN1 (2p16) and TNS1 (2q35), never reported associated and/or interacting with ETS exposure for asthma, eczema or more generally for allergic diseases. TNS1 is a promising candidate gene because of its link to lung and skin diseases with possible interactive effect with tobacco smoke exposure.

CONCLUSIONS

This first GWIS of asthma-plus-eczema with ETS exposure underlines the importance of studying sub-phenotypes such as co-morbidities as well as G-E interactions to detect new susceptibility genes.

Auditory oddball hypoactivation in schizophrenia

Individuals with schizophrenia (SZ) show aberrant activations, assessed via functional magnetic resonance imaging (fMRI), during auditory oddball tasks. However, associations with cognitive performance and genetic contributions remain unknown. This study compares individuals with SZ to healthy volunteers (HVs) using two cross-sectional data sets from multi-center brain imaging studies. It examines brain activation to auditory oddball targets, and their associations with cognitive domain performance, schizophrenia polygenic risk scores (PRS), and genetic variation (loci). Both sample 1 (137 SZ vs. 147 HV) and sample 2 (91 SZ vs. 98 HV), showed hypoactivation in SZ in the left-frontal pole, and right frontal orbital, frontal pole, paracingulate, intracalcarine, precuneus, supramarginal and hippocampal cortices, and right thalamus. In SZ, precuneus activity was positively related to cognitive performance. Schizophrenia PRS showed a negative correlation with brain activity in the right-supramarginal cortex. GWA analyses revealed significant single-nucleotide polymorphisms associated with right-supramarginal gyrus activity. RPL36 also predicted right-supramarginal gyrus activity. In addition to replicating hypoactivation for oddball targets in SZ, this study identifies novel relationships between regional activity, cognitive performance, and genetic loci that warrant replication, emphasizing the need for continued data sharing and collaborative efforts.

Molecular Factors in PAD2 (PADI2) and PAD4 (PADI4) Are Associated with Interstitial Lung Disease Susceptibility in Rheumatoid Arthritis Patients

Around 50% of rheumatoid arthritis (RA) patients show some extra-articular manifestation, with the lung a usually affected organ; in addition, the presence of anti-citrullinated protein antibodies (ACPA) is a common feature, which is caused by protein citrullination modifications, catalyzed by the peptidyl arginine deiminases (PAD) enzymes. We aimed to identify single nucleotide variants (SNV) in PADI2 and PADI4 genes (PAD2 and PAD4 proteins, respectively) associated with susceptibility to interstitial lung disease (ILD) in RA patients and the PAD2 and PAD4 levels. Material and methods: 867 subjects were included: 118 RA-ILD patients, 133 RA patients, and 616 clinically healthy subjects (CHS). Allelic discrimination was performed in eight SNVs using qPCR, four in PADI2 and four in PADI4. The ELISA technique determined PAD2 and PAD4 levels in serum and bronchoalveolar lavage (BAL) samples, and the population structure was evaluated using 14 informative ancestry markers. Results: The rs1005753-GG (OR = 4.9) in PADI2 and rs11203366-AA (OR = 3.08), rs11203367-GG (OR = 2.4) in PADI4 are associated with genetic susceptibility to RA-ILD as well as the ACTC haplotype (OR = 2.64). In addition, the PAD4 protein is increased in RA-ILD individuals harboring the minor allele homozygous genotype in PADI4 SNVs. Moreover, rs1748033 in PADI4, rs2057094, and rs2076615 in PADI2 are associated with RA susceptibility. In conclusion, in RA patients, single nucleotide variants in PADI4 and PADI2 are associated with ILD susceptibility. The rs1748033 in PADI4 and two different SNVs in PADI2 are associated with RA development but not ILD. PAD4 serum levels are increased in RA-ILD patients.

Gene-Smoking Interaction Analysis for the Identification of Novel Asthma-Associated Genetic Factors

Asthma is a complex heterogeneous disease caused by gene-environment interactions. Although numerous genome-wide association studies have been conducted, these interactions have not been systemically investigated. We sought to identify genetic factors associated with the asthma phenotype in 66,857 subjects from the Health Examination Study, Cardiovascular Disease Association Study, and Korea Association Resource Study cohorts. We investigated asthma-associated gene-environment (smoking status) interactions at the level of single nucleotide polymorphisms, genes, and gene sets. We identified two potentially novel (SETDB1 and ZNF8) and five previously reported (DM4C, DOCK8, MMP20, MYL7, and ADCY9) genes associated with increased asthma risk. Numerous gene ontology processes, including regulation of T cell differentiation in the thymus (GO:0033081), were significantly enriched for asthma risk. Functional annotation analysis confirmed the causal relationship between five genes (two potentially novel and three previously reported genes) and asthma through genome-wide functional prediction scores (combined annotation-dependent depletion, deleterious annotation of genetic variants using neural networks, and RegulomeDB). Our findings elucidate the genetic architecture of asthma and improve the understanding of its biological mechanisms. However, further studies are necessary for developing preventive treatments based on environmental factors and understanding the immune system mechanisms that contribute to the etiology of asthma.

Fibroblast growth factor 10 attenuates chronic obstructive pulmonary disease by protecting against glycocalyx impairment and endothelial apoptosis

Background

The defects and imbalance in lung repair and structural maintenance contribute to the pathogenesis of chronic obstructive pulmonary diseases (COPD), yet the molecular mechanisms that regulate lung repair process are so far incompletely understood. We hypothesized that cigarette smoking causes glycocalyx impairment and endothelial apoptosis in COPD, which could be repaired by the stimulation of fibroblast growth factor 10 (FGF10)/FGF receptor 1 (FGFR1) signaling.

Methods

We used immunostaining (immunohistochemical [IHC] and immunofluorescence [IF]) and enzyme-linked immunosorbent assay (ELISA) to detect the levels of glycocalyx components and endothelial apoptosis in animal models and in patients with COPD. We used the murine emphysema model and in vitro studies to determine the protective and reparative role of FGF10/FGFR1.

Results

Exposure to cigarette smoke caused endothelial glycocalyx impairment and emphysematous changes in murine models and human specimens. Pretreatment of FGF10 attenuated the development of emphysema and the shedding of glycocalyx components induced by CSE in vivo. However, FGF10 did not attenuate the emphysema induced by endothelial-specific killing peptide CGSPGWVRC-GG-_D(KLAKLAK)₂. Mechanistically, FGF10 alleviated smoke-induced endothelial apoptosis and glycocalyx repair through FGFR1/ERK/SOX9/HS6ST1 signaling in vitro. FGF10 was shown to repair pulmonary glycocalyx injury and endothelial apoptosis, and attenuate smoke-induced COPD through FGFR1 signaling.

Conclusions

Our results suggest that FGF10 may serve as a potential therapeutic strategy against COPD via endothelial repair and glycocalyx reconstitution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-022-02193-5.

It is the first time to prove the confirm the endothelial glycocalyx impairment in COPD.

FGF10 attenuates the development of emphysema and the shedding of glycocalyx induced by CSE in vivo.

FGF10 alleviates smoke-induced endothelial apoptosis and glycocalyx repair through FGFR1/ERK/SOX9/HS6ST1 signaling.

The missing heritability in type 1 diabetes

Type 1 diabetes (T1D) is a complex autoimmune disease characterized by an absolute deficiency of insulin. It affects more than 20 million people worldwide and imposes an enormous financial burden on patients. The underlying pathogenic mechanisms of T1D are still obscure, but it is widely accepted that both genetics and the environment play an important role in its onset and development. Previous studies have identified more than 60 susceptible loci associated with T1D, explaining approximately 80%-85% of the heritability. However, most identified variants confer only small increases in risk, which restricts their potential clinical application. In addition, there is still a so-called 'missing heritability' phenomenon. While the gap between known heritability and true heritability in T1D is small compared with that in other complex traits and disorders, further elucidation of T1D genetics has the potential to bring novel insights into its aetiology and provide new therapeutic targets. Many hypotheses have been proposed to explain the missing heritability, including variants remaining to be found (variants with small effect sizes, rare variants and structural variants) and interactions (gene-gene and gene-environment interactions; e.g. epigenetic effects). In the following review, we introduce the possible sources of missing heritability and discuss the existing related knowledge in the context of T1D.

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.

Sex-specific genetic modifiers identified susceptibility of cold stored red blood cells to osmotic hemolysis

Background

Genetic variants have been found to influence red blood cell (RBC) susceptibility to hemolytic stress and affect transfusion outcomes and the severity of blood diseases. Males have a higher susceptibility to hemolysis than females, but little is known about the genetic mechanism contributing to the difference.

Results

To investigate the sex differences in RBC susceptibility to hemolysis, we conducted a sex-stratified genome-wide association study and a genome-wide gene-by-sex interaction scan in a multi-ethnic dataset with 12,231 blood donors who have in vitro osmotic hemolysis measurements during routine blood storage. The estimated SNP-based heritability for osmotic hemolysis was found to be significantly higher in males than in females (0.46 vs. 0.41). We identified SNPs associated with sex-specific susceptibility to osmotic hemolysis in five loci (SPTA1, KCNA6, SLC4A1, SUMO1P1, and PAX8) that impact RBC function and hemolysis.

Conclusion

Our study established a best practice to identify sex-specific genetic modifiers for sexually dimorphic traits in datasets with mixed ancestries, providing evidence of different genetic regulations of RBC susceptibility to hemolysis between sexes. These and other variants may help explain observed sex differences in the severity of hemolytic diseases, such as sickle cell and malaria, as well as the viability of red cell storage and recovery.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08461-4.

Genome-wide gene-air pollution interaction analysis of lung function in 300,000 individuals

BACKGROUND

Impaired lung function is predictive of mortality and is a key component of chronic obstructive pulmonary disease. Lung function has a strong genetic component but is also affected by environmental factors such as increased exposure to air pollution, but the effect of their interactions is not well understood.

OBJECTIVES

To identify interactions between genetic variants and air pollution measures which affect COPD risk and lung function. Additionally, to determine whether previously identified lung function genetic association signals showed evidence of interaction with air pollution, considering both individual effects and combined effects using a genetic risk score (GRS).

METHODS

We conducted a genome-wide gene-air pollution interaction analysis of spirometry measures with three measures of air pollution at home address: particulate matter (PM2.5 & PM10) and nitrogen dioxide (NO2), in approximately 300,000 unrelated European individuals from UK Biobank. We explored air pollution interactions with previously identified lung function signals and determined their combined interaction effect using a GRS.

RESULTS

We identified seven new genome-wide interaction signals (P<5×10-8), and a further ten suggestive interaction signals (P<5×10-7). Additionally, we found statistical evidence of interaction for FEV1/FVC between PM2.5 and previously identified lung function signal, rs10841302, near AEBP2, suggesting increased susceptibility as copies of the G allele increased (but size of the impact was small - interaction beta: -0.363 percentage points, 95% CI: -0.523, -0.203 per 5 µg/m3). There was no observed interaction between air pollutants and the weighted GRS.

DISCUSSION

We carried out the largest genome-wide gene-air pollution interaction study of lung function and identified potential effects of clinically relevant size and significance. We observed up to 440 ml lower lung function for certain genotypes when exposed to mean levels of outdoor air pollution, which is approximately equivalent to nine years of average normal loss of lung function in adults.

Genetic Modifying Factors of Cystic Fibrosis Phenotype: A Challenge for Modern Medicine

Cystic fibrosis (CF) is a monogenic autosomal recessive disease caused by cystic fibrosis transmembrane conductance regulator (CFTR) gene mutations. CF is characterized by a high phenotypic variability present even in patients with the same genotype. This is due to the intervention of modifier genes that interact with both the CFTR gene and environmental factors. The purpose of this review is to highlight the role of non-CFTR genetic factors (modifier genes) that contribute to phenotypic variability in CF. We analyzed literature data starting with candidate gene studies and continuing with extensive studies, such as genome-wide association studies (GWAS) and whole exome sequencing (WES). The results of both types of studies revealed that the number of modifier genes in CF patients is impressive. Their identification offers a new perspective on the pathophysiological mechanisms of the disease, paving the way for the understanding of other genetic disorders. In conclusion, in the future, genetic analysis, such as GWAS and WES, should be performed routinely. A challenge for future research is to integrate their results in the process of developing new classes of drugs, with a goal to improve the prognosis, increase life expectancy, and enhance quality of life among CF patients.

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.

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.

Genetic Influences on Disease Subtypes

Disease classification, or nosology, was historically driven by careful examination of clinical features of patients. As technologies to measure and understand human phenotypes advanced, so too did classifications of disease, and the advent of genetic data has led to a surge in genetic subtyping in the past decades. Although the fundamental process of refining disease definitions and subtypes is shared across diverse fields, each field is driven by its own goals and technological expertise, leading to inconsistent and conflicting definitions of disease subtypes. Here, we review several classical and recent subtypes and subtyping approaches and provide concrete definitions to delineate subtypes. In particular, we focus on subtypes with distinct causal disease biology, which are of primary interest to scientists, and subtypes with pragmatic medical benefits, which are of primary interest to physicians. We propose genetic heterogeneity as a gold standard for establishing biologically distinct subtypes of complex polygenic disease. We focus especially on methods to find and validate genetic subtypes, emphasizing common pitfalls and how to avoid them.

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.

Cystic Fibrosis Lung Disease Modifiers and Their Relevance in the New Era of Precision Medicine

Our understanding of cystic fibrosis (CF) has grown exponentially since the discovery of the cystic fibrosis transmembrane conductance regulator (CFTR) gene in 1989. With evolving genetic and genomic tools, we have come to better understand the role of CFTR genotypes in the pathophysiology of the disease. This, in turn, has paved the way for the development of modulator therapies targeted at mutations in the CFTR, which are arguably one of the greatest advances in the treatment of CF. These modulator therapies, however, do not target all the mutations in CFTR that are seen in patients with CF and, furthermore, a variation in response is seen in patients with the same genotype who are taking modulator therapies. There is growing evidence to support the role of non-CFTR modifiers, both genetic and environmental, in determining the variation seen in CF morbidity and mortality and also in the response to existing therapies. This review focusses on key findings from studies using candidate gene and genome-wide approaches to identify CF modifier genes of lung disease in cystic fibrosis and considers the interaction between modifiers and the response to modulator therapies. As the use of modulator therapies expands and we gain data around outcomes, it will be of great interest to investigate this interaction further. Going forward, it will also be crucial to better understand the relative influence of genomic versus environmental factors. With this understanding, we can truly begin to deliver personalised care by better profiling the likely disease phenotype for each patient and their response to treatment.

Circulating chemerin levels are determined through circulating platelet counts in nondiabetic Taiwanese people: A bidirectional Mendelian randomization study

BACKGROUND AND AIMS

Platelet count (PLT) is a predictor of metabolic and inflammation-related disorders. Platelets can release prochemerin, which acts as a link between coagulation and inflammation and between innate and adaptive immunity. The causal effect between PLT and circulating chemerin level has not been elucidated.

METHODS

Nondiabetic participants with samples in the Taiwan Biobank were recruited for a genome-wide association study (GWAS) based on PLT (17,037 participants) and chemerin levels (3887 participants). A bidirectional Mendelian randomization (MR) study was conducted to determine the association between circulating PLT and chemerin levels.

RESULTS

For a GWAS of PLT, 11 gene loci were found to have genome-wide significance. For a GWAS of chemerin levels, two gene loci, RARRES2 and HLADQA2-HLADQB1, were found to have genome-wide significance. Age, sex, body mass index, leukocyte count, hemoglobin, mean blood pressure, hemoglobin A1C, serum total bilirubin, aspartate aminotransferase, triglyceride, and low-density-lipoprotein cholesterol levels, estimated glomerular filtration rate, and circulating chemerin level were found to be independently associated with PLT through a stepwise regression analysis. A bidirectional MR study revealed weighted genetic risk scores (WGRSs) for PLT were significantly associated with chemerin levels by using a two-stage least-square method in a multivariate analysis (p = 0.0031), and no significant association between chemerin level WGRSs and PLT was noted. Sensitivity analysis further revealed no violation of the exclusion-restriction assumption with PLT-determining genotypes on chemerin levels.

CONCLUSIONS

Through a bidirectional MR analysis, our data revealed that chemerin levels were determined based on circulating PLT. Circulating chemerin levels can be intermediates between PLT and future metabolic and inflammation-related disorders.

Sex-Specific Genetic Associations for Barrett's Esophagus and Esophageal Adenocarcinoma

BACKGROUND & AIMS

Esophageal adenocarcinoma (EA) and its premalignant lesion, Barrett's esophagus (BE), are characterized by a strong and yet unexplained male predominance (with a male-to-female ratio in EA incidence of up to 6:1). Genome-wide association studies (GWAS) have identified more than 20 susceptibility loci for these conditions. However, potential sex differences in genetic associations with BE/EA remain largely unexplored.

METHODS

Given strong genetic overlap, BE and EA cases were combined into a single case group for analysis. These were compared with population-based controls. We performed sex-specific GWAS of BE/EA in 3 separate studies and then used fixed-effects meta-analysis to provide summary estimates for >9 million variants for male and female individuals. A series of downstream analyses were conducted separately in male and female individuals to identify genes associated with BE/EA and the genetic correlations between BE/EA and other traits.

RESULTS

We included 6758 male BE/EA cases, 7489 male controls, 1670 female BE/EA cases, and 6174 female controls. After Bonferroni correction, our meta-analysis of sex-specific GWAS identified 1 variant at chromosome 6q11.1 (rs112894788, KHDRBS2-MTRNR2L9, PBONF = .039) that was statistically significantly associated with BE/EA risk in male individuals only, and 1 variant at chromosome 8p23.1 (rs13259457, PRSS55-RP1L1, PBONF = 0.057) associated, at borderline significance, with BE/EA risk in female individuals only. We also observed strong genetic correlations of BE/EA with gastroesophageal reflux disease in male individuals and obesity in female individuals.

CONCLUSIONS

The identified novel sex-specific variants associated with BE/EA could improve the understanding of the genetic architecture of the disease and the reasons for the male predominance.

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.

Genome-Wide Gene-by-Smoking Interaction Study of Chronic Obstructive Pulmonary Disease

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008-2010) and SpiroMeta Consortium (multiple countries, 1947-2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic β4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.

Whole genome sequence analysis of pulmonary function and COPD in 19,996 multi-ethnic participants

Chronic obstructive pulmonary disease (COPD), diagnosed by reduced lung function, is a leading cause of morbidity and mortality. We performed whole genome sequence (WGS) analysis of lung function and COPD in a multi-ethnic sample of 11,497 participants from population- and family-based studies, and 8499 individuals from COPD-enriched studies in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We identify at genome-wide significance 10 known GWAS loci and 22 distinct, previously unreported loci, including two common variant signals from stratified analysis of African Americans. Four novel common variants within the regions of PIAS1, RGN (two variants) and FTO show evidence of replication in the UK Biobank (European ancestry n ~ 320,000), while colocalization analyses leveraging multi-omic data from GTEx and TOPMed identify potential molecular mechanisms underlying four of the 22 novel loci. Our study demonstrates the value of performing WGS analyses and multi-omic follow-up in cohorts of diverse ancestry.

Effect of 6p21 region on lung function is modified by smoking: a genome-wide interaction study

Smoking is a major risk factor for chronic obstructive pulmonary disease (COPD); however, more than 25% of COPD patients are non-smokers, and gene-by-smoking interactions are expected to affect COPD onset. We aimed to identify the common genetic variants interacting with pack-years of smoking on FEV₁/FVC ratios in individuals with normal lung function. A genome-wide interaction study (GWIS) on FEV₁/FVC was performed for individuals with FEV₁/FVC ratio ≥ 70 in the Korea Associated Resource cohort data, and significant SNPs were validated using data from two other Korean cohorts. The GWIS revealed that rs10947231 and rs8192575 met genome-wide significant levels; For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{H}}_{0} :\beta_{SNP} = \beta_{SNP*pack - years} = 0{\text{ vs H}}_{1} : not {\text{H}}_{0} ,$$\end{document}H0:βSNP=βSNP∗pack-years=0vs H1:notH0, the likelihood ratio (LR) test was conducted, and its P values, P_LR, for rs10947231 and rs8192575 were 2.23 × 10^–12 and 1.18 × 10^–8, respectively. Interaction between rs8192575 and smoking is significantly replicated with two additional data (P_INT = 0.0454, 0.0131). Expression quantitative trait loci, topologically associated domains, and PrediXcan analyses revealed that rs8192575 is significantly associated with AGER expression. SNPs on the 6p21 region are associated with FEV₁/FVC, and the effect of smoking on FEV₁/FVC differs among the associated genotypes.

A Genome-Wide Association Study Identifies Two Novel Susceptible Regions for Squamous Cell Carcinoma of the Head and Neck

To identify genetic variants for risk of squamous cell carcinoma of the head and neck (SCCHN), we conducted a two-phase genome-wide association study consisting of 7,858,089 SNPs in 2,171 cases and 4,493 controls of non-Hispanic white, of which, 434,839 typed and 7,423,250 imputed SNPs were used as the discovery. SNPs with P < 1 × 10-3 were further validated in the OncoArray study of oral and pharynx cancer (5,205 cases and 3,232 controls of European ancestry) from databases of Genotypes and Phenotypes. Meta-analysis of the discovery and replication studies identified one novel locus 6p22.1 (P = 2.96 × 10-9 for the leading rs259919) and two cancer susceptibility loci 6p21.32 (rs3135001, HLA-DQB1) and 6p21.33 (rs1265081, CCHCR1) associated with SCCHN risk. Further stratification by tumor site revealed four known cancer loci (5p15.33, 6p21.32, 6p21.33, and 2p23.1) associated with oral cavity cancer risk and oropharyngeal cancer risk, respectively. In addition, one novel locus 18q22.2 (P = 2.54 × 10-9 for the leading SNP rs142021700) was identified for hypopharynx and larynx cancer risk. For SNPs in those reported or novel loci, we also performed functional annotations by bioinformatics prediction and expression quantitative trait loci analysis. Collectively, our identification of four reported loci (2p23.1, 5p15.33, 6p21.32, and 6p21.33) and two novel loci (6p22.1 and 18q22.2) for SCCHN risk highlight the importance of human leukocyte antigen loci for oropharyngeal cancer risk, suggesting that immunologic mechanisms are implicated in the etiology of this subset of SCCHN. SIGNIFICANCE: Two novel risk loci for SCCHN in non-Hispanic white individuals highlight the importance of immunologic mechanism in the disease etiology.

Dentate gyrus volume deficit in schizophrenia

BACKGROUND

Schizophrenia is associated with robust hippocampal volume deficits but subregion volume deficits, their associations with cognition, and contributing genes remain to be determined.

METHODS

Hippocampal formation (HF) subregion volumes were obtained using FreeSurfer 6.0 from individuals with schizophrenia (n = 176, mean age ± s.d. = 39.0 ± 11.5, 132 males) and healthy volunteers (n = 173, mean age ± s.d. = 37.6 ± 11.3, 123 males) with similar mean age, gender, handedness, and race distributions. Relationships between the HF subregion volume with the largest between group difference, neuropsychological performance, and single-nucleotide polymorphisms were assessed.

RESULTS

This study found a significant group by region interaction on hippocampal subregion volumes. Compared to healthy volunteers, individuals with schizophrenia had significantly smaller dentate gyrus (DG) (Cohen's d = -0.57), Cornu Ammonis (CA) 4, molecular layer of the hippocampus, hippocampal tail, and CA 1 volumes, when statistically controlling for intracranial volume; DG (d = -0.43) and CA 4 volumes remained significantly smaller when statistically controlling for mean hippocampal volume. DG volume showed the largest between group difference and significant positive associations with visual memory and speed of processing in the overall sample. Genome-wide association analysis with DG volume as the quantitative phenotype identified rs56055643 (β = 10.8, p < 5 × 10-8, 95% CI 7.0-14.5) on chromosome 3 in high linkage disequilibrium with MOBP. Gene-based analyses identified associations between SLC25A38 and RPSA and DG volume.

CONCLUSIONS

This study suggests that DG dysfunction is fundamentally involved in schizophrenia pathophysiology, that it may contribute to cognitive abnormalities in schizophrenia, and that underlying biological mechanisms may involve contributions from MOBP, SLC25A38, and RPSA.

Spirometric traits show quantile-dependent heritability, which may contribute to their gene-environment interactions with smoking and pollution

BACKGROUND

"Quantile-dependent expressivity" refers to a genetic effect that is dependent upon whether the phenotype (e.g., spirometric data) is high or low relative to its population distribution. Forced vital capacity (FVC), forced expiratory volume in 1 second (FEV1), and the FEV1/FVC ratio are moderately heritable spirometric traits. The aim of the analyses is to test whether their heritability (h2 ) is constant over all quantiles of their distribution.

METHODS

Quantile regression was applied to the mean age, sex, height and smoking-adjusted spirometric data over multiple visits in 9,993 offspring-parent pairs and 1,930 sibships from the Framingham Heart Study to obtain robust estimates of offspring-parent (βOP), offspring-midparent (βOM), and full-sib regression slopes (βFS). Nonparametric significance levels were obtained from 1,000 bootstrap samples. βOPs were used as simple indicators of quantile-specific heritability (i.e., h 2 = 2βOP/(1+rspouse), where rspouse was the correlation between spouses).

RESULTS

βOP ± standard error (SE) decreased by 0.0009 ± 0.0003 (P = 0.003) with every one-percent increment in the population distribution of FEV1/FVC, i.e., βOP ± SE were: 0.182 ± 0.031, 0.152 ± 0.015; 0.136 ± 0.011; 0.121 ± 0.013; and 0.099 ± 0.013 at the 10th, 25th, 50th, 75th, and 90th percentiles of the FEV1/FVC distribution, respectively. These correspond to h2  ± SEs of 0.350 ± 0.060 at the 10th, 0.292 ± 0.029 at the 25th, 0.262 ± 0.020 at the 50th, 0.234 ± 0.025 at the 75th, and 0.191 ± 0.025 at the 90th percentiles of the FEV1/FVC ratio. Maximum mid-expiratory flow (MMEF) h2  ± SEs increased 0.0025 ± 0.0007 (P = 0.0004) with every one-percent increment in its distribution, i.e.: 0.467 ± 0.046, 0.467 ± 0.033, 0.554 ± 0.038, 0.615 ± 0.042, and 0.675 ± 0.060 at the 10th, 25th, 50th, 75th, and 90th percentiles of its distribution. This was due to forced expiratory flow at 75% of FVC (FEF75%), whose quantile-specific h2 increased an average of 0.0042 ± 0.0008 for every one-percent increment in its distribution. It is speculated that previously reported gene-environment interactions may be partially attributable to quantile-specific h2 , i.e., greater heritability in individuals with lower FEV1/FVC due to smoking or airborne particles exposure vs. nonsmoking, unexposed individuals.

CONCLUSION

Heritabilities of FEV1/FVC, MMEF, and FEF75% from quantile-regression of offspring-parent and sibling spirometric data suggest their quantile-dependent expressivity.

Smoking-by-genotype interaction in type 2 diabetes risk and fasting glucose

Smoking is a potentially causal behavioral risk factor for type 2 diabetes (T2D), but not all smokers develop T2D. It is unknown whether genetic factors partially explain this variation. We performed genome-environment-wide interaction studies to identify loci exhibiting potential interaction with baseline smoking status (ever vs. never) on incident T2D and fasting glucose (FG). Analyses were performed in participants of European (EA) and African ancestry (AA) separately. Discovery analyses were conducted using genotype data from the 50,000-single-nucleotide polymorphism (SNP) ITMAT-Broad-CARe (IBC) array in 5 cohorts from from the Candidate Gene Association Resource Consortium (n = 23,189). Replication was performed in up to 16 studies from the Cohorts for Heart Aging Research in Genomic Epidemiology Consortium (n = 74,584). In meta-analysis of discovery and replication estimates, 5 SNPs met at least one criterion for potential interaction with smoking on incident T2D at p<1x10-7 (adjusted for multiple hypothesis-testing with the IBC array). Two SNPs had significant joint effects in the overall model and significant main effects only in one smoking stratum: rs140637 (FBN1) in AA individuals had a significant main effect only among smokers, and rs1444261 (closest gene C2orf63) in EA individuals had a significant main effect only among nonsmokers. Three additional SNPs were identified as having potential interaction by exhibiting a significant main effects only in smokers: rs1801232 (CUBN) in AA individuals, rs12243326 (TCF7L2) in EA individuals, and rs4132670 (TCF7L2) in EA individuals. No SNP met significance for potential interaction with smoking on baseline FG. The identification of these loci provides evidence for genetic interactions with smoking exposure that may explain some of the heterogeneity in the association between smoking and T2D.

Smoking Modifies the Genetic Risk for Early-Onset Periodontitis

Periodontitis has low-prevalence, highly severe disease manifestations with an early onset and rapid progression. The diagnosis is based on severe destruction of the alveolar bone in adolescents and young adults. Genetic susceptibility variants and smoking are well-established risk factors, but their interactions in modifying disease susceptibility have not been studied. We aimed to identify genetic risk variants of early-onset periodontitis that unmask their effects on tobacco smoke exposure. To this end, we analyzed 79,780,573 common variants in 741 northwest Europeans diagnosed to have >30% bone loss at >2 teeth before 35 y of age, using imputed genotypes of the OmniExpress BeadChip. Never versus ever smokers were compared in a logistic regression analysis via a case-only approach. To explore the effect of tobacco smoke on the expression of the G×S-associated genes, cultures of primary gingival fibroblasts ( n = 9) were exposed to cigarette smoke extract, and transcripts were quantified by reverse transcription polymerase chain reaction. We identified 16 loci for which our analysis suggested an association with G×S increased disease risk ( P < 5 × 10−5). Nine loci had previously been reported to be associated with spirometric measures of pulmonary function by an earlier G×S genome-wide association study. Genome-wide significant cis expression quantitative trait loci were reported for G×S-associated single-nucleotide polymorphisms at ST8SIA1 and SOST, indicating a causal role of these genes in tobacco-related etiopathology. Notably, SOST is a negative regulator of bone growth, and ST8SIA1 has a role in tissue remodeling. Cigarette smoke extract significantly altered the expression of 2 associated genes: SSH1 ( P = 5 × 10−07), which is required for NF-κB activation and innate immune responses to bacterial invasion, and ST8SIA1 ( P = 0.0048). We conclude that the genetic predisposition to early-onset periodontitis is in part triggered by smoking and that tobacco smoke directly affects the expression of genes involved in bone homeostasis, tissue repair, and immune response.

Mitochondrial biology in airway pathogenesis and the role of NRF2

A constant improvement in understanding of mitochondrial biology has provided new insights into mitochondrial dysfunction in human disease pathogenesis. Impaired mitochondrial dynamics caused by various stressors are characterized by structural abnormalities and leakage, compromised turnover, and reactive oxygen species overproduction in mitochondria as well as increased mitochondrial DNA mutation frequency, which leads to modified energy production and mitochondria-derived cell signaling. The mitochondrial dysfunction in airway epithelial, smooth muscle, and endothelial cells has been implicated in diseases including chronic obstructive lung diseases and acute lung injury. Increasing evidence indicates that the NRF2-antioxidant response element (ARE) pathway not only enhances redox defense but also facilitates mitochondrial homeostasis and bioenergetics. Identification of functional or potential AREs further supports the role for Nrf2 in mitochondrial dysfunction-associated airway disorders. While clinical reports indicate mixed efficacy, NRF2 agonists acting on respiratory mitochondrial dynamics are potentially beneficial. In lung cancer, growth advantage provided by sustained NRF2 activation is suggested to be through increased cellular antioxidant defense as well as mitochondria reinforcement and metabolic reprogramming to the preferred pathways to meet the increased energy demands of uncontrolled cell proliferation. Further studies are warranted to better understand NRF2 regulation of mitochondrial functions as therapeutic targets in airway disorders.

Genetic underpinnings of lung function and COPD

Spirometry based measurement of lung function is a global initiative for chronic obstructive lung disease (GOLD) standard to diagnose chronic obstructive pulmonary disease (COPD), one of the leading causes of mortality worldwide. The environmental and behavioural risk factors for COPD includes tobacco smoking, air pollutants and biomass fuel exposure, which can induce one or more abnormal lung function patterns. While smoking remains the primary risk factor, only 15-20% smokers develop COPD, indicating that the genetic factors are also likely to play a role. According to the study of Global Burden of Disease 2015, ∼174 million people across the world have COPD. From a comprehensive literature search conducted using the 'PubMed' and 'GWAS Catalogue' databases, and reviewing the literature available, only a limited number of studies were identified which had attempted to investigate the genetics of COPD and lung volumes, implying a huge research gap. With the advent of genomewide association studies several genetic variants linked to lung function and COPD, like HHIP, HTR4, ADAM19 and GSTCD etc., have been found and validated in different population groups, suggesting their potential role in determining lung volume and risk for COPD. This article aims at reviewing the present knowledge of the genetics of lung function and COPD.

Defining a role for lung function associated gene GSTCD in cell homeostasis

Genome wide association (GWA) studies have reproducibly identified signals on chromosome 4q24 associated with lung function and COPD. GSTCD (Glutathione S-transferase C-terminal domain containing) represents a candidate causal gene in this locus, however little is currently known about the function of this protein. We set out to further our understanding of the role of GSTCD in cell functions and homeostasis using multiple molecular and cellular approaches in airway relevant cells. Recombinant expression of human GSTCD in conjunction with a GST activity assay did not identify any enzymatic activity for two GSTCD isoforms questioning the assignment of this protein to this family of enzymes. Protein structure analyses identified a potential methyltransferase domain contained within GSTCD, with these enzymes linked to cell viability and apoptosis. Targeted knockdown (siRNA) of GSTCD in bronchial epithelial cells identified a role for GSTCD in cell viability as proliferation rates were not altered. To provide greater insight we completed transcriptomic analyses on cells with GSTCD expression knocked down and identified several differentially expressed genes including those implicated in airway biology; fibrosis e.g. TGFBR1 and inflammation e.g. IL6R. Pathway based transcriptomic analyses identified an over-representation of genes related to adipogenesis which may suggest additional functions for GSTCD. These findings identify potential additional functions for GSTCD in the context of airway biology beyond the hypothesised GST activity and warrant further investigation.

The Notch ligand DNER regulates macrophage IFNγ release in chronic obstructive pulmonary disease

BACKGROUND

Chronic Obstructive Pulmonary Disease (COPD) is the third leading cause of death worldwide with no curative therapy. A non-canonical Notch ligand, DNER, has been recently identified in GWAS to associate with COPD severity, but its function and contribution to COPD is unknown.

METHODS

DNER localisation was assessed in lung tissue from healthy and COPD patients, and cigarette smoke (CS) exposed mice. Microarray analysis was performed on WT and DNER deficient M1 and M2 bone marrow-derived macrophages (BMDM), and gene set enrichment undertaken. WT and DNER deficient mice were exposed to CS or filtered air for 3 day and 2 months to assess IFNγ-expressing macrophages and emphysema development. Notch and NFKB active subunits were quantified in WT and DNER deficient LPS-treated and untreated BMDM.

FINDINGS

Immunofluorescence staining revealed DNER localised to macrophages in lung tissue from COPD patients and mice. Human and murine macrophages showed enhanced DNER expression in response to inflammation. Interestingly, pro-inflammatory DNER deficient BMDMs exhibited impaired NICD1/NFKB dependent IFNγ signalling and reduced nuclear NICD1/NFKB translocation. Furthermore, decreased IFNγ production and Notch1 activation in recruited macrophages from CS exposed DNER deficient mice were observed, protecting against emphysema and lung dysfunction.

INTERPRETATION

DNER is a novel protein induced in COPD patients and 6 months CS-exposed mice that regulates IFNγ secretion via non-canonical Notch in pro-inflammatory recruited macrophages. These results provide a new pathway involved in COPD immunity that could contribute to the discovery of innovative therapeutic targets.

FUNDING

This work was supported from the Helmholtz Alliance 'Aging and Metabolic Programming, AMPro'.

Family History, Tobacco Smoking, and Risk of Ischemic Stroke

BACKGROUND AND PURPOSE

Both genetic factors and smoking are associated with ischemic stroke (IS) risk. However, little is known about the potential interaction of these factors. We aimed to assess whether smoking and a positive family history interact to increase the risk of IS.

METHODS

The nationwide prospective study recruited 210,000 men and 300,000 women in 2004 to 2008 at ages 30 to 79 years. During 9.7 years of follow-up, we documented 16,923 and 20,656 incident IS cases in men and women without major chronic diseases at baseline, respectively. Multivariable Cox regression models were used to examine associations between family history and IS. Likelihood ratio tests were used to test the smoking-family history interactions on IS.

RESULTS

About 67.8% (n=135,168) of men ever smoked regularly compared with 2.7% (n=7,775) of women. Among men, a significant interaction between family history and smoking on IS was observed (P for interaction=0.03), with more pronounced association between family history and IS among ever-regular smokers (hazard ratio [HR], 1.21; 95% confidence interval [CI], 1.16 to 1.27) than among never-smokers (HR, 1.11; 95% CI, 1.01 to 1.23). The association between family history and IS among ex-smokers after more than 10 years of cessation (HR, 1.01; 95% CI, 0.85 to 1.20) appeared similar to that among never-smokers. Among women, a similar but not significant interaction between family history and smoking on IS was observed. Ever-regular smokers who had a family history of stroke had the highest risk of IS.

CONCLUSIONS

Among Chinese men, the association of family history with IS was accentuated by smoking, and such accentuation tended to be lowered by cessation.

A powerful and data-adaptive test for rare-variant-based gene-environment interaction analysis

As whole-exome/genome sequencing data become increasingly available in genetic epidemiology research consortia, there is emerging interest in testing the interactions between rare genetic variants and environmental exposures that modify the risk of complex diseases. However, testing rare-variant-based gene-by-environment interactions (GxE) is more challenging than testing the genetic main effects due to the difficulty in correctly estimating the latter under the null hypothesis of no GxE effects and the presence of neutral variants. In response, we have developed a family of powerful and data-adaptive GxE tests, called "aGE" tests, in the framework of the adaptive powered score test, originally proposed for testing the genetic main effects. Using extensive simulations, we show that aGE tests can control the type I error rate in the presence of a large number of neutral variants or a nonlinear environmental main effect, and the power is more resilient to the inclusion of neutral variants than that of existing methods. We demonstrate the performance of the proposed aGE tests using Pancreatic Cancer Case-Control Consortium Exome Chip data. An R package "aGE" is available at http://github.com/ytzhong/projects/.

Omega-3 Fatty Acids and Genome-Wide Interaction Analyses Reveal DPP10-Pulmonary Function Association

RATIONALE

Omega-3 polyunsaturated fatty acids (n-3 PUFAs) have anti-inflammatory properties that could benefit adults with comprised pulmonary health.

OBJECTIVE

To investigate n-3 PUFA associations with spirometric measures of pulmonary function tests (PFTs) and determine underlying genetic susceptibility.

METHODS

Associations of n-3 PUFA biomarkers (α-linolenic acid, eicosapentaenoic acid, docosapentaenoic acid [DPA], and docosahexaenoic acid [DHA]) were evaluated with PFTs (FEV1, FVC, and FEV1/FVC) in meta-analyses across seven cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (N = 16,134 of European or African ancestry). PFT-associated n-3 PUFAs were carried forward to genome-wide interaction analyses in the four largest cohorts (N = 11,962) and replicated in one cohort (N = 1,687). Cohort-specific results were combined using joint 2 degree-of-freedom (2df) meta-analyses of SNP associations and their interactions with n-3 PUFAs.

RESULTS

DPA and DHA were positively associated with FEV1 and FVC (P < 0.025), with evidence for effect modification by smoking and by sex. Genome-wide analyses identified a novel association of rs11693320-an intronic DPP10 SNP-with FVC when incorporating an interaction with DHA, and the finding was replicated (P2df = 9.4 × 10-9 across discovery and replication cohorts). The rs11693320-A allele (frequency, ∼80%) was associated with lower FVC (PSNP = 2.1 × 10-9; βSNP = -161.0 ml), and the association was attenuated by higher DHA levels (PSNP×DHA interaction = 2.1 × 10-7; βSNP×DHA interaction = 36.2 ml).

CONCLUSIONS

We corroborated beneficial effects of n-3 PUFAs on pulmonary function. By modeling genome-wide n-3 PUFA interactions, we identified a novel DPP10 SNP association with FVC that was not detectable in much larger studies ignoring this interaction.

Joint Analysis of Multiple Interaction Parameters in Genetic Association Studies

With growing human genetic and epidemiologic data, there has been increased interest for the study of gene-by-environment (G-E) interaction effects. Still, major questions remain on how to test jointly a large number of interactions between multiple SNPs and multiple exposures. In this study, we first compared the relative performance of four fixed-effect joint analysis approaches using simulated data, considering up to 10 exposures and 300 SNPs: (1) omnibus test, (2) multi-exposure and genetic risk score (GRS) test, (3) multi-SNP and environmental risk score (ERS) test, and (4) GRS-ERS test. Our simulations explored both linear and logistic regression while considering three statistics: the Wald test, the Score test, and the likelihood ratio test (LRT). We further applied the approaches to three large sets of human cohort data (n = 37,664), focusing on type 2 diabetes (T2D), obesity, hypertension, and coronary heart disease with smoking, physical activity, diets, and total energy intake. Overall, GRS-based approaches were the most robust, and had the highest power, especially when the G-E interaction effects were correlated with the marginal genetic and environmental effects. We also observed severe miscalibration of joint statistics in logistic models when the number of events per variable was too low when using either the Wald test or LRT test. Finally, our real data application detected nominally significant interaction effects for three outcomes (T2D, obesity, and hypertension), mainly from the GRS-ERS approach. In conclusion, this study provides guidelines for testing multiple interaction parameters in modern human cohorts including extensive genetic and environmental data.

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

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N = 60,552), African (N = 8429), Asian (N = 9959), and Hispanic/Latino (N = 11,775) ethnicities. We identify over 50 additional loci at genome-wide significance in ancestry-specific or multiethnic meta-analyses. Using recent fine-mapping methods incorporating functional annotation, gene expression, and differences in linkage disequilibrium between ethnicities, we further shed light on potential causal variants and genes at known and newly identified loci. Several of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12. Our study highlights the utility of multiethnic and integrative genomics approaches to extend existing knowledge of the genetics of lung function and clinical relevance of implicated loci.

A Genome-Wide Association Study in Hispanics/Latinos Identifies Novel Signals for Lung Function. The Hispanic Community Health Study/Study of Latinos

RATIONALE

Lung function and chronic obstructive pulmonary disease (COPD) are heritable traits. Genome-wide association studies (GWAS) have identified numerous pulmonary function and COPD loci, primarily in cohorts of European ancestry.

OBJECTIVES

Perform a GWAS of COPD phenotypes in Hispanic/Latino populations to identify loci not previously detected in European populations.

METHODS

GWAS of lung function and COPD in Hispanic/Latino participants from a population-based cohort. We performed replication studies of novel loci in independent studies.

MEASUREMENTS AND MAIN RESULTS

Among 11,822 Hispanic/Latino participants, we identified eight novel signals; three replicated in independent populations of European Ancestry. A novel locus for FEV1 in ZSWIM7 (rs4791658; P = 4.99 × 10-9) replicated. A rare variant (minor allele frequency = 0.002) in HAL (rs145174011) was associated with FEV1/FVC (P = 9.59 × 10-9) in a region previously identified for COPD-related phenotypes; it remained significant in conditional analyses but did not replicate. Admixture mapping identified a novel region, with a variant in AGMO (rs41331850), associated with Amerindian ancestry and FEV1, which replicated. A novel locus for FEV1 identified among ever smokers (rs291231; P = 1.92 × 10-8) approached statistical significance for replication in admixed populations of African ancestry, and a novel SNP for COPD in PDZD2 (rs7709630; P = 1.56 × 10-8) regionally replicated. In addition, loci previously identified for lung function in European samples were associated in Hispanic/Latino participants in the Hispanic Community Health Study/Study of Latinos at the genome-wide significance level.

CONCLUSIONS

We identified novel signals for lung function and COPD in a Hispanic/Latino cohort. Including admixed populations when performing genetic studies may identify variants contributing to genetic etiologies of COPD.

Genome-wide assessment of gene-by-smoking interactions in COPD

Cigarette smoke exposure is a major risk factor in chronic obstructive pulmonary disease (COPD) and its interactions with genetic variants could affect lung function. However, few gene-smoking interactions have been reported. In this report, we evaluated the effects of gene-smoking interactions on lung function using Korea Associated Resource (KARE) data with the spirometric variables—forced expiratory volume in 1 s (FEV₁). We found that variations in FEV₁ were different among smoking status. Thus, we considered a linear mixed model for association analysis under heteroscedasticity according to smoking status. We found a previously identified locus near SOX9 on chromosome 17 to be the most significant based on a joint test of the main and interaction effects of smoking. Smoking interactions were replicated with Gene-Environment of Interaction and phenotype (GENIE), Multi-Ethnic Study of Atherosclerosis-Lung (MESA-Lung), and COPDGene studies. We found that individuals with minor alleles, rs17765644, rs17178251, rs11870732, and rs4793541, tended to have lower FEV₁ values, and lung function decreased much faster with age for smokers. There have been very few reports to replicate a common variant gene-smoking interaction, and our results revealed that statistical models for gene-smoking interaction analyses should be carefully selected.

Established and Emerging Environmental Contributors to Disparities in Asthma and Chronic Obstructive Pulmonary Disease

PURPOSE OF REVIEW

Multiple respiratory diseases, including asthma and chronic obstructive pulmonary disease (COPD), display significant socioeconomic and racial/ethnic disparities. The objective of this review is to evaluate the evidence supporting a link between disproportionate environmental exposures and these health disparities.

RECENT FINDINGS

Studies suggest that various co-occurring factors related to the home environment, neighborhood environment, non-modifiable individual factors, and individual behaviors and attributes can increase or modify the risk of adverse respiratory outcomes among socioeconomically-disadvantaged and racially/ethnically diverse populations. Pollutants in the home environment, including particulate matter, nitrogen dioxide, and pesticides, are elevated among lower socioeconomic status populations and have been implicated in the development or exacerbation of respiratory-related conditions. Neighborhood crime and green space are socioeconomically patterned and linked with asthma outcomes through psychosocial pathways. Non-modifiable individual factors such as genetic predisposition cannot explain environmental health disparities but can increase susceptibility to air pollution and other stressors. Individual behaviors and attributes, including obesity and physical activity, contribute to worse outcomes among those with asthma or COPD.

SUMMARY

The root causes of these multifactorial exposures are complex, but many likely stem from economic forces and racial/ethnic and economic segregation that influence the home environment, neighborhood environment, and access to healthy foods and consumer products. Critical research needs include investigations that characterize exposure to and health implications of numerous stressors simultaneously, both to guard against potential confounding in epidemiological investigations and to consider the cumulative impact of multiple elevated environmental exposures and sociodemographic stressors on health disparities.

Surgical necrotizing enterocolitis in extremely premature neonates is associated with genetic variations in an intergenic region of chromosome 8

BackgroundTwin studies suggest that genetic factors may account for up to 50% increased risk for necrotizing enterocolitis (NEC), but genome-wide association studies for NEC are lacking.MethodsGenotyping was done on Illumina BeadChip, followed by analysis using PLINK with logistic regression under an additive model.ResultsAmong 751 extremely-low-birth-weight (<1,000 g, >401 g) neonates, 30 had surgical NEC. Two hundred and sixty-one single-nucleotide polymorphisms (SNPs) showed association with NEC at P<0.05, of which 35 were significant at P<10-7. Minor allele(s) in a cluster of SNPs spanning a 43-kb region of chromosome 8 (8q23.3) conferred an odds ratio of 4.72 (95% confidence interval (CI): 2.51-8.88) for elevated risk of NEC. Two smaller clusters on chromosome 14 and chromosome 11 exhibited P values of 10-7-10-8. The chromosome 8 cluster is in an intergenic region between CUB and Sushi multiple domains 3 (-1.43 Mb) and trichorhinophalangeal syndrome I (+542 kb). RNA sequencing in this region identified a potential novel open-reading frame corresponding to a long interspersed element-1 retrotransposable element.ConclusionGenetic variation in an intergenic region of chromosome 8 is associated with increased risk for NEC with a mechanism that is yet to be identified.

Association between Genetic Variations of MERTK and Chronic Obstructive Pulmonary Disease in Koreans

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a debilitating lung disease. To date, a large number of clinical studies have been conducted to investigate the association between genetic variations and COPD. However, little is known regarding the genetic susceptibility of Koreans to this disease. MER receptor tyrosine kinase (MERTK) plays important roles in the inhibition of inflammation and in the clearance of apoptotic cells. Here, we investigated the association between genetic variations in MERTK and the development of COPD in Koreans.

METHODS

We conducted genetic analysis of MERTK using genomic DNA samples from 87 patients with COPD and 88 healthy controls and compared the frequency of each variation or haplotype between the patient and control groups. Subsequently, the effect of each variation was evaluated using in vitro assays.

RESULTS

Ten variations were identified in this study, four of them for the first time. In addition, we found that the frequency of each variation or haplotype was comparable between the patient and control groups. However, we observed that the frequency for the wild-type haplotype was higher in the control group, compared to that in the group of patients with COPD, in the subgroup analysis of current smokers, although the difference was not statistically significant (P = 0.080). In in vitro assays, we observed that none of the variations affected the activity of the promoter or the expression of MERTK.

CONCLUSION

Our findings indicate that the susceptibility to COPD is not related to the genetic variations or haplotypes of MERTK in Koreans.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV ₁), forced vital capacity (FVC) and the ratio of FEV ₁ to FVC (FEV ₁/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 ^-7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.