Lacticaseibacillus rhamnosus GG Improves Periodontal Bone Repair via Gut-Blood Axis in Hyperlipidemia

Periodontal bone regeneration remains a clinical challenge, and hyperlipidemia can aggravate alveolar bone resorption. Probiotics have recently been reported to improve bone mass. We aimed to determine the role of Lacticaseibacillus rhamnosus GG (LGG) in periodontal bone regeneration improvement within the context of periodontitis with hyperlipidemia. A Sprague Dawley rat model for periodontitis, hyperlipidemia, and periodontal fenestration defect was constructed (n = 36) and administered LGG gavage for 6 wk (the rats were subsequently sacrificed). Fecal microbiota from donor rats 3 wk after LGG gavage was transplanted into recipient rats to evaluate the role of LGG-modulated gut microbiota in periodontal bone regeneration. Regenerated bone mass was detected using micro-computerized tomography and hematoxylin and eosin stain. Gut microbiota was analyzed using 16S ribosomal RNA sequencing. Serum metabolites were detected by liquid chromatography-mass spectrometry (6 wk after LGG gavage). The pro-osteogenic effects of screened serum metabolite were verified in vitro on bone marrow mesenchymal stem cells (BMMSCs). We found that the bone mineral density, bone volume (BV), trabecular bone volume fraction (BV/TV), and trabecular thickness of the regenerated periodontal bone increased after LGG gavage (P < 0.05) but had little effect on oral flora. After LGG gavage, Staphylococcus, Corynebacterium, and Collinsella in the gut of donors were significantly changed, and these differences were maintained in recipients, who also showed increased trabecular thickness of the regenerated periodontal bone (P < 0.05). These key genera were correlated with BV/TV and BV (P < 0.05). In addition, LGG gavage significantly regulated bone-related blood metabolites, of which selenomethionine promoted BMMSC osteogenesis. Notably, selenomethionine was associated with key gut genera (P < 0.05). Collectively, LGG improved periodontal bone regeneration in the context of periodontitis with hyperlipidemia by modulating gut microbiota and increasing pro-osteogenic metabolites in the blood. These results reveal new insights into the use of probiotics to promote periodontal bone regeneration via the gut-blood-bone axis.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

[Structural characteristics of lower respiratory tract microflora in patients with pneumoconiosis]

Objective: To explore the composition of bacteria in lower respiratory tract of patients with pneumoconiosis and dust exposure, and to compare and analyze the difference and correlation between them. Methods: From May 2020 to January 2021, a prospective multicenter cross-sectional study was conducted to select patients with pneumoconiosis who underwent bronchoalveolar lavage treatment at the Respiratory and Critical Care Medical Department of the 920th Hospital of the Joint Support Force and the Respiratory Department of Tongren Hospital in Kunming, as well as the population of dust recipients. A total of 24 patients with pneumoconiosis (pneumoconiosis group) were included, and 16 dust exposed individuals (dust exposed group) were used as controls. Two groups of patients' alveolar lavage fluid were collected. The 16SrRNA gene V3-V4 sequencing technology and bioinformatics analysis platform were used to measure and analyze the differences in microbial structure composition and associations between bacterial communities. Results: Compared with the dust exposed group, the top 5 bacterial phyla in the alveolar lavage fluid level of patients with pneumoconiosis were the same, followed by Proteobacteria, Firmicutes, Bacteroidetes, Fusobacteria, and Actinobacteria. Compared with the dust exposure group, the pneumoconiosis group patients belong to the top 5 genera of horizontal flora abundance, which are different. The dust exposure group is respectively: Pseudomonas, Proctor, Streptococcus, Achromobacter, and Neisseria. The pneumoconiosis group is respectively: Pseudomonas, Achromobacter, Streptococcus, Ralstonia, and Proctor. The Alpha diversity analysis results showed that compared with the dust exposed group, the level of bacterial diversity in the pneumoconiosis group was difference (P<0.05), and there was no statistically significant difference in bacterial evenness (P>0.05) ; Beta diversity showed differences in microbial community structure between the two groups (P<0.05 ). Single factor microbial association network analysis showed that there was a high correlation between Firmicutes and Bacteroidetes in the pneumoconiosis and dust exposed groups and other species, showing a positive correlation; The correlation between Proteobacteria and other species is high, showing a negative correlation. Conclusion: The structure and relative abundance of bacteria in lower respiratory tract were different between patients with pneumoconiosis and dust exposure, and the diversity of bacteria in lower respiratory tract increased in patients with pneumoconiosis, which may be related to disease status.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Assessing the contribution of rare genetic variants to phenotypes of chronic obstructive pulmonary disease using whole-genome sequence data

RATIONALE

Genetic variation has a substantial contribution to chronic obstructive pulmonary disease (COPD) and lung function measurements. Heritability estimates using genome-wide genotyping data can be biased if analyses do not appropriately account for the nonuniform distribution of genetic effects across the allele frequency and linkage disequilibrium (LD) spectrum. In addition, the contribution of rare variants has been unclear.

OBJECTIVES

We sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program.

METHODS

Using the genome-based restricted maximum likelihood method, we partitioned the genome into bins based on minor allele frequency and LD scores and estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio in 11 051 European ancestry and 5853 African-American participants.

MEASUREMENTS AND MAIN RESULTS

In European ancestry participants, the estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio were 35.5%, 55.6% and 32.5%, of which 18.8%, 19.7%, 17.8% were from common variants, and 16.6%, 35.8%, and 14.6% were from rare variants. These estimates had wide confidence intervals, with common variants and some sets of rare variants showing a statistically significant contribution (P-value < 0.05). In African-Americans, common variant heritability was similar to European ancestry participants, but lower sample size precluded calculation of rare variant heritability.

CONCLUSIONS

Our study provides updated and unbiased estimates of heritability for COPD and lung function, and suggests an important contribution of rare variants. Larger studies of more diverse ancestry will improve accuracy of these estimates.

Rare genetic variants explain missing heritability in smoking

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

A statistical physics approach for disease module detection

Extensive evidence indicates that the pathobiological processes of a complex disease are associated with perturbation in specific neighborhoods of the human protein-protein interaction (PPI) network (also known as the interactome), often referred to as the disease module. Many computational methods have been developed to integrate the interactome and omics profiles to extract context-dependent disease modules. Yet, existing methods all have fundamental limitations in terms of rigor and/or efficiency. Here, we developed a statistical physics approach based on the random-field Ising model (RFIM) for disease module detection, which is both mathematically rigorous and computationally efficient. We applied our RFIM approach to genome-wide association studies (GWAS) of ten complex diseases to examine its performance for disease module detection. We found that our RFIM approach outperforms existing methods in terms of computational efficiency, connectivity of disease modules, and robustness to the interactome incompleteness.

A polygenic risk score and age of diagnosis of COPD

BACKGROUND

Genetic susceptibility may be associated with earlier onset of chronic obstructive pulmonary disease (COPD). We hypothesised that a polygenic risk score (PRS) for COPD would be associated with earlier age of diagnosis of COPD.

METHODS

In 6647 non-Hispanic White (NHW) and 2464 African American (AA) participants from COPDGene, and 6812 participants from the Framingham Heart Study (FHS), we tested the relationship of the PRS and age of COPD diagnosis. Age at diagnosis was determined by: 1) self-reported age at COPD diagnosis or 2) age at visits when moderate-to-severe airflow limitation (Global Initiative for Chronic Obstructive Lung Disease (GOLD) grade 2-4) was observed on spirometry. We used Cox regression to examine the overall and time-dependent effects of the PRS on incident COPD. In the COPDGene study, we also examined the PRS's predictive value for COPD at age <50 years (COPD50) using logistic regression and area under the curve (AUC) analyses, with and without the addition of other risk factors present at early life (e.g. childhood asthma).

RESULTS

In Cox models, the PRS demonstrated age-dependent associations with incident COPD, with larger effects at younger ages in both cohorts. The PRS was associated with COPD50 (OR 1.55 (95% CI 1.41-1.71) for NHW, OR 1.23 (95% CI 1.05-1.43) for AA and OR 2.47 (95% CI 2.12-2.88) for FHS participants). In COPDGene, adding the PRS to known early-life risk factors improved prediction of COPD50 in NHW (AUC 0.69 versus 0.74; p<0.0001) and AA (AUC 0.61 versus 0.64; p=0.04) participants.

CONCLUSIONS

A COPD PRS is associated with earlier age of diagnosis of COPD and retains predictive value when added to known early-life risk factors.

Interstitial lung abnormalities are associated with decreased mean telomere length

BACKGROUND

Interstitial lung abnormalities (ILA) share many features with idiopathic pulmonary fibrosis; however, it is not known if ILA are associated with decreased mean telomere length (MTL).

METHODS

Telomere length was measured with quantitative PCR in the Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) and Age Gene/Environment Susceptibility Reykjavik (AGES-Reykjavik) cohorts and Southern blot analysis was used in the Framingham Heart Study (FHS). Logistic and linear regression were used to assess the association between ILA and MTL; Cox proportional hazards models were used to assess the association between MTL and mortality.

RESULTS

In all three cohorts, ILA were associated with decreased MTL. In the COPDGene and AGES-Reykjavik cohorts, after adjustment there was greater than twofold increase in the odds of ILA when comparing the shortest quartile of telomere length to the longest quartile (OR 2.2, 95% CI 1.5-3.4, p=0.0001, and OR 2.6, 95% CI 1.4-4.9, p=0.003, respectively). In the FHS, those with ILA had shorter telomeres than those without ILA (-767 bp, 95% CI 76-1584 bp, p=0.03). Although decreased MTL was associated with chronic obstructive pulmonary disease (OR 1.3, 95% CI 1.1-1.6, p=0.01) in COPDGene, the effect estimate was less than that noted with ILA. There was no consistent association between MTL and risk of death when comparing the shortest quartile of telomere length in COPDGene and AGES-Reykjavik (HR 0.82, 95% CI 0.4-1.7, p=0.6, and HR 1.2, 95% CI 0.6-2.2, p=0.5, respectively).

CONCLUSION

ILA are associated with decreased MTL.

C FTR variants are associated with chronic bronchitis in smokers

INTRODUCTION

Loss-of-function variants in both copies of the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF); however, there is evidence that reduction in CFTR function due to the presence of one deleterious variant can have clinical consequences. Here, we hypothesise that CFTR variants in individuals with a history of smoking are associated with chronic obstructive pulmonary disease (COPD) and related phenotypes.

METHODS

Whole-genome sequencing was performed through the National Heart, Lung, and Blood Institute TOPMed (TransOmics in Precision Medicine) programme in 8597 subjects from the COPDGene (Genetic Epidemiology of COPD) study, an observational study of current and former smokers. We extracted clinically annotated CFTR variants and performed single-variant and variant-set testing for COPD and related phenotypes. Replication was performed in 2118 subjects from the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) study.

RESULTS

We identified 301 coding variants within the CFTR gene boundary: 147 of these have been reported in individuals with CF, including 36 CF-causing variants. We found that CF-causing variants were associated with chronic bronchitis in variant-set testing in COPDGene (one-sided p=0.0025; OR 1.53) and in meta-analysis of COPDGene and ECLIPSE (one-sided p=0.0060; OR 1.52). Single-variant testing revealed that the F508del variant was associated with chronic bronchitis in COPDGene (one-sided p=0.015; OR 1.47). In addition, we identified 32 subjects with two or more CFTR variants on separate alleles and these subjects were enriched for COPD cases (p=0.010).

CONCLUSIONS

Cigarette smokers who carry one deleterious CFTR variant have higher rates of chronic bronchitis, while presence of two CFTR variants may be associated with COPD. These results indicate that genetically mediated reduction in CFTR function contributes to COPD related phenotypes, in particular chronic bronchitis.

Machine Learning Prediction of Progression in Forced Expiratory Volume in 1 Second in the COPDGene® Study

BACKGROUND

The heterogeneous nature of chronic obstructive pulmonary disease (COPD) complicates the identification of the predictors of disease progression. We aimed to improve the prediction of disease progression in COPD by using machine learning and incorporating a rich dataset of phenotypic features.

METHODS

We included 4496 smokers with available data from their enrollment and 5-year follow-up visits in the COPD Genetic Epidemiology (COPDGene^®) study. We constructed linear regression (LR) and supervised random forest models to predict 5-year progression in forced expiratory in 1 second (FEV₁) from 46 baseline features. Using cross-validation, we randomly partitioned participants into training and testing samples. We also validated the results in the COPDGene 10-year follow-up visit.

RESULTS

Predicting the change in FEV₁ over time is more challenging than simply predicting the future absolute FEV₁ level. For random forest, R-squared was 0.15 and the area under the receiver operator characteristic (ROC) curves for the prediction of participants in the top quartile of observed progression was 0.71 (testing) and respectively, 0.10 and 0.70 (validation). Random forest provided slightly better performance than LR. The accuracy was best for Global initiative for chronic Obstructive Lung Disease (GOLD) grades 1-2 participants, and it was harder to achieve accurate prediction in advanced stages of the disease. Predictive variables differed in their relative importance as well as for the predictions by GOLD.

CONCLUSION

Random forest, along with deep phenotyping, predicts FEV₁ progression with reasonable accuracy. There is significant room for improvement in future models. This prediction model facilitates the identification of smokers at increased risk for rapid disease progression. Such findings may be useful in the selection of patient populations for targeted clinical trials.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV₁] and its ratio to forced vital capacity [FEV₁/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV₁ and FEV₁/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10^-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Mendelian randomization supports bidirectional causality between telomere length and clonal hematopoiesis of indeterminate potential

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Association of clonal hematopoiesis with chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is associated with age and smoking, but other determinants of the disease are incompletely understood. Clonal hematopoiesis of indeterminate potential (CHIP) is a common, age-related state in which somatic mutations in clonal blood populations induce aberrant inflammatory responses. Patients with CHIP have an elevated risk for cardiovascular disease, but the association of CHIP with COPD remains unclear. We analyzed whole-genome sequencing and whole-exome sequencing data to detect CHIP in 48 835 patients, of whom 8444 had moderate to very severe COPD, from four separate cohorts with COPD phenotyping and smoking history. We measured emphysema in murine models in which Tet2 was deleted in hematopoietic cells. In the COPDGene cohort, individuals with CHIP had risks of moderate-to-severe, severe, or very severe COPD that were 1.6 (adjusted 95% confidence interval [CI], 1.1-2.2) and 2.2 (adjusted 95% CI, 1.5-3.2) times greater than those for noncarriers. These findings were consistently observed in three additional cohorts and meta-analyses of all patients. CHIP was also associated with decreased FEV1% predicted in the COPDGene cohort (mean between-group differences, -5.7%; adjusted 95% CI, -8.8% to -2.6%), a finding replicated in additional cohorts. Smoke exposure was associated with a small but significant increased risk of having CHIP (odds ratio, 1.03 per 10 pack-years; 95% CI, 1.01-1.05 per 10 pack-years) in the meta-analysis of all patients. Inactivation of Tet2 in mouse hematopoietic cells exacerbated the development of emphysema and inflammation in models of cigarette smoke exposure. Somatic mutations in blood cells are associated with the development and severity of COPD, independent of age and cumulative smoke exposure.

Negative associations of morning serum cortisol levels with obesity: the Henan rural cohort study

AIMS

To evaluate the associations of morning serum cortisol levels with obesity defined by different indices in Chinese rural populations.

MATERIALS AND METHODS

A cross-sectional study was performed including 6198 participants (2566 males and 3632 females). Serum cortisol was collected in morning and quantified by liquid chromatography-tandem mass spectrometry. Obesity was defined by body mass index (BMI), body fat percentage (BFP), waist-to-height ratio (WHtR), waist circumference (WC), visceral fat index (VFI) and waist-to-hip ratio (WHR). Both multivariable liner regression, logistic regression and restrictive cubic splines models were used to estimate the gender-specific relationships between cortisol levels and obesity defined by different indices, respectively.

RESULTS

After adjusting for potential confounders, serum cortisol was negatively associated with different obesity measures, except obese females defined by BFP (for instance, overall obesity defined by BMI, Quartile 4 vs. Quartile 1, odds ratio (OR) = 0.25, 95% confidence interval (CI):0.15, 0.41 in males, and OR = 0.58, 95% CI: 0.42,0.80 in females, central obesity defined by WC, OR = 0.52, 95% CI:0.39,0.69 in males and OR = 0.63, 95% CI:0.51,0.77 in females). Similarly, restrictive cubic splines showed the nonlinear relationship between high levels of cortisol and different obesity indices. Furthermore, ROC curve analysis indicated that cortisol could improve the discrimination of model with common biomarkers.

CONCLUSION

Morning serum cortisol were negatively related to obesity defined by different indices in Chinese rural populations. In addition, cortisol could be as a biomarker for prediction of obesity in males.

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.

Heterozygosity of the Alpha 1-Antitrypsin Pi*Z Allele and Risk of Liver Disease

The serpin family A member 1 (SERPINA1) Z allele is present in approximately one in 25 individuals of European ancestry. Z allele homozygosity (Pi*ZZ) is the most common cause of alpha 1-antitrypsin deficiency and is a proven risk factor for cirrhosis. We examined whether heterozygous Z allele (Pi*Z) carriers in United Kingdom (UK) Biobank, a population-based cohort, are at increased risk of liver disease. We replicated findings in Massachusetts General Brigham Biobank, a hospital-based cohort. We also examined variants associated with liver disease and assessed for gene-gene and gene-environment interactions. In UK Biobank, we identified 1,493 cases of cirrhosis, 12,603 Z allele heterozygotes, and 129 Z allele homozygotes among 312,671 unrelated white British participants. Heterozygous carriage of the Z allele was associated with cirrhosis compared to noncarriage (odds ratio [OR], 1.53; P = 1.1×10^-04); homozygosity of the Z allele also increased the risk of cirrhosis (OR, 11.8; P = 1.8 × 10^-09). The OR for cirrhosis of the Z allele was comparable to that of well-established genetic variants, including patatin-like phospholipase domain containing 3 (PNPLA3) I148M (OR, 1.48; P = 1.1 × 10^-22) and transmembrane 6 superfamily member 2 (TM6SF2) E167K (OR, 1.34; P = 2.6 × 10^-06). In heterozygotes compared to noncarriers, the Z allele was associated with higher alanine aminotransferase (ALT; P = = 4.6 × 10^-46), aspartate aminotransferase (AST; P = 2.2 × 10^-27), alkaline phosphatase (P = 3.3 × 10^-43), gamma-glutamyltransferase (P = 1.2 × 10^-05), and total bilirubin (P = 6.4 × 10^-06); Z allele homozygotes had even greater elevations in liver biochemistries. Body mass index (BMI) amplified the association of the Z allele for ALT (P interaction = 0.021) and AST (P interaction = 0.0040), suggesting a gene-environment interaction. Finally, we demonstrated genetic interactions between variants in PNPLA3, TM6SF2, and hydroxysteroid 17-beta dehydrogenase 13 (HSD17B13); there was no evidence of epistasis between the Z allele and these variants. Conclusion: SERPINA1 Z allele heterozygosity is an important risk factor for liver disease; this risk is amplified by increasing BMI.

[Analysis of Dimethyltryptamine in Captured Suspicious Substances：Three Case Reports]

Abstract.

Tempo-spatial regulation of the Wnt pathway by FAM13A modulates the stemness of alveolar epithelial progenitors

BACKGROUND

Family with Sequence Similarity 13, Member A (FAM13A) gene has been consistently associated with COPD by Genome-wide association studies (GWAS). Our previous study demonstrated that FAM13A was mainly expressed in the lung epithelial progenitors including Club cells and alveolar type II epithelial (ATII) cells. Fam13a^-/- mice were resistant to cigarette smoke (CS)-induced emphysema through promoting β-catenin/Wnt activation. Given the important roles of β-catenin/Wnt activation in alveolar regeneration during injury, it is unclear when and where FAM13A regulates the Wnt pathway, the requisite pathway for alveolar epithelial repair, in vivo during CS exposure in lung epithelial progenitors.

METHODS

Fam13a^+/+ or Fam13a^-/- mice were crossed with TCF/Lef:H2B-GFP Wnt-signaling reporter mouse line to indicate β-catenin/Wnt-activated cells labeled with GFP followed by acute (1 month) or chronic (7 months) CS exposure. Fluorescence-activated flow cytometry analysis, immunofluorescence and organoid culture system were performed to identify the β-catenin/Wnt-activated cells in Fam13a^+/+ or Fam13a^-/- mice exposed to CS. Fam13a;SftpcCreERT2;Rosa26RmTmG mouse line, where GFP labels ATII cells, was generated for alveolar organoid culture followed by analyses of organoid number, immunofluorescence and gene expression. Single cell RNA-seq data from COPD ever smokers and nonsmoker control lungs were further analyzed.

FINDINGS

We found that FAM13A-deficiency significantly increased Wnt activation mainly in lung epithelial cells. Consistently, after long-term CS exposure in vivo, FAM13A deficiency bestows alveolar epithelial progenitor cells with enhanced proliferation and differentiation in the ex vivo organoid model. Importantly, expression of FAM13A is significantly increased in human COPD-derived ATII cells compared to healthy ATII cells as suggested by single cell RNA-sequencing data.

INTERPRETATION

Our findings suggest that FAM13A-deficiency promotes the Wnt pathway-mediated ATII cell repair/regeneration, and thereby possibly mitigating CS-induced alveolar destruction. FUND: This project is funded by the National Institutes of Health of United States of America (NIH) grants R01HL127200, R01HL137927, R01HL148667 and R01HL147148 (XZ).

Connecting COPD GWAS genes: FAM13A controls TGFβ2 secretion by modulating AP-3 transport

Chronic Obstructive Pulmonary Disease (COPD) is a common, complex disease and a major cause of morbidity and mortality. Although multiple genetic determinants of COPD have been implicated by genome-wide association studies (GWAS), the pathophysiologic significance of these associations remains largely unknown. From a COPD protein-protein interaction network module, we selected a network path between two COPD GWAS genes for validation studies: FAM13A-AP3D1-CTGF-TGFB2. We find that TGFβ2, FAM13A, and AP3D1 (but not CTGF) form a cellular protein complex. Functional characterization suggests that this complex mediates the secretion of TGFβ2 through an AP-3-dependent pathway, with FAM13A acting as a negative regulator by targeting a late stage of this transport that involves the dissociation of coat-cargo interaction. Moreover, we find that TGFβ2 is a transmembrane protein that engages the AP-3 complex for delivery to the late endosomal compartments for subsequent secretion through exosomes. These results identify a pathophysiologic context that unifies the biological network role of two COPD GWAS proteins and reveal novel mechanisms of cargo transport through an intracellular pathway.

Genome-wide association analysis of COVID-19 mortality risk in SARS-CoV-2 genomes identifies mutation in the SARS-CoV-2 spike protein that colocalizes with P.1 of the Brazilian strain

SARS‐CoV‐2 mortality has been extensively studied in relation to host susceptibility. How sequence variations in the SARS‐CoV‐2 genome affect pathogenicity is poorly understood. Starting in October 2020, using the methodology of genome‐wide association studies (GWAS), we looked at the association between whole‐genome sequencing (WGS) data of the virus and COVID‐19 mortality as a potential method of early identification of highly pathogenic strains to target for containment. Although continuously updating our analysis, in December 2020, we analyzed 7548 single‐stranded SARS‐CoV‐2 genomes of COVID‐19 patients in the GISAID database and associated variants with mortality using a logistic regression. In total, evaluating 29,891 sequenced loci of the viral genome for association with patient/host mortality, two loci, at 12,053 and 25,088 bp, achieved genome‐wide significance (p values of 4.09e−09 and 4.41e−23, respectively), though only 25,088 bp remained significant in follow‐up analyses. Our association findings were exclusively driven by the samples that were submitted from Brazil (p value of 4.90e−13 for 25,088 bp). The mutation frequency of 25,088 bp in the Brazilian samples on GISAID has rapidly increased from about 0.4 in October/December 2020 to 0.77 in March 2021. Although GWAS methodology is suitable for samples in which mutation frequencies varies between geographical regions, it cannot account for mutation frequencies that change rapidly overtime, rendering a GWAS follow‐up analysis of the GISAID samples that have been submitted after December 2020 as invalid. The locus at 25,088 bp is located in the P.1 strain, which later (April 2021) became one of the distinguishing loci (precisely, substitution V1176F) of the Brazilian strain as defined by the Centers for Disease Control. Specifically, the mutations at 25,088 bp occur in the S2 subunit of the SARS‐CoV‐2 spike protein, which plays a key role in viral entry of target host cells. Since the mutations alter amino acid coding sequences, they potentially imposing structural changes that could enhance viral infectivity and symptom severity. Our analysis suggests that GWAS methodology can provide suitable analysis tools for the real‐time detection of new more transmissible and pathogenic viral strains in databases such as GISAID, though new approaches are needed to accommodate rapidly changing mutation frequencies over time, in the presence of simultaneously changing case/control ratios. Improvements of the associated metadata/patient information in terms of quality and availability will also be important to fully utilize the potential of GWAS methodology in this field.

An algorithm to simulate nonstationary and non-Gaussian stochastic processes

We proposed a new iterative power and amplitude correction (IPAC) algorithm to simulate nonstationary and non-Gaussian processes. The proposed algorithm is rooted in the concept of defining the stochastic processes in the transform domain, which is elaborated and extend. The algorithm extends the iterative amplitude adjusted Fourier transform algorithm for generating surrogate and the spectral correction algorithm for simulating stationary non-Gaussian process. The IPAC algorithm can be used with different popular transforms, such as the Fourier transform, S-transform, and continuous wavelet transforms. The targets for the simulation are the marginal probability distribution function of the process and the power spectral density function of the process that is defined based on the variables in the transform domain for the adopted transform. The algorithm is versatile and efficient. Its application is illustrated using several numerical examples.

Robust, flexible, and scalable tests for Hardy-Weinberg equilibrium across diverse ancestries

Traditional Hardy-Weinberg equilibrium (HWE) tests (the χ2 test and the exact test) have long been used as a metric for evaluating genotype quality, as technical artifacts leading to incorrect genotype calls often can be identified as deviations from HWE. However, in data sets composed of individuals from diverse ancestries, HWE can be violated even without genotyping error, complicating the use of HWE testing to assess genotype data quality. In this manuscript, we present the Robust Unified Test for HWE (RUTH) to test for HWE while accounting for population structure and genotype uncertainty, and to evaluate the impact of population heterogeneity and genotype uncertainty on the standard HWE tests and alternative methods using simulated and real sequence data sets. Our results demonstrate that ignoring population structure or genotype uncertainty in HWE tests can inflate false-positive rates by many orders of magnitude. Our evaluations demonstrate different tradeoffs between false positives and statistical power across the methods, with RUTH consistently among the best across all evaluations. RUTH is implemented as a practical and scalable software tool to rapidly perform HWE tests across millions of markers and hundreds of thousands of individuals while supporting standard VCF/BCF formats. RUTH is publicly available at https://www.github.com/statgen/ruth.

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

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

A linear prognostic score based on the ratio of interleukin-6 to interleukin-10 predicts outcomes in COVID-19

BACKGROUND

Prognostic tools are required to guide clinical decision-making in COVID-19.

METHODS

We studied the relationship between the ratio of interleukin (IL)-6 to IL-10 and clinical outcome in 80 patients hospitalized for COVID-19, and created a simple 5-point linear score predictor of clinical outcome, the Dublin-Boston score. Clinical outcome was analysed as a three-level ordinal variable ("Improved", "Unchanged", or "Declined"). For both IL-6:IL-10 ratio and IL-6 alone, we associated clinical outcome with a) baseline biomarker levels, b) change in biomarker level from day 0 to day 2, c) change in biomarker from day 0 to day 4, and d) slope of biomarker change throughout the study. The associations between ordinal clinical outcome and each of the different predictors were performed with proportional odds logistic regression. Associations were run both "unadjusted" and adjusted for age and sex. Nested cross-validation was used to identify the model for incorporation into the Dublin-Boston score.

FINDINGS

The 4-day change in IL-6:IL-10 ratio was chosen to derive the Dublin-Boston score. Each 1 point increase in the score was associated with a 5.6 times increased odds for a more severe outcome (OR 5.62, 95% CI -3.22-9.81, P = 1.2 × 10^-9). Both the Dublin-Boston score and the 4-day change in IL-6:IL-10 significantly outperformed IL-6 alone in predicting clinical outcome at day 7.

INTERPRETATION

The Dublin-Boston score is easily calculated and can be applied to a spectrum of hospitalized COVID-19 patients. More informed prognosis could help determine when to escalate care, institute or remove mechanical ventilation, or drive considerations for therapies.

FUNDING

Funding was received from the Elaine Galwey Research Fellowship, American Thoracic Society, National Institutes of Health and the Parker B Francis Research Opportunity Award.

Machine Learning and Prediction of All-Cause Mortality in COPD

BACKGROUND

COPD is a leading cause of mortality.

RESEARCH QUESTION

We hypothesized that applying machine learning to clinical and quantitative CT imaging features would improve mortality prediction in COPD.

STUDY DESIGN AND METHODS

We selected 30 clinical, spirometric, and imaging features as inputs for a random survival forest. We used top features in a Cox regression to create a machine learning mortality prediction (MLMP) in COPD model and also assessed the performance of other statistical and machine learning models. We trained the models in subjects with moderate to severe COPD from a subset of subjects in Genetic Epidemiology of COPD (COPDGene) and tested prediction performance in the remainder of individuals with moderate to severe COPD in COPDGene and Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE). We compared our model with the BMI, airflow obstruction, dyspnea, exercise capacity (BODE) index; BODE modifications; and the age, dyspnea, and airflow obstruction index.

RESULTS

We included 2,632 participants from COPDGene and 1,268 participants from ECLIPSE. The top predictors of mortality were 6-min walk distance, FEV1 % predicted, and age. The top imaging predictor was pulmonary artery-to-aorta ratio. The MLMP-COPD model resulted in a C index ≥ 0.7 in both COPDGene and ECLIPSE (6.4- and 7.2-year median follow-ups, respectively), significantly better than all tested mortality indexes (P < .05). The MLMP-COPD model had fewer predictors but similar performance to that of other models. The group with the highest BODE scores (7-10) had 64% mortality, whereas the highest mortality group defined by the MLMP-COPD model had 77% mortality (P = .012).

INTERPRETATION

An MLMP-COPD model outperformed four existing models for predicting all-cause mortality across two COPD cohorts. Performance of machine learning was similar to that of traditional statistical methods. The model is available online at: https://cdnm.shinyapps.io/cgmortalityapp/.

Statistical considerations for the analysis of massively parallel reporter assays data

Noncoding DNA contains gene regulatory elements that alter gene expression, and the function of these elements can be modified by genetic variation. Massively parallel reporter assays (MPRA) enable high-throughput identification and characterization of functional genetic variants, but the statistical methods to identify allelic effects in MPRA data have not been fully developed. In this study, we demonstrate how the baseline allelic imbalance in MPRA libraries can produce biased results, and we propose a novel, nonparametric, adaptive testing method that is robust to this bias. We compare the performance of this method with other commonly used methods, and we demonstrate that our novel adaptive method controls Type I error in a wide range of scenarios while maintaining excellent power. We have implemented these tests along with routines for simulating MPRA data in the Analysis Toolset for MPRA (@MPRA), an R package for the design and analyses of MPRA experiments. It is publicly available at http://github.com/redaq/atMPRA.

FAM13A Represses AMPK Activity and Regulates Hepatic Glucose and Lipid Metabolism

Obesity commonly co-exists with fatty liver disease with increasing health burden worldwide. Family with Sequence Similarity 13, Member A (FAM13A) has been associated with lipid levels and fat mass by genome-wide association studies (GWAS). However, the function of FAM13A in maintaining metabolic homeostasis in vivo remains unclear. Here, we demonstrated that rs2276936 in this locus has allelic-enhancer activity in massively parallel reporter assays (MPRA) and reporter assay. The DNA region containing rs2276936 regulates expression of endogenous FAM13A in HepG2 cells. In vivo, Fam13a-/- mice are protected from high-fat diet (HFD)-induced fatty liver accompanied by increased insulin sensitivity and reduced glucose production in liver. Mechanistically, loss of Fam13a led to the activation of AMP-activated protein kinase (AMPK) and increased mitochondrial respiration in primary hepatocytes. These findings demonstrate that FAM13A mediates obesity-related dysregulation of lipid and glucose homeostasis. Targeting FAM13A might be a promising treatment of obesity and fatty liver disease.

Association of obesity with bone mineral density and osteoporosis in adults: a systematic review and meta-analysis

OBJECTIVES

Results from previous studies investigating the association of obesity with bone mineral density (BMD) and osteoporosis (OP) are inconsistent. The purpose of the present study was to examine whether obesity is associated with BMD and the risk of OP.

STUDY DESIGN

This is a meta-analysis.

METHODS

A computerised literature search was conducted using the PubMed, Web of Science, China National Knowledge Internet (CNKI) (Chinese) and Wanfang (Chinese) databases for relevant articles published in English or Chinese before the end of December 2017. Means with standard deviations and odds ratios (ORs) with corresponding 95% confidence intervals (95% CIs) were computed to estimate the association between obesity and BMD and the risk of OP by using a random-effects model. In addition, a heterogeneity test and sensitivity analysis were conducted.

RESULTS

Adults with obesity had higher BMD than healthy-weight subjects in both the lumbar spine and femoral neck measurement sites. Obesity was positively related to BMD in the two sites, with merged weighted mean difference (WMD) of 0.07 g/cm² in the lumbar spine (95% CI: 0.045, 0.095; P < 0.001; P_{heterogeneity} <0.001; I² = 89.1%) and WMD of 0.087 g/cm² in the femoral neck (95% CI: 0.063, 0.112; P < 0.001; P_{heterogeneity} <0.001; I² = 92.8%). General obesity was negatively related to femoral neck OP, indicating that obesity is a protective factor for OP, with a merged OR of 0.169 (95% CI: 0.132, 0.216; P < 0.001; P_{heterogeneity} = 0.716; I² = 0.0%).

CONCLUSIONS

Adults with obesity had significantly higher BMD than healthy-weight adults. Obesity was positively associated with BMD and negatively correlated with OP.

Whole Genome Sequencing Identifies CRISPLD2 as a Lung Function Gene in Children With Asthma

BACKGROUND

Asthma is a common respiratory disorder with a highly heterogeneous nature that remains poorly understood. The objective was to use whole genome sequencing (WGS) data to identify regions of common genetic variation contributing to lung function in individuals with a diagnosis of asthma.

METHODS

WGS data were generated for 1,053 individuals from trios and extended pedigrees participating in the family-based Genetic Epidemiology of Asthma in Costa Rica study. Asthma affection status was defined through a physician's diagnosis of asthma, and most participants with asthma also had airway hyperresponsiveness (AHR) to methacholine. Family-based association tests for single variants were performed to assess the associations with lung function phenotypes.

RESULTS

A genome-wide significant association was identified between baseline FEV1/FVC ratio and a single-nucleotide polymorphism in the top hit cysteine-rich secretory protein LCCL domain-containing 2 (CRISPLD2) (rs12051168; P = 3.6 × 10-8 in the unadjusted model) that retained suggestive significance in the covariate-adjusted model (P = 5.6 × 10-6). Rs12051168 was also nominally associated with other related phenotypes: baseline FEV1 (P = 3.3 × 10-3), postbronchodilator (PB) FEV1 (7.3 × 10-3), and PB FEV1/FVC ratio (P = 2.7 × 10-3). The identified baseline FEV1/FVC ratio and rs12051168 association was meta-analyzed and replicated in three independent cohorts in which most participants with asthma also had confirmed AHR (combined weighted z-score P = .015) but not in cohorts without information about AHR.

CONCLUSIONS

These findings suggest that using specific asthma characteristics, such as AHR, can help identify more genetically homogeneous asthma subgroups with genotype-phenotype associations that may not be observed in all children with asthma. CRISPLD2 also may be important for baseline lung function in individuals with asthma who also may have AHR.

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.

Lipid profiles as potential mediators linking body mass index to osteoporosis among Chinese adults: the Henan Rural Cohort Study

This study is to examine the relationship between body mass index (BMI) and osteoporosis in rural areas of China, and then explore whether this association was mediated by lipid profiles. Results showed that with the increasing of BMI, the risk of osteoporosis was reduced, and lipid profiles mediate this association.

INTRODUCTION

To examine the relationship between BMI and osteoporosis, and explore whether this association was mediated by lipid profiles.

METHODS

A total of 8272 participants (18-79 years) were enrolled from the Henan Rural Cohort Study. The bone mineral density of the calcaneus was measured using an ultrasonic bone density apparatus. Logistic regression and restricted cubic splines were used to evaluate the odds ratio (OR) and 95% confidence intervals (95% CI). Mediation analysis using bootstrap was performed to examine the contribution of lipid profiles to BMI-related osteoporosis.

RESULTS

The crude and age-standardized prevalence of osteoporosis were 15.93% and 11.77%, respectively. The mean BMIs were 24.12 kg/m² for participants with osteoporosis and 25.06 kg/m² for non-osteoporosis participants (P < 0.001). After adjusting for potential confounders, subjects with obesity had a lower OR of osteoporosis (0.493 [95% CI: 0.405-0.600], P_trend < 0.001) compared with normal-weight individuals. Mediation analysis showed that lipid profile partly mediated the relationship between BMI and osteoporosis with indirect effect OR (95% CI) of 0.985 (0.978-0.992), and the proportion explained of BMI was 15.48% for lipid profile.

CONCLUSION

With the increasing of BMI, the risk of osteoporosis was reduced in the Chinese adult population, and lipid profiles may be a potential mediator linking reduced risk of osteoporosis. Elucidating the underlying mechanisms will facilitate developing feasible preventive and therapeutic measures for osteoporosis. Chinese clinical trial register: ChiCTR-OOC-15006699.

metaFARVAT: An Efficient Tool for Meta-Analysis of Family-Based, Case-Control, and Population-Based Rare Variant Association Studies

Family-based designs have been shown to be powerful in detecting the significant rare variants associated with human diseases. However, very few significant results have been found owing to relatively small sample sizes and the fact that statistical analyses often suffer from high false-negative error rates. These limitations can be avoided by combining results from multiple studies via meta-analysis. However, statistical methods for meta-analysis with rare variants are limited for family-based samples. In this report, we propose a tool for the meta-analysis of family-based rare variant associations, metaFARVAT. metaFARVAT is based on a quasi-likelihood score for each variant. These scores are combined to generate burden test, variable-threshold test, sequence kernel association test (SKAT), and optimal SKAT statistics. The proposed method tests homogeneous and heterogeneous effects of variants among different studies and can be applied to both quantitative and dichotomous phenotypes. Simulation results demonstrated the robustness and efficiency of the proposed method in different scenarios. By applying metaFARVAT to data from a family-based study and a case-control study, we identified a few promising candidate genes, including DLEC1, which is associated with chronic obstructive pulmonary disease.

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.

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.

Genetic landscape of chronic obstructive pulmonary disease identifies heterogeneous cell-type and phenotype associations

Chronic obstructive pulmonary disease (COPD) is the leading cause of respiratory mortality worldwide. Genetic risk loci provide new insights into disease pathogenesis. We performed a genome-wide association study in 35,735 cases and 222,076 controls from the UK Biobank and additional studies from the International COPD Genetics Consortium. We identified 82 loci associated with P < 5 × 10-8; 47 of these were previously described in association with either COPD or population-based measures of lung function. Of the remaining 35 new loci, 13 were associated with lung function in 79,055 individuals from the SpiroMeta consortium. Using gene expression and regulation data, we identified functional enrichment of COPD risk loci in lung tissue, smooth muscle, and several lung cell types. We found 14 COPD loci shared with either asthma or pulmonary fibrosis. COPD genetic risk loci clustered into groups based on associations with quantitative imaging features and comorbidities. Our analyses provide further support for the genetic susceptibility and heterogeneity of COPD.

Identification of Functional Variants in the FAM13A Chronic Obstructive Pulmonary Disease Genome-Wide Association Study Locus by Massively Parallel Reporter Assays

RATIONALE

The identification of causal variants responsible for disease associations from genome-wide association studies (GWASs) facilitates functional understanding of the biological mechanisms by which those genetic variants influence disease susceptibility.

OBJECTIVE

We aim to identify causal variants in or near the FAM13A (family with sequence similarity member 13A) GWAS locus associated with chronic obstructive pulmonary disease (COPD).

METHODS

We used an integrated approach featuring conditional genetic analysis, massively parallel reporter assays (MPRAs), traditional reporter assays, chromatin conformation capture assays, and clustered regularly interspaced short palindromic repeats (CRISPR)-based gene editing to characterize COPD-associated regulatory variants in the FAM13A region in human bronchial epithelial cell lines.

MEASUREMENTS AND MAIN RESULTS

Conditional genetic association suggests the presence of two independent COPD association signals in FAM13A. MPRAs identified 45 regulatory variants within FAM13A, among which six variants were prioritized for further investigation. Three COPD-associated variants demonstrated significant allele-specific activity in reporter assays. One of three variants, rs2013701, was tested in the endogenous genomic context by CRISPR-based genome editing that confirmed its allele-specific effects on FAM13A expression and on cell proliferation, providing functional characterization for this COPD-associated variant.

CONCLUSIONS

The human GWAS association near FAM13A may contain independent association signals. MPRAs identified multiple functional variants in this region, including rs2013701, a putative COPD-causing variant with allele-specific regulatory activity.

Whole exome sequencing analysis in severe chronic obstructive pulmonary disease

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

Rare Variants in Known Susceptibility Loci and Their Contribution to Risk of Lung Cancer

BACKGROUND

Genome-wide association studies are widely used to map genomic regions contributing to lung cancer (LC) susceptibility, but they typically do not identify the precise disease-causing genes/variants. To unveil the inherited genetic variants that cause LC, we performed focused exome-sequencing analyses on genes located in 121 genome-wide association study-identified loci previously implicated in the risk of LC, chronic obstructive pulmonary disease, pulmonary function level, and smoking behavior.

METHODS

Germline DNA from 260 case patients with LC and 318 controls were sequenced by utilizing VCRome 2.1 exome capture. Filtering was based on enrichment of rare and potential deleterious variants in cases (risk alleles) or controls (protective alleles). Allelic association analyses of single-variant and gene-based burden tests of multiple variants were performed. Promising candidates were tested in two independent validation studies with a total of 1773 case patients and 1123 controls.

RESULTS

We identified 48 rare variants with deleterious effects in the discovery analysis and validated 12 of the 43 candidates that were covered in the validation platforms. The top validated candidates included one well-established truncating variant, namely, BRCA2, DNA repair associated gene (BRCA2) K3326X (OR = 2.36, 95% confidence interval [CI]: 1.38-3.99), and three newly identified variations, namely, lymphotoxin beta gene (LTB) p.Leu87Phe (OR = 7.52, 95% CI: 1.01-16.56), prolyl 3-hydroxylase 2 gene (P3H2) p.Gln185His (OR = 5.39, 95% CI: 0.75-15.43), and dishevelled associated activator of morphogenesis 2 gene (DAAM2) p.Asp762Gly (OR = 0.25, 95% CI: 0.10-0.79). Burden tests revealed strong associations between zinc finger protein 93 gene (ZNF93), DAAM2, bromodomain containing 9 gene (BRD9), and the gene LTB and LC susceptibility.

CONCLUSION

Our results extend the catalogue of regions associated with LC and highlight the importance of germline rare coding variants in LC susceptibility.

WISARD: workbench for integrated superfast association studies for related datasets

BACKGROUND

A Mendelian transmission produces phenotypic and genetic relatedness between family members, giving family-based analytical methods an important role in genetic epidemiological studies-from heritability estimations to genetic association analyses. With the advance in genotyping technologies, whole-genome sequence data can be utilized for genetic epidemiological studies, and family-based samples may become more useful for detecting de novo mutations. However, genetic analyses employing family-based samples usually suffer from the complexity of the computational/statistical algorithms, and certain types of family designs, such as incorporating data from extended families, have rarely been used.

RESULTS

We present a Workbench for Integrated Superfast Association studies for Related Data (WISARD) programmed in C/C++. WISARD enables the fast and a comprehensive analysis of SNP-chip and next-generation sequencing data on extended families, with applications from designing genetic studies to summarizing analysis results. In addition, WISARD can automatically be run in a fully multithreaded manner, and the integration of R software for visualization makes it more accessible to non-experts.

CONCLUSIONS

Comparison with existing toolsets showed that WISARD is computationally suitable for integrated analysis of related subjects, and demonstrated that WISARD outperforms existing toolsets. WISARD has also been successfully utilized to analyze the large-scale massive sequencing dataset of chronic obstructive pulmonary disease data (COPD), and we identified multiple genes associated with COPD, which demonstrates its practical value.

A Genome-Wide Linkage Study for Chronic Obstructive Pulmonary Disease in a Dutch Genetic Isolate Identifies Novel Rare Candidate Variants

Chronic obstructive pulmonary disease (COPD) is a complex and heritable disease, associated with multiple genetic variants. Specific familial types of COPD may be explained by rare variants, which have not been widely studied. We aimed to discover rare genetic variants underlying COPD through a genome-wide linkage scan. Affected-only analysis was performed using the 6K Illumina Linkage IV Panel in 142 cases clustered in 27 families from a genetic isolate, the Erasmus Rucphen Family (ERF) study. Potential causal variants were identified by searching for shared rare variants in the exome-sequence data of the affected members of the families contributing most to the linkage peak. The identified rare variants were then tested for association with COPD in a large meta-analysis of several cohorts. Significant evidence for linkage was observed on chromosomes 15q14-15q25 [logarithm of the odds (LOD) score = 5.52], 11p15.4-11q14.1 (LOD = 3.71) and 5q14.3-5q33.2 (LOD = 3.49). In the chromosome 15 peak, that harbors the known COPD locus for nicotinic receptors, and in the chromosome 5 peak we could not identify shared variants. In the chromosome 11 locus, we identified four rare (minor allele frequency (MAF) <0.02), predicted pathogenic, missense variants. These were shared among the affected family members. The identified variants localize to genes including neuroblast differentiation-associated protein (AHNAK), previously associated with blood biomarkers in COPD, phospholipase C Beta 3 (PLCB3), shown to increase airway hyper-responsiveness, solute carrier family 22-A11 (SLC22A11), involved in amino acid metabolism and ion transport, and metallothionein-like protein 5 (MTL5), involved in nicotinate and nicotinamide metabolism. Association of SLC22A11 and MTL5 variants were confirmed in the meta-analysis of 9,888 cases and 27,060 controls. In conclusion, we have identified novel rare variants in plausible genes related to COPD. Further studies utilizing large sample whole-genome sequencing should further confirm the associations at chromosome 11 and investigate the chromosome 15 and 5 linked regions.

Selecting cases and controls for DNA sequencing studies using family histories of disease

Recent improvements in sequencing technology have enabled the investigation of so-called missing heritability, and a large number of affected subjects have been sequenced in order to detect significant associations between human diseases and rare variants. However, the cost of genome sequencing is still high, and a statistically powerful strategy for selecting informative subjects would be useful. Therefore, in this report, we propose a new statistical method for selecting cases and controls for sequencing studies based on family history. We assume that disease status is determined by unobserved liability scores. Our method consists of two steps: first, the conditional means of liability are estimated with the liability threshold model given the individual's disease status and those of their relatives. Second, the informative subjects are selected with the estimated conditional means. Our simulation studies showed that statistical power is substantially affected by the subject selection strategy chosen, and power is maximized when affected (unaffected) subjects with high (low) risks are selected as cases (controls). The proposed method was successfully applied to genome-wide association studies for type 2 diabetes, and our analysis results reveal the practical value of the proposed methods. Copyright © 2017 John Wiley & Sons, Ltd.

Exome Sequencing Analysis in Severe, Early-Onset Chronic Obstructive Pulmonary Disease

RATIONALE

Genomic regions identified by genome-wide association studies explain only a small fraction of heritability for chronic obstructive pulmonary disease (COPD). Alpha-1 antitrypsin deficiency shows that rare coding variants of large effect also influence COPD susceptibility. We hypothesized that exome sequencing in families identified through a proband with severe, early-onset COPD would identify additional rare genetic determinants of large effect.

OBJECTIVES

To identify rare genetic determinants of severe COPD.

METHODS

We applied filtering approaches to identify potential causal variants for COPD in whole exomes from 347 subjects in 49 extended pedigrees from the Boston Early-Onset COPD Study. We assessed the power of this approach under different levels of genetic heterogeneity using simulations. We tested genes identified in these families using gene-based association tests in exomes of 204 cases with severe COPD and 195 resistant smokers from the COPDGene study. In addition, we examined previously described loci associated with COPD using these datasets.

MEASUREMENTS AND MAIN RESULTS

We identified 69 genes with predicted deleterious nonsynonymous, stop, or splice variants that segregated with severe COPD in at least two pedigrees. Four genes (DNAH8, ALCAM, RARS, and GBF1) also demonstrated an increase in rare nonsynonymous, stop, and/or splice mutations in cases compared with resistant smokers from the COPDGene study; however, these results were not statistically significant. We demonstrate the limitations of the power of this approach under genetic heterogeneity through simulation.

CONCLUSIONS

Rare deleterious coding variants may increase risk for COPD, but multiple genes likely contribute to COPD susceptibility.

Genetic loci associated with chronic obstructive pulmonary disease overlap with loci for lung function and pulmonary fibrosis

Chronic obstructive pulmonary disease (COPD) is a leading cause of mortality worldwide. We performed a genetic association study in 15,256 cases and 47,936 controls, with replication of select top results (P < 5 × 10^-6) in 9,498 cases and 9,748 controls. In the combined meta-analysis, we identified 22 loci associated at genome-wide significance, including 13 new associations with COPD. Nine of these 13 loci have been associated with lung function in general population samples, while 4 (EEFSEC, DSP, MTCL1, and SFTPD) are new. We noted two loci shared with pulmonary fibrosis (FAM13A and DSP) but that had opposite risk alleles for COPD. None of our loci overlapped with genome-wide associations for asthma, although one locus has been implicated in joint susceptibility to asthma and obesity. We also identified genetic correlation between COPD and asthma. Our findings highlight new loci associated with COPD, demonstrate the importance of specific loci associated with lung function to COPD, and identify potential regions of genetic overlap between COPD and other respiratory diseases.

Gene-based segregation method for identifying rare variants in family-based sequencing studies

Whole-exome sequencing using family data has identified rare coding variants in Mendelian diseases or complex diseases with Mendelian subtypes, using filters based on variant novelty, functionality, and segregation with the phenotype within families. However, formal statistical approaches are limited. We propose a gene-based segregation test (GESE) that quantifies the uncertainty of the filtering approach. It is constructed using the probability of segregation events under the null hypothesis of Mendelian transmission. This test takes into account different degrees of relatedness in families, the number of functional rare variants in the gene, and their minor allele frequencies in the corresponding population. In addition, a weighted version of this test allows incorporating additional subject phenotypes to improve statistical power. We show via simulations that the GESE and weighted GESE tests maintain appropriate type I error rate, and have greater power than several commonly used region-based methods. We apply our method to whole-exome sequencing data from 49 extended pedigrees with severe, early-onset chronic obstructive pulmonary disease (COPD) in the Boston Early-Onset COPD study (BEOCOPD) and identify several promising candidate genes. Our proposed methods show great potential for identifying rare coding variants of large effect and high penetrance for family-based sequencing data. The proposed tests are implemented in an R package that is available on CRAN (https://cran.r-project.org/web/packages/GESE/).