Rare loss of function variants in the hepatokine gene INHBE protect from abdominal obesity

Identifying genetic variants associated with lower waist-to-hip ratio can reveal new therapeutic targets for abdominal obesity. We use exome sequences from 362,679 individuals to identify genes associated with waist-to-hip ratio adjusted for BMI (WHRadjBMI), a surrogate for abdominal fat that is causally linked to type 2 diabetes and coronary heart disease. Predicted loss of function (pLOF) variants in INHBE associate with lower WHRadjBMI and this association replicates in data from AMP-T2D-GENES. INHBE encodes a secreted protein, the hepatokine activin E. In vitro characterization of the most common INHBE pLOF variant in our study, indicates an in-frame deletion resulting in a 90% reduction in secreted protein levels. We detect associations with lower WHRadjBMI for variants in ACVR1C, encoding an activin receptor, further highlighting the involvement of activins in regulating fat distribution. These findings highlight activin E as a potential therapeutic target for abdominal obesity, a phenotype linked to cardiometabolic disease.

Rare coding variants in DNA damage repair genes associated with timing of natural menopause

The age of menopause is associated with fertility and disease risk, and its genetic control is of great interest. We use whole-exome sequences from 132,370 women in the UK Biobank to test for associations between rare damaging variants and age at natural menopause. Rare damaging variants in five genes are significantly associated with menopause: CHEK2 (p = 3.3 × 10⁻⁵¹), DCLRE1A (p = 8.4 × 10⁻¹³), and HELB (p = 5.7 × 10⁻⁷) with later menopause and TOP3A (p = 7.6 × 10⁻⁸) and CLPB (p = 8.1 × 10⁻⁷) with earlier menopause. Two additional genes are suggestive: RAD54L (p = 2.4 × 10⁻⁶) with later menopause and HROB (p = 2.9 × 10⁻⁶) with earlier menopause. In a follow-up analysis of repeated questionnaires in women who were initially premenopausal, CHEK2, TOP3A, and RAD54L genotypes are associated with subsequent menopause. Consistent with previous genome-wide association studies (GWASs), six of the seven genes are involved in the DNA damage repair pathway. Phenome-wide scans across 398,569 men and women revealed that in addition to known associations with cancers and blood cell counts, rare variants in CHEK2 are also associated with increased risk for uterine fibroids, polycystic ovary syndrome, and prostate hypertrophy; these associations are not shared with higher-penetrance breast cancer genes. Causal mediation analysis suggests that approximately 8% of the breast cancer risk conferred by CHEK2 pathogenic variants after menopause is mediated through delayed menopause.

Gene-level analysis of rare variants in 379,066 whole exome sequences identifies an association of GIGYF1 loss of function with type 2 diabetes

Sequencing of large cohorts offers an unprecedented opportunity to identify rare genetic variants and to find novel contributors to human disease. We used gene-based collapsing tests to identify genes associated with glucose, HbA1c and type 2 diabetes (T2D) diagnosis in 379,066 exome-sequenced participants in the UK Biobank. We identified associations for variants in GCK, HNF1A and PDX1, which are known to be involved in Mendelian forms of diabetes. Notably, we uncovered novel associations for GIGYF1, a gene not previously implicated by human genetics in diabetes. GIGYF1 predicted loss of function (pLOF) variants associated with increased levels of glucose (0.77 mmol/L increase, p = 4.42 × 10^–12) and HbA1c (4.33 mmol/mol, p = 1.28 × 10^–14) as well as T2D diagnosis (OR = 4.15, p = 6.14 × 10^–11). Multiple rare variants contributed to these associations, including singleton variants. GIGYF1 pLOF also associated with decreased cholesterol levels as well as an increased risk of hypothyroidism. The association of GIGYF1 pLOF with T2D diagnosis replicated in an independent cohort from the Geisinger Health System. In addition, a common variant association for glucose and T2D was identified at the GIGYF1 locus. Our results highlight the role of GIGYF1 in regulating insulin signaling and protecting from diabetes.

GWAS of serum ALT and AST reveals an association of SLC30A10 Thr95Ile with hypermanganesemia symptoms

Understanding mechanisms of hepatocellular damage may lead to new treatments for liver disease, and genome-wide association studies (GWAS) of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) serum activities have proven useful for investigating liver biology. Here we report 100 loci associating with both enzymes, using GWAS across 411,048 subjects in the UK Biobank. The rare missense variant SLC30A10 Thr95Ile (rs188273166) associates with the largest elevation of both enzymes, and this association replicates in the DiscovEHR study. SLC30A10 excretes manganese from the liver to the bile duct, and rare homozygous loss of function causes the syndrome hypermanganesemia with dystonia-1 (HMNDYT1) which involves cirrhosis. Consistent with hematological symptoms of hypermanganesemia, SLC30A10 Thr95Ile carriers have increased hematocrit and risk of iron deficiency anemia. Carriers also have increased risk of extrahepatic bile duct cancer. These results suggest that genetic variation in SLC30A10 adversely affects more individuals than patients with diagnosed HMNDYT1.

A catalog of associations between rare coding variants and COVID-19 outcomes

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease-19 (COVID-19), a respiratory illness that can result in hospitalization or death. We investigated associations between rare genetic variants and seven COVID-19 outcomes in 543,213 individuals, including 8,248 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome-wide or when specifically focusing on (i) 14 interferon pathway genes in which rare deleterious variants have been reported in severe COVID-19 patients; (ii) 167 genes located in COVID-19 GWAS risk loci; or (iii) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, with results publicly browsable at https://rgc-covid19.regeneron.com.

Surviving sepsis campaign international guidelines for the management of septic shock and sepsis-associated organ dysfunction in children

Objectives

To develop evidence-based recommendations for clinicians caring for children (including infants, school-aged children, and adolescents) with septic shock and other sepsis-associated organ dysfunction.

Design

A panel of 49 international experts, representing 12 international organizations, as well as three methodologists and three public members was convened. Panel members assembled at key international meetings (for those panel members attending the conference), and a stand-alone meeting was held for all panel members in November 2018. A formal conflict-of-interest policy was developed at the onset of the process and enforced throughout. Teleconferences and electronic-based discussion among the chairs, co-chairs, methodologists, and group heads, as well as within subgroups, served as an integral part of the guideline development process.

Methods

The panel consisted of six subgroups: recognition and management of infection, hemodynamics and resuscitation, ventilation, endocrine and metabolic therapies, adjunctive therapies, and research priorities. We conducted a systematic review for each Population, Intervention, Control, and Outcomes question to identify the best available evidence, statistically summarized the evidence, and then assessed the quality of evidence using the Grading of Recommendations Assessment, Development, and Evaluation approach. We used the evidence-to-decision framework to formulate recommendations as strong or weak, or as a best practice statement. In addition, “in our practice” statements were included when evidence was inconclusive to issue a recommendation, but the panel felt that some guidance based on practice patterns may be appropriate.

Results

The panel provided 77 statements on the management and resuscitation of children with septic shock and other sepsis-associated organ dysfunction. Overall, six were strong recommendations, 49 were weak recommendations, and nine were best-practice statements. For 13 questions, no recommendations could be made; but, for 10 of these, “in our practice” statements were provided. In addition, 52 research priorities were identified.

Conclusions

A large cohort of international experts was able to achieve consensus regarding many recommendations for the best care of children with sepsis, acknowledging that most aspects of care had relatively low quality of evidence resulting in the frequent issuance of weak recommendations. Despite this challenge, these recommendations regarding the management of children with septic shock and other sepsis-associated organ dysfunction provide a foundation for consistent care to improve outcomes and inform future research.

Transthyretin-stabilising mutation T119M is not associated with protection against vascular disease or death in the UK Biobank

Background: Destabilised transthyretin (TTR) can result in the progressive, fatal disease transthyretin-mediated (ATTR) amyloidosis. A stabilising TTR mutation, T119M, is the basis for a therapeutic strategy to reduce destabilised TTR. Recently, T119M was associated with extended lifespan and lower risk of cerebrovascular disease in a Danish cohort. We aimed to determine whether this finding could be replicated in the UK Biobank.Methods: TTR T119M carriers were identified in the UK Biobank, a large prospective cohort of ∼500,000 individuals. Association between T119M genotype and inpatient diagnosis of vascular disease, cardiovascular disease, cerebrovascular disease, and mortality was analysed.Results: Frequency of T119M within the white UK Biobank population (n = 337,148) was 0.4%. Logistic regression comparing T119M carriers to non-carriers found no association between T119M and vascular disease (odds ratio [OR] = 1.08; p = .27), cardiovascular disease (OR = 1.08; p = .31), cerebrovascular disease (OR = 1.1; p = .42), or death (OR = 1.2; p = .06). Cox proportional hazards regression showed similar results (hazard ratio >1, p>.05). Age at death and vascular disease diagnosis were similar between T119M carriers and non-carriers (p = .12 and p = .38, respectively).Conclusions: There was no association between the TTR T119M genotype and risk of vascular disease or death in a large prospective cohort study, indicating that TTR tetramer stabilisation through T119M is not protective in this setting.

Unsupervised discovery of phenotype-specific multi-omics networks

MOTIVATION

Complex diseases often involve a wide spectrum of phenotypic traits. Better understanding of the biological mechanisms relevant to each trait promotes understanding of the etiology of the disease and the potential for targeted and effective treatment plans. There have been many efforts towards omics data integration and network reconstruction, but limited work has examined the incorporation of relevant (quantitative) phenotypic traits.

RESULTS

We propose a novel technique, sparse multiple canonical correlation network analysis (SmCCNet), for integrating multiple omics data types along with a quantitative phenotype of interest, and for constructing multi-omics networks that are specific to the phenotype. As a case study, we focus on miRNA-mRNA networks. Through simulations, we demonstrate that SmCCNet has better overall prediction performance compared to popular gene expression network construction and integration approaches under realistic settings. Applying SmCCNet to studies on chronic obstructive pulmonary disease (COPD) and breast cancer, we found enrichment of known relevant pathways (e.g. the Cadherin pathway for COPD and the interferon-gamma signaling pathway for breast cancer) as well as less known omics features that may be important to the diseases. Although those applications focus on miRNA-mRNA co-expression networks, SmCCNet is applicable to a variety of omics and other data types. It can also be easily generalized to incorporate multiple quantitative phenotype simultaneously. The versatility of SmCCNet suggests great potential of the approach in many areas.

AVAILABILITY AND IMPLEMENTATION

The SmCCNet algorithm is written in R, and is freely available on the web at https://cran.r-project.org/web/packages/SmCCNet/index.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mixed-model admixture mapping identifies smoking-dependent loci of lung function in African Americans

Admixture mapping has led to the discovery of many genes associated with differential disease risk by ancestry, highlighting the importance of ancestry-based approaches to association studies. However, the potential of admixture mapping in deciphering the interplay between genes and environment exposures has been seldom explored. Here we performed a genome-wide screening of local ancestry-smoking interactions for five spirometric lung function phenotypes in 3300 African Americans from the COPDGene study. To account for population structure and outcome heterogeneity across exposure groups, we developed a multi-component linear mixed model for mapping gene-environment interactions and empirically showed its robustness and increased power. When applied to the COPDGene study, our approach identified two 11p15.2-3 and 2q37 loci, exhibiting local ancestry-smoking interactions at genome-wide significant level, which would have been missed by standard single-nucleotide polymorphism analyses. These two loci harbor the PARVA and RAB17 genes previously recognized to be involved in smoking behavior. Overall, our study provides the first evidence for potential synergistic effects between African ancestry and smoking on pulmonary function, and underlines the importance of ethnic diversity in genetic studies.

Heme metabolism genes Downregulated in COPD Cachexia

INTRODUCTION

Cachexia contributes to increased mortality and reduced quality of life in Chronic Obstructive Pulmonary Disease (COPD) and may be associated with underlying gene expression changes. Our goal was to identify differential gene expression signatures associated with COPD cachexia in current and former smokers.

METHODS

We analyzed whole-blood gene expression data from participants with COPD in a discovery cohort (COPDGene, N = 400) and assessed replication (ECLIPSE, N = 114). To approximate the consensus definition using available criteria, cachexia was defined as weight-loss > 5% in the past 12 months or low body mass index (BMI) (< 20 kg/m2) and 1/3 criteria: decreased muscle strength (six-minute walk distance < 350 m), anemia (hemoglobin < 12 g/dl), and low fat-free mass index (FFMI) (< 15 kg/m2 among women and < 17 kg/m2 among men) in COPDGene. In ECLIPSE, cachexia was defined as weight-loss > 5% in the past 12 months or low BMI and 3/5 criteria: decreased muscle strength, anorexia, abnormal biochemistry (anemia or high c-reactive protein (> 5 mg/l)), fatigue, and low FFMI. Differential gene expression was assessed between cachectic and non-cachectic subjects, adjusting for age, sex, white blood cell counts, and technical covariates. Gene set enrichment analysis was performed using MSigDB.

RESULTS

The prevalence of COPD cachexia was 13.7% in COPDGene and 7.9% in ECLIPSE. Fourteen genes were differentially downregulated in cachectic versus non-cachectic COPD patients in COPDGene (FDR < 0.05) and ECLIPSE (FDR < 0.05).

DISCUSSION

Several replicated genes regulating heme metabolism were downregulated among participants with COPD cachexia. Impaired heme biosynthesis may contribute to cachexia development through free-iron buildup and oxidative tissue damage.

Common and Rare Variants Genetic Association Analysis of Cigarettes per Day Among Ever-Smokers in Chronic Obstructive Pulmonary Disease Cases and Controls

INTRODUCTION

Cigarette smoking is a major environmental risk factor for many diseases, including chronic obstructive pulmonary disease (COPD). There are shared genetic influences on cigarette smoking and COPD. Genetic risk factors for cigarette smoking in cohorts enriched for COPD are largely unknown.

METHODS

We performed genome-wide association analyses for average cigarettes per day (CPD) across the Genetic Epidemiology of COPD (COPDGene) non-Hispanic white (NHW) (n = 6659) and African American (AA) (n = 3260), GenKOLS (the Genetics of Chronic Obstructive Lung Disease) (n = 1671), and ECLIPSE (the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) (n = 1942) cohorts. In addition, we performed exome array association analyses across the COPDGene NHW and AA cohorts. We considered analyses across the entire cohort and stratified by COPD case-control status.

RESULTS

We identified genome-wide significant associations for CPD on chromosome 15q25 across all cohorts (lowest p = 1.78 × 10-15), except in the COPDGene AA cohort alone. Previously reported associations on chromosome 19 had suggestive and directionally consistent associations (RAB4, p = 1.95 × 10-6; CYP2A7, p = 7.50 × 10-5; CYP2B6, p = 4.04 × 10-4). When we stratified by COPD case-control status, single nucleotide polymorphisms on chromosome 15q25 were nominally associated with both NHW COPD cases (β = 0.11, p = 5.58 × 10-4) and controls (β = 0.12, p = 3.86 × 10-5) For the gene-based exome array association analysis of rare variants, there were no exome-wide significant associations. For these previously replicated associations, the most significant results were among COPDGene NHW subjects for CYP2A7 (p = 5.2 × 10-4).

CONCLUSIONS

In a large genome-wide association study of both common variants and a gene-based association of rare coding variants in ever-smokers, we found genome-wide significant associations on chromosome 15q25 with CPD for common variants, but not for rare coding variants. These results were directionally consistent among COPD cases and controls.

IMPLICATIONS

We examined both common and rare coding variants associated with CPD in a large population of heavy smokers with and without COPD of NHW and AA descent. We replicated genome-wide significant associations on chromosome 15q25 with CPD for common variants among NHW subjects, but not for rare variants. We demonstrated for the first time that common variants on chromosome 15q25 associated with CPD are similar among COPD cases and controls. Previously reported associations on chromosome 19 showed suggestive and directionally consistent associations among common variants (RAB4, CYP2A7, and CYP2B6) and for rare variants (CYP2A7) among COPDGene NHW subjects. Although the genetic effect sizes for these single nucleotide polymorphisms on chromosome 15q25 are modest, we show that this creates a substantial smoking burden over the lifetime of a smoker.

Surviving Sepsis Campaign International Guidelines for the Management of Septic Shock and Sepsis-Associated Organ Dysfunction in Children

OBJECTIVES

To develop evidence-based recommendations for clinicians caring for children (including infants, school-aged children, and adolescents) with septic shock and other sepsis-associated organ dysfunction.

DESIGN

A panel of 49 international experts, representing 12 international organizations, as well as three methodologists and three public members was convened. Panel members assembled at key international meetings (for those panel members attending the conference), and a stand-alone meeting was held for all panel members in November 2018. A formal conflict-of-interest policy was developed at the onset of the process and enforced throughout. Teleconferences and electronic-based discussion among the chairs, co-chairs, methodologists, and group heads, as well as within subgroups, served as an integral part of the guideline development process.

METHODS

The panel consisted of six subgroups: recognition and management of infection, hemodynamics and resuscitation, ventilation, endocrine and metabolic therapies, adjunctive therapies, and research priorities. We conducted a systematic review for each Population, Intervention, Control, and Outcomes question to identify the best available evidence, statistically summarized the evidence, and then assessed the quality of evidence using the Grading of Recommendations Assessment, Development, and Evaluation approach. We used the evidence-to-decision framework to formulate recommendations as strong or weak, or as a best practice statement. In addition, "in our practice" statements were included when evidence was inconclusive to issue a recommendation, but the panel felt that some guidance based on practice patterns may be appropriate.

RESULTS

The panel provided 77 statements on the management and resuscitation of children with septic shock and other sepsis-associated organ dysfunction. Overall, six were strong recommendations, 52 were weak recommendations, and nine were best-practice statements. For 13 questions, no recommendations could be made; but, for 10 of these, "in our practice" statements were provided. In addition, 49 research priorities were identified.

CONCLUSIONS

A large cohort of international experts was able to achieve consensus regarding many recommendations for the best care of children with sepsis, acknowledging that most aspects of care had relatively low quality of evidence resulting in the frequent issuance of weak recommendations. Despite this challenge, these recommendations regarding the management of children with septic shock and other sepsis-associated organ dysfunction provide a foundation for consistent care to improve outcomes and inform future research.

Integrative Genomics Analysis Identifies ACVR1B as a Candidate Causal Gene of Emphysema Distribution

Genome-wide association studies (GWAS) have identified multiple associations with emphysema apicobasal distribution (EABD), but the biological functions of these variants are unknown. To characterize the functions of EABD-associated variants, we integrated GWAS results with 1) expression quantitative trait loci (eQTL) from the Genotype Tissue Expression (GTEx) project and subjects in the COPDGene (Genetic Epidemiology of COPD) study and 2) cell type epigenomic marks from the Roadmap Epigenomics project. On the basis of these analyses, we selected a variant near ACVR1B (activin A receptor type 1B) for functional validation. SNPs from 168 loci with P values less than 5 × 10^-5 in the largest GWAS meta-analysis of EABD were analyzed. Eighty-four loci overlapped eQTL, with 12 of these loci showing greater than 80% likelihood of harboring a single, shared GWAS and eQTL causal variant. Seventeen cell types were enriched for overlap between EABD loci and Roadmap Epigenomics marks (permutation P < 0.05), with the strongest enrichment observed in CD4⁺, CD8⁺, and regulatory T cells. We selected a putative causal variant, rs7962469, associated with ACVR1B expression in lung tissue for additional functional investigation, and reporter assays confirmed allele-specific regulatory activity for this variant in human bronchial epithelial and Jurkat immune cell lines. ACVR1B expression levels exhibit a nominally significant association with emphysema distribution. EABD-associated loci are preferentially enriched in regulatory elements of multiple cell types, most notably T-cell subsets. Multiple EABD loci colocalize to regulatory elements that are active across multiple tissues and cell types, and functional analyses confirm the presence of an EABD-associated functional variant that regulates ACVR1B expression, indicating that transforming growth factor-β signaling plays a role in the EABD phenotype. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Identification of an emphysema-associated genetic variant near TGFB2 with regulatory effects in lung fibroblasts

Murine studies have linked TGF-β signaling to emphysema, and human genome-wide association studies (GWAS) studies of lung function and COPD have identified associated regions near genes in the TGF-β superfamily. However, the functional regulatory mechanisms at these loci have not been identified. We performed the largest GWAS of emphysema patterns to date, identifying 10 GWAS loci including an association peak spanning a 200 kb region downstream from TGFB2. Integrative analysis of publicly available eQTL, DNaseI, and chromatin conformation data identified a putative functional variant, rs1690789, that may regulate TGFB2 expression in human fibroblasts. Using chromatin conformation capture, we confirmed that the region containing rs1690789 contacts the TGFB2 promoter in fibroblasts, and CRISPR/Cas-9 targeted deletion of a ~ 100 bp region containing rs1690789 resulted in decreased TGFB2 expression in primary human lung fibroblasts. These data provide novel mechanistic evidence linking genetic variation affecting the TGF-β pathway to emphysema in humans.

Systemic Markers of Adaptive and Innate Immunity Are Associated with Chronic Obstructive Pulmonary Disease Severity and Spirometric Disease Progression

The progression of chronic obstructive pulmonary disease (COPD) is associated with marked alterations in circulating immune cell populations, but no studies have characterized alterations in these cell types across the full spectrum of lung function impairment in current and former smokers. In 6,299 subjects from the COPDGene and ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) studies, we related Coulter blood counts and proportions to cross-sectional forced expiratory volume in 1 second (FEV₁), adjusting for current smoking status. We also related cell count measures to 3-year change in FEV₁ in ECLIPSE subjects. In a subset of subjects with blood gene expression data, we used cell type deconvolution methods to infer the proportions of immune cell subpopulations, and we related these to COPD clinical status. We observed that FEV₁ levels are positively correlated with lymphocytes and negatively correlated with myeloid populations, such as neutrophils and monocytes. In multivariate models, absolute cell counts and proportions were associated with cross-sectional FEV₁, and lymphocytes, monocytes, and eosinophil counts were predictive of 3-year change in lung function. Using cell type deconvolution to study immune cell subpopulations, we observed that subjects with COPD had a lower proportion of CD4⁺ resting memory cells and naive B cells compared with smokers without COPD. Alterations in circulating immune cells in COPD support a mixed pattern of lymphocyte suppression and an enhanced myeloid cell immune response. Cell counts and proportions contribute independent information to models predicting lung function, suggesting a critical role for immune response in long-term COPD outcomes. Cell type deconvolution is a promising method for immunophenotyping in large cohorts.

Whole-Genome Sequencing in Severe Chronic Obstructive Pulmonary Disease

Genome-wide association studies have identified common variants associated with chronic obstructive pulmonary disease (COPD). Whole-genome sequencing (WGS) offers comprehensive coverage of the entire genome, as compared with genotyping arrays or exome sequencing. We hypothesized that WGS in subjects with severe COPD and smoking control subjects with normal pulmonary function would allow us to identify novel genetic determinants of COPD. We sequenced 821 patients with severe COPD and 973 control subjects from the COPDGene and Boston Early-Onset COPD studies, including both non-Hispanic white and African American individuals. We performed single-variant and grouped-variant analyses, and in addition, we assessed the overlap of variants between sequencing- and array-based imputation. Our most significantly associated variant was in a known region near HHIP (combined P = 1.6 × 10^-9); additional variants approaching genome-wide significance included previously described regions in CHRNA5, TNS1, and SERPINA6/SERPINA1 (the latter in African American individuals). None of our associations were clearly driven by rare variants, and we found minimal evidence of replication of genes identified by previously reported smaller sequencing studies. With WGS, we identified more than 20 million new variants, not seen with imputation, including more than 10,000 of potential importance in previously identified COPD genome-wide association study regions. WGS in severe COPD identifies a large number of potentially important functional variants, with the strongest associations being in known COPD risk loci, including HHIP and SERPINA1. Larger sample sizes will be needed to identify associated variants in novel regions of the genome.

Author Correction: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

RNA-sequencing across three matched tissues reveals shared and tissue-specific gene expression and pathway signatures of COPD

BACKGROUND

Multiple gene expression studies have been performed separately in peripheral blood, lung, and airway tissues to study COPD. We performed RNA-sequencing gene expression profiling of large-airway epithelium, alveolar macrophage and peripheral blood samples from the same subset of COPD cases and controls from the COPDGene study who underwent bronchoscopy at a single center. Using statistical and gene set enrichment approaches, we sought to improve the understanding of COPD by studying gene sets and pathways across these tissues, beyond the individual genomic determinants.

METHODS

We performed differential expression analysis using RNA-seq data obtained from 63 samples from 21 COPD cases and controls (includes four non-smokers) via the R package DESeq2. We tested associations between gene expression and variables related to lung function, smoking history, and CT scan measures of emphysema and airway disease. We examined the correlation of differential gene expression across the tissues and phenotypes, hypothesizing that this would reveal preserved and private gene expression signatures. We performed gene set enrichment analyses using curated databases and findings from prior COPD studies to provide biological and disease relevance.

RESULTS

The known smoking-related genes CYP1B1 and AHRR were among the top differential expression results for smoking status in the large-airway epithelium data. We observed a significant overlap of genes primarily across large-airway and macrophage results for smoking and airway disease phenotypes. We did not observe specific genes differentially expressed in all three tissues for any of the phenotypes. However, we did observe hemostasis and immune signaling pathways in the overlaps across all three tissues for emphysema, and amyloid and telomere-related pathways for smoking. In peripheral blood, the emphysema results were enriched for B cell related genes previously identified in lung tissue studies.

CONCLUSIONS

Our integrative analyses across COPD-relevant tissues and prior studies revealed shared and tissue-specific disease biology. These replicated and novel findings in the airway and peripheral blood have highlighted candidate genes and pathways for COPD pathogenesis.

New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function-associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD.

Assessing pleiotropy and mediation in genetic loci associated with chronic obstructive pulmonary disease

Genetic association studies have increasingly recognized variant effects on multiple phenotypes. Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease with environmental and genetic causes. Multiple genetic variants have been associated with COPD, many of which show significant associations to additional phenotypes. However, it is unknown if these associations represent biological pleiotropy or if they exist through correlation of related phenotypes ("mediated pleiotropy"). Using 6,670 subjects from the COPDGene study, we describe the association of known COPD susceptibility loci with other COPD-related phenotypes and distinguish if these act directly on the phenotypes (i.e., biological pleiotropy) or if the association is due to correlation (i.e., mediated pleiotropy). We identified additional associated phenotypes for 13 of 25 known COPD loci. Tests for pleiotropy between genotype and associated outcomes were significant for all loci. In cases of significant pleiotropy, we performed mediation analysis to test if SNPs had a direct association to phenotype. Most loci showed a mediated effect through the hypothesized causal pathway. However, many loci also had direct associations, suggesting causal explanations (i.e., emphysema leading to reduced lung function) are incomplete. Our results highlight the high degree of pleiotropy in complex disease-associated loci and provide novel insights into the mechanisms underlying COPD.

Inferring disease risk genes from sequencing data in multiplex pedigrees through sharing of rare variants

We previously demonstrated how sharing of rare variants (RVs) in distant affected relatives can be used to identify variants causing a complex and heterogeneous disease. This approach tested whether single RVs were shared by all sequenced affected family members. However, as with other study designs, joint analysis of several RVs (e.g., within genes) is sometimes required to obtain sufficient statistical power. Further, phenocopies can lead to false negatives for some causal RVs if complete sharing among affected is required. Here, we extend our methodology (Rare Variant Sharing, RVS) to address these issues. Specifically, we introduce gene-based analyses, a partial sharing test based on RV sharing probabilities for subsets of affected relatives and a haplotype-based RV definition. RVS also has the desirable feature of not requiring external estimates of variant frequency or control samples, provides functionality to assess and address violations of key assumptions, and is available as open source software for genome-wide analysis. Simulations including phenocopies, based on the families of an oral cleft study, revealed the partial and complete sharing versions of RVS achieved similar statistical power compared with alternative methods (RareIBD and the Gene-Based Segregation Test), and had superior power compared with the pedigree Variant Annotation, Analysis, and Search Tool (pVAAST) linkage statistic. In studies of multiplex cleft families, analysis of rare single nucleotide variants in the exome of 151 affected relatives from 54 families revealed no significant excess sharing in any one gene, but highlighted different patterns of sharing revealed by the complete and partial sharing tests.

Elevated circulating MMP-9 is linked to increased COPD exacerbation risk in SPIROMICS and COPDGene

BACKGROUND

Matrix metalloprotease 9 (MMP-9) is associated with inflammation and lung remodeling in chronic obstructive pulmonary disease (COPD). We hypothesized that elevated circulating MMP-9 represents a potentially novel biomarker that identifies a subset of individuals with COPD with an inflammatory phenotype who are at increased risk for acute exacerbation (AECOPD).

METHODS

We analyzed Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) and Genetic Epidemiology of COPD (COPDGene) cohorts for which baseline and prospective data were available. Elevated MMP-9 was defined based on >95th percentile plasma values from control (non-COPD) sample in SPIROMICS. COPD subjects were classified as having elevated or nonelevated MMP-9. Logistic, Poisson, and Kaplan-Meier analyses were used to identify associations with prospective AECOPD in both cohorts.

RESULTS

Elevated MMP-9 was present in 95/1,053 (9%) of SPIROMICS and 41/140 (29%) of COPDGene participants with COPD. COPD subjects with elevated MMP-9 had a 13%-16% increased absolute risk for AECOPD and a higher median (interquartile range; IQR) annual AECOPD rate (0.33 [0-0.74] versus 0 [0-0.80] events/year and 0.9 [0.5-2] versus 0.5 [0-1.4] events/year for SPIROMICS and COPDGene, respectively). In adjusted models within each cohort, elevated MMP-9 was associated with increased odds (odds ratio [OR], 1.71; 95%CI, 1.00-2.90; and OR, 3.03; 95%CI, 1.02-9.01), frequency (incidence rate ratio [IRR], 1.45; 95%CI, 1.23-1.7; and IRR, 1.24; 95%CI, 1.03-1.49), and shorter time-to-first AECOPD (21.7 versus 31.7 months and 14 versus 21 months) in SPIROMICS and COPDGene, respectively.

CONCLUSIONS

Elevated MMP-9 was independently associated with AECOPD risk in 2 well-characterized COPD cohorts. These findings provide evidence for MMP-9 as a prognostic biomarker and potential therapeutic target in COPD.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT01969344 (SPIROMICS) and NCT00608764 (COPDGene).

FUNDING

This work was funded by K08 HL123940 to JMW; R01HL124233 to PJC; Merit Review I01 CX000911 to JLC; R01 (R01HL102371, R01HL126596) and VA Merit (I01BX001756) to AG. SPIROMICS (Subpopulations and Intermediate Outcomes in COPD Study) is funded by contracts from the NHLBI (HHSN268200900013C, HHSN268200900014C,HHSN268200900015C HHSN268200900016C, HHSN268200900017C, HHSN268200900018C, HHSN268200900019C, and HHSN268200900020C) and a grant from the NIH/NHLBI (U01 HL137880), and supplemented by contributions made through the Foundation for the NIH and the COPD Foundation from AstraZeneca/MedImmune; Bayer; Bellerophon Therapeutics; Boehringer-Ingelheim Pharmaceuticals Inc.; Chiesi Farmaceutici; Forest Research Institute Inc.; GlaxoSmithKline; Grifols Therapeutics Inc.; Ikaria Inc.; Novartis Pharmaceuticals Corporation; Nycomed GmbH; ProterixBio; Regeneron Pharmaceuticals Inc.; Sanofi; Sunovion; Takeda Pharmaceutical Company; and Theravance Biopharma and Mylan. COPDGene is funded by the NHLBI (R01 HL089897 and R01 HL089856) and by the COPD Foundation through contributions made to an Industry Advisory Board composed of AstraZeneca, Boehringer Ingelheim, GlaxoSmithKline, Novartis, Pfizer, Siemens, and Sunovion.

Whole exome sequencing analysis in severe chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD), one of the leading causes of death worldwide, is substantially influenced by genetic factors. Alpha-1 antitrypsin deficiency demonstrates that rare coding variants of large effect can influence COPD susceptibility. To identify additional rare coding variants in patients with severe COPD, we conducted whole exome sequencing analysis in 2543 subjects from two family-based studies (Boston Early-Onset COPD Study and International COPD Genetics Network) and one case-control study (COPDGene). Applying a gene-based segregation test in the family-based data, we identified significant segregation of rare loss of function variants in TBC1D10A and RFPL1 (P-value < 2x10-6), but were unable to find similar variants in the case-control study. In single-variant, gene-based and pathway association analyses, we were unable to find significant findings that replicated or were significant in meta-analysis. However, we found that the top results in the two datasets were in proximity to each other in the protein-protein interaction network (P-value = 0.014), suggesting enrichment of these results for similar biological processes. A network of these association results and their neighbors was significantly enriched in the transforming growth factor beta-receptor binding and cilia-related pathways. Finally, in a more detailed examination of candidate genes, we identified individuals with putative high-risk variants, including patients harboring homozygous mutations in genes associated with cutis laxa and Niemann-Pick Disease Type C. Our results likely reflect heterogeneity of genetic risk for COPD along with limitations of statistical power and functional annotation, and highlight the potential of network analysis to gain insight into genetic association studies.

Longitudinal Modeling of Lung Function Trajectories in Smokers with and without Chronic Obstructive Pulmonary Disease

RATIONALE

The relationship between longitudinal lung function trajectories, chest computed tomography (CT) imaging, and genetic predisposition to chronic obstructive pulmonary disease (COPD) has not been explored.

OBJECTIVES

1) To model trajectories using a data-driven approach applied to longitudinal data spanning adulthood in the Normative Aging Study (NAS), and 2) to apply these models to demographically similar subjects in the COPDGene (Genetic Epidemiology of COPD) Study with detailed phenotypic characterization including chest CT.

METHODS

We modeled lung function trajectories in 1,060 subjects in NAS with a median follow-up time of 29 years. We assigned 3,546 non-Hispanic white males in COPDGene to these trajectories for further analysis. We assessed phenotypic and genetic differences between trajectories and across age strata.

MEASUREMENTS AND MAIN RESULTS

We identified four trajectories in NAS with differing levels of maximum lung function and rate of decline. In COPDGene, 617 subjects (17%) were assigned to the lowest trajectory and had the greatest radiologic burden of disease (P < 0.01); 1,283 subjects (36%) were assigned to a low trajectory with evidence of airway disease preceding emphysema on CT; 1,411 subjects (40%) and 237 subjects (7%) were assigned to the remaining two trajectories and tended to have preserved lung function and negligible emphysema. The genetic contribution to these trajectories was as high as 83% (P = 0.02), and membership in lower lung function trajectories was associated with greater parental histories of COPD, decreased exercise capacity, greater dyspnea, and more frequent COPD exacerbations.

CONCLUSIONS

Data-driven analysis identifies four lung function trajectories. Trajectory membership has a genetic basis and is associated with distinct lung structural abnormalities.

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.

Association Between the New York Sepsis Care Mandate and In-Hospital Mortality for Pediatric Sepsis

IMPORTANCE

The death of a pediatric patient with sepsis motivated New York to mandate statewide sepsis treatment in 2013. The mandate included a 1-hour bundle of blood cultures, broad-spectrum antibiotics, and a 20-mL/kg intravenous fluid bolus. Whether completing the bundle elements within 1 hour improves outcomes is unclear.

OBJECTIVE

To determine the risk-adjusted association between completing the 1-hour pediatric sepsis bundle and individual bundle elements with in-hospital mortality.

DESIGN, SETTINGS, AND PARTICIPANTS

Statewide cohort study conducted from April 1, 2014, to December 31, 2016, in emergency departments, inpatient units, and intensive care units across New York State. A total of 1179 patients aged 18 years and younger with sepsis and septic shock reported to the New York State Department of Health who had a sepsis protocol initiated were included.

EXPOSURES

Completion of a 1-hour sepsis bundle within 1 hour compared with not completing the 1-hour sepsis bundle within 1 hour.

MAIN OUTCOMES AND MEASURES

Risk-adjusted in-hospital mortality.

RESULTS

Of 1179 patients with sepsis reported at 54 hospitals (mean [SD] age, 7.2 [6.2] years; male, 54.2%; previously healthy, 44.5%; diagnosed as having shock, 68.8%), 139 (11.8%) died. The entire sepsis bundle was completed in 1 hour in 294 patients (24.9%). Antibiotics were administered to 798 patients (67.7%), blood cultures were obtained in 740 patients (62.8%), and the fluid bolus was completed in 548 patients (46.5%) within 1 hour. Completion of the entire bundle within 1 hour was associated with lower risk-adjusted odds of in-hospital mortality (odds ratio [OR], 0.59 [95% CI, 0.38 to 0.93], P = .02; predicted risk difference [RD], 4.0% [95% CI, 0.9% to 7.0%]). However, completion of each individual bundle element within 1 hour was not significantly associated with lower risk-adjusted mortality (blood culture: OR, 0.73 [95% CI, 0.51 to 1.06], P = .10; RD, 2.6% [95% CI, -0.5% to 5.7%]; antibiotics: OR, 0.78 [95% CI, 0.55 to 1.12], P = .18; RD, 2.1% [95% CI, -1.1% to 5.2%], and fluid bolus: OR, 0.88 [95% CI, 0.56 to 1.37], P = .56; RD, 1.1% [95% CI, -2.6% to 4.8%]).

CONCLUSIONS AND RELEVANCE

In New York State following a mandate for sepsis care, completion of a sepsis bundle within 1 hour compared with not completing the 1-hour sepsis bundle within 1 hour was associated with lower risk-adjusted in-hospital mortality among patients with pediatric sepsis and septic shock.

Correction to: Comparative Analysis of the Expression of Glutathione Peroxidase and Glutathione Reductase Genes in Human Sperm after Cryopreservation

The second author's name should read A. Yu. Romanov.

Blood eosinophil count thresholds and exacerbations in patients with chronic obstructive pulmonary disease

BACKGROUND

Eosinophilic airway inflammation in patients with chronic obstructive pulmonary disease (COPD) is associated with exacerbations and responsivity to steroids, suggesting potential shared mechanisms with eosinophilic asthma. However, there is no consistent blood eosinophil count that has been used to define the increased exacerbation risk.

OBJECTIVE

We sought to investigate blood eosinophil counts associated with exacerbation risk in patients with COPD.

METHODS

Blood eosinophil counts and exacerbation risk were analyzed in patients with moderate-to-severe COPD by using 2 independent studies of former and current smokers with longitudinal data. The Genetic Epidemiology of COPD (COPDGene) study was analyzed for discovery (n = 1,553), and the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) study was analyzed for validation (n = 1,895). A subset of the ECLIPSE study subjects were used to assess the stability of blood eosinophil counts over time.

RESULTS

COPD exacerbation risk increased with higher eosinophil counts. An eosinophil count threshold of 300 cells/μL or greater showed adjusted incidence rate ratios for exacerbations of 1.32 in the COPDGene study (95% CI, 1.10-1.63). The cutoff of 300 cells/μL or greater was validated for prospective risk of exacerbation in the ECLIPSE study, with adjusted incidence rate ratios of 1.22 (95% CI, 1.06-1.41) using 3-year follow-up data. Stratified analysis confirmed that the increased exacerbation risk associated with an eosinophil count of 300 cells/μL or greater was driven by subjects with a history of frequent exacerbations in both the COPDGene and ECLIPSE studies.

CONCLUSIONS

Patients with moderate-to-severe COPD and blood eosinophil counts of 300 cells/μL or greater had an increased risk exacerbations in the COPDGene study, which was prospectively validated in the ECLIPSE study.

Genotype imputation performance of three reference panels using African ancestry individuals

Genotype imputation estimates unobserved genotypes from genome-wide makers, to increase genome coverage and power for genome-wide association studies. Imputation has been successful for European ancestry populations in which very large reference panels are available. Smaller subsets of African descent populations are available in 1000 Genomes (1000G), the Consortium on Asthma among African ancestry Populations in the Americas (CAAPA) and the Haplotype Reference Consortium (HRC). We compared the performance of these reference panels when imputing variation in 3747 African Americans (AA) from two cohorts (HCV and COPDGene) genotyped using Illumina Omni microarrays. The haplotypes of 2504 (1000G), 883 (CAAPA) and 32,470 individuals (HRC) were used as reference. We compared the number of variants, imputation quality, imputation accuracy and coverage between panels. In both cohorts, 1000G imputed 1.5-1.6× more variants than CAAPA and 1.2× more than HRC. Similar findings were observed for variants with imputation R2 > 0.5 and for rare, low-frequency, and common variants. When merging imputed variants of the three panels, the total number was 62-63 M with 20 M overlapping variants imputed by all three panels, and a range of 5-15 M variants imputed exclusively with one of them. For overlapping variants, imputation quality was highest for HRC, followed by 1000G, then CAAPA, and improved as the minor allele frequency increased. 1000G, HRC and CAAPA provided high performance and accuracy for imputation of African American individuals, increasing the number of variants available for subsequent analyses. These panels are complementary and would benefit from the development of an integrated African reference panel.

RNA sequencing identifies novel non-coding RNA and exon-specific effects associated with cigarette smoking

BACKGROUND

Cigarette smoking is the leading modifiable risk factor for disease and death worldwide. Previous studies quantifying gene-level expression have documented the effect of smoking on mRNA levels. Using RNA sequencing, it is possible to analyze the impact of smoking on complex regulatory phenomena (e.g. alternative splicing, differential isoform usage) leading to a more detailed understanding of the biology underlying smoking-related disease.

METHODS

We used whole-blood RNA sequencing to describe gene and exon-level expression differences between 229 current and 286 former smokers in the COPDGene study. We performed differential gene expression and differential exon usage analyses using the voom/limma and DEXseq R packages. Samples from current and former smokers were compared while controlling for age, gender, race, lifetime smoke exposure, cell counts, and technical covariates.

RESULTS

At an adjusted p-value <0.05, 171 genes were differentially expressed between current and former smokers. Differentially expressed genes included 7 long non-coding RNAs that have not been previously associated with smoking: LINC00599, LINC01362, LINC00824, LINC01624, RP11-563D10.1, RP11-98G13.1, AC004791.2. Secondary analysis of acute smoking (having smoked within 2-h) revealed 5 of the 171 smoking genes demonstrated an acute response above the baseline effect of chronic smoking. Exon-level analyses identified 9 exons from 8 genes with significant differential usage by smoking status, suggesting smoking-induced changes in isoform expression.

CONCLUSIONS

Transcriptomic changes at the gene and exon levels from whole blood can refine our understanding of the molecular mechanisms underlying the response to smoking.

Analysis of sequence data to identify potential risk variants for oral clefts in multiplex families

BACKGROUND

Nonsyndromic oral clefts are craniofacial malformations, which include cleft lip with or without cleft palate. The etiology for oral clefts is complex with both genetic and environmental factors contributing to risk. Previous genome-wide association (GWAS) studies have identified multiple loci with small effects; however, many causal variants remain elusive.

METHODS

In this study, we address this by specifically looking for rare, potentially damaging variants in family-based data. We analyzed both whole exome sequence (WES) data and whole genome sequence (WGS) data in multiplex cleft families to identify variants shared by affected individuals.

RESULTS

Here we present the results from these analyses. Our most interesting finding was from a single Syrian family, which showed enrichment of nonsynonymous and potentially damaging rare variants in two genes: CASP9 and FAT4.

CONCLUSION

Neither of these candidate genes has previously been associated with oral clefts and, if confirmed as contributing to disease risk, may indicate novel biological pathways in the genetic etiology for oral clefts.

Whole exome association of rare deletions in multiplex oral cleft families

By sequencing the exomes of distantly related individuals in multiplex families, rare mutational and structural changes to coding DNA can be characterized and their relationship to disease risk can be assessed. Recently, several rare single nucleotide variants (SNVs) were associated with an increased risk of nonsyndromic oral cleft, highlighting the importance of rare sequence variants in oral clefts and illustrating the strength of family-based study designs. However, the extent to which rare deletions in coding regions of the genome occur and contribute to risk of nonsyndromic clefts is not well understood. To identify putative structural variants underlying risk, we developed a pipeline for rare hemizygous deletions in families from whole exome sequencing and statistical inference based on rare variant sharing. Among 56 multiplex families with 115 individuals, we identified 53 regions with one or more rare hemizygous deletions. We found 45 of the 53 regions contained rare deletions occurring in only one family member. Members of the same family shared a rare deletion in only eight regions. We also devised a scalable global test for enrichment of shared rare deletions.

Hemizygous Deletion on Chromosome 3p26.1 Is Associated with Heavy Smoking among African American Subjects in the COPDGene Study

Many well-powered genome-wide association studies have identified genetic determinants of self-reported smoking behaviors and measures of nicotine dependence, but most have not considered the role of structural variants, such as copy number variation (CNVs), influencing these phenotypes. Here, we included 2,889 African American and 6,187 non-Hispanic White subjects from the COPDGene cohort (http://www.copdgene.org) to carefully investigate the role of polymorphic CNVs across the genome on various measures of smoking behavior. We identified a CNV component (a hemizygous deletion) on chromosome 3p26.1 associated with two quantitative phenotypes related to smoking behavior among African Americans. This polymorphic hemizygous deletion is significantly associated with pack-years and cigarettes smoked per day among African American subjects in the COPDGene study. We sought evidence of replication in African Americans from the population based Atherosclerosis Risk in Communities (ARIC) study. While we observed similar CNV counts, the extent of exposure to cigarette smoking among ARIC subjects was quite different and the smaller sample size of heavy smokers in ARIC severely limited statistical power, so we were unable to replicate our findings from the COPDGene cohort. But meta-analyses of COPDGene and ARIC study subjects strengthened our association signal. However, a few linkage studies have reported suggestive linkage to the 3p26.1 region, and a few genome-wide association studies (GWAS) have reported markers in the gene (GRM7) nearest to this 3p26.1 area of polymorphic deletions are associated with measures of nicotine dependence among subjects of European ancestry.

Exome Array Analysis Identifies a Common Variant in IL27 Associated with Chronic Obstructive Pulmonary Disease

RATIONALE

Chronic obstructive pulmonary disease (COPD) susceptibility is in part related to genetic variants. Most genetic studies have been focused on genome-wide common variants without a specific focus on coding variants, but common and rare coding variants may also affect COPD susceptibility.

OBJECTIVES

To identify coding variants associated with COPD.

METHODS

We tested nonsynonymous, splice, and stop variants derived from the Illumina HumanExome array for association with COPD in five study populations enriched for COPD. We evaluated single variants with a minor allele frequency greater than 0.5% using logistic regression. Results were combined using a fixed effects meta-analysis. We replicated novel single-variant associations in three additional COPD cohorts.

MEASUREMENTS AND MAIN RESULTS

We included 6,004 control subjects and 6,161 COPD cases across five cohorts for analysis. Our top result was rs16969968 (P = 1.7 × 10(-14)) in CHRNA5, a locus previously associated with COPD susceptibility and nicotine dependence. Additional top results were found in AGER, MMP3, and SERPINA1. A nonsynonymous variant, rs181206, in IL27 (P = 4.7 × 10(-6)) was just below the level of exome-wide significance but attained exome-wide significance (P = 5.7 × 10(-8)) when combined with results from other cohorts. Gene expression datasets revealed an association of rs181206 and the surrounding locus with expression of multiple genes; several were differentially expressed in COPD lung tissue, including TUFM.

CONCLUSIONS

In an exome array analysis of COPD, we identified nonsynonymous variants at previously described loci and a novel exome-wide significant variant in IL27. This variant is at a locus previously described in genome-wide associations with diabetes, inflammatory bowel disease, and obesity and appears to affect genes potentially related to COPD pathogenesis.

Family Presence During Pediatric Tracheal Intubations

IMPORTANCE

Family-centered care, which supports family presence (FP) during procedures, is now a widely accepted standard at health care facilities that care for children. However, there is a paucity of data regarding the practice of FP during tracheal intubation (TI) in pediatric intensive care units (PICUs). Family presence during procedures in PICUs has been advocated.

OBJECTIVE

To describe the current practice of FP during TI and evaluate the association with procedural and clinician (including physician, respiratory therapist, and nurse practitioner) outcomes across multiple PICUs.

DESIGN, SETTING, AND PARTICIPANTS

Prospective cohort study in which all TIs from July 2010 to March 2014 in the multicenter TI database (National Emergency Airway Registry for Children [NEAR4KIDS]) were analyzed. Family presence was defined as a family member present during TI. This study included all TIs in patients younger than 18 years in 22 international PICUs.

EXPOSURES

Family presence and no FP during TI in the PICU.

MAIN OUTCOMES AND MEASURES

The percentage of FP during TIs. First attempt success rate, adverse TI-associated events, multiple attempts (≥ 3), oxygen desaturation (oxygen saturation as measured by pulse oximetry <80%), and self-reported team stress level.

RESULTS

A total of 4969 TI encounters were reported. Among those, 81% (n = 4030) of TIs had documented FP status (with/without). The median age of participants with FP was 2 years and 1 year for those without FP. The average percentage of TIs with FP was 19% and varied widely across sites (0%-43%; P < .001). Tracheal intubations with FP (vs without FP) were associated with older patients (median, 2 years vs 1 year; P = .04), lower Paediatric Index of Mortality 2 score, and pediatric resident as the first airway clinician (23%, n = 179 vs 18%, n = 584; odds ratio [OR], 1.4; 95% CI, 1.2-1.7). Tracheal intubations with FP and without FP were no different in the first attempt success rate (OR, 1.00; 95% CI, 0.85-1.18), adverse TI-associated events (any events: OR, 1.06; 95% CI, 0.85-1.30 and severe events: OR, 1.04; 95% CI, 0.75-1.43), multiple attempts (≥ 3) (OR, 1.03; 95% CI, 0.82-1.28), oxygen desaturation (oxygen saturation <80%) (OR, 0.97; 95% CI, 0.80-1.18), or self-reported team stress level (OR, 1.09; 95% CI, 0.92-1.31). This result persisted after adjusting for patient and clinician confounders.

CONCLUSIONS AND RELEVANCE

Wide variability exists in FP during TIs across PICUs. Family presence was not associated with first attempt success, adverse TI-associated events, oxygen desaturation (<80%), or higher team stress level. Our data suggest that FP during TI can safely be implemented as part of a family-centered care model in the PICU.

Identifying a Deletion Affecting Total Lung Capacity Among Subjects in the COPDGene Study Cohort

Chronic obstructive pulmonary disease (COPD) is a progressive disease with both environmental and genetic risk factors. Genome-wide association studies (GWAS) have identified multiple genomic regions influencing risk of COPD. To thoroughly investigate the genetic etiology of COPD, however, it is also important to explore the role of copy number variants (CNVs) because the presence of structural variants can alter gene expression and can be causal for some diseases. Here, we investigated effects of polymorphic CNVs on quantitative measures of pulmonary function and chest computed tomography (CT) phenotypes among subjects enrolled in COPDGene, a multisite study. COPDGene subjects consist of roughly one-third African American (AA) and two-thirds non-Hispanic white adult smokers (with or without COPD). We estimated CNVs using PennCNV on 9,076 COPDGene subjects using Illumina's Omni-Express genome-wide marker array. We tested for association between polymorphic CNV components (defined as disjoint intervals of copy number regions) for several quantitative phenotypes associated with COPD within each racial group. Among the AAs, we identified a polymorphic CNV on chromosome 5q35.2 located between two genes (FAM153B and SIMK1, but also harboring several pseudo-genes) giving genome-wide significance in tests of association with total lung capacity (TLCCT ) as measured by chest CT scans. This is the first study of genome-wide association tests of polymorphic CNVs and TLCCT . Although the ARIC cohort did not have the phenotype of TLCCT , we found similar counts of CNV deletions and amplifications among AA and European subjects in this second cohort.

A genome-wide study of inherited deletions identified two regions associated with nonsyndromic isolated oral clefts

BACKGROUND

DNA copy number variants play an important part in the development of common birth defects such as oral clefts. Individual patients with multiple birth defects (including oral clefts) have been shown to carry small and large chromosomal deletions.

METHODS

We investigated the role of polymorphic copy number deletions by comparing transmission rates of deletions from parents to offspring in case-parent trios of European ancestry ascertained through a cleft proband with trios ascertained through a normal offspring. DNA copy numbers in trios were called using the joint hidden Markov model in the freely available PennCNV software. All statistical analyses were performed using Bioconductor tools in the open source environment R.

RESULTS

We identified a 67 kb region in the gene MGAM on chromosome 7q34, and a 206 kb region overlapping genes ADAM3A and ADAM5 on chromosome 8p11, where deletions are more frequently transmitted to cleft offspring than control offspring.

CONCLUSIONS

These genes or nearby regulatory elements may be involved in the etiology of oral clefts.

Identification of functional variants for cleft lip with or without cleft palate in or near PAX7, FGFR2, and NOG by targeted sequencing of GWAS loci

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.

Admixture mapping identifies a quantitative trait locus associated with FEV1/FVC in the COPDGene Study

African Americans are admixed with genetic contributions from European and African ancestral populations. Admixture mapping leverages this information to map genes influencing differential disease risk across populations. We performed admixture and association mapping in 3,300 African American current or former smokers from the COPDGene Study. We analyzed estimated local ancestry and SNP genotype information to identify regions associated with FEV1 /FVC, the ratio of forced expiratory volume in one second to forced vital capacity, measured by spirometry performed after bronchodilator administration. Global African ancestry inversely associated with FEV1 /FVC (P = 0.035). Genome-wide admixture analysis, controlling for age, gender, body mass index, current smoking status, pack-years smoked, and four principal components summarizing the genetic background of African Americans in the COPDGene Study, identified a region on chromosome 12q14.1 associated with FEV1 /FVC (P = 2.1 × 10(-6) ) when regressed on local ancestry. Allelic association in this region of chromosome 12 identified an intronic variant in FAM19A2 (rs348644) as associated with FEV1 /FVC (P = 1.76 × 10(-6) ). By combining admixture and association mapping, a marker on chromosome 12q14.1 was identified as being associated with reduced FEV1 /FVC ratio among African Americans in the COPDGene Study.