Genome-Wide Association Study Identifies Pharmacogenomic Variants Associated With Metformin Glycemic Response in African American Patients With Type 2 Diabetes

OBJECTIVE

Metformin is the most common treatment for type 2 diabetes (T2D). However, there have been no pharmacogenomic studies for T2D in which a population of color was used in the discovery analysis. This study sought to identify genomic variants associated with metformin response in African American patients with diabetes.

RESEARCH DESIGN AND METHODS

Patients in the discovery set were adult, African American participants from the Diabetes Multi-omic Investigation of Drug Response (DIAMOND), a cohort study of patients with T2D from a health system serving southeast Michigan. DIAMOND participants had genome-wide genotype data and longitudinal electronic records of laboratory results and medication fills. The genome-wide discovery analysis identified polymorphisms correlated to changes in glycated hemoglobin (HbA1c) levels among individuals on metformin monotherapy. Lead associations were assessed for replication in an independent cohort of African American participants from Kaiser Permanente Northern California (KPNC) and in European American participants from DIAMOND.

RESULTS

The discovery set consisted of 447 African American participants, whereas the replication sets included 353 African American KPNC participants and 466 European American DIAMOND participants. The primary analysis identified a variant, rs143276236, in the gene ARFGEF3, which met the threshold for genome-wide significance, replicated in KPNC African Americans, and was still significant in the meta-analysis (P = 1.17 × 10-9). None of the significant discovery variants replicated in European Americans DIAMOND participants.

CONCLUSIONS

We identified a novel and biologically plausible genetic variant associated with a change in HbA1c levels among African American patients on metformin monotherapy. These results highlight the importance of diversity in pharmacogenomic studies.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Evaluation of a New Aptamer-Based Array for Soluble Suppressor of Tumorgenicity (ST2) and N-terminal Pro-B-Type Natriuretic Peptide (NTproBNP) in Heart Failure Patients

BACKGROUND

Recent advances in multi-marker platforms offer faster data generation, but the fidelity of these methods compared to the ELISA is not established. We tested the correlation and predictive performance of SOMAscan vs. ELISA methods for NTproBNP and ST2.

METHODS

Patients ≥ 18 years with heart failure and ejection fraction < 50% were enrolled. We tested the correlation between SOMA and ELISA for each biomarker and their association with outcomes.

RESULTS

There was good correlation of SOMA vs. ELISA for ST2 (ρ = 0.71) and excellent correlation for NTproBNP (ρ = 0.94). The two versions of both markers were not significantly different regarding survival association. The two ST2 assays and NTproBNP assays were similarly associated with all-cause mortality and cardiovascular mortality. These associations remained statistically significant when adjusted for MAGGIC risk score (all p < 0.05).

CONCLUSION

SOMAscan quantifications of ST2 and NTproBNP correlate to ELISA versions and carry similar prognosis.

Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction

In patients with heart failure with reduced ejection fraction (HFrEF), individual responses to beta-blockers vary. Candidate gene pharmacogenetic studies yielded significant but inconsistent results, and they may have missed important associations. Our objective was to use an unbiased genome-wide association study (GWAS) to identify loci influencing beta-blocker survival benefit in HFrEF patients. Genetic variant × beta-blocker exposure interactions were tested in Cox proportional hazards models for all-cause mortality stratified by self-identified race. The models were adjusted for clinical risk factors and propensity scores. A prospective HFrEF registry (469 black and 459 white patients) was used for discovery, and linkage disequilibrium (LD) clumped variants with a beta-blocker interaction of p < 5 × 10-5, were tested for Bonferroni-corrected validation in a multicenter HFrEF clinical trial (288 black and 579 white patients). A total of 229 and 18 variants in black and white HFrEF patients, respectively, had interactions with beta-blocker exposure at p < 5 × 10-5 upon discovery. After LD-clumping, 100 variants and 4 variants in the black and white patients, respectively, remained for validation but none reached statistical significance. In conclusion, genetic variants of potential interest were identified in a discovery-based GWAS of beta-blocker survival benefit in HFrEF patients, but none were validated in an independent dataset. Larger cohorts or alternative approaches, such as polygenic scores, are needed.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Common Variants on FGD5 Increase Hazard of Mortality or Rehospitalization in Patients With Heart Failure From the ASCEND-HF Trial

BACKGROUND

Heart failure remains a global health burden, and patients hospitalized are particularly at risk, but genetic associates for subsequent death or rehospitalization are still lacking.

METHODS

The genetic substudy of the ASCEND-HF trial (Acute Study of Clinical Effectiveness of Nesiritide in Decompensated Heart Failure) was used to perform genome-wide association study and transethnic meta-analysis. The overall trial included the patients of self-reported European ancestry (n=2173) and African ancestry (n=507). The end point was death or heart failure rehospitalization within 180 days. Cox models adjusted for 11 a priori predictors of rehospitalization and 5 genetic principal components were used to test the association between single-nucleotide polymorphisms and outcome. Summary statistics from the 2 populations were combined via meta-analysis with the significance threshold considered P<5×10-8.

RESULTS

Common variants (rs2342882 and rs35850039 in complete linkage disequilibrium) located in FGD5 were significantly associated with the primary outcome in both ancestry groups (European Americans: hazard ratio [HR], 1.38; P=2.42×10-6; African ancestry: HR, 1.51; P=4.43×10-3; HR in meta-analysis, 1.41; P=4.25×10-8). FGD5 encodes a regulator of VEGF (vascular endothelial growth factor)-mediated angiogenesis, and in silico investigation revealed several previous genome-wide association study hits in this gene, among which rs748431 was associated with our outcome (HR, 1.20; meta P<0.01). Sensitivity analysis proved FGD5 common variants survival association did not appear to operate via coronary artery disease or nesiritide treatment (P>0.05); and the signal was still significant when changing the censoring time from 180 to 30 days (HR, 1.39; P=1.59×10-5).

CONCLUSIONS

In this multiethnic genome-wide association study of ASCEND-HF, single-nucleotide polymorphisms in FGD5 were associated with increased risk of death or rehospitalization. Additional investigation is required to examine biological mechanisms and whether FGD5 could be a therapeutic target.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT00475852.

Identification of Genetic Variation Influencing Metformin Response in a Multiancestry Genome-Wide Association Study in the Diabetes Prevention Program (DPP)

Genome-wide significant loci for metformin response in type 2 diabetes reported elsewhere have not been replicated in the Diabetes Prevention Program (DPP). To assess pharmacogenetic interactions in prediabetes, we conducted a genome-wide association study (GWAS) in the DPP. Cox proportional hazards models tested associations with diabetes incidence in the metformin (MET; n = 876) and placebo (PBO; n = 887) arms. Multiple linear regression assessed association with 1-year change in metformin-related quantitative traits, adjusted for baseline trait, age, sex, and 10 ancestry principal components. We tested for gene-by-treatment interaction. No significant associations emerged for diabetes incidence. We identified four genome-wide significant variants after correcting for correlated traits (P < 9 × 10-9). In the MET arm, rs144322333 near ENOSF1 (minor allele frequency [MAF]AFR = 0.07; MAFEUR = 0.002) was associated with an increase in percentage of glycated hemoglobin (per minor allele, β = 0.39 [95% CI 0.28, 0.50]; P = 2.8 × 10-12). rs145591055 near OMSR (MAF = 0.10 in American Indians) was associated with weight loss (kilograms) (per G allele, β = -7.55 [95% CI -9.88, -5.22]; P = 3.2 × 10-10) in the MET arm. Neither variant was significant in PBO; gene-by-treatment interaction was significant for both variants [P(G×T) < 1.0 × 10-4]. Replication in individuals with diabetes did not yield significant findings. A GWAS for metformin response in prediabetes revealed novel ethnic-specific associations that require further investigation but may have implications for tailored therapy.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

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.

Differences in Self-Reported Food Allergy and Food-Associated Anaphylaxis by Race and Ethnicity Among SAPPHIRE Cohort Participants

BACKGROUND

Although food allergies are considered common, relatively little is known about disparities in food allergy by race in the United States.

OBJECTIVE

To evaluate differences in reported food allergy and food-associated anaphylaxis among individuals enrolled in a longitudinal cohort study from metropolitan Detroit, Michigan.

METHODS

Participants in the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity (SAPPHIRE) were asked about food allergies, including the inciting food and associated symptoms. Individuals were considered to have food-associated anaphylaxis if symptoms coincided with established clinical criteria. Logistic regression was used to assess whether race difference persisted after adjusting for and stratifying by potential confounders. African genetic ancestry was individually estimated among African American SAPPHIRE participants to assess whether ancestry was associated with food allergy.

RESULTS

Within the SAPPHIRE cohort, African American participants were significantly more likely to report food allergy (26.1% vs 17%; P = 3.47 × 10-18) and have food-associated anaphylactic symptoms (12.7% vs 7%; P = 4.65 × 10-14) when compared with European American participants. Allergy to seafood accounted for the largest difference (13.1% vs 4.6%; P = 1.38 × 10-31). Differences in food allergy by race persisted after adjusting for potential confounders including asthma status. Among African American participants, the proportion of African ancestry was not associated with any outcome evaluated.

CONCLUSION

Compared with European Americans, African Americans appear to be at higher risk for developing food allergy and food-associated anaphylaxis, particularly with regard to seafood allergy. The lack of association with genetic ancestry suggests that socioenvironmental determinants may play a role in these disparities.

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.

Understanding TBI as a Risk Factor Versus a Means of Suicide Death Using Electronic Health Record Data

OBJECTIVE

The aim of this research was to examine predictors and characterize causes of suicide death in people with traumatic brain injury (TBI) and conduct sensitivity analyses with and without people whose first diagnosis of TBI occurred within 3 days of their suicide death.

METHODS

This case-control study examined suicide risk for people with TBI in eight Mental Health Research Network-affiliated healthcare systems. Sample 1 included 61 persons with TBI who died by suicide and their 75 matched controls with TBI who did not die by suicide between January 1, 2000, and December 31, 2013. Sample 2 excluded the 34 persons with TBI whose first TBI diagnosis occurred within 3 days of their suicide death and their 46 matched controls. Descriptive statistics characterized the sample stratified by cases and controls, while conditional logistic regression models estimated the adjusted odds of suicide.

RESULTS

Over half of suicide deaths occurred within 3 days of a person's first diagnosis of TBI in the larger sample. After excluding these persons, people with TBI were 2.84 (95% confidence interval [CI]: 2.15-2.73) times more likely to die by suicide than were people without TBI. Among those with TBI, men were 16.39 times (95% CI: 1.89-142.15) more likely to die by suicide than were women.

CONCLUSIONS

Accounting for TBI as a potential consequence of suicide attenuates the association between TBI and suicide, but a robust association persists-especially among men. Ultimately, all people with TBI should be carefully screened and monitored for suicide risk.HIGHLIGHTSPeople with traumatic brain injury (TBI) were at considerably elevated risk for suicide deathMen with TBI had significantly increased risk of suicide death compared to women with TBITBI timing suggests confusion of risk factors for and consequences of suicide.

Admixture mapping of severe asthma exacerbations in Hispanic/Latino children and youth

BACKGROUND

In the USA, genetically admixed populations have the highest asthma prevalence and severe asthma exacerbations rates. This could be explained not only by environmental factors but also by genetic variants that exert ethnic-specific effects. However, no admixture mapping has been performed for severe asthma exacerbations.

OBJECTIVE

We sought to identify genetic variants associated with severe asthma exacerbations in Hispanic/Latino subgroups by means of admixture mapping analyses and fine mapping, and to assess their transferability to other populations and potential functional roles.

METHODS

We performed an admixture mapping in 1124 Puerto Rican and 625 Mexican American children with asthma. Fine-mapping of the significant peaks was performed via allelic testing of common and rare variants. We performed replication across Hispanic/Latino subgroups, and the transferability to non-Hispanic/Latino populations was assessed in 1001 African Americans, 1250 Singaporeans and 941 Europeans with asthma. The effects of the variants on gene expression and DNA methylation from whole blood were also evaluated in participants with asthma and in silico with data obtained through public databases.

RESULTS

Genomewide significant associations of Indigenous American ancestry with severe asthma exacerbations were found at 5q32 in Mexican Americans as well as at 13q13-q13.2 and 3p13 in Puerto Ricans. The single nucleotide polymorphism (SNP) rs1144986 (C5orf46) showed consistent effects for severe asthma exacerbations across Hispanic/Latino subgroups, but it was not validated in non-Hispanics/Latinos. This SNP was associated with DPYSL3 DNA methylation and SCGB3A2 gene expression levels.

CONCLUSIONS

Admixture mapping study of asthma exacerbations revealed a novel locus that exhibited Hispanic/Latino-specific effects and regulated DPYSL3 and SCGB3A2.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Genome-wide association study in minority children with asthma implicates DNAH5 in bronchodilator responsiveness

Variability in response to short-acting β2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10-8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P < 0.05). No genetic variant was associated with BDR at the genome-wide significant threshold in Mexicans and African Americans. Our findings help inform genetic underpinnings of BDR for understudied minority patients with asthma, but the limited availability of genetic data for racial/ethnic minority children with asthma remains a paramount challenge.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Integrative omics of schizophrenia: from genetic determinants to clinical classification and risk prediction

Schizophrenia (SCZ) is a debilitating neuropsychiatric disorder with high heritability and complex inheritance. In the past decade, successful identification of numerous susceptibility loci has provided useful insights into the molecular etiology of SCZ. However, applications of these findings to clinical classification and diagnosis, risk prediction, or intervention for SCZ have been limited, and elucidating the underlying genomic and molecular mechanisms of SCZ is still challenging. More recently, multiple Omics technologies - genomics, transcriptomics, epigenomics, proteomics, metabolomics, connectomics, and gut microbiomics - have all been applied to examine different aspects of SCZ pathogenesis. Integration of multi-Omics data has thus emerged as an approach to provide a more comprehensive view of biological complexity, which is vital to enable translation into assessments and interventions of clinical benefit to individuals with SCZ. In this review, we provide a broad survey of the single-omics studies of SCZ, summarize the advantages and challenges of different Omics technologies, and then focus on studies in which multiple omics data are integrated to unravel the complex pathophysiology of SCZ. We believe that integration of multi-Omics technologies would provide a roadmap to create a more comprehensive picture of interactions involved in the complex pathogenesis of SCZ, constitute a rich resource for elucidating the potential molecular mechanisms of the illness, and eventually improve clinical assessments and interventions of SCZ to address clinical translational questions from bench to bedside.

Predicting death from COVID-19 using pre-existing conditions: implications for vaccination triage

Introduction

Global shortages in the supply of SARS-CoV-2 vaccines have resulted in campaigns to first inoculate individuals at highest risk for death from COVID-19. Here, we develop a predictive model of COVID-19-related death using longitudinal clinical data from patients in metropolitan Detroit.

Methods

All individuals included in the analysis had a laboratory-confirmed SARS-CoV-2 infection. Thirty-six pre-existing conditions with a false discovery rate p<0.05 were combined with other demographic variables to develop a parsimonious prediction model using least absolute shrinkage and selection operator regression. The model was then prospectively validated in a separate set of individuals with confirmed COVID-19.

Results

The study population consisted of 15 502 individuals with laboratory-confirmed SARS-CoV-2. The main prediction model was developed using data from 11 635 individuals with 709 reported deaths (case fatality ratio 6.1%). The final prediction model consisted of 14 variables with 11 comorbidities. This model was then prospectively assessed among the remaining 3867 individuals (185 deaths; case fatality ratio 4.8%). When compared with using an age threshold of 65 years, the 14-variable model detected 6% more of the individuals who would die from COVID-19. However, below age 45 years and its risk equivalent, there was no benefit to using the prediction model over age alone.

Discussion

Using a prediction model, such as the one described here, may help identify individuals who would most benefit from COVID-19 inoculation, and thereby may produce more dramatic initial drops in deaths through targeted vaccination.

Serious Suicide Attempts and Risk of Suicide Death

Background: In the US, more than one million people attempt suicide each year. History of suicide attempt is a significant risk factor for death by suicide; however, there is a paucity of data from the US general population on this relationship. Aim: The objective of this study was to examine suicide attempts needing medical attention as a risk for suicide death. Method: We conducted a case-control study involving eight US healthcare systems. A total of 2,674 individuals who died by suicide from 2000 to 2013 were matched to 267,400 individuals by year and location. Results: Prior suicide attempt associated with a medical visit increases risk for suicide death by 39.1 times, particularly for women (OR = 79.2). However, only 11.3% of suicide deaths were associated with an attempt that required medical attention. The association was the strongest for children 10-14 years old (OR = 98.0). Most suicide attempts were recorded during the 20-week period prior to death. Limitations: Our study is limited to suicide attempts for which individuals sought medical care. Conclusion: In the US, prior suicide attempt is associated with an increased risk of suicide death; the risk is high especially during the period immediately following a nonlethal attempt.

Plasma Proteomic Profile Predicts Survival in Heart Failure With Reduced Ejection Fraction

BACKGROUND

It remains unclear whether the plasma proteome adds value to established predictors in heart failure (HF) with reduced ejection fraction (HFrEF). We sought to derive and validate a plasma proteomic risk score (PRS) for survival in patients with HFrEF (HFrEF-PRS).

METHODS

Patients meeting Framingham criteria for HF with EF<50% were enrolled (N=1017) and plasma underwent SOMAscan profiling (4453 targets). Patients were randomly divided 2:1 into derivation and validation cohorts. The HFrEF-PRS was derived using Cox regression of all-cause mortality adjusted for clinical score and NT-proBNP (N-terminal pro-B-type natriuretic peptide), then was tested in the validation cohort. Risk stratification improvement was evaluated by C statistic, integrated discrimination index, continuous net reclassification index, and median improvement in risk score for 1-year and 3-year mortality.

RESULTS

Participants' mean age was 68 years, 48% identified as Black, 35% were female, and 296 deaths occurred. In derivation (n=681), 128 proteins associated with mortality, 8 comprising the optimized HFrEF-PRS. In validation (n=336) the HFrEF-PRS associated with mortality (hazard ratio, 2.27 [95% CI, 1.84-2.82], P=6.3×10^-14), Kaplan-Meier curves differed significantly between HFrEF-PRS quartiles (P=2.2×10^-6), and it remained significant after adjustment for clinical score and NT-proBNP (hazard ratio, 1.37 [95% CI, 1.05-1.79], P=0.021). The HFrEF-PRS improved metrics of risk stratification (C statistic change, 0.009, P=0.612; integrated discrimination index, 0.041, P=0.010; net reclassification index=0.391, P=0.078; median improvement in risk score=0.039, P=0.016) and associated with cardiovascular death and HF phenotypes (eg, 6-minute walk distance, EF change). Most HFrEF-PRS proteins had little known connection to HFrEF.

CONCLUSIONS

A plasma multiprotein score improved risk stratification in patients with HFrEF and identified novel candidates.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Mapping the 17q12-21.1 Locus for Variants Associated with Early-Onset Asthma in African Americans

Rationale: The 17q12-21.1 locus is one of the most highly replicated genetic associations with asthma. Individuals of African descent have lower linkage disequilibrium in this region, which could facilitate identifying causal variants.Objectives: To identify functional variants at 17q12-21.1 associated with early-onset asthma among African American individuals.Methods: We evaluated African American participants from SAPPHIRE (Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity) (n = 1,940), SAGE II (Study of African Americans, Asthma, Genes and Environment) (n = 885), and GCPD-A (Study of the Genetic Causes of Complex Pediatric Disorders-Asthma) (n = 2,805). Associations with asthma onset at ages under 5 years were meta-analyzed across cohorts. The lead signal was reevaluated considering haplotypes informed by genetic ancestry (i.e., African vs. European). Both an expression-quantitative trait locus analysis and a phenome-wide association study were performed on the lead variant.Measurements and Main Results: The meta-analyzed results from SAPPHIRE, SAGE II, and the GCPD-A identified rs11078928 as the top association for early-onset asthma. A haplotype analysis suggested that the asthma association partitioned most closely with the rs11078928 genotype. Genetic ancestry did not appear to influence the effect of this variant. In the expression-quantitative trait locus analysis, rs11078928 was related to alternative splicing of GSDMB (gasdermin-B) transcripts. The phenome-wide association study of rs11078928 suggested that this variant was predominantly associated with asthma and asthma-associated symptoms.Conclusions: A splice-acceptor polymorphism appears to be a causal variant for asthma at the 17q12-21.1 locus. This variant appears to have the same magnitude of effect in individuals of African and European descent.

Asthma and its relationship to mitochondrial copy number: Results from the Asthma Translational Genomics Collaborative (ATGC) of the Trans-Omics for Precision Medicine (TOPMed) program

Background

Mitochondria support critical cellular functions, such as energy production through oxidative phosphorylation, regulation of reactive oxygen species, apoptosis, and calcium homeostasis.

Objective

Given the heightened level of cellular activity in patients with asthma, we sought to determine whether mitochondrial DNA (mtDNA) copy number measured in peripheral blood differed between individuals with and without asthma.

Methods

Whole genome sequence data was generated as part of the Trans-Omics for Precision Medicine (TOPMed) Program on participants from the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-ethnicity (SAPPHIRE) and the Study of African Americans, Asthma, Genes, & Environment II (SAGE II). We restricted our analysis to individuals who self-identified as African American (3,651 asthma cases and 1,344 controls). Mitochondrial copy number was estimated using the sequencing read depth ratio for the mitochondrial and nuclear genomes. Respiratory complex expression was assessed using RNA-sequencing.

Results

Average mitochondrial copy number was significantly higher among individuals with asthma when compared with controls (SAPPHIRE: 218.60 vs. 200.47, P<0.001; SAGE II: 235.99 vs. 223.07, P<0.001). Asthma status was significantly associated with mitochondrial copy number after accounting for potential explanatory variables, such as participant age, sex, leukocyte counts, and mitochondrial haplogroup. Despite the consistent relationship between asthma status and mitochondrial copy number, the latter was not associated with time-to-exacerbation or patient-reported asthma control. Mitochondrial respiratory complex gene expression was disproportionately lower in individuals with asthma when compared with individuals without asthma and other protein-encoding genes.

Conclusions

We observed a robust association between asthma and higher mitochondrial copy number. Asthma having an effect on mitochondria function was also supported by lower respiratory complex gene expression in this group.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.¹ Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations leading to clonal expansion was associated with both hematologic cancer ^2–4 and coronary heart disease⁵, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).⁶ Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer resulting in increased hematopoietic stem cell self-renewal. Overall, we observe that germline genetic variation shapes hematopoietic stem cell function leading to CHIP through mechanisms that are both specific to clonal hematopoiesis and shared mechanisms leading to somatic mutations across tissues.

Polygenic Score for β-Blocker Survival Benefit in European Ancestry Patients With Reduced Ejection Fraction Heart Failure

BACKGROUND

β-Blockers (BBs) are mainstay therapy for heart failure with reduced ejection fraction. However, individual patient responses to BB vary, which may be partially due to genetic variation. The goal of this study was to derive and validate the first polygenic response predictor (PRP) for BB survival benefit in heart failure with reduced ejection fraction patients.

METHODS

Derivation and validation analyses were performed in n=1436 total HF patients of European descent and with ejection fraction <50%. The PRP was derived in a random subset of the Henry Ford Heart Failure Pharmacogenomic Registry (n=248) and then validated in a meta-analysis of the remaining patients from Henry Ford Heart Failure Pharmacogenomic Registry (n=247), the TIME-CHF (Trial of Intensified Versus Standard Medical Therapy in Elderly Patients With Congestive Heart Failure; n=431), and HF-ACTION trial (Heart Failure: a Controlled Trial Investigating Outcomes of Exercise Training; n=510). The PRP was constructed from a genome-wide analysis of BB×genotype interaction predicting time to all-cause mortality, adjusted for Meta-Analysis Global Group in Chronic Heart Failure score, genotype, level of BB exposure, and BB propensity score.

RESULTS

Five-fold cross-validation summaries out to 1000 single-nucleotide polymorphisms identified optimal prediction with a 44 single-nucleotide polymorphism score and cutoff at the 30th percentile. In validation testing (n=1188), greater BB exposure was associated with reduced all-cause mortality in patients with low PRP score (n=251; hazard ratio, 0.19 [95% CI, 0.04-0.51]; P=0.0075) but not high PRP score (n=937; hazard ratio, 0.84 [95% CI, 0.53-1.3]; P=0.448)-a difference that was statistically significant (P interaction, 0.0235). Results were consistent regardless of atrial fibrillation, ejection fraction (≤40% versus 41%-50%), or when examining cardiovascular death.

CONCLUSIONS

Among patients of European ancestry with heart failure with reduced ejection fraction, a PRP distinguished patients who derived substantial survival benefit from BB exposure from a larger group that did not. Additional work is needed to prospectively test clinical utility and to develop PRPs for other population groups and other medications.

Suppression tumorigenicity 2 (ST2) turbidimetric immunoassay compared to enzyme-linked immunosorbent assay in predicting survival in heart failure patients with reduced ejection fraction

BACKGROUND

Suppressor of tumorigenicity 2 (ST2) is a powerful marker of prognosis and treatment response in heart failure (HF), however, it is an enzyme-linked immunosorbent assay (ELISA) which may be cumbersome and costly. A turbidimetric immunoassay (TIA) that can run on common chemistry analyzers could overcome this. We studied a novel TIA for ST2, comparing it to commercial ST2 (ELISA).

METHODS

Patients age ≥ 18 years meeting Framingham definition for HF were enrolled in a prospective registry (Oct 2007 - March 2015) at Henry Ford Hospital and donated blood samples. Participants with reduced ejection fraction (<50%) and available plasma samples were included and valid ST2 measurements were obtained on the same sample using both TIA and ELISA (N = 721). The primary endpoint was all cause death. Correlation between the methods was quantified. The association with survival was tested using unadjusted and adjusted (for MAGGIC score and NTproBNP) Cox models and comparing the Area Under the Curve (AUC).

RESULTS

The inter-assay Spearman correlation coefficient was 0.77. Nonparametric regression showed no significant proportional difference (slope = 0.97) and a very small systematic difference (3.2 ng/mL). In univariate analyses, both TIA and ELISA ST2 were significant associates of survival with similar effect sizes (HR 4.46 and 3.50, respectively, both p < 0.001). In models adjusted for MAGGIC score, both ST2 remained significant in Cox models and incrementally improved AUC vs. MAGGIC alone (MAGGIC AUC = 0.757; TIA + MAGGIC AUC = 0.786, p = 0.025; ELISA + MAGGIC AUC = 0.793, p = 0.033). In models with both MAGGIC and NTproBNP included, both ST2 still remained significant but did not improve AUC.

CONCLUSIONS

A novel TIA method for ST2 quantification correlates highly with ELISA and offers similarly powerful risk-stratification.

Lung Function in African American Children with Asthma Is Associated with Novel Regulatory Variants of the KIT Ligand KITLG/SCF and Gene-By-Air-Pollution Interaction

Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV1), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole-genome sequencing data from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine project, we identified a novel genetic association with FEV1 on chromosome 12 in 867 African American children with asthma (P = 1.26 × 10-8, β = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded nine variants as the most likely candidates responsible for the association with FEV1 Hi-C data and expression QTL analysis demonstrated that these variants physically interacted with KITLG (KIT ligand, also known as SCF), and their minor alleles were associated with increased expression of the KITLG gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year ambient sulfur dioxide exposure (P = 0.003, β = 0.32). This study identified a novel protective genetic association with FEV1, possibly mediated through KITLG, in African American children with asthma. This is the first study that has identified a genetic association between lung function and KITLG, which has established a role in orchestrating allergic inflammation in asthma.

The Hyperechoic Appearance of the Deltoid Muscle on Shoulder Ultrasound Imaging as a Predictor of Diabetes and Prediabetes

OBJECTIVES

To evaluate whether the ultrasound appearance of the deltoid muscle in diabetic patients differs from that in obese nondiabetic patients.

METHODS

Ultrasound images of the deltoid muscle from 137 type 2 diabetic patients (including 13 prediabetic patients) and 49 obese nondiabetic patients were blindly reviewed by 2 musculoskeletal radiologists, and by a third when arbitration was needed, to determine whether the appearance was "normal," "suspected diabetes," or "definite diabetes." Age, sex, race, body mass index (BMI), insulin use, and hemoglobin A_1c were analyzed. This retrospective study included patients presenting between October 2005 and November 2017. Statistical analyses included a 2-sided sample t test or Wilcoxon rank sum test and a χ² or Fisher exact test. Statistical significance was defined as P < .05.

RESULTS

The type 2 diabetic patients included 98 women and 39 men aged 29 to 92 years, and the nondiabetic patients included 19 women and 30 men aged 18 to 75 years. A consensus diagnosis of definite diabetes by the musculoskeletal radiologists based on a hyperechoic deltoid was a powerful predictor of diabetes, with a positive predictive value of 89%. A hyperechoic deltoid was also a powerful predictor of prediabetes. Of the 13 prediabetic patients, all had the same hyperechoic appearance of the diabetic deltoid, regardless of BMI. Although obese diabetic patients more often had a diagnosis of definite diabetes, the BMI alone could not explain the increased echogenicity, as obese nondiabetic patients' deltoid muscles did not appear as hyperechoic and were correctly categorized as not having definite diabetes with 82% specificity.

CONCLUSIONS

The characteristic hyperechoic deltoid appearance is a strong predictor of both diabetes and prediabetes and differs from that of obese nondiabetic patients.

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.

Association of Regulatory Genetic Variants for Protein Kinase Cα with Mortality and Drug Efficacy in Patients with Heart Failure

PURPOSE

Protein kinase C alpha (gene: PRKCA) is a key regulator of cardiac contractility. Two genetic variants have recently been discovered to regulate PRKCA expression in failing human heart tissue (rs9909004 [T → C] and rs9303504 [C → G]). The association of those variants with clinical outcomes in patients with heart failure (HF), and their interaction with HF drug efficacy, is unknown.

METHODS

Patients with HF in a prospective registry starting in 2007 were genotyped by whole genome array (n = 951). The primary outcome was all-cause mortality. Cox proportional hazards models adjusted for established clinical risk factors and genomic ancestry tested the independent association of rs9909004 or rs9303504 and the variant interactions with cornerstone HF pharmacotherapies (beta-blockers or angiotensin-converting enzyme inhibitors/angiotensin receptor blockers) in additive genetic models.

RESULTS

The minor allele of rs9909004, but not of rs9303504, was independently associated with a decreased risk for all-cause mortality: adjusted HR = 0.81 (95% CI = 0.67-0.98), p = 0.032. The variants did not significantly interact with mortality benefit associated with cornerstone HF pharmacotherapies (p > 0.1 for all).

CONCLUSIONS

A recently discovered cardiac-specific regulatory variant for PRKCA (rs9909004) was independently associated with a decreased risk for all-cause mortality in patients with HF. The variant did not interact with mortality benefit associated with cornerstone HF pharmacotherapies.

How does race and ethnicity effect the precision treatment of asthma?

Introduction

Asthma is a common condition that affects large numbers of children and adults, yet the burden of disease is not equally distributed amongst groups. In the United States, African Americans and Puerto Ricans have higher rates of asthma and its complications when compared with European Americans. However, clinical trials and genetic studies have largely focused on the latter group.

Areas covered

Here we examine what is known regarding differences in asthma treatment response by race-ethnicity. We also review existing genetic studies related to the use of asthma medications, paying special attention to studies that included substantial numbers of non-white population groups. Publicly accessible search engines of the medical literature were queried using combinations of the terms asthma, race, ethnicity, pharmacogenomics, and pharmacogenetics, as well as the names of individual asthma medication classes. The list of articles reviewed was supplemented by bibliographies and expert knowledge.

Expert opinion

A substantial and coordinated effort is still needed to both identify and validate genetic biomarkers of asthma medication response, as currently there are no clinically actionable genetic markers available for this purpose. The path to identifying such markers in non-white populations is even more formidable, since these groups are underrepresented in existing data.

P6336Plasma proteomic profile predicts survival in heart failure with reduced ejection fraction

Background

Whether the plasma proteome can predict the course of heart failure (HF) and has incremental value to established predictors is uncertain.

Methods

Patients meeting Framingham HF criteria with history of reduced ejection fraction (n=1017) were prospectively enrolled in a registry and donated fasting blood samples. Plasma underwent analysis on the SOMAscan proteomic discovery platform, quantifying 4789 proteins using standard assay and quality controls. Patients were randomly divided into derivation (n=681) and validation (n=336) cohorts. We derived a proteomic risk score (PRS) in the derivation cohort using Lasso-penalized Cox regression and then tested it in the validation cohort. Both models were adjusted for an establish HF clinical risk score (MAGGIC) and NTproBNP. We assessed risk stratification improvement in the validation cohort by comparing models with and without PRS using the model C statistic, continuous net reclassification index (NRI), integrated discrimination index (IDI), and the median improvement in risk score (MIRS).

Results

Overall 47.5% of patients were African American, 35.2% were female, mean ejection fraction was 34.8%, and average age was 67.9 years. After median follow-up of 3.6 years, there were 296 deaths (194 in derivation and 102 in validation). Optimized modeling defined a 32 protein PRS (hazard ratio [HR] 2.33, p<2.00E-16) which was also statistically significant when tested in the validation cohort (PRS HR=1.19, p=4.87E-02) and showed some improvement in risk stratification (Table).

Methods Variables Estimate 95% CI P Validation Testing MAGGIC 1.06 1.027, 1.092 2.84E-04 NTproBNP 1.84 1.430, 2.359 1.88E-06 PRS 1.19 1.001, 1.408 4.87E-02 Risk Stratification Assessment C-statistic improvement 0.012 −0.076, 0.101 8.30E-01 IDI 0.034 0.007, 0.095 <2.00E-16 Continuous-NRI 0.286 −0.062, 0.475 9.00E-02 Median Improvement in Risk Score 0.015 0.001, 0.078 2.00E-02

Conclusion

A plasma multi-protein predictive score can improve risk stratification in HF patients on top of a validated clinical score and NTproBNP. Additional investigation is warranted to define mechanisms underlying individual proteins and explore proteomic clinical applications.

Diagnosed Mental Health Conditions and Risk of Suicide Mortality

OBJECTIVE

Although mental health conditions are risk factors for suicide, limited data are available on suicide mortality associated with specific mental health conditions in the U.S. population. This study aimed to fill this gap.

METHODS

This study used a case-control design. Patients in the case group were those who died by suicide between 2000 and 2013 and who were patients in eight health care systems in the Mental Health Research Network (N=2,674). Each was matched with 100 general population patients from the same system (N=267,400). Diagnostic codes for five mental health conditions in the year before death were obtained from medical records: anxiety disorders, attention deficit-hyperactivity disorder (ADHD), bipolar disorder, depressive disorders, and schizophrenia spectrum disorder.

RESULTS

Among patients in the case group, 51.3% had a recorded psychiatric diagnosis in the year before death, compared with 12.7% of control group patients. Risk of suicide mortality was highest among those with schizophrenia spectrum disorder, after adjustment for age and sociodemographic characteristics (adjusted odds ratio [AOR]=15.0) followed by bipolar disorder (AOR=13.2), depressive disorders (AOR=7.2), anxiety disorders (AOR=5.8), and ADHD (AOR=2.4). The risk of suicide death among those with a diagnosed bipolar disorder was higher in women than men.

CONCLUSIONS

Half of those who died by suicide had at least one diagnosed mental health condition in the year before death, and most mental health conditions were associated with an increased risk of suicide. Findings suggest the importance of suicide screening and providing an approach to improve awareness of mental health conditions.

Performance of the Meta-Analysis Global Group in Chronic Heart Failure Score in Black Patients Compared With Whites

BACKGROUND

Risk stratification is critical in heart failure (HF) and the Meta-Analysis Global Group in Chronic HF (MAGGIC) score is a validated tool derived from ~40,000 patients. However, few of these patients self-identified as black, raising uncertainty regarding performance in blacks with HF.

METHODS AND RESULTS

This study analyzed a racially diverse group of 4046 patients (1646 black and 2400 white) from a single center from 2007 to 2015. Baseline characteristics were collected to tabulate MAGGIC score and test its discrimination and calibration within race groups. The primary end point was all-cause mortality. Death was detected using system records and the social security death master file. Discrimination was tested using Cox models of MAGGIC score stratified by race, and combined analysis including MAGGIC, race, and MAGGIC×race. Calibration was assessed using linear regression models and plots of observed versus predicted data. Overall, 901 (21%) patients died during 1-year follow-up. MAGGIC score discrimination was similar in both race groups in terms of C statistic (0.707±0.027 versus 0.725±0.014, for black versus white; P=0.556) and the hazard ratio (HR) per MAGGIC point was 1.12 in black patients (95% CI, 1.10-1.14) and 1.13 in white patients (95% CI, 1.12-1.14). Race was a significant correlate of survival, with better survival in black patients compared with white (HR, 0.66; 95% CI, 0.56-0.78), but the interaction of MAGGIC×race was not significant (β=-0.013; P=0.16), and adding race to the model did not improve discrimination (C statistic for MAGGIC versus MAGGIC+race, 0.721 versus 0.722; P=0.79). In calibration testing, the slope was not significantly different from 1 in either group, but the groups differed from each other, and it was closer to unity among black patients (0.94 versus 1.4; P=0.004).

CONCLUSIONS

These data support the use of the MAGGIC score to risk stratify black patients with HF.

An admixture mapping meta-analysis implicates genetic variation at 18q21 with asthma susceptibility in Latinos

BACKGROUND

Asthma is a common but complex disease with racial/ethnic differences in prevalence, morbidity, and response to therapies.

OBJECTIVE

We sought to perform an analysis of genetic ancestry to identify new loci that contribute to asthma susceptibility.

METHODS

We leveraged the mixed ancestry of 3902 Latinos and performed an admixture mapping meta-analysis for asthma susceptibility. We replicated associations in an independent study of 3774 Latinos, performed targeted sequencing for fine mapping, and tested for disease correlations with gene expression in the whole blood of more than 500 subjects from 3 racial/ethnic groups.

RESULTS

We identified a genome-wide significant admixture mapping peak at 18q21 in Latinos (P = 6.8 × 10-6), where Native American ancestry was associated with increased risk of asthma (odds ratio [OR], 1.20; 95% CI, 1.07-1.34; P = .002) and European ancestry was associated with protection (OR, 0.86; 95% CI, 0.77-0.96; P = .008). Our findings were replicated in an independent childhood asthma study in Latinos (P = 5.3 × 10-3, combined P = 2.6 × 10-7). Fine mapping of 18q21 in 1978 Latinos identified a significant association with multiple variants 5' of SMAD family member 2 (SMAD2) in Mexicans, whereas a single rare variant in the same window was the top association in Puerto Ricans. Low versus high SMAD2 blood expression was correlated with case status (13.4% lower expression; OR, 3.93; 95% CI, 2.12-7.28; P < .001). In addition, lower expression of SMAD2 was associated with more frequent exacerbations among Puerto Ricans with asthma.

CONCLUSION

Ancestry at 18q21 was significantly associated with asthma in Latinos and implicated multiple ancestry-informative noncoding variants upstream of SMAD2 with asthma susceptibility. Furthermore, decreased SMAD2 expression in blood was strongly associated with increased asthma risk and increased exacerbations.

Author Correction: Assembly of a pan-genome from deep sequencing of 910 humans of African descent

In the version of this article initially published, the statement "there are no pan-genomes for any other animal or plant species" was incorrect. The statement has been corrected to "there are no reported pan-genomes for any other animal species, to our knowledge." We thank David Edwards for bringing this error to our attention. The error has been corrected in the HTML and PDF versions of the article.

Secondhand smoke exposure and asthma outcomes among African-American and Latino children with asthma

BACKGROUND

Secondhand smoke (SHS) exposures have been linked to asthma-related outcomes but quantitative dose-responses using biomarkers of exposure have not been widely reported.

OBJECTIVES

Assess dose-response relationships between plasma cotinine-determined SHS exposure and asthma outcomes in minority children, a vulnerable population exposed to higher levels of SHS and under-represented in the literature.

METHODS

We performed analyses in 1172 Latino and African-American children with asthma from the mainland USA and Puerto Rico. We used logistic regression to assess relationships of cotinine levels ≥0.05 ng/mL with asthma exacerbations (defined as asthma-related hospitalisations, emergency room visits or oral steroid prescription) in the previous year and asthma control. The shape of dose-response relationships was assessed using a continuous exposure variable in generalised additive logistic models with penalised splines.

RESULTS

The OR for experiencing asthma exacerbations in the previous year for cotinine levels ≥0.05 ng/mL, compared with <0.05 ng/mL, was 1.40 (95% CI 1.03 to 1.89), while the OR for poor asthma control was 1.53 (95% CI 1.12 to 2.13). Analyses for dose-response relationships indicated increasing odds of asthma outcomes related with increasing exposure, even at cotinine levels associated with light SHS exposures.

CONCLUSIONS

Exposure to SHS was associated with higher odds of asthma exacerbations and having poorly controlled asthma with an increasing dose-response even at low levels of exposure. Our results support the conclusion that there are no safe levels of SHS exposures.

Reliably colonising broiler chickens with Campylobacter spp. using a litter-based method

1. Chicken-associated Campylobacter spp. are the cause of most food poisoning cases in Europe. In order to study the host-pathogen interactions, a reliable and reproducible method of colonising chickens with the bacteria is required. 2. This study aimed to identify a more appropriate and less invasive method of colonisation (cf. gavaging) by seeding bedding material (litter) that commercial chickens are kept on with a mixture of Campylobacter spp., broth and faeces. 3. The first phase of the study tested the longevity of Campylobacter spp. recovery in seeded litter over 24 h: significantly more Campylobacter spp. was recovered at 0 or 3 h post-seeding than at 6 and 24 h post-seeding, indicating that the pathogen can survive to detectable levels for at least 3 h in this environment. 4. In the second phase, three groups of 10 broiler chickens (negative for Campylobacter spp. prior to exposure) were exposed at 21 days of age to one of three different Campylobacter jejuni and C. coli mixes (A, B, C), using the method above. At 28 days of age, birds were euthanised by overdose of barbiturate or cervical dislocation, and livers and caeca removed for Campylobacter spp. assessment. 5. All liver and 28/30 caeca samples tested positive for Campylobacter spp., with mix A and C giving higher counts in the caeca than mix B. The method of euthanasia did not affect Campylobacter spp. counts. 6. In conclusion, a successful method for reliably colonising broiler chickens with Campylobacter spp. has been developed which negates the need for gavaging and is more representative of how contamination occurs in the field.

A genome-wide association and admixture mapping study of bronchodilator drug response in African Americans with asthma

Short-acting β₂-adrenergic receptor agonists (SABAs) are the most commonly prescribed asthma medications worldwide. Response to SABAs is measured as bronchodilator drug response (BDR), which varies among racial/ethnic groups in the U.S^{1, 2}. However, the genetic variation that contributes to BDR is largely undefined in African Americans with asthma³. To identify genetic variants that may contribute to differences in BDR in African Americans with asthma, we performed a genome-wide association study (GWAS) of BDR in 949 African American children with asthma, genotyped with the Axiom World Array 4 (Affymetrix, Santa Clara, CA) followed by imputation using 1000 Genomes phase III genotypes. We used linear regression models adjusting for age, sex, body mass index (BMI) and genetic ancestry to test for an association between BDR and genotype at single nucleotide polymorphisms (SNPs). To increase power and distinguish between shared vs. population-specific associations with BDR in children with asthma, we performed a meta-analysis across 949 African Americans and 1,830 Latinos (Total=2,779). Lastly, we performed genome-wide admixture mapping to identify regions whereby local African or European ancestry is associated with BDR in African Americans. We identified a population-specific association with an intergenic SNP on chromosome 9q21 that was significantly associated with BDR (rs73650726, p=7.69×10⁻⁹). A trans-ethnic meta-analysis across African Americans and Latinos identified three additional SNPs within the intron of PRKG1 that were significantly associated with BDR (rs7903366, rs7070958, and rs7081864, p≤5×10⁻⁸). Our results failed to replicate in three additional populations of 416 Latinos and 1,615 African Americans. Our findings indicate that both population specific and shared genetic variation contributes to differences in BDR in minority children with asthma, and that the genetic underpinnings of BDR may differ between racial/ethnic groups.

Author Correction: Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis

In the version of this article initially published, in Fig. 3, the y-axis numbering did not match the log scale indicated in the axis label. The error has been corrected in the HTML and PDF version of the article.