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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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 (Padjusted<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 Padjusted<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.

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.

The combined finding of HPV 16, 18, or 45 and cytologic Atypical Glandular Cells (AGC) indicates a greatly elevated risk of in situ and invasive cervical adenocarcinoma.

BACKGROUND: Cervical screening has not effectively controlled cervical adenocarcinoma (AC). Human papillomavirus (HPV) testing is recommended for cervical screening but the optimal management of HPV-positive individuals to prevent AC remains a question. Cytology and HPV typing are two triage options to predict the risk of AC. We combined two potential biomarkers (atypical glandular cell, AGC, cytology and HPV-types 16, 18, or 45) to assess their joint effect on detecting AC. METHODS: Kaiser Permanente Northern California (KPNC) used triennial co-testing with cytology and HPV testing (positive/negative) for routine cervical screening between 2003 and 2020. HPV typing of a sample of residual HPV test specimens was performed on a separate cohort selected from KPNC (Persistence and Progression, PaP, cohort). We compared risk of prevalent and incident histologic AC/AIS (adenocarcinoma in situ) associated with preceding combinations of cytologic results and HPV typing. Risk of squamous cell cancer (SCC)/cervical intraepithelial neoplasia grade 3 (CIN3) (SCC/CIN3) was also included for comparison. RESULTS: Among HPV-positive individuals in PaP cohort, 99% of prevalent AC and 96% of AIS were linked to HPV-types 16, 18, or 45 (denoted HPV 16/18/45). Although rare (0.09% of screening population), the concurrent detection of HPV 16/18/45 with AGC cytology predicted a highly elevated relative risk of underlying histologic AC/AIS; the absolute risk of diagnosing AC/AIS was 12% and odds ratio (OR) was 1341 (95%CI:495-3630) compared to patients with other high-risk HPV types and normal cytology. Cumulatively (allowing non-concurrent results), approximately one-third of the AC/AIS cases ever had HPV 16/18/45 and AGC cytology (OR = 1785; 95%CI:872-3656). AGC was not as strongly associated with SCC/CIN3. CONCLUSION: Detection of HPV 16/18/45 positivity elevates risk of adenocarcinoma, particularly if AGC cytology is also found.

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.

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.

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.

Analysis of genetically driven alternative splicing identifies FBXO38 as a novel COPD susceptibility gene.

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p<1x10-6) in a published GWAS have at least one shared causal variant between the GWAS and sQTL studies. Among the genes identified to have splice sites associated with top GWAS SNPs was FBXO38, in which a novel exon was discovered to be protective against COPD. Importantly, the sQTL in this locus was validated by qPCR in both blood and lung tissue, demonstrating that splice variants relevant to lung tissue can be identified in blood. Other identified genes included CDK11A and SULT1A2. Overall, these data indicate that analysis of alternative splicing can provide novel insights into disease mechanisms. In particular, we demonstrated that SNPs in a known COPD GWAS locus on chromosome 5q32 influence alternative splicing in the gene FBXO38.

Multiethnic genome-wide and HLA association study of total serum IgE level.

BACKGROUND: Total serum IgE (tIgE) is an important intermediate phenotype of allergic disease. Whole genome genetic association studies across ancestries may identify important determinants of IgE. OBJECTIVE: We aimed to increase understanding of genetic variants affecting tIgE production across the ancestry and allergic disease spectrum by leveraging data from the National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine program; the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA); and the Atopic Dermatitis Research Network (N = 21,901). METHODS: We performed genome-wide association within strata of study, disease, and ancestry groups, and we combined results via a meta-regression approach that models heterogeneity attributable to ancestry. We also tested for association between HLA alleles called from whole genome sequence data and tIgE, assessing replication of associations in HLA alleles called from genotype array data. RESULTS: We identified 6 loci at genome-wide significance (P < 5 × 10-9), including 4 loci previously reported as genome-wide significant for tIgE, as well as new regions in chr11q13.5 and chr15q22.2, which were also identified in prior genome-wide association studies of atopic dermatitis and asthma. In the HLA allele association study, HLA-A∗02:01 was associated with decreased tIgE level (Pdiscovery = 2 × 10-4; Preplication = 5 × 10-4; Pdiscovery+replication = 4 × 10-7), and HLA-DQB1∗03:02 was strongly associated with decreased tIgE level in Hispanic/Latino ancestry populations (PHispanic/Latino discovery+replication = 8 × 10-8). CONCLUSION: We performed the largest genome-wide association study and HLA association study of tIgE focused on ancestrally diverse populations and found several known tIgE and allergic disease loci that are relevant in non-European ancestry populations.

Damaging Mutations in AFDN 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.