Phenome-wide Association Study Research Articles

Synaptic proteome perturbations after maternal immune activation: Identification of embryonic and adult hippocampal changes

BackgroundMaternal immune activation (MIA) triggers neurobiological changes in offspring, potentially reshaping the molecular synaptic landscape, with the hippocampus being particularly vulnerable. However, critical details regarding developmental timing of these changes and whether they differ between males and females remain unclear. MethodsWe induced MIA in C57BL/6J mice on gestational day nine using the viral mimetic poly(I:C) and performed mass spectrometry-based proteomic analyses on hippocampal synaptoneurosomes of embryonic (E18) and adult (20 ± 1 weeks) MIA offspring. ResultsIn the embryonic synaptoneurosomes, MIA led to lipid, polysaccharide, and glycoprotein metabolism pathway disruptions. In the adult synaptic proteome, we observed a dynamic shift toward transmembrane trafficking, intracellular signalling cascades, including cell death and growth, and cytoskeletal organisation. In adults, many associated pathways overlapped between males and females. However, we found distinct sex-specific enrichment of dopaminergic and glutamatergic pathways. We identified 50 proteins altered by MIA in both embryonic and adult samples (28 with the same directionality), mainly involved in presynaptic structure and synaptic vesicle function. We probed human phenome-wide association study data in the cognitive and psychiatric domains, and 49 of the 50 genes encoding these proteins were significantly associated with the investigated phenotypes. ConclusionsOur data emphasise the dynamic effects of viral-like MIA on developing and mature hippocampi and provide novel targets for study following prenatal immune challenges. The 22 proteins that changed directionality from the embryonic to adult hippocampus, suggestive of compensatory over-adaptions, are particularly attractive for future investigations.

Rare damaging CCR2 variants are associated with lower lifetime cardiovascular risk.

Previous work has shown a role of CCL2, a key chemokine governing monocyte trafficking, in atherosclerosis. However, it remains unknown whether targeting CCR2, the cognate receptor of CCL2, provides protection against human atherosclerotic cardiovascular disease. Computationally predicted damaging or loss-of-function (REVEL>0.5) variants within CCR2 were detected in whole-exome-sequencing data from 454,775 UK Biobank participants and tested for association with cardiovascular endpoints in gene-burden tests. Given the key role of CCR2 in monocyte mobilization, variants associated with lower monocyte count were prioritized for experimental validation. The response to CCL2 of human cells transfected with these variants was tested in migration and cAMP assays. Validated damaging variants were tested for association with cardiovascular endpoints, atherosclerosis burden, and vascular risk factors. Significant associations were replicated in six independent datasets (n=1,062,595). Carriers of 45 predicted damaging or loss-of-function CCR2 variants (n=787 individuals) were at lower risk of myocardial infarction and coronary artery disease. One of these variants (M249K, n=585, 0.15% of European ancestry individuals) was associated with lower monocyte count and with both decreased downstream signaling and chemoattraction in response to CCL2. While M249K showed no association with conventional vascular risk factors, it was consistently associated with a lower risk of myocardial infarction (Odds Ratio [OR]: 0.66 95% Confidence Interval [CI]: 0.54-0.81, p=6.1×10-5) and coronary artery disease (OR: 0.74 95%CI: 0.63-0.87, p=2.9×10-4) in the UK Biobank and in six replication cohorts. In a phenome-wide association study, there was no evidence of a higher risk of infections among M249K carriers. Carriers of an experimentally confirmed damaging CCR2 variant are at a lower lifetime risk of myocardial infarction and coronary artery disease without carrying a higher risk of infections. Our findings provide genetic support for the translational potential of CCR2-targeting as an atheroprotective approach.

Bridging the gap: exploring the causal relationship between metformin and tumors.

Numerous studies have reported that metformin can reduce the risk of tumor development. However, some of the results of these studies are conflicting, necessitating a more reliable evaluation. We conducted a Mendelian randomization phenome-wide association study (MR-PheWAS) of tumors to explore the causal relationship between metformin and tumors. Two cohorts of patients taking metformin were obtained from the UK Biobank. Complete phenotype data of the tumors were obtained from FinnGen_R10. We elucidated the causal relationship using a two-sample Mendelian randomization (MR) analysis. More importantly, we conducted a meta-analysis to ensure relatively unbiased results. In the MR analysis, we used the inverse-variance weighted (IVW) method as the main outcome indicator. Subsequently, two cohorts were integrated for the meta-analysis. Finally, we investigated the mechanisms through mediational MR analysis. MR analysis revealed that metformin might have a causal relationship with 13 tumor-associated phenotypes in the training cohort. Four phenotypes were validated in the testing cohort. In the training and testing cohorts, metformin exhibited a protective effect against brain meningiomas and malignant neoplasms of the breast (HER-positive), oral cavity, tonsils, and the base of the tongue. Intriguingly, after integrating the results of the two cohorts for the meta-analysis, 12 results were statistically significant. Mediational MR analysis suggested that the effects of metformin on brain meningiomas may be weakened by the presence of the family Oxalobacteraceae. Metformin exhibits potential preventive and therapeutic effects on four types of tumors: brain meningioma, malignant neoplasms of the breast (HER-positive), oral cavity and tonsils, and the base of the tongue. Large randomized controlled trials are required to confirm these findings.

Genome-wide association studies of lifetime and frequency cannabis use in 131,895 individuals.

Cannabis is one of the most widely used drugs globally. Decriminalization of cannabis is further increasing cannabis consumption. We performed genome-wide association studies (GWASs) of lifetime (N=131,895) and frequency (N=73,374) of cannabis use. Lifetime cannabis use GWAS identified two loci, one near CADM2 (rs11922956, p=2.40E-11) and another near GRM3 (rs12673181, p=6.90E-09). Frequency of use GWAS identified one locus near CADM2 (rs4856591, p=8.10E-09; r2 =0.76 with rs11922956). Both traits were heritable and genetically correlated with previous GWASs of lifetime use and cannabis use disorder (CUD), as well as other substance use and cognitive traits. Polygenic scores (PGSs) for lifetime and frequency of cannabis use associated cannabis use phenotypes in AllofUs participants. Phenome-wide association study of lifetime cannabis use PGS in a hospital cohort replicated associations with substance use and mood disorders, and uncovered associations with celiac and infectious diseases. This work demonstrates the value of GWASs of CUD transition risk factors.

Genome-wide association studies of coffee intake in UK/US participants of European ancestry uncover cohort-specific genetic associations.

Coffee is one of the most widely consumed beverages. We performed a genome-wide association study (GWAS) of coffee intake in US-based 23andMe participants (N = 130,153) and identified 7 significant loci, with many replicating in three multi-ancestral cohorts. We examined genetic correlations and performed a phenome-wide association study across hundreds of biomarkers, health, and lifestyle traits, then compared our results to the largest available GWAS of coffee intake from the UK Biobank (UKB; N = 334,659). We observed consistent positive genetic correlations with substance use and obesity in both cohorts. Other genetic correlations were discrepant, including positive genetic correlations between coffee intake and psychiatric illnesses, pain, and gastrointestinal traits in 23andMe that were absent or negative in the UKB, and genetic correlations with cognition that were negative in 23andMe but positive in the UKB. Phenome-wide association study using polygenic scores of coffee intake derived from 23andMe or UKB summary statistics also revealed consistent associations with increased odds of obesity- and red blood cell-related traits, but all other associations were cohort-specific. Our study shows thatthe genetics of coffee intake associate with substance use and obesity across cohorts, but also that GWAS performed in different populations could capture cultural differences in the relationship between behavior and genetics.

Novel Protein Biomarkers and Therapeutic Targets for Type 1 Diabetes and Its Complications: Insights from Summary-Data-Based Mendelian Randomization and Colocalization Analysis.

Millions of patients suffer from type 1 diabetes (T1D) and its associated complications. Nevertheless, the pursuit of a cure for T1D has encountered significant challenges, with a crucial impediment being the lack of biomarkers that can accurately predict the progression of T1D and reliable therapeutic targets for T1D. Hence, there is an urgent need to discover novel protein biomarkers and therapeutic targets, which holds promise for targeted therapy for T1D. In this study, we extracted summary-level data on 4907 plasma proteins from 35,559 Icelanders and 2923 plasma proteins from 54,219 UK participants as exposures. The genome-wide association study (GWAS) summary statistics on T1D and T1D with complications were obtained from the R9 release results from the FinnGen consortium. Summary-data-based Mendelian randomization (SMR) analysis was employed to evaluate the causal associations between the genetically predicted levels of plasma proteins and T1D-associated outcomes. Colocalization analysis was utilized to investigate the shared genetic variants between the exposure and outcome. Moreover, transcriptome analysis and a protein-protein interaction (PPI) network further illustrated the expression patterns of the identified protein targets and their interactions with the established targets of T1D. Finally, a Mendelian randomization phenome-wide association study evaluated the potential side effects of the identified core protein targets. In the primary SMR analysis, we identified 72 potential protein targets for T1D and its complications, and nine of them were considered crucial protein targets. Within the group were five risk targets and four protective targets. Backed by evidence from the colocalization analysis, the protein targets were classified into four tiers, with MANSC4, CTRB1, SIGLEC5 and MST1 being categorized as tier 1 targets. Delving into the DrugBank database, we retrieved 11 existing medications for T1D along with their therapeutic targets. The PPI network clarified the interactions among the identified potential protein targets and established ones. Finally, the Mendelian randomization phenome-wide association study corroborated MANSC4 as a reliable target capable of mitigating the risk of various forms of diabetes, and it revealed the absence of adverse effects linked to CTRB1, SIGLEC5 and MST1. This study unveiled many protein biomarkers and therapeutic targets for T1D and its complications. Such advancements hold great promise for the progression of drug development and targeted therapy for T1D.

A Mendelian randomization study of the entire phenome to explore the causal links between epilepsy.

The causes and triggering factors of epilepsy are still unknown. The results of genome-wide association studies can be utilized for a phenome-wide association study using Mendelian randomization (MR) to identify potential risk factors for epilepsy. This study utilizes two-sample MR analysis to investigate whether 316 phenotypes, including lifestyle, environmental factors, blood biomarker, and more, are causally associated with the occurrence of epilepsy. The primary analysis employed the inverse variance weighted (IVW) model, while complementary MR analysis methods (MR Egger, Wald ratio) were also employed. Sensitivity analyses were also conducted to evaluate heterogeneity and pleiotropy. There was no evidence of a statistically significant causal association between the examined phenotypes and epilepsy following Bonferroni correction (p<1.58 × 10-4) or false discovery rate correction. The results of the MR analysis indicate that the frequency of tiredness or lethargy in the last 2weeks (p=0.042), blood uridine (p=0.003), blood propionylcarnitine (p=0.041), and free cholesterol (p=0.044) are suggestive causal risks for epilepsy. Lifestyle choices, such as sleep duration and alcohol consumption, as well as biomarkers including steroid hormone levels, hippocampal volume, and amygdala volume were not identified as causal factors for developing epilepsy (p>0.05). Our study provides additional insights into the underlying causes of epilepsy, which will serve as evidence for the prevention and control of epilepsy. The associations observed in epidemiological studies may be partially attributed to shared biological factors or lifestyle confounders.

Comorbidity Profiles of Posttraumatic Stress Disorder Across the Medical Phenome

BackgroundPrevious epidemiological research has linked posttraumatic stress disorder (PTSD) with specific physical health problems, but the comprehensive landscape of medical conditions associated with PTSD remains uncharacterized. Electronic health records provide an opportunity to overcome clinical knowledge gaps and uncover associations with biological relevance that potentially vary by sex. MethodsPTSD was defined among biobank participants (N = 145,959) in 3 major healthcare systems using 2 ICD code-based definitions: broad (≥1 PTSD or acute stress codes vs. 0; ncases = 16,706) and narrow (≥2 PTSD codes vs. 0; ncases = 3325). Using a phenome-wide association study design, we tested associations between each PTSD definition and all prevalent disease umbrella categories, i.e., phecodes. We also conducted sex-stratified phenome-wide association study analyses including a sex × diagnosis interaction term in each logistic regression. ResultsA substantial number of phecodes were significantly associated with PTSDNarrow (61%) and PTSDBroad (83%). While the strongest associations were shared between the 2 definitions, PTSDBroad captured 334 additional phecodes not significantly associated with PTSDNarrow and exhibited a wider range of significantly associated phecodes across various categories, including respiratory, genitourinary, and circulatory conditions. Sex differences were observed in that PTSDBroad was more strongly associated with osteoporosis, respiratory failure, hemorrhage, and pulmonary heart disease among male patients and with urinary tract infection, acute pharyngitis, respiratory infections, and overweight among female patients. ConclusionsThis study provides valuable insights into a diverse range of comorbidities associated with PTSD, including both known and novel associations, while highlighting the influence of sex differences and the impact of defining PTSD using electronic health records.

#1880 Disease outcomes associated with decreased estimated glomerular filtration rate: a Phenome-wide association study

Abstract Background and Aims This study aimed to identify associations between decreased creatinine-based estimated glomerular filtration rate (eGFR) and a broad range of related outcomes. Method Polygenic risk score (PRS) was calculated using PLINK1.9 software according to the tutorial with summing of risk alleles weighted by the effect sizes of 347308 white participants in UK Biobank. Genome-wide association study (GWAS) summary statistical data of creatinine-based eGFR from CKDGen consortium was utilized as base data. The phenome codes of disease outcomes were derived from the primary and secondary diagnoses of ICD-10 and ICD-9 codes in the HSE records in the U.K. Biobank. Then, the association between genetically predicted eGFR and following disease outcomes were evaluated using logistic regression adjusting for age (continuous variables), sex (male/female), assessment cancer (category variables), and the first 10 genetic principal components (PCs) (continuous variables). PRS of eGFR was scaled to a unit reflecting one SD (equals to 1) decrease of PRS. We applied the Bonferroni correction (P threshold = 0.05/665 phe codes = 7.52 × 10−5) to correct for multiple testing. Results Totally, 19857 SNPs (pruned by r2 &amp;lt; 0.1 within the 250 kb region, P &amp;lt; 0.01) from base data were retained for PRS calculation. Our Phe was identified that one standard deviation (SD) decrease in the PRS of eGFR significantly associated with chronic renal failure (OR = 1.27, 95% CI: 1.25-1.30, P = 1.66 × 10−158), intestinal malabsorption (OR = 1.18, 95% CI: 1.13-1.22, P = 1.76 × 10−16), other hypothyroidism (OR = 1.06, 95% CI: 1.04-1.07, P = 4.98 × 10−13), calculus of kidney and ureter (OR = 0.92, 95% CI: 0.89-0.94, P = 4.91 × 10−11), gout (OR = 1.07, 95% CI: 1.05-1.10, P = 2.34 × 10−8), chronic nephritic syndrome (OR = 1.21, 95% CI: 1.12-1.30, P = 3.20 × 10−7), Hypertensive renal disease (OR = 1.15, 95% CI: 1.09-1.21, P = 4.72 × 10−7), acute renal failure (OR = 1.04, 95% CI: 1.02-1.06, P = 4.58 × 10−6), glaucoma (OR = 0.95, 95% CI: 0.93-0.98, P = 2.16 × 10−5) and unspecified renal failure (OR = 1.10, 95% CI: 1.05-1.16, P = 2.70 × 10−5). Conclusion Genetically predicted deceased eGFR associated with a series of disease outcomes, providing evidence for prevention of following diseases in patients with kidney function impairment.

Proteome-wide Mendelian randomization identifies potential therapeutic targets for nonalcoholic fatty liver diseases

Nonalcoholic fatty liver disease (NAFLD) is the predominant cause of liver pathology. Current evidence highlights plasma proteins as potential therapeutic targets. However, their mechanistic roles in NAFLD remain unclear. This study investigated the involvement of specific plasma proteins and intermediate risk factors in NAFLD progression. Two-sample Mendelian randomization (MR) analysis was conducted to examine the association between plasma proteins and NAFLD. Colocalization analysis determined the shared causal variants between the identified proteins and NAFLD. The MR analysis was applied separately to proteins, risk factors, and NAFLD. Mediator shares were computed by detecting the correlations among these elements. Phenome-wide association studies (phewas) were utilized to assess the safety implications of targeting these proteins. Among 1,834 cis-protein quantitative trait loci (cis-pQTLs), after-FDR correction revealed correlations between the plasma levels of four gene-predicted proteins (CSPG3, CILP2, Apo-E, and GCKR) and NAFLD. Colocalization analysis indicated shared causal variants for CSPG3 and GCKR in NAFLD (posterior probability > 0.8). Out of the 22 risk factors screened for MR analysis, only 8 showed associations with NAFLD (p ≤ 0.05), while 4 linked to CSPG3 and GCKR. The mediator shares for these associations were calculated separately. Additionally, reverse MR analysis was performed on the pQTLs, risk factors, and NAFLD, which exhibited a causal relationship with forward MR analysis. Finally, phewas summarized the potential side effects of associated-targeting proteins, including CSPG3 and GCKR. Our research emphasized the potential therapeutic targets for NAFLD and provided modifiable risk factors for preventing NAFLD.

#985 Disease outcomes of chronic kidney disease: a Phenome-wide association study

Abstract Background and Aims Chronic kidney disease (CKD) is one of the prevalent diseases worldwide leading to increasing global disease burden, however, comprehensive studies of the subsequently developed disease outcomes are lacking. We aimed to explore the associations between CKD and a broad range of following disease conditions. Method Polygenic risk score (PRS) of 347308 white participants in UK Biobank were calculated using PLINK1.9 software with summing of risk alleles weighted by the effect sizes of. Genome-wide association study (GWAS) summary statistical data of CKD from CKDGen consortium was utilized as base data. The phenome codes of disease outcomes were derived from the primary and secondary diagnoses of ICD-10 and ICD-9 codes in the HSE records in the U.K. Biobank. Then, the association between genetically predicted eGFR and following disease outcomes were evaluated using logistic regression adjusting for age (continuous variables), sex (male/female), assessment cancer (category variables), and the first 10 genetic principal components (PCs) (continuous variables). PRS of eGFR was scaled to a unit reflecting one SD (equals to 1) decrease of PRS. We applied the Bonferroni correction (P threshold = 0.05/665 phecodes = 7.52 × 10−5) to correct for multiple testing. Results Totally, 261,254 SNPs (pruned by r2 &amp;lt; 0.1 within the 250 kb region, P &amp;lt; 0.4) from base data were retained for PRS calculation. Our PheWas identified that one standard deviation (SD) increase in the PRS of CKD significantly associated with chronic renal failure (OR = 1.16, 95% CI: 1.14-1.18, P = 2.01 × 10−66), essential (primary) hypertension (OR = 1.05, 95% CI: 1.04-1.06, P = 2.12 × 10−29), hypertensive renal disease (OR = 1.22, 95% CI: 1.16-1.29, P = 3.81 × 10−14), gout (OR = 1.09, 95% CI: 1.07-1.12, P = 6.57 × 10−13), intestinal malabsorption (OR = 1.12, 95% CI: 1.08-1.16, P = 6.71 × 10−9), acute renal failure (OR = 1.05, 95% CI: 1.03-1.07, P = 2.08 × 10−8), disorders of lipoprotein metabolism and other lipidaemias (OR = 1.03, 95% CI: 1.01-1.04, P = 9.11 × 10−7), non-insulin-dependent diabetes mellitus (OR = 1.03, 95% CI: 1.02-1.04, P = 3.11 × 10−6), obesity (OR = 1.03, 95% CI: 1.02-1.05, P = 3.13 × 10−6), unspecified renal failure (OR = 1.11, 95% CI: 1.06-1.16, P = 1.18 × 10−5), chronic ischaemic heart disease (OR = 1.03, 95% CI: 1.01-1.04, P = 1.52 × 10−5). Conclusion Increased risk of CKD significantly associated with a series of subsequent conditions, providing evidence for prevention of following diseases in patients with CKD.

The landscape of rare genetic variation associated with inflammatory bowel disease and Parkinson’s disease comorbidity

BackgroundInflammatory bowel disease (IBD) and Parkinson’s disease (PD) are chronic disorders that have been suggested to share common pathophysiological processes. LRRK2 has been implicated as playing a role in both diseases. Exploring the genetic basis of the IBD-PD comorbidity through studying high-impact rare genetic variants can facilitate the identification of the novel shared genetic factors underlying this comorbidity.MethodsWe analyzed whole exomes from the BioMe BioBank and UK Biobank, and whole genomes from a cohort of 67 European patients diagnosed with both IBD and PD to examine the effects of LRRK2 missense variants on IBD, PD and their co-occurrence (IBD-PD). We performed optimized sequence kernel association test (SKAT-O) and network-based heterogeneity clustering (NHC) analyses using high-impact rare variants in the IBD-PD cohort to identify novel candidate genes, which we further prioritized by biological relatedness approaches. We conducted phenome-wide association studies (PheWAS) employing BioMe BioBank and UK Biobank whole exomes to estimate the genetic relevance of the 14 prioritized genes to IBD-PD.ResultsThe analysis of LRRK2 missense variants revealed significant associations of the G2019S and N2081D variants with IBD-PD in addition to several other variants as potential contributors to increased or decreased IBD-PD risk. SKAT-O identified two significant genes, LRRK2 and IL10RA, and NHC identified 6 significant gene clusters that are biologically relevant to IBD-PD. We observed prominent overlaps between the enriched pathways in the known IBD, PD, and candidate IBD-PD gene sets. Additionally, we detected significantly enriched pathways unique to the IBD-PD, including MAPK signaling, LPS/IL-1 mediated inhibition of RXR function, and NAD signaling. Fourteen final candidate IBD-PD genes were prioritized by biological relatedness methods. The biological importance scores estimated by protein–protein interaction networks and pathway and ontology enrichment analyses indicated the involvement of genes related to immunity, inflammation, and autophagy in IBD-PD. Additionally, PheWAS provided support for the associations of candidate genes with IBD and PD.ConclusionsOur study confirms and uncovers new LRRK2 associations in IBD-PD. The identification of novel inflammation and autophagy-related genes supports and expands previous findings related to IBD-PD pathogenesis, and underscores the significance of therapeutic interventions for reducing systemic inflammation.

To weight or not to weight? The effect of selection bias in 3 large electronic health record-linked biobanks and recommendations for practice.

To develop recommendations regarding the use of weights to reduce selection bias for commonly performed analyses using electronic health record (EHR)-linked biobank data. We mapped diagnosis (ICD code) data to standardized phecodes from 3 EHR-linked biobanks with varying recruitment strategies: All of Us (AOU; n = 244071), Michigan Genomics Initiative (MGI; n = 81243), and UK Biobank (UKB; n = 401167). Using 2019 National Health Interview Survey data, we constructed selection weights for AOU and MGI to represent the US adult population more. We used weights previously developed for UKB to represent the UKB-eligible population. We conducted 4 common analyses comparing unweighted and weighted results. For AOU and MGI, estimated phecode prevalences decreased after weighting (weighted-unweighted median phecode prevalence ratio [MPR]: 0.82 and 0.61), while UKB estimates increased (MPR: 1.06). Weighting minimally impacted latent phenome dimensionality estimation. Comparing weighted versus unweighted phenome-wide association study for colorectal cancer, the strongest associations remained unaltered, with considerable overlap in significant hits. Weighting affected the estimated log-odds ratio for sex and colorectal cancer to align more closely with national registry-based estimates. Weighting had a limited impact on dimensionality estimation and large-scale hypothesis testing but impacted prevalence and association estimation. When interested in estimating effect size, specific signals from untargeted association analyses should be followed up by weighted analysis. EHR-linked biobanks should report recruitment and selection mechanisms and provide selection weights with defined target populations. Researchers should consider their intended estimands, specify source and target populations, and weight EHR-linked biobank analyses accordingly.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses.

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

Phenome-wide identification of therapeutic genetic targets, leveraging knowledge graphs, graph neural networks, and UK Biobank data.

The ongoing expansion of human genomic datasets propels therapeutic target identification; however, extracting gene-disease associations from gene annotations remains challenging. Here, we introduce Mantis-ML 2.0, a framework integrating AstraZeneca's Biological Insights Knowledge Graph and numerous tabular datasets, to assess gene-disease probabilities throughout the phenome. We use graph neural networks, capturing the graph's holistic structure, and train them on hundreds of balanced datasets via a robust semi-supervised learning framework to provide gene-disease probabilities across the human exome. Mantis-ML 2.0 incorporates natural language processing to automate disease-relevant feature selection for thousands of diseases. The enhanced models demonstrate a 6.9% average classification power boost, achieving a median receiver operating characteristic (ROC) area under curve (AUC) score of 0.90 across 5220 diseases from Human Phenotype Ontology, OpenTargets, and Genomics England. Notably, Mantis-ML 2.0 prioritizes associations from an independent UK Biobank phenome-wide association study (PheWAS), providing a stronger form of triaging and mitigating against underpowered PheWAS associations. Results are exposed through an interactive web resource.

Distinguishing vulnerability and resilience to posttraumatic stress disorder evaluating traumatic experiences, genetic risk and electronic health records

What distinguishes vulnerability and resilience to posttraumatic stress disorder (PTSD) remains unclear. Levering traumatic experiences reporting, genetic data, and electronic health records (EHR), we investigated and predicted the clinical comorbidities (co-phenome) of PTSD vulnerability and resilience in the UK Biobank (UKB) and All of Us Research Program (AoU), respectively. In 60,354 trauma-exposed UKB participants, we defined PTSD vulnerability and resilience considering PTSD symptoms, trauma burden, and polygenic risk scores. EHR-based phenome-wide association studies (PheWAS) were conducted to dissect the co-phenomes of PTSD vulnerability and resilience. Significant diagnostic endpoints were applied as weights, yielding a phenotypic risk score (PheRS) to conduct PheWAS of PTSD vulnerability and resilience PheRS in up to 95,761 AoU participants. EHR-based PheWAS revealed three significant phenotypes positively associated with PTSD vulnerability (top association "Sleep disorders") and five outcomes inversely associated with PTSD resilience (top association "Irritable Bowel Syndrome"). In the AoU cohort, PheRS analysis showed a partial inverse relationship between vulnerability and resilience with distinct comorbid associations. While PheRSvulnerability associations were linked to multiple phenotypes, PheRSresilience showed inverse relationships with eye conditions. Our study unveils phenotypic differences in PTSD vulnerability and resilience, highlighting that these concepts are not simply the absence and presence of PTSD.

Phenomic landscape and pharmacogenomic implications for HLA region in a Taiwan Han Chinese population

BackgroundThe human leukocyte antigen (HLA) genes, exhibiting significant genetic diversity, are associated with susceptibility to various clinical diseases and diverse in drug responses. High costs of HLA sequencing and the population-specific architecture of this genetic region necessitate the establishment of a population-specific HLA imputation reference panel. Moreover, there is a lack of understanding about the genetic and phenotypic landscape of HLA variations within the Taiwanese population.MethodsWe created models for a Taiwanese-specific HLA imputation reference panel. These models were trained with the array genotype data and HLA sequencing data from 845 Taiwanese subjects. HLA imputation was applied for 59,448 Taiwanese subjects to characterize the HLA allele and haplotype frequencies. Additionally, a phenome-wide association study (PheWAS) was conducted to identify the phenotypes associated with HLA variations. The association of the biallelic HLA variants with the binary and quantitative traits were evaluated with additive logistic and linear regression models, respectively. Furthermore, an omnibus test with likelihood-ratio test was applied for each HLA amino acid position in the multiallelic HLA amino acid polymorphisms to compare the difference between a fitted model and a null model following a χ2 distribution of n-1 degree of freedom at a position with n residues. Finally, we estimated the prevalence of adverse drug reactions (ADR)-related HLA alleles in the Taiwanese population.ResultsIn this study, the reference panel models displayed remarkable accuracy, with averages of 99.3%, 98.9%, and 99.1% for 2-, 4-, 6-digit alleles of the eight classical HLA genes, respectively. For PheWAS, a total of 18,136 significant associations with HLA variants across 26 phenotypes are identified (p < 5×10-8), highlighting the pleiotropy feature of the HLA region. Among the independent signals, 15 are novel, including the association of HLA-B pos 138 variation with ankylosing spondylitis (AS), and rs9266290 and rs9266292 with allergy. Through an analysis spanning the entire HLA region, we identified clusters of phenotype correlations. Finally, the carriers of pharmacogenomic related HLA alleles, including HLA-C*01:02 (35.86%), HLA-B*58:01 (20.9%), and HLA-B*15:02 (8.38%), were characterized in the Taiwanese general population.ConclusionsWe successfully delivered the HLA imputation for 59,448 Taiwanese subjects and characterized the genetic and phenotypic landscapes of the HLA variations. In addition, we quantified the estimated prevalence of the ADR-related HLA alleles in the Taiwanese population. The developed HLA imputation reference panel could be used for estimation of population HLA allele frequencies, which can facilitate further studies in the role of HLA variants in a wider range of phenotypes in the population.

A clustering approach to improve our understanding of the genetic and phenotypic complexity of chronic kidney disease

Chronic kidney disease (CKD) is a complex disorder that causes a gradual loss of kidney function, affecting approximately 9.1% of the world's population. Here, we use a soft-clustering algorithm to deconstruct its genetic heterogeneity. First, we selected 322 CKD-associated independent genetic variants from published genome-wide association studies (GWAS) and added association results for 229 traits from the GWAS catalog. We then applied nonnegative matrix factorization (NMF) to discover overlapping clusters of related traits and variants. We computed cluster-specific polygenic scores and validated each cluster with a phenome-wide association study (PheWAS) on the BioMe biobank (n = 31,701). NMF identified nine clusters that reflect different aspects of CKD, with the top-weighted traits signifying areas such as kidney function, type 2 diabetes (T2D), and body weight. For most clusters, the top-weighted traits were confirmed in the PheWAS analysis. Results were found to be more significant in the cross-ancestry analysis, although significant ancestry-specific associations were also identified. While all alleles were associated with a decreased kidney function, associations with CKD-related diseases (e.g., T2D) were found only for a smaller subset of variants and differed across genetic ancestry groups. Our findings leverage genetics to gain insights into the underlying biology of CKD and investigate population-specific associations.

Genome-wide association studies reveal differences in genetic susceptibility between single events vs. recurrent events of atrial fibrillation and myocardial infarction: the HUNT study.

Genetic research into atrial fibrillation (AF) and myocardial infarction (MI) has predominantly focused on comparing afflicted individuals with their healthy counterparts. However, this approach lacks granularity, thus overlooking subtleties within patient populations. In this study, we explore the distinction between AF and MI patients who experience only a single disease event and those experiencing recurrent events. Integrating hospital records, questionnaire data, clinical measurements, and genetic data from more than 500,000 HUNT and United Kingdom Biobank participants, we compare both clinical and genetic characteristics between the two groups using genome-wide association studies (GWAS) meta-analyses, phenome-wide association studies (PheWAS) analyses, and gene co-expression networks. We found that the two groups of patients differ in both clinical characteristics and genetic risks. More specifically, recurrent AF patients are significantly younger and have better baseline health, in terms of reduced cholesterol and blood pressure, than single AF patients. Also, the results of the GWAS meta-analysis indicate that recurrent AF patients seem to be at greater genetic risk for recurrent events. The PheWAS and gene co-expression network analyses highlight differences in the functions associated with the sets of single nucleotide polymorphisms (SNPs) and genes for the two groups. However, for MI patients, we found that those experiencing single events are significantly younger and have better baseline health than those with recurrent MI, yet they exhibit higher genetic risk. The GWAS meta-analysis mostly identifies genetic regions uniquely associated with single MI, and the PheWAS analysis and gene co-expression networks support the genetic differences between the single MI and recurrent MI groups. In conclusion, this work has identified novel genetic regions uniquely associated with single MI and related PheWAS analyses, as well as gene co-expression networks that support the genetic differences between the patient subgroups of single and recurrent occurrence for both MI and AF.

Large-scale cross-ancestry genome-wide meta-analysis of serum urate

Hyperuricemia is an essential causal risk factor for gout and is associated with cardiometabolic diseases. Given the limited contribution of East Asian ancestry to genome-wide association studies of serum urate, the genetic architecture of serum urate requires exploration. A large-scale cross-ancestry genome-wide association meta-analysis of 1,029,323 individuals and ancestry-specific meta-analysis identifies a total of 351 loci, including 17 previously unreported loci. The genetic architecture of serum urate control is similar between European and East Asian populations. A transcriptome-wide association study, enrichment analysis, and colocalization analysis in relevant tissues identify candidate serum urate-associated genes, including CTBP1, SKIV2L, and WWP2. A phenome-wide association study using polygenic risk scores identifies serum urate-correlated diseases including heart failure and hypertension. Mendelian randomization and mediation analyses show that serum urate-associated genes might have a causal relationship with serum urate-correlated diseases via mediation effects. This study elucidates our understanding of the genetic architecture of serum urate control.