Damaging Mutations in Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.

We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.

We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.

This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

Detecting Gene-Gene Interaction among DNA Repair Genes in Chinese non-Syndromic Cleft lip with or Without Palate Trios

OBJECTIVE

The objective of this study is to investigate the gene-gene interactions associated with NSCL/P among DNA repair genes.

DESIGN

This study included 806 NSCL/P case-parent trios from China. Quality control process was conducted for genotyped single nucleotide polymorphisms (SNPs) located in six DNA repair genes (ATR, ERCC4, RFC1, TYMS, XRCC1 and XRCC3). We tested gene-gene interactions with Cordell's method using statistical package TRIO in R software. Bonferroni corrected significance level was set as P = 4.24 × 10-4. We also test the robustness of the interactions by permutation tests.

SETTING

Not applicable.

PATIENTS/PARTICIPANTS

A total of 806 NSCL/P case-parent trios (complete trios: 682, incomplete trios: 124) with Chinese ancestry.

INTERVENTIONS

Not applicable.

MAIN OUTCOME MEASURE(S)

Not applicable.

RESULTS

A total of 118 SNPs were extracted for the interaction tests. Fourteen pairs of significant interactions were identified after Bonferroni correction, which were confirmed in permutation tests. Twelve pairs were between ATR and ERCC4 or XRCC3. The most significant interaction occurred between rs2244500 in TYMS and rs3213403 in XRCC1(P = 8.16 × 10-15).

CONCLUSIONS

The current study identified gene-gene interactions among DNA repair genes in 806 Chinese NSCL/P trios, providing additional evidence for the complicated genetic structure underlying NSCL/P. ATR, ERCC4, XRCC3, TYMS and RFC1 were suggested to be possible candidate genes for NSCL/P.

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

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

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate

Cleft palate (CP) is one of the most common craniofacial birth defects; however, there are relatively few established genetic risk factors associated with its occurrence despite high heritability. Historically, CP has been studied as a single phenotype, although it manifests across a spectrum of defects involving the hard and/or soft palate. We performed a genome-wide association study using transmission disequilibrium tests of 435 case-parent trios to evaluate broad risks for any cleft palate (ACP) (n = 435), and subtype-specific risks for any cleft soft palate (CSP), (n = 259) and any cleft hard palate (CHP) (n = 125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p = 4.24 × 10-8) associated with CHP. One gene at this locus, angiopoietin-like 2 (ANGPTL2), plays a role in osteoblast differentiation. It is expressed both in craniofacial tissue of human embryos and developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p < 5 × 10-6), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios for the 20 loci were most similar across all 3 groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either subtype. We also found nominal evidence of replication (p < 0.05) for 22 SNPs previously associated with orofacial clefts. Our study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting

PURPOSE

Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls.

METHODS

We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria.

RESULTS

9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance.

CONCLUSION

These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

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

BACKGROUND

COPD diagnosis is tightly linked to the fixed-ratio spirometry criteria of FEV1/FVC < 0.7. African-Americans are less often diagnosed with COPD.

OBJECTIVE

Compare COPD diagnosis by fixed-ratio with findings and outcomes by race.

DESIGN

Genetic Epidemiology of COPD (COPDGene) (2007-present), cross-sectional comparing non-Hispanic white (NHW) and African-American (AA) participants for COPD diagnosis, manifestations, and outcomes.

SETTING

Multicenter, longitudinal US cohort study.

PARTICIPANTS

Current or former smokers with ≥ 10-pack-year smoking history enrolled at 21 clinical centers including over-sampling of participants with known COPD and AA. Exclusions were pre-existing non-COPD lung disease, except for a history of asthma.

MEASUREMENTS

Subject diagnosis by conventional criteria. Mortality, imaging, respiratory symptoms, function, and socioeconomic characteristics, including area deprivation index (ADI). Matched analysis (age, sex, and smoking status) of AA vs. NHW within participants without diagnosed COPD (GOLD 0; FEV1 ≥ 80% predicted and FEV1/FVC ≥ 0.7).

RESULTS

Using the fixed ratio, 70% of AA (n = 3366) were classified as non-COPD, versus 49% of NHW (n = 6766). AA smokers were younger (55 vs. 62 years), more often current smoking (80% vs. 39%), with fewer pack-years but similar 12-year mortality. Density distribution plots for FEV1 and FVC raw spirometry values showed disproportionate reductions in FVC relative to FEV1 in AA that systematically led to higher ratios. The matched analysis demonstrated GOLD 0 AA had greater symptoms, worse DLCO, spirometry, BODE scores (1.03 vs 0.54, p < 0.0001), and greater deprivation than NHW.

LIMITATIONS

Lack of an alternative diagnostic metric for comparison.

CONCLUSIONS

The fixed-ratio spirometric criteria for COPD underdiagnosed potential COPD in AA participants when compared to broader diagnostic criteria. Disproportionate reductions in FVC relative to FEV1 leading to higher FEV1/FVC were identified in these participants and associated with deprivation. Broader diagnostic criteria for COPD are needed to identify the disease across all populations.

Causes of and Clinical Features Associated with Death in Tobacco Cigarette Users by Lung Function Impairment

Rationale: Cigarette smoking contributes to the risk of death through different mechanisms. Objectives: To determine how causes of and clinical features associated with death vary in tobacco cigarette users by lung function impairment. Methods: We stratified current and former tobacco cigarette users enrolled in Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) into normal spirometry, PRISm (Preserved Ratio Impaired Spirometry), Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1-2 COPD, and GOLD 3-4 COPD. Deaths were identified via longitudinal follow-up and Social Security Death Index search. Causes of death were adjudicated after a review of death certificates, medical records, and next-of-kin interviews. We tested associations between baseline clinical variables and all-cause mortality using multivariable Cox proportional hazards models. Measurements and Main Results: Over a 10.1-year median follow-up, 2,200 deaths occurred among 10,132 participants (age 59.5 ± 9.0 yr; 46.6% women). Death from cardiovascular disease was most frequent in PRISm (31% of deaths). Lung cancer deaths were most frequent in GOLD 1-2 (18% of deaths vs. 9-11% in other groups). Respiratory deaths outpaced competing causes of death in GOLD 3-4, particularly when BODE index ⩾7. St. George's Respiratory Questionnaire score ⩾25 was associated with higher mortality in all groups: Hazard ratio (HR), 1.48 (1.20-1.84) normal spirometry; HR, 1.40 (1.05-1.87) PRISm; HR, 1.80 (1.49-2.17) GOLD 1-2; HR, 1.65 (1.26-2.17) GOLD 3-4. History of respiratory exacerbations was associated with higher mortality in GOLD 1-2 and GOLD 3-4, quantitative emphysema in GOLD 1-2, and airway wall thickness in PRISm and GOLD 3-4. Conclusions: Leading causes of death vary by lung function impairment in tobacco cigarette users. Worse respiratory-related quality of life is associated with all-cause mortality regardless of lung function.

Xanthine oxidoreductase gene polymorphisms are associated with high risk of sepsis and organ failure

BACKGROUND

Sepsis and associated organ failures confer substantial morbidity and mortality. Xanthine oxidoreductase (XOR) is implicated in the development of tissue oxidative damage in a wide variety of respiratory and cardiovascular disorders including sepsis and sepsis-associated acute respiratory distress syndrome (ARDS). We examined whether single nucleotide polymorphisms (SNPs) in the XDH gene (encoding XOR) might influence susceptibility to and outcome in patients with sepsis.

METHODS

We genotyped 28 tag SNPs in XDH gene in the CELEG cohort, including 621 European American (EA) and 353 African American (AA) sepsis patients. Serum XOR activity was measured in a subset of CELEG subjects. Additionally, we assessed the functional effects of XDH variants utilizing empirical data from different integrated software tools and datasets.

RESULTS

Among AA patients, six intronic variants (rs206805, rs513311, rs185925, rs561525, rs2163059, rs13387204), in a region enriched with regulatory elements, were associated with risk of sepsis (P < 0.008-0.049). Two out of six SNPs (rs561525 and rs2163059) were associated with risk of sepsis-associated ARDS in an independent validation cohort (GEN-SEP) of 590 sepsis patients of European descent. Two common SNPs (rs1884725 and rs4952085) in tight linkage disequilibrium (LD) provided strong evidence for association with increased levels of serum creatinine (P_adjusted<0.0005 and 0.0006, respectively), suggesting a role in increased risk of renal dysfunction. In contrast, among EA ARDS patients, the missense variant rs17011368 (I703V) was associated with enhanced mortality at 60-days (P < 0.038). We found higher serum XOR activity in 143 sepsis patients (54.5 ± 57.1 mU/mL) compared to 31 controls (20.9 ± 12.4 mU/mL, P = 1.96 × 10^- 13). XOR activity was associated with the lead variant rs185925 among AA sepsis patients with ARDS (P < 0.005 and P_adjusted<0.01). Multifaceted functions of prioritized XDH variants, as suggested by various functional annotation tools, support their potential causality in sepsis.

CONCLUSIONS

Our findings suggest that XOR is a novel combined genetic and biochemical marker for risk and outcome in patients with sepsis and ARDS.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate

Orofacial clefts (OFCs) are the most common craniofacial birth defects and are often categorized into two etiologically distinct groups: cleft lip with or without cleft palate (CL/P) and isolated cleft palate (CP). CP is highly heritable, but there are still relatively few established genetic risk factors associated with its occurrence compared to CL/P. Historically, CP has been studied as a single phenotype despite manifesting across a spectrum of defects involving the hard and/or soft palate. We performed GWAS using transmission disequilibrium tests using 435 case-parent trios to evaluate broad risks for any cleft palate (ACP, n=435), as well as subtype-specific risks for any cleft soft palate (CSP, n=259) and any cleft hard palate (CHP, n=125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p=4.24×10 -8 ) associated with CHP. One gene at this locus, angiopoietin-like 2 ( ANGPTL2 ), plays a role in osteoblast differentiation. It is expressed in craniofacial tissue of human embryos, as well as in the developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10 -6 ), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios (ORs) for each of the 20 loci were most similar across all three groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either the CSP or CHP groups. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with cleft palate (including CL/P). Interestingly, most SNPs associated with CL/P cases were found to convey the opposite effect in those replicated in our dataset for CP only. Ours is the first study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

Evidence of the folate-mediated one-carbon metabolism pathway genes in controlling the non-syndromic oral clefts risks

The folate-mediated one-carbon metabolism pathway is thought to play an important role in the etiology of non-syndromic oral clefts (NSOFC), although none of the genes in this pathway has shown significant signals in genome-wide association studies (GWAS). Recent evidence indicated that enhanced understanding could be gained by aggregating multiple SNPs effect simultaneously into polygenic risk score (PRS) to assess its association with disease risks. This study is aimed to assess the association between the genetic effect of folate-mediated one-carbon metabolism pathway and NSOFC risks using PRS based on a case-parent trio design. A total of 297 SNPs mapped from 18 genes in the folate-mediated one-carbon metabolism pathway were aggregated from a GWAS of 2458 case-parent trios recruited from an international consortium. We found a PRS based on the folate-mediated one-carbon metabolism pathway was significant among all NSOFC trios (OR = 1.95, 95% CI: 1.66-2.28, p = 2.39 × 10^-16 ), as well as two major subtypes, non-syndromic cleft lip with or without cleft palate (NSCL/P) trios (OR = 1.71, 95% CI: 1.50-1.96, p = 7.66 × 10^-15 ) and non-syndromic cleft palate only (NSCPO) trios (OR = 1.51, 95% CI: 1.36-1.68, p = 2.1 × 10^-14 ). Similar results were also observed in further subgroup analyses stratified into Asian and European trios. The averaged PRS of the folate-mediated one-carbon metabolism pathway varied between the NSOFC case group and its comparison group (p < 0.05) with higher average PRS in the cases. Moreover, the top 5% pathway PRS group had 2.25 (95% CI: 1.85-2.73) times increased NSOFC risk, also 3.09 (95% CI: 2.50-3.81) and 2.06 (95% CI: 1.39-3.02) times increased risk of NSCL/P and NSCPO compared to the remainder of the distribution. The results of our study confirmed the folate-mediated one-carbon metabolism pathway was important in controlling risk to NSOFC and this study enhanced evidence towards understanding the genetic risks of NSOFC.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10-5) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

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

BACKGROUND

The association between genetic variants on the X chromosome to risk of COPD has not been fully explored. We hypothesize that the X chromosome harbors variants important in determining risk of COPD related phenotypes and may drive sex differences in COPD manifestations.

METHODS

Using X chromosome data from three COPD-enriched cohorts of adult smokers, we performed X chromosome specific quality control, imputation, and testing for association with COPD case-control status, lung function, and quantitative emphysema. Analyses were performed among all subjects, then stratified by sex, and subsequently combined in meta-analyses.

RESULTS

Among 10,193 subjects of non-Hispanic white or European ancestry, a variant near TMSB4X, rs5979771, reached genome-wide significance for association with lung function measured by FEV₁/FVC ([Formula: see text] 0.020, SE 0.004, p 4.97 × 10^-08), with suggestive evidence of association with FEV₁ ([Formula: see text] 0.092, SE 0.018, p 3.40 × 10^-07). Sex-stratified analyses revealed X chromosome variants that were differentially trending in one sex, with significantly different effect sizes or directions.

CONCLUSIONS

This investigation identified loci influencing lung function, COPD, and emphysema in a comprehensive genetic association meta-analysis of X chromosome genetic markers from multiple COPD-related datasets. Sex differences play an important role in the pathobiology of complex lung disease, including X chromosome variants that demonstrate differential effects by sex and variants that may be relevant through escape from X chromosome inactivation. Comprehensive interrogation of the X chromosome to better understand genetic control of COPD and lung function is important to further understanding of disease pathology. Trial registration Genetic Epidemiology of COPD Study (COPDGene) is registered at ClinicalTrials.gov, NCT00608764 (Active since January 28, 2008). Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints Study (ECLIPSE), GlaxoSmithKline study code SCO104960, is registered at ClinicalTrials.gov, NCT00292552 (Active since February 16, 2006). Genetics of COPD in Norway Study (GenKOLS) holds GlaxoSmithKline study code RES11080, Genetics of Chronic Obstructive Lung Disease.

Benchmarking statistical methods for analyzing parent-child dyads in genetic association studies

Genetic association studies of child health outcomes often employ family-based study designs. One of the most popular family-based designs is the case-parent trio design that considers the smallest possible nuclear family consisting of two parents and their affected child. This trio design is particularly advantageous for studying relatively rare disorders because it is less prone to type 1 error inflation due to population stratification compared to population-based study designs (e.g., case-control studies). However, obtaining genetic data from both parents is difficult, from a practical perspective, and many large studies predominantly measure genetic variants in mother-child dyads. While some statistical methods for analyzing parent-child dyad data (most commonly involving mother-child pairs) exist, it is not clear if they provide the same advantage as trio methods in protecting against population stratification, or if a specific dyad design (e.g., case-mother dyads vs. case-mother/control-mother dyads) is more advantageous. In this article, we review existing statistical methods for analyzing genome-wide marker data on dyads and perform extensive simulation experiments to benchmark their type I errors and statistical power under different scenarios. We extend our evaluation to existing methods for analyzing a combination of case-parent trios and dyads together. We apply these methods on genotyped and imputed data from multiethnic mother-child pairs only, case-parent trios only or combinations of both dyads and trios from the Gene, Environment Association Studies consortium (GENEVA), where each family was ascertained through a child affected by nonsyndromic cleft lip with or without cleft palate. Results from the GENEVA study corroborate the findings from our simulation experiments. Finally, we provide recommendations for using statistical genetic association methods for dyads.

A genome-wide association study of bronchodilator response in participants of European and African ancestry from six independent cohorts

Introduction

Bronchodilator response (BDR) is a measurement of acute bronchodilation in response to short-acting β2-agonists, with a heritability between 10 and 40%. Identifying genetic variants associated with BDR may lead to a better understanding of its complex pathophysiology.

Methods

We performed a genome-wide association study (GWAS) of BDR in six adult cohorts with participants of European ancestry (EA) and African ancestry (AA) including community cohorts and cohorts ascertained on the basis of obstructive pulmonary disease. Validation analysis was carried out in two paediatric asthma cohorts.

Results

A total of 10 623 EA and 3597 AA participants were included in the analyses. No single nucleotide polymorphism (SNP) was associated with BDR at the conventional genome-wide significance threshold (p<5×10-8). Performing fine mapping and using a threshold of p<5×10-6 to identify suggestive variants of interest, we identified three SNPs with possible biological relevance: rs35870000 (within FREM1), which may be involved in IgE- and IL5-induced changes in airway smooth muscle cell responsiveness; rs10426116 (within ZNF284), a zinc finger protein, which has been implicated in asthma and BDR previously; and rs4782614 (near ATP2C2), involved in calcium transmembrane transport. Validation in paediatric cohorts yielded no significant SNPs, possibly due to age-genotype interaction effects.

Conclusion

Ancestry-stratified and ancestry-combined GWAS meta-analyses of over 14 000 participants did not identify genetic variants associated with BDR at the genome-wide significance threshold, although a less stringent threshold identified three variants showing suggestive evidence of association. A common definition and protocol for measuring BDR in research may improve future efforts to identify variants associated with BDR.

Individual and Combined Association Between Prenatal Polysubstance Exposure and Childhood Risk of Attention-Deficit/Hyperactivity Disorder

IMPORTANCE

Polysubstance use among pregnant women has increased because of the opioid epidemic and the increasing legalization of cannabis along with persistent tobacco and alcohol consumption. Previous research on prenatal substance use and the child's risk of attention-deficit/hyperactivity disorder (ADHD) has mostly focused on single-substance exposures; simultaneous examination of multiple substance use and assessment of their synergistic health consequences is needed.

OBJECTIVES

To assess the consequences of the use of specific substances during pregnancy, investigate whether the interaction of multiple prenatal substance exposures is associated with increases in the risk of childhood ADHD, and estimate the aggregate burden of polysubstance exposure during gestation.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study analyzed data from the Boston Birth Cohort from 1998 to 2019. The sample of the present study comprised a multiethnic urban cohort of mother-child pairs who were predominantly low income. A total of 3138 children who were enrolled shortly after birth at Boston Medical Center were included and followed up from age 6 months to 21 years.

EXPOSURES

Substance use during pregnancy was identified based on self-reported tobacco smoking, alcohol consumption, and use of cannabis, cocaine, or opioids in any trimester of pregnancy. Diagnostic codes for neonatal opioid withdrawal syndrome or neonatal abstinence syndrome from the International Classification of Diseases, Ninth Revision, and the International Classification of Diseases, Tenth Revision, were also used to identify opioid exposure during gestation.

MAIN OUTCOMES AND MEASURES

ADHD diagnosis in the child's electronic medical record.

RESULTS

Among 3138 children (1583 boys [50.4%]; median age, 12 years [IQR, 9-14 years]; median follow-up, 10 years [IQR, 7-12 years]) in the final analytic sample, 486 (15.5%) had an ADHD diagnosis and 2652 (84.5%) were neurotypical. The median postnatal follow-up duration was 12 years (IQR, 9-14 years). Among mothers, 46 women (1.5%) self-identified as Asian (non-Pacific Islander), 701 (22.3%) as Hispanic, 1838 (58.6%) as non-Hispanic Black, 227 (7.2%) as non-Hispanic White, and 326 (10.4%) as other races and/or ethnicities (including American Indian or Indigenous, Cape Verdean, Pacific Islander, multiracial, other, or unknown). A total of 759 women (24.2%) reported the use of at least 1 substance during pregnancy, with tobacco being the most frequently reported (580 women [18.5%]). Cox proportional hazards models revealed that opioid exposure (60 children) had the highest adjusted hazard ratio (HR) for ADHD (2.19; 95% CI, 1.10-4.37). After including main statistical effects of all individual substances in an elastic net regression model, the HR of opioids was reduced to 1.60, and evidence of a statistical interaction between opioids and both cannabis and alcohol was found, producing 1.42 and 1.15 times higher risk of ADHD, respectively. The interaction between opioids and smoking was also associated with a higher risk of ADHD (HR, 1.17).

CONCLUSIONS AND RELEVANCE

The findings of this study suggest that it is important to consider prenatal concurrent exposure to multiple substances and their possible interactions when counseling women regarding substance use during pregnancy, the future risk of ADHD for their children, and strategies for cessation and treatment programs.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

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

IMPORTANCE

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

OBJECTIVE

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

DESIGN, SETTING, AND PARTICIPANTS

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

EXPOSURES

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

MAIN OUTCOMES AND MEASURES

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

RESULTS

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

CONCLUSIONS AND RELEVANCE

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

Accurate assignment of disease liability to genetic variants using only population data

PURPOSE

The growing size of public variant repositories prompted us to test the accuracy of pathogenicity prediction of DNA variants using population data alone.

METHODS

Under the a priori assumption that the ratio of the prevalence of variants in healthy population vs that in affected populations form 2 distinct distributions (pathogenic and benign), we used a Bayesian method to assign probability to a variant belonging to either distribution.

RESULTS

The approach, termed Bayesian prevalence ratio (BayPR), accurately parsed 300 of 313 expertly curated CFTR variants: 284 of 296 pathogenic/likely pathogenic variants in 1 distribution and 16 of 17 benign/likely benign variants in another. BayPR produced an area under the receiver operating characteristic curve of 0.99 for 103 functionally confirmed missense CFTR variants, which is equal to or exceeds 10 commonly used algorithms (area under the receiver operating characteristic curve range = 0.54-0.99). Application of BayPR to expertly curated variants in 8 genes associated with 7 Mendelian conditions led to the assignment of a disease-causing probability of ≥80% to 1350 of 1374 (98.3%) pathogenic/likely pathogenic variants and of ≤20% to 22 of 23 (95.7%) benign/likely benign variants.

CONCLUSION

Irrespective of the variant type or functional effect, the BayPR approach provides probabilities of pathogenicity for DNA variants responsible for Mendelian disorders using only the variant counts in affected and unaffected population samples.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Pleiotropy method reveals genetic overlap between orofacial clefts at multiple novel loci from GWAS of multi-ethnic trios

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts- the most common craniofacial birth defects in humans- are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10-8, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10-6, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

FAT4 identified as a potential modifier of orofacial cleft laterality

Orofacial clefts (OFCs) are common (1 in 700 births) congenital malformations that include a cleft lip (CL) and cleft lip and palate (CLP). These OFC subtypes are also heterogeneous themselves, with the CL occurring on the left, right, or both sides of the upper lip. Unilateral CL and CLP have a 2:1 bias towards left-sided clefts, suggesting a nonrandom process. Here, we performed a study of left- and right-sided clefts within the CL and CLP subtypes to better understand the genetic factors controlling cleft laterality. We conducted genome-wide modifier analyses by comparing cases that had right unilateral CL (RCL; N = 130), left unilateral CL (LCL; N = 216), right unilateral CLP (RCLP; N = 416), or left unilateral CLP (LCLP; N = 638), and identified a candidate region on 4q28, 400 kb downstream from FAT4, that approached genome-wide significance for LCL versus RCL (p = 8.4 × 10^-8 ). Consistent with its potential involvement as a genetic modifier of CL, we found that Fat4 exhibits a specific domain of expression in the mesenchyme of the medial nasal processes that form the median upper lip. Overall, these results suggest that the epidemiological similarities in left- to right-sided clefts in CL and CLP are not reflected in the genetic association results.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Detecting Gene-Environment Interaction for Maternal Exposures Using Case-Parent Trios Ascertained Through a Case With Non-Syndromic Orofacial Cleft

Two large studies of case-parent trios ascertained through a proband with a non-syndromic orofacial cleft (OFC, which includes cleft lip and palate, cleft lip alone, or cleft palate alone) were used to test for possible gene-environment (G × E) interaction between genome-wide markers (both observed and imputed) and self-reported maternal exposure to smoking, alcohol consumption, and multivitamin supplementation during pregnancy. The parent studies were as follows: GENEVA, which included 1,939 case-parent trios recruited largely through treatment centers in Europe, the United States, and Asia, and 1,443 case-parent trios from the Pittsburgh Orofacial Cleft Study (POFC) also ascertained through a proband with an OFC including three major racial/ethnic groups (European, Asian, and Latin American). Exposure rates to these environmental risk factors (maternal smoking, alcohol consumption, and multivitamin supplementation) varied across studies and among racial/ethnic groups, creating substantial differences in power to detect G × E interaction, but the trio design should minimize spurious results due to population stratification. The GENEVA and POFC studies were analyzed separately, and a meta-analysis was conducted across both studies to test for G × E interaction using the 2 df test of gene and G × E interaction and the 1 df test for G × E interaction alone. The 2 df test confirmed effects for several recognized risk genes, suggesting modest G × E effects. This analysis did reveal suggestive evidence for G × Vitamin interaction for CASP9 on 1p36 located about 3 Mb from PAX7, a recognized risk gene. Several regions gave suggestive evidence of G × E interaction in the 1 df test. For example, for G × Smoking interaction, the 1 df test suggested markers in MUSK on 9q31.3 from meta-analysis. Markers near SLCO3A1 also showed suggestive evidence in the 1 df test for G × Alcohol interaction, and rs41117 near RETREG1 (a.k.a. FAM134B) also gave suggestive significance in the meta-analysis of the 1 df test for G × Vitamin interaction. While it remains quite difficult to obtain definitive evidence for G × E interaction in genome-wide studies, perhaps due to small effect sizes of individual genes combined with low exposure rates, this analysis of two large case-parent trio studies argues for considering possible G × E interaction in any comprehensive study of complex and heterogeneous disorders such as OFC.

The PAX1 locus at 20p11 is a potential genetic modifier for bilateral cleft lip

Nonsyndromic orofacial clefts (OFCs) are a common birth defect and are phenotypically heterogenous in the structure affected by the cleft - cleft lip (CL) and cleft lip and palate (CLP) - as well as other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity, because the genetic architecture of these subtypes is not well understood. We tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using a mixed-model approach. While no novel loci were found, the genetic architecture of UCL was distinct compared to BCL, with 44.03% of suggestive loci having different effects between the two subtypes. To further understand the subtype-specific genetic risk factors, we performed a genome-wide scan for modifiers and found a significant modifier locus on 20p11 (p=7.53×10-9), 300kb downstream of PAX1, that associated with higher odds of BCL vs. UCL, and replicated in an independent cohort (p=0.0018) with no effect in BCLP (p>0.05). We further found that this locus was associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have different genetic architectures. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of modifiers for OFC subtypes and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs.

Relative contributions of family history and a polygenic risk score on COPD and related outcomes: COPDGene and ECLIPSE studies

INTRODUCTION

Family history is a risk factor for chronic obstructive pulmonary disease (COPD). We previously developed a COPD risk score from genome-wide genetic markers (Polygenic Risk Score, PRS). Whether the PRS and family history provide complementary or redundant information for predicting COPD and related outcomes is unknown.

METHODS

We assessed the predictive capacity of family history and PRS on COPD and COPD-related outcomes in non-Hispanic white (NHW) and African American (AA) subjects from COPDGene and ECLIPSE studies. We also performed interaction and mediation analyses.

RESULTS

In COPDGene, family history and PRS were significantly associated with COPD in a single model (PFamHx <0.0001; PPRS<0.0001). Similar trends were seen in ECLIPSE. The area under the receiver operator characteristic curve for a model containing family history and PRS was significantly higher than a model with PRS (p=0.00035) in NHWs and a model with family history (p<0.0001) alone in NHWs and AAs. Both family history and PRS were significantly associated with measures of quantitative emphysema and airway thickness. There was a weakly positive interaction between family history and the PRS under the additive, but not multiplicative scale in NHWs (relative excess risk due to interaction=0.48, p=0.04). Mediation analyses found that a significant proportion of the effect of family history on COPD was mediated through PRS in NHWs (16.5%, 95% CI 9.4% to 24.3%), but not AAs.

CONCLUSION

Family history and the PRS provide complementary information for predicting COPD and related outcomes. Future studies can address the impact of obtaining both measures in clinical practice.

Co-Morbidity Patterns Identified Using Latent Class Analysis of Medications Predict All-Cause Mortality Independent of Other Known Risk Factors: The COPDGene® Study

PURPOSE

Medication patterns include all medications in an individual's clinical profile. We aimed to identify chronic co-morbidity treatment patterns through medication use among COPDGene participants and determine whether these patterns were associated with mortality, acute exacerbations of chronic obstructive pulmonary disease (AECOPD) and quality of life.

MATERIALS AND METHODS

Participants analyzed here completed Phase 1 (P1) and/or Phase 2 (P2) of COPDGene. Latent class analysis (LCA) was used to identify medication patterns and assign individuals into unobserved LCA classes. Mortality, AECOPD, and the St. George's Respiratory Questionnaire (SGRQ) health status were compared in different LCA classes through survival analysis, logistic regression, and Kruskal-Wallis test, respectively.

RESULTS

LCA identified 8 medication patterns from 32 classes of chronic comorbid medications. A total of 8110 out of 10,127 participants with complete covariate information were included. Survival analysis adjusted for covariates showed, compared to a low medication use class, mortality was highest in participants with hypertension+diabetes+statin+antiplatelet medication group. Participants in hypertension+SSRI+statin medication group had the highest odds of AECOPD and the highest SGRQ score at both P1 and P2.

CONCLUSION

Medication pattern can serve as a good indicator of an individual's comorbidities profile and improves models predicting clinical outcomes.

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

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

Protective effect of club cell secretory protein (CC-16) on COPD risk and progression: a Mendelian randomisation study

BACKGROUND

The anti-inflammatory pneumoprotein club cell secretory protein-16 (CC-16) is associated with the clinical expression of chronic obstructive pulmonary disease (COPD). We aimed to determine if there is a causal effect of serum CC-16 level on the risk of having COPD and/or its progression using Mendelian randomisation (MR) analysis.

METHODS

We performed a genome-wide association meta-analysis for serum CC-16 in two COPD cohorts (Lung Health Study (LHS), n=3850 and ECLIPSE, n=1702). We then used the CC-16-associated single-nucleotide polymorphisms (SNPs) as instrumental variables in MR analysis to identify a causal effect of serum CC-16 on 'COPD risk' (ie, case status in the International COPD Genetics Consortium/UK-Biobank dataset; n=35 735 COPD cases, n=222 076 controls) and 'COPD progression' (ie, annual change in forced expiratory volume in 1 s in LHS and ECLIPSE). We also determined the associations between SNPs associated with CC-16 and gene expression using n=1111 lung tissue samples from the Lung Expression Quantitative Trait Locus Study.

RESULTS

We identified seven SNPs independently associated (p<5×10^-8) with serum CC-16 levels; six of these were novel. MR analysis suggested a protective causal effect of increased serum CC-16 on COPD risk (MR estimate (SE) -0.11 (0.04), p=0.008) and progression (LHS only, MR estimate (SE) 7.40 (3.28), p=0.02). Five of the SNPs were also associated with gene expression in lung tissue (at false discovery rate <0.1) of several genes, including the CC-16-encoding gene SCGB1A1.

CONCLUSION

We have identified several novel genetic variants associated with serum CC-16 level in COPD cohorts. These genetic associations suggest a potential causal effect of serum CC-16 on the risk of having COPD and its progression, the biological basis of which warrants further investigation.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

A pseudolikelihood approach for assessing genetic association in case-control studies with unmeasured population structure

The case-control study design is one of the main tools for detecting associations between genetic markers and diseases. It is well known that population substructure can lead to spurious association between disease status and a genetic marker if the prevalence of disease and the marker allele frequency vary across subpopulations. In this paper, we propose a novel statistical method to estimate the association in case-control studies with unmeasured population substructure. The proposed method takes two steps. First, the information on genomic markers and disease status is used to infer the population substructure; second, the association between the disease and the test marker adjusting for the population substructure is modeled and estimated parametrically through polytomous logistic regression. The performance of the proposed method, relative to the existing methods, on bias, coverage probability and computational time, is assessed through simulations. The method is applied to an end-stage renal disease study in African Americans population.

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.

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.

Genetic Advances in Chronic Obstructive Pulmonary Disease. Insights from COPDGene

Chronic obstructive pulmonary disease (COPD) is a common and progressive disease that is influenced by both genetic and environmental factors. For many years, knowledge of the genetic basis of COPD was limited to Mendelian syndromes, such as alpha-1 antitrypsin deficiency and cutis laxa, caused by rare genetic variants. Over the past decade, the proliferation of genome-wide association studies, the accessibility of whole-genome sequencing, and the development of novel methods for analyzing genetic variation data have led to a substantial increase in the understanding of genetic variants that play a role in COPD susceptibility and COPD-related phenotypes. COPDGene (Genetic Epidemiology of COPD), a multicenter, longitudinal study of over 10,000 current and former cigarette smokers, has been pivotal to these breakthroughs in understanding the genetic basis of COPD. To date, over 20 genetic loci have been convincingly associated with COPD affection status, with additional loci demonstrating association with COPD-related phenotypes such as emphysema, chronic bronchitis, and hypoxemia. In this review, we discuss the contributions of the COPDGene study to the discovery of these genetic associations as well as the ongoing genetic investigations of COPD subtypes, protein biomarkers, and post-genome-wide association study analysis.

Whole genome sequencing of orofacial cleft trios from the Gabriella Miller Kids First Pediatric Research Consortium identifies a new locus on chromosome 21

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic, and the genetic etiology of non-syndromic OFCs is only partially determined. Here, we analyze whole genome sequence (WGS) data for association with risk of OFCs in European and Colombian families selected from a multicenter family-based OFC study. This is the first large-scale WGS study of OFC in parent-offspring trios, and a part of the Gabriella Miller Kids First Pediatric Research Program created for the study of childhood cancers and structural birth defects. WGS provides deeper and more specific genetic data than using imputation on present-day single nucleotide polymorphic (SNP) marker panels. Genotypes of case-parent trios at single nucleotide variants (SNV) and short insertions and deletions (indels) spanning the entire genome were called from their sequences using human GRCh38 genome assembly, and analyzed for association using the transmission disequilibrium test. Among genome-wide significant associations, we identified a new locus on chromosome 21 in Colombian families, not previously observed in other larger OFC samples of Latin American ancestry. This locus is situated within a region known to be expressed during craniofacial development. Based on deeper investigation of this locus, we concluded that it contributed risk for OFCs exclusively in the Colombians. This study reinforces the ancestry differences seen in the genetic etiology of OFCs, and underscores the need for larger samples when studying for OFCs and other birth defects in populations with diverse ancestry.

The Evolving Field of Genetic Epidemiology: From Familial Aggregation to Genomic Sequencing

The field of genetic epidemiology is relatively young and brings together genetics, epidemiology, and biostatistics to identify and implement the best study designs and statistical analyses for identifying genes controlling risk for complex and heterogeneous diseases (i.e., those where genes and environmental risk factors both contribute to etiology). The field has moved quickly over the past 40 years partly because the technology of genotyping and sequencing has forced it to adapt while adhering to the fundamental principles of genetics. In the last two decades, the available tools for genetic epidemiology have expanded from a genetic focus (considering 1 gene at a time) to a genomic focus (considering the entire genome), and now they must further expand to integrate information from other “-omics” (e.g., epigenomics, transcriptomics as measured by RNA expression) at both the individual and the population levels. Additionally, we can now also evaluate gene and environment interactions across populations to better understand exposure and the heterogeneity in disease risk. The future challenges facing genetic epidemiology are considerable both in scale and techniques, but the importance of the field will not diminish because by design it ties scientific goals with public health applications.

The pharmacogenomics of inhaled corticosteroids and lung function decline in COPD

Inhaled corticosteroids (ICS) are widely prescribed for patients with chronic obstructive pulmonary disease (COPD), yet have variable outcomes and adverse reactions, which may be genetically determined. The primary aim of the study was to identify the genetic determinants for forced expiratory volume in 1 s (FEV₁) changes related to ICS therapy.In the Lung Health Study (LHS)-2, 1116 COPD patients were randomised to the ICS triamcinolone acetonide (n=559) or placebo (n=557) with spirometry performed every 6 months for 3 years. We performed a pharmacogenomic genome-wide association study for the genotype-by-ICS treatment effect on 3 years of FEV₁ changes (estimated as slope) in 802 genotyped LHS-2 participants. Replication was performed in 199 COPD patients randomised to the ICS, fluticasone or placebo.A total of five loci showed genotype-by-ICS interaction at p<5×10^-6; of these, single nucleotide polymorphism (SNP) rs111720447 on chromosome 7 was replicated (discovery p=4.8×10^-6, replication p=5.9×10^-5) with the same direction of interaction effect. ENCODE (Encyclopedia of DNA Elements) data revealed that in glucocorticoid-treated (dexamethasone) A549 alveolar cell line, glucocorticoid receptor binding sites were located near SNP rs111720447. In stratified analyses of LHS-2, genotype at SNP rs111720447 was significantly associated with rate of FEV₁ decline in patients taking ICS (C allele β 56.36 mL·year^-1, 95% CI 29.96-82.76 mL·year^-1) and in patients who were assigned to placebo, although the relationship was weaker and in the opposite direction to that in the ICS group (C allele β -27.57 mL·year^-1, 95% CI -53.27- -1.87 mL·year^-1).The study uncovered genetic factors associated with FEV₁ changes related to ICS in COPD patients, which may provide new insight on the potential biology of steroid responsiveness in COPD.

COPDGene® 2019: Redefining the Diagnosis of Chronic Obstructive Pulmonary Disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) remains a major cause of morbidity and mortality. Present-day diagnostic criteria are largely based solely on spirometric criteria. Accumulating evidence has identified a substantial number of individuals without spirometric evidence of COPD who suffer from respiratory symptoms and/or increased morbidity and mortality. There is a clear need for an expanded definition of COPD that is linked to physiologic, structural (computed tomography [CT]) and clinical evidence of disease. Using data from the COPD Genetic Epidemiology study (COPDGene®), we hypothesized that an integrated approach that includes environmental exposure, clinical symptoms, chest CT imaging and spirometry better defines disease and captures the likelihood of progression of respiratory obstruction and mortality.

METHODS

Four key disease characteristics - environmental exposure (cigarette smoking), clinical symptoms (dyspnea and/or chronic bronchitis), chest CT imaging abnormalities (emphysema, gas trapping and/or airway wall thickening), and abnormal spirometry - were evaluated in a group of 8784 current and former smokers who were participants in COPDGene® Phase 1. Using these 4 disease characteristics, 8 categories of participants were identified and evaluated for odds of spirometric disease progression (FEV1 > 350 ml loss over 5 years), and the hazard ratio for all-cause mortality was examined.

RESULTS

Using smokers without symptoms, CT imaging abnormalities or airflow obstruction as the reference population, individuals were classified as Possible COPD, Probable COPD and Definite COPD. Current Global initiative for obstructive Lung Disease (GOLD) criteria would diagnose 4062 (46%) of the 8784 study participants with COPD. The proposed COPDGene® 2019 diagnostic criteria would add an additional 3144 participants. Under the new criteria, 82% of the 8784 study participants would be diagnosed with Possible, Probable or Definite COPD. These COPD groups showed increased risk of disease progression and mortality. Mortality increased in patients as the number of their COPD characteristics increased, with a maximum hazard ratio for all cause-mortality of 5.18 (95% confidence interval [CI]: 4.15-6.48) in those with all 4 disease characteristics.

CONCLUSIONS

A substantial portion of smokers with respiratory symptoms and imaging abnormalities do not manifest spirometric obstruction as defined by population normals. These individuals are at significant risk of death and spirometric disease progression. We propose to redefine the diagnosis of COPD through an integrated approach using environmental exposure, clinical symptoms, CT imaging and spirometric criteria. These expanded criteria offer the potential to stimulate both current and future interventions that could slow or halt disease progression in patients before disability or irreversible lung structural changes develop.

A systematic genetic analysis and visualization of phenotypic heterogeneity among orofacial cleft GWAS signals

Phenotypic heterogeneity is a hallmark of complex traits, and genetic studies of such traits may focus on them as a single diagnostic entity or by analyzing specific components. For example, in orofacial clefting (OFC), three subtypes-cleft lip (CL), cleft lip and palate (CLP), and cleft palate (CP) have been studied separately and in combination. To further dissect the genetic architecture of OFCs and how a given associated locus may be contributing to distinct subtypes of a trait we developed a framework for quantifying and interpreting evidence of subtype-specific or shared genetic effects in complex traits. We applied this technique to create a "cleft map" of the association of 30 genetic loci with three OFC subtypes. In addition to new associations, we found loci with subtype-specific effects (e.g., GRHL3 [CP], WNT5A [CLP]), as well as loci associated with two or all three subtypes. We cross-referenced these results with mouse craniofacial gene expression datasets, which identified additional promising candidate genes. However, we found no strong correlation between OFC subtypes and expression patterns. In aggregate, the cleft map revealed that neither subtype-specific nor shared genetic effects operate in isolation in OFC architecture. Our approach can be easily applied to any complex trait with distinct phenotypic subgroups.

Gene-gene interaction among cell adhesion genes and risk of nonsyndromic cleft lip with or without cleft palate in Chinese case-parent trios

Background

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with complex etiology. One strategy for studying the genetic risk factors of NSCL/P is to consider gene–gene interaction (G × G) among gene pathways having a role in craniofacial development. The present study aimed to investigate the G × G among cell adhesion gene pathway.

Methods

We carried out an interaction analysis of eight genes involved in cell adherens junctions among 806 NSCL/P Chinese case‐parent trios originally recruited for a genome‐wide association study (GWAS). Regression‐based approach was used to test for two‐way G × G interaction, while machine learning algorithm was run for exploring both two‐way and multi‐way interaction that may affect the risk of NSCL/P.

Results

A two‐way ACTN1 × CTNNB1 interaction reached the adjusted significance level. The single nucleotide polymorphisms pair composed of rs17252114 (CTNNB1) and rs1274944 (ACTN1) yielded a p value of .0002, and this interaction was also supported by the logic regression algorithm. Higher order interactions involving ACTN1, CTNNB1, and CDH1 were picked out by logic regression, suggesting a potential role in NSCL/P risk.

Conclusion

This study suggests for the first time evidence of both two‐way and multi‐way G × G interactions among cell adhesion genes contributing to the NSCL/P risk.

Exploring the interaction between FGF Genes and T-box genes among chinese nonsyndromic cleft lip with or without cleft palate case-parent trios

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect. Genetic variants causing syndromic orofacial clefts can also contribute to the etiology of NSCL/P. The purpose of the present study was to explore gene-gene (G × G) interaction using common single nucleotide polymorphic (SNP) markers in fibroblast growth factor (FGF) family and its receptors and T-box genes, which were associated with syndromic orofacial clefts. Our study was conducted in 806 Chinese NSCL/P case-parent trios drawn from an international consortium. A total of 252 SNPs in FGF8, FGF10, FGFR1, FGFR2, and TBX5 passed the quality control criteria and were included in the analysis. The interactions between SNPs in different genes were assessed using Cordell's method, which fitted a conditional logistic regression model. The analysis was performed using the R-package trio (Version 3.8.0). Bonferroni correction was used to adjust for multiple comparisons, and the overall significance threshold was set as P = 1.98 × 10^-4 (0.05/252). Conditional logistic regression revealed the most significant interaction between rs2330542 in FGF10 and rs1946295 in TBX5, which remained significant (P = 9.63 × 10^-6 ) after Bonferroni correction. The relative risk of allele C in rs2330542 (FGF10) was 1.02 (95%CI 0.81-1.28), while the relative risk was 1.42 (95%CI 1.03-1.97) when the exposure was a combination of allele C in rs2330542 and allele A in rs1946295 (TBX5). Our findings confirmed the importance of considering G × G interaction when exploring the genetic risk factors of NSCL/P. Further investigations are warranted to validate the potential interaction and reveal the biological function of FGF10/FGFR2/TBX5. Environ. Mol. Mutagen. 2019. © 2019 Wiley Periodicals, Inc.

DSP variants may be associated with longitudinal change in quantitative emphysema

Background

Emphysema, characterized by lung destruction, is a key component of Chronic Obstructive Pulmonary Disease (COPD) and is associated with increased morbidity and mortality. Genome-wide association studies (GWAS) have identified multiple genetic factors associated with cross-sectional measures of quantitative emphysema, but the genetic determinants of longitudinal change in quantitative measures of emphysema remain largely unknown. Our study aims to identify genetic variants associated with longitudinal change in quantitative emphysema measured by computed tomography (CT) imaging.

Methods

We included current and ex-smokers from two longitudinal cohorts: COPDGene, a study of Non-Hispanic Whites (NHW) and African Americans (AA), and the Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points (ECLIPSE). We calculated annual change in two quantitative measures of emphysema based on chest CT imaging: percent low attenuation area (≤ − 950HU) (%LAA-950) and adjusted lung density (ALD). We conducted GWAS, separately in 3030 NHW and 1158 AA from COPDGene and 1397 Whites from ECLIPSE. We further explored effects of 360 previously reported variants and a lung function based polygenic risk score on annual change in quantitative emphysema.

Results

In the genome-wide association analysis, no variants achieved genome-wide significance (P < 5e-08). However, in the candidate region analysis, rs2076295 in the DSP gene, previously associated with COPD, lung function and idiopathic pulmonary fibrosis, was associated with change in %LAA-950 (β (SE) = 0.09 (0.02), P = 3.79e-05) and in ALD (β (SE) = − 0.06 (0.02), P = 2.88e-03). A lung function based polygenic risk score was associated with annual change in %LAA-950 (P = 4.03e-02) and with baseline measures of quantitative emphysema (P < 1e-03) and showed a trend toward association with annual change in ALD (P = 7.31e-02).

Conclusions

DSP variants may be associated with longitudinal change in quantitative emphysema. Additional investigation of the DSP gene are likely to provide further insights into the disease progression in emphysema and COPD.

Trial registration

Clinicaltrials.gov Identifier: NCT00608764, NCT00292552.

Electronic supplementary material

The online version of this article (10.1186/s12931-019-1097-8) contains supplementary material, which is available to authorized users.

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.

Predictors of dental care utilization in north-central Appalachia in the USA

OBJECTIVES

Dental utilization is an important determinant of oral health and well-being. The aim of this study was to evaluate potential associations between a variety of biopsychosocial factors and dental utilization in north-central Appalachia, USA, a region where oral health disparities are profound.

METHODS

This study used household-based data from the Center for Oral Health Research in Appalachia (COHRA1) study in north-central Appalachia, including 449 families with 868 adults. The generalized estimating equation (GEE) approach was used to determine the best-fitting predictor model for dental utilization among adult family members.

RESULTS

On average across West Virginia and Pennsylvania, having dental insurance was associated with greater dental utilization over a 3-year time period (OR = 2.20, 95% CI = 1.54, 3.14). When stratified by state, the association held for only West Virginia (OR = 2.41, 95% CI = 1.54, 3.79) and was nonsignificant for Pennsylvania residents (OR = 1.50, 95% CI = 0.80, 2.79). Individuals from Pennsylvania were more likely to utilize dental care and participants from West Virginia less so (2.31, 95% CI = 1.57, 3.40). Females from Pennsylvania were more likely than males to regularly seek dental care (OR = 1.44, 95% CI = 1.00, 2.05), and a higher income was associated with greater frequency of regular dental visits (OR = 1.21, 95% CI = 1.09, 1.34) in West Virginia. Individuals from Pennsylvania who scored higher on the Physiological Arousal subscale of the Dental Fear Survey were more likely to attend routine care visits (OR = 1.18, 95% CI = 1.03, 1.35). Across both states, more fatalistic beliefs related to oral health care also predicted less routine care (OR = 0.87, 95% CI = 0.81, 0.94), and more investment in or more positive attitudes towards one's oral health also was associated with higher utilization (OR = 1.18, 95% CI = 1.13, 1.23).

CONCLUSIONS

Overall, the findings of this study suggest state residency, sex, insurance, income, fatalistic beliefs, health values, and aspects of dental care-related anxiety and fear predicted dental care utilization in north-central Appalachia. These findings reinforce the need to address insurance and other economic factors affecting utilization and to consider how individual-level fatalistic beliefs and oral health values may affect utilization of routine oral health care.

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.