Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) lncRNA differentially regulates gene and protein expression, signaling and morphology of human ocular cells.

Pseudoexfoliation glaucoma (PEXG) is characterized by dysregulated extracellular matrix (ECM) homeostasis that disrupts conventional outflow function and increases intraocular pressure (IOP). Prolonged IOP elevation results in optic nerve head damage and vision loss. Uniquely, PEXG is a form of open angle glaucoma that has variable penetrance, is difficult to treat and does not respond well to common IOP-lowering pharmaceuticals. Therefore, understanding modulators of disease severity will aid in targeted therapies for PEXG. Genome-wide association studies have identified polymorphisms in the long non-coding RNA lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) as a risk factor for PEXG. Risk alleles, oxidative stress and mechanical stretch all alter LOXL1-AS1 expression. As a long non-coding RNA, LOXL1-AS1 binds hnRNPL and regulates global gene expression. In this study, we focus on the role of LOXL1-AS1 in the ocular cells (trabecular meshwork and Schlemm's canal) that regulate IOP. We show that selective knockdown of LOXL1-AS1 leads to cell-type-specific changes in gene expression, ECM homeostasis, signaling and morphology. These results implicate LOXL1-AS1 as a modulator of cellular homeostasis, altering cell contractility and ECM turnover, both of which are well-known contributors to PEXG. These findings support LOXL1-AS1 as a key target for modifying the disease.

Nontargeted Plasma Proteomic Analysis of Renal Disease and Pulmonary Hypertension in Patients with Sickle Cell Disease.

Sickle cell disease (SCD) is characterized by red blood cell sickling, vaso-occlusion, hemolytic anemia, damage to multiple organ systems, and, as a result, shortened life expectancy. Sickle cell disease nephropathy (SCDN) and pulmonary hypertension (pHTN) are common and frequently co-occurring complications of SCD; both are associated with markedly accelerated mortality. To identify candidate circulating biomarkers of SCDN and pHTN, we used mass spectrometry to quantify the relative abundance of >1000 proteins in plasma samples from 189 adults with SCD from the Outcome Modifying Genes in SCD (OMG-SCD) cohort (ProteomeXchange identifier PXD048716). Forty-four proteins were differentially abundant in SCDN, most significantly cystatin-C and collagen α-1(XVIII) chain (COIA1), and 55 proteins were dysregulated in patients with SCDN and pHTN, most significantly insulin-like growth factor-binding protein 6 (IBP6). Network analysis identified a module of 133 coregulated proteins significantly associated with SCDN, that was enriched for extracellular matrix proteins, insulin-like growth factor binding proteins, cell adhesion proteins, EGF-like calcium binding proteins, and several cadherin family members. Collectively, these data provide a comprehensive understanding of plasma protein changes in SCDN and pHTN which validate numerous studies of chronic kidney disease and suggest shared profiles of protein disruption in kidney dysfunction and pHTN among SCD patients.

Rare and de novo coding variants in chromodomain genes in Chiari I malformation.

Chiari I malformation (CM1), the displacement of the cerebellum through the foramen magnum into the spinal canal, is one of the most common pediatric neurological conditions. Individuals with CM1 can present with neurological symptoms, including severe headaches and sensory or motor deficits, often as a consequence of brainstem compression or syringomyelia (SM). We conducted whole-exome sequencing (WES) on 668 CM1 probands and 232 family members and performed gene-burden and de novo enrichment analyses. A significant enrichment of rare and de novo non-synonymous variants in chromodomain (CHD) genes was observed among individuals with CM1 (combined p = 2.4 × 10-10), including 3 de novo loss-of-function variants in CHD8 (LOF enrichment p = 1.9 × 10-10) and a significant burden of rare transmitted variants in CHD3 (p = 1.8 × 10-6). Overall, individuals with CM1 were found to have significantly increased head circumference (p = 2.6 × 10-9), with many harboring CHD rare variants having macrocephaly. Finally, haploinsufficiency for chd8 in zebrafish led to macrocephaly and posterior hindbrain displacement reminiscent of CM1. These results implicate chromodomain genes and excessive brain growth in CM1 pathogenesis.

A genome-wide association study of alloimmunization in the TOPMed OMG-SCD cohort identifies a locus on chromosome 12.

BACKGROUND: Red cell alloimmunization after exposure to donor red cells is a very common complication of transfusion for patients with sickle cell disease (SCD), resulting frequently in accelerated donor red blood cell destruction. Patients show substantial differences in their predisposition to alloimmunization, and genetic variability is one proposed component. Although several genetic association studies have been conducted for alloimmunization, the results have been inconsistent, and the genetic determinants of alloimmunization remain largely unknown. STUDY DESIGN AND METHODS: We performed a genome-wide association study (GWAS) in 236 African American (AA) SCD patients from the Outcome Modifying Genes in Sickle Cell Disease (OMG-SCD) cohort, which is part of Trans-Omics for Precision Medicine (TOPMed), with whole-genome sequencing data available. We also performed sensitivity analyses adjusting for different sets of covariates and applied different sample grouping strategies based on the number of alloantibodies patients developed. RESULTS: We identified one genome-wide significant locus on chr12 (p = 3.1e-9) with no evidence of genomic inflation (lambda = 1.003). Further leveraging QTL evidence from GTEx whole blood and/or Jackson Heart Study PBMC RNA-Seq data, we identified a number of potential genes, such as ARHGAP9, STAT6, and ATP23, that may be driving the association signal. We also discovered some suggestive loci using different analysis strategies. DISCUSSION: We call for the community to collect additional alloantibody information within SCD cohorts to further the understanding of the genetic basis of alloimmunization in order to improve transfusion outcomes.

Variation and impact of polygenic hematologic traits in monogenic sickle cell disease.

Several of the complications observed in sickle cell disease (SCD) are influenced by variation in hematologic traits (HT), such as fetal hemoglobin (HbF) level and neutrophil count. Previous large-scale genome-wide association studies carried out in largely healthy individuals have identified thousands of variants associated with HT, which have then been used to develop multi-ancestry polygenic trait scores (PTS). Here, we tested whether these PTS associate with HT in SCD patients and if they can improve statistical models associated with SCD-related complications. In 2,056 SCD patients, we found that the PTS predicted less HT variance than in non-SCD individuals of African ancestry. This was particularly striking at the Duffy/DARC locus, where we observed an epistatic interaction between the SCD genotype and the Duffy null variant (rs2814778) that led to a two-fold weaker effect on neutrophil count. PTS for these HT which are measured as part of routine practice were not associated with complications in SCD. In contrast, we found that a simple PTS for HbF that includes only six variants explained a large fraction of the phenotypic variation (20.5-27.1%), associated with acute chest syndrome and stroke risk, and improved the statistical modeling of the vaso-occlusive crisis rate. Using Mendelian randomization, we found that increasing HbF by 4.8% reduces stroke risk by 39% (P=0.0006). Taken together, our results highlight the importance of validating PTS in large diseased populations before proposing their implementation in the context of precision medicine initiatives.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci.

To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.

Evaluating Associations between Average Pain Intensity and Genetic Variation in People with Sickle Cell Disease: An Exploratory Study.

BACKGROUND: Pain is one of the most common and deleterious symptoms experienced by individuals with sickle cell disease (SCD). There is a paucity of studies identifying potential genetic mechanisms of pain in this population. AIM: Examine associations between 11 functional single nucleotide polymorphisms in 9 candidate genes with reports of average pain intensity in individuals with sickle cell disease. METHOD: Cross-sectional analyses were performed on data and blood samples collected through the Duke SCD Implementation Consortium Registry. Participants were asked to rate their pain "on the average" using an 11-point numeric rating scale (0 = no pain; 10 = pain as bad as you can imagine). We genotyped 11 single nucleotide polymorphisms in 9 pain-related genes using TaqMan® Genotyping Assays. Associations between each polymorphism and reports of average pain were evaluated. RESULTS: The 86 participants (mean age: 28.7 years; 64% female) included in this study reported moderate pain on average (Mean = 4, Standard Deviation = 2.4). ICAM1 rs1799969 was the only genetic polymorphism that was significantly associated with pain (p = .01). Individuals with one or more minor alleles had lower average pain (Mean = 1.25, Standard Deviation = 1.50) than individuals without a minor allele (Mean = 4.13, Standard Deviation = 2.25). The effect size for ICAM1 rs1799969 was 1.30, which is considered large. The effect sizes for all other single nucleotide polymorphisms ranged from small to medium (range: 0-0.3). CONCLUSIONS: Our findings provide preliminary evidence that the minor allele in ICAM1 rs1799969 had protective effects against experiencing more severe pain in sickle cell disease.

A Genocentric Approach to Discovery of Mendelian Disorders.

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

The TFAP2A-IRF6-GRHL3 genetic pathway is conserved in neurulation.

Mutations in IRF6, TFAP2A and GRHL3 cause orofacial clefting syndromes in humans. However, Tfap2a and Grhl3 are also required for neurulation in mice. Here, we found that homeostasis of Irf6 is also required for development of the neural tube and associated structures. Over-expression of Irf6 caused exencephaly, a rostral neural tube defect, through suppression of Tfap2a and Grhl3 expression. Conversely, loss of Irf6 function caused a curly tail and coincided with a reduction of Tfap2a and Grhl3 expression in tail tissues. To test whether Irf6 function in neurulation was conserved, we sequenced samples obtained from human cases of spina bifida and anencephaly. We found two likely disease-causing variants in two samples from patients with spina bifida. Overall, these data suggest that the Tfap2a-Irf6-Grhl3 genetic pathway is shared by two embryologically distinct morphogenetic events that previously were considered independent during mammalian development. In addition, these data suggest new candidates to delineate the genetic architecture of neural tube defects and new therapeutic targets to prevent this common birth defect.

Potential causal role of l-glutamine in sickle cell disease painful crises: A Mendelian randomization analysis.

In a recent clinical trial, the metabolite l-glutamine was shown to reduce painful crises in sickle cell disease (SCD) patients. To support this observation and identify other metabolites implicated in SCD clinical heterogeneity, we profiled 129 metabolites in the plasma of 705 SCD patients. We tested correlations between metabolite levels and six SCD-related complications (painful crises, cholecystectomy, retinopathy, leg ulcer, priapism, aseptic necrosis) or estimated glomerular filtration rate (eGFR), and used Mendelian randomization (MR) to assess causality. We found a potential causal relationship between l-glutamine levels and painful crises (N = 1278, odds ratio (OR) [95% confidence interval] = 0.68 [0.52-0.89], P = 0.0048). In two smaller SCD cohorts (N = 299 and 406), the protective effect of l-glutamine was observed (OR = 0.82 [0.50-1.34]), although the MR result was not significant (P = 0.44). We identified 66 significant correlations between the levels of other metabolites and SCD-related complications or eGFR. We tested these correlations for causality using MR analyses and found no significant causal relationship. The baseline levels of quinolinic acid were associated with prospectively ascertained survival in SCD patients, and this effect was dependent on eGFR. Metabolomics provide a promising approach to prioritize small molecules that may serve as biomarkers or drug targets in SCD.

Rapid decline in estimated glomerular filtration rate in sickle cell anemia: results of a multicenter pooled analysis.

Not available.

Diversity of variant alleles encoding Kidd, Duffy, and Kell antigens in individuals with sickle cell disease using whole genome sequencing data from the NHLBI TOPMed Program.

BACKGROUND: Genetic variants in the SLC14A1, ACKR1, and KEL genes, which encode Kidd, Duffy, and Kell red blood cell antigens, respectively, may result in weakened expression of antigens or a null phenotype. These variants are of particular interest to individuals with sickle cell disease (SCD), who frequently undergo chronic transfusion therapy with antigen-matched units. The goal was to describe the diversity and the frequency of variants in SLC14A1, ACKR1, and KEL genes among individuals with SCD using whole genome sequencing (WGS) data. STUDY DESIGN AND METHODS: Two large SCD cohorts were studied: the Recipient Epidemiology and Donor Evaluation Study III (REDS-III) (n = 2634) and the Outcome Modifying Gene in SCD (OMG) (n = 640). Most of the studied individuals were of mixed origin. WGS was performed as part of the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program. RESULTS: In SLC14A1, variants included four encoding a weak Jka phenotype and five null alleles (JKnull ). JKA*01N.09 was the most common JKnull . One possible JKnull mutation was novel: c.812G>T. In ACKR1, identified variants included two that predicted Fyx (FY*X) and one corresponding to the c.-67T>C GATA mutation. The c.-67T>C mutation was associated with FY*A (FY*01N.01) in four participants. FY*X was identified in 49 individuals. In KEL, identified variants included three null alleles (KEL*02N.17, KEL*02N.26, and KEL*02N.04) and one allele predicting Kmod phenotype, all in heterozygosity. CONCLUSIONS: We described the diversity and distribution of SLC14A1, ACKR1, and KEL variants in two large SCD cohorts, comprising mostly individuals of mixed ancestry. This information may be useful for planning the transfusion support of patients with SCD.

Genetic variation in dopamine neurotransmission and motor development of infants born extremely-low-birthweight.

AIM: To determine if genetic variation associated with decreased dopamine neurotransmission predicts a decrease in motor development in a convenience cohort study of infants born extremely-low-birthweight (ELBW). METHOD: Four hundred and ninety-eight infants born ELBW had genome-wide genotyping and a neurodevelopmental evaluation at 18 to 22 months of age, corrected for preterm birth. A polygenic risk score (PRS) was created to combine into one predictor variable the hypothesized influences on motor development of alleles at seven independent single nucleotide polymorphisms previously associated with relative decreases in both dopamine neurotransmission and motor learning, by summing the number of alleles present in each infant (range=0-14). The motor development outcome was the Psychomotor Development Index (PDI) of the Bayley Scales of Infant Development, Second Edition. The linear regression models were adjusted for seven clinical and four genetic ancestry covariates. The mean PRS of infants with cerebral palsy (CP) was compared to those without CP. RESULTS: PRS was inversely related to PDI (p=0.011). Each 1-point increase in PRS resulted in an average decrease in PDI of 1.37 points. Patients with CP did not have a greater mean PRS than those without (p=0.67), both with and without adjustment for covariates. INTERPRETATION: Genetic variation that favors a decrease in dopamine neurotransmission predisposes to a decrease in motor development in infants born ELBW, but not to the diagnosis of CP. WHAT THIS PAPER ADDS: Genetic variation in dopamine neurotransmission was associated with a decrease in motor development in infants born at an extremely-low-birthweight. It does not predispose to the diagnosis of cerebral palsy.

Transcriptome analysis of adult and fetal trabecular meshwork, cornea, and ciliary body tissues by RNA sequencing.

PURPOSE: To characterize the transcriptional landscape of human adult and fetal trabecular meshwork (TM), cornea, and ciliary body (CB) tissues, and to evaluate the expression level of candidate genes selected from genetic association studies of primary-open angle glaucoma, central corneal thickness, intraocular pressure, vertical cup to disc ratio, and optic nerve parameters. METHODS: Deep RNA sequencing was performed on the selected human tissues. Transcriptome analyses were performed to 1) characterize the total number of expressed genes, 2) identify the most highly expressed genes, 3) estimate the number of novel transcripts, and 4) evaluate the expression of candidate genes in each tissue. Finally, a differential gene expression analysis was conducted to compare the adult and fetal ocular tissues. RESULTS: There was an average of 12,362 protein coding genes and 3725 novel transcripts expressed in each tissue. The top most expressed genes in each tissue included SPARC (fetal cornea and TM), APOD (adult TM), CLU (adult cornea), and PTGDS (adult and fetal CB). Twenty-nine candidate genes selected from genetic association studies primarily showed high expression levels in the trabecular meshwork and cornea. Comparison of adult and fetal samples identified 2012 and 1261 differentially expressed protein-coding genes within the cornea and trabecular meshwork, respectively. CONCLUSIONS: This study has provided an unbiased glimpse into the transcriptome of three essential anterior ocular tissues, resulting in the development of several novel hypotheses. These data can be used in the future to better guide ocular research questions.

Clinical and metabolomic risk factors associated with rapid renal function decline in sickle cell disease.

Sickle cell disease (SCD) nephropathy and lower estimated glomerular filtration rate (eGFR) are risk factors for early mortality. Furthermore, rate of eGFR decline predicts progression to end-stage renal disease in many clinical settings. However, factors predicting renal function decline in SCD are poorly documented. Using clinical, laboratory, genetic, and metabolomic data, we evaluated predictors of renal function decline in a longitudinal cohort of 288 adults (mean age 33.0 years). In 193 subjects with 5-year follow-up data, mean rate of eGFR decline was 2.35 mL/min/1.73 m2 /year, nearly twice that of African American adults overall. Hyperfiltration was prevalent at baseline (61.1%), and 36.8% of subjects experienced rapid eGFR decline (≥3 mL/min/1.73 m2 /year). Severe Hb genotype; proteinuria; higher platelet and reticulocyte counts, and systolic BP; and lower Hb level and BMI were associated with rapid decline. A risk scoring system was created using these 7 variables and was highly predictive of rapid eGFR decline, with odds of rapid decline increasing 1.635-fold for every point increment (P < 0.0001). Rapid eGFR decline was also associated with higher organ system severity score and peak creatinine. Additionally, two metabolites (asymmetric dimethylarginine and quinolinic acid) were associated with rapid decline. Further investigation into longitudinal SCD nephropathy (SCDN) trajectory, early markers of SCDN, and tools for risk stratification should inform interventional studies targeted to slowing GFR decline and improving SCD outcomes.