Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.

DNA methylation is globally disrupted and associated with expression changes in chronic obstructive pulmonary disease small airways

DNA methylation is an epigenetic modification that is highly disrupted in response to cigarette smoke and involved in a wide spectrum of malignant and nonmalignant diseases, but surprisingly not previously assessed in small airways of patients with chronic obstructive pulmonary disease (COPD). Small airways are the primary sites of airflow obstruction in COPD. We sought to determine whether DNA methylation patterns are disrupted in small airway epithelia of patients with COPD, and evaluate whether changes in gene expression are associated with these disruptions. Genome-wide methylation and gene expression analysis were performed on small airway epithelial DNA and RNA obtained from the same patient during bronchoscopy, using Illumina's Infinium HM27 and Affymetrix's Genechip Human Gene 1.0 ST arrays. To control for known effects of cigarette smoking on DNA methylation, methylation and gene expression profiles were compared between former smokers with and without COPD matched for age, pack-years, and years of smoking cessation. Our results indicate that aberrant DNA methylation is (1) a genome-wide phenomenon in small airways of patients with COPD, and (2) associated with altered expression of genes and pathways important to COPD, such as the NF-E2-related factor 2 oxidative response pathway. DNA methylation is likely an important mechanism contributing to modulation of genes important to COPD pathology. Because these methylation events may underlie disease-specific gene expression changes, their characterization is a critical first step toward the development of epigenetic markers and an opportunity for developing novel epigenetic therapeutic interventions for COPD.

Longer telomere length in COPD patients with α1-antitrypsin deficiency independent of lung function

Oxidative stress is involved in the pathogenesis of airway obstruction in α₁-antitrypsin deficient patients. This may result in a shortening of telomere length, resulting in cellular senescence. To test whether telomere length differs in α₁-antitrypsin deficient patients compared with controls, we measured telomere length in DNA from peripheral blood cells of 217 α₁-antitrypsin deficient patients and 217 control COPD patients. We also tested for differences in telomere length between DNA from blood and DNA from lung tissue in a subset of 51 controls. We found that telomere length in the blood was significantly longer in α₁-antitrypsin deficient COPD patients compared with control COPD patients (p = 1×10⁻²⁹). Telomere length was not related to lung function in α₁-antitrypsin deficient patients (p = 0.3122) or in COPD controls (p = 0.1430). Although mean telomere length was significantly shorter in the blood when compared with the lungs (p = 0.0078), telomere length was correlated between the two tissue types (p = 0.0122). Our results indicate that telomere length is better preserved in α₁-antitrypsin deficient COPD patients than in non-deficient patients. In addition, measurement of telomere length in the blood may be a suitable surrogate for measurement in the lung.

Susceptibility to chronic mucus hypersecretion, a genome wide association study

Background

Chronic mucus hypersecretion (CMH) is associated with an increased frequency of respiratory infections, excess lung function decline, and increased hospitalisation and mortality rates in the general population. It is associated with smoking, but it is unknown why only a minority of smokers develops CMH. A plausible explanation for this phenomenon is a predisposing genetic constitution. Therefore, we performed a genome wide association (GWA) study of CMH in Caucasian populations.

Methods

GWA analysis was performed in the NELSON-study using the Illumina 610 array, followed by replication and meta-analysis in 11 additional cohorts. In total 2,704 subjects with, and 7,624 subjects without CMH were included, all current or former heavy smokers (≥20 pack-years). Additional studies were performed to test the functional relevance of the most significant single nucleotide polymorphism (SNP).

Results

A strong association with CMH, consistent across all cohorts, was observed with rs6577641 (p = 4.25×10⁻⁶, OR = 1.17), located in intron 9 of the special AT-rich sequence-binding protein 1 locus (SATB1) on chromosome 3. The risk allele (G) was associated with higher mRNA expression of SATB1 (4.3×10⁻⁹) in lung tissue. Presence of CMH was associated with increased SATB1 mRNA expression in bronchial biopsies from COPD patients. SATB1 expression was induced during differentiation of primary human bronchial epithelial cells in culture.

Conclusions

Our findings, that SNP rs6577641 is associated with CMH in multiple cohorts and is a cis-eQTL for SATB1, together with our additional observation that SATB1 expression increases during epithelial differentiation provide suggestive evidence that SATB1 is a gene that affects CMH.

APOM and high-density lipoprotein cholesterol are associated with lung function and per cent emphysema

Chronic obstructive pulmonary disease (COPD) is linked to cardiovascular disease; however, there are few studies on the associations of cardiovascular genes with COPD. We assessed the association of lung function with 2100 genes selected for cardiovascular diseases among 20 077 European-Americans and 6900 African-Americans. We performed replication of significant loci in the other racial group and an independent consortium of Europeans, tested the associations of significant loci with per cent emphysema and examined gene expression in an independent sample. We then tested the association of a related lipid biomarker with forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) ratio and per cent emphysema. We identified one new polymorphism for FEV1/FVC (rs805301) in European-Americans (p=1.3×10(-6)) and a second (rs707974) in the combined European-American and African-American analysis (p=1.38×10(-7)). Both single-nucleotide polymorphisms (SNPs) flank the gene for apolipoprotein M (APOM), a component of high-density lipoprotein (HDL) cholesterol. Both were replicated in an independent cohort. SNPs in a second gene related to apolipoprotein M and HDL, PCSK9, were associated with FEV1/FVC ratio among African-Americans. rs707974 was associated with per cent emphysema among European-Americans and African-Americans and APOM expression was related to FEV1/FVC ratio and per cent emphysema. Higher HDL levels were associated with lower FEV1/FVC ratio and greater per cent emphysema. These findings suggest a novel role for the apolipoprotein M/HDL pathway in the pathogenesis of COPD and emphysema.

Admixture fine-mapping in African Americans implicates XAF1 as a possible sarcoidosis risk gene

Sarcoidosis is a complex, multi-organ granulomatous disease with a likely genetic component. West African ancestry confers a higher risk for sarcoidosis than European ancestry. Admixture mapping provides the most direct method to locate genes that underlie such ethnic variation in disease risk. We sought to identify genetic risk variants within four previously-identified ancestry-associated regions—6p24.3–p12.1, 17p13.3–13.1, 2p13.3–q12.1, and 6q23.3–q25.2—in a sample of 2,727 African Americans. We used logistic regression fit by generalized estimating equations and the MIX score statistic to determine which variants within ancestry-associated regions were associated with risk and responsible for the admixture signal. Fine mapping was performed by imputation, based on a previous genome-wide association study; significant variants were validated by direct genotyping. Within the 6p24.3–p12.1 locus, the most significant ancestry-adjusted SNP was rs74318745 (p = 9.4*10⁻¹¹), an intronic SNP within the HLA-DRA gene that did not solely explain the admixture signal, indicating the presence of more than a single risk variant within this well-established sarcoidosis risk region. The locus on chromosome 17p13.3–13.1 revealed a novel sarcoidosis risk SNP, rs6502976 (p = 9.5*10⁻⁶), within intron 5 of the gene X-linked Inhibitor of Apoptosis Associated Factor 1 (XAF1) that accounted for the majority of the admixture linkage signal. Immunohistochemical expression studies demonstrated lack of expression of XAF1 and a corresponding high level of expression of its downstream target, X-linked Inhibitor of Apoptosis (XIAP) in sarcoidosis granulomas. In conclusion, ancestry and association fine mapping revealed a novel sarcoidosis susceptibility gene, XAF1, which has not been identified by previous genome-wide association studies. Based on the known biology of the XIAP/XAF1 apoptosis pathway and the differential expression patterns of XAF1 and XIAP in sarcoidosis granulomas, we suggest that this pathway may play a role in the maintenance of sarcoidosis granulomas.

A draft de novo genome assembly for the northern bobwhite (Colinus virginianus) reveals evidence for a rapid decline in effective population size beginning in the Late Pleistocene

Wild populations of northern bobwhites (Colinus virginianus; hereafter bobwhite) have declined across nearly all of their U.S. range, and despite their importance as an experimental wildlife model for ecotoxicology studies, no bobwhite draft genome assembly currently exists. Herein, we present a bobwhite draft de novo genome assembly with annotation, comparative analyses including genome-wide analyses of divergence with the chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genomes, and coalescent modeling to reconstruct the demographic history of the bobwhite for comparison to other birds currently in decline (i.e., scarlet macaw; Ara macao). More than 90% of the assembled bobwhite genome was captured within <40,000 final scaffolds (N50 = 45.4 Kb) despite evidence for approximately 3.22 heterozygous polymorphisms per Kb, and three annotation analyses produced evidence for >14,000 unique genes and proteins. Bobwhite analyses of divergence with the chicken and zebra finch genomes revealed many extremely conserved gene sequences, and evidence for lineage-specific divergence of noncoding regions. Coalescent models for reconstructing the demographic history of the bobwhite and the scarlet macaw provided evidence for population bottlenecks which were temporally coincident with human colonization of the New World, the late Pleistocene collapse of the megafauna, and the last glacial maximum. Demographic trends predicted for the bobwhite and the scarlet macaw also were concordant with how opposing natural selection strategies (i.e., skewness in the r-/K-selection continuum) would be expected to shape genome diversity and the effective population sizes in these species, which is directly relevant to future conservation efforts.

Risk loci for chronic obstructive pulmonary disease: a genome-wide association study and meta-analysis

BACKGROUND

The genetic risk factors for susceptibility to chronic obstructive pulmonary disease (COPD) are still largely unknown. Additional genetic variants are likely to be identified by genome-wide association studies in larger cohorts or specific subgroups. We sought to identify risk loci for moderate to severe and severe COPD with data from several cohort studies.

METHODS

We combined genome-wide association analysis data from participants in the COPDGene study (non-Hispanic white and African-American ethnic origin) and the ECLIPSE, NETT/NAS, and Norway GenKOLS studies (self-described white ethnic origin). We did analyses comparing control individuals with individuals with moderate to severe COPD and with a subset of individuals with severe COPD. Single nucleotide polymorphisms yielding a p value of less than 5 × 10(-7) in the meta-analysis at loci not previously described were genotyped in individuals from the family-based ICGN study. We combined results in a joint meta-analysis (threshold for significance p<5 × 10(-8)).

FINDINGS

Analysis of 6633 individuals with moderate to severe COPD and 5704 control individuals confirmed association at three known loci: CHRNA3 (p=6·38 × 10(-14)), FAM13A (p=1·12 × 10(-14)), and HHIP (p=1·57 × 10(-12)). We also showed significant evidence of association at a novel locus near RIN3 (p=5·25 × 10(-9)). In the overall meta-analysis (ie, including data from 2859 ICGN participants), the association with RIN3 remained significant (p=5·4 × 10(-9)). 3497 individuals were included in our analysis of severe COPD. The effect estimates for the loci near HHIP and CHRNA3 were significantly stronger in severe disease than in moderate to severe disease (p<0·01). We also identified associations at two additional loci: MMP12 (overall joint meta-analysis p=2·6 × 10(-9)) and TGFB2 (overall joint meta-analysis p=8·3 × 10(-9)).

INTERPRETATION

We have confirmed associations with COPD at three known loci and identified three new genome-wide significant associations. Genetic variants other than in α-1 antitrypsin increase the risk of COPD.

FUNDING

US National Heart, Lung, and Blood Institute; the Alpha-1 Foundation; the COPD Foundation through contributions from AstraZeneca, Boehringer Ingelheim, Novartis, and Sepracor; GlaxoSmithKline; Centers for Medicare and Medicaid Services; Agency for Healthcare Research and Quality; and US Department of Veterans Affairs.

Follow-up on genome-wide main effects: do polymorphisms modify the air pollution effect on lung function decline in adults?

Improved air quality has been found associated with attenuated age-related decline in lung function. But whether genetic polymorphisms strongly associated with lung function play a modifying role in this attenuation process has so far not been investigated. We selected ten single nucleotide polymorphisms derived from the largest genome-wide association studies on lung function and examined whether they modified the association between the change in exposure to particulate matter ≤10μm (ΔPM10) and lung function decline. 4310 participants from the SAPALDIA cohort provided valid spirometry measurements, a detailed pulmonary health questionnaire both at baseline and 11years later as well as blood samples for genetic testing. Spatially and temporally resolved air pollution exposures were assigned on an individual level based on participants' residences. Statistically significant interactions of moderate strength with ΔPM10 were detected for rs2284746. Individuals with the CC genotype had a 21ml slower annual decline of the mid expiratory flow per 10μg/m(3) PM10 reduction over an 10-year period, while the benefits of CG and GG carriers were smaller (14 and 7ml per year, respectively; Pinteraction=0.04). The attenuated annual decline in the percentage of the forced expiratory volume in one second relative to the forced vital capacity (FEV1/FVC) was also increased with the presence of each C-allele (Pinteraction=0.009). We observed further suggestive interactions of similar magnitude in never-smokers, but none of the results would remain statistically significant after correction for multiple testing. We could not find strong evidence that lung function benefits from improved air quality are modified by polymorphisms associated with lung function level in large meta-analyzed genome-wide association studies.

Genome-wide study of percent emphysema on computed tomography in the general population. The Multi-Ethnic Study of Atherosclerosis Lung/SNP Health Association Resource Study

RATIONALE

Pulmonary emphysema overlaps partially with spirometrically defined chronic obstructive pulmonary disease and is heritable, with moderately high familial clustering.

OBJECTIVES

To complete a genome-wide association study (GWAS) for the percentage of emphysema-like lung on computed tomography in the Multi-Ethnic Study of Atherosclerosis (MESA) Lung/SNP Health Association Resource (SHARe) Study, a large, population-based cohort in the United States.

METHODS

We determined percent emphysema and upper-lower lobe ratio in emphysema defined by lung regions less than -950 HU on cardiac scans. Genetic analyses were reported combined across four race/ethnic groups: non-Hispanic white (n = 2,587), African American (n = 2,510), Hispanic (n = 2,113), and Chinese (n = 704) and stratified by race and ethnicity.

MEASUREMENTS AND MAIN RESULTS

Among 7,914 participants, we identified regions at genome-wide significance for percent emphysema in or near SNRPF (rs7957346; P = 2.2 × 10(-8)) and PPT2 (rs10947233; P = 3.2 × 10(-8)), both of which replicated in an additional 6,023 individuals of European ancestry. Both single-nucleotide polymorphisms were previously implicated as genes influencing lung function, and analyses including lung function revealed independent associations for percent emphysema. Among Hispanics, we identified a genetic locus for upper-lower lobe ratio near the α-mannosidase-related gene MAN2B1 (rs10411619; P = 1.1 × 10(-9); minor allele frequency [MAF], 4.4%). Among Chinese, we identified single-nucleotide polymorphisms associated with upper-lower lobe ratio near DHX15 (rs7698250; P = 1.8 × 10(-10); MAF, 2.7%) and MGAT5B (rs7221059; P = 2.7 × 10(-8); MAF, 2.6%), which acts on α-linked mannose. Among African Americans, a locus near a third α-mannosidase-related gene, MAN1C1 (rs12130495; P = 9.9 × 10(-6); MAF, 13.3%) was associated with percent emphysema.

CONCLUSIONS

Our results suggest that some genes previously identified as influencing lung function are independently associated with emphysema rather than lung function, and that genes related to α-mannosidase may influence risk of emphysema.

Gene set enrichment analysis of the NF-κB/Snail/YY1/RKIP circuitry in multiple myeloma

The presence of a dysregulated NF-κB/Snail/YY1/RKIP loop was recently established in metastatic prostate cancer cells and non-Hodgkin's lymphoma; however, its involvement in multiple myeloma (MM) has yet to be investigated. Aim of the study was to investigate the role of the NF-κB/Snail/YY1/RKIP circuitry in MM and how each gene is correlated with the remaining genes of the loop. Using gene set enrichment analysis and gene neighbours analysis in data received from four datasets included in the Multiple Myeloma Genomics Portal of the Multiple Myeloma Research Consortium, we identified various enriched gene sets associated with each member of the NF-κB/Snail/YY1/RKIP circuitry. In each dataset, the 20 most co-expressed genes with the circuitry genes were isolated subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment. Among many, we highlighted on FNDC3B, TPD52, BBX, MBNL1 and MFAP2. Many co-expressed genes participated in the regulation of metabolic processes and nucleic acid binding, or were transcription factor binding genes and genes with metallopeptidase activity. The transcription factors FOXO4, GATA binding factor, Sp1 and AP4 most likely affect the expression of the NF-κB/Snail/YY1/RKIP circuitry genes. Computational analysis of various GEO datasets revealed elevated YY1 and RKIP levels in MM vs. the normal plasma cells, as well as elevated RKIP levels in MM vs. normal B lymphocytes. The present study highlights the relationships of the NF-κB/Snail/YY1/RKIP circuitry genes with specific cancer-related gene sets in multiple myeloma.

Biology of the KCNQ1 Potassium Channel

Ion channels are essential for basic cellular function and for processes including sensory perception and intercellular communication in multicellular organisms. Voltage-gated potassium (Kv) channels facilitate dynamic cellular repolarization during an action potential, opening in response to membrane depolarization to facilitate K+ efflux. In both excitable and nonexcitable cells other, constitutively active, K+ channels provide a relatively constant repolarizing force to control membrane potential, ion homeostasis, and secretory processes. Of the forty known human Kv channel pore-forming α subunits that coassemble in various combinations to form the fundamental tetrameric channel pore and voltage sensor module, KCNQ1 is unique. KCNQ1 stands alone in having the capacity to form either channels that are voltage-dependent and require membrane depolarization for activation, or constitutively active channels. In mammals, KCNQ1 regulates processes including gastric acid secretion, thyroid hormone biosynthesis, salt and glucose homeostasis, and cell volume and in some species is required for rhythmic beating of the heart. In this review, the author discusses the unique functional properties, regulation, cell biology, diverse physiological roles, and involvement in human disease states of this chameleonic K+ channel.

Arrhythmogenic KCNE gene variants: current knowledge and future challenges

There are twenty-five known inherited cardiac arrhythmia susceptibility genes, all of which encode either ion channel pore-forming subunits or proteins that regulate aspects of ion channel biology such as function, trafficking, and localization. The human KCNE gene family comprises five potassium channel regulatory subunits, sequence variants in each of which are associated with cardiac arrhythmias. KCNE gene products exhibit promiscuous partnering and in some cases ubiquitous expression, hampering efforts to unequivocally correlate each gene to specific native potassium currents. Likewise, deducing the molecular etiology of cardiac arrhythmias in individuals harboring rare KCNE gene variants, or more common KCNE polymorphisms, can be challenging. In this review we provide an update on putative arrhythmia-causing KCNE gene variants, and discuss current thinking and future challenges in the study of molecular mechanisms of KCNE-associated cardiac rhythm disturbances.

Targeting EGFR signalling in chronic lung disease: therapeutic challenges and opportunities

Chronic respiratory diseases, including pulmonary fibrosis, chronic obstructive pulmonary disease (COPD) and lung cancer, are the second leading cause of death among Europeans. Despite this, there have been only a few therapeutic advances in these conditions over the past 20 years. In this review we provide evidence that targeting the epidermal growth factor receptor (EGFR) signalling pathway may represent a novel therapeutic panacea for treating chronic lung disease. Using evidence from human patient samples, transgenic animal models, and cell and molecular biology studies we highlight the roles of this signalling pathway in lung development, homeostasis, repair, and disease ontogeny. We identify mechanisms underlying lung EGFR pathway regulation and suggest how targeting these mechanisms using new and existing therapies has the potential to improve future lung cancer, COPD and pulmonary fibrosis patient outcomes.

Genome-wide association studies identify locus on 6p21 influencing lung function in the Korean population

BACKGROUND AND OBJECTIVE

Loss of lung function is an important chronic obstructive pulmonary disease phenotype and decreased forced expiratory volume in 1 s (FEV₁) is an independent risk factor of morbidity and mortality. Genome-wide association studies (GWAS) identifying genetic variants underlying lung function have been performed mostly in Caucasian populations. In this study, we aimed to identify genetic variants influencing lung function in a Korean population.

METHODS

GWAS on lung function (FEV₁ and FEV₁/forced vital capacity (FVC) ratio) were performed in two cohort studies. A population-based cohort, the Korean Association Resource phase 3 (KARE3) (6223 subjects), served as a discovery set. The replication analysis was performed in a family-based cohort, the Healthy Twin Study (HTS; 2730 subjects). Dense single-nucleotide polymorphism array data from each study were imputed and used for genetic analysis.

RESULTS

At the discovery phase, variants in 6p21 and 17q24 showed the strongest association with FEV₁/FVC ratio and FEV₁. Several variants in FAM13A on 4q22 locus exhibited positive association with FEV₁/FVC ratio. In the replication set, PPT2 in the 6p21 region showed significant association with lung function in the HTS, although the 4q22 locus and the 17q24 locus were not replicated.

CONCLUSIONS

We identified that PPT2 on chromosome 6p21 is associated with loss of lung function in the Korean population.

Copy number variation of the beta-defensin genes in europeans: no supporting evidence for association with lung function, chronic obstructive pulmonary disease or asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02-1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72-1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.

Heritability of chronic obstructive pulmonary disease and related phenotypes in smokers

RATIONALE

Previous studies of chronic obstructive pulmonary disease (COPD) have suggested that genetic factors play an important role in the development of disease. However, single-nucleotide polymorphisms that are associated with COPD in genome-wide association studies have been shown to account for only a small percentage of the genetic variance in phenotypes of COPD, such as spirometry and imaging variables. These phenotypes are highly predictive of disease, and family studies have shown that spirometric phenotypes are heritable.

OBJECTIVES

To assess the heritability and coheritability of four major COPD-related phenotypes (measurements of FEV1, FEV1/FVC, percent emphysema, and percent gas trapping), and COPD affection status in smokers of non-Hispanic white and African American descent using a population design.

METHODS

Single-nucleotide polymorphisms from genome-wide association studies chips were used to calculate the relatedness of pairs of individuals and a mixed model was adopted to estimate genetic variance and covariance.

MEASUREMENTS AND MAIN RESULTS

In the non-Hispanic whites, estimated heritabilities of FEV1 and FEV1/FVC were both about 37%, consistent with estimates in the literature from family-based studies. For chest computed tomography scan phenotypes, estimated heritabilities were both close to 25%. Heritability of COPD affection status was estimated as 37.7% in both populations.

CONCLUSIONS

This study suggests that a large portion of the genetic risk of COPD is yet to be discovered and gives rationale for additional genetic studies of COPD. The estimates of coheritability (genetic covariance) for pairs of the phenotypes suggest considerable overlap of causal genetic loci.

Systemic soluble receptor for advanced glycation endproducts is a biomarker of emphysema and associated with AGER genetic variants in patients with chronic obstructive pulmonary disease

RATIONALE

Emphysema in chronic obstructive pulmonary disease (COPD) can be characterized by high-resolution chest computed tomography (HRCT); however, the repeated use of HRCT is limited because of concerns regarding radiation exposure and cost.

OBJECTIVES

To evaluate biomarkers associated with emphysema and COPD-related clinical characteristics, and to assess the relationships of soluble receptor for advanced glycation endproducts (sRAGE), a candidate systemic biomarker identified in this study, with single-nucleotide polymorphisms (SNPs) in the gene coding for RAGE (AGER locus) and with clinical characteristics.

METHODS

Circulating levels of 111 biomarkers were analyzed for association with clinical characteristics in 410 patients with COPD enrolled in the TESRA study. sRAGE was also measured in the ECLIPSE cohort in 1,847 patients with COPD, 298 smokers and 204 nonsmokers. The association between 21 SNPs in the AGER locus with sRAGE levels and clinical characteristics was also investigated.

MEASUREMENTS AND MAIN RESULTS

sRAGE was identified as a biomarker of diffusing capacity of carbon monoxide and lung density in the TESRA cohort. In the ECLIPSE cohort, lower sRAGE levels were associated with increased emphysema, increased Global Initiative for Chronic Obstructive Lung Disease stage, and COPD disease status. The associations with emphysema in both cohorts remained significant after covariate adjustment (P < 0.0001). One SNP in the AGER locus, rs2070600, was associated with circulating sRAGE levels both in TESRA (P = 0.0014) and ECLIPSE (7.07 × 10(-16)), which exceeded genome-wide significance threshold. Another SNP (rs2071288) was also associated with sRAGE levels (P = 0.01) and diffusing capacity of carbon monoxide (P = 0.01) in the TESRA study.

CONCLUSIONS

Lower circulating sRAGE levels are associated with emphysema severity and genetic polymorphisms in the AGER locus are associated with systemic sRAGE levels. Clinical trial registered with www.clinicaltrials.gov (NCT 00413205 and NCT 00292552).

Gene-environment interaction effects on lung function- a genome-wide association study within the Framingham heart study

BACKGROUND

Previous studies in occupational exposure and lung function have focused only on the main effect of occupational exposure or genetics on lung function. Some disease-susceptible genes may be missed due to their low marginal effects, despite potential involvement in the disease process through interactions with the environment. Through comprehensive genome-wide gene-environment interaction studies, we can uncover these susceptibility genes. Our objective in this study was to explore gene by occupational exposure interaction effects on lung function using both the individual SNPs approach and the genetic network approach.

METHODS

The study population comprised the Offspring Cohort and the Third Generation from the Framingham Heart Study. We used forced expiratory volume in one second (FEV1) and ratio of FEV1 to forced vital capacity (FVC) as outcomes. Occupational exposures were classified using a population-specific job exposure matrix. We performed genome-wide gene-environment interaction analysis, using the Affymetrix 550 K mapping array for genotyping. A linear regression-based generalized estimating equation was applied to account for within-family relatedness. Network analysis was conducted using results from single-nucleotide polymorphism (SNP)-level analyses and from gene expression study results.

RESULTS

There were 4,785 participants in total. SNP-level analysis and network analysis identified SNP rs9931086 (P(interaction) =1.16 × 10(-7)) in gene SLC38A8, which may significantly modify the effects of occupational exposure on FEV1. Genes identified from the network analysis included CTLA-4, HDAC, and PPAR-alpha.

CONCLUSIONS

Our study implies that SNP rs9931086 in SLC38A8 and genes CTLA-4, HDAC, and PPAR-alpha, which are related to inflammatory processes, may modify the effect of occupational exposure on lung function.

Vitamin D-responsive SGPP2 variants associated with lung cell expression and lung function

Background

Vitamin D is associated with lung health in epidemiologic studies, but mechanisms mediating observed associations are poorly understood. This study explores mechanisms for an effect of vitamin D in lung through an in vivo gene expression study, an expression quantitative trait loci (eQTL) analysis in lung tissue, and a population-based cohort study of sequence variants.

Methods

Microarray analysis investigated the association of gene expression in small airway epithelial cells with serum 25(OH)D in adult non-smokers. Sequence variants in candidate genes identified by the microarray were investigated in a lung tissue eQTL database, and also in relation to cross-sectional pulmonary function in the Health, Aging, and Body Composition (Health ABC) study, stratified by race, with replication in the Framingham Heart Study (FHS).

Results

13 candidate genes had significant differences in expression by serum 25(OH)D (nominal p < 0.05), and a genome-wide significant eQTL association was detected for SGPP2. In Health ABC, SGPP2 SNPs were associated with FEV₁ in both European- and African-Americans, and the gene-level association was replicated in European-American FHS participants. SNPs in 5 additional candidate genes (DAPK1, FSTL1, KAL1, KCNS3, and RSAD2) were associated with FEV₁ in Health ABC participants.

Conclusions

SGPP2, a sphingosine-1-phosphate phosphatase, is a novel vitamin D-responsive gene associated with lung function. The identified associations will need to be followed up in further studies.

The Rotterdam Study: 2014 objectives and design update

The Rotterdam Study is a prospective cohort study ongoing since 1990 in the city of Rotterdam in The Netherlands. The study targets cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, oncological, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. The findings of the Rotterdam Study have been presented in over a 1,000 research articles and reports (see www.erasmus-epidemiology.nl/rotterdamstudy ). This article gives the rationale of the study and its design. It also presents a summary of the major findings and an update of the objectives and methods.

On the simultaneous association analysis of large genomic regions: a massive multi-locus association test

MOTIVATION

For samples of unrelated individuals, we propose a general analysis framework in which hundred thousands of genetic loci can be tested simultaneously for association with complex phenotypes. The approach is built on spatial-clustering methodology, assuming that genetic loci that are associated with the target phenotype cluster in certain genomic regions. In contrast to standard methodology for multilocus analysis, which has focused on the dimension reduction of the data, our multilocus association-clustering test profits from the availability of large numbers of genetic loci by detecting clusters of loci that are associated with the phenotype.

RESULTS

The approach is computationally fast and powerful, enabling the simultaneous association testing of large genomic regions. Even the entire genome or certain chromosomes can be tested simultaneously. Using simulation studies, the properties of the approach are evaluated. In an application to a genome-wide association study for chronic obstructive pulmonary disease, we illustrate the practical relevance of the proposed method by simultaneously testing all genotyped loci of the genome-wide association study and by testing each chromosome individually. Our findings suggest that statistical methodology that incorporates spatial-clustering information will be especially useful in whole-genome sequencing studies in which millions or billions of base pairs are recorded and grouped by genomic regions or genes, and are tested jointly for association.

AVAILABILITY AND IMPLEMENTATION

Implementation of the approach is available upon request.

Personalizing and targeting therapy for COPD: the role of molecular and clinical biomarkers

Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease characterized by persistent airflow limitation. It is the third leading cause of death worldwide, and there are currently no curative strategies for this disease. Many factors contribute to COPD susceptibility, progression and exacerbations. These include cigarette smoking, environmental and occupational pollutants, respiratory infections and comorbidities. As the clinical phenotypes of COPD are so variable, it has been difficult to devise an individualized treatment plan for patients with this complex chronic disease. This review will highlight how potential clinical, inflammatory, genomic and epigenomic biomarkers for COPD could be used to personalize treatment, leading to improved disease management and prevention for our patients.

Genetic variation in antioxidant enzymes, cigarette smoking, and longitudinal change in lung function

Antioxidant enzymes play an important role in the defense against oxidative stress in the lung and in the pathogenesis of chronic obstructive pulmonary disease (COPD). Sequence variation in genes encoding antioxidant enzymes may alter susceptibility to COPD by affecting longitudinal change in lung function in adults. We genotyped 384 sequence variants in 56 candidate genes in 1281 African American and 1794 European American elderly adults in the Health, Aging, and Body Composition study. Single-marker associations and gene-by-smoking interactions with rate of change in FEV₁ and FEV₁/FVC were evaluated using linear mixed-effects models, stratified by race/ethnicity. In European Americans, rs17883901 in GCLC was statistically significantly associated with rate of change in FEV₁/FVC; the recessive genotype (TT) was associated with a 0.9% per year steeper decline (P = 4.50 × 10(-5)). Statistically significant gene-by-smoking interactions were observed for variants in two genes in European Americans: the minor allele of rs2297765 in mGST3 attenuated the accelerated decline in FEV₁/FVC in smokers by 0.45% per year (P = 1.13 × 10(-4)); for participants with greater baseline smoking pack-years, the minor allele of rs2073192 in IDH3B was associated with an accelerated decline in FEV₁/FVC (P = 2.10 × 10(-4)). For both genes, nominally significant interactions (P < 0.01) were observed at the gene level in African Americans (P = 0.007 and 4.60 × 10(-4), respectively). Nominally significant evidence of association was observed for variants in SOD3 and GLRX2 in multiple analyses. This study identifies two novel genes associated with longitudinal lung function phenotypes in both African and European Americans and confirms a prior finding for GCLC. These findings suggest novel mechanisms and molecular targets for future research and advance the understanding of genetic determinants of lung function and COPD risk.

GSTCD and INTS12 regulation and expression in the human lung

Genome-Wide Association Study (GWAS) meta-analyses have identified a strong association signal for lung function, which maps to a region on 4q24 containing two oppositely transcribed genes: glutathione S-transferase, C-terminal domain containing (GSTCD) and integrator complex subunit 12 (INTS12). Both genes were found to be expressed in a range of human airway cell types. The promoter regions and transcription start sites were determined in mRNA from human lung and a novel splice variant was identified for each gene. We obtained the following evidence for GSTCD and INTS12 co-regulation and expression: (i) correlated mRNA expression was observed both via Q-PCR and in a lung expression quantitative trait loci (eQTL) study, (ii) induction of both GSTCD and INTS12 mRNA expression in human airway smooth muscle cells was seen in response to TGFβ1, (iii) a lung eQTL study revealed that both GSTCD and INTS12 mRNA levels positively correlate with percent predicted FEV1, and (iv) FEV1 GWAS associated SNPs in 4q24 were found to act as an eQTL for INTS12 in a number of tissues. In fixed sections of human lung tissue, GSTCD protein expression was ubiquitous, whereas INTS12 expression was predominantly in epithelial cells and pneumocytes. During human fetal lung development, GSTCD protein expression was observed to be highest at the earlier pseudoglandular stage (10-12 weeks) compared with the later canalicular stage (17-19 weeks), whereas INTS12 expression levels did not alter throughout these stages. Knowledge of the transcriptional and translational regulation and expression of GSTCD and INTS12 provides important insights into the potential role of these genes in determining lung function. Future work is warranted to fully define the functions of INTS12 and GSTCD.

Heritability of pulmonary function estimated from pedigree and whole-genome markers

Asthma and chronic obstructive pulmonary disease (COPD) are major worldwide health problems. Pulmonary function testing is a useful diagnostic tool for these diseases, and is known to be influenced by genetic and environmental factors. Previous studies have demonstrated that a substantial proportion of the variation in pulmonary function phenotypes can be explained by familial relationships. The availability of whole-genome single nucleotide polymorphism (SNP) data enables us to further evaluate the extent to which genetic factors account for variation in pulmonary function and to compare pedigree- to SNP-based estimates of heritability. Here, we employ methods developed in the animal breeding field to estimate the heritability of forced expiratory volume in one second (FEV1), forced vital capacity (FVC), and the ratio of these two measures (FEV1/FVC) among subjects in the Framingham Heart Study dataset. We compare heritability estimates based on pedigree-based relationships to those based on genome-wide SNPs. We find that, in a family-based study, estimates of heritability using SNP data are nearly identical to estimates based on pedigree information, and range from 0.50 for FEV1 to 0.66 for FEV1/FVC. Therefore, we conclude that genetic factors account for a sizable proportion of inter-individual differences in pulmonary function, and that estimates of heritability based on SNP data are nearly identical to estimates based on pedigree data. Finally, our findings suggest a higher heritability for FEV1/FVC compared to either FEV1 or FVC.

Gene network analysis in a pediatric cohort identifies novel lung function genes

Lung function is a heritable trait and serves as an important clinical predictor of morbidity and mortality for pulmonary conditions in adults, however, despite its importance, no studies have focused on uncovering pediatric-specific loci influencing lung function. To identify novel genetic determinants of pediatric lung function, we conducted a genome-wide association study (GWAS) of four pulmonary function traits, including FVC, FEV₁, FEV₁/FVC and FEF_25–75% in 1556 children. Further, we carried out gene network analyses for each trait including all SNPs with a P-value of <1.0×10⁻³ from the individual GWAS. The GWAS identified SNPs with notable trends towards association with the pulmonary function measures, including the previously described INTS12 locus association with FEV1 (p_meta = 1.41×10⁻⁷). The gene network analyses identified 34 networks of genes associated with pulmonary function variables in Caucasians. Of those, the glycoprotein gene network reached genome-wide significance for all four variables. P-value range p_meta = 6.29×10⁻⁴ - 2.80×10⁻⁸ on meta-analysis. In this study, we report on specific pathways that are significantly associated with pediatric lung function at genome-wide significance. In addition, we report the first loci associated with lung function in both pediatric Caucasian and African American populations.

Causal and synthetic associations of variants in the SERPINA gene cluster with alpha1-antitrypsin serum levels

Several infrequent genetic polymorphisms in the SERPINA1 gene are known to substantially reduce concentration of alpha1-antitrypsin (AAT) in the blood. Since low AAT serum levels fail to protect pulmonary tissue from enzymatic degradation, these polymorphisms also increase the risk for early onset chronic obstructive pulmonary disease (COPD). The role of more common SERPINA1 single nucleotide polymorphisms (SNPs) in respiratory health remains poorly understood.

We present here an agnostic investigation of genetic determinants of circulating AAT levels in a general population sample by performing a genome-wide association study (GWAS) in 1392 individuals of the SAPALDIA cohort.

Five common SNPs, defined by showing minor allele frequencies (MAFs) >5%, reached genome-wide significance, all located in the SERPINA gene cluster at 14q32.13. The top-ranking genotyped SNP rs4905179 was associated with an estimated effect of β = −0.068 g/L per minor allele (P = 1.20*10⁻¹²). But denser SERPINA1 locus genotyping in 5569 participants with subsequent stepwise conditional analysis, as well as exon-sequencing in a subsample (N = 410), suggested that AAT serum level is causally determined at this locus by rare (MAF<1%) and low-frequent (MAF 1–5%) variants only, in particular by the well-documented protein inhibitor S and Z (PI S, PI Z) variants. Replication of the association of rs4905179 with AAT serum levels in the Copenhagen City Heart Study (N = 8273) was successful (P<0.0001), as was the replication of its synthetic nature (the effect disappeared after adjusting for PI S and Z, P = 0.57). Extending the analysis to lung function revealed a more complex situation. Only in individuals with severely compromised pulmonary health (N = 397), associations of common SNPs at this locus with lung function were driven by rarer PI S or Z variants. Overall, our meta-analysis of lung function in ever-smokers does not support a functional role of common SNPs in the SERPINA gene cluster in the general population.

Low levels of alpha1-antitrypsin (AAT) in the blood are a well-established risk factor for accelerated loss in lung function and chronic obstructive pulmonary disease. While a few infrequent genetic polymorphisms are known to influence the serum levels of this enzyme, the role of common genetic variants has not been examined so far. The present genome-wide scan for associated variants in approximately 1400 Swiss inhabitants revealed a chromosomal locus containing the functionally established variants of AAT deficiency and variants previously associated with lung function and emphysema. We used dense genotyping of this genetic region in more than 5500 individuals and subsequent conditional analyses to unravel which of these associated variants contribute independently to the phenotype's variability. All associations of common variants could be attributed to the rarer functionally established variants, a result which was then replicated in an independent population-based Danish cohort. Hence, this locus represents a textbook example of how a large part of a trait's heritability can be hidden in infrequent genetic polymorphisms. The attempt to transfer these results to lung function furthermore suggests that effects of common variants in this genetic region in ever-smokers may also be explained by rarer variants, but only in individuals with hampered pulmonary health.

Genome-wide association study of lung function phenotypes in a founder population

BACKGROUND

Lung function is a long-term predictor of mortality and morbidity.

OBJECTIVE

We sought to identify single nucleotide polymorphisms (SNPs) associated with lung function.

METHODS

We performed a genome-wide association study (GWAS) of FEV1, forced vital capacity (FVC), and FEV1/FVC in 1144 Hutterites aged 6 to 89 years, who are members of a founder population of European descent. We performed least absolute shrinkage and selection operation regression to select the minimum set of SNPs that best predict FEV1/FVC in the Hutterites and used the GRAIL algorithm to mine the Gene Ontology database for evidence of functional connections between genes near the predictive SNPs.

RESULTS

Our GWAS identified significant associations between FEV1/FVC and SNPs at the THSD4-UACA-TLE3 locus on chromosome 15q23 (P = 5.7 × 10(-8) to 3.4 × 10(-9)). Nine SNPs at or near 4 additional loci had P < 10(-5) with FEV1/FVC. Only 2 SNPs were found with P < 10(-5) for FEV1 or FVC. We found nominal levels of significance with SNPs at 9 of the 27 previously reported loci associated with lung function measures. Among a predictive set of 80 SNPs, 6 loci were identified that had a significant degree of functional connectivity (GRAIL P < .05), including 3 clusters of β-defensin genes, 2 chemokine genes (CCL18 and CXCL12), and TNFRSF13B.

CONCLUSION

This study identifies genome-wide significant associations and replicates results of previous GWASs. Multimarker modeling implicated for the first time common variation in genes involved in antimicrobial immunity in airway mucosa that influences lung function.

HTR4 gene structure and altered expression in the developing lung

BACKGROUND

Meta-analyses of genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) spanning the 5-hydroxytryptamine receptor 4 (5-HT₄R) gene (HTR4) associated with lung function. The aims of this study were to i) investigate the expression profile of HTR4 in adult and fetal lung tissue and cultured airway cells, ii) further define HTR4 gene structure and iii) explore the potential functional implications of key SNPs using a bioinformatic approach.

METHODS

Following reverse transcription (RT)-PCR in human brain, 5' rapid amplification of cDNA ends (5' RACE) was used to examine the exonic structure of HTR4 at the 5' end. Quantitative (Q)-PCR was used to quantify HTR4 mRNA expression in total RNA from cultured airway cells and whole lung tissue. Publically available gene microarray data on fetal samples of estimated gestational age 7-22 weeks were mined for HTR4 expression. Immunohistochemistry (IHC; in adult and fetal lung tissue) and a radioligand binding assay (in cultured airway cells) were used to analyze 5-HT₄R protein expression.

RESULTS

IHC in adult lung, irrespective of the presence of chronic obstructive pulmonary disease (COPD), suggested low level expression of 5-HT₄R protein, which was most prominent in alveolar pneumocytes. There was evidence of differential 5-HT₄R protein levels during gestation in fetal lung, which was also evident in gene expression microarray data. HTR4 mRNA expression, assessed by Q-PCR, was <0.5% relative to brain in total adult lung tissue and in human airway smooth muscle (HASM) and bronchial epithelial cells (HBEC) derived from adult donors. Radioligand binding experiments also indicated that HBEC and HASM cells did not express a significant 5-HT₄R population. 5' RACE in brain identified a novel N-terminal variant, containing an extended N-terminal sequence. The functional significance of key HTR4 SNPs was investigated using the encyclopedia of DNA elements consortium (ENCODE) dataset. These analyses identified multiple alterations in regulatory motifs for transcription factors implicated in lung development, including Foxp1.

CONCLUSIONS

Taken together, these data suggest a role for HTR4 in lung development, which may at least in part explain the genetic association with lung function.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Translation initiation requires cell division cycle 123 (Cdc123) to facilitate biogenesis of the eukaryotic initiation factor 2 (eIF2)

The eukaryotic translation initiation factor 2 (eIF2) is central to the onset of protein synthesis and its modulation in response to physiological demands. eIF2, a heterotrimeric G-protein, is activated by guanine nucleotide exchange to deliver the initiator methionyl-tRNA to the ribosome. Here we report that assembly of the eIF2 complex in vivo depends on Cdc123, a cell proliferation protein conserved among eukaryotes. Mutations of CDC123 in budding yeast reduced the association of eIF2 subunits, diminished polysome levels, and increased GCN4 expression indicating that Cdc123 is critical for eIF2 activity. Cdc123 bound the unassembled eIF2γ subunit, but not the eIF2 complex, and the C-terminal domain III region of eIF2γ was both necessary and sufficient for Cdc123 binding. Alterations of the binding site revealed a strict correlation between Cdc123 binding, the biological function of eIF2γ, and its ability to assemble with eIF2α and eIF2β. Interestingly, high levels of Cdc123 neutralized the assembly defect and restored the biological function of an eIF2γ mutant. Moreover, the combined overexpression of eIF2 subunits rescued an otherwise inviable cdc123 deletion mutant. Thus, Cdc123 is a specific eIF2 assembly factor indispensable for the onset of protein synthesis. Human Cdc123 is encoded by a disease risk locus, and, therefore, eIF2 biogenesis control by Cdc123 may prove relevant for normal cell physiology and human health. This work identifies a novel step in the eukaryotic translation initiation pathway and assigns a biochemical function to a protein that is essential for growth and viability of eukaryotic cells.

Investigation of novel genes for lung function in children and their interaction with tobacco smoke exposure: a preliminary report

AIM

To replicate newly reported single nucleotide polymorphism (SNP) associations with adult lung function in a cohort of children and to investigate interactions with tobacco smoke exposure on lung function.

METHODS

Of 37 reported SNPs in adults, 21 were available in our genome-wide association study data set, either directly genotyped or as proxy SNPs. We tested for association with lung function (FEV1 , FVC, FEV1 /FVC%) in 8-year-old children and analysed interaction with tobacco smoke exposure both prenatally and/or during infancy, as well as at the age of 8 years.

RESULTS

SNPs in TNS1, ADAM19, THSD4 and ADCY2 showed significant association with lung function in children, although ADCY2 variants not in the expected direction. For example, variant allele A in THSD4 rs12899618 was associated with 50.2 mL higher FEV1 (p = 0.007) and variant allele C in ADAM19 rs2277027 with 1% unit lower FEV1 /FVC% (p = 0.03) per allele. DAAM2 rs2395730 showed interaction with current tobacco smoke exposure, with variant C allele associated with 1.3% unit decrease in FEV1 /FVC% in exposed children, but not in unexposed (p for interaction 0.04).

CONCLUSION

Our data replicate in children an association for TNS1, ADAM19, THSD4 and ADCY2 gene variants with lung function and suggest that interaction with tobacco smoke exposure may be of importance.

Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis

We performed a genome-wide association study in non-Hispanic white subjects with fibrotic idiopathic interstitial pneumonias (N=1616) and controls (N=4683); replication was assessed in 876 cases and 1890 controls. We confirmed association with TERT and MUC5B on chromosomes 5p15 and 11p15, respectively, the chromosome 3q26 region near TERC, and identified 7 novel loci (P_Meta = 2.4×10⁻⁸ to P_Meta = 1.1×10⁻¹⁹). The novel loci include FAM13A (4q22), DSP (6p24), OBFC1 (10q24), ATP11A (13q34), DPP9 (19p13), and chromosomal regions 7q22 and 15q14-15. Our results demonstrate that genes involved in host defense, cell-cell adhesion, and DNA repair contribute to the risk of fibrotic IIP.

Normal spirometry values in healthy elderly: the Rotterdam Study

Although many different reference values for spirometry are available from various studies, the elderly are usually underrepresented. Therefore, our objective was to assess reference values in a sample of healthy participants from a prospective population-based cohort study, including a large proportion of elderly. We included spirometry measurements of healthy, never smokers, from the Rotterdam Study and excluded participants with respiratory symptoms or prescriptions for respiratory medication. Age- and height-specific curves for the 5th (lower limit of normal) and the 50th (median) percentile of Forced Expiratory Volume in 1 s (FEV1), Forced Vital Capacity (FVC), and the ratio (FEV1/FVC) were calculated by quantile regression models. The group of healthy elderly study subjects consisted of 1,125 individuals, with a mean age of 68 years, ranging from 47 to 96 years of age. Sex stratified equations for the median and the lower limit of normal were calculated adjusted for age and height. In this study, we report age- and height-dependent reference limits for FEV1, FVC, and FEV1/FVC in a large population, and prediction equations for the lower limit of normal and median values for a sample containing a large proportion of healthy elderly.

Genome-wide association study identifies TH1 pathway genes associated with lung function in asthmatic patients

BACKGROUND

Recent meta-analyses of genome-wide association studies in general populations of European descent have identified 28 loci for lung function.

OBJECTIVE

We sought to identify novel lung function loci specifically for asthma and to confirm lung function loci identified in general populations.

METHODS

Genome-wide association studies of lung function (percent predicted FEV1 [ppFEV1], percent predicted forced vital capacity, and FEV1/forced vital capacity ratio) were performed in 4 white populations of European descent (n = 1544), followed by meta-analyses.

RESULTS

Seven of 28 previously identified lung function loci (HHIP, FAM13A, THSD4, GSTCD, NOTCH4-AGER, RARB, and ZNF323) identified in general populations were confirmed at single nucleotide polymorphism (SNP) levels (P < .05). Four of 32 loci (IL12A, IL12RB1, STAT4, and IRF2) associated with ppFEV1 (P < 10(-4)) belong to the TH1 or IL-12 cytokine family pathway. By using a linear additive model, these 4 TH1 pathway SNPs cumulatively explained 2.9% to 7.8% of the variance in ppFEV1 values in 4 populations (P = 3 × 10(-11)). Genetic scores of these 4 SNPs were associated with ppFEV1 values (P = 2 × 10(-7)) and the American Thoracic Society severe asthma classification (P = .005) in the Severe Asthma Research Program population. TH2 pathway genes (IL13, TSLP, IL33, and IL1RL1) conferring asthma susceptibility were not associated with lung function.

CONCLUSION

Genes involved in airway structure/remodeling are associated with lung function in both general populations and asthmatic subjects. TH1 pathway genes involved in anti-virus/bacterial infection and inflammation modify lung function in asthmatic subjects. Genes associated with lung function that might affect asthma severity are distinct from those genes associated with asthma susceptibility.

Pedigree and genotyping quality analyses of over 10,000 DNA samples from the Generation Scotland: Scottish Family Health Study

Background

Generation Scotland: Scottish Family Health Study (GS:SFHS) is a family-based biobank of 24,000 participants with rich phenotype and DNA available for genetic research. This paper describes the laboratory results from genotyping 32 single nucleotide polymorphisms (SNPs) on DNA from over 10,000 participants who attended GS:SFHS research clinics. The analysis described here was undertaken to test the quality of genetic information available to researchers. The success rate of each marker genotyped (call rate), minor allele frequency and adherence to Mendelian inheritance are presented. The few deviations in marker transmission in the 925 parent-child trios analysed were assessed as to whether they were likely to be miscalled genotypes, data or sample handling errors, or pedigree inaccuracies including non-paternity.

Methods

The first 10,450 GS:SFHS clinic participants who had spirometry and smoking data available and DNA extracted were selected. 32 SNPs were assayed, chosen as part of a replication experiment from a Genome-Wide Association Study meta-analysis of lung function.

Results

In total 325,336 genotypes were returned. The overall project pass rate (32 SNPs on 10,450 samples) was 97.29%. A total of 925 parent-child trios were assessed for transmission of the SNP markers, with 16 trios indicating evidence of inconsistency in the recorded pedigrees.

Conclusions

The Generation Scotland: Scottish Family Health Study used well-validated study methods and can produce good quality genetic data, with a low error rate. The GS:SFHS DNA samples are of high quality and the family groups were recorded and processed with accuracy during collection of the cohort.

A large scale gene-centric association study of lung function in newly-hired female cotton textile workers with endotoxin exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P_trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.

Respiratory diseases call for special attention from clinical and translational science

Respiratory diseases will become one of the top 3 leading causes of estimated mortality in 2020 and become about one third of total causes of estimated mortality. The journal of Translational Respiratory Medicine is a truly international, peer-reviewed journal devoted to the publication of articles on outstanding work with translational potentials between basic research and clinical application to understanding respiratory disease. Translational respiratory medicine will more focus on biomarker identification and validation in pulmonary diseases in combination with clinical informatics, targeted proteomics, bioinformatics, systems medicine, or mathematical science; on different translational strategies of cell-based therapy to clinical application to treat lung diseases; on targeted therapies in combination with personalized medicine; and on distant electronic medicine to monitor a large population of people’s health. Translational Respiratory Medicine is an additional but unique opportunity for scientists and clinicians who work on pulmonary diseases to publish their outstanding findings, initiative results, and critical and perceptive opinions in the journal.

Meta-analysis identifies a MECOM gene as a novel predisposing factor of osteoporotic fracture

BACKGROUND

Osteoporotic fracture (OF) as a clinical endpoint is a major complication of osteoporosis. To screen for OF susceptibility genes, we performed a genome-wide association study and carried out de novo replication analysis of an East Asian population.

METHODS

Association was tested using a logistic regression analysis. A meta-analysis was performed on the combined results using effect size and standard errors estimated for each study.

RESULTS

In a combined meta-analysis of a discovery cohort (288 cases and 1139 controls), three hospital based sets in replication stage I (462 cases and 1745 controls), and an independent ethnic group in replication stage II (369 cases and 560 for controls), we identified a new locus associated with OF (rs784288 in the MECOM gene) that showed genome-wide significance (p=3.59×10(-8); OR 1.39). RNA interference revealed that a MECOM knockdown suppresses osteoclastogenesis.

CONCLUSIONS

Our findings provide new insights into the genetic architecture underlying OF in East Asians.

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

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV(1)), and its ratio to forced vital capacity (FEV(1)/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV(1) and FEV(1)/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P(JMA = )5.00×10(-11)), HLA-DQB1 and HLA-DQA2 (smallest P(JMA = )4.35×10(-9)), and KCNJ2 and SOX9 (smallest P(JMA = )1.28×10(-8)) were associated with FEV(1)/FVC or FEV(1) in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

What can genetics tell us about the cause of fixed airflow obstruction?

Chronic obstructive pulmonary disease (COPD) is a major cause of chronic morbidity and mortality worldwide with smoking being the most important risk factor of the disease. However, lung function and COPD are known to also have a genetic component and a deeper knowledge of the genetic architecture of the disease could lead to further understanding of predisposition to COPD and also to development of new therapeutic interventions. Genetic linkage studies and candidate gene association studies have not provided evidence to convincingly identify the genes underlying lung function or COPD. However, recent large genome-wide association studies (GWAS) including tens of thousands of individuals have identified 26 variants at different loci in the human genome that show robust association with quantitative lung function measures in the general population. A growing number of these variants are being shown to be associated with COPD. Following the identification of these new lung function loci, the challenge now lies in refining the signals to identify the causative variants underlying the association signals and relating these signals to the molecular pathways that underlie lung function.

Opportunities and challenges of disease biomarkers: a new section in the Journal of Translational Medicine

Disease biomarkers are defined to diagnose various phases of diseases, monitor severities of diseases and responses to therapies, or predict prognosis of patients. Disease-specific biomarkers should benefit drug discovery and development, integrate multidisciplinary sciences, be validated by molecular imaging. The opportunities and challenges in biomarker development are emphasized and considered. The Journal of Translational Medicine opens a new Section of Disease Biomarkers to bridge identification and validation of gene or protein-based biomarkers, network biomarkers, dynamic network biomarkers in human diseases, patient phenotypes, and clinical applications. Disease biomarkers are also important for determining drug effects, target specificities and binding, dynamic metabolism and pharmacological kinetics, or toxicity profiles.

Identification of the BCAR1-CFDP1-TMEM170A locus as a determinant of carotid intima-media thickness and coronary artery disease risk

BACKGROUND

Carotid intima-media thickness (cIMT) is a widely accepted marker of subclinical atherosclerosis. To date, large-scale investigations of genetic determinants of cIMT are sparse.

METHODS AND RESULTS

To identify cIMT-associated genes and genetic variants, a discovery analysis using the Illumina 200K CardioMetabochip was conducted in 3430 subjects with detailed ultrasonographic determinations of cIMT from the IMPROVE (Carotid Intima Media Thickness [IMT] and IMT-Progression as Predictors of Vascular Events in a High Risk European Population) study. Segment-specific IMT measurements of common carotid, bifurcation, and internal carotid arteries, and composite IMT variables considering the whole carotid tree (IMT(mean), IMT(max), and IMT(mean-max)), were analyzed. A replication stage investigating 42 single-nucleotide polymorphisms for association with common carotid IMT was undertaken in 5 independent European cohorts (total n=11,590). A locus on chromosome 16 (lead single-nucleotide polymorphism rs4888378, intronic in CFDP1) was associated with cIMT at significance levels passing multiple testing correction at both stages (array-wide significant discovery P=6.75 × 10(-7) for IMT(max); replication P=7.24×10(-6) for common cIMT; adjustments for sex, age, and population substructure where applicable; minor allele frequency 0.43 and 0.41, respectively). The protective minor allele was associated with lower carotid plaque score in a replication cohort (P=0.04, n=2120) and lower coronary artery disease risk in 2 case-control studies of subjects with European ancestry (odds ratio [95% confidence interval] 0.83 [0.77-0.90], P=6.53 × 10(-6), n=13 591; and 0.95 [0.92-0.98], P=1.83 × 10(-4), n=82 297, respectively). Queries of human biobank data sets revealed associations of rs4888378 with nearby gene expression in vascular tissues (n=126-138).

CONCLUSIONS

This study identified rs4888378 in the BCAR1-CFDP1-TMEM170A locus as a novel genetic determinant of cIMT and coronary artery disease risk in individuals of European descent.

Gene-environment interactions and obesity traits among postmenopausal African-American and Hispanic women in the Women's Health Initiative SHARe Study

Genome-wide association studies (GWAS) of obesity measures have identified associations with single nucleotide polymorphisms (SNPs). However, no large-scale evaluation of gene-environment interactions has been performed. We conducted a search of gene-environment (G × E) interactions in post-menopausal African-American and Hispanic women from the Women's Health Initiative SNP Health Association Resource GWAS study. Single SNP linear regression on body mass index (BMI) and waist-to-hip circumference ratio (WHR) adjusted for multidimensional-scaling-derived axes of ancestry and age was run in race-stratified data with 871,512 SNPs available from African-Americans (N = 8,203) and 786,776 SNPs from Hispanics (N = 3,484). Tests of G × E interaction at all SNPs for recreational physical activity (m h/week), dietary energy intake (kcal/day), alcohol intake (categorical), cigarette smoking years, and cigarette smoking (ever vs. never) were run in African-Americans and Hispanics adjusted for ancestry and age at interview, followed by meta-analysis of G × E interaction terms. The strongest evidence for concordant G × E interactions in African-Americans and Hispanics was for smoking and marker rs10133840 (Q statistic P = 0.70, beta = -0.01, P = 3.81 × 10(-7)) with BMI as the outcome. The strongest evidence for G × E interaction within a cohort was in African-Americans with WHR as outcome for dietary energy intake and rs9557704 (SNP × kcal = -0.04, P = 2.17 × 10(-7)). No results exceeded the Bonferroni-corrected statistical significance threshold.