R PheWAS: data analysis and plotting tools for phenome-wide association studies in the R environment

Genetic insights into smoking behaviours in 10,558 men of African ancestry from continental Africa and the UK

Smoking is a leading risk factor for many of the top ten causes of death worldwide. Of the 1.3 billion smokers globally, 80% live in low- and middle-income countries, where the number of deaths due to tobacco use is expected to double in the next decade according to the World Health Organization. Genetic studies have helped to identify biological pathways for smoking behaviours, but have mostly focussed on individuals of European ancestry or living in either North America or Europe. We performed a genome-wide association study of two smoking behaviour traits in 10,558 men of African ancestry living in five African countries and the UK. Eight independent variants were associated with either smoking initiation or cessation at P-value < 5 × 10-6, four being monomorphic or rare in European populations. Gene prioritisation strategy highlighted five genes, including SEMA6D, previously described as associated with several smoking behaviour traits. These results confirm the importance of analysing underrepresented populations in genetic epidemiology, and the urgent need for larger genomic studies to boost discovery power to better understand smoking behaviours, as well as many other traits.

Stroke genetics informs drug discovery and risk prediction across ancestries

Previous genome-wide association studies (GWASs) of stroke - the second leading cause of death worldwide - were conducted predominantly in populations of European ancestry1,2. Here, in cross-ancestry GWAS meta-analyses of 110,182 patients who have had a stroke (five ancestries, 33% non-European) and 1,503,898 control individuals, we identify association signals for stroke and its subtypes at 89 (61 new) independent loci: 60 in primary inverse-variance-weighted analyses and 29 in secondary meta-regression and multitrait analyses. On the basis of internal cross-ancestry validation and an independent follow-up in 89,084 additional cases of stroke (30% non-European) and 1,013,843 control individuals, 87% of the primary stroke risk loci and 60% of the secondary stroke risk loci were replicated (P < 0.05). Effect sizes were highly correlated across ancestries. Cross-ancestry fine-mapping, in silico mutagenesis analysis3, and transcriptome-wide and proteome-wide association analyses revealed putative causal genes (such as SH3PXD2A and FURIN) and variants (such as at GRK5 and NOS3). Using a three-pronged approach4, we provide genetic evidence for putative drug effects, highlighting F11, KLKB1, PROC, GP1BA, LAMC2 and VCAM1 as possible targets, with drugs already under investigation for stroke for F11 and PROC. A polygenic score integrating cross-ancestry and ancestry-specific stroke GWASs with vascular-risk factor GWASs (integrative polygenic scores) strongly predicted ischaemic stroke in populations of European, East Asian and African ancestry5. Stroke genetic risk scores were predictive of ischaemic stroke independent of clinical risk factors in 52,600 clinical-trial participants with cardiometabolic disease. Our results provide insights to inform biology, reveal potential drug targets and derive genetic risk prediction tools across ancestries.

Multi-omics analysis reveals expression complexity and functional diversity of mouse kinome

Protein kinases are a crucial component of signaling pathways involved in a wide range of cellular responses, including growth, proliferation, differentiation, and migration. Systematic investigation of protein kinases is critical to better understand phosphorylation-mediated signaling pathways and may provide insights into the development of potential therapeutic drug targets. Here we perform a systems-level analysis of the mouse kinome by analyzing multi-omics data. We used bulk and single-cell transcriptomic data from the C57BL/6J mouse strain to define tissue- and cell-type-specific expression of protein kinases, followed by investigating variations in sequence and expression between C57BL/6J and DBA/2J strains. We then profiled a deep brain phosphoproteome from C57BL/6J and DBA/2J strains as well as their reciprocal hybrids to infer the activity of the mouse kinome. Finally, we performed phenome-wide association analysis using the BXD recombinant inbred (RI) mice (a cross between C57BL/6J and DBA/2J strains) to identify any associations between variants in protein kinases and phenotypes. Collectively, our study provides a comprehensive analysis of the mouse kinome by investigating genetic sequence variation, tissue-specific expression patterns, and associations with downstream phenotypes.

Multiancestral polygenic risk score for pediatric asthma

BACKGROUND

Asthma is the most common chronic condition in children and the third leading cause of hospitalization in pediatrics. The genome-wide association study catalog reports 140 studies with genome-wide significance. A polygenic risk score (PRS) with predictive value across ancestries has not been evaluated for this important trait.

OBJECTIVES

This study aimed to train and validate a PRS relying on genetic determinants for asthma to provide predictions for disease occurrence in pediatric cohorts of diverse ancestries.

METHODS

This study applied a Bayesian regression framework method using the Trans-National Asthma Genetic Consortium genome-wide association study summary statistics to derive a multiancestral PRS score, used one Electronic Medical Records and Genomics (eMERGE) cohort as a training set, used a second independent eMERGE cohort to validate the score, and used the UK Biobank data to replicate the findings. A phenome-wide association study was performed using the PRS to identify shared genetic etiology with other phenotypes.

RESULTS

The multiancestral asthma PRS was associated with asthma in the 2 pediatric validation datasets. Overall, the multiancestral asthma PRS has an area under the curve (AUC) of 0.70 (95% CI, 0.69-0.72) in the pediatric validation 1 and AUC of 0.66 (0.65-0.66) in the pediatric validation 2 datasets. We found significant discrimination across pediatric subcohorts of European (AUC, 95% CI, 0.60 and 0.66), African (AUC, 95% CI, 0.61 and 0.66), admixed American (AUC, 0.64 and 0.70), Southeast Asian (AUC, 0.65), and East Asian (AUC, 0.73) ancestry. Pediatric participants with the top 5% PRS had 2.80 to 5.82 increased odds of asthma compared to the bottom 5% across the training, validation 1, and validation 2 cohorts when adjusted for ancestry. Phenome-wide association study analysis confirmed the strong association of the identified PRS with asthma (odds ratio, 2.71, PFDR = 3.71 × 10-65) and related phenotypes.

CONCLUSIONS

A multiancestral PRS for asthma based on Bayesian posterior genomic effect sizes identifies increased odds of pediatric asthma.

Comorbidities in lichen planus by phenome-wide association study in two biobank population cohorts

BACKGROUND

Lichen planus (LP) is a relatively frequent mucocutaneous inflammatory disease affecting the skin, skin appendages and mucosae, including oral mucosae, and less frequently the anogenital area, conjunctivae, oesophagus or larynx.

OBJECTIVES

To estimate the association of LP, with emphasis on dermatological and gastrointestinal conditions, in two large independent population cohorts.

MATERIALS AND METHODS

We performed a phenome-wide association study (PheWAS) and examined conditions associated with LP in two unrelated cohorts, i.e. the multicentre, community-based UK Biobank (UKB: 501 381 controls; 1130 LP subjects) and the healthcare-associated Penn Medicine BioBank (PMBB; 42 702 controls; 764 LP subjects). The data were analysed in 2021. The 'PheWAS' R package was used to perform the PheWAS analyses and Bonferroni correction was used to adjust for multiple testing. Odds ratios (ORs) were adjusted for age, sex and body mass index.

RESULTS

In the UKB, PheWAS revealed 133 phenome codes (PheCodes) significantly associated with LP and most of them were confirmed in PMBB. Dermatological and digestive PheCodes were the most abundant: 29 and 34 of these disorders, respectively, were significantly overrepresented in LP individuals from both cohorts. The 29 dermatological and 12 oral disorders were often highly enriched, whereas hepatic, gastric, oesophageal and intestinal PheCodes displayed ORs in the range of 1·6-4·5. Several autoimmune disorders also exhibited OR > 5 in both cohorts.

CONCLUSIONS

PheWAS in two large unrelated cohorts identified previously unknown comorbidities and may support clinical counselling of patients with LP. What is already known about this topic? Lichen planus (LP) is known to affect the skin, skin appendages and mucosae, including oral mucosae, and less frequently the anogenital area, conjunctivae, oesophagus or larynx. What does this study add? Our data provide the most comprehensive collection of associated dermatological, digestive and autoimmune disorders to date. Our findings are expected to be useful for the evaluation and management of patients with LP.

Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores

Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real-data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, owing to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched linkage disequilibrium (LD) references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power.

Genetic contribution to cancer risk in patients with tooth loss: a genetic association study

Early-stage cancer diagnosis is critical for higher survival rates. Because early cancers can be difficult to detect, our focus is on the identification of cancer risk markers such as pleiotropic genes involved in the etiology of both craniofacial conditions and cancers. In this study we aimed to test if our previously detected association between ERN1 rs196929 marker and oral health outcomes would be detected in individuals diagnosed with cancer as well as in a subpopulation of individuals who also had one or more teeth missing due to dental caries, periodontal disease, or periapical lesions. We genotyped a total of 1,671 subjects and selected a subset of 1,421 subjects for stratified analysis of cancer types; three hundred and twelve self-reported a diagnosis of various cancer types and 1,109 reported never receiving a diagnosis of cancer. Our results showed a statistically significant association between the rs196929 in ERN1, and cancer overall in both the additive and dominant models (OR = 1.37, 95% C.I. 1.06-1.79, p = 0.014). When we stratified the analysis for each cancer type, our results show that the rs196929 ERN1 variant is associated with skin cancer (OR = 2.07, 95% C.I. 1.27-3.37, p = 0.003) and breast cancer (OR = 1.83, 95% C.I. 1.13-2.99, p = 0.013) in the subset of patients that had tooth loss. An additional nominal association between the rs196929 in ERN1 and male's reproductive system cancers (OR = 1.96, 95% C.I. 1.07-3.59, p = 0.028) was identified. We hope that our study helps guide future genetic studies on these cancers and this specific genetic variant as well as drive attention to the potential for oral health outcomes to serve as indicators for cancer risk. The early identification of genetic markers and/or oral conditions that indicate increased cancer risk could positively impact cancer outcomes and survival rates with timely implementation of preventive and diagnostic measures. In conclusion, our results suggest that the genetic variant in ERN1 (rs196929) is associated with increased risk of skin and breast cancers.

Genome-wide association analyses of common infections in a large practice-based biobank

INTRODUCTION

Infectious diseases are common causes of morbidity and mortality worldwide. Susceptibility to infection is highly heritable; however, little has been done to identify the genetic determinants underlying common infectious diseases. One GWAS was performed using 23andMe information about self-reported infections; we set out to confirm previous loci and identify new ones using medically diagnosed infections.

METHODS

We used the electronic health record (EHR)-based biobank at Vanderbilt and diagnosis codes to identify cases of 12 infectious diseases in white patients: urinary tract infection, pneumonia, chronic sinus infections, otitis media, candidiasis, streptococcal pharyngitis, herpes zoster, herpes labialis, hepatitis B, infectious mononucleosis, tuberculosis (TB) or a positive TB test, and hepatitis C. We selected controls from patients with no diagnosis code for the candidate disease and matched by year of birth, sex, and calendar year at first and last EHR visits. We conducted GWAS using SAIGE and transcriptome-wide analysis (TWAS) using S-PrediXcan. We also conducted phenome-wide association study to understand associations between identified genetic variants and clinical phenotypes.

RESULTS

We replicated three 23andMe loci (p ≤ 0.05): herpes zoster and rs7047299-A (p = 2.6 × 10^-3) and rs2808290-C (p = 9.6 × 10^-3;); otitis media and rs114947103-C (p = 0.04). We also identified 2 novel regions (p ≤ 5 × 10^-8): rs113235453-G for otitis media (p = 3.04 × 10^-8), and rs10422015-T for candidiasis (p = 3.11 × 10^-8). In TWAS, four gene-disease associations were significant: SLC30A9 for otitis media (p = 8.06 × 10^-7); LRP3 and WDR88 for candidiasis (p = 3.91 × 10^-7 and p = 1.95 × 10^-6); and AAMDC for hepatitis B (p = 1.51 × 10^-6).

CONCLUSION

We conducted GWAS and TWAS for 12 infectious diseases and identified novel genetic contributors to the susceptibility of infectious diseases.

Leveraging genomic diversity for discovery in an electronic health record linked biobank: the UCLA ATLAS Community Health Initiative

BACKGROUND

Large medical centers in urban areas, like Los Angeles, care for a diverse patient population and offer the potential to study the interplay between genetic ancestry and social determinants of health. Here, we explore the implications of genetic ancestry within the University of California, Los Angeles (UCLA) ATLAS Community Health Initiative-an ancestrally diverse biobank of genomic data linked with de-identified electronic health records (EHRs) of UCLA Health patients (N=36,736).

METHODS

We quantify the extensive continental and subcontinental genetic diversity within the ATLAS data through principal component analysis, identity-by-descent, and genetic admixture. We assess the relationship between genetically inferred ancestry (GIA) and >1500 EHR-derived phenotypes (phecodes). Finally, we demonstrate the utility of genetic data linked with EHR to perform ancestry-specific and multi-ancestry genome and phenome-wide scans across a broad set of disease phenotypes.

RESULTS

We identify 5 continental-scale GIA clusters including European American (EA), African American (AA), Hispanic Latino American (HL), South Asian American (SAA) and East Asian American (EAA) individuals and 7 subcontinental GIA clusters within the EAA GIA corresponding to Chinese American, Vietnamese American, and Japanese American individuals. Although we broadly find that self-identified race/ethnicity (SIRE) is highly correlated with GIA, we still observe marked differences between the two, emphasizing that the populations defined by these two criteria are not analogous. We find a total of 259 significant associations between continental GIA and phecodes even after accounting for individuals' SIRE, demonstrating that for some phenotypes, GIA provides information not already captured by SIRE. GWAS identifies significant associations for liver disease in the 22q13.31 locus across the HL and EAA GIA groups (HL p-value=2.32×10-16, EAA p-value=6.73×10-11). A subsequent PheWAS at the top SNP reveals significant associations with neurologic and neoplastic phenotypes specifically within the HL GIA group.

CONCLUSIONS

Overall, our results explore the interplay between SIRE and GIA within a disease context and underscore the utility of studying the genomes of diverse individuals through biobank-scale genotyping linked with EHR-based phenotyping.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Phenome of coeliac disease vs. inflammatory bowel disease

Coeliac disease (CeD) is characterized by gliadin-induced intestinal inflammation appearing in genetically susceptible individuals, such as HLA-DQ2.5 carriers. CeD, as well as other chronic intestinal disorders, such as Crohn's disease (CD) and ulcerative colitis, has been associated with increased morbidity and mortality, but the causes are unknown. We systematically analysed CeD-associated diagnoses and compared them to conditions enriched in subjects with CD/UC as well as in HLA-DQ2.5 carriers without CeD. We compared the overall and cause-specific mortality and morbidity of 3,001 patients with CeD, 2,020 with CD, 4,399 with UC and 492,200 controls in the community-based UK Biobank. Disease-specific phenotypes were assessed with the multivariable Phenome Wide Association Study (PheWAS) method. Associations were adjusted for age, sex and body mass index. All disease groups displayed higher overall mortality than controls (CD: aHR = 1.91[1.70-2.17]; UC: aHR = 1.32 [1.20-1.46]; CeD: aHR = 1.38 [1.22-1.55]). Cardiovascular and cancer-related deaths were responsible for the majority of fatalities. PheWAS analysis revealed 166 Phecodes overrepresented in all three disorders, whereas only ~ 20% of enriched Phecodes were disease specific. Seven of the 58 identified CeD-specific Phecodes were enriched in individuals homozygous for HLA-DQ2.5 without diagnosed CeD. Four out of these seven Phecodes and eight out of 19 HLA-DQ2.5 specific Phecodes were more common in homozygous HLA-DQ2.5 subjects with vs. without CeD, highlighting the interplay between genetics and diagnosis-related factors. Our study illustrates that the morbidity and mortality in CeD share similarities with CD/UC, while the CeD-restricted conditions might be driven by both inherited and acquired factors.

The All of Us Research Program: Data quality, utility, and diversity

The All of Us Research Program seeks to engage at least one million diverse participants to advance precision medicine and improve human health. We describe here the cloud-based Researcher Workbench that uses a data passport model to democratize access to analytical tools and participant information including survey, physical measurement, and electronic health record (EHR) data. We also present validation study findings for several common complex diseases to demonstrate use of this novel platform in 315,000 participants, 78% of whom are from groups historically underrepresented in biomedical research, including 49% self-reporting non-White races. Replication findings include medication usage pattern differences by race in depression and type 2 diabetes, validation of known cancer associations with smoking, and calculation of cardiovascular risk scores by reported race effects. The cloud-based Researcher Workbench represents an important advance in enabling secure access for a broad range of researchers to this large resource and analytical tools.

Exploring the clinical and genetic associations of adult weight trajectories using electronic health records in a racially diverse biobank: a phenome-wide and polygenic risk study

BACKGROUND

Weight trajectories might reflect individual health status. In this study, we aimed to examine the clinical and genetic associations of adult weight trajectories using electronic health records (EHRs) in the BioMe Biobank.

METHODS

We constructed four weight trajectories based on a-priori definitions of weight changes (5% or 10%) using annual weight in EHRs (stable weight, weight gain, weight loss, and weight cycle); the final weight dataset included 21 487 participants with 162 783 annual weight measures. To confirm accurate assignment of weight trajectories, we manually reviewed weight trajectory plots for 100 random individuals. We then did a hypothesis-free phenome-wide association study (PheWAS) to identify diseases associated with each weight trajectory. Next, we estimated the single-nucleotide polymorphism-based heritability (hSNP2) of weight trajectories using GCTA-GREML, and we did a hypothesis-driven analysis of anorexia nervosa and depression polygenic risk scores (PRS) on these weight trajectories, given both diseases are associated with weight changes. We extended our analyses to the UK Biobank to replicate findings from a patient population to a generally healthy population.

FINDINGS

We found high concordance between manually assigned weight trajectories and those assigned by the algorithm (accuracy ≥98%). Stable weight was consistently associated with lower disease risks among those passing Bonferroni-corrected p value in our PheWAS (p≤4·4 × 10-5). Additionally, we identified an association between depression and weight cycle (odds ratio [OR] 1·42, 95% CI 1·31-1·55, p≤7·7 × 10-16). The adult weight trajectories were heritable (using 5% weight change as the cutoff: hSNP2 of 2·1%, 95% CI 0·9-3·3, for stable weight; 4·1%, 1·4-6·8, for weight gain; 5·5%, 2·8-8·2, for weight loss; and 4·7%, 2·3-7·1%, for weight cycle). Anorexia nervosa PRS was positively associated with weight loss trajectory among individuals without eating disorder diagnoses (OR1SD 1·16, 95% CI 1·07-1·26, per 1 SD higher PRS, p=0·011), and the association was not attenuated by obesity PRS. No association was found between depression PRS and weight trajectories after permutation tests. All main findings were replicated in the UK Biobank (p<0·05).

INTERPRETATION

Our findings suggest the importance of considering weight from a longitudinal aspect for its association with health and highlight a crucial role of weight management during disease development and progression.

FUNDING

Klarman Family Foundation, US National Institute of Mental Health (NIMH).

Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease

More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.

Microbiome-associated human genetic variants impact phenome-wide disease risk

Human genetic variation associates with the composition of the gut microbiome, yet its influence on clinical traits remains largely unknown. We analyzed the consequences of nearly a thousand gut microbiome-associated variants (MAVs) on phenotypes reported in electronic health records from tens of thousands of individuals. We discovered and replicated associations of MAVs with neurological, metabolic, digestive, and circulatory diseases. Five significant MAVs in these categories correlate with the relative abundance of microbes down to the strain level. We also demonstrate that these relationships are independently observed and concordant with microbe by disease associations reported in case-control studies. Moreover, a selective sweep and population differentiation impacted some disease-linked MAVs. Combined, these findings establish triad relationships among the human genome, microbiome, and disease. Consequently, human genetic influences may offer opportunities for precision diagnostics of microbiome-associated diseases but also highlight the relevance of genetic background for microbiome modulation and therapeutics.

Genetic Architecture of Plasma Alpha-Aminoadipic Acid Reveals a Relationship With High-Density Lipoprotein Cholesterol

Background Elevated plasma levels of alpha-aminoadipic acid (2-AAA) have been associated with the development of type 2 diabetes and atherosclerosis. However, the nature of the association remains unknown. Methods and Results We identified genetic determinants of plasma 2-AAA through meta-analysis of genome-wide association study data in 5456 individuals of European, African, and Asian ancestry from the Framingham Heart Study, Diabetes Prevention Program, Jackson Heart Study, and Shanghai Women's and Men's Health Studies. No single nucleotide polymorphisms reached genome-wide significance across all samples. However, the top associations from the meta-analysis included single-nucleotide polymorphisms in the known 2-AAA pathway gene DHTKD1, and single-nucleotide polymorphisms in genes involved in mitochondrial respiration (NDUFS4) and macrophage function (MSR1). We used a Mendelian randomization instrumental variable approach to evaluate relationships between 2-AAA and cardiometabolic phenotypes in large disease genome-wide association studies. Mendelian randomization identified a suggestive inverse association between increased 2-AAA and lower high-density lipoprotein cholesterol (P=0.005). We further characterized the genetically predicted relationship through measurement of plasma 2-AAA and high-density lipoprotein cholesterol in 2 separate samples of individuals with and without cardiometabolic disease (N=98), and confirmed a significant negative correlation between 2-AAA and high-density lipoprotein (rs=-0.53, P<0.0001). Conclusions 2-AAA levels in plasma may be regulated, in part, by common variants in genes involved in mitochondrial and macrophage function. Elevated plasma 2-AAA associates with reduced levels of high-density lipoprotein cholesterol. Further mechanistic studies are required to probe this as a possible mechanism linking 2-AAA to future cardiometabolic risk.

Impacts of ADH1B rs1229984 and ALDH2 rs671 polymorphisms on risks of alcohol-related disorder and cancer

BACKGROUND

ADH1B rs1229984 and ALDH2 rs671 are the specifically prevalent functional variants in the East Asians. These variants, which result in a dramatic change in enzyme activity, are highly associated with alcohol-related disorders and cancer. Previous studies focusing on the additive and synergic effects of the variants are few and inconsistent. The aim of the research was to evaluate the associations of ADH1B rs1229984 and ALDH2 rs671 with the risks of alcohol-related disorder and cancer.

METHODS

This cohort study enrolled 42,665 participants from the Taiwan Precision Medicine Initiative database, including 19,522 and 20,534, ADH1B and ALDH2 carriers, respectively. The associations between the two variants and cancer risk were analyzed by univariable and multivariable logistic regression.

RESULTS

Compared with the noncarriers, the ADH1B rs1229984 variant had a stronger effect on alcohol-related disorders and was related to an increased risk of alcohol-related cancers. The CC genotype of ADH1B rs1229984 was significantly associated with cancer of the larynx, pharynx, and nasal cavities [odds ratio (OR) = 1.56, p = 0.0009], cancer of the pancreas (OR = 1.66, p = 0.018), and cancer of the esophagus (OR = 4.10, p < 0.001). Participants who carried the rs1229984 TC/CC and rs671 GG genotypes were at higher risk of esophageal cancer (OR = 3.02, p < 0.001). The risk of esophageal cancer was increased by 381% (OR = 4.81, p < 0.001) in those carrying the rs1229984 TC/CC and rs671 GA/AA genotypes.

CONCLUSION

rs1229984 and rs671 are common and functionally important genetic variants in the Taiwanese population. Our findings provide strong evidence of additive and synergic risks of ADH1B and ALDH2 variants for alcohol-related disorders and cancer. The results suggested that are reduction in alcohol consumption should be advised as a preventive measure for high-risk patients carrying ADH1B rs1229984 C or the ALDH2 rs671 A allele.

Comprehensive genetic analysis of the human lipidome identifies loci associated with lipid homeostasis with links to coronary artery disease

We integrated lipidomics and genomics to unravel the genetic architecture of lipid metabolism and identify genetic variants associated with lipid species putatively in the mechanistic pathway for coronary artery disease (CAD). We quantified 596 lipid species in serum from 4,492 individuals from the Busselton Health Study. The discovery GWAS identified 3,361 independent lipid-loci associations, involving 667 genomic regions (479 previously unreported), with validation in two independent cohorts. A meta-analysis revealed an additional 70 independent genomic regions associated with lipid species. We identified 134 lipid endophenotypes for CAD associated with 186 genomic loci. Associations between independent lipid-loci with coronary atherosclerosis were assessed in ∼456,000 individuals from the UK Biobank. Of the 53 lipid-loci that showed evidence of association (P < 1 × 10-3), 43 loci were associated with at least one lipid endophenotype. These findings illustrate the value of integrative biology to investigate the aetiology of atherosclerosis and CAD, with implications for other complex diseases.

Analysis of MRI-derived spleen iron in the UK Biobank identifies genetic variation linked to iron homeostasis and hemolysis

The spleen plays a key role in iron homeostasis. It is the largest filter of the blood and performs iron reuptake from old or damaged erythrocytes. Despite this role, spleen iron concentration has not been measured in a large, population-based cohort. In this study, we quantify spleen iron in 41,764 participants of the UK Biobank by using magnetic resonance imaging and provide a reference range for spleen iron in an unselected population. Through genome-wide association study, we identify associations between spleen iron and regulatory variation at two hereditary spherocytosis genes, ANK1 and SPTA1. Spherocytosis-causing coding mutations in these genes are associated with lower reticulocyte volume and increased reticulocyte percentage, while these common alleles are associated with increased expression of ANK1 and SPTA1 in blood and with larger reticulocyte volume and reduced reticulocyte percentage. As genetic modifiers, these common alleles may explain mild spherocytosis phenotypes that have been observed clinically. Our genetic study also identifies a signal that co-localizes with a splicing quantitative trait locus for MS4A7, and we show this gene is abundantly expressed in the spleen and in macrophages. The combination of deep learning and efficient image processing enables non-invasive measurement of spleen iron and, in turn, characterization of genetic factors related to the lytic phase of the erythrocyte life cycle and iron reuptake in the spleen.

A multi-population phenome-wide association study of genetically-predicted height in the Million Veteran Program

BACKGROUND

Height has been associated with many clinical traits but whether such associations are causal versus secondary to confounding remains unclear in many cases. To systematically examine this question, we performed a Mendelian Randomization-Phenome-wide association study (MR-PheWAS) using clinical and genetic data from a national healthcare system biobank.

METHODS AND FINDINGS

Analyses were performed using data from the US Veterans Affairs (VA) Million Veteran Program in non-Hispanic White (EA, n = 222,300) and non-Hispanic Black (AA, n = 58,151) adults in the US. We estimated height genetic risk based on 3290 height-associated variants from a recent European-ancestry genome-wide meta-analysis. We compared associations of measured and genetically-predicted height with phenome-wide traits derived from the VA electronic health record, adjusting for age, sex, and genetic principal components. We found 345 clinical traits associated with measured height in EA and an additional 17 in AA. Of these, 127 were associated with genetically-predicted height at phenome-wide significance in EA and 2 in AA. These associations were largely independent from body mass index. We confirmed several previously described MR associations between height and cardiovascular disease traits such as hypertension, hyperlipidemia, coronary heart disease (CHD), and atrial fibrillation, and further uncovered MR associations with venous circulatory disorders and peripheral neuropathy in the presence and absence of diabetes. As a number of traits associated with genetically-predicted height frequently co-occur with CHD, we evaluated effect modification by CHD status of genetically-predicted height associations with risk factors for and complications of CHD. We found modification of effects of MR associations by CHD status for atrial fibrillation/flutter but not for hypertension, hyperlipidemia, or venous circulatory disorders.

CONCLUSIONS

We conclude that height may be an unrecognized but biologically plausible risk factor for several common conditions in adults. However, more studies are needed to reliably exclude horizontal pleiotropy as a driving force behind at least some of the MR associations observed in this study.

Genome-wide and phenome-wide analysis of ideal cardiovascular health in the VA Million Veteran Program

BACKGROUND

Genetic studies may help identify causal pathways; therefore, we sought to identify genetic determinants of ideal CVH and their association with CVD outcomes in the multi-population Veteran Administration Million Veteran Program.

METHODS

An ideal health score (IHS) was calculated from 3 clinical factors (blood pressure, total cholesterol, and blood glucose levels) and 3 behavioral factors (smoking status, physical activity, and BMI), ascertained at baseline. Multi-population genome-wide association study (GWAS) was performed on IHS and binary ideal health using linear and logistic regression, respectively. Using the genome-wide significant SNPs from the IHS GWAS, we created a weighted IHS polygenic risk score (PRSIHS) which was used (i) to conduct a phenome-wide association study (PheWAS) of associations between PRSIHS and ICD-9 phenotypes and (ii) to further test for associations with mortality and selected CVD outcomes using logistic and Cox regression and, as an instrumental variable, in Mendelian Randomization.

RESULTS

The discovery and replication cohorts consisted of 142,404 (119,129 European American (EUR); 16,495 African American (AFR)), and 45,766 (37,646 EUR; 5,366 AFR) participants, respectively. The mean age was 65.8 years (SD = 11.2) and 92.7% were male. Overall, 4.2% exhibited ideal CVH based on the clinical and behavioral factors. In the multi-population meta-analysis, variants at 17 loci were associated with IHS and each had known GWAS associations with multiple components of the IHS. PheWAS analysis in 456,026 participants showed that increased PRSIHS was associated with a lower odds ratio for many CVD outcomes and risk factors. Both IHS and PRSIHS measures of ideal CVH were associated with significantly less CVD outcomes and CVD mortality.

CONCLUSION

A set of high interest genetic variants contribute to the presence of ideal CVH in a multi-ethnic cohort of US Veterans. Genetically influenced ideal CVH is associated with lower odds of CVD outcomes and mortality.

Psychiatric manifestations of rare variation in medically actionable genes: a PheWAS approach

BACKGROUND

As genomic sequencing moves closer to clinical implementation, there has been an increasing acceptance of returning incidental findings to research participants and patients for mutations in highly penetrant, medically actionable genes. A curated list of genes has been recommended by the American College of Medical Genetics and Genomics (ACMG) for return of incidental findings. However, the pleiotropic effects of these genes are not fully known. Such effects could complicate genetic counseling when returning incidental findings. In particular, there has been no systematic evaluation of psychiatric manifestations associated with rare variation in these genes.

RESULTS

Here, we leveraged a targeted sequence panel and real-world electronic health records from the eMERGE network to assess the burden of rare variation in the ACMG-56 genes and two psychiatric-associated genes (CACNA1C  and TCF4) across common mental health conditions in 15,181 individuals of European descent. As a positive control, we showed that this approach replicated the established association between rare mutations in LDLR and hypercholesterolemia with no visible inflation from population stratification. However, we did not identify any genes significantly enriched with rare deleterious variants that confer risk for common psychiatric disorders after correction for multiple testing. Suggestive associations were observed between depression and rare coding variation in PTEN (P = 1.5 × 10-4), LDLR (P = 3.6 × 10-4), and CACNA1S (P = 5.8 × 10-4). We also observed nominal associations between rare variants in KCNQ1 and substance use disorders (P = 2.4 × 10-4), and APOB and tobacco use disorder (P = 1.1 × 10-3).

CONCLUSIONS

Our results do not support an association between psychiatric disorders and incidental findings in medically actionable gene mutations, but power was limited with the available sample sizes. Given the phenotypic and genetic complexity of psychiatric phenotypes, future work will require a much larger sequencing dataset to determine whether incidental findings in these genes have implications for risk of psychopathology.

GWAS of Hematuria

BACKGROUND AND OBJECTIVES

Glomerular hematuria has varied causes but can have a genetic basis, including Alport syndrome and IgA nephropathy.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

We used summary statistics to identify genetic variants associated with hematuria in White British UK Biobank participants. Individuals with glomerular hematuria were enriched by excluding participants with genitourinary conditions. A strongly associated locus on chromosome 2 (COL4A4-COL4A3) was identified. The region was reimputed using the Trans-Omics for Precision Medicine Program followed by sequential rounds of regional conditional analysis, conditioning on previous genetic signals. Similarly, we applied conditional analysis to identify independent variants in the MHC region on chromosome 6 using imputed HLA haplotypes.

RESULTS

In total, 16,866 hematuria cases and 391,420 controls were included. Cases had higher urinary albumin-creatinine compared with controls (women: 13.01 mg/g [8.05-21.33] versus 12.12 mg/g [7.61-19.29]; P<0.001; men: 8.85 mg/g [5.66-16.19] versus 7.52 mg/g [5.04-12.39]; P<0.001) and lower eGFR (women: 88±14 versus 90±13 ml/min per 1.72 m²; P<0.001; men: 87±15 versus 90±13 ml/min per 1.72 m²; P<0.001), supporting enrichment of glomerular hematuria. Variants at six loci (PDPN, COL4A4-COL4A3, HLA-B, SORL1, PLLP, and TGFB1) met genome-wide significance (P<5E-8). At chromosome 2, COL4A4 p.Ser969X (rs35138315; minor allele frequency=0.00035; P<7.95E-35; odds ratio, 87.3; 95% confidence interval, 47.9 to 159.0) had the most significant association, and two variants in the locus remained associated with hematuria after conditioning for this variant: COL4A3 p.Gly695Arg (rs200287952; minor allele frequency=0.00021; P<2.16E-7; odds ratio, 45.5; 95% confidence interval, 11.8 to 168.0) and a common COL4A4 intron 25 variant (not previously reported; rs58261427; minor allele frequency=0.214; P<2.00E-9; odds ratio, 1.09; 95% confidence interval, 1.06 to 1.12). Of the HLA haplotypes, HLA-B (*0801; minor allele frequency=0.14; P<4.41E-24; odds ratio, 0.84; 95% confidence interval, 0.82 to 0.88) displayed the most statistically significant association. For remaining loci, we identified three novel associations, which were replicated in the deCODE dataset for dipstick hematuria (nearest genes: PDPN, SORL1, and PLLP).

CONCLUSIONS

Our study identifies six loci associated with hematuria, including independent variants in COL4A4-COL4A3 and HLA-B. Additionally, three novel loci are reported, including an association with an intronic variant in PDPN expressed in the podocyte.

PODCAST

This article contains a podcast at https://www.asn-online.org/media/podcast/CJASN/2022_04_26_CJN13711021.mp3.

Genetic analysis of over half a million people characterises C-reactive protein loci

Chronic low-grade inflammation is linked to a multitude of chronic diseases. We report the largest genome-wide association study (GWAS) on C-reactive protein (CRP), a marker of systemic inflammation, in UK Biobank participants (N = 427,367, European descent) and the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (total N = 575,531 European descent). We identify 266 independent loci, of which 211 are not previously reported. Gene-set analysis highlighted 42 gene sets associated with CRP levels (p ≤ 3.2 ×10-6) and tissue expression analysis indicated a strong association of CRP related genes with liver and whole blood gene expression. Phenome-wide association study identified 27 clinical outcomes associated with genetically determined CRP and subsequent Mendelian randomisation analyses supported a causal association with schizophrenia, chronic airway obstruction and prostate cancer. Our findings identified genetic loci and functional properties of chronic low-grade inflammation and provided evidence for causal associations with a range of diseases.

Large-scale real-world data analysis identifies comorbidity patterns in schizophrenia

Schizophrenia affects >3.2 million people in the USA. However, its comorbidity patterns have not been systematically characterized in real-world populations. To address this gap, we conducted an observational study using a cohort of 86 million patients in a nationwide health insurance dataset. We identified participants with schizophrenia and those without schizophrenia matched by age, sex, and the first three digits of zip code. For each phenotype encoded in phecodes, we compared their prevalence in schizophrenia patients and the matched non-schizophrenic participants, and we performed subgroup analyses stratified by age and sex. Results show that anxiety, posttraumatic stress disorder, and substance abuse commonly occur in adolescents and young adults prior to schizophrenia diagnoses. Patients aged 60 and above are at higher risks of developing delirium, alcoholism, dementia, pelvic fracture, and osteomyelitis than their matched controls. Type 2 diabetes, sleep apnea, and eating disorders were more prevalent in women prior to schizophrenia diagnosis, whereas acute renal failure, rhabdomyolysis, and developmental delays were found at higher rates in men. Anxiety and obesity are more commonly seen in patients with schizoaffective disorders compared to patients with other types of schizophrenia. Leveraging a large-scale insurance claims dataset, this study identified less-known comorbidity patterns of schizophrenia and confirmed known ones. These comorbidity profiles can guide clinicians and researchers to take heed of early signs of co-occurring diseases.

pyPheWAS: A Phenome-Disease Association Tool for Electronic Medical Record Analysis

Along with the increasing availability of electronic medical record (EMR) data, phenome-wide association studies (PheWAS) and phenome-disease association studies (PheDAS) have become a prominent, first-line method of analysis for uncovering the secrets of EMR. Despite this recent growth, there is a lack of approachable software tools for conducting these analyses on large-scale EMR cohorts. In this article, we introduce pyPheWAS, an open-source python package for conducting PheDAS and related analyses. This toolkit includes 1) data preparation, such as cohort censoring and age-matching; 2) traditional PheDAS analysis of ICD-9 and ICD-10 billing codes; 3) PheDAS analysis applied to a novel EMR phenotype mapping: current procedural terminology (CPT) codes; and 4) novelty analysis of significant disease-phenotype associations found through PheDAS. The pyPheWAS toolkit is approachable and comprehensive, encapsulating data prep through result visualization all within a simple command-line interface. The toolkit is designed for the ever-growing scale of available EMR data, with the ability to analyze cohorts of 100,000 + patients in less than 2 h. Through a case study of Down Syndrome and other intellectual developmental disabilities, we demonstrate the ability of pyPheWAS to discover both known and potentially novel disease-phenotype associations across different experiment designs and disease groups. The software and user documentation are available in open source at https://github.com/MASILab/pyPheWAS .

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (nrs495828 = 53 and nrs505922 = 59); strongest association with venous embolism, odds ratio (ORrs495828 1.33 (p = 1.32 x 10-199), and thrombosis ORrs505922 1.33, p = 2.2 x10-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10-191; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10-12. The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10-23, lupus OR 0.84, p = 3.97 x 10-06. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10-13 among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.

Phenome-Wide Association Study of Polygenic Risk Score for Alzheimer's Disease in Electronic Health Records

Alzheimer's disease (AD) is the most common form of dementia and a growing public health burden in the United States. Significant progress has been made in identifying genetic risk for AD, but limited studies have investigated how AD genetic risk may be associated with other disease conditions in an unbiased fashion. In this study, we conducted a phenome-wide association study (PheWAS) by genetic ancestry groups within a large academic health system using the polygenic risk score (PRS) for AD. PRS was calculated using LDpred2 with genome-wide association study (GWAS) summary statistics. Phenotypes were extracted from electronic health record (EHR) diagnosis codes and mapped to more clinically meaningful phecodes. Logistic regression with Firth's bias correction was used for PRS phenotype analyses. Mendelian randomization was used to examine causality in significant PheWAS associations. Our results showed a strong association between AD PRS and AD phenotype in European ancestry (OR = 1.26, 95% CI: 1.13, 1.40). Among a total of 1,515 PheWAS tests within the European sample, we observed strong associations of AD PRS with AD and related phenotypes, which include mild cognitive impairment (MCI), memory loss, and dementias. We observed a phenome-wide significant association between AD PRS and gouty arthropathy (OR = 0.90, adjusted p = 0.05). Further causal inference tests with Mendelian randomization showed that gout was not causally associated with AD. We concluded that genetic predisposition of AD was negatively associated with gout, but gout was not a causal risk factor for AD. Our study evaluated AD PRS in a real-world EHR setting and provided evidence that AD PRS may help to identify individuals who are genetically at risk of AD and other related phenotypes. We identified non-neurodegenerative diseases associated with AD PRS, which is essential to understand the genetic architecture of AD and potential side effects of drugs targeting genetic risk factors of AD. Together, these findings expand our understanding of AD genetic and clinical risk factors, which provide a framework for continued research in aging with the growing number of real-world EHR linked with genetic data.

Association of Telomere Length With Risk of Disease and Mortality

IMPORTANCE

Telomeres protect DNA from damage. Because they shorten with each mitotic cycle, leukocyte telomere length (LTL) serves as a mitotic clock. Reduced LTL has been associated with multiple human disorders.

OBJECTIVE

To determine the association between LTL and overall as well as disease-specific mortality and morbidity.

DESIGN, SETTING, AND PARTICIPANTS

This multicenter, community-based cohort study conducted from March 2006 to December 2010 included longitudinal follow-up (mean [SD], 12 [2] years) for 472 432 English participants from the United Kingdom Biobank (UK Biobank) and analyzed morbidity and mortality. The data were analyzed in 2021.

MAIN OUTCOMES AND MEASURES

Hazard ratios (HRs) and odds ratios for mortality and morbidity associated with a standard deviation change in LTL, adjusted for age, sex, body mass index (calculated as weight in kilograms divided by height in meters squared), and ethnicity.

RESULTS

This study included a total of 472 432 English participants, of whom 54% were women (mean age, 57 years). Reduced LTL was associated with increased overall (HR, 1.08; 95% CI, 1.07-1.09), cardiovascular (HR, 1.09; 95% CI, 1.06-1.12), respiratory (HR, 1.40; 95% CI, 1.34-1.45), digestive (HR, 1.26; 95% CI, 1.19-1.33), musculoskeletal (HR, 1.51; 95% CI, 1.35-1.92), and COVID-19 (HR, 1.15; 95% CI, 1.07-1.23) mortality, but not cancer-related mortality. A total of 214 disorders were significantly overrepresented and 37 underrepresented in participants with shorter LTL. Respiratory (11%), digestive/liver-related (14%), circulatory (18%), and musculoskeletal conditions (6%), together with infections (5%), accounted for most positive associations, whereas (benign) neoplasms and endocrinologic/metabolic disorders were the most underrepresented entities. Malignant tumors, esophageal cancer, and lymphoid and myeloid leukemia were significantly more common in participants with shorter LTL, whereas brain cancer and melanoma were less prevalent. While smoking and alcohol consumption were associated with shorter LTL, additional adjustment for both factors, as well as cognitive function/major comorbid conditions, did not significantly alter the results.

CONCLUSIONS AND RELEVANCE

This cohort study found that shorter LTL was associated with a small risk increase of overall mortality, but a higher risk of mortality was associated with specific organs and diseases.

Association of Preinjury Medical Diagnoses With Pediatric Persistent Postconcussion Symptoms in Electronic Health Records

OBJECTIVE

To identify risk factors and generate hypotheses for pediatric persistent postconcussion symptoms (PPCS).

SETTING

A regional healthcare system in the Southeastern United States.

PARTICIPANTS

An electronic health record-based algorithm was developed and validated to identify PPCS cases and controls from an institutional database of more than 2.8 million patients. PPCS cases (n = 274) were patients aged 5 to 18 years with PPCS-related diagnostic codes or with PPCS key words identified by natural language processing of clinical notes. Age, sex, and year of index event-matched controls (n = 1096) were patients with mild traumatic brain injury codes only. Patients with moderate or severe traumatic brain injury were excluded. All patients used our healthcare system at least 3 times 180 days before their injury.

DESIGN

Case-control study.

MAIN MEASURES

The outcome was algorithmic classification of PPCS. Exposures were all preinjury medical diagnoses assigned at least 180 days before the injury.

RESULTS

Cases and controls both had a mean of more than 9 years of healthcare system use preinjury. Of 221 preinjury medical diagnoses, headache disorder was associated with PPCS after accounting for multiple testing (odds ratio [OR] = 2.9; 95% confidence interval [CI]: 1.6-5.0; P = 2.1e-4). Six diagnoses were associated with PPCS at a suggestive threshold for statistical significance (false discovery rate P < .10): gastritis/duodenitis (OR = 2.8; 95% CI: 1.6-5.1; P = 5.0e-4), sleep disorders (OR = 2.3; 95% CI: 1.4-3.7; P = 7.4e-4), abdominal pain (OR = 1.6; 95% CI: 1.2-2.2; P = 9.2e-4), chronic sinusitis (OR = 2.8; 95% CI: 1.5-5.2; P = 1.3e-3), congenital anomalies of the skin (OR = 2.9; 95% CI: 1.5-5.5; P = 1.9e-3), and chronic pharyngitis/nasopharyngitis (OR = 2.4; 95% CI: 1.4-4.3; P = 2.5e-3).

CONCLUSIONS

These results support the strong association of preinjury headache disorders with PPCS. An association of PPCS with prior gastritis/duodenitis, sinusitis, and pharyngitis/nasopharyngitis suggests a role for chronic inflammation in PPCS pathophysiology and risk, although results could equally be attributable to a higher likelihood of somatization among PPCS cases. Identified risk factors should be investigated further and potentially considered during the management of pediatric mild traumatic brain injury cases.

Identifying Potential Therapeutic Applications and Diagnostic Harms of Increased Bilirubin Concentrations: A Clinical and Genetic Approach

Bilirubin has antioxidant and anti-inflammatory properties in vitro and in animal studies and protects against inflammatory, cardiovascular, and other diseases in observational studies; therefore, bilirubin has potential as a therapeutic agent. However, observational studies could be confounded by many factors. We used a genetic (n = 61,281) and clinical (n = 234,670) approach to define the association between bilirubin and 19 conditions with a putative protective signal in observational studies. We also tested if individuals with genetically higher bilirubin levels underwent more diagnostic tests. We used a common variant in UGT1A1 (rs6742078) associated with an 26% increase in bilirubin levels in the genetic studies. Carriers of the variant had higher bilirubin levels (P = 2.2 × 10-16 ) but there was no significant association with any of the 19 conditions. In a phenome-wide association study (pheWAS) to seek undiscovered genetic associations, the only significant finding was increased risk of "jaundice-not of newborn." Carriers of the variant allele were more likely to undergo an abdominal ultrasound (odds ratio = 1.04, [1.00-1.08], P = 0.03). In contrast, clinically measured bilirubin levels were significantly associated with 15 of the 19 conditions (P < 0.003) and with 431 clinical diagnoses in the pheWAS (P < 1 × 10-5 adjusted for sex, age, and follow-up). With additional adjustment for smoking and body mass index, 7 of 19 conditions and 260 pheWAS diagnoses remained significantly associated with bilirubin levels. In conclusion, bilirubin does not protect against inflammatory or other diseases using a genetic approach; the many putative beneficial associations reported clinically are likely due to confounding.

Health effects of high serum calcium levels: Updated phenome-wide Mendelian randomisation investigation and review of Mendelian randomisation studies

Background

Calcium plays a role in a wide range of biological functions. Here we conducted a phenome-wide Mendelian randomisation (MR-PheWAS) analysis and a systematic review for MR studies to comprehensively investigate the health effects of serum calcium.

Methods

One-hundred and thirty genetic variants strongly associated with serum calcium levels were used as instrumental variables. A phenome-wide association analysis (PheWAS) was conducted to examine the associations of genetically predicted serum calcium with 1473 distinct phenotypes in the UK Biobank including 339,197 individuals. Observed associations in PheWAS were further tested for replication in two-sample MR replication analysis. A systematic review for MR studies on serum calcium was performed to synthesize the published evidence and compare with the current MR-PheWAS findings.

Findings

Higher genetically predicted calcium levels were associated with decreased risk of 5 diseases in dermatologic and musculoskeletal systems and increased risk of 17 diseases in circulatory, digestive, endocrine, genitourinary and immune systems. Eight associations were replicated in two-sample MR analysis. These included decreased risk of osteoarthritis and increased risk of coronary artery disease, myocardial infarction, coronary atherosclerosis, hyperparathyroidism, disorder of parathyroid gland, gout, and calculus of kidney and ureter with increased serum calcium. Systematic review of 25 MR studies provided supporting evidence on five out of the eight disease outcomes, while the increased risk of gout, hyperparathyroidism and disorder of parathyroid gland were novel findings.

Interpretation

This study found wide-ranged health effects of high serum calcium, which suggests that the benefits and adversities of strategies promoting calcium intake should be assessed.

Funding

ET is supported by a CRUK Career Development Fellowship (C31250/A22804). XL is supported by the Natural Science Fund for Distinguished Young Scholars of Zhejiang Province. SCL acknowledges research funding from the Swedish Heart Lung Foundation (Hjärt-Lungfonden, 20210351), the Swedish Research Council (Vetenskapsrådet, 2019-00977), and the Swedish Cancer Society (Cancerfonden).

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.

Associations between genetic loci, environment factors and mental disorders: a genome-wide survival analysis using the UK Biobank data

It is well-accepted that both environment and genetic factors contribute to the development of mental disorders (MD). However, few genetic studies used time-to-event data analysis to identify the susceptibility genetic variants associated with MD and explore the role of environment factors in these associations. In order to detect novel genetic loci associated with MD based on the time-to-event data and identify the role of environmental factors in them, this study recruited 376,806 participants from the UK Biobank cohort. The MD outcomes (including overall MD status, anxiety, depression and substance use disorders (SUD)) were defined based on in-patient hospital, self-reported and death registry data collected in the UK Biobank. SPACOX approach was used to identify the susceptibility loci for MD using the time-to-event data of the UK Biobank cohort. And then we estimated the associations between identified candidate loci, fourteen environment factors and MD through a phenome-wide association study and mediation analysis. SPACOX identified multiple candidate loci for overall MD status, depression and SUD, such as rs139813674 (P value = 8.39 × 10^-9, ZNF684) for overall MD status, rs7231178 (DCC, P value = 2.11 × 10^-9) for depression, and rs10228494 (FOXP2, P value = 6.58 × 10^-10) for SUD. Multiple environment factors could influence the associations between identified loci and MD, such as confide in others and felt hated. Our study identified novel candidate loci for MD, highlighting the strength of time-to-event data based genetic association studies. We also observed that multiple environment factors could influence the association between susceptibility loci and MD.

Portability of 245 polygenic scores when derived from the UK Biobank and applied to 9 ancestry groups from the same cohort

The low portability of polygenic scores (PGSs) across global populations is a major concern that must be addressed before PGSs can be used for everyone in the clinic. Indeed, prediction accuracy has been shown to decay as a function of the genetic distance between the training and test cohorts. However, such cohorts differ not only in their genetic distance but also in their geographical distance and their data collection and assaying, conflating multiple factors. In this study, we examine the extent to which PGSs are transferable between ancestries by deriving polygenic scores for 245 curated traits from the UK Biobank data and applying them in nine ancestry groups from the same cohort. By restricting both training and testing to the UK Biobank data, we reduce the risk of environmental and genotyping confounding from using different cohorts. We define the nine ancestry groups at a sub-continental level, based on a simple, robust, and effective method that we introduce here. We then apply two different predictive methods to derive polygenic scores for all 245 phenotypes and show a systematic and dramatic reduction in portability of PGSs trained using Northwestern European individuals and applied to nine ancestry groups. These analyses demonstrate that prediction already drops off within European ancestries and reduces globally in proportion to genetic distance. Altogether, our study provides unique and robust insights into the PGS portability problem.

Maturation and application of phenome-wide association studies

In the past 10 years since its introduction, phenome-wide association studies (PheWAS) have uncovered novel genotype-phenotype relationships. Along the way, PheWAS have evolved in many aspects as a study design with the expanded availability of large data repositories with genome-wide data linked to detailed phenotypic data. Advancement in methods, including algorithms, software, and publicly available integrated resources, makes it feasible to more fully realize the potential of PheWAS, overcoming the previous computational and analytical limitations. We review here the most recent improvements and notable applications of PheWAS since the second half of the decade from its inception. We also note the challenges that remain embedded along the entire PheWAS analytical pipeline that necessitate further development of tools and resources to further advance the understanding of the complex genetic architecture underlying human diseases and traits.

Using phenome-wide association studies and the SF-12 quality of life metric to identify profound consequences of adverse childhood experiences on adult mental and physical health in a Northern Nevadan population

In this research, we examine and identify the implications of Adverse Childhood Experiences (ACEs) on a range of health outcomes, with particular focus on a number of mental health disorders. Many previous studies observed that traumatic childhood events are linked to long-term adult diseases using the standard Adverse Childhood Experience Questionnaire. The study cohort was derived from the Healthy Nevada Project, a volunteer-based population health study in which each adult participant is invited to take a retrospective questionnaire that includes the Adverse Childhood Experience Questionnaire, the 12-item Short Form Survey measuring quality of life, and self-reported incidence of nine mental disorders. Using participant's cross-referenced electronic health records, a phenome-wide association analysis of 1,703 phenotypes and the incidence of ACEs examined links between traumatic events in childhood and adult disease. These analyses showed that many mental disorders were significantly associated with ACEs in a dose-response manner. Similarly, a dose response between ACEs and obesity, chronic pain, migraine, and other physical phenotypes was identified. An examination of the prevalence of self-reported mental disorders and incidence of ACEs showed a positive relationship. Furthermore, participants with less adverse childhood events experienced a higher quality of life, both physically and mentally. The whole-phenotype approach confirms that ACEs are linked with many negative adult physical and mental health outcomes. With the nationwide prevalence of ACEs as high as 67%, these findings suggest a need for new public health resources: ACE-specific interventions and early childhood screenings.

Systems Approach to Integrating Preclinical Apolipoprotein E-Knockout Investigations Reveals Novel Etiologic Pathways and Master Atherosclerosis Network in Humans

OBJECTIVE

Animal models of atherosclerosis are used extensively to interrogate molecular mechanisms in serial fashion. We tested whether a novel systems biology approach to integration of preclinical data identifies novel pathways and regulators in human disease. Approach and Results: Of 716 articles published in ATVB from 1995 to 2019 using the apolipoprotein E knockout mouse to study atherosclerosis, data were extracted from 360 unique studies in which a gene was experimentally perturbed to impact plaque size or composition and analyzed using Ingenuity Pathway Analysis software. TREM1 (triggering receptor expressed on myeloid cells) signaling and LXR/RXR (liver X receptor/retinoid X receptor) activation were identified as the top atherosclerosis-associated pathways in mice (both P<1.93×10-4, TREM1 implicated early and LXR/RXR in late atherogenesis). The top upstream regulatory network in mice (sc-58125, a COX2 inhibitor) linked 64.0% of the genes into a single network. The pathways and networks identified in mice were interrogated by testing for associations between the genetically predicted gene expression of each mouse pathway-identified human homolog with clinical atherosclerosis in a cohort of 88 660 human subjects. Homologous human pathways and networks were significantly enriched for gene-atherosclerosis associations (empirical P<0.01 for TREM1 and LXR/RXR pathways and COX2 network). This included 12(60.0%) TREM1 pathway genes, 15(53.6%) LXR/RXR pathway genes, and 67(49.3%) COX2 network genes. Mouse analyses predicted, and human study validated, the strong association of COX2 expression (PTGS2) with increased likelihood of atherosclerosis (odds ratio, 1.68 per SD of genetically predicted gene expression; P=1.07×10-6).

CONCLUSIONS

PRESCIANT (Preclinical Science Integration and Translation) leverages published preclinical investigations to identify high-confidence pathways, networks, and regulators of human disease.

Item-Level Genome-Wide Association Study of the Alcohol Use Disorders Identification Test in Three Population-Based Cohorts

OBJECTIVE

Genome-wide association studies (GWASs) of the Alcohol Use Disorders Identification Test (AUDIT), a 10-item screen for alcohol use disorder (AUD), have elucidated novel loci for alcohol consumption and misuse. However, these studies also revealed that GWASs can be influenced by numerous biases (e.g., measurement error, selection bias), which may have led to inconsistent genetic correlations between alcohol involvement and AUD, as well as paradoxically negative genetic correlations between alcohol involvement and psychiatric disorders and/or medical conditions. The authors used genomic structural equation modeling to elucidate the genetics of alcohol consumption and problematic consequences of alcohol use as measured by AUDIT.

METHODS

To explore these unexpected differences in genetic correlations, the authors conducted the first item-level and the largest GWAS of AUDIT items (N=160,824) and applied a multivariate framework to mitigate previous biases.

RESULTS

The authors identified novel patterns of similarity (and dissimilarity) among the AUDIT items and found evidence of a correlated two-factor structure at the genetic level ("consumption" and "problems," rg=0.80). Moreover, by applying empirically derived weights to each of the AUDIT items, the authors constructed an aggregate measure of alcohol consumption that was strongly associated with alcohol dependence (rg=0.67), moderately associated with several other psychiatric disorders, and no longer positively associated with health and positive socioeconomic outcomes. Lastly, by conducting polygenic analyses in three independent cohorts that differed in their ascertainment and prevalence of AUD, the authors identified novel genetic associations between alcohol consumption, alcohol misuse, and health.

CONCLUSIONS

This work further emphasizes the value of AUDIT for both clinical and genetic studies of AUD and the importance of using multivariate methods to study genetic associations that are more closely related to AUD.

Genomic and pleiotropic analyses of resting QT interval identifies novel loci and overlap with atrial electrical disorders

The resting QT interval, an electrocardiographic (ECG) measure of ventricular myocardial repolarization, is a heritable risk marker of cardiovascular mortality, but the mechanisms remain incompletely understood. Previously reported candidate genes have provided insights into the regulatory mechanisms of the QT interval. However, there are still important knowledge gaps. We aimed to gain new insights by (i) providing new candidate genes, (ii) identifying pleiotropic associations with other cardiovascular traits, and (iii) scanning for sexually dimorphic genetic effects. We conducted a genome-wide association analysis for resting QT interval with ~9.8 million variants in 52 107 individuals of European ancestry without known cardiovascular disease from the UK Biobank. We identified 40 loci, 13 of which were novel, including 2 potential sex-specific loci, explaining ~11% of the trait variance. Candidate genes at novel loci were involved in myocardial structure and arrhythmogenic cardiomyopathy. Investigation of pleiotropic effects of QT interval variants using phenome-wide association analyses in 302 000 unrelated individuals from the UK Biobank and pairwise genome-wide comparisons with other ECG and cardiac imaging traits revealed genetic overlap with atrial electrical pathology. These findings provide novel insights into how abnormal myocardial repolarization and increased cardiovascular mortality may be linked.

Childhood Obesity and Risk of Stroke: A Mendelian Randomisation Analysis

Background: The causal relationship between childhood obesity and stroke remains unclear. Our objective was to elucidate the causal relationship between childhood obesity and the risk of stroke and its subtypes by performing Mendelian randomisation (MR) analyses. Methods: Genetic instruments for childhood obesity were obtained from a genome-wide association study (GWAS) of 13,848 European participants. Summary level data for stroke, intracerebral haemorrhage, ischaemic stroke (IS), and its subtypes were evaluated using the MEGASTROKE GWAS dataset, which included 446,696 European adults. Inverse-variance weighting, weighted-median analysis, MR-Egger regression, MR Pleiotropy RESidual Sum and Outlier test (MR-PRESSO), and MR-Robust Adjusted Profile Score were applied in this MR analysis. The leave-one-out sensitivity test, MR-PRESSO Global test, and Cochran's Q test were conducted to confirm the accuracy and robustness of our results. Results: Genetic evaluations revealed that childhood obesity was associated with a higher risk of stroke (OR = 1.04, 95%CI: 1.01-1.07, p = 0.005) and IS (OR = 1.05, 95%CI: 1.02-1.08, p = 0.003), but not with intracerebral haemorrhage (ICH, OR = 0.93, 95%CI: 0.80-1.09, p = 0.39). In the subtype analysis, childhood obesity was also associated with large artery stroke (LAS, OR = 1.12, 95%CI: 1.02-1.22, p = 0.016) but not with cardioembolic stroke (OR = 1.06, 95%CI: 0.96-1.18, p = 0.21) and small vessel stroke (OR = 1.06, 95%CI: 0.98-1.15, p = 0.17). These results were stable in the sensitivity analysis and remained significant after Bonferroni correction. Conclusion: Our study provides evidence that childhood obesity is associated with a higher risk of stroke, IS, and LAS. The prevention of stroke, especially IS and LAS, should be promoted in populations with childhood obesity.

Comparing Genetic and Socioenvironmental Contributions to Ethnic Differences in C-Reactive Protein

C-reactive protein (CRP) is a routinely measured blood biomarker for inflammation. Elevated levels of circulating CRP are associated with response to infection, risk for a number of complex common diseases, and psychosocial stress. The objective of this study was to compare the contributions of genetic ancestry, socioenvironmental factors, and inflammation-related health conditions to ethnic differences in C-reactive protein levels. We used multivariable regression to compare CRP blood serum levels between Black and White ethnic groups from the United Kingdom Biobank (UKBB) prospective cohort study. CRP serum levels are significantly associated with ethnicity in an age and sex adjusted model. Study participants who identify as Black have higher average CRP than those who identify as White, CRP increases with age, and females have higher average CRP than males. Ethnicity and sex show a significant interaction effect on CRP. Black females have higher average CRP levels than White females, whereas White males have higher average CRP than Black males. Significant associations between CRP, ethnicity, and genetic ancestry are almost completely attenuated in a fully adjusted model that includes socioenvironmental factors and inflammation-related health conditions. BMI, smoking, and socioeconomic deprivation all have high relative effects on CRP. These results indicate that socioenvironmental factors contribute more to CRP ethnic differences than genetics. Differences in CRP are associated with ethnic disparities for a number of chronic diseases, including type 2 diabetes, essential hypertension, sarcoidosis, and lupus erythematosus. Our results indicate that ethnic differences in CRP are linked to both socioenvironmental factors and numerous ethnic health disparities.

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.

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n=35) or hospitalization (n=42) due to severe COVID-19 using genome-wide association summary from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (n_rs495828=53 and n_rs505922=59); strongest association with venous embolism, odds ratio (ORrs495828 1.33 (p=1.32 × 10^-199), and thrombosis OR_rs505922 1.33, p=2.2 × 10^-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p=4.12 × 10^-191; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p=2.26 × 10^-12. The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p=6.48 × 10^-23, lupus OR 0.84, p=3.97 × 10^-06. PheWAS stratified by genetic ancestry demonstrated differences in genotype-phenotype associations across ancestry. LMNA (rs581342) associated with neutropenia OR 1.29 p=4.1 × 10^-13 among Veterans of African ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.

HIF-1α Pulmonary Phenotype Wide Association Study Unveils a Link to Inflammatory Airway Conditions

Despite experimental data linking HIF-1α dysfunction to inflammatory airway conditions, the effect of single nucleotide polymorphisms within the HIF1A gene on these conditions remains poorly understood. In the current study, we complete a phenotype wide association study to assess the link between SNPs with known disease associations and respiratory phenotypes. We report two SNPs of the HIF1A gene, the intronic rs79865957 and the missense rs41508050. In these positions the A and the T allele are significantly associated with allergic rhinitis and acute bronchitis and bronchiolitis, respectively. These findings further support the role of HIF-1α in inflammatory pulmonary conditions and may serve as a basis to refine our understanding of other HIF-1α associated phenotypes.

A genome-first approach to rare variants in hypertrophic cardiomyopathy genes MYBPC3 and MYH7 in a medical biobank

'Genome-first' approaches to analyzing rare variants can reveal new insights into human biology and disease. Because pathogenic variants are often rare, new discovery requires aggregating rare coding variants into 'gene burdens' for sufficient power. However, a major challenge is deciding which variants to include in gene burden tests. Pathogenic variants in MYBPC3 and MYH7 are well-known causes of hypertrophic cardiomyopathy (HCM), and focusing on these 'positive control' genes in a genome-first approach could help inform variant selection methods and gene burdening strategies for other genes and diseases. Integrating exome sequences with electronic health records among 41 759 participants in the Penn Medicine BioBank, we evaluated the performance of aggregating predicted loss-of-function (pLOF) and/or predicted deleterious missense (pDM) variants in MYBPC3 and MYH7 for gene burden phenome-wide association studies (PheWAS). The approach to grouping rare variants for these two genes produced very different results: pLOFs but not pDM variants in MYBPC3 were strongly associated with HCM, whereas the opposite was true for MYH7. Detailed review of clinical charts revealed that only 38.5% of patients with HCM diagnoses carrying an HCM-associated variant in MYBPC3 or MYH7 had a clinical genetic test result. Additionally, 26.7% of MYBPC3 pLOF carriers without HCM diagnoses had clear evidence of left atrial enlargement and/or septal/LV hypertrophy on echocardiography. Our study shows the importance of evaluating both pLOF and pDM variants for gene burden testing in future studies to uncover novel gene-disease relationships and identify new pathogenic loss-of-function variants across the human genome through genome-first analyses of healthcare-based populations.

Mammographic features are associated with cardiometabolic disease risk and mortality

AIMS

In recent years, microcalcifications identified in routine mammograms were found to be associated with cardiometabolic disease in women. Here, we aimed to systematically evaluate the association of microcalcifications and other mammographic features with cardiometabolic disease risk and mortality in a large screening cohort and to understand a potential genetic contribution.

METHODS AND RESULTS

This study included 57 867 women from a prospective mammographic screening cohort in Sweden (KARMA) and 49 583 sisters. Cardiometabolic disease diagnoses and mortality and medication were extracted by linkage to Swedish population registries with virtually no missing data. In the cardiometabolic phenome-wide association study, we found that a higher number of microcalcifications were associated with increased risk for multiple cardiometabolic diseases, particularly in women with pre-existing cardiometabolic diseases. In contrast, dense breasts were associated with a lower incidence of cardiometabolic diseases. Importantly, we observed similar associations in sisters of KARMA women, indicating a potential genetic overlap between mammographic features and cardiometabolic traits. Finally, we observed that the presence of microcalcifications was associated with increased cardiometabolic mortality in women with pre-existing cardiometabolic diseases (hazard ratio and 95% confidence interval: 1.79 [1.24-2.58], P = 0.002) while we did not find such effects in women without cardiometabolic diseases.

CONCLUSIONS

We found that mammographic features are associated with cardiometabolic risk and mortality. Our results strengthen the notion that a combination of mammographic features and other breast cancer risk factors could be a novel and affordable tool to assess cardiometabolic health in women attending mammographic screening.

COVID-19 susceptibility variants associate with blood clots, thrombophlebitis and circulatory diseases

Epidemiological studies suggest that individuals with comorbid conditions including diabetes, chronic lung, inflammatory and vascular disease, are at higher risk of adverse COVID-19 outcomes. Genome-wide association studies have identified several loci associated with increased susceptibility and severity for COVID-19. However, it is not clear whether these associations are genetically determined or not. We used a Phenome-Wide Association (PheWAS) approach to investigate the role of genetically determined COVID-19 susceptibility on disease related outcomes. PheWAS analyses were performed in order to identify traits and diseases related to COVID-19 susceptibility and severity, evaluated through a predictive COVID-19 risk score. We utilised phenotypic data in up to 400,000 individuals from the UK Biobank, including Hospital Episode Statistics and General Practice data. We identified a spectrum of associations between both genetically determined COVID-19 susceptibility and severity with a number of traits. COVID-19 risk was associated with increased risk for phlebitis and thrombophlebitis (OR = 1.11, p = 5.36e-08). We also identified significant signals between COVID-19 susceptibility with blood clots in the leg (OR = 1.1, p = 1.66e-16) and with increased risk for blood clots in the lung (OR = 1.12, p = 1.45 e-10). Our study identifies significant association of genetically determined COVID-19 with increased blood clot events in leg and lungs. The reported associations between both COVID-19 susceptibility and severity and other diseases adds to the identification and stratification of individuals at increased risk, adverse outcomes and long-term effects.

The genetic architecture of plasma kynurenine includes cardiometabolic disease mechanisms associated with the SH2B3 gene

Inflammation increases the risk of cardiometabolic disease. Delineating specific inflammatory pathways and biomarkers of their activity could identify the mechanistic underpinnings of the increased risk. Plasma levels of kynurenine, a metabolite involved in inflammation, associates with cardiometabolic disease risk. We used genetic approaches to identify inflammatory mechanisms associated with kynurenine variability and their relationship to cardiometabolic disease. We identified single-nucleotide polymorphisms (SNPs) previously associated with plasma kynurenine, including a missense-variant (rs3184504) in the inflammatory gene SH2B3/LNK. We examined the association between rs3184504 and plasma kynurenine in independent human samples, and measured kynurenine levels in SH2B3-knock-out mice and during human LPS-evoked endotoxemia. We conducted phenome scanning to identify clinical phenotypes associated with each kynurenine-related SNP and with a kynurenine polygenic score using the UK-Biobank (n = 456,422), BioVU (n = 62,303), and Electronic Medical Records and Genetics (n = 32,324) databases. The SH2B3 missense variant associated with plasma kynurenine levels and SH2B3-/- mice had significant tissue-specific differences in kynurenine levels.LPS, an acute inflammatory stimulus, increased plasma kynurenine in humans. Mendelian randomization showed increased waist-circumference, a marker of central obesity, associated with increased kynurenine, and increased kynurenine associated with C-reactive protein (CRP). We found 30 diagnoses associated (FDR q < 0.05) with the SH2B3 variant, but not with SNPs mapping to genes known to regulate tryptophan-kynurenine metabolism. Plasma kynurenine may be a biomarker of acute and chronic inflammation involving the SH2B3 pathways. Its regulation lies upstream of CRP, suggesting that kynurenine may be a biomarker of one inflammatory mechanism contributing to increased cardiometabolic disease risk.

Associations of Genetically Predicted Lp(a) (Lipoprotein [a]) Levels With Cardiovascular Traits in Individuals of European and African Ancestry

BACKGROUND

Lp(a) (lipoprotein [a]) levels are higher in individuals of African ancestry (AA) than in individuals of European ancestry (EA). We examined associations of genetically predicted Lp(a) levels with (1) atherosclerotic cardiovascular disease subtypes: coronary heart disease, cerebrovascular disease, peripheral artery disease, and abdominal aortic aneurysm and (2) nonatherosclerotic cardiovascular disease phenotypes, stratified by ancestry.

METHODS

We performed (1) Mendelian randomization analyses for previously reported cardiovascular associations and (2) Mendelian randomization-phenome-wide association analyses for novel associations. Analyses were stratified by ancestry in electronic Medical Records and Genomics, United Kingdom Biobank, and Million Veteran Program cohorts separately and in a combined cohort of 804 507 EA and 103 580 AA participants.

RESULTS

In Mendelian randomization analyses using the combined cohort, a 1-SD genetic increase in Lp(a) level was associated with atherosclerotic cardiovascular disease subtypes in EA-odds ratio and 95% CI for coronary heart disease 1.28 (1.16-1.41); cerebrovascular disease 1.14 (1.07-1.21); peripheral artery disease 1.22 (1.11-1.34); abdominal aortic aneurysm 1.28 (1.17-1.40); in AA, the effect estimate was lower than in EA and nonsignificant for coronary heart disease 1.11 (0.99-1.24) and cerebrovascular disease 1.06 (0.99-1.14) but similar for peripheral artery disease 1.16 (1.01-1.33) and abdominal aortic aneurysm 1.34 (1.11-1.62). In EA, a 1-SD genetic increase in Lp(a) level was associated with aortic valve disorders 1.34 (1.10-1.62), mitral valve disorders 1.18 (1.09-1.27), congestive heart failure 1.12 (1.05-1.19), and chronic kidney disease 1.07 (1.01-1.14). In AA, no significant associations were noted for aortic valve disorders 1.08 (0.94-1.25), mitral valve disorders 1.02 (0.89-1.16), congestive heart failure 1.02 (0.95-1.10), or chronic kidney disease 1.05 (0.99-1.12). Mendelian randomization-phenome-wide association analyses identified novel associations in EA with arterial thromboembolic disease, nonaortic aneurysmal disease, atrial fibrillation, cardiac conduction disorders, and hypertension.

CONCLUSIONS

Many cardiovascular associations of genetically increased Lp(a) that were significant in EA were not significant in AA. Lp(a) was associated with atherosclerotic cardiovascular disease in four major arterial beds in EA but only with peripheral artery disease and abdominal aortic aneurysm in AA. Additionally, novel cardiovascular associations were detected in EA.