Global Biobank Meta-analysis Initiative: Powering genetic discovery across human disease

Principles and methods for transferring polygenic risk scores across global populations

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Open Science Practices in Psychiatric Genetics: A Primer

Open science ensures that research is transparently reported and freely accessible for all to assess and collaboratively build on. Psychiatric genetics has led among the health sciences in implementing some open science practices in common study designs, such as replication as part of genome-wide association studies. However, thorough open science implementation guidelines are limited and largely not specific to data, privacy, and research conduct challenges in psychiatric genetics. Here, we present a primer of open science practices, including selection of a research topic with patients/nonacademic collaborators, equitable authorship and citation practices, design of replicable, reproducible studies, preregistrations, open data, and privacy issues. We provide tips for informative figures and inclusive, precise reporting. We discuss considerations in working with nonacademic collaborators and distributing research through preprints, blogs, social media, and accessible lecture materials. Finally, we provide extra resources to support every step of the research process.

Prediction and stratification of longitudinal risk for chronic obstructive pulmonary disease across smoking behaviors

Smoking is the leading risk factor for chronic obstructive pulmonary disease (COPD) worldwide, yet many people who never smoke develop COPD. We perform a longitudinal analysis of COPD in the UK Biobank to derive and validate the Socioeconomic and Environmental Risk Score which captures additive and cumulative environmental, behavioral, and socioeconomic exposure risks beyond tobacco smoking. The Socioeconomic and Environmental Risk Score is more predictive of COPD than smoking status and pack-years. Individuals in the highest decile of the risk score have a greater risk for incident COPD compared to the remaining population. Never smokers in the highest decile of exposure risk are more likely to develop COPD than previous and current smokers in the lowest decile. In general, the prediction accuracy of the Social and Environmental Risk Score is lower in non-European populations. While smoking status is often considered in screening COPD, our finding highlights the importance of other non-smoking environmental and socioeconomic variables.

Analysis across Taiwan Biobank, Biobank Japan, and UK Biobank identifies hundreds of novel loci for 36 quantitative traits

Genome-wide association studies (GWASs) have identified tens of thousands of genetic loci associated with human complex traits. However, the majority of GWASs were conducted in individuals of European ancestries. Failure to capture global genetic diversity has limited genomic discovery and has impeded equitable delivery of genomic knowledge to diverse populations. Here we report findings from 102,900 individuals across 36 human quantitative traits in the Taiwan Biobank (TWB), a major biobank effort that broadens the population diversity of genetic studies in East Asia. We identified 968 novel genetic loci, pinpointed novel causal variants through statistical fine-mapping, compared the genetic architecture across TWB, Biobank Japan, and UK Biobank, and evaluated the utility of cross-phenotype, cross-population polygenic risk scores in disease risk prediction. These results demonstrated the potential to advance discovery through diversifying GWAS populations and provided insights into the common genetic basis of human complex traits in East Asia.

Russian Biodiversity Collections: A Professional Opinion Survey

Biodiversity collections are important vehicles for protecting endangered wildlife in situations of adverse anthropogenic influence. In Russia, there are currently a number of institution- and museum-based biological collections, but there are no nation-wide centres of biodiversity collections. In this paper, we report on the results of our survey of 324 bioconservation, big-data, and ecology specialists from different regions of Russia in regard to the necessity to create several large national biodiversity centres of wildlife protection. The survey revealed specific goals that have to be fulfilled during the development of these centres for the protection and restoration of endangered wildlife species. The top three problems/tasks (topics) are the following: (1) the necessity to create large national centres for different types of specimens; (2) the full sequencing and creation of different "omic" (genomic, proteomic, transcriptomic, etc.) databases; (3) full digitisation of a biodiversity collection/centre. These goals may constitute a guideline for the future of biodiversity collections in Russia that would be targeted at protecting and restoring endangered species. With the due network service level, the translation of the website into English, and permission from the regulator (Ministry of Science and Higher Education of Russian Federation), it can also become an international project.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Mendelian randomization suggests a causal relationship between gut dysbiosis and thyroid cancer

Background

Alterations in gut microbiota composition and function have been linked to the development and progression of thyroid cancer (TC). However, the exact nature of the causal relationship between them remains uncertain.

Methods

A bidirectional two-sample Mendelian randomization (TSMR) analysis was conducted to assess the causal connection between gut microbiota (18,340 individuals) and TC (6,699 cases combined with 1,613,655 controls) using data from a genome-wide association study (GWAS). The primary analysis used the inverse-variance weighted (IVW) method to estimate the causal effect, with supplementary approaches including the weighted median, weighted mode, simple mode, and MR-Egger. Heterogeneity and pleiotropy were assessed using the Cochrane Q test, MR-Egger intercept test, and MR-PRESSO global test. A reverse TSMR analysis was performed to explore reverse causality.

Results

This study identified seven microbial taxa with significant associations with TC. Specifically, the genus Butyrivibrio (OR: 1.127, 95% CI: 1.008-1.260, p = 0.036), Fusicatenibacter (OR: 1.313, 95% CI: 1.066-1.618, p = 0.011), Oscillospira (OR: 1.240, 95% CI: 1.001-1.536, p = 0.049), Ruminococcus2 (OR: 1.408, 95% CI: 1.158-1.711, p < 0.001), Terrisporobacter (OR: 1.241, 95% CI: 1.018-1.513, p = 0.032) were identified as risk factors for TC, while The genus Olsenella (OR: 0.882, 95% CI: 0.787-0.989, p = 0.031) and Ruminococcaceae UCG004 (OR: 0.719, 95% CI: 0.566-0.914, p = 0.007) were associated with reduced TC risk. The reverse MR analysis found no evidence of reverse causality and suggested that TC may lead to increased levels of the genus Holdemanella (β: 0.053, 95% CI: 0.012~0.094, p = 0.011) and decreased levels of the order Bacillales (β: -0.075, 95% CI: -0.143~-0.006, p = 0.033). No significant bias, heterogeneity, or pleiotropy was detected in this study.

Conclusion

This study suggests a potential causal relationship between gut microbiota and TC, providing new insights into the role of gut microbiota in TC. Further research is needed to explore the underlying biological mechanisms.

Genetic risk prediction in Hispanics/Latinos: milestones, challenges, and social-ethical considerations

Genome-wide association studies (GWAS) have allowed the identification of disease-associated variants, which can be leveraged to build polygenic scores (PGSs). Even though PGSs can be a valuable tool in personalized medicine, their predictive power is limited in populations of non-European ancestry, particularly in admixed populations. Recent efforts have focused on increasing racial and ethnic diversity in GWAS, thus, addressing some of the limitations of genetic risk prediction in these populations. Even with these efforts, few studies focus exclusively on Hispanics/Latinos. Additionally, Hispanic/Latino populations are often considered a single population despite varying admixture proportions between and within ethnic groups, diverse genetic heterogeneity, and demographic history. Combined with highly heterogeneous environmental and socioeconomic exposures, this diversity can reduce the transferability of genetic risk prediction models. Given the recent increase of genomic studies that include Hispanics/Latinos, we review the milestones and efforts that focus on genetic risk prediction, summarize the potential for improving PGS transferability, and highlight the challenges yet to be addressed. Additionally, we summarize social-ethical considerations and provide ideas to promote genetic risk prediction models that can be implemented equitably.

Interactive Association Between Gut Microbiota and Thyroid Cancer

CONTEXT

The association between the gut microbiota and thyroid cancer remains controversial.

OBJECTIVE

We aimed to systematically investigate the interactive causal relationships between the abundance and metabolism pathways of gut microbiota and thyroid cancer.

METHODS

We leveraged genome-wide association studies for the abundance of 211 microbiota taxa from the MiBioGen study (N = 18 340), 205 microbiota metabolism pathways from the Dutch Microbiome Project (N = 7738), and thyroid cancer from the Global Biobank Meta-analysis Initiative (N cases = 6699 and N participants = 1 620 354). We performed a bidirectional Mendelian randomization (MR) to investigate the causality from microbiota taxa and metabolism pathways to thyroid cancer and vice versa. We performed a systematic review of previous observational studies and compared MR results with observational findings.

RESULTS

Eight taxa and 12 metabolism pathways had causal effects on thyroid cancer, where RuminococcaceaeUCG004 genus (P = .001), Streptococcaceae family (P = .016), Olsenella genus (P = .029), ketogluconate metabolism pathway (P = .003), pentose phosphate pathway (P = .016), and L-arginine degradation II in the AST pathway (P = .0007) were supported by sensitivity analyses. Conversely, thyroid cancer had causal effects on 3 taxa and 2 metabolism pathways, where the Holdemanella genus (P = .015) was supported by sensitivity analyses. The Proteobacteria phylum, Streptococcaceae family, Ruminococcus2 genus, and Holdemanella genus were significantly associated with thyroid cancer in both the systematic review and MR, whereas the other 121 significant taxa in observational results were not supported by MR.

DISCUSSIONS

These findings implicated the potential role of host-microbiota crosstalk in thyroid cancer, while the discrepancy among observational studies calls for further investigations.

New Insights into Polygenic Score-Lifestyle Interactions for Cardiometabolic Risk Factors from Genome-Wide Interaction Analyses

The relationship between lifestyles and cardiometabolic outcomes varies between individuals. In 382,275 UK Biobank Europeans, we investigate how lifestyles interact with polygenic scores (PGS) of cardiometabolic risk factors. We identify six interactions (PGS for body mass index with meat diet, physical activity, sedentary behaviour and insomnia; PGS for high-density lipoprotein cholesterol with sedentary behaviour; PGS for triglycerides with meat diet) in multivariable linear regression models including an interaction term and show stronger associations between lifestyles and cardiometabolic risk factors among individuals with high PGSs than those with low PGSs. Genome-wide interaction analyses pinpoint three genetic variants (FTO rs72805613 for BMI; CETP rs56228609 for high-density lipoprotein cholesterol; TRIB2 rs4336630 for triglycerides; PInteraction < 5 × 10-8). The associations between lifestyles and cardiometabolic risk factors differ between individuals grouped by the genotype of these variants, with the degree of differences being similar to that between individuals with high and low values for the corresponding PGSs. This study demonstrates that associations between lifestyles and cardiometabolic risk factors can differ between individuals based upon their genetic profiles. It further suggests that genetic variants with interaction effects contribute more to such differences compared to those without interaction effects, which has potential implications for developing PGSs for personalised intervention.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Next-generation phenotyping: introducing phecodeX for enhanced discovery research in medical phenomics

MOTIVATION

Phecodes are widely used and easily adapted phenotypes based on International Classification of Diseases codes. The current version of phecodes (v1.2) was designed primarily to study common/complex diseases diagnosed in adults; however, there are numerous limitations in the codes and their structure.

RESULTS

Here, we present phecodeX, an expanded version of phecodes with a revised structure and 1,761 new codes. PhecodeX adds granularity to phenotypes in key disease domains that are under-represented in the current phecode structure-including infectious disease, pregnancy, congenital anomalies, and neonatology-and is a more robust representation of the medical phenome for global use in discovery research.

AVAILABILITY AND IMPLEMENTATION

phecodeX is available at https://github.com/PheWAS/phecodeX.

Polygenic risk scores for disease risk prediction in Africa: current challenges and future directions

Early identification of genetic risk factors for complex diseases can enable timely interventions and prevent serious outcomes, including mortality. While the genetics underlying many Mendelian diseases have been elucidated, it is harder to predict risk for complex diseases arising from the combined effects of many genetic variants with smaller individual effects on disease aetiology. Polygenic risk scores (PRS), which combine multiple contributing variants to predict disease risk, have the potential to influence the implementation for precision medicine. However, the majority of existing PRS were developed from European data with limited transferability to African populations. Notably, African populations have diverse genetic backgrounds, and a genomic architecture with smaller haplotype blocks compared to European genomes. Subsequently, growing evidence shows that using large-scale African ancestry cohorts as discovery for PRS development may generate more generalizable findings. Here, we (1) discuss the factors contributing to the poor transferability of PRS in African populations, (2) showcase the novel Africa genomic datasets for PRS development, (3) explore the potential clinical utility of PRS in African populations, and (4) provide insight into the future of PRS in Africa.

Decoding Genetics, Ancestry, and Geospatial Context for Precision Health

Mass General Brigham, an integrated healthcare system based in the Greater Boston area of Massachusetts, annually serves 1.5 million patients. We established the Mass General Brigham Biobank (MGBB), encompassing 142,238 participants, to unravel the intricate relationships among genomic profiles, environmental context, and disease manifestations within clinical practice. In this study, we highlight the impact of ancestral diversity in the MGBB by employing population genetics, geospatial assessment, and association analyses of rare and common genetic variants. The population structures captured by the genetics mirror the sequential immigration to the Greater Boston area throughout American history, highlighting communities tied to shared genetic and environmental factors. Our investigation underscores the potency of unbiased, large-scale analyses in a healthcare-affiliated biobank, elucidating the dynamic interplay across genetics, immigration, structural geospatial factors, and health outcomes in one of the earliest American sites of European colonization.

Causal inference of sex hormone-binding globulin on venous thromboembolism: evidence from Mendelian randomisation

BACKGROUND

Previous cohort studies have shown that exogenous sex hormone use, such as testosterone replacement therapy and oestrogen-containing contraceptives, can increase the risk of venous thromboembolism (VTE). However, the relationship between endogenous sex hormone levels and VTE remains unclear. The goal of the present study was to explore the causal roles of endogenous sex hormones, including hormone-binding globulin (SHBG), bioactive testosterone (BT), and total testosterone (TT), in VTE and its two subgroups, deep vein thrombosis (DVT) and pulmonary embolism (PE).

METHODS

We used a genome-wide association study of sex hormones as exposure data and Finnish VTE data as the outcome. Inverse variance weighting, MR-Egger, and weighted median were used for two-sample Mendelian randomisation (MR). Sensitivity analyses included MR-Egger, MR-PRESSO, Cochrane Q test, MR Steiger, leave-one-out analysis, and funnel plot, combined with multivariate MR and replicated MR analyses using larger VTE data from the global biobank meta-analysis initiative. Linkage disequilibrium score regression (LDSC) was used to determine genetic associations and estimate sample overlap.

RESULTS

Our findings genetically predicted that an increase in serum SHBG levels by one standard deviation (SD) caused 25% higher odds for VTE (OR: 1.25, 95% CI: 1.01-1.55) and 58% higher odds for PE (OR: 1.58, 95% CI: 1.20-2.08). LDSC supported the genetic correlation between these two traits and replicated analyses confirm SHBG's genetic effect on VTE in both sexes (OR: 1.46, 95% CI: 1.20-1.78) and in females (OR: 1.49, 95% CI: 1.17-1.91). In addition, an increase in serum TT levels by one SD caused 32% higher odds for VTE (OR: 1.32, 95% CI: 1.08-1.62) and 31% higher odds for DVT (OR: 1.31, 95% CI: 1.01-1.69); however, LDSC and replicated analyses did not find a genetic correlation between TT and VTE or its subtypes. No significant correlation was observed between BT and all three outcome traits.

CONCLUSION

Our study provides evidence that elevated serum SHBG levels, as predicted by genetics, increase VTE risk. However, the causal effect of testosterone levels on VTE requires further investigation.

Revealing polygenic pleiotropy using genetic risk scores for asthma

In this study we examined how genetic risk for asthma associates with different features of the disease and with other medical conditions and traits. Using summary statistics from two multi-ancestry genome-wide association studies of asthma, we modeled polygenic risk scores (PRSs) and validated their predictive performance in the UK Biobank. We then performed phenome-wide association studies of the asthma PRSs with 371 heritable traits in the UK Biobank. We identified 228 total significant associations across a variety of organ systems, including associations that varied by PRS model, sex, age of asthma onset, ancestry, and human leukocyte antigen region alleles. Our results highlight pervasive pleiotropy between asthma and numerous other traits and conditions and elucidate pathways that contribute to asthma and its comorbidities.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Detecting biased validation of predictive models in the positive-unlabeled setting: disease gene prioritization case study

Motivation

Positive-unlabeled data consists of points with either positive or unknown labels. It is widespread in medical, genetic, and biological settings, creating a high demand for predictive positive-unlabeled models. The performance of such models is usually estimated using validation sets, assumed to be selected completely at random (SCAR) from known positive examples. For certain metrics, this assumption enables unbiased performance estimation when treating positive-unlabeled data as positive/negative. However, the SCAR assumption is often adopted without proper justifications, simply for the sake of convenience.

Results

We provide an algorithm that under the weak assumptions of a lower bound on the number of positive examples can test for the violation of the SCAR assumption. Applying it to the problem of gene prioritization for complex genetic traits, we illustrate that the SCAR assumption is often violated there, causing the inflation of performance estimates, which we refer to as validation bias. We estimate the potential impact of validation bias on performance estimation. Our analysis reveals that validation bias is widespread in gene prioritization data and can significantly overestimate the performance of models. This finding elucidates the discrepancy between the reported good performance of models and their limited practical applications.

Availability and implementation

Python code with examples of application of the validation bias detection algorithm is available at github.com/ArtomovLab/ValidationBias.

Demonstrating paths for unlocking the value of cloud genomics through cross cohort analysis

Recently, large scale genomic projects such as All of Us and the UK Biobank have introduced a new research paradigm where data are stored centrally in cloud-based Trusted Research Environments (TREs). To characterize the advantages and drawbacks of different TRE attributes in facilitating cross-cohort analysis, we conduct a Genome-Wide Association Study of standard lipid measures using two approaches: meta-analysis and pooled analysis. Comparison of full summary data from both approaches with an external study shows strong correlation of known loci with lipid levels (R2 ~ 83-97%). Importantly, 90 variants meet the significance threshold only in the meta-analysis and 64 variants are significant only in pooled analysis, with approximately 20% of variants in each of those groups being most prevalent in non-European, non-Asian ancestry individuals. These findings have important implications, as technical and policy choices lead to cross-cohort analyses generating similar, but not identical results, particularly for non-European ancestral populations.

Extremely sparse models of linkage disequilibrium in ancestrally diverse association studies

Linkage disequilibrium (LD) is the correlation among nearby genetic variants. In genetic association studies, LD is often modeled using large correlation matrices, but this approach is inefficient, especially in ancestrally diverse studies. In the present study, we introduce LD graphical models (LDGMs), which are an extremely sparse and efficient representation of LD. LDGMs are derived from genome-wide genealogies; statistical relationships among alleles in the LDGM correspond to genealogical relationships among haplotypes. We published LDGMs and ancestry-specific LDGM precision matrices for 18 million common variants (minor allele frequency >1%) in five ancestry groups, validated their accuracy and demonstrated order-of-magnitude improvements in runtime for commonly used LD matrix computations. We implemented an extremely fast multiancestry polygenic prediction method, BLUPx-ldgm, which performs better than a similar method based on the reference LD correlation matrix. LDGMs will enable sophisticated methods that scale to ancestrally diverse genetic association data across millions of variants and individuals.

Interplay of chronic obstructive pulmonary disease and colorectal cancer development: unravelling the mediating role of fatty acids through a comprehensive multi-omics analysis

BACKGROUND

Chronic obstructive pulmonary disease (COPD) patients often exhibit gastrointestinal symptoms, A potential association between COPD and Colorectal Cancer (CRC) has been indicated, warranting further examination.

METHODS

In this study, we collected COPD and CRC data from the National Health and Nutrition Examination Survey, genome-wide association studies, and RNA sequence for a comprehensive analysis. We used weighted logistic regression to explore the association between COPD and CRC incidence risk. Mendelian randomization analysis was performed to assess the causal relationship between COPD and CRC, and cross-phenotype meta-analysis was conducted to pinpoint crucial loci. Multivariable mendelian randomization was used to uncover mediating factors connecting the two diseases. Our results were validated using both NHANES and GEO databases.

RESULTS

In our analysis of the NHANES dataset, we identified COPD as a significant contributing factor to CRC development. MR analysis revealed that COPD increased the risk of CRC onset and progression (OR: 1.16, 95% CI 1.01-1.36). Cross-phenotype meta-analysis identified four critical genes associated with both CRC and COPD. Multivariable Mendelian randomization suggested body fat percentage, omega-3, omega-6, and the omega-3 to omega-6 ratio as potential mediating factors for both diseases, a finding consistent with the NHANES dataset. Further, the interrelation between fatty acid-related modules in COPD and CRC was demonstrated via weighted gene co-expression network analysis and Kyoto Encyclopedia of Genes and Genomes enrichment results using RNA expression data.

CONCLUSIONS

This study provides novel insights into the interplay between COPD and CRC, highlighting the potential impact of COPD on the development of CRC. The identification of shared genes and mediating factors related to fatty acid metabolism deepens our understanding of the underlying mechanisms connecting these two diseases.

Whole Person Modeling: a transdisciplinary approach to mental health research

The growing global burden of mental illness has prompted calls for innovative research strategies. Theoretical models of mental health include complex contributions of biological, psychosocial, experiential, and other environmental influences. Accordingly, neuropsychiatric research has self-organized into largely isolated disciplines working to decode each individual contribution. However, research directly modeling objective biological measurements in combination with cognitive, psychological, demographic, or other environmental measurements is only now beginning to proliferate. This review aims to (1) to describe the landscape of modern mental health research and current movement towards integrative study, (2) to provide a concrete framework for quantitative integrative research, which we call Whole Person Modeling, (3) to explore existing and emerging techniques and methods used in Whole Person Modeling, and (4) to discuss our observations about the scarcity, potential value, and untested aspects of highly transdisciplinary research in general. Whole Person Modeling studies have the potential to provide a better understanding of multilevel phenomena, deliver more accurate diagnostic and prognostic tests to aid in clinical decision making, and test long standing theoretical models of mental illness. Some current barriers to progress include challenges with interdisciplinary communication and collaboration, systemic cultural barriers to transdisciplinary career paths, technical challenges in model specification, bias, and data harmonization, and gaps in transdisciplinary educational programs. We hope to ease anxiety in the field surrounding the often mysterious and intimidating world of transdisciplinary, data-driven mental health research and provide a useful orientation for students or highly specialized researchers who are new to this area.

The All of Us Data and Research Center: Creating a Secure, Scalable, and Sustainable Ecosystem for Biomedical Research

The All of Us Research Program's Data and Research Center (DRC) was established to help acquire, curate, and provide access to one of the world's largest and most diverse datasets for precision medicine research. Already, over 500,000 participants are enrolled in All of Us, 80% of whom are underrepresented in biomedical research, and data are being analyzed by a community of over 2,300 researchers. The DRC created this thriving data ecosystem by collaborating with engaged participants, innovative program partners, and empowered researchers. In this review, we first describe how the DRC is organized to meet the needs of this broad group of stakeholders. We then outline guiding principles, common challenges, and innovative approaches used to build the All of Us data ecosystem. Finally, we share lessons learned to help others navigate important decisions and trade-offs in building a modern biomedical data platform.

Addressing the Challenge of Biomedical Data Inequality: An Artificial Intelligence Perspective

Artificial intelligence (AI) and other data-driven technologies hold great promise to transform healthcare and confer the predictive power essential to precision medicine. However, the existing biomedical data, which are a vital resource and foundation for developing medical AI models, do not reflect the diversity of the human population. The low representation in biomedical data has become a significant health risk for non-European populations, and the growing application of AI opens a new pathway for this health risk to manifest and amplify. Here we review the current status of biomedical data inequality and present a conceptual framework for understanding its impacts on machine learning. We also discuss the recent advances in algorithmic interventions for mitigating health disparities arising from biomedical data inequality. Finally, we briefly discuss the newly identified disparity in data quality among ethnic groups and its potential impacts on machine learning.

Importance of Diversity in Precision Medicine: Generalizability of Genetic Associations Across Ancestry Groups Toward Better Identification of Disease Susceptibility Variants

Genome-wide association studies (GWAS) revolutionized our understanding of common genetic variation and its impact on common human disease and traits. Developed and adopted in the mid-2000s, GWAS led to searchable genotype-phenotype catalogs and genome-wide datasets available for further data mining and analysis for the eventual development of translational applications. The GWAS revolution was swift and specific, including almost exclusively populations of European descent, to the neglect of the majority of the world's genetic diversity. In this narrative review, we recount the GWAS landscape of the early years that established a genotype-phenotype catalog that is now universally understood to be inadequate for a complete understanding of complex human genetics. We then describe approaches taken to augment the genotype-phenotype catalog, including the study populations, collaborative consortia, and study design approaches aimed to generalize and then ultimately discover genome-wide associations in non-European descent populations. The collaborations and data resources established in the efforts to diversify genomic findings undoubtedly provide the foundations of the next chapters of genetic association studies with the advent of budget-friendly whole-genome sequencing.

A FinnGen pilot clinical recall study for Alzheimer's disease

Successful development of novel therapies requires that clinical trials are conducted in patient cohorts with the highest benefit-to-risk ratio. Population-based biobanks with comprehensive health and genetic data from large numbers of individuals hold promise to facilitate identification of trial participants, particularly when interventions need to start while symptoms are still mild, such as for Alzheimer's disease (AD). This study describes a process for clinical recall studies from FinnGen. We demonstrate the feasibility to systematically ascertain customized clinical data from FinnGen participants with ICD10 diagnosis of AD or mild cognitive disorder (MCD) in a single-center cross-sectional study testing blood-based biomarkers and cognitive functioning in-person, computer-based and remote. As a result, 19% (27/140) of a pre-specified FinnGen subcohort were successfully recalled and completed the study. Hospital records largely validated registry entries. For 8/12 MCD patients, other reasons than AD were identified as underlying diagnosis. Cognitive measures correlated across platforms, with highest consistencies for dementia screening (r = 0.818) and semantic fluency (r = 0.764), respectively, for in-person versus telephone-administered tests. Glial fibrillary acidic protein (GFAP) (p < 0.002) and phosphorylated-tau 181 (pTau-181) (p < 0.020) most reliably differentiated AD from MCD participants. We conclude that informative, customized clinical recall studies from FinnGen are feasible.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

The Impact of Genomic Variation on Function (IGVF) Consortium

Our genomes influence nearly every aspect of human biology from molecular and cellular functions to phenotypes in health and disease. Human genetics studies have now associated hundreds of thousands of differences in our DNA sequence ("genomic variation") with disease risk and other phenotypes, many of which could reveal novel mechanisms of human biology and uncover the basis of genetic predispositions to diseases, thereby guiding the development of new diagnostics and therapeutics. Yet, understanding how genomic variation alters genome function to influence phenotype has proven challenging. To unlock these insights, we need a systematic and comprehensive catalog of genome function and the molecular and cellular effects of genomic variants. Toward this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations, and predictive modeling to investigate the relationships among genomic variation, genome function, and phenotypes. Through systematic comparisons and benchmarking of experimental and computational methods, we aim to create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how both coding and noncoding variants may connect through gene regulatory and protein interaction networks. These experimental data, computational predictions, and accompanying standards and pipelines will be integrated into an open resource that will catalyze community efforts to explore genome function and the impact of genetic variation on human biology and disease across populations.

Genome-wide association analysis and Mendelian randomization proteomics identify drug targets for heart failure

We conduct a large-scale meta-analysis of heart failure genome-wide association studies (GWAS) consisting of over 90,000 heart failure cases and more than 1 million control individuals of European ancestry to uncover novel genetic determinants for heart failure. Using the GWAS results and blood protein quantitative loci, we perform Mendelian randomization and colocalization analyses on human proteins to provide putative causal evidence for the role of druggable proteins in the genesis of heart failure. We identify 39 genome-wide significant heart failure risk variants, of which 18 are previously unreported. Using a combination of Mendelian randomization proteomics and genetic cis-only colocalization analyses, we identify 10 additional putatively causal genes for heart failure. Findings from GWAS and Mendelian randomization-proteomics identify seven (CAMK2D, PRKD1, PRKD3, MAPK3, TNFSF12, APOC3 and NAE1) proteins as potential targets for interventions to be used in primary prevention of heart failure.

Fine-mapping across diverse ancestries drives the discovery of putative causal variants underlying human complex traits and diseases

Genome-wide association studies (GWAS) of human complex traits or diseases often implicate genetic loci that span hundreds or thousands of genetic variants, many of which have similar statistical significance. While statistical fine-mapping in individuals of European ancestries has made important discoveries, cross-population fine-mapping has the potential to improve power and resolution by capitalizing on the genomic diversity across ancestries. Here we present SuSiEx, an accurate and computationally efficient method for cross-population fine-mapping, which builds on the single-population fine-mapping framework, Sum of Single Effects (SuSiE). SuSiEx integrates data from an arbitrary number of ancestries, explicitly models population-specific allele frequencies and LD patterns, accounts for multiple causal variants in a genomic region, and can be applied to GWAS summary statistics. We comprehensively evaluated SuSiEx using simulations, a range of quantitative traits measured in both UK Biobank and Taiwan Biobank, and schizophrenia GWAS across East Asian and European ancestries. In all evaluations, SuSiEx fine-mapped more association signals, produced smaller credible sets and higher posterior inclusion probability (PIP) for putative causal variants, and captured population-specific causal variants.

A multi-ancestry polygenic risk score improves risk prediction for coronary artery disease

Identification of individuals at highest risk of coronary artery disease (CAD)-ideally before onset-remains an important public health need. Prior studies have developed genome-wide polygenic scores to enable risk stratification, reflecting the substantial inherited component to CAD risk. Here we develop a new and significantly improved polygenic score for CAD, termed GPSMult, that incorporates genome-wide association data across five ancestries for CAD (>269,000 cases and >1,178,000 controls) and ten CAD risk factors. GPSMult strongly associated with prevalent CAD (odds ratio per standard deviation 2.14, 95% confidence interval 2.10-2.19, P < 0.001) in UK Biobank participants of European ancestry, identifying 20.0% of the population with 3-fold increased risk and conversely 13.9% with 3-fold decreased risk as compared with those in the middle quintile. GPSMult was also associated with incident CAD events (hazard ratio per standard deviation 1.73, 95% confidence interval 1.70-1.76, P < 0.001), identifying 3% of healthy individuals with risk of future CAD events equivalent to those with existing disease and significantly improving risk discrimination and reclassification. Across multiethnic, external validation datasets inclusive of 33,096, 124,467, 16,433 and 16,874 participants of African, European, Hispanic and South Asian ancestry, respectively, GPSMult demonstrated increased strength of associations across all ancestries and outperformed all available previously published CAD polygenic scores. These data contribute a new GPSMult for CAD to the field and provide a generalizable framework for how large-scale integration of genetic association data for CAD and related traits from diverse populations can meaningfully improve polygenic risk prediction.

Diversity and Scale: Genetic Architecture of 2,068 Traits in the VA Million Veteran Program

Genome-wide association studies (GWAS) have underrepresented individuals from non-European populations, impeding progress in characterizing the genetic architecture and consequences of health and disease traits. To address this, we present a population-stratified phenome-wide GWAS followed by a multi-population meta-analysis for 2,068 traits derived from electronic health records of 635,969 participants in the Million Veteran Program (MVP), a longitudinal cohort study of diverse U.S. Veterans genetically similar to the respective African (121,177), Admixed American (59,048), East Asian (6,702), and European (449,042) superpopulations defined by the 1000 Genomes Project. We identified 38,270 independent variants associating with one or more traits at experiment-wide P < 4.6 × 10 - 11 significance; fine-mapping 6,318 signals identified from 613 traits to single-variant resolution. Among these, a third (2,069) of the associations were found only among participants genetically similar to non-European reference populations, demonstrating the importance of expanding diversity in genetic studies. Our work provides a comprehensive atlas of phenome-wide genetic associations for future studies dissecting the architecture of complex traits in diverse populations.

The associations between gut microbiota and chronic respiratory diseases: a Mendelian randomization study

Introduction

Growing evidence indicates that variations in the composition of the gut microbiota are linked to the onset and progression of chronic respiratory diseases (CRDs), albeit the causal relationship between the two remains unclear.

Methods

We conducted a comprehensive two-sample Mendelian randomization (MR) analysis to investigate the relationship between gut microbiota and five main CRDs, including chronic obstructive pulmonary disease (COPD), asthma, idiopathic pulmonary fibrosis (IPF), sarcoidosis, and pneumoconiosis. For MR analysis, the inverse variance weighted (IVW) method was utilized as the primary method. The MR-Egger, weighted median, and MR-PRESSO statistical methods were used as a supplement. To detect heterogeneity and pleiotropy, the Cochrane and Rucker Q test, MR-Egger intercept test, and MR-PRESSO global test were then implemented. The leave-one-out strategy was also applied to assess the consistency of the MR results.

Results

Based on substantial genetic data obtained from genome-wide association studies (GWAS) comprising 3,504,473 European participants, our study offers evidence that several gut microbial taxa, including 14 probable microbial taxa (specifically, 5, 3, 2, 3 and 1 for COPD, asthma, IPF, sarcoidosis, and pneumoconiosis, respectively) and 33 possible microbial taxa (specifically, 6, 7, 8, 7 and 5 for COPD, asthma, IPF, sarcoidosis, and pneumoconiosis, respectively) play significant roles in the formation of CRDs.

Discussion

This work implies causal relationships between the gut microbiota and CRDs, thereby shedding new light on the gut microbiota-mediated prevention of CRDs.

Impact of cross-ancestry genetic architecture on GWASs in admixed populations

Genome-wide association studies (GWASs) have identified thousands of variants for disease risk. These studies have predominantly been conducted in individuals of European ancestries, which raises questions about their transferability to individuals of other ancestries. Of particular interest are admixed populations, usually defined as populations with recent ancestry from two or more continental sources. Admixed genomes contain segments of distinct ancestries that vary in composition across individuals in the population, allowing for the same allele to induce risk for disease on different ancestral backgrounds. This mosaicism raises unique challenges for GWASs in admixed populations, such as the need to correctly adjust for population stratification. In this work we quantify the impact of differences in estimated allelic effect sizes for risk variants between ancestry backgrounds on association statistics. Specifically, while the possibility of estimated allelic effect-size heterogeneity by ancestry (HetLanc) can be modeled when performing a GWAS in admixed populations, the extent of HetLanc needed to overcome the penalty from an additional degree of freedom in the association statistic has not been thoroughly quantified. Using extensive simulations of admixed genotypes and phenotypes, we find that controlling for and conditioning effect sizes on local ancestry can reduce statistical power by up to 72%. This finding is especially pronounced in the presence of allele frequency differentiation. We replicate simulation results using 4,327 African-European admixed genomes from the UK Biobank for 12 traits to find that for most significant SNPs, HetLanc is not large enough for GWASs to benefit from modeling heterogeneity in this way.

Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens

Most variants associated with complex traits and diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown effects. Using ancestrally diverse, biobank-scale GWAS data, massively parallel CRISPR screens, and single-cell transcriptomic and proteomic sequencing, we discovered 124 cis-target genes of 91 noncoding blood trait GWAS loci. Using precise variant insertion through base editing, we connected specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs. Networks were themselves enriched for GWAS variants and demonstrated polygenic contributions to complex traits. This platform enables massively parallel characterization of the target genes and mechanisms of human noncoding variants in both cis and trans.

Family Matters: Germline Testing in Thoracic Cancers

Most thoracic cancers arise via a series of stepwise somatic alterations driven by a well-defined carcinogen (ie, tobacco or asbestos for lung cancer and mesothelioma, respectively). A small proportion can emerge on a background of pathogenic germline variants (PGVs), which have the property of heritability. In general, PGVs may be initially suspected on the basis of the presence of specific clinical features. Such gene × environment interactions significantly increase the risk of developing lung cancer (1.5- to 3.2-fold). PGVs have been discovered involving the actionable driver oncogene, epidermal growth factor receptor (EGFR), with an EGFR T790M PGV rate of 0.3%-0.9% in the nonsquamous non-small-cell lung cancer subtype. Its appearance during routine somatic DNA sequencing in those patients who have not had a previous tyrosine kinase inhibitor should raise suspicion. In patients with sporadic mesothelioma, BAP1 is the most frequently mutated tumor driver, with a PGV rate between 2.8% and 8%, associated with a favorable prognosis. BAP1 PGVs accelerate mesothelioma tumorigenesis after asbestos exposure in preclinical models and may be partly predicted by clinical criteria. At present, routine germline genetic testing for thoracic cancers is not a standard practice. Expert genetic counseling is, therefore, required for patients who carry a PGV. Ongoing studies aim to better understand the natural history of patients harboring PGVs to underpin future cancer prevention, precise counseling, and cancer management with the goal of improving the quality and length of life.

Inference of chronic obstructive pulmonary disease with deep learning on raw spirograms identifies new genetic loci and improves risk models

Chronic obstructive pulmonary disease (COPD), the third leading cause of death worldwide, is highly heritable. While COPD is clinically defined by applying thresholds to summary measures of lung function, a quantitative liability score has more power to identify genetic signals. Here we train a deep convolutional neural network on noisy self-reported and International Classification of Diseases labels to predict COPD case-control status from high-dimensional raw spirograms and use the model's predictions as a liability score. The machine-learning-based (ML-based) liability score accurately discriminates COPD cases and controls, and predicts COPD-related hospitalization without any domain-specific knowledge. Moreover, the ML-based liability score is associated with overall survival and exacerbation events. A genome-wide association study on the ML-based liability score replicates existing COPD and lung function loci and also identifies 67 new loci. Lastly, our method provides a general framework to use ML methods and medical-record-based labels that does not require domain knowledge or expert curation to improve disease prediction and genomic discovery for drug design.

Prediction and stratification of longitudinal risk for chronic obstructive pulmonary disease across smoking behaviors

Smoking is the leading risk factor for chronic obstructive pulmonary disease (COPD) worldwide, yet many people who never smoke develop COPD. We hypothesize that considering other socioeconomic and environmental factors can better predict and stratify the risk of COPD in both non-smokers and smokers. We performed longitudinal analysis of COPD in the UK Biobank to develop the Socioeconomic and Environmental Risk Score (SERS) which captures additive and cumulative environmental, behavioral, and socioeconomic exposure risks beyond tobacco smoking. We tested the ability of SERS to predict and stratify the risk of COPD in current, previous, and never smokers of European and non-European ancestries in comparison to a composite genome-wide polygenic risk score (PGS). We tested associations using Cox regression models and assessed the predictive performance of models using Harrell's C index. SERS (C index = 0.770, 95% CI 0.756 to 0.784) was more predictive of COPD than smoking status (C index = 0.738, 95% CI 0.724 to 0.752), pack-years (C index = 0.742, 95% CI 0.727 to 0.756). Compared to the remaining population, individuals in the highest decile of the SERS had hazard ratios (HR) = 7.24 (95% CI 6.51 to 8.05, P < 0.0001) for incident COPD. Never smokers in the highest decile of exposure risk were more likely to develop COPD than previous and current smokers in the lowest decile with HR=4.95 (95% CI 1.56 to 15.69, P=6.65×10-3) and 2.92 (95%CI 1.51 to 5.61, P=1.38×10-3), respectively. In general, the prediction accuracy of SERS was lower in the non-European populations compared to the European evaluation set. In addition to genetic factors, socioeconomic and environmental factors beyond smoking can predict and stratify COPD risk for both non- and smoking individuals. Smoking status is often considered in screening; other non-smoking environmental and non-genetic variables should be evaluated prospectively for their clinical utility.

Causal effects on complex traits are similar for common variants across segments of different continental ancestries within admixed individuals

Individuals of admixed ancestries (for example, African Americans) inherit a mosaic of ancestry segments (local ancestry) originating from multiple continental ancestral populations. This offers the unique opportunity of investigating the similarity of genetic effects on traits across ancestries within the same population. Here we introduce an approach to estimate correlation of causal genetic effects (radmix) across local ancestries and analyze 38 complex traits in African-European admixed individuals (N = 53,001) to observe very high correlations (meta-analysis radmix = 0.95, 95% credible interval 0.93-0.97), much higher than correlation of causal effects across continental ancestries. We replicate our results using regression-based methods from marginal genome-wide association study summary statistics. We also report realistic scenarios where regression-based methods yield inflated heterogeneity-by-ancestry due to ancestry-specific tagging of causal effects, and/or polygenicity. Our results motivate genetic analyses that assume minimal heterogeneity in causal effects by ancestry, with implications for the inclusion of ancestry-diverse individuals in studies.

Dissecting the causal effect between gut microbiota, DHA, and urate metabolism: A large-scale bidirectional Mendelian randomization

Objectives

Our aim was to investigate the interactive causal effects between gut microbiota and host urate metabolism and explore the underlying mechanism using genetic methods.

Methods

We extracted summary statistics from the abundance of 211 microbiota taxa from the MiBioGen (N =18,340), 205 microbiota metabolism pathways from the Dutch Microbiome Project (N =7738), gout from the Global Biobank Meta-analysis Initiative (N =1,448,128), urate from CKDGen (N =288,649), and replication datasets from the Global Urate Genetics Consortium (N gout =69,374; N urate =110,347). We used linkage disequilibrium score regression and bidirectional Mendelian randomization (MR) to detect genetic causality between microbiota and gout/urate. Mediation MR and colocalization were performed to investigate potential mediators in the association between microbiota and urate metabolism.

Results

Two taxa had a common causal effect on both gout and urate, whereas the Victivallaceae family was replicable. Six taxa were commonly affected by both gout and urate, whereas the Ruminococcus gnavus group genus was replicable. Genetic correlation supported significant results in MR. Two microbiota metabolic pathways were commonly affected by gout and urate. Mediation analysis indicated that the Bifidobacteriales order and Bifidobacteriaceae family had protective effects on urate mediated by increasing docosahexaenoic acid. These two bacteria shared a common causal variant rs182549 with both docosahexaenoic acid and urate, which was located within MCM6/LCT locus.

Conclusions

Gut microbiota and host urate metabolism had a bidirectional causal association, implicating the critical role of host-microbiota crosstalk in hyperuricemic patients. Changes in gut microbiota can not only ameliorate host urate metabolism but also become a foreboding indicator of urate metabolic diseases.

The Michigan Genomics Initiative: A biobank linking genotypes and electronic clinical records in Michigan Medicine patients

Biobanks of linked clinical patient histories and biological samples are an efficient strategy to generate large cohorts for modern genetics research. Biobank recruitment varies by factors such as geographic catchment and sampling strategy, which affect biobank demographics and research utility. Here, we describe the Michigan Genomics Initiative (MGI), a single-health-system biobank currently consisting of >91,000 participants recruited primarily during surgical encounters at Michigan Medicine. The surgical enrollment results in a biobank enriched for many diseases and ideally suited for a disease genetics cohort. Compared with the much larger population-based UK Biobank, MGI has higher prevalence for nearly all diagnosis-code-based phenotypes and larger absolute case counts for many phenotypes. Genome-wide association study (GWAS) results replicate known findings, thereby validating the genetic and clinical data. Our results illustrate that opportunistic biobank sampling within single health systems provides a unique and complementary resource for exploring the genetics of complex diseases.

From 'Omics to Multi-omics Technologies: the Discovery of Novel Causal Mediators

PURPOSE OF REVIEW

'Omics studies provide a comprehensive characterisation of a biological entity, such as the genome, epigenome, transcriptome, proteome, metabolome, or microbiome. This review covers the unique properties of these types of 'omics and their roles as causal mediators in cardiovascular disease. Moreover, applications and challenges of integrating multiple types of 'omics data to increase predictive power, improve causal inference, and elucidate biological mechanisms are discussed.

RECENT FINDINGS

Multi-omics approaches are growing in adoption as they provide orthogonal evidence and overcome the limitations of individual types of 'omics data. Studies with multiple types of 'omics data have improved the diagnosis and prediction of disease states and afforded a deeper understanding of underlying pathophysiological mechanisms, beyond any single type of 'omics data. For instance, disease-associated loci in the genome can be supplemented with other 'omics to prioritise causal genes and understand the function of non-coding variants. Alternatively, techniques, such as Mendelian randomisation, can leverage genetics to provide evidence supporting a causal role for disease-associated molecules, and elucidate their role in disease pathogenesis. As technologies improve, costs for 'omics studies will continue to fall and datasets will become increasingly accessible to researchers. The intrinsically unbiased nature of 'omics data is well-suited to exploratory analyses that discover causal mediators of disease, and multi-omics is an emerging discipline that leverages the strengths of each type of 'omics data to provide insights greater than the sum of its parts.

Pharmacogenomics: current status and future perspectives

Inter-individual variability in drug response, be it efficacy or safety, is common and likely to become an increasing problem globally given the growing elderly population requiring treatment. Reasons for this inter-individual variability include genomic factors, an area of study called pharmacogenomics. With genotyping technologies now widely available and decreasing in cost, implementing pharmacogenomics into clinical practice - widely regarded as one of the initial steps in mainstreaming genomic medicine - is currently a focus in many countries worldwide. However, major challenges of implementation lie at the point of delivery into health-care systems, including the modification of current clinical pathways coupled with a massive knowledge gap in pharmacogenomics in the health-care workforce. Pharmacogenomics can also be used in a broader sense for drug discovery and development, with increasing evidence suggesting that genomically defined targets have an increased success rate during clinical development.

Impact of cross-ancestry genetic architecture on GWAS in admixed populations

Genome-wide association studies (GWAS) have identified thousands of variants for disease risk. These studies have predominantly been conducted in individuals of European ancestries, which raises questions about their transferability to individuals of other ancestries. Of particular interest are admixed populations, usually defined as populations with recent ancestry from two or more continental sources. Admixed genomes contain segments of distinct ancestries that vary in composition across individuals in the population, allowing for the same allele to induce risk for disease on different ancestral backgrounds. This mosaicism raises unique challenges for GWAS in admixed populations, such as the need to correctly adjust for population stratification to balance type I error with statistical power. In this work we quantify the impact of differences in estimated allelic effect sizes for risk variants between ancestry backgrounds on association statistics. Specifically, while the possibility of estimated allelic effect-size heterogeneity by ancestry (HetLanc) can be modeled when performing GWAS in admixed populations, the extent of HetLanc needed to overcome the penalty from an additional degree of freedom in the association statistic has not been thoroughly quantified. Using extensive simulations of admixed genotypes and phenotypes we find that modeling HetLanc in its absence reduces statistical power by up to 72%. This finding is especially pronounced in the presence of allele frequency differentiation. We replicate simulation results using 4,327 African-European admixed genomes from the UK Biobank for 12 traits to find that for most significant SNPs HetLanc is not large enough for GWAS to benefit from modeling heterogeneity.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.