Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Genetic examination of hematological parameters in SARS-CoV-2 infection and COVID-19.

People hospitalized with COVID-19 often exhibit altered hematological traits associated with disease prognosis (e.g., lower lymphocyte and platelet counts). We investigated whether inter-individual variability in baseline hematological traits influences risk of acute SARS-CoV-2 infection or progression to severe COVID-19. We report inconsistent associations between blood cell traits with incident SARS-CoV-2 infection and severe COVID-19 in UK Biobank and the Vanderbilt University Medical Center Synthetic Derivative (VUMC SD). Since genetically determined blood cell measures better represent cell abundance across the lifecourse, we also assessed the shared genetic architecture of baseline blood cell traits on COVID-19 related outcomes by Mendelian randomization (MR) analyses. We found significant relationships between COVID-19 severity and mean sphered cell volume after adjusting for multiple testing. However, MR results differed significantly across different freezes of COVID-19 summary statistics and genetic correlation between these traits was modest (0.1), decreasing our confidence in these results. We observed overlapping genetic association signals between other hematological and COVID-19 traits at specific loci such as MAPT and TYK2. In conclusion, we did not find convincing evidence of relationships between the genetic architecture of blood cell traits and either SARS-CoV-2 infection or COVID-19 hospitalization, though we do see evidence of shared signals at specific loci.

Transcriptional control of parallel-acting pathways that remove specific presynaptic proteins in remodeling neurons.

Synapses are actively dismantled to mediate circuit refinement, but the developmental pathways that regulate synaptic disassembly are largely unknown. We have previously shown that the epithelial sodium channel ENaC/UNC-8 triggers an activity-dependent mechanism that drives the removal of presynaptic proteins liprin-α/SYD-2, Synaptobrevin/SNB-1, RAB-3 and Endophilin/UNC-57 in remodeling GABAergic neurons in C. elegans (Miller-Fleming et al., 2016). Here, we report that the conserved transcription factor Iroquois/IRX-1 regulates UNC-8 expression as well as an additional pathway, independent of UNC-8, that functions in parallel to dismantle functional presynaptic terminals. We show that the additional IRX-1-regulated pathway is selectively required for the removal of the presynaptic proteins, Munc13/UNC-13 and ELKS, which normally mediate synaptic vesicle fusion and neurotransmitter release. Our findings are notable because they highlight the key role of transcriptional regulation in synapse elimination during development and reveal parallel-acting pathways that coordinate synaptic disassembly by removing specific active zone proteins.SIGNIFICANT STATEMENT:Synaptic pruning is a conserved feature of developing neural circuits but the mechanisms that dismantle the presynaptic apparatus are largely unknown. We have determined that synaptic disassembly is orchestrated by parallel-acting mechanisms that target distinct components of the active zone. Thus, our finding suggests that synaptic disassembly is not accomplished by en masse destruction but depends on mechanisms that dismantle the structure in an organized process.

Analyses of biomarker traits in diverse UK biobank participants identify associations missed by European-centric analysis strategies.

Despite the dramatic underrepresentation of non-European populations in human genetics studies, researchers continue to exclude participants of non-European ancestry, as well as variants rare in European populations, even when these data are available. This practice perpetuates existing research disparities and can lead to important and large effect size associations being missed. Here, we conducted genome-wide association studies (GWAS) of 31 serum and urine biomarker quantitative traits in African (n = 9354), East Asian (n = 2559), and South Asian (n = 9823) ancestry UK Biobank (UKBB) participants. We adjusted for all known GWAS catalog variants for each trait, as well as novel signals identified in a recent European ancestry-focused analysis of UKBB participants. We identify 7 novel signals in African ancestry and 2 novel signals in South Asian ancestry participants (p < 1.61E-10). Many of these signals are highly plausible, including a cis pQTL for the gene encoding gamma-glutamyl transferase and PIEZO1 and G6PD variants with impacts on HbA1c through likely erythrocytic mechanisms. This work illustrates the importance of using the genetic data we already have in diverse populations, with novel discoveries possible in even modest sample sizes.

Ca2+ permeability and Na+ conductance in cellular toxicity caused by hyperactive DEG/ENaC channels.

Hyperactivated DEG/ENaC channels cause neuronal death mediated by intracellular Ca2+ overload. Mammalian ASIC1a channels and MEC-4(d) neurotoxic channels in Caenorhabditis elegans both conduct Na+ and Ca2+, raising the possibility that direct Ca2+ influx through these channels contributes to intracellular Ca2+ overload. However, we showed that the homologous C. elegans DEG/ENaC channel UNC-8(d) is not Ca2+ permeable, yet it is neurotoxic, suggesting that Na+ influx is sufficient to induce cell death. Interestingly, UNC-8(d) shows small currents due to extracellular Ca2+ block in the Xenopus oocyte expression system. Thus, MEC-4(d) and UNC-8(d) differ both in current amplitude and Ca2+ permeability. Given that these two channels show a striking difference in toxicity, we wondered how Na+ conductance vs. Ca2+ permeability contributes to cell death. To address this question, we built an UNC-8/MEC-4 chimeric channel that retains the calcium permeability of MEC-4 and characterized its properties in Xenopus oocytes. Our data support the hypothesis that for Ca2+-permeable DEG/ENaC channels, both Ca2+ permeability and Na+ conductance contribute to toxicity. However, for Ca2+-impermeable DEG/ENaCs (e.g., UNC-8), our evidence shows that constitutive Na+ conductance is sufficient to induce toxicity, and that this effect is enhanced as current amplitude increases. Our work further refines the contribution of different channel properties to cellular toxicity induced by hyperactive DEG/ENaC channels.

Developing a Phenotype Risk Score for Tic Disorders in a Large, Clinical Biobank.

Importance: Tics are a common feature of early-onset neurodevelopmental disorders, characterized by involuntary and repetitive movements or sounds. Despite affecting up to 2% of young children and having a genetic contribution, the underlying causes remain poorly understood, likely due to the complex phenotypic and genetic heterogeneity among affected individuals. Objective: In this study, we leverage dense phenotype information from electronic health records to identify the disease features associated with tic disorders within the context of a clinical biobank. These disease features are then used to generate a phenotype risk score for tic disorder. Design: Using de-identified electronic health records from a tertiary care center, we extracted individuals with tic disorder diagnosis codes. We performed a phenome-wide association study to identify the features enriched in tic cases versus controls (N=1,406 and 7,030; respectively). These disease features were then used to generate a phenotype risk score for tic disorder, which was applied across an independent set of 90,051 individuals. A previously curated set of tic disorder cases from an electronic health record algorithm followed by clinician chart review was used to validate the tic disorder phenotype risk score. Main Outcomes and Measures: Phenotypic patterns associated with a tic disorder diagnosis in the electronic health record. Results: Our tic disorder phenome-wide association study revealed 69 significantly associated phenotypes, predominantly neuropsychiatric conditions, including obsessive compulsive disorder, attention-deficit hyperactivity disorder, autism, and anxiety. The phenotype risk score constructed from these 69 phenotypes in an independent population was significantly higher among clinician-validated tic cases versus non-cases. Conclusions and Relevance: Our findings provide support for the use of large-scale medical databases to better understand phenotypically complex diseases, such as tic disorders. The tic disorder phenotype risk score provides a quantitative measure of disease risk that can be leveraged for the assignment of individuals in case-control studies or for additional downstream analyses.

Transcriptome-Wide Structural Equation Modeling of 13 Major Psychiatric Disorders for Cross-Disorder Risk and Drug Repurposing.

Importance: Psychiatric disorders display high levels of comorbidity and genetic overlap, necessitating multivariate approaches for parsing convergent and divergent psychiatric risk pathways. Identifying gene expression patterns underlying cross-disorder risk also stands to propel drug discovery and repurposing in the face of rising levels of polypharmacy. Objective: To identify gene expression patterns underlying genetic convergence and divergence across psychiatric disorders along with existing pharmacological interventions that target these genes. Design, Setting, and Participants: This genomic study applied a multivariate transcriptomic method, transcriptome-wide structural equation modeling (T-SEM), to investigate gene expression patterns associated with 5 genomic factors indexing shared risk across 13 major psychiatric disorders. Follow-up tests, including overlap with gene sets for other outcomes and phenome-wide association studies, were conducted to better characterize T-SEM results. The Broad Institute Connectivity Map Drug Repurposing Database and Drug-Gene Interaction Database public databases of drug-gene pairs were used to identify drugs that could be repurposed to target genes found to be associated with cross-disorder risk. Data were collected from database inception up to February 20, 2023. Main Outcomes and Measures: Gene expression patterns associated with genomic factors or disorder-specific risk and existing drugs that target these genes. Results: In total, T-SEM identified 466 genes whose expression was significantly associated (z ≥ 5.02) with genomic factors and 36 genes with disorder-specific effects. Most associated genes were found for a thought disorders factor, defined by bipolar disorder and schizophrenia. Several existing pharmacological interventions were identified that could be repurposed to target genes whose expression was associated with the thought disorders factor or a transdiagnostic p factor defined by all 13 disorders. Conclusions and Relevance: The findings from this study shed light on patterns of gene expression associated with genetic overlap and uniqueness across psychiatric disorders. Future versions of the multivariate drug repurposing framework outlined here have the potential to identify novel pharmacological interventions for increasingly common, comorbid psychiatric presentations.

Association Between Genetically Predicted Expression of TPMT and Azathioprine Adverse Events.

Polymorphisms thiopurine-S-methyltransferase (TPMT) and nudix hydrolase 15 (NUDT15) can increase the risk of azathioprine myelotoxicity, but little is known about other genetic factors that increase risk for azathioprine-associated side effects. PrediXcan is a gene-based association method that estimates the expression of individuals' genes and examines their correlation to specified phenotypes. As proof of concept for using PrediXcan as a tool to define the association between genetic factors and azathioprine side effects, we aimed to determine whether the genetically predicted expression of TPMT or NUDT15 was associated with leukopenia or other known side effects. In a retrospective cohort of 1364 new users of azathioprine with EHR-reported White race, we used PrediXcan to impute expression in liver tissue, tested its association with pre-specified phecodes representing known side effects (e.g., skin cancer), and completed chart review to confirm cases. Among confirmed cases, patients in the lowest tertile (i.e., lowest predicted) of TPMT expression had significantly higher odds of developing leukopenia (OR=3.30, 95%CI: 1.07-10.20, p=0.04) versus those in the highest tertile; no other side effects were significant. The results suggest that this methodology could be deployed on a larger scale to uncover associations between genetic factors and drug side effects for more personalized care.

Circulating extracellular vesicles in human cardiorenal syndrome promote renal injury in a kidney-on-chip system.

BACKGROUNDCardiorenal syndrome (CRS) - renal injury during heart failure (HF) - is linked to high morbidity. Whether circulating extracellular vesicles (EVs) and their RNA cargo directly impact its pathogenesis remains unclear.METHODSWe investigated the role of circulating EVs from patients with CRS on renal epithelial/endothelial cells using a microfluidic kidney-on-chip (KOC) model. The small RNA cargo of circulating EVs was regressed against serum creatinine to prioritize subsets of functionally relevant EV-miRNAs and their mRNA targets investigated using in silico pathway analysis, human genetics, and interrogation of expression in the KOC model and in renal tissue. The functional effects of EV-RNAs on kidney epithelial cells were experimentally validated.RESULTSRenal epithelial and endothelial cells in the KOC model exhibited uptake of EVs from patients with HF. HF-CRS EVs led to higher expression of renal injury markers (IL18, LCN2, HAVCR1) relative to non-CRS EVs. A total of 15 EV-miRNAs were associated with creatinine, targeting 1,143 gene targets specifying pathways relevant to renal injury, including TGF-ß and AMPK signaling. We observed directionally consistent changes in the expression of TGF-ß pathway members (BMP6, FST, TIMP3) in the KOC model exposed to CRS EVs, which were validated in epithelial cells treated with corresponding inhibitors and mimics of miRNAs. A similar trend was observed in renal tissue with kidney injury. Mendelian randomization suggested a role for FST in renal function.CONCLUSIONPlasma EVs in patients with CRS elicit adverse transcriptional and phenotypic responses in a KOC model by regulating biologically relevant pathways, suggesting a role for EVs in CRS.TRIAL REGISTRATIONClinicalTrials.gov NCT03345446.FUNDINGAmerican Heart Association (AHA) (SFRN16SFRN31280008); National Heart, Lung, and Blood Institute (1R35HL150807-01); National Center for Advancing Translational Sciences (UH3 TR002878); and AHA (23CDA1045944).

A phenome-wide comparative analysis of genetic discordance between obesity and type 2 diabetes.

Obesity and type 2 diabetes are causally related, yet there is considerable heterogeneity in the consequences of both conditions and the mechanisms of action are poorly defined. Here we show a genetic-driven approach defining two obesity profiles that convey highly concordant and discordant diabetogenic effects. We annotate and then compare association signals for these profiles across clinical and molecular phenotypic layers. Key differences are identified in a wide range of traits, including cardiovascular mortality, fat distribution, liver metabolism, blood pressure, specific lipid fractions and blood levels of proteins involved in extracellular matrix remodelling. We find marginal differences in abundance of Bacteroidetes and Firmicutes bacteria in the gut. Instrumental analyses reveal prominent causal roles for waist-to-hip ratio, blood pressure and cholesterol content of high-density lipoprotein particles in the development of diabetes in obesity. We prioritize 17 genes from the discordant signature that convey protection against type 2 diabetes in obesity, which may represent logical targets for precision medicine approaches.

Predicted expression of genes involved in the thiopurine metabolic pathway and azathioprine discontinuation due to myelotoxicity.

TPMT and NUDT15 variants explain less than 25% of azathioprine-associated myelotoxicity. There are 25 additional genes in the thiopurine pathway that could also contribute to azathioprine myelotoxicity. We hypothesized that among TPMT and NUDT15 normal metabolizers, a score combining the genetically predicted expression of other proteins in the thiopurine pathway would be associated with a higher risk for azathioprine discontinuation due to myelotoxicity. We conducted a retrospective cohort study of new users of azathioprine who were normal TPMT and NUDT15 metabolizers. In 1201 White patients receiving azathioprine for an inflammatory disease, we used relaxed Least Absolute Shrinkage and Selection Operator (LASSO) regression to select genes that built a score for discontinuing azathioprine due to myelotoxicity. The score incorporated the predicted expression of AOX1 and NME1. Patients in the highest score tertile had a higher risk of discontinuing azathioprine compared to those in the lowest tertile (hazard ratio [HR] = 2.15, 95% confidence interval [CI] = 1.11-4.19, p = 0.024). Results remained significant after adjusting for a propensity score, including sex, tertile of calendar year at initial dose, initial dose, age at baseline, indication, prior TPMT testing, and the first 10 principal components of the genetic data (HR = 2.11, 95% CI = 1.08-4.13, p = 0.030). We validated the results in a cohort (N = 517 non-White patients and those receiving azathioprine to prevent transplant rejection) that included all other patients receiving azathioprine (HR = 2.00, (95% CI = 1.09-3.65, p = 0.024). In conclusion, among patients who were TPMT and NUDT15 normal metabolizers, a score combining the predicted expression of AOX1 and NME1 was associated with an increased risk for discontinuing azathioprine due to myelotoxicity.

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.

Integration of genetic, transcriptomic, and clinical data provides insight into 16p11.2 and 22q11.2 CNV genes.

BACKGROUND: Deletions and duplications of the multigenic 16p11.2 and 22q11.2 copy number variant (CNV) regions are associated with brain-related disorders including schizophrenia, intellectual disability, obesity, bipolar disorder, and autism spectrum disorder (ASD). The contribution of individual CNV genes to each of these identified phenotypes is unknown, as well as the contribution of these CNV genes to other potentially subtler health implications for carriers. Hypothesizing that DNA copy number exerts most effects via impacts on RNA expression, we attempted a novel in silico fine-mapping approach in non-CNV carriers using both GWAS and biobank data. METHODS: We first asked whether gene expression level in any individual gene in the CNV region alters risk for a known CNV-associated behavioral phenotype(s). Using transcriptomic imputation, we performed association testing for CNV genes within large genotyped cohorts for schizophrenia, IQ, BMI, bipolar disorder, and ASD. Second, we used a biobank containing electronic health data to compare the medical phenome of CNV carriers to controls within 700,000 individuals in order to investigate the full spectrum of health effects of the CNVs. Third, we used genotypes for over 48,000 individuals within the biobank to perform phenome-wide association studies between imputed expressions of individual 16p11.2 and 22q11.2 genes and over 1500 health traits. RESULTS: Using large genotyped cohorts, we found individual genes within 16p11.2 associated with schizophrenia (TMEM219, INO80E, YPEL3), BMI (TMEM219, SPN, TAOK2, INO80E), and IQ (SPN), using conditional analysis to identify upregulation of INO80E as the driver of schizophrenia, and downregulation of SPN and INO80E as increasing BMI. We identified both novel and previously observed over-represented traits within the electronic health records of 16p11.2 and 22q11.2 CNV carriers. In the phenome-wide association study, we found seventeen significant gene-trait pairs, including psychosis (NPIPB11, SLX1B) and mood disorders (SCARF2), and overall enrichment of mental traits. CONCLUSIONS: Our results demonstrate how integration of genetic and clinical data aids in understanding CNV gene function and implicates pleiotropy and multigenicity in CNV biology.

Integrative genomic analyses identify susceptibility genes underlying COVID-19 hospitalization.

Despite rapid progress in characterizing the role of host genetics in SARS-Cov-2 infection, there is limited understanding of genes and pathways that contribute to COVID-19. Here, we integrate a genome-wide association study of COVID-19 hospitalization (7,885 cases and 961,804 controls from COVID-19 Host Genetics Initiative) with mRNA expression, splicing, and protein levels (n = 18,502). We identify 27 genes related to inflammation and coagulation pathways whose genetically predicted expression was associated with COVID-19 hospitalization. We functionally characterize the 27 genes using phenome- and laboratory-wide association scans in Vanderbilt Biobank (n = 85,460) and identified coagulation-related clinical symptoms, immunologic, and blood-cell-related biomarkers. We replicate these findings across trans-ethnic studies and observed consistent effects in individuals of diverse ancestral backgrounds in Vanderbilt Biobank, pan-UK Biobank, and Biobank Japan. Our study highlights and reconfirms putative causal genes impacting COVID-19 severity and symptomology through the host inflammatory response.

Soluble Urokinase Plasminogen Activator Receptor: Genetic Variation and Cardiovascular Disease Risk in Black Adults.

BACKGROUND: suPAR (Soluble urokinase plasminogen activator receptor) has emerged as an important biomarker of coagulation, inflammation, and cardiovascular disease (CVD) risk. The contribution of suPAR to CVD risk and its genetic influence in Black populations have not been evaluated. METHODS: We measured suPAR in 3492 Black adults from the prospective, community-based JHS (Jackson Heart Study). Cross-sectional associations of suPAR with lifestyle and CVD risk factors were assessed, whole-genome sequence data were used to evaluate genetic associations of suPAR, and relationships of suPAR with incident CVD outcomes and overall mortality were estimated over follow-up. RESULTS: In Cox models adjusted for traditional CVD risk factors, estimated glomerular filtration rate, and CRP (C-reactive protein), each 1-SD higher suPAR was associated with a 21% to 31% increased risk of incident coronary heart disease, heart failure, stroke, and mortality. In the genome-wide association study, 2 missense (rs399145 encoding p.Thr86Ala, rs4760 encoding p.Phe272Leu) and 2 noncoding regulatory variants (rs73935023 within an enhancer element and rs4251805 within the promoter) of PLAUR on chromosome 19 were each independently associated with suPAR and together explained 14% of suPAR phenotypic variation. The allele frequencies of each of the four suPAR-associated genetic variants differ considerably across African and European populations. We further show that PLAUR rs73935023 can alter transcriptional activity in vitro. We did not find any association between genetically determined suPAR and CVD in JHS or a larger electronic medical record-based analyses of Blacks or Whites. CONCLUSIONS: Our results demonstrate the importance of ancestry-differentiated genetic variation on suPAR levels and indicate suPAR is a CVD biomarker in Black adults.

Integrative analyses identify susceptibility genes underlying COVID-19 hospitalization.

Despite rapid progress in characterizing the role of host genetics in SARS-Cov-2 infection, there is limited understanding of genes and pathways that contribute to COVID-19. Here, we integrated a genome-wide association study of COVID-19 hospitalization (7,885 cases and 961,804 controls from COVID-19 Host Genetics Initiative) with mRNA expression, splicing, and protein levels (n=18,502). We identified 27 genes related to inflammation and coagulation pathways whose genetically predicted expression was associated with COVID-19 hospitalization. We functionally characterized the 27 genes using phenome- and laboratory-wide association scans in Vanderbilt Biobank (BioVU; n=85,460) and identified coagulation-related clinical symptoms, immunologic, and blood-cell-related biomarkers. We replicated these findings across trans-ethnic studies and observed consistent effects in individuals of diverse ancestral backgrounds in BioVU, pan-UK Biobank, and Biobank Japan. Our study highlights putative causal genes impacting COVID-19 severity and symptomology through the host inflammatory response.

The DEG/ENaC cation channel protein UNC-8 drives activity-dependent synapse removal in remodeling GABAergic neurons.

Genetic programming and neural activity drive synaptic remodeling in developing neural circuits, but the molecular components that link these pathways are poorly understood. Here we show that the C. elegans Degenerin/Epithelial Sodium Channel (DEG/ENaC) protein, UNC-8, is transcriptionally controlled to function as a trigger in an activity-dependent mechanism that removes synapses in remodeling GABAergic neurons. UNC-8 cation channel activity promotes disassembly of presynaptic domains in DD type GABA neurons, but not in VD class GABA neurons where unc-8 expression is blocked by the COUP/TF transcription factor, UNC-55. We propose that the depolarizing effect of UNC-8-dependent sodium import elevates intracellular calcium in a positive feedback loop involving the voltage-gated calcium channel UNC-2 and the calcium-activated phosphatase TAX-6/calcineurin to initiate a caspase-dependent mechanism that disassembles the presynaptic apparatus. Thus, UNC-8 serves as a link between genetic and activity-dependent pathways that function together to promote the elimination of GABA synapses in remodeling neurons.