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.

Three-dimensional mitochondria reconstructions of murine cardiac muscle changes in size across aging.

With sparse treatment options, cardiac disease remains a significant cause of death among humans. As a person ages, mitochondria breakdown and the heart becomes less efficient. Heart failure is linked to many mitochondria-associated processes, including endoplasmic reticulum stress, mitochondrial bioenergetics, insulin signaling, autophagy, and oxidative stress. The roles of key mitochondrial complexes that dictate the ultrastructure, such as the mitochondrial contact site and cristae organizing system (MICOS), in aging cardiac muscle are poorly understood. To better understand the cause of age-related alteration in mitochondrial structure in cardiac muscle, we used transmission electron microscopy (TEM) and serial block facing-scanning electron microscopy (SBF-SEM) to quantitatively analyze the three-dimensional (3-D) networks in cardiac muscle samples of male mice at aging intervals of 3 mo, 1 yr, and 2 yr. Here, we present the loss of cristae morphology, the inner folds of the mitochondria, across age. In conjunction with this, the three-dimensional (3-D) volume of mitochondria decreased. These findings mimicked observed phenotypes in murine cardiac fibroblasts with CRISPR/Cas9 knockout of Mitofilin, Chchd3, Chchd6 (some members of the MICOS complex), and Opa1, which showed poorer oxidative consumption rate and mitochondria with decreased mitochondrial length and volume. In combination, these data show the need to explore if loss of the MICOS complex in the heart may be involved in age-associated mitochondrial and cristae structural changes.NEW & NOTEWORTHY This article shows how mitochondria in murine cardiac changes, importantly elucidating age-related changes. It also is the first to show that the MICOS complex may play a role in outer membrane mitochondrial structure.

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.

The physical and psychiatric health conditions related to autism genetic scores, across genetic ancestries, sexes and age-groups in electronic health records.

BACKGROUND: Although polygenic scores (PGS) for autism have been related to various psychiatric and medical conditions, most studies to date have been conducted in research ascertained populations. We aimed to identify the psychiatric and physical conditions associated with autism PGS in a health care setting. METHODS: We computed PGS for 12,383 unrelated participants of African genetic ancestry (AF) and 65,363 unrelated participants of European genetic ancestry (EU) from Vanderbilt's de-identified biobank. Next, we performed phenome wide association studies of the autism PGS within these two genetic ancestries. RESULTS: Seven associations surpassed the Bonferroni adjusted threshold for statistical significance (p = 0.05/1374 = 3.6 × 10-5) in the EU participants, including mood disorders (OR (95%CI) = 1.08(1.05 to 1.10), p = 1.0 × 10-10), autism (OR (95%CI) = 1.34(1.24 to 1.43), p = 1.2 × 10-9), and breast cancer (OR (95%CI) = 1.09(1.05 to 1.14), 2.6 × 10-5). There was no statistical evidence for PGS-phenotype associations in the AF participants. Conditioning on the diagnosis of autism or on median body mass index (BMI) did not impact the strength of the reported associations. Although we observed some sex differences in the pattern of associations, there was no significant interaction between sex and autism PGS. Finally, the associations between autism PGS and autism diagnosis were stronger in childhood and adolescence, while the associations with mood disorders and breast cancer were stronger in adulthood. DISCUSSION: Our findings indicate that autism PGS is not only related to autism diagnosis but may also be related to adult-onset conditions, including mood disorders and some cancers. CONCLUSIONS: Our study raises the hypothesis that genes associated with autism may also increase the risk for cancers later in life. Future studies are necessary to replicate and extend our findings.

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.

Genetic Variants Associated With Hidradenitis Suppurativa.

Importance: Hidradenitis suppurativa (HS) is a common and severely morbid chronic inflammatory skin disease that is reported to be highly heritable. However, the genetic understanding of HS is insufficient, and limited genome-wide association studies (GWASs) have been performed for HS, which have not identified significant risk loci. Objective: To identify genetic variants associated with HS and to shed light on the underlying genes and genetic mechanisms. Design, Setting, and Participants: This genetic association study recruited 753 patients with HS in the HS Program for Research and Care Excellence (HS ProCARE) at the University of North Carolina Department of Dermatology from August 2018 to July 2021. A GWAS was performed for 720 patients (after quality control) with controls from the Add Health study and then meta-analyzed with 2 large biobanks, UK Biobank (247 cases) and FinnGen (673 cases). Variants at 3 loci were tested for replication in the BioVU biobank (290 cases). Data analysis was performed from September 2021 to December 2022. Main Outcomes and Measures: Main outcome measures are loci identified, with association of P < 1 × 10-8 considered significant. Results: A total of 753 patients were recruited, with 720 included in the analysis. Mean (SD) age at symptom onset was 20.3 (10.57) years and at enrollment was 35.3 (13.52) years; 360 (50.0%) patients were Black, and 575 (79.7%) were female. In a meta-analysis of the 4 studies, 2 HS-associated loci were identified and replicated, with lead variants rs10512572 (P = 2.3 × 10-11) and rs17090189 (P = 2.1 × 10-8) near the SOX9 and KLF5 genes, respectively. Variants at these loci are located in enhancer regulatory elements detected in skin tissue. Conclusions and Relevance: In this genetic association study, common variants associated with HS located near the SOX9 and KLF5 genes were associated with risk of HS. These or other nearby genes may be associated with genetic risk of disease and the development of clinical features, such as cysts, comedones, and inflammatory tunnels, that are unique to HS. New insights into disease pathogenesis related to these genes may help predict disease progression and novel treatment approaches in the future.

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.

Polygenic risk score-based phenome-wide association study identifies novel associations for Tourette syndrome.

Tourette Syndrome (TS) is a complex neurodevelopmental disorder characterized by vocal and motor tics lasting more than a year. It is highly polygenic in nature with both rare and common previously associated variants. Epidemiological studies have shown TS to be correlated with other phenotypes, but large-scale phenome wide analyses in biobank level data have not been performed to date. In this study, we used the summary statistics from the latest meta-analysis of TS to calculate the polygenic risk score (PRS) of individuals in the UK Biobank data and applied a Phenome Wide Association Study (PheWAS) approach to determine the association of disease risk with a wide range of phenotypes. A total of 57 traits were found to be significantly associated with TS polygenic risk, including multiple psychosocial factors and mental health conditions such as anxiety disorder and depression. Additional associations were observed with complex non-psychiatric disorders such as Type 2 diabetes, heart palpitations, and respiratory conditions. Cross-disorder comparisons of phenotypic associations with genetic risk for other childhood-onset disorders (e.g.: attention deficit hyperactivity disorder [ADHD], autism spectrum disorder [ASD], and obsessive-compulsive disorder [OCD]) indicated an overlap in associations between TS and these disorders. ADHD and ASD had a similar direction of effect with TS while OCD had an opposite direction of effect for all traits except mental health factors. Sex-specific PheWAS analysis identified differences in the associations with TS genetic risk between males and females. Type 2 diabetes and heart palpitations were significantly associated with TS risk in males but not in females, whereas diseases of the respiratory system were associated with TS risk in females but not in males. This analysis provides further evidence of shared genetic and phenotypic architecture of different complex disorders.

Genetic Examination of Hematological Parameters in SARS-CoV-2 Infection and COVID-19.

Background: People hospitalized with COVID-19 often exhibit hematological alterations, such as lower lymphocyte and platelet counts, which have been reported to associate with disease prognosis. It is unclear whether inter-individual variability in baseline hematological parameters prior to acute infection influences risk of SARS-CoV-2 infection and progression to severe COVID-19. Methods: We assessed the association of blood cell counts and indices 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 used summary statistics from genome-wide association studies to assess the shared genetic architecture of baseline blood cell counts and indices on COVID-19 outcomes. Results: We observed inconsistent associations between measured blood cell indices and both SARS-CoV-2 infection and COVID-19 hospitalization in UK Biobank and VUMC SD. In Mendelian randomization analyses using genetic summary statistics, no putative causal relationships were identified between COVID-19 related outcomes and hematological indices after adjusting for multiple testing. We observed overlapping genetic association signals between hematological parameters and COVID-19 traits. For example, we observed overlap between infection susceptibility-associated variants at PPP1R15A and red blood cell parameters, and between disease severity-associated variants at TYK2 and lymphocyte and platelet phenotypes. Conclusions: We did not find convincing evidence of a relationship between baseline hematological parameters and susceptibility to SARS-CoV-2 infection or COVID-19 severity, though this relationship should be re-examined as larger and better-powered genetic analyses of SARS-CoV-2 infection and severe COVID-19 become available.

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.