Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Novel insight into the etiology of ischemic stroke gained by integrative multiome-wide association study.

Stroke, characterized by sudden neurological deficits, is the second leading cause of death worldwide. Although genome-wide association studies (GWAS) have successfully identified many genomic regions associated with ischemic stroke (IS), the genes underlying risk and their regulatory mechanisms remain elusive. Here, we integrate a large-scale GWAS (N = 1 296 908) for IS together with molecular QTLs data, including mRNA, splicing, enhancer RNA (eRNA), and protein expression data from up to 50 tissues (total N = 11 588). We identify 136 genes/eRNA/proteins associated with IS risk across 60 independent genomic regions and find IS risk is most enriched for eQTLs in arterial and brain-related tissues. Focusing on IS-relevant tissues, we prioritize 9 genes/proteins using probabilistic fine-mapping TWAS analyses. In addition, we discover that blood cell traits, particularly reticulocyte cells, have shared genetic contributions with IS using TWAS-based pheWAS and genetic correlation analysis. Lastly, we integrate our findings with a large-scale pharmacological database and identify a secondary bile acid, deoxycholic acid, as a potential therapeutic component. Our work highlights IS risk genes/splicing-sites/enhancer activity/proteins with their phenotypic consequences using relevant tissues as well as identify potential therapeutic candidates for IS.

Are phenome-wide association studies feasible in a developing country?

Phenome-wide association studies (PheWASs), a powerful approach that examines phenotypes associated with a genetic marker, have been used extensively in highly developed countries. Although there may be a clear need for PheWAS in a developing country such as the Philippines, limitations related to resources and practicality would make conducting them a challenge.

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.

Combining evidence from Mendelian randomization and colocalization: Review and comparison of approaches.

Mendelian randomization and colocalization are two statistical approaches that can be applied to summarized data from genome-wide association studies (GWASs) to understand relationships between traits and diseases. However, despite similarities in scope, they are different in their objectives, implementation, and interpretation, in part because they were developed to serve different scientific communities. Mendelian randomization assesses whether genetic predictors of an exposure are associated with the outcome and interprets an association as evidence that the exposure has a causal effect on the outcome, whereas colocalization assesses whether two traits are affected by the same or distinct causal variants. When considering genetic variants in a single genetic region, both approaches can be performed. While a positive colocalization finding typically implies a non-zero Mendelian randomization estimate, the reverse is not generally true: there are several scenarios which would lead to a non-zero Mendelian randomization estimate but lack evidence for colocalization. These include the existence of distinct but correlated causal variants for the exposure and outcome, which would violate the Mendelian randomization assumptions, and a lack of strong associations with the outcome. As colocalization was developed in the GWAS tradition, typically evidence for colocalization is concluded only when there is strong evidence for associations with both traits. In contrast, a non-zero estimate from Mendelian randomization can be obtained despite only nominally significant genetic associations with the outcome at the locus. In this review, we discuss how the two approaches can provide complementary information on potential therapeutic targets.

Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study.

BACKGROUND: Male-pattern baldness (MPB) is the most common cause of hair loss in men. It can be categorized into three types: type 2 (T2), type 3 (T3), and type 4 (T4), with type 1 (T1) being considered normal. Although various MPB-associated genetic variants have been suggested, a comprehensive study for linking these variants to gene expression regulation has not been performed to the best of our knowledge. RESULTS: In this study, we prioritized MPB-related tissue panels using tissue-specific enrichment analysis and utilized single-tissue panels from genotype-tissue expression version 8, as well as cross-tissue panels from context-specific genetics. Through a transcriptome-wide association study and colocalization analysis, we identified 52, 75, and 144 MPB associations for T2, T3, and T4, respectively. To assess the causality of MPB genes, we performed a conditional and joint analysis, which revealed 10, 11, and 54 putative causality genes for T2, T3, and T4, respectively. Finally, we conducted drug repositioning and identified potential drug candidates that are connected to MPB-associated genes. CONCLUSIONS: Overall, through an integrative analysis of gene expression and genotype data, we have identified robust MPB susceptibility genes that may help uncover the underlying molecular mechanisms and the novel drug candidates that may alleviate MPB.

Stress combinations and their interactions in plants database: a one-stop resource on combined stress responses in plants.

We have developed a compendium and interactive platform, named Stress Combinations and their Interactions in Plants Database (SCIPDb; http://www.nipgr.ac.in/scipdb.php), which offers information on morpho-physio-biochemical (phenome) and molecular (transcriptome and metabolome) responses of plants to different stress combinations. SCIPDb is a plant stress informatics hub for data mining on phenome, transcriptome, trait-gene ontology, and data-driven research for advancing mechanistic understanding of combined stress biology. We analyzed global phenome data from 939 studies to delineate the effects of various stress combinations on yield in major crops and found that yield was substantially affected under abiotic-abiotic stresses. Transcriptome datasets from 36 studies hosted in SCIPDb identified novel genes, whose roles have not been earlier established in combined stress. Integretome analysis under combined drought-heat stress pinpointed carbohydrate, amino acid, and energy metabolism pathways as the crucial metabolic, proteomic, and transcriptional components in plant tolerance to combined stress. These examples illustrate the application of SCIPDb in identifying novel genes and pathways involved in combined stress tolerance. Further, we showed the application of this database in identifying novel candidate genes and pathways for combined drought and pathogen stress tolerance. To our knowledge, SCIPDb is the only publicly available platform offering combined stress-specific omics big data visualization tools, such as an interactive scrollbar, stress matrix, radial tree, global distribution map, meta-phenome analysis, search, BLAST, transcript expression pattern table, Manhattan plot, and co-expression network. These tools facilitate a better understanding of the mechanisms underlying plant responses to combined stresses.

Genome-scale metabolic network model and phenome of solvent-tolerant Pseudomonas putida S12.

BACKGROUND: Pseudomonas putida S12 is a gram-negative bacterium renowned for its high tolerance to organic solvents and metabolic versatility, making it attractive for various applications, including bioremediation and the production of aromatic compounds, bioplastics, biofuels, and value-added compounds. However, a metabolic model of S12 has yet to be developed. RESULTS: In this study, we present a comprehensive and highly curated genome-scale metabolic network model of S12 (iSH1474), containing 1,474 genes, 1,436 unique metabolites, and 2,938 metabolic reactions. The model was constructed by leveraging existing metabolic models and conducting comparative analyses of genomes and phenomes. Approximately 2,000 different phenotypes were measured for S12 and its closely related KT2440 strain under various nutritional and environmental conditions. These phenotypic data, combined with the reported experimental data, were used to refine and validate the reconstruction. Model predictions quantitatively agreed well with in vivo flux measurements and the batch cultivation of S12, which demonstrated that iSH1474 accurately represents the metabolic capabilities of S12. Furthermore, the model was simulated to investigate the maximum theoretical metabolic capacity of S12 growing on toxic organic solvents. CONCLUSIONS: iSH1474 represents a significant advancement in our understanding of the cellular metabolism of P. putida S12. The combined results of metabolic simulation and comparative genome and phenome analyses identified the genetic and metabolic determinants of the characteristic phenotypes of S12. This study could accelerate the development of this versatile organism as an efficient cell factory for various biotechnological applications.

The Role of Electronic Health Records in Advancing Genomic Medicine.

Recent advances in genomic technology and widespread adoption of electronic health records (EHRs) have accelerated the development of genomic medicine, bringing promising research findings from genome science into clinical practice. Genomic and phenomic data, accrued across large populations through biobanks linked to EHRs, have enabled the study of genetic variation at a phenome-wide scale. Through new quantitative techniques, pleiotropy can be explored with phenome-wide association studies, the occurrence of common complex diseases can be predicted using the cumulative influence of many genetic variants (polygenic risk scores), and undiagnosed Mendelian syndromes can be identified using EHR-based phenotypic signatures (phenotype risk scores). In this review, we trace the role of EHRs from the development of genome-wide analytic techniques to translational efforts to test these new interventions to the clinic. Throughout, we describe the challenges that remain when combining EHRs with genetics to improve clinical care.

A powerful subset-based method identifies gene set associations and improves interpretation in UK Biobank.

Tests of association between a phenotype and a set of genes in a biological pathway can provide insights into the genetic architecture of complex phenotypes beyond those obtained from single-variant or single-gene association analysis. However, most existing gene set tests have limited power to detect gene set-phenotype association when a small fraction of the genes are associated with the phenotype and cannot identify the potentially "active" genes that might drive a gene set-based association. To address these issues, we have developed Gene set analysis Association Using Sparse Signals (GAUSS), a method for gene set association analysis that requires only GWAS summary statistics. For each significantly associated gene set, GAUSS identifies the subset of genes that have the maximal evidence of association and can best account for the gene set association. Using pre-computed correlation structure among test statistics from a reference panel, our p value calculation is substantially faster than other permutation- or simulation-based approaches. In simulations with varying proportions of causal genes, we find that GAUSS effectively controls type 1 error rate and has greater power than several existing methods, particularly when a small proportion of genes account for the gene set signal. Using GAUSS, we analyzed UK Biobank GWAS summary statistics for 10,679 gene sets and 1,403 binary phenotypes. We found that GAUSS is scalable and identified 13,466 phenotype and gene set association pairs. Within these gene sets, we identify an average of 17.2 (max = 405) genes that underlie these gene set associations.

Leveraging health systems data to characterize a large effect variant conferring risk for liver disease in Puerto Ricans.

The integration of genomic data into health systems offers opportunities to identify genomic factors underlying the continuum of rare and common disease. We applied a population-scale haplotype association approach based on identity-by-descent (IBD) in a large multi-ethnic biobank to a spectrum of disease outcomes derived from electronic health records (EHRs) and uncovered a risk locus for liver disease. We used genome sequencing and in silico approaches to fine-map the signal to a non-coding variant (c.2784-12T>C) in the gene ABCB4. In vitro analysis confirmed the variant disrupted splicing of the ABCB4 pre-mRNA. Four of five homozygotes had evidence of advanced liver disease, and there was a significant association with liver disease among heterozygotes, suggesting the variant is linked to increased risk of liver disease in an allele dose-dependent manner. Population-level screening revealed the variant to be at a carrier rate of 1.95% in Puerto Rican individuals, likely as the result of a Puerto Rican founder effect. This work demonstrates that integrating EHR and genomic data at a population scale can facilitate strategies for understanding the continuum of genomic risk for common diseases, particularly in populations underrepresented in genomic medicine.

Efficient mixed model approach for large-scale genome-wide association studies of ordinal categorical phenotypes.

In genome-wide association studies, ordinal categorical phenotypes are widely used to measure human behaviors, satisfaction, and preferences. However, because of the lack of analysis tools, methods designed for binary or quantitative traits are commonly used inappropriately to analyze categorical phenotypes. To accurately model the dependence of an ordinal categorical phenotype on covariates, we propose an efficient mixed model association test, proportional odds logistic mixed model (POLMM). POLMM is computationally efficient to analyze large datasets with hundreds of thousands of samples, can control type I error rates at a stringent significance level regardless of the phenotypic distribution, and is more powerful than alternative methods. In contrast, the standard linear mixed model approaches cannot control type I error rates for rare variants when the phenotypic distribution is unbalanced, although they performed well when testing common variants. We applied POLMM to 258 ordinal categorical phenotypes on array genotypes and imputed samples from 408,961 individuals in UK Biobank. In total, we identified 5,885 genome-wide significant variants, of which, 424 variants (7.2%) are rare variants with MAF < 0.01.

Genetic investigation into the broad health implications of caffeine: evidence from phenome-wide, proteome-wide and metabolome-wide Mendelian randomization.

BACKGROUND: Caffeine is one of the most utilized drugs in the world, yet its clinical effects are not fully understood. Circulating caffeine levels are influenced by the interplay between consumption behaviour and metabolism. This study aimed to investigate the effects of circulating caffeine levels by considering genetically predicted variation in caffeine metabolism. METHODS: Leveraging genetic variants related to caffeine metabolism that affect its circulating levels, we investigated the clinical effects of plasma caffeine in a phenome-wide association study (PheWAS). We validated novel findings using a two-sample Mendelian randomization framework and explored the potential mechanisms underlying these effects in proteome-wide and metabolome-wide Mendelian randomization. RESULTS: Higher levels of genetically predicted circulating caffeine among caffeine consumers were associated with a lower risk of obesity (odds ratio (OR) per standard deviation increase in caffeine = 0.97, 95% confidence interval (CI) CI: 0.95-0.98, p = 2.47 × 10-4), osteoarthrosis (OR = 0.97, 95% CI: 0.96-0.98, P=1.10 × 10-8) and osteoarthritis (OR: 0.97, 95% CI: 0.96 to 0.98, P = 1.09 × 10-6). Approximately one third of the protective effect of plasma caffeine on osteoarthritis risk was estimated to be mediated through lower bodyweight. Proteomic and metabolomic perturbations indicated lower chronic inflammation, improved lipid profiles, and altered protein and glycogen metabolism as potential biological mechanisms underlying these effects. CONCLUSIONS: We report novel evidence suggesting that long-term increases in circulating caffeine may reduce bodyweight and the risk of osteoarthrosis and osteoarthritis. We confirm prior genetic evidence of a protective effect of plasma caffeine on risk of overweight and obesity. Further clinical study is warranted to understand the translational relevance of these findings before clinical practice or lifestyle interventions related to caffeine consumption are introduced.

Polygenic risk score-based phenome-wide association study of head and neck cancer across two large biobanks.

BACKGROUND: Numerous observational studies have highlighted associations of genetic predisposition of head and neck squamous cell carcinoma (HNSCC) with diverse risk factors, but these findings are constrained by design limitations of observational studies. In this study, we utilized a phenome-wide association study (PheWAS) approach, incorporating a polygenic risk score (PRS) derived from a wide array of genomic variants, to systematically investigate phenotypes associated with genetic predisposition to HNSCC. Furthermore, we validated our findings across heterogeneous cohorts, enhancing the robustness and generalizability of our results. METHODS: We derived PRSs for HNSCC and its subgroups, oropharyngeal cancer and oral cancer, using large-scale genome-wide association study summary statistics from the Genetic Associations and Mechanisms in Oncology Network. We conducted a comprehensive investigation, leveraging genotyping data and electronic health records from 308,492 individuals in the UK Biobank and 38,401 individuals in the Penn Medicine Biobank (PMBB), and subsequently performed PheWAS to elucidate the associations between PRS and a wide spectrum of phenotypes. RESULTS: We revealed the HNSCC PRS showed significant association with phenotypes related to tobacco use disorder (OR, 1.06; 95% CI, 1.05-1.08; P = 3.50 × 10-15), alcoholism (OR, 1.06; 95% CI, 1.04-1.09; P = 6.14 × 10-9), alcohol-related disorders (OR, 1.08; 95% CI, 1.05-1.11; P = 1.09 × 10-8), emphysema (OR, 1.11; 95% CI, 1.06-1.16; P = 5.48 × 10-6), chronic airway obstruction (OR, 1.05; 95% CI, 1.03-1.07; P = 2.64 × 10-5), and cancer of bronchus (OR, 1.08; 95% CI, 1.04-1.13; P = 4.68 × 10-5). These findings were replicated in the PMBB cohort, and sensitivity analyses, including the exclusion of HNSCC cases and the major histocompatibility complex locus, confirmed the robustness of these associations. Additionally, we identified significant associations between HNSCC PRS and lifestyle factors related to smoking and alcohol consumption. CONCLUSIONS: The study demonstrated the potential of PRS-based PheWAS in revealing associations between genetic risk factors for HNSCC and various phenotypic traits. The findings emphasized the importance of considering genetic susceptibility in understanding HNSCC and highlighted shared genetic bases between HNSCC and other health conditions and lifestyles.

A quality control portal for sequencing data deposited at the European genome-phenome archive.

Since its launch in 2008, the European Genome-Phenome Archive (EGA) has been leading the archiving and distribution of human identifiable genomic data. In this regard, one of the community concerns is the potential usability of the stored data, as of now, data submitters are not mandated to perform any quality control (QC) before uploading their data and associated metadata information. Here, we present a new File QC Portal developed at EGA, along with QC reports performed and created for 1 694 442 files [Fastq, sequence alignment map (SAM)/binary alignment map (BAM)/CRAM and variant call format (VCF)] submitted at EGA. QC reports allow anonymous EGA users to view summary-level information regarding the files within a specific dataset, such as quality of reads, alignment quality, number and type of variants and other features. Researchers benefit from being able to assess the quality of data prior to the data access decision and thereby, increasing the reusability of data (https://ega-archive.org/blog/data-upcycling-powered-by-ega/).

An early-onset specific polygenic risk score optimizes age-based risk estimate and stratification of prostate cancer: population-based cohort study.

BACKGROUND: Early-onset prostate cancer (EOPC, ≤ 55 years) has a unique clinical entity harboring high genetic risk, but the majority of EOPC patients still substantial opportunity to be early-detected thus suffering an unfavorable prognosis. A refined understanding of age-based polygenic risk score (PRS) for prostate cancer (PCa) would be essential for personalized risk stratification. METHODS: We included 167,517 male participants [4882 cases including 205 EOPC and 4677 late-onset PCa (LOPC)] from UK Biobank. A General-, an EOPC- and an LOPC-PRS were derived from age-specific genome-wide association studies. Weighted Cox proportional hazard models were applied to estimate the risk of PCa associated with PRSs. The discriminatory capability of PRSs were validated using time-dependent receiver operating characteristic (ROC) curves with additional 4238 males from PLCO and TCGA. Phenome-wide association studies underlying Mendelian Randomization were conducted to discover EOPC linking phenotypes. RESULTS: The 269-PRS calculated via well-established risk variants was more strongly associated with risk of EOPC [hazard ratio (HR) = 2.35, 95% confidence interval (CI) 1.99-2.78] than LOPC (HR = 1.95, 95% CI 1.89-2.01; I2 = 79%). EOPC-PRS was dramatically related to EOPC risk (HR = 4.70, 95% CI 3.98-5.54) but not to LOPC (HR = 0.98, 95% CI 0.96-1.01), while LOPC-PRS had similar risk estimates for EOPC and LOPC (I2 = 0%). Particularly, EOPC-PRS performed optimal discriminatory capability for EOPC (area under the ROC = 0.613). Among the phenomic factors to PCa deposited in the platform of ProAP (Prostate cancer Age-based PheWAS; https://mulongdu.shinyapps.io/proap ), EOPC was preferentially associated with PCa family history while LOPC was prone to environmental and lifestyles exposures. CONCLUSIONS: This study comprehensively profiled the distinct genetic and phenotypic architecture of EOPC. The EOPC-PRS may optimize risk estimate of PCa in young males, particularly those without family history, thus providing guidance for precision population stratification.

Phenome-wide association study on miRNA-related sequence variants: the UK Biobank.

BACKGROUND: Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown. METHODS: Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis. RESULTS: We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia. CONCLUSIONS: Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Associations between attention-deficit hyperactivity disorder genetic liability and ICD-10 medical conditions in adults: utilizing electronic health records in a Phenome-Wide Association Study.

BACKGROUND: Attention-deficit hyperactivity disorder (ADHD) is often comorbid with other medical conditions in adult patients. However, ADHD is extremely underdiagnosed in adults and little is known about the medical comorbidities in undiagnosed adult individuals with high ADHD liability. In this study we investigated associations between ADHD genetic liability and electronic health record (EHR)-based ICD-10 diagnoses across all diagnostic categories, in individuals without ADHD diagnosis history. METHODS: We used data from the Estonian Biobank cohort (N = 111 261) and generated polygenic risk scores (PRS) for ADHD (PRSADHD) based on the ADHD genome-wide association study. We performed a phenome-wide association study (PheWAS) to test for associations between standardized PRSADHD and 1515 EHR-based ICD-10 diagnoses in the full and sex-stratified sample. We compared the observed significant ICD-10 associations to associations with (1) ADHD diagnosis and (2) questionnaire-based high ADHD risk analyses. RESULTS: After Bonferroni correction (p = 3.3 × 10-5) we identified 80 medical conditions associated with PRSADHD. The strongest evidence was seen with chronic obstructive pulmonary disease (OR 1.15, CI 1.11-1.18), obesity (OR 1.13, CI 1.11-1.15), and type 2 diabetes (OR 1.11, CI 1.09-1.14). Sex-stratified analysis generally showed similar associations in males and females. Out of all identified associations, 40% and 78% were also observed using ADHD diagnosis or questionnaire-based ADHD, respectively, as the predictor. CONCLUSIONS: Overall our findings indicate that ADHD genetic liability is associated with an increased risk of a substantial number of medical conditions in undiagnosed individuals. These results highlight the need for timely detection and improved management of ADHD symptoms in adults.

Disease trajectory of high neuroticism and the relevance to psychiatric disorders: A retro-prospective cohort study.

BACKGROUND: Neuroticism is a psychological personality trait that has a significant impact on public health and is also a potential predisposing factor for adverse disease outcomes; however, comprehensive studies of the subsequently developed conditions are lacking. The starting point of disease trajectory in terms of genetic variation remains unclear. METHOD: Our study included 344,609 adult participants from the UK Biobank cohort who were virtually followed up from January 1, 1997. Neuroticism levels were assessed using 12 items from the Eysenck Personality Questionnaire. We performed a phenome-wide association analysis of neuroticism and subsequent diseases. Binomial tests and logistic regression models were used to test the temporal directionality and association between disease pairs to construct disease trajectories. We also investigated the association between polygenic risk scores (PRSs) for five psychiatric traits and high neuroticism. RESULTS: The risk for 59 diseases was significantly associated with high neuroticism. Depression, anxiety, irritable bowel syndrome, migraine, spondylosis, and sleep disorders were the most likely to develop, with hazard ratios of 6.13, 3.66, 2.28, 1.74, 1.74, and 1.71, respectively. The disease trajectory network revealed two major disease clusters: cardiometabolic and chronic inflammatory diseases. Medium/high genetic risk groups stratified by the PRSs of four psychiatric traits were associated with an elevated risk of high neuroticism. We further identified eight complete phenotypic trajectory clusters of medium or high genetic risk for psychotic, anxiety-, depression-, and stress-related disorders. CONCLUSION: Neuroticism plays an important role in the development of somatic and mental disorders. The full picture of disease trajectories from the genetic risk of psychiatric traits and neuroticism in early life to a series of diseases later provides evidence for future research to explore the etiological mechanisms and precision management.

A phenome-wide association and factorial Mendelian randomization study on the repurposing of uric acid-lowering drugs for cardiovascular outcomes.

Uric acid has been linked to various disease outcomes. However, it remains unclear whether uric acid-lowering therapy could be repurposed as a treatment for conditions other than gout. We first performed both observational phenome-wide association study (Obs-PheWAS) and polygenic risk score PheWAS (PRS-PheWAS) to identify associations of uric acid levels with a wide range of disease outcomes. Then, trajectory analysis was conducted to explore temporal progression patterns of the observed disease outcomes. Finally, we investigated whether uric acid-lowering drugs could be repurposed using a factorial Mendelian randomization (MR) study design. A total of 41 overlapping phenotypes associated with uric acid levels were identified by both Obs- and PRS- PheWASs, primarily cardiometabolic diseases. The trajectory analysis illustrated how elevated uric acid levels contribute to cardiometabolic diseases, and finally death. Meanwhile, we found that uric acid-lowering drugs exerted a protective role in reducing the risk of coronary atherosclerosis (OR = 0.96, 95%CI: 0.93, 1.00, P = 0.049), congestive heart failure (OR = 0.64, 95%CI: 0.42, 0.99, P = 0.043), occlusion of cerebral arteries (OR = 0.93, 95%CI: 0.87, 1.00, P = 0.044) and peripheral vascular disease (OR = 0.60, 95%CI: 0.38, 0.94, P = 0.025). Furthermore, the combination of uric acid-lowering therapy (e.g. xanthine oxidase inhibitors) with antihypertensive treatment (e.g. calcium channel blockers) exerted additive effects and was associated with a 6%, 8%, 8%, 10% reduction in risk of coronary atherosclerosis, heart failure, occlusion of cerebral arteries and peripheral vascular disease, respectively. Our findings support a role of elevated uric acid levels in advancing cardiovascular dysfunction and identify potential repurposing opportunities for uric acid-lowering drugs in cardiovascular treatment.