Polygenic transcriptome risk scores (PTRS) can improve portability of polygenic risk scores across ancestries

Susceptibility gene identification and risk evaluation model construction by transcriptome-wide association analysis for salt sensitivity of blood pressure

BACKGROUND

Salt sensitivity of blood pressure (SSBP) is an intermediate phenotype of hypertension and is a predictor of long-term cardiovascular events and death. However, the genetic structures of SSBP are uncertain, and it is difficult to precisely diagnose SSBP in population. So, we aimed to identify genes related to susceptibility to the SSBP, construct a risk evaluation model, and explore the potential functions of these genes.

METHODS AND RESULTS

A genome-wide association study of the systemic epidemiology of salt sensitivity (EpiSS) cohort was performed to obtain summary statistics for SSBP. Then, we conducted a transcriptome-wide association study (TWAS) of 12 tissues using FUSION software to predict the genes associated with SSBP and verified the genes with an mRNA microarray. The potential roles of the genes were explored. Risk evaluation models of SSBP were constructed based on the serial P value thresholds of polygenetic risk scores (PRSs), polygenic transcriptome risk scores (PTRSs) and their combinations of the identified genes and genetic variants from the TWAS. The TWAS revealed that 2605 genes were significantly associated with SSBP. Among these genes, 69 were differentially expressed according to the microarray analysis. The functional analysis showed that the genes identified in the TWAS were enriched in metabolic process pathways. The PRSs were correlated with PTRSs in the heart atrial appendage, adrenal gland, EBV-transformed lymphocytes, pituitary, artery coronary, artery tibial and whole blood. Multiple logistic regression models revealed that a PRS of P < 0.05 had the best predictive ability compared with other PRSs and PTRSs. The combinations of PRSs and PTRSs did not significantly increase the prediction accuracy of SSBP in the training and validation datasets.

CONCLUSIONS

Several known and novel susceptibility genes for SSBP were identified via multitissue TWAS analysis. The risk evaluation model constructed with the PRS of susceptibility genes showed better diagnostic performance than the transcript levels, which could be applied to screen for SSBP high-risk individuals.

Transcriptome-wide association analysis identifies candidate susceptibility genes for prostate-specific antigen levels in men without prostate cancer

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility for prostate cancer (PCa) screening. Using genome-wide association study (GWAS) summary statistics from 95,768 PCa-free men, we conducted a transcriptome-wide association study (TWAS) to examine impacts of genetically predicted gene expression on PSA. Analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10 × 10-6) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61 × 10-6) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses identified 155 statistically significantly (p < 0.05/22,249 = 2.25 × 10-6) genes. Out of 173 unique PSA-associated genes across analyses, we replicated 151 (87.3%) in a TWAS of 209,318 PCa-free individuals from the Million Veteran Program. Based on conditional analyses, we found 20 genes (11 single tissue, nine cross-tissue) that were associated with PSA levels in the discovery TWAS that were not attributable to a lead variant from a GWAS. Ten of these 20 genes replicated, and two of the replicated genes had colocalization probability of >0.5: CCNA2 and HIST1H2BN. Six of the 20 identified genes are not known to impact PCa risk. Fine-mapping based on whole blood and prostate tissue revealed five protein-coding genes with evidence of causal relationships with PSA levels. Of these five genes, four exhibited evidence of colocalization and one was conditionally independent of previous GWAS findings. These results yield hypotheses that should be further explored to improve understanding of genetic factors underlying PSA levels.

A multi-ancestry cerebral cortex transcriptome-wide association study identifies genes associated with smoking behaviors

Transcriptome-wide association studies (TWAS) have provided valuable insight in identifying genes that may impact cigarette smoking. Most of previous studies, however, mainly focused on European ancestry. Limited TWAS studies have been conducted across multiple ancestries to explore genes that may impact smoking behaviors. In this study, we used cis-eQTL data of cerebral cortex from multiple ancestries in MetaBrain, including European, East Asian, and African samples, as reference panels to perform multi-ancestry TWAS analyses on ancestry-matched GWASs of four smoking behaviors including smoking initiation, smoking cessation, age of smoking initiation, and number of cigarettes per day in GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN). Multiple-ancestry fine-mapping approach was conducted to identify credible gene sets associated with these four traits. Enrichment and module network analyses were further performed to explore the potential roles of these identified gene sets. A total of 719 unique genes were identified to be associated with at least one of the four smoking traits across ancestries. Among those, 249 genes were further prioritized as putative causal genes in multiple ancestry-based fine-mapping approach. Several well-known smoking-related genes, including PSMA4, IREB2, and CHRNA3, showed high confidence across ancestries. Some novel genes, e.g., TSPAN3 and ANK2, were also identified in the credible sets. The enrichment analysis identified a series of critical pathways related to smoking such as synaptic transmission and glutamate receptor activity. Leveraging the power of the latest multi-ancestry GWAS and eQTL data sources, this study revealed hundreds of genes and relevant biological processes related to smoking behaviors. These findings provide new insights for future functional studies on smoking behaviors.

Haplotype function score improves biological interpretation and cross-ancestry polygenic prediction of human complex traits

We propose a new framework for human genetic association studies: at each locus, a deep learning model (in this study, Sei) is used to calculate the functional genomic activity score for two haplotypes per individual. This score, defined as the Haplotype Function Score (HFS), replaces the original genotype in association studies. Applying the HFS framework to 14 complex traits in the UK Biobank, we identified 3619 independent HFS-trait associations with a significance of p < 5 × 10-8. Fine-mapping revealed 2699 causal associations, corresponding to a median increase of 63 causal findings per trait compared with single-nucleotide polymorphism (SNP)-based analysis. HFS-based enrichment analysis uncovered 727 pathway-trait associations and 153 tissue-trait associations with strong biological interpretability, including 'circadian pathway-chronotype' and 'arachidonic acid-intelligence'. Lastly, we applied least absolute shrinkage and selection operator (LASSO) regression to integrate HFS prediction score with SNP-based polygenic risk scores, which showed an improvement of 16.1-39.8% in cross-ancestry polygenic prediction. We concluded that HFS is a promising strategy for understanding the genetic basis of human complex traits.

A perspective on translating genomic discoveries into targets for brain-machine interface and deep brain stimulation devices

Mental illnesses have a huge impact on individuals, families, and society, so there is a growing need for more efficient treatments. In this context, brain-computer interface (BCI) technology has the potential to revolutionize the options for neuropsychiatric therapies. However, the development of BCI-based therapies faces enormous challenges, such as power dissipation constraints, lack of credible feedback mechanisms, uncertainty of which brain areas and frequencies to target, and even which patients to treat. Some of these setbacks are due to the large gap in our understanding of brain function. In recent years, large-scale genomic analyses uncovered an unprecedented amount of information regarding the biology of the altered brain function observed across the psychopathology spectrum. We believe findings from genetic studies can be useful to refine BCI technology to develop novel treatment options for mental illnesses. Here, we assess the latest advancements in both fields, the possibilities that can be generated from their intersection, and the challenges that these research areas will need to address to ensure that translational efforts can lead to effective and reliable interventions. Specifically, starting from highlighting the overlap between mechanisms uncovered by large-scale genetic studies and the current targets of deep brain stimulation treatments, we describe the steps that could help to translate genomic discoveries into BCI targets. Because these two research areas have not been previously presented together, the present article can provide a novel perspective for scientists with different research backgrounds. This article is categorized under: Neurological Diseases > Genetics/Genomics/Epigenetics Neurological Diseases > Biomedical Engineering.

Integration of risk factor polygenic risk score with disease polygenic risk score for disease prediction

Polygenic risk score (PRS) is useful for capturing an individual's genetic susceptibility. However, previous studies have not fully exploited the potential of the risk factor PRS (RFPRS) for disease prediction. We explored the potential of integrating disease-related RFPRSs with disease PRS to enhance disease prediction performance. We constructed 112 RFPRSs and analyzed the association of RFPRSs with diseases to identify disease-related RFPRSs in 700 diseases, using the UK Biobank dataset. We uncovered 6157 statistically significant associations between 247 diseases and 109 RFPRSs. We estimated the disease PRSs of 70 diseases that exhibited statistically significant heritability, to generate RFDiseasemetaPRS-a combined PRS integrating RFPRSs and disease PRS-and compare the prediction performance metrics between RFDiseasemetaPRS and disease PRS. RFDiseasemetaPRS showed better performance for Nagelkerke's pseudo-R2, odds ratio (OR) per 1 SD, net reclassification improvement (NRI) values and difference of R2 considered by variance of R2 in 31 out of 70 diseases. Additionally, we assessed risk classification between two models by examining OR between the top 10% and remaining 90% individuals for the 31 diseases; RFDiseasemetaPRS exhibited better R2, NRI and OR than disease PRS. These findings highlight the importance of utilizing RFDiseasemetaPRS, which can provide personalized healthcare and tailored prevention strategies.

The declining insulinogenic index correlates with inflammation and metabolic dysregulation in non-obese individuals assessed by blood gene expression

AIMS

Diabetes onset is difficult to predict. Since decreased insulinogenic index (IGI) is observed in prediabetes, and blood gene expression correlates with insulin secretion, candidate biomarkers can be identified.

METHODS

We collected blood from 96 participants (54 males, 42 females) in 2008 (age: 52.5 years) and 2016 for clinical and gene expression analyses. IGI was derived from values of insulin and glucose at fasting and at 30 min post-OGTT. Two subgroups were identified based on IGI variation: "Minor change in IGI" group with absolute value variation between -0.05 and +0.05, and "Decrease in IGI" group with a variation between -20 and -0.05.

RESULTS

Following the comparison of "Minor change in IGI" and "Decrease in IGI" groups at time 0 (2008), we identified 77 genes correlating with declining IGI, related to response to lipid, carbohydrate, and hormone metabolism, response to stress and DNA metabolic processes. Over the eight years, genes correlating to declining IGI were related to inflammation, metabolic and hormonal dysregulation. Individuals with minor change in IGI, instead, featured homeostatic and regenerative responses.

CONCLUSIONS

By blood gene expression analysis of non-obese individuals, we identified potential gene biomarkers correlating to declining IGI, associated to a pathophysiology of inflammation and metabolic dysregulation.

Whole Blood-based Transcriptional Risk Score for Nonobese Type 2 Diabetes Predicts Dynamic Changes in Glucose Metabolism

CONTEXT

The performance of peripheral blood transcriptional markers in evaluating risk of type 2 diabetes (T2D) with normal body mass index (BMI) is unknown.

OBJECTIVE

We developed a whole blood-based transcriptional risk score (wb-TRS) for nonobese T2D and assessed its contributions on disease risk and dynamic changes in glucose metabolism.

METHODS

Using a community-based cohort with blood transcriptome data, we developed the wb-TRS in 1105 participants aged ≥40 years who maintained a normal BMI for up to 10 years, and we validated the wb-TRS in an external dataset. Potential biological significance was explored.

RESULTS

The wb-TRS included 144 gene transcripts. Compared to the lowest tertile, wb-TRS in tertile 3 was associated with 8.91-fold (95% CI, 3.53-22.5) higher risk and each 1-unit increment was associated with 2.63-fold (95% CI, 1.87-3.68) higher risk of nonobese T2D. Furthermore, baseline wb-TRS significantly associated with dynamic changes in average, daytime, nighttime, and 24-hour glucose, HbA1c values, and area under the curve of glucose measured by continuous glucose monitoring over 6 months of intervention. The wb-TRS improved the prediction performance for nonobese T2D, combined with fasting glucose, triglycerides, and demographic and anthropometric parameters. Multi-contrast gene set enrichment (Mitch) analysis implicated oxidative phosphorylation, mTORC1 signaling, and cholesterol metabolism involved in nonobese T2D pathogenesis.

CONCLUSION

A whole blood-based nonobese T2D-associated transcriptional risk score was validated to predict dynamic changes in glucose metabolism. These findings suggested several biological pathways involved in the pathogenesis of nonobese T2D.

Unraveling the causal genes and transcriptomic determinants of human telomere length

Telomere length (TL) shortening is a pivotal indicator of biological aging and is associated with many human diseases. The genetic determinates of human TL have been widely investigated, however, most existing studies were conducted based on adult tissues which are heavily influenced by lifetime exposure. Based on the analyses of terminal restriction fragment (TRF) length of telomere, individual genotypes, and gene expressions on 166 healthy placental tissues, we systematically interrogate TL-modulated genes and their potential functions. We discover that the TL in the placenta is comparatively longer than in other adult tissues, but exhibiting an intra-tissue homogeneity. Trans-ancestral TL genome-wide association studies (GWASs) on 644,553 individuals identify 20 newly discovered genetic associations and provide increased polygenic determination of human TL. Next, we integrate the powerful TL GWAS with placental expression quantitative trait locus (eQTL) mapping to prioritize 23 likely causal genes, among which 4 are functionally validated, including MMUT, RRM1, KIAA1429, and YWHAZ. Finally, modeling transcriptomic signatures and TRF-based TL improve the prediction performance of human TL. This study deepens our understanding of causal genes and transcriptomic determinants of human TL, promoting the mechanistic research on fine-grained TL regulation.

Phenotype integration improves power and preserves specificity in biobank-based genetic studies of major depressive disorder

Biobanks often contain several phenotypes relevant to diseases such as major depressive disorder (MDD), with partly distinct genetic architectures. Researchers face complex tradeoffs between shallow (large sample size, low specificity/sensitivity) and deep (small sample size, high specificity/sensitivity) phenotypes, and the optimal choices are often unclear. Here we propose to integrate these phenotypes to combine the benefits of each. We use phenotype imputation to integrate information across hundreds of MDD-relevant phenotypes, which significantly increases genome-wide association study (GWAS) power and polygenic risk score (PRS) prediction accuracy of the deepest available MDD phenotype in UK Biobank, LifetimeMDD. We demonstrate that imputation preserves specificity in its genetic architecture using a novel PRS-based pleiotropy metric. We further find that integration via summary statistics also enhances GWAS power and PRS predictions, but can introduce nonspecific genetic effects depending on input. Our work provides a simple and scalable approach to improve genetic studies in large biobanks by integrating shallow and deep phenotypes.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

An integrated gene-to-outcome multimodal database for metabolic dysfunction-associated steatotic liver disease

Metabolic dysfunction-associated steatotic liver disease (MASLD) is the commonest cause of chronic liver disease worldwide and represents an unmet precision medicine challenge. We established a retrospective national cohort of 940 histologically defined patients (55.4% men, 44.6% women; median body mass index 31.3; 32% with type 2 diabetes) covering the complete MASLD severity spectrum, and created a secure, searchable, open resource (SteatoSITE). In 668 cases and 39 controls, we generated hepatic bulk RNA sequencing data and performed differential gene expression and pathway analysis, including exploration of gender-specific differences. A web-based gene browser was also developed. We integrated histopathological assessments, transcriptomic data and 5.67 million days of time-stamped longitudinal electronic health record data to define disease-stage-specific gene expression signatures, pathogenic hepatic cell subpopulations and master regulator networks associated with adverse outcomes in MASLD. We constructed a 15-gene transcriptional risk score to predict future hepatic decompensation events (area under the receiver operating characteristic curve 0.86, 0.81 and 0.83 for 1-, 3- and 5-year risk, respectively). Additionally, thyroid hormone receptor beta regulon activity was identified as a critical suppressor of disease progression. SteatoSITE supports rational biomarker and drug development and facilitates precision medicine approaches for patients with MASLD.

The functional impact of rare variation across the regulatory cascade

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Projecting genetic associations through gene expression patterns highlights disease etiology and drug mechanisms

Genes act in concert with each other in specific contexts to perform their functions. Determining how these genes influence complex traits requires a mechanistic understanding of expression regulation across different conditions. It has been shown that this insight is critical for developing new therapies. Transcriptome-wide association studies have helped uncover the role of individual genes in disease-relevant mechanisms. However, modern models of the architecture of complex traits predict that gene-gene interactions play a crucial role in disease origin and progression. Here we introduce PhenoPLIER, a computational approach that maps gene-trait associations and pharmacological perturbation data into a common latent representation for a joint analysis. This representation is based on modules of genes with similar expression patterns across the same conditions. We observe that diseases are significantly associated with gene modules expressed in relevant cell types, and our approach is accurate in predicting known drug-disease pairs and inferring mechanisms of action. Furthermore, using a CRISPR screen to analyze lipid regulation, we find that functionally important players lack associations but are prioritized in trait-associated modules by PhenoPLIER. By incorporating groups of co-expressed genes, PhenoPLIER can contextualize genetic associations and reveal potential targets missed by single-gene strategies.

Evaluating the Classification Accuracy of Expression Quantitative Trait Loci Calculated Polygenic Risk Scores in Alzheimer's Disease

Polygenic risk scores (PRS) hold promise for the early identification of those at risk from neurodegenerative disorders such as Alzheimer's Disease (AD), allowing for intervention to occur prior to neuronal damage. The current selection of informative single nucleotide polymorphisms (SNPs) to generate the risk scores is based on the modelling of large genome-wide association data using significance thresholds. However, the biological relevance of these SNPs is largely unknown. This study, in contrast, aims to identify SNPs with biological relevance to AD and then assess them for their ability to accurately classify cases and controls. Samples selected from the Brains for Dementia Research (BDR) were used to produce gene expression data to identify potential expression quantitative trait loci (eQTLs) relevant to AD. These SNPs were then incorporated into a PRS model to classify AD and controls in the full BDR cohort. Models derived from these eQTLs demonstrate modest classification potential with an accuracy between 61% and 67%. Although the model accuracy is not as high as some values in the literature based on significance thresholds from genome-wide association studies, these models may reflect a more biologically relevant model, which may provide novel targets for therapeutic intervention.

Biobank-scale methods and projections for sparse polygenic prediction from machine learning

In this paper we characterize the performance of linear models trained via widely-used sparse machine learning algorithms. We build polygenic scores and examine performance as a function of training set size, genetic ancestral background, and training method. We show that predictor performance is most strongly dependent on size of training data, with smaller gains from algorithmic improvements. We find that LASSO generally performs as well as the best methods, judged by a variety of metrics. We also investigate performance characteristics of predictors trained on one genetic ancestry group when applied to another. Using LASSO, we develop a novel method for projecting AUC and correlation as a function of data size (i.e., for new biobanks) and characterize the asymptotic limit of performance. Additionally, for LASSO (compressed sensing) we show that performance metrics and predictor sparsity are in agreement with theoretical predictions from the Donoho-Tanner phase transition. Specifically, a future predictor trained in the Taiwan Precision Medicine Initiative for asthma can achieve an AUC of [Formula: see text] and for height a correlation of [Formula: see text] for a Taiwanese population. This is above the measured values of [Formula: see text] and [Formula: see text], respectively, for UK Biobank trained predictors applied to a European population.

From function to translation: Decoding genetic susceptibility to human diseases via artificial intelligence

While genome-wide association studies (GWAS) have discovered thousands of disease-associated loci, molecular mechanisms for a considerable fraction of the loci remain to be explored. The logical next steps for post-GWAS are interpreting these genetic associations to understand disease etiology (GWAS functional studies) and translating this knowledge into clinical benefits for the patients (GWAS translational studies). Although various datasets and approaches using functional genomics have been developed to facilitate these studies, significant challenges remain due to data heterogeneity, multiplicity, and high dimensionality. To address these challenges, artificial intelligence (AI) technology has demonstrated considerable promise in decoding complex functional datasets and providing novel biological insights into GWAS findings. This perspective first describes the landmark progress driven by AI in interpreting and translating GWAS findings and then outlines specific challenges followed by actionable recommendations related to data availability, model optimization, and interpretation, as well as ethical concerns.

Combinations of genes at the 16p11.2 and 22q11.2 CNVs contribute to neurobehavioral traits

The 16p11.2 and 22q11.2 copy number variants (CNVs) are associated with neurobehavioral traits including autism spectrum disorder (ASD), schizophrenia, bipolar disorder, obesity, and intellectual disability. Identifying specific genes contributing to each disorder and dissecting the architecture of CNV-trait association has been difficult, inspiring hypotheses of more complex models, such as multiple genes acting together. Using multi-tissue data from the GTEx consortium, we generated pairwise expression imputation models for CNV genes and then applied these elastic net models to GWAS for: ASD, bipolar disorder, schizophrenia, BMI (obesity), and IQ (intellectual disability). We compared the variance in these five traits explained by gene pairs with the variance explained by single genes and by traditional interaction models. We also modeled polygene region-wide effects using summed predicted expression ranks across many genes to create a regionwide score. We found that in all CNV-trait pairs except for bipolar disorder at 22q11.2, pairwise effects explain more variance than single genes. Pairwise model superiority was specific to the CNV region for all 16p11.2 traits and ASD at 22q11.2. We identified novel individual genes over-represented in top pairs that did not show single-gene signal. We also found that BMI and IQ have significant regionwide association with both CNV regions. Overall, we observe that genetic architecture differs by trait and region, but 9/10 CNV-trait combinations demonstrate evidence for multigene contribution, and for most of these, the importance of combinatorial models appears unique to CNV regions. Our results suggest that mechanistic insights for CNV pathology may require combinational models.

Transcriptome-Wide Association Analysis Identifies Novel Candidate Susceptibility Genes for Prostate-Specific Antigen Levels in Men Without Prostate Cancer

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility to screen for prostate cancer (PCa). We thus conducted a transcriptome-wide association study (TWAS) of PSA levels using genome-wide summary statistics from 95,768 PCa-free men, the MetaXcan framework, and gene prediction models trained in Genotype-Tissue Expression (GTEx) project data. Tissue-specific analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10e-6) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61e-6) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses that combined associations across 45 tissues identified 155 genes that were statistically significantly (p < 0.05/22,249 = 2.25e-6) associated with PSA levels. Based on conditional analyses that assessed whether TWAS associations were attributable to a lead GWAS variant, we found 20 novel genes (11 single-tissue, 9 cross-tissue) that were associated with PSA levels in the TWAS. Of these novel genes, five showed evidence of colocalization (colocalization probability > 0.5): EXOSC9, CCNA2, HIST1H2BN, RP11-182L21.6, and RP11-327J17.2. Six of the 20 novel genes are not known to impact PCa risk. These findings yield new hypotheses for genetic factors underlying PSA levels that should be further explored toward improving our understanding of PSA biology.

Examining the biological mechanisms of human mental disorders resulting from gene-environment interdependence using novel functional genomic approaches

We explore how functional genomics approaches that integrate datasets from human and non-human model systems can improve our understanding of the effect of gene-environment interplay on the risk for mental disorders. We start by briefly defining the G-E paradigm and its challenges and then discuss the different levels of regulation of gene expression and the corresponding data existing in humans (genome wide genotyping, transcriptomics, DNA methylation, chromatin modifications, chromosome conformational changes, non-coding RNAs, proteomics and metabolomics), discussing novel approaches to the application of these data in the study of the origins of mental health. Finally, we discuss the multilevel integration of diverse types of data. Advance in the use of functional genomics in the context of a G-E perspective improves the detection of vulnerabilities, informing the development of preventive and therapeutic interventions.

Detecting and Adjusting for Hidden Biases due to Phenotype Misclassification in Genome-Wide Association Studies

With the advent of healthcare-based genotyped biobanks, genome-wide association studies (GWAS) leverage larger sample sizes, incorporate patients with diverse ancestries and introduce noisier phenotypic definitions. Yet the extent and impact of phenotypic misclassification on large-scale datasets is not currently well understood due to a lack of statistical methods to estimate relevant parameters from empirical data. Here, we develop a statistical method and scalable software, PheMED, Phenotypic Measurement of Effective Dilution, to quantify phenotypic misclassification across GWAS using only summary statistics. We illustrate how the parameters estimated by PheMED relate to the negative and positive predictive value of the labeled phenotype, compared to ground truth, and how misclassification of the phenotype yields diluted effect-sizes of variant-phenotype associations. Furthermore, we apply our methodology to detect multiple instances of statistically significant dilution in real-world data. We demonstrate how effective dilution biases downstream GWAS replication and heritability analyses despite utilizing current best practices, and provide a dilution-aware meta-analysis approach that outperforms existing methods. Consequently, we anticipate that PheMED will be a valuable tool for researchers to address phenotypic data quality issues both within and across cohorts.

The Genetically Informed Neurobiology of Addiction (GINA) model

Addictions are heritable and unfold dynamically across the lifespan. One prominent neurobiological theory proposes that substance-induced changes in neural circuitry promote the progression of addiction. Genome-wide association studies have begun to characterize the polygenic architecture undergirding addiction liability and revealed that genetic loci associated with risk can be divided into those associated with a general broad-spectrum liability to addiction and those associated with drug-specific addiction risk. In this Perspective, we integrate these genomic findings with our current understanding of the neurobiology of addiction to propose a new Genetically Informed Neurobiology of Addiction (GINA) model.

Best practices for multi-ancestry, meta-analytic transcriptome-wide association studies: Lessons from the Global Biobank Meta-analysis Initiative

The Global Biobank Meta-analysis Initiative (GBMI), through its diversity, provides a valuable opportunity to study population-wide and ancestry-specific genetic associations. However, with multiple ascertainment strategies and multi-ancestry study populations across biobanks, GBMI presents unique challenges in implementing statistical genetics methods. Transcriptome-wide association studies (TWASs) boost detection power for and provide biological context to genetic associations by integrating genetic variant-to-trait associations from genome-wide association studies (GWASs) with predictive models of gene expression. TWASs present unique challenges beyond GWASs, especially in a multi-biobank, meta-analytic setting. Here, we present the GBMI TWAS pipeline, outlining practical considerations for ancestry and tissue specificity, meta-analytic strategies, and open challenges at every step of the framework. We advise conducting ancestry-stratified TWASs using ancestry-specific expression models and meta-analyzing results using inverse-variance weighting, showing the least test statistic inflation. Our work provides a foundation for adding transcriptomic context to biobank-linked GWASs, allowing for ancestry-aware discovery to accelerate genomic medicine.

The role of parental genotype in the intergenerational transmission of externalizing behavior: Evidence for genetic nurturance

The purpose of this study was to examine possible pathways by which genetic risk associated with externalizing is transmitted in families. We used molecular data to disentangle the genetic and environmental pathways contributing to adolescent externalizing behavior in a sample of 1,111 adolescents (50% female; 719 European and 392 African ancestry) and their parents from the Collaborative Study on the Genetics of Alcoholism. We found evidence for genetic nurture such that parental externalizing polygenic scores were associated with adolescent externalizing behavior, over and above the effect of adolescents' own externalizing polygenic scores. Mediation analysis indicated that parental externalizing psychopathology partly explained the effect of parental genotype on children's externalizing behavior. We also found evidence for evocative gene-environment correlation, whereby adolescent externalizing polygenic scores were associated with lower parent-child communication, less parent-child closeness, and lower parental knowledge, controlling for parental genotype. These effects were observed among participants of European ancestry but not African ancestry, likely due to the limited predictive power of polygenic scores across ancestral background. These results demonstrate that in addition to genetic transmission, genes influence offspring behavior through the influence of parental genotypes on their children's environmental experiences, and the role of children's genotypes in shaping parent-child relationships.

Multi-ancestry fine-mapping improves precision to identify causal genes in transcriptome-wide association studies

Transcriptome-wide association studies (TWASs) are a powerful approach to identify genes whose expression is associated with complex disease risk. However, non-causal genes can exhibit association signals due to confounding by linkage disequilibrium (LD) patterns and eQTL pleiotropy at genomic risk regions, which necessitates fine-mapping of TWAS signals. Here, we present MA-FOCUS, a multi-ancestry framework for the improved identification of genes underlying traits of interest. We demonstrate that by leveraging differences in ancestry-specific patterns of LD and eQTL signals, MA-FOCUS consistently outperforms single-ancestry fine-mapping approaches with equivalent total sample sizes across multiple metrics. We perform TWASs for 15 blood traits using genome-wide summary statistics (average n_EA = 511 k, n_AA = 13 k) and lymphoblastoid cell line eQTL data from cohorts of primarily European and African continental ancestries. We recapitulate evidence demonstrating shared genetic architectures for eQTL and blood traits between the two ancestry groups and observe that gene-level effects correlate 20% more strongly across ancestries than SNP-level effects. Lastly, we perform fine-mapping using MA-FOCUS and find evidence that genes at TWAS risk regions are more likely to be shared across ancestries than they are to be ancestry specific. Using multiple lines of evidence to validate our findings, we find that gene sets produced by MA-FOCUS are more enriched in hematopoietic categories than alternative approaches (p = 2.36 × 10^-15). Our work demonstrates that including and appropriately accounting for genetic diversity can drive more profound insights into the genetic architecture of complex traits.

Novel functional genomics approaches bridging neuroscience and psychiatry

The possibility of establishing a metric of individual genetic risk for a particular disease or trait has sparked the interest of the clinical and research communities, with many groups developing and validating genomic profiling methodologies for their potential application in clinical care. Current approaches for calculating genetic risk to specific psychiatric conditions consist of aggregating genome-wide association studies-derived estimates into polygenic risk scores, which broadly represent the number of inherited risk alleles for an individual. While the traditional approach for polygenic risk score calculation aggregates estimates of gene-disease associations, novel alternative approaches have started to consider functional molecular phenotypes that are closer to genetic variation and are less penalized by the multiple testing required in genome-wide association studies. Moving the focus from genotype-disease to genotype-gene regulation frameworks, these novel approaches incorporate prior knowledge regarding biological processes involved in disease and aggregate estimates for the association of genotypes and phenotypes using multi-omics data modalities. In this review, we discuss and list different functional genomics tools that can be used and integrated to inform researchers and clinicians for a better understanding and diagnosis of psychopathology. We suggest that these novel approaches can help generate biologically driven hypotheses for polygenic signals that can ultimately serve the clinical community as potential biomarkers of psychiatric disease susceptibility.

Can genetics guide exercise prescriptions in osteoarthritis?

Osteoarthritis (OA) is the most common form of arthritis and has a multifactorial etiology. Current management for OA focuses on minimizing pain and functional loss, typically involving pharmacological, physical, psychosocial, and mind-body interventions. However, there remain challenges in determining which patients will benefit most from which interventions. Although exercise-based interventions are recommended as first-line treatments and are known to be beneficial for managing both the disease and illness of OA, the optimal exercise “prescription” is unknown, due in part to our limited understanding of the precise mechanisms underlying its action. Here we present our perspective on the potential role of genetics in guiding exercise prescription for persons with OA. We describe key publications in the areas of exercise and OA, genetics and OA, and exercise and genetics, and point to a paucity of knowledge at the intersection of exercise, genetics, and OA. We suggest there is emerging evidence to support the use of genetics and epigenetics to explain the beneficial effects of exercise for OA. We identify missing links in the existing research relating to exercise, genetics, and OA, and highlight epigenetics as a promising mechanism through which environmental exposures such as exercise may impact OA outcomes. We anticipate future studies will improve our understanding of how genetic and epigenetic factors mediate exercise-based interventions to support implementation and ultimately improve OA patient care.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV₁] and its ratio to forced vital capacity [FEV₁/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV₁ and FEV₁/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10^-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.

Integrative transcriptomic, evolutionary, and causal inference framework for region-level analysis: Application to COVID-19

We developed an integrative transcriptomic, evolutionary, and causal inference framework for a deep region-level analysis, which integrates several published approaches and a new summary-statistics-based methodology. To illustrate the framework, we applied it to understanding the host genetics of COVID-19 severity. We identified putative causal genes, including SLC6A20, CXCR6, CCR9, and CCR5 in the locus on 3p21.31, quantifying their effect on mediating expression and on severe COVID-19. We confirmed that individuals who carry the introgressed archaic segment in the locus have a substantially higher risk of developing the severe disease phenotype, estimating its contribution to expression-mediated heritability using a new summary-statistics-based approach we developed here. Through a large-scale phenome-wide scan for the genes in the locus, several potential complications, including inflammatory, immunity, olfactory, and gustatory traits, were identified. Notably, the introgressed segment showed a much higher concentration of expression-mediated causal effect on severity (0.9–11.5 times) than the entire locus, explaining, on average, 15.7% of the causal effect. The region-level framework (implemented in publicly available software, SEGMENT-SCAN) has important implications for the elucidation of molecular mechanisms of disease and the rational design of potentially novel therapeutics.