10 Years of GWAS Discovery: Biology, Function, and Translation

Building machine learning prediction models for well-being using predictors from the exposome and genome in a population cohort

Effective personalized well-being interventions require the ability to predict who will thrive or not, and the understanding of underlying mechanisms. Here, using longitudinal data of a large population cohort (the Netherlands Twin Register, collected 1991-2022), we aim to build machine learning prediction models for adult well-being from the exposome and genome, and identify the most predictive factors (N between 702 and 5874). The specific exposome was captured by parent and self-reports of psychosocial factors from childhood to adulthood, the genome was described by polygenic scores, and the general exposome was captured by linkage of participants' postal codes to objective, registry-based exposures. Not the genome (R 2 = -0.007 [-0.026-0.010]), but the general exposome (R 2 = 0.047 [0.015-0.076]) and especially the specific exposome (R 2 = 0.702 [0.637-0.753]) were predictive of well-being in an independent test set. Adding the genome (P = 0.334) and general exposome (P = 0.695) independently or jointly (P = 0.029) beyond the specific exposome did not improve prediction. Risk/protective factors such as optimism, personality, social support and neighborhood housing characteristics were most predictive. Our findings highlight the importance of longitudinal monitoring and promises of different data modalities for well-being prediction.

Single nucleotide polymorphism profile for quantitative trait nucleotide in populations with small effective size and its impact on mapping and genomic predictions

Increasing SNP density by incorporating sequence information only marginally increases prediction accuracies of breeding values in livestock. To find out why, we used statistical models and simulations to investigate the shape of distribution of estimated SNP effects (a profile) around quantitative trait nucleotides (QTNs) in populations with a small effective population size (Ne). A QTN profile created by averaging SNP effects around each QTN was similar to the shape of expected pairwise linkage disequilibrium (PLD) based on Ne and genetic distance between SNP, with a distinct peak for the QTN. Populations with smaller Ne showed lower but wider QTN profiles. However, adding more genotyped individuals with phenotypes dragged the profile closer to the QTN. The QTN profile was higher and narrower for populations with larger compared to smaller Ne. Assuming the PLD curve for the QTN profile, 80% of the additive genetic variance explained by each QTN was contained in ± 1/Ne Morgan interval around the QTN, corresponding to 2 Mb in cattle and 5 Mb in pigs and chickens. With such large intervals, identifying QTN is difficult even if all of them are in the data and the assumed genetic architecture is simplistic. Additional complexity in QTN detection arises from confounding of QTN profiles with signals due to relationships, overlapping profiles with closely spaced QTN, and spurious signals. However, small Ne allows for accurate predictions with large data even without QTN identification because QTNs are accounted for by QTN profiles if SNP density is sufficient to saturate the segments.

Genome-wide fine-mapping improves identification of causal variants

Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 600 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.

Insights into preeclampsia: a bioinformatics approach to deciphering genetic and immune contributions

Background

Preeclampsia (PE) is a global pregnancy concern, characterized by hypertension with an unclear etiology. This study employs Mendelian randomization (MR) and single-cell RNA sequencing (scRNA-seq) to clarify its genetic and molecular roots, offering insights into diagnosis and treatment avenues.

Methods

We integrated PE-specific genome-wide association study (GWAS) data, expression and protein quantitative trait loci (eQTL and pQTL) data, and single-cell data from peripheral blood mononuclear cells (PBMCs). We identified highly variable genes using single-cell information and employed MR to determine potential causality. We also combined pQTL and GWAS data, discerned genes positively associated with PE through scRNA-seq, and leveraged the Enrichr platform to unearth drug-gene interactions.

Results

Our scRNA-seq pinpointed notable cell type distribution variances, especially in T helper cells (Th cells), between PE and control groups. We unveiled 591 highly variable genes and 6 directly PE-associated genes. Although MR revealed correlations with PE risk, pQTL analysis was inconclusive due to data constraints. Using DSigDB, 93 potential therapeutic agents, like Retinoic acid targeting core genes (IFITM3, NINJ1, COTL1, CD69, and YWHAZ), emerged as prospective multi-target treatments.

Conclusion

Utilizing MR and scRNA-seq, this study underscores significant cellular disparities, particularly in Th cells, and identifies crucial genes related to PE. Despite some limitations, these genes have been revealed in PE's underlying mechanism. Potential therapeutic agents, such as Retinoic acid, suggest promising treatment pathways.

All of Us diversity and scale improve polygenic prediction contextually with greatest improvements for under-represented populations

Recent studies have demonstrated that polygenic risk scores (PRS) trained on multi-ancestry data can improve prediction accuracy in groups historically underrepresented in genomic studies, but the availability of linked health and genetic data from large-scale diverse cohorts representative of a wide spectrum of human diversity remains limited. To address this need, the All of Us research program (AoU) generated whole-genome sequences of 245,388 individuals who collectively reflect the diversity of the USA. Leveraging this resource and another widely-used population-scale biobank, the UK Biobank (UKB) with a half million participants, we developed PRS trained on multi-ancestry and multi-biobank data with up to ~750,000 participants for 32 common, complex traits and diseases across a range of genetic architectures. We then compared effects of ancestry, PRS methodology, and genetic architecture on PRS accuracy across a held out subset of ancestrally diverse AoU participants. Due to the more heterogeneous study design of AoU, we found lower heritability on average compared to UKB (0.075 vs 0.165), which limited the maximal achievable PRS accuracy in AoU. Overall, we found that the increased diversity of AoU significantly improved PRS performance in some participants in AoU, especially underrepresented individuals, across multiple phenotypes. Notably, maximizing sample size by combining discovery data across AoU and UKB is not the optimal approach for predicting some phenotypes in African ancestry populations; rather, using data from only AoU for these traits resulted in the greatest accuracy. This was especially true for less polygenic traits with large ancestry-enriched effects, such as neutrophil count (R 2: 0.055 vs. 0.035 using AoU vs. cross-biobank meta-analysis, respectively, because of e.g. DARC). Lastly, we calculated individual-level PRS accuracies rather than grouping by continental ancestry, a critical step towards interpretability in precision medicine. Individualized PRS accuracy decays linearly as a function of ancestry divergence, but the slope was smaller using multi-ancestry GWAS compared to using European GWAS. Our results highlight the potential of biobanks with more balanced representations of human diversity to facilitate more accurate PRS for the individuals least represented in genomic studies.

Exposome Profiling of Environmental Pollutants in Seminal Plasma and Novel Associations with Semen Parameters

Indicators of male fertility are in decline globally, but the underlying causes, including the role of environmental exposures, are unclear. This study aimed to examine organic chemical pollutants in seminal plasma, including both known priority environmental chemicals and less studied chemicals, to identify uncharacterized male reproductive environmental toxicants. Semen samples were collected from 100 individuals and assessed for sperm concentration, percent motility, and total motile sperm. Targeted and nontargeted organic pollutant exposures were measured from seminal plasma using gas chromatography, which showed widespread detection of organic pollutants in seminal plasma across all exposure classes. We used principal component pursuit (PCP) on our targeted panel and derived one component (driven by etriadizole) associated with total motile sperm (p < 0.001) and concentration (p = 0.03). This was confirmed by the exposome-wide association models using individual chemicals, where etriadizole was negatively associated with total motile sperm (FDR q = 0.01) and concentration (q = 0.07). Using PCP on 814 nontargeted spectral peaks identified a component that was associated with total motile sperm (p = 0.001). Bayesian kernel machine regression identified one principal driver of this association, which was analytically confirmed to be N-nitrosodiethylamine. These findings are promising and consistent with experimental evidence showing that etridiazole and N-nitrosodiethylamine may be reproductive toxicants.

Genome-wide fine-mapping improves identification of causal variants

Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 599 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.

Identification of Potential Drug Targets for Myopia Through Mendelian Randomization

Purpose

The purpose of this study was to identify potential drug targets for myopia and explore underlying mechanisms.

Methods

Mendelian randomization (MR) was implemented to assess the effect of 2684 pharmacologically targetable genes in the blood and retina on the risk of myopia from a genomewide association study (GWAS) for age-at-onset of spectacle wearing-inferred mean spherical equivalent (MSE; discovery cohort, N = 287,448, European), which was further validated in a GWAS for autorefraction-measured MSE (replication cohort, N = 95,619, European). The reliability of the identified significant potential targets was strengthened by colocalization analysis. Additionally, enrichment analysis, protein-protein interaction network, and molecular docking were performed to explore the functional roles and the druggability of these targets.

Results

This systematic drug target identification has unveiled 6 putative genetically causal targets for myopia-CD34, CD55, Wnt3, LCAT, BTN3A1, and TSSK6-each backed by colocalization evidence in adult blood eQTL datasets. Functional analysis found that dopaminergic neuron differentiation, cell adhesion, Wnt signaling pathway, and plasma lipoprotein-associated pathways may be involved in myopia pathogenesis. Finally, drug prediction and molecular docking corroborated the pharmacological value of these targets with LCAT demonstrating the strongest binding affinity.

Conclusions

Our study not only opens new avenues for the development of therapeutic interventions for myopia but may also help to understand the underlying mechanisms of myopia.

Cellular and molecular mechanisms of skin wound healing

Wound healing is a complex process that involves the coordinated actions of many different tissues and cell lineages. It requires tight orchestration of cell migration, proliferation, matrix deposition and remodelling, alongside inflammation and angiogenesis. Whereas small skin wounds heal in days, larger injuries resulting from trauma, acute illness or major surgery can take several weeks to heal, generally leaving behind a fibrotic scar that can impact tissue function. Development of therapeutics to prevent scarring and successfully repair chronic wounds requires a fuller knowledge of the cellular and molecular mechanisms driving wound healing. In this Review, we discuss the current understanding of the different phases of wound healing, from clot formation through re-epithelialization, angiogenesis and subsequent scar deposition. We highlight the contribution of different cell types to skin repair, with emphasis on how both innate and adaptive immune cells in the wound inflammatory response influence classically studied wound cell lineages, including keratinocytes, fibroblasts and endothelial cells, but also some of the less-studied cell lineages such as adipocytes, melanocytes and cutaneous nerves. Finally, we discuss newer approaches and research directions that have the potential to further our understanding of the mechanisms underpinning tissue repair.

Schoolhouse risk: Can we mitigate the polygenic Pygmalion effect?

BACKGROUND

Although limited in predictive accuracy, polygenic scores (PGS) for educational outcomes are currently available to the public via direct-to-consumer genetic testing companies. Further, there is a growing movement to apply PGS in educational settings via 'precision education.' Prior scholarship highlights the potentially negative impacts of such applications, as disappointing results may give rise a "polygenic Pygmalion effect." In this paper two studies were conducted to identify factors that may mitigate or exacerbate negative impacts of PGS.

METHODS

Two studies were conducted. In each, 1188 students were randomized to one of four conditions: Low-percentile polygenic score for educational attainment (EA-PGS), Low EA-PGS + Mitigating information, Low EA-PGS + Exacerbating information, or Control. Regression analyses were used to examine differences between conditions.

RESULTS

In Study 1, participants randomized to Control reported significantly higher on the Rosenberg Self-Esteem Scale (RSES), Competence Scale (CS), Academic Efficacy Scale (AES) and Educational Potential Scale (EPS). CS was significantly higher in the Low EA-PGS + Mitigating information condition. CS and AES were significantly lower in the Low EA-PGS + Exacerbating information condition compared to the Low EA-PGS + Mitigating information condition. In Study 2, participants randomized to Control reported significantly higher CS and AES. Pairwise comparisons did not show significant differences in CS and AES. Follow-up pairwise comparisons using Tukey P-value correction did not find significant associations between non-control conditions.

CONCLUSION

These studies replicated the polygenic Pygmalion effect yet were insufficiently powered to detect significant effects of mitigating contextual information. Regardless of contextual information, disappointing EA-PGS results were significantly associated with lower assessments of self-esteem, competence, academic efficacy, and educational potential.

Sources of gene expression variation in a globally diverse human cohort

Genetic variation that influences gene expression and splicing is a key source of phenotypic diversity1-5. Although invaluable, studies investigating these links in humans have been strongly biased towards participants of European ancestries, which constrains generalizability and hinders evolutionary research. Here to address these limitations, we developed MAGE, an open-access RNA sequencing dataset of lymphoblastoid cell lines from 731 individuals from the 1000 Genomes Project6, spread across 5 continental groups and 26 populations. Most variation in gene expression (92%) and splicing (95%) was distributed within versus between populations, which mirrored the variation in DNA sequence. We mapped associations between genetic variants and expression and splicing of nearby genes (cis-expression quantitative trait loci (eQTLs) and cis-splicing QTLs (sQTLs), respectively). We identified more than 15,000 putatively causal eQTLs and more than 16,000 putatively causal sQTLs that are enriched for relevant epigenomic signatures. These include 1,310 eQTLs and 1,657 sQTLs that are largely private to underrepresented populations. Our data further indicate that the magnitude and direction of causal eQTL effects are highly consistent across populations. Moreover, the apparent 'population-specific' effects observed in previous studies were largely driven by low resolution or additional independent eQTLs of the same genes that were not detected. Together, our study expands our understanding of human gene expression diversity and provides an inclusive resource for studying the evolution and function of human genomes.

Finemap-MiXeR: A variational Bayesian approach for genetic finemapping

Genome-wide association studies (GWAS) implicate broad genomic loci containing clusters of highly correlated genetic variants. Finemapping techniques can select and prioritize variants within each GWAS locus which are more likely to have a functional influence on the trait. Here, we present a novel method, Finemap-MiXeR, for finemapping causal variants from GWAS summary statistics, controlling for correlation among variants due to linkage disequilibrium. Our method is based on a variational Bayesian approach and direct optimization of the Evidence Lower Bound (ELBO) of the likelihood function derived from the MiXeR model. After obtaining the analytical expression for ELBO's gradient, we apply Adaptive Moment Estimation (ADAM) algorithm for optimization, allowing us to obtain the posterior causal probability of each variant. Using these posterior causal probabilities, we validated Finemap-MiXeR across a wide range of scenarios using both synthetic data, and real data on height from the UK Biobank. Comparison of Finemap-MiXeR with two existing methods, FINEMAP and SuSiE RSS, demonstrated similar or improved accuracy. Furthermore, our method is computationally efficient in several aspects. For example, unlike many other methods in the literature, its computational complexity does not increase with the number of true causal variants in a locus and it does not require any matrix inversion operation. The mathematical framework of Finemap-MiXeR is flexible and may also be applied to other problems including cross-trait and cross-ancestry finemapping.

Genome-wide association study for morphological and hunting-behavior traits in Braque Français Type Pyrénées dogs: A preliminary study

High-throughput genotyping offers great potential to increase our understanding of the genomic basis of canid variation. Braque Français Type Pyrénées (BRA) are smart, agile, and friendly dogs originally developed for tracking, hunting, and retrieving feathered game. On a population of 44 unrelated BRA dogs, single nucleotide polymorphism (SNP) genotype data from the CanineHD Whole-Genome Genotyping BeadChip and evaluation scores for 12 traits related to morphology and hunting performance were available. After quality filtering, 95,859 SNPs on the 38 dog autosomes (CFA) were retained. Phenotypic scores were expressed on a scale from 1 (worst) to 6 (best) and were mostly poorly to moderately correlated except for some morphological traits (e.g. r = 0.81 between the conformation of the head and that of the eye). From GWAS, a total of 378 SNP-phenotype associations with posterior odds of association > 1 have been detected. The strongest associations were found for the eye conformation, for the skull/muzzle ratio, and for connection to the hunter. These included both new and previously identified markers and genes potentially involved with type and behavior traits in BRA. Six of the significant markers mapped within SETDB2, a gene known to be related to pointing behavior in dogs. These results advance our understanding of the genetic basis for morphology and hunting behavior in dogs and identify new variants which are potential targets for further research.

Common genetic variants associated with urinary phthalate levels in children: A genome-wide study

INTRODUCTION

Phthalates, or dieters of phthalic acid, are a ubiquitous type of plasticizer used in a variety of common consumer and industrial products. They act as endocrine disruptors and are associated with increased risk for several diseases. Once in the body, phthalates are metabolized through partially known mechanisms, involving phase I and phase II enzymes.

OBJECTIVE

In this study we aimed to identify common single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) associated with the metabolism of phthalate compounds in children through genome-wide association studies (GWAS).

METHODS

The study used data from 1,044 children with European ancestry from the Human Early Life Exposome (HELIX) cohort. Ten phthalate metabolites were assessed in a two-void pooled urine collected at the mean age of 8 years. Six ratios between secondary and primary phthalate metabolites were calculated. Genome-wide genotyping was done with the Infinium Global Screening Array (GSA) and imputation with the Haplotype Reference Consortium (HRC) panel. PennCNV was used to estimate copy number variants (CNVs) and CNVRanger to identify consensus regions. GWAS of SNPs and CNVs were conducted using PLINK and SNPassoc, respectively. Subsequently, functional annotation of suggestive SNPs (p-value < 1E-05) was done with the FUMA web-tool.

RESULTS

We identified four genome-wide significant (p-value < 5E-08) loci at chromosome (chr) 3 (FECHP1 for oxo-MiNP_oh-MiNP ratio), chr6 (SLC17A1 for MECPP_MEHHP ratio), chr9 (RAPGEF1 for MBzP), and chr10 (CYP2C9 for MECPP_MEHHP ratio). Moreover, 115 additional loci were found at suggestive significance (p-value < 1E-05). Two CNVs located at chr11 (MRGPRX1 for oh-MiNP and SLC35F2 for MEP) were also identified. Functional annotation pointed to genes involved in phase I and phase II detoxification, molecular transfer across membranes, and renal excretion.

CONCLUSION

Through genome-wide screenings we identified known and novel loci implicated in phthalate metabolism in children. Genes annotated to these loci participate in detoxification, transmembrane transfer, and renal excretion.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Polygenic Indices (a.k.a. Polygenic Scores) in Social Science: A Guide for Interpretation and Evaluation

Polygenic indices (PGI)-the new recommended label for polygenic scores (PGS) in social science-are genetic summary scales often used to represent an individual's liability for a disease, trait, or behavior based on the additive effects of measured genetic variants. Enthusiasm for linking genetic data with social outcomes and the inclusion of premade PGIs in social science datasets have facilitated increased uptake of PGIs in social science research-a trend that will likely continue. Yet, most social scientists lack the expertise to interpret and evaluate PGIs in social science research. Here, we provide a primer on PGIs for social scientists focusing on key concepts, unique statistical genetic considerations, and best practices in calculation, estimation, reporting, and interpretation. We summarize our recommended best practices as a checklist to aid social scientists in evaluating and interpreting studies with PGIs. We conclude by discussing the similarities between PGIs and standard social science scales and unique interpretative considerations.

Genome-Wide Association Study as an Efficacious Approach to Discover Candidate Genes Associated with Body Linear Type Traits in Dairy Cattle

Body shape traits are very important and play a crucial role in the economic development of dairy farming. By improving the accuracy of selection for body size traits, we can enhance economic returns across the dairy industry and on farms, contributing to the future profitability of the dairy sector. Registered body conformation traits are reliable and cost-effective tools for use in national cattle breeding selection programs. These traits are significantly related to the production, longevity, mobility, health, fertility, and environmental adaptation of dairy cows. Therefore, they can be considered indirect indicators of economically important traits in dairy cows. Utilizing efficacious genetic methods, such as genome-wide association studies (GWASs), allows for a deeper understanding of the genetic architecture of complex traits through the identification and application of genetic markers. In the current review, we summarize information on candidate genes and genomic regions associated with body conformation traits in dairy cattle worldwide. The manuscript also reviews the importance of body conformation, the relationship between body conformation traits and other traits, heritability, influencing factors, and the genetics of body conformation traits. The information on candidate genes related to body conformation traits provided in this review may be helpful in selecting potential genetic markers for the genetic improvement of body conformation traits in dairy cattle.

A high-dimensional omnibus test for set-based association analysis

Set-based association analysis is a valuable tool in studying the etiology of complex diseases in genome-wide association studies, as it allows for the joint testing of variants in a region or group. Two common types of single nucleotide polymorphism (SNP)-disease functional models are recognized when evaluating the joint function of a set of SNP: the cumulative weak signal model, in which multiple functional variants with small effects contribute to disease risk, and the dominating strong signal model, in which a few functional variants with large effects contribute to disease risk. However, existing methods have two main limitations that reduce their power. Firstly, they typically only consider one disease-SNP association model, which can result in significant power loss if the model is misspecified. Secondly, they do not account for the high-dimensional nature of SNPs, leading to low power or high false positives. In this study, we propose a solution to these challenges by using a high-dimensional inference procedure that involves simultaneously fitting many SNPs in a regression model. We also propose an omnibus testing procedure that employs a robust and powerful P-value combination method to enhance the power of SNP-set association. Our results from extensive simulation studies and a real data analysis demonstrate that our set-based high-dimensional inference strategy is both flexible and computationally efficient and can substantially improve the power of SNP-set association analysis. Application to a real dataset further demonstrates the utility of the testing strategy.

Unsupervised clustering identified clinically relevant metabolic syndrome endotypes in UK and Taiwan Biobanks

Metabolic syndrome (MetS) is a collection of cardiovascular risk factors; however, the high prevalence and heterogeneity impede effective clinical management. We conducted unsupervised clustering on individuals from UK Biobank to reveal endotypes. Five MetS subgroups were identified: Cluster 1 (C1): non-descriptive, Cluster 2 (C2): hypertensive, Cluster 3 (C3): obese, Cluster 4 (C4): lipodystrophy-like, and Cluster 5 (C5): hyperglycemic. For all of the endotypes, we identified the corresponding cardiometabolic traits and their associations with clinical outcomes. Genome-wide association studies (GWASs) were conducted to identify associated genotypic traits. We then determined endotype-specific genotypic traits and constructed polygenic risk score (PRS) models specific to each endotype. GWAS of each MetS clusters revealed different genotypic traits. C1 GWAS revealed novel findings of TRIM63, MYBPC3, MYLPF, and RAPSN. Intriguingly, C1, C3, and C4 were associated with genes highly expressed in brain tissues. MetS clusters with comparable phenotypic and genotypic traits were identified in Taiwan Biobank.

Trait selection strategy in multi-trait GWAS: Boosting SNP discoverability

Since the first genome-wide association studies (GWASs), thousands of variant-trait associations have been discovered. However, comprehensively mapping the genetic determinant of complex traits through univariate testing can require prohibitive sample sizes. Multi-trait GWAS can circumvent this issue and improve statistical power by leveraging the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been solved, the strategy to select traits has been overlooked. In this study, we conducted multi-trait GWAS on approximately 20,000 combinations of 72 traits using an omnibus test as implemented in the Joint Analysis of Summary Statistics. We assessed which genetic features of the sets of traits analyzed were associated with an increased detection of variants compared with univariate screening. Several features of the set of traits, including the heritability, the number of traits, and the genetic correlation, drive the multi-trait test gain. Using these features jointly in predictive models captures a large fraction of the power gain of the multi-trait test (Pearson's r between the observed and predicted gain equals 0.43, p < 1.6 × 10-60). Applying an alternative multi-trait approach (Multi-Trait Analysis of GWAS), we identified similar features of interest, but with an overall 70% lower number of new associations. Finally, selecting sets based on our data-driven models systematically outperformed the common strategy of selecting clinically similar traits. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outlines practical strategies for multi-trait testing.

An externally validated resting-state brain connectivity signature of pain-related learning

Pain can be conceptualized as a precision signal for reinforcement learning in the brain and alterations in these processes are a hallmark of chronic pain conditions. Investigating individual differences in pain-related learning therefore holds important clinical and translational relevance. Here, we developed and externally validated a novel resting-state brain connectivity-based predictive model of pain-related learning. The pre-registered external validation indicates that the proposed model explains 8-12% of the inter-individual variance in pain-related learning. Model predictions are driven by connections of the amygdala, posterior insula, sensorimotor, frontoparietal, and cerebellar regions, outlining a network commonly described in aversive learning and pain. We propose the resulting model as a robust and highly accessible biomarker candidate for clinical and translational pain research, with promising implications for personalized treatment approaches and with a high potential to advance our understanding of the neural mechanisms of pain-related learning.

Gut microbiota, plasma metabolites, and osteoporosis: unraveling links via Mendelian randomization

Objective

Osteoporosis, characterized by reduced bone density and heightened fracture risk, is influenced by genetic and environmental factors. This study investigates the interplay between gut microbiota, plasma metabolomics, and osteoporosis, identifying potential causal relationships mediated by plasma metabolites.

Methods

Utilizing aggregated genome-wide association studies (GWAS) data, a comprehensive two-sample Mendelian Randomization (MR) analysis was performed involving 196 gut microbiota taxa, 1,400 plasma metabolites, and osteoporosis indicators. Causal relationships between gut microbiota, plasma metabolites, and osteoporosis were explored.

Results

The MR analyses revealed ten gut microbiota taxa associated with osteoporosis, with five taxa positively linked to increased risk and five negatively associated. Additionally, 96 plasma metabolites exhibited potential causal relationships with osteoporosis, with 49 showing positive associations and 47 displaying negative associations. Mediation analyses identified six causal pathways connecting gut microbiota to osteoporosis through ten mediating relationships involving seven distinct plasma metabolites, two of which demonstrated suppression effects.

Conclusion

This study provides suggestive evidence of genetic correlations and causal links between gut microbiota, plasma metabolites, and osteoporosis. The findings underscore the complex, multifactorial nature of osteoporosis and suggest the potential of gut microbiota and plasma metabolite profiles as biomarkers or therapeutic targets in the management of osteoporosis.

Post-GWAS Validation of Target Genes Associated with HbF and HbA2 Levels

Genome-Wide Association Studies (GWASs) have identified a huge number of variants associated with different traits. However, their validation through in vitro and in vivo studies often lags well behind their identification. For variants associated with traits or diseases of biomedical interest, this gap delays the development of possible therapies. This issue also impacts beta-hemoglobinopathies, such as beta-thalassemia and sickle cell disease (SCD). The definitive cures for these diseases are currently bone marrow transplantation and gene therapy. However, limitations regarding their effective use restrict their worldwide application. Great efforts have been made to identify whether modulators of fetal hemoglobin (HbF) and, to a lesser extent, hemoglobin A2 (HbA2) are possible therapeutic targets. Herein, we performed the post-GWAS in vivo validation of two genes, cyclin D3 (CCND3) and nuclear factor I X (NFIX), previously associated with HbF and HbA2 levels. The absence of Ccnd3 expression in vivo significantly increased g (HbF) and d (HbA2) globin gene expression. Our data suggest that CCND3 is a possible therapeutic target in sickle cell disease. We also confirmed the association of Nfix with γ-globin gene expression and present data suggesting a possible role for Nfix in regulating Kruppel-like transcription factor 1 (Klf1), a master regulator of hemoglobin switching. This study contributes to filling the gap between GWAS variant identification and target validation for beta-hemoglobinopathies.

Dare to be resilient: the key to future pesticide-free orchards?

Considering the urgent need for more sustainable fruit tree production, it is high time to find durable alternatives to the systematic use of phytosanitary products in orchards. To this end, resilience can deliver a number of benefits. Relying on a combination of tolerance, resistance, and recovery traits, disease resilience appears as a cornerstone to cope with the multiple pest and disease challenges over an orchard's lifetime. Here, we describe resilience as the capacity of a tree to be minimally affected by external disturbances or to rapidly bounce back to normal functioning after being exposed to these disturbances. Based on a literature survey largely inspired from research on livestock, we highlight different approaches for dissecting phenotypic and genotypic components of resilience. In particular, multisite experimental designs and longitudinal measures of so-called 'resilience biomarkers' are required. We identified a list of promising biomarkers relying on ecophysiological and digital measurements. Recent advances in high-throughput phenotyping and genomics tools will likely facilitate fine scale temporal monitoring of tree health, allowing identification of resilient genotypes with the calculation of specific resilience indicators. Although resilience could be considered as a 'black box' trait, we demonstrate how it could become a realistic breeding goal.

Whole exome sequencing analysis identifies genes for alcohol consumption

Alcohol consumption is a heritable behavior seriously endangers human health. However, genetic studies on alcohol consumption primarily focuses on common variants, while insights from rare coding variants are lacking. Here we leverage whole exome sequencing data across 304,119 white British individuals from UK Biobank to identify protein-coding variants associated with alcohol consumption. Twenty-five variants are associated with alcohol consumption through single variant analysis and thirteen genes through gene-based analysis, ten of which have not been reported previously. Notably, the two unreported alcohol consumption-related genes GIGYF1 and ANKRD12 show enrichment in brain function-related pathways including glial cell differentiation and are strongly expressed in the cerebellum. Phenome-wide association analyses reveal that alcohol consumption-related genes are associated with brain white matter integrity and risk of digestive and neuropsychiatric diseases. In summary, this study enhances the comprehension of the genetic architecture of alcohol consumption and implies biological mechanisms underlying alcohol-related adverse outcomes.

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

How group structure impacts the numbers at risk for coronary artery disease: polygenic risk scores and nongenetic risk factors in the UK Biobank cohort

The UK Biobank (UKB) is a large cohort study that recruited over 500,000 British participants aged 40-69 in 2006-2010 at 22 assessment centers from across the United Kingdom. Self-reported health outcomes and hospital admission data are 2 types of records that include participants' disease status. Coronary artery disease (CAD) is the most common cause of death in the UKB cohort. After distinguishing between prevalence and incidence CAD events for all UKB participants, we identified geographical variations in age-standardized rates of CAD between assessment centers. Significant distributional differences were found between the pooled cohort equation scores of UKB participants from England and Scotland using the Mann-Whitney test. Polygenic risk scores of UKB participants from England and Scotland and from different assessment centers differed significantly using permutation tests. Our aim was to discriminate between assessment centers with different disease rates by collecting data on disease-related risk factors. However, relying solely on individual-level predictions and averaging them to obtain group-level predictions proved ineffective, particularly due to the presence of correlated covariates resulting from participation bias. By using the Mundlak model, which estimates a random effects regression by including the group means of the independent variables in the model, we effectively addressed these issues. In addition, we designed a simulation experiment to demonstrate the functionality of the Mundlak model. Our findings have applications in public health funding and strategy, as our approach can be used to predict case rates in the future, as both population structure and lifestyle changes are uncertain.

Functional Genomics in Psoriasis

Psoriasis is an autoimmune cutaneous condition that significantly impacts quality of life and represents a burden on society due to its prevalence. Genome-wide association studies (GWASs) have pinpointed several psoriasis-related risk loci, underlining the disease's complexity. Functional genomics is paramount to unveiling the role of such loci in psoriasis and disentangling its complex nature. In this review, we aim to elucidate the main findings in this field and integrate our discussion with gold-standard techniques in molecular biology-i.e., Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-and high-throughput technologies. These tools are vital to understanding how disease risk loci affect gene expression in psoriasis, which is crucial in identifying new targets for personalized treatments in advanced precision medicine.

PathGPS: discover shared genetic architecture using GWAS summary data

The increasing availability and scale of biobanks and "omic" datasets bring new horizons for understanding biological mechanisms. PathGPS is an exploratory data analysis tool to discover genetic architectures using Genome Wide Association Studies (GWAS) summary data. PathGPS is based on a linear structural equation model where traits are regulated by both genetic and environmental pathways. PathGPS decouples the genetic and environmental components by contrasting the GWAS associations of "signal" genes with those of "noise" genes. From the estimated genetic component, PathGPS then extracts genetic pathways via principal component and factor analysis, leveraging the low-rank and sparse properties. In addition, we provide a bootstrap aggregating ("bagging") algorithm to improve stability under data perturbation and hyperparameter tuning. When applied to a metabolomics dataset and the UK Biobank, PathGPS confirms several known gene-trait clusters and suggests multiple new hypotheses for future investigations.

Distinct genetic liability profiles define clinically relevant patient strata across common diseases

Stratified medicine holds great promise to tailor treatment to the needs of individual patients. While genetics holds great potential to aid patient stratification, it remains a major challenge to operationalize complex genetic risk factor profiles to deconstruct clinical heterogeneity. Contemporary approaches to this problem rely on polygenic risk scores (PRS), which provide only limited clinical utility and lack a clear biological foundation. To overcome these limitations, we develop the CASTom-iGEx approach to stratify individuals based on the aggregated impact of their genetic risk factor profiles on tissue specific gene expression levels. The paradigmatic application of this approach to coronary artery disease or schizophrenia patient cohorts identified diverse strata or biotypes. These biotypes are characterized by distinct endophenotype profiles as well as clinical parameters and are fundamentally distinct from PRS based groupings. In stark contrast to the latter, the CASTom-iGEx strategy discovers biologically meaningful and clinically actionable patient subgroups, where complex genetic liabilities are not randomly distributed across individuals but rather converge onto distinct disease relevant biological processes. These results support the notion of different patient biotypes characterized by partially distinct pathomechanisms. Thus, the universally applicable approach presented here has the potential to constitute an important component of future personalized medicine paradigms.

Towards a Novel Frontier in the Use of Epigenetic Clocks in Epidemiology

Health problems associated with aging are a major public health concern for the future. Aging is a complex process with wide intervariability among individuals. Therefore, there is a need for innovative public health strategies that target factors associated with aging and the development of tools to assess the effectiveness of these strategies accurately. Novel approaches to measure biological age, such as epigenetic clocks, have become relevant. These clocks use non-sequential variable information from the genome and employ mathematical algorithms to estimate biological age based on DNA methylation levels. Therefore, in the present study, we comprehensively review the current status of the epigenetic clocks and their associations across the human phenome. We emphasize the potential utility of these tools in an epidemiological context, particularly in evaluating the impact of public health interventions focused on promoting healthy aging. Our review describes associations between epigenetic clocks and multiple traits across the life and health span. Additionally, we highlighted the evolution of studies beyond mere associations to establish causal mechanisms between epigenetic age and disease. We explored the application of epigenetic clocks to measure the efficacy of interventions focusing on rejuvenation.

Social determinants of health, the microbiome, and surgical injury

ABSTRACT

Postinjury infection continues to plague trauma and emergency surgery patients fortunate enough to survive the initial injury. Rapid response systems, massive transfusion protocols, and the development of level 1 trauma centers, among others, have improved the outcome for millions of patients worldwide. Nonetheless, despite this excellent initial care, patients still remain vulnerable to postinjury infections that can result in organ failure, prolonged critical illness, and even death. While risk factors have been identified (degree of injury, blood loss, time to definitive care, immunocompromise, etc.), they remain probabilistic, not deterministic, and do not explain outcome variability at the individual case level. Here, we assert that analysis of the social determinants of health, as reflected in the patient's microbiome composition (i.e., community structure, membership) and function (metabolomic output), may offer a "window" with which to define individual variability following traumatic injury. Given emerging knowledge in the field, a more comprehensive evaluation of biomarkers within the patient's microbiome, from stool-based microbial metabolites to those in plasma and those present in exhaled breath, when coupled with clinical metadata and machine learning, could lead to a more deterministic assessment of an individual's risk for a poor outcome and those factors that are modifiable. The aim of this piece is to examine how measurable elements of the social determinants of health and the life history of the patient may be buried within the ecologic memory of the gut microbiome. Here we posit that interrogation of the gut microbiome in this manner may be used to inform novel approaches to drive recovery following a surgical injury.

Dual-trait genomic analysis in highly stratified Arabidopsis thaliana populations using genome-wide association summary statistics

Genome-wide association study (GWAS) is a powerful tool to identify genomic loci underlying complex traits. However, the application in natural populations comes with challenges, especially power loss due to population stratification. Here, we introduce a bivariate analysis approach to a GWAS dataset of Arabidopsis thaliana. We demonstrate the efficiency of dual-phenotype analysis to uncover hidden genetic loci masked by population structure via a series of simulations. In real data analysis, a common allele, strongly confounded with population structure, is discovered to be associated with late flowering and slow maturation of the plant. The discovered genetic effect on flowering time is further replicated in independent datasets. Using Mendelian randomization analysis based on summary statistics from our GWAS and expression QTL scans, we predicted and replicated a candidate gene AT1G11560 that potentially causes this association. Further analysis indicates that this locus is co-selected with flowering-time-related genes. The discovered pleiotropic genotype-phenotype map provides new insights into understanding the genetic correlation of complex traits.

Genomic, molecular, and cellular divergence of the human brain

While many core biological processes are conserved across species, the human brain has evolved with unique capacities. Current understanding of the neurobiological mechanisms that endow human traits as well as associated vulnerabilities remains limited. However, emerging data have illuminated species divergence in DNA elements and genome organization, in molecular, morphological, and functional features of conserved neural cell types, as well as temporal differences in brain development. Here, we summarize recent data on unique features of the human brain and their complex implications for the study and treatment of brain diseases. We also consider key outstanding questions in the field and discuss the technologies and foundational knowledge that will be required to accelerate understanding of human neurobiology.

Genome-wide association study identifies high-impact susceptibility loci for HCC in North America

BACKGROUND AND AIMS

Despite the substantial impact of environmental factors, individuals with a family history of liver cancer have an increased risk for HCC. However, genetic factors have not been studied systematically by genome-wide approaches in large numbers of individuals from European descent populations (EDP).

APPROACH AND RESULTS

We conducted a 2-stage genome-wide association study (GWAS) on HCC not affected by HBV infections. A total of 1872 HCC cases and 2907 controls were included in the discovery stage, and 1200 HCC cases and 1832 controls in the validation. We analyzed the discovery and validation samples separately and then conducted a meta-analysis. All analyses were conducted in the presence and absence of HCV. The liability-scale heritability was 24.4% for overall HCC. Five regions with significant ORs (95% CI) were identified for nonviral HCC: 3p22.1, MOBP , rs9842969, (0.51, [0.40-0.65]); 5p15.33, TERT , rs2242652, (0.70, (0.62-0.79]); 19q13.11, TM6SF2 , rs58542926, (1.49, [1.29-1.72]); 19p13.11 MAU2 , rs58489806, (1.53, (1.33-1.75]); and 22q13.31, PNPLA3 , rs738409, (1.66, [1.51-1.83]). One region was identified for HCV-induced HCC: 6p21.31, human leukocyte antigen DQ beta 1, rs9275224, (0.79, [0.74-0.84]). A combination of homozygous variants of PNPLA3 and TERT showing a 6.5-fold higher risk for nonviral-related HCC compared to individuals lacking these genotypes. This observation suggests that gene-gene interactions may identify individuals at elevated risk for developing HCC.

CONCLUSIONS

Our GWAS highlights novel genetic susceptibility of nonviral HCC among European descent populations from North America with substantial heritability. Selected genetic influences were observed for HCV-positive HCC. Our findings indicate the importance of genetic susceptibility to HCC development.

A genomic association study revealing subphenotypes of childhood steroid-sensitive nephrotic syndrome in a larger genomic sequencing cohort

Dissecting the genetic components that contribute to the two main subphenotypes of steroid-sensitive nephrotic syndrome (SSNS) using genome-wide association studies (GWAS) strategy is important for understanding the disease. We conducted a multicenter cohort study (360 patients and 1835 controls) combined with a GWAS strategy to identify susceptibility variants associated with the following two subphenotypes of SSNS: steroid-sensitive nephrotic syndrome without relapse (SSNSWR, 181 patients) and steroid-dependent/frequent relapse nephrotic syndrome (SDNS/FRNS, 179 patients). The distribution of two single-nucleotide polymorphisms (SNPs) in ANKRD36 and ALPG was significant between SSNSWR and healthy controls, and that of two SNPs in GAD1 and HLA-DQA1 was significant between SDNS/FRNS and healthy controls. Interestingly, rs1047989 in HLA-DQA1 was a candidate locus for SDNS/FRNS but not for SSNSWR. No significant SNPs were observed between SSNSWR and SDNS/FRNS. Meanwhile, chromosome 2:171713702 in GAD1 was associated with a greater steroid dose (>0.75 mg/kg/d) upon relapse to first remission in patients with SDNS/FRNS (odds ratio = 3.14; 95% confidence interval, 0.97-9.87; P = 0.034). rs117014418 in APOL4 was significantly associated with a decrease in eGFR of greater than 20% compared with the baseline in SDNS/FRNS patients (P = 0.0001). Protein-protein intersection network construction suggested that HLA-DQA1 and HLA-DQB1 function together through GSDMA. Thus, SSNSWR belongs to non-HLA region-dependent nephropathy, and the HLA-DQA/DQB region is likely strongly associated with disease relapse, especially in SDNS/FRNS. The study provides a novel approach for the GWAS strategy of SSNS and contributes to our understanding of the pathological mechanisms of SSNSWR and SDNS/FRNS.

Association between gut microbiota and hepatocellular carcinoma and biliary tract cancer: A mendelian randomization study

BACKGROUND

An increasing number of studies have begun to discuss the relationship between gut microbiota and diseases, yet there is currently a lack of corresponding articles describing the association between gut microbiota and hepatocellular carcinoma (HCC) and biliary tract cancer (BTC). This study aims to explore the relationship between them using Mendelian randomization (MR) analysis method.

AIM

To assess the relationship between gut microbiota and HCC and BTC.

METHODS

We obtained Genome-wide association study (GWAS) data for the gut microbiome from the intestinal microbiota genomic library (MiBioGen, https://mibiogen.gcc.rug.nl/). Additionally, we accessed data pertaining to HCC and BTC from the IEU open GWAS platform (https://gwas.mrcieu.ac.uk/). Our analysis employed fundamental instrumental variable analysis methods, including inverse-variance weighted, MR and Egger. To ensure the dependability of the results, we subjected the results to tests for multiple biases and heterogeneity.

RESULTS

During our investigation, we discovered 11 gut microbiota linked to an increased risk to BTC and HCC. The former included the genus Eubacterium hallii group (P = 0.017), Candidatus Soleaferrea (P = 0.034), Flavonifractor (P = 0.021), Lachnospiraceae FCS020 (P = 0.034), the order Victivallales (P = 0.018), and the class Lentisphaeria (P = 0.0.18). The latter included the genus Desulfovibrio (P = 0.042), Oscillibacter (P = 0.023), the family Coriobacteriaceae (P = 0.048), the order Coriobacteriales (P = 0.048), and the class Coriobacteriia (P = 0.048). Furthermore, in BTC, we observed 2 protective gut microbiota namely the genus Dorea (P = 0.041) and Lachnospiraceae ND3007 group (P = 0.045). All results showed no evidence of multiplicity or heterogeneity.

CONCLUSION

This study explores a causal link between gut microbiota and HCC and BTC. These insights may enhance the mechanistic knowledge of microbiota-related HCC and BTC pathways, potentially informing therapeutic strategies.

Distributed transformer for high order epistasis detection in large-scale datasets

Understanding the genetic basis of complex diseases is one of the most important challenges in current precision medicine. To this end, Genome-Wide Association Studies aim to correlate Single Nucleotide Polymorphisms (SNPs) to the presence or absence of certain traits. However, these studies do not consider interactions between several SNPs, known as epistasis, which explain most genetic diseases. Analyzing SNP combinations to detect epistasis is a major computational task, due to the enormous search space. A possible solution is to employ deep learning strategies for genomic prediction, but the lack of explainability derived from the black-box nature of neural networks is a challenge yet to be addressed. Herein, a novel, flexible, portable, and scalable framework for network interpretation based on transformers is proposed to tackle any-order epistasis. The results on various epistasis scenarios show that the proposed framework outperforms state-of-the-art methods for explainability, while being scalable to large datasets and portable to various deep learning accelerators. The proposed framework is validated on three WTCCC datasets, identifying SNPs related to genes known in the literature that have direct relationships with the studied diseases.

Fast and flexible joint fine-mapping of multiple traits via the Sum of Single Effects model

We introduce mvSuSiE, a multi-trait fine-mapping method for identifying putative causal variants from genetic association data (individual-level or summary data). mvSuSiE learns patterns of shared genetic effects from data, and exploits these patterns to improve power to identify causal SNPs. Comparisons on simulated data show that mvSuSiE is competitive in speed, power and precision with existing multi-trait methods, and uniformly improves on single-trait fine-mapping (SuSiE) in each trait separately. We applied mvSuSiE to jointly fine-map 16 blood cell traits using data from the UK Biobank. By jointly analyzing the traits and modeling heterogeneous effect sharing patterns, we discovered a much larger number of causal SNPs (>3,000) compared with single-trait fine-mapping, and with narrower credible sets. mvSuSiE also more comprehensively characterized the ways in which the genetic variants affect one or more blood cell traits; 68% of causal SNPs showed significant effects in more than one blood cell type.

Machine learning integrative approaches to advance computational immunology

The study of immunology, traditionally reliant on proteomics to evaluate individual immune cells, has been revolutionized by single-cell RNA sequencing. Computational immunologists play a crucial role in analysing these datasets, moving beyond traditional protein marker identification to encompass a more detailed view of cellular phenotypes and their functional roles. Recent technological advancements allow the simultaneous measurements of multiple cellular components-transcriptome, proteome, chromatin, epigenetic modifications and metabolites-within single cells, including in spatial contexts within tissues. This has led to the generation of complex multiscale datasets that can include multimodal measurements from the same cells or a mix of paired and unpaired modalities. Modern machine learning (ML) techniques allow for the integration of multiple "omics" data without the need for extensive independent modelling of each modality. This review focuses on recent advancements in ML integrative approaches applied to immunological studies. We highlight the importance of these methods in creating a unified representation of multiscale data collections, particularly for single-cell and spatial profiling technologies. Finally, we discuss the challenges of these holistic approaches and how they will be instrumental in the development of a common coordinate framework for multiscale studies, thereby accelerating research and enabling discoveries in the computational immunology field.

Epistolution: a new principle necessary to a learning-first theory of life

Biological theory assumes the organized appearance of life and the reliable recurrence of traits are due to inheritance. Natural selection acting on blind variations produces phenotypes with heritable traits, one of which may be natural learning. The aim of learning, then, is solving problems related to survival and reproduction. But what if these views confuse cause with effect? Perhaps a learning algorithm is required for any phenotype at all to arise. If so, evolution proceeds learning-first, with individuals pursuing another telos entirely. I argue that this aim may be epistemological, the drive to understand the world through an umwelt. By "understand" I mean neither association nor prediction but Karl Popper's concept of explanation through conjecture and refutation. I propose that if only genetic materials are truly heritable, not traits, then testing a successful physical theory of life will depend on building abiotic machines which can perform natural learning without the presence of any inherited materials or conditions. I name this process "epistolution," combining "epistemology" and "evolution," to distinguish it from other concepts. Epistolution is an integral consequence of any learning-first view of life, such as the Cellular Basis of Consciousness theory. This type of theory suggests that in all cells during the history of life full-blown agency, involving beliefs, intentions, and desires, generated all the phenotypes that have then been winnowed by natural selection. Unlike in other versions, I posit that the aim of agential living systems is the explanation of reality rather than inductive prediction or survival/reproduction.

A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

tstrait: a quantitative trait simulator for ancestral recombination graphs

SUMMARY

Ancestral recombination graphs (ARGs) encode the ensemble of correlated genealogical trees arising from recombination in a compact and efficient structure and are of fundamental importance in population and statistical genetics. Recent breakthroughs have made it possible to simulate and infer ARGs at biobank scale, and there is now intense interest in using ARG-based methods across a broad range of applications, particularly in genome-wide association studies (GWAS). Sophisticated methods exist to simulate ARGs using population genetics models, but there is currently no software to simulate quantitative traits directly from these ARGs. To apply existing quantitative trait simulators users must export genotype data, losing important information about ancestral processes and producing prohibitively large files when applied to the biobank-scale datasets currently of interest in GWAS. We present tstrait, an open-source Python library to simulate quantitative traits on ARGs, and show how this user-friendly software can quickly simulate phenotypes for biobank-scale datasets on a laptop computer.

AVAILABILITY AND IMPLEMENTATION

tstrait is available for download on the Python Package Index. Full documentation with examples and workflow templates is available on https://tskit.dev/tstrait/docs/, and the development version is maintained on GitHub (https://github.com/tskit-dev/tstrait).

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

A multi-ethnic reference panel to impute HLA classical and non-classical class I alleles in admixed samples: Testing imputation accuracy in an admixed sample from Brazil

The MHC class I region contains crucial genes for the innate and adaptive immune response, playing a key role in susceptibility to many autoimmune and infectious diseases. Genome-wide association studies have identified numerous disease-associated SNPs within this region. However, these associations do not fully capture the immune-biological relevance of specific HLA alleles. HLA imputation techniques may leverage available SNP arrays by predicting allele genotypes based on the linkage disequilibrium between SNPs and specific HLA alleles. Successful imputation requires diverse and large reference panels, especially for admixed populations. This study employed a bioinformatics approach to call SNPs and HLA alleles in multi-ethnic samples from the 1000 genomes (1KG) dataset and admixed individuals from Brazil (SABE), utilising 30X whole-genome sequencing data. Using HIBAG, we created three reference panels: 1KG (n = 2504), SABE (n = 1171), and the full model (n = 3675) encompassing all samples. In extensive cross-validation of these reference panels, the multi-ethnic 1KG reference exhibited overall superior performance than the reference with only Brazilian samples. However, the best results were achieved with the full model. Additionally, we expanded the scope of imputation by developing reference panels for non-classical, MICA, MICB and HLA-H genes, previously unavailable for multi-ethnic populations. Validation in an independent Brazilian dataset showcased the superiority of our reference panels over the Michigan Imputation Server, particularly in predicting HLA-B alleles among Brazilians. Our investigations underscored the need to enhance or adapt reference panels to encompass the target population's genetic diversity, emphasising the significance of multiethnic references for accurate imputation across different populations.

From GWASs toward Mechanistic Understanding with Case Studies in Dermatogenetics

Many human skin diseases result from the complex interplay of genetic and environmental mechanisms that are largely unknown. GWASs have yielded insight into the genetic aspect of complex disease by highlighting regions of the genome or specific genetic variants associated with disease. Leveraging this information to identify causal genes and cell types will provide insight into fundamental biology, inform diagnostics, and aid drug discovery. However, the etiological mechanisms from genetic variant to disease are still unestablished in most cases. There now exists an unprecedented wealth of data and computational methods for variant interpretation in a functional context. It can be challenging to decide where to start owing to a lack of consensus on the best way to identify causal genetic mechanisms. This article highlights 3 key aspects of genetic variant interpretation: prioritizing causal genes, cell types, and pathways. We provide a practical overview of the main methods and datasets, giving examples from recent atopic dermatitis studies to provide a blueprint for variant interpretation. A collection of resources, including brief description and links to the packages and web tools, is provided for researchers looking to start in silico follow-up genetic analysis of associated genetic variants.

Exploring the genetic prediction of academic underachievement and overachievement

Academic underachievement refers to school performance which falls below expectations. Focusing on the pivotal first stage of education, we explored a quantitative measure of underachievement using genomically predicted achievement delta (GPAΔ), which reflects the difference between observed and expected achievement predicted by genome-wide polygenic scores. We analyzed the relationship between GPAΔ at age 7 and achievement trajectories from ages 7 to 16, using longitudinal data from 4175 participants in the Twins Early Development Study to assess empirically the extent to which students regress to their genomically predicted levels by age 16. We found that the achievement of underachievers and overachievers who deviated from their genomic predictions at age 7 regressed on average by one-third towards their genomically predicted levels. We also found that GPAΔ at age 7 was as predictive of achievement trajectories as a traditional ability-based index of underachievement. Targeting GPAΔ underachievers might prove cost-effective because such interventions seem more likely to succeed by going with the genetic flow rather than swimming upstream, helping GPAΔ underachievers reach their genetic potential as predicted by their GPS. However, this is a hypothesis that needs to be tested in intervention research investigating whether GPAΔ underachievers respond better to the intervention than other underachievers. We discuss the practicality of genomic indices in assessing underachievement.

Robust use of phenotypic heterogeneity at drug target genes for mechanistic insights: Application of cis-multivariable Mendelian randomization to GLP1R gene region

Phenotypic heterogeneity at genomic loci encoding drug targets can be exploited by multivariable Mendelian randomization to provide insight into the pathways by which pharmacological interventions may affect disease risk. However, statistical inference in such investigations may be poor if overdispersion heterogeneity in measured genetic associations is unaccounted for. In this work, we first develop conditional F statistics for dimension-reduced genetic associations that enable more accurate measurement of phenotypic heterogeneity. We then develop a novel extension for two-sample multivariable Mendelian randomization that accounts for overdispersion heterogeneity in dimension-reduced genetic associations. Our empirical focus is to use genetic variants in the GLP1R gene region to understand the mechanism by which GLP1R agonism affects coronary artery disease (CAD) risk. Colocalization analyses indicate that distinct variants in the GLP1R gene region are associated with body mass index and type 2 diabetes (T2D). Multivariable Mendelian randomization analyses that were corrected for overdispersion heterogeneity suggest that bodyweight lowering rather than T2D liability lowering effects of GLP1R agonism are more likely contributing to reduced CAD risk. Tissue-specific analyses prioritized brain tissue as the most likely to be relevant for CAD risk, of the tissues considered. We hope the multivariable Mendelian randomization approach illustrated here is widely applicable to better understand mechanisms linking drug targets to diseases outcomes, and hence to guide drug development efforts.

Large-scale whole-exome sequencing of neuropsychiatric diseases and traits in 350,770 adults

While numerous genomic loci have been identified for neuropsychiatric conditions, the contribution of protein-coding variants has yet to be determined. Here we conducted a large-scale whole-exome-sequencing study to interrogate the impact of protein-coding variants on 46 neuropsychiatric diseases and 23 traits in 350,770 adults from the UK Biobank. Twenty new genes were associated with neuropsychiatric diseases through coding variants, among which 16 genes had impacts on the longitudinal risks of diseases. Thirty new genes were associated with neuropsychiatric traits, with SYNGAP1 showing pleiotropic effects across cognitive function domains. Pairwise estimation of genetic correlations at the coding-variant level highlighted shared genetic associations among pairs of neurodegenerative diseases and mental disorders. Lastly, a comprehensive multi-omics analysis suggested that alterations in brain structures, blood proteins and inflammation potentially contribute to the gene-phenotype linkages. Overall, our findings characterized a compendium of protein-coding variants for future research on the biology and therapeutics of neuropsychiatric phenotypes.