Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

SharePro: an accurate and efficient genetic colocalization method accounting for multiple causal signals.

MOTIVATION: Colocalization analysis is commonly used to assess whether two or more traits share the same genetic signals identified in genome-wide association studies (GWAS), and is important for prioritizing targets for functional follow-up of GWAS results. Existing colocalization methods can have suboptimal performance when there are multiple causal variants in one genomic locus. RESULTS: We propose SharePro to extend the COLOC framework for colocalization analysis. SharePro integrates linkage disequilibrium (LD) modeling and colocalization assessment by grouping correlated variants into effect groups. With an efficient variational inference algorithm, posterior colocalization probabilities can be accurately estimated. In simulation studies, SharePro demonstrated increased power with a well-controlled false positive rate at a low computational cost. Compared to existing methods, SharePro provided stronger and more consistent colocalization evidence for known lipid-lowering drug target proteins and their corresponding lipid traits. Through an additional challenging case of the colocalization analysis of the circulating abundance of R-spondin 3 GWAS and estimated bone mineral density GWAS, we demonstrated the utility of SharePro in identifying biologically plausible colocalized signals. AVAILABILITY AND IMPLEMENTATION: SharePro for colocalization analysis is written in Python and openly available at https://github.com/zhwm/SharePro_coloc.

Genetically Predicted Body Mass Index and Mortality in COPD.

RATIONALE: BMI is associated with COPD mortality, but the underlying mechanisms are unclear. The effect of genetic variants aggregated into a polygenic score may elucidate causal mechanisms and predict risk. OBJECTIVES: To examine the associations of genetically predicted BMI with all-cause and cause-specific mortality in COPD. METHODS: We developed a polygenic score for BMI (PGSBMI) and tested for associations of the PGSBMI with all-cause, respiratory, and cardiovascular mortality in participants with COPD from the COPDGene, ECLIPSE, and Framingham Heart studies. We calculated the difference between measured BMI and PGS-predicted BMI (BMIdiff) and categorized participants into groups of discordantly low (BMIdiff < 20th percentile), concordant (BMIdiff between 20th - 80th percentile), and discordantly high (BMIdiff > 80th percentile) BMI. We applied Cox models, examined potential non-linear associations of the PGSBMI and BMIdiff with mortality, and summarized results with meta-analysis. MEASUREMENTS AND MAIN RESULTS: We observed significant non-linear associations of measured BMI and BMIdiff, but not PGSBMI, with all-cause mortality. In meta-analyses, a one standard deviation increase in the PGSBMI was associated with an increased hazard for cardiovascular mortality (HR=1.29, 95% CI=1.12-1.49), but not with respiratory or all-cause mortality. Compared to participants with concordant measured and genetically predicted BMI, those with discordantly low BMI had higher mortality risk for all-cause (HR=1.57, CI=1.41-1.74) and respiratory death (HR=2.01, CI=1.61-2.51). CONCLUSIONS: In people with COPD, higher genetically predicted BMI is associated with higher cardiovascular mortality but not respiratory mortality. Individuals with discordantly low BMI have higher all-cause and respiratory mortality compared to those with concordant BMI.

Genome-wide association study of preserved ratio impaired spirometry (PRISm).

BACKGROUND: Preserved ratio impaired spirometry (PRISm) is defined as a forced expiratory volume in 1 s (FEV1) <80% predicted and FEV1/forced vital capacity ≥0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities. METHODS: We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed. RESULTS: 22 signals reached significance in the joint meta-analysis, including four signals novel for lung function. A strong genome-wide genetic correlation (rg) between PRISm and spirometric COPD (rg=0.62, p<0.001) was observed, and genetic correlation with type 2 diabetes (rg=0.12, p=0.007). Phenome-wide association studies showed that 18 of 22 signals were associated with diabetic traits and seven with blood pressure traits. CONCLUSION: This is the first GWAS to successfully identify SNPs associated with PRISm. Four of the signals, rs7652391 (nearest gene MECOM), rs9431040 (HLX), rs62018863 (TMEM114) and rs185937162 (HLA-B), have not been described in association with lung function before, demonstrating the utility of using different lung function phenotypes in GWAS. Genetic factors associated with PRISm are strongly correlated with risk of both other lung diseases and extrapulmonary comorbidity.

Association of Dynamic Trajectories of Time-Series Data and Life-Threatening Mass Effect in Large Middle Cerebral Artery Stroke.

BACKGROUND: Life-threatening, space-occupying mass effect due to cerebral edema and/or hemorrhagic transformation is an early complication of patients with middle cerebral artery stroke. Little is known about longitudinal trajectories of laboratory and vital signs leading up to radiographic and clinical deterioration related to this mass effect. METHODS: We curated a retrospective data set of 635 patients with large middle cerebral artery stroke totaling 95,463 data points for 10 longitudinal covariates and 40 time-independent covariates. We assessed trajectories of the 10 longitudinal variables during the 72 h preceding three outcomes representative of life-threatening mass effect: midline shift ≥ 5 mm, pineal gland shift (PGS) > 4 mm, and decompressive hemicraniectomy (DHC). We used a "backward-looking" trajectory approach. Patients were aligned based on outcome occurrence time and the trajectory of each variable was assessed before that outcome by accounting for cases and noncases, adjusting for confounders. We evaluated longitudinal trajectories with Cox proportional time-dependent regression. RESULTS: Of 635 patients, 49.0% were female, and the mean age was 69 years. Thirty five percent of patients had midline shift ≥ 5 mm, 24.3% of patients had PGS > 4 mm, and 10.7% of patients underwent DHC. Backward-looking trajectories showed mild increases in white blood cell count (10-11 K/UL within 72 h), temperature (up to half a degree within 24 h), and sodium levels (1-3 mEq/L within 24 h) before the three outcomes of interest. We also observed a decrease in heart rate (75-65 beats per minute) 24 h before DHC. We found a significant association between increased white blood cell count with PGS > 4 mm (hazard ratio 1.05, p value 0.007). CONCLUSIONS: Longitudinal profiling adjusted for confounders demonstrated that white blood cell count, temperature, and sodium levels appear to increase before radiographic and clinical indicators of space-occupying mass effect. These findings will inform the development of multivariable dynamic risk models to aid prediction of life-threatening, space-occupying mass effect.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data.

INTRODUCTION: Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci. METHODS: We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants. RESULTS: Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses. DISCUSSION: This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function.

BACKGROUND: Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function. OBJECTIVES: We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions. METHODS: We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings. RESULTS: Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers. CONCLUSION: Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Association of large core middle cerebral artery stroke and hemorrhagic transformation with hospitalization outcomes.

Historically, investigators have not differentiated between patients with and without hemorrhagic transformation (HT) in large core ischemic stroke at risk for life-threatening mass effect (LTME) from cerebral edema. Our objective was to determine whether LTME occurs faster in those with HT compared to those without. We conducted a two-center retrospective study of patients with ≥ 1/2 MCA territory infarct between 2006 and 2021. We tested the association of time-to-LTME and HT subtype (parenchymal, petechial) using Cox regression, controlling for age, mean arterial pressure, glucose, tissue plasminogen activator, mechanical thrombectomy, National Institute of Health Stroke Scale, antiplatelets, anticoagulation, temperature, and stroke side. Secondary and exploratory outcomes included mass effect-related death, all-cause death, disposition, and decompressive hemicraniectomy. Of 840 patients, 358 (42.6%) had no HT, 403 (48.0%) patients had petechial HT, and 79 (9.4%) patients had parenchymal HT. LTME occurred in 317 (37.7%) and 100 (11.9%) had mass effect-related deaths. Parenchymal (HR 8.24, 95% CI 5.46-12.42, p < 0.01) and petechial HT (HR 2.47, 95% CI 1.92-3.17, p < 0.01) were significantly associated with time-to-LTME and mass effect-related death. Understanding different risk factors and sequelae of mass effect with and without HT is critical for informed clinical decisions.

Later Midline Shift Is Associated with Better Outcomes after Large Middle Cerebral Artery Stroke.

Background/Objective: Space occupying cerebral edema is the most feared early complication after large ischemic stroke, occurring in up to 30% of patients with middle cerebral artery (MCA) occlusion, and is reported to peak 2-4 days after injury. Little is known about the factors and outcomes associated with peak edema timing, especially when it occurs after 96 hours. We aimed to characterize differences between patients who experienced maximum midline shift (MLS) or decompressive hemicraniectomy (DHC) in the acute (<48 hours), average (48-96 hours), and subacute (>96 hours) groups and determine whether patients with subacute peak edema timing have improved discharge dispositions. Methods: We performed a two-center, retrospective study of patients with ≥1/2 MCA territory infarct and MLS. We constructed a multivariable model to test the association of subacute peak edema and favorable discharge disposition, adjusting for age, admission Alberta Stroke Program Early CT Score (ASPECTS), National Institute of Health Stroke Scale (NIHSS), acute thrombolytic intervention, cerebral atrophy, maximum MLS, parenchymal hemorrhagic transformation, DHC, and osmotic therapy receipt. Results: Of 321 eligible patients with MLS, 32%, 36%, and 32% experienced acute, average, and subacute peak edema. Subacute peak edema was significantly associated with higher odds of favorable discharge than non-subacute swelling, adjusting for confounders (aOR, 1.85; 95% CI, 1.05-3.31). Conclusions: Subacute peak edema after large MCA stroke is associated with better discharge disposition compared to earlier peak edema courses. Understanding how the timing of cerebral edema affects risk of unfavorable discharge has important implications for treatment decisions and prognostication.

Follow-up ASPECTS improves prediction of potentially lethal malignant edema in patients with large middle cerebral artery stroke.

BACKGROUND: Recent studies have shown that follow-up head CT is a strong predictor of functional outcomes in patients with middle cerebral artery stroke and mechanical thrombectomy. We sought to determine whether total and/or regional follow-up Alberta Stroke Program Early CT Score (ASPECTSfu) are associated with important clinical outcomes during hospitalization and improve the performance of clinical prediction models of potentially lethal malignant edema (PLME). METHODS: We conducted a retrospective study of patients at three medical centers in a major North American metropolitan area with baseline and follow-up head CTs after large middle cerebral artery stroke between 2006 and 2022. We used multivariable logistic regression to test the association of total and regional ASPECTSfu with PLME (cerebral edema related death or surgery), adjusting for total baseline ASPECTS, age, sex, admission glucose, tissue plasminogen activator, and mechanical thrombectomy. We compared existing clinical risk models with and without total or regional ASPECTSfu using area under the curve. RESULTS: In our 560 patient cohort, lower total ASPECTSfu was significantly associated with higher odds of PLME when adjusting for confounders (OR 1.69, 95% CI 1.49 to 2.0), and improved model discrimination compared with existing models and models using baseline ASPECTS. Deep territory involvement (OR 2.46, 95% CI 1.53 to 4.01) and anterior territory involvement (OR 3.23, 95% CI 1.88 to 5.71) were significantly associated with PLME. CONCLUSIONS: Lower ASPECTSfu and certain locations on regional ASPECTSfu, including deep and anterior areas, were significantly associated with PLME. Including ASPECTSfu information improved discrimination of established edema prediction models and could be used immediately to help facilitate clinical management decisions and prognostication.

Bik promotes proteasomal degradation to control low-grade inflammation.

Although chronic low-grade inflammation does not cause immediate clinical symptoms, over the longer term, it can enhance other insults or age-dependent damage to organ systems and thereby contribute to age-related disorders, such as respiratory disorders, heart disease, metabolic disorders, autoimmunity, and cancer. However, the molecular mechanisms governing low-level inflammation are largely unknown. We discovered that Bcl-2-interacting killer (Bik) deficiency causes low-level inflammation even at baseline and the development of spontaneous emphysema in female but not male mice. Similarly, a single nucleotide polymorphism that reduced Bik levels was associated with increased inflammation and enhanced decline in lung function in humans. Transgenic expression of Bik in the airways of Bik-deficient mice inhibited allergen- or LPS-induced lung inflammation and reversed emphysema in female mice. Bik deficiency increased nuclear but not cytosolic p65 levels because Bik, by modifying the BH4 domain of Bcl-2, interacted with regulatory particle non-ATPase 1 (RPN1) and RPN2 and enhanced proteasomal degradation of nuclear proteins. Bik deficiency increased inflammation primarily in females because Bcl-2 and Bik levels were reduced in lung tissues and airway cells of female compared with male mice. Therefore, controlling low-grade inflammation by modifying the unappreciated role of Bik and Bcl-2 in facilitating proteasomal degradation of nuclear proteins may be crucial in treating chronic age-related diseases.