Controlling for polygenic genetic confounding in epidemiologic association studies.

Epidemiologic associations estimated from observational data are often confounded by genetics due to pervasive pleiotropy among complex traits. Many studies either neglect genetic confounding altogether or rely on adjusting for polygenic scores (PGS) in regression analysis. In this study, we unveil that the commonly employed PGS approach is inadequate for removing genetic confounding due to measurement error and model misspecification. To tackle this challenge, we introduce PENGUIN, a principled framework for polygenic genetic confounding control based on variance component estimation. In addition, we present extensions of this approach that can estimate genetically unconfounded associations using GWAS summary statistics alone as input and between multiple generations of study samples. Through simulations, we demonstrate superior statistical properties of PENGUIN compared to the existing approaches. Applying our method to multiple population cohorts, we reveal and remove substantial genetic confounding in the associations of educational attainment with various complex traits and between parental and offspring education. Our results show that PENGUIN is an effective solution for genetic confounding control in observational data analysis with broad applications in future epidemiologic association studies.

Decomposing heritability and genetic covariance by direct and indirect effect paths.

Estimation of heritability and genetic covariance is crucial for quantifying and understanding complex trait genetic architecture and is employed in almost all recent genome-wide association studies (GWAS). However, many existing approaches for heritability estimation and almost all methods for estimating genetic correlation ignore the presence of indirect genetic effects, i.e., genotype-phenotype associations confounded by the parental genome and family environment, and may thus lead to incorrect interpretation especially for human sociobehavioral phenotypes. In this work, we introduce a statistical framework to decompose heritability and genetic covariance into multiple components representing direct and indirect effect paths. Applied to five traits in UK Biobank, we found substantial involvement of indirect genetic components in shared genetic architecture across traits. These results demonstrate the effectiveness of our approach and highlight the importance of accounting for indirect effects in variance component analysis of complex traits.

GWAS on birth year infant mortality rates provides evidence of recent natural selection.

Following more than a century of phenotypic measurement of natural selection processes, much recent work explores relationships between molecular genetic measurements and realized fitness in the next generation. We take an innovative approach to the study of contemporary selective pressure by examining which genetic variants are "sustained" in populations as mortality exposure increases. Specifically, we deploy a so-called "regional GWAS" (genome-wide association study) that links the infant mortality rate (IMR) by place and year in the United Kingdom with common genetic variants among birth cohorts in the UK Biobank. These cohorts (born between 1936 and 1970) saw a decline in IMR from above 65 to under 20 deaths per 1,000 live births, with substantial subnational variations and spikes alongside wartime exposures. Our results show several genome-wide significant loci, including LCT and TLR10/1/6, related to area-level cohort IMR exposure during gestation and infancy. Genetic correlations are found across multiple domains, including fertility, cognition, health behaviors, and health outcomes, suggesting an important role for cohort selection in modern populations.

A quantile integral linear model to quantify genetic effects on phenotypic variability.

Detecting genetic variants associated with the variance of complex traits, that is, variance quantitative trait loci (vQTLs), can provide crucial insights into the interplay between genes and environments and how they jointly shape human phenotypes in the population. We propose a quantile integral linear model (QUAIL) to estimate genetic effects on trait variability. Through extensive simulations and analyses of real data, we demonstrate that QUAIL provides computationally efficient and statistically powerful vQTL mapping that is robust to non-Gaussian phenotypes and confounding effects on phenotypic variability. Applied to UK Biobank (n = 375,791), QUAIL identified 11 vQTLs for body mass index (BMI) that have not been previously reported. Top vQTL findings showed substantial enrichment for interactions with physical activities and sedentary behavior. Furthermore, variance polygenic scores (vPGSs) based on QUAIL effect estimates showed superior predictive performance on both population-level and within-individual BMI variability compared to existing approaches. Overall, QUAIL is a unified framework to quantify genetic effects on the phenotypic variability at both single-variant and vPGS levels. It addresses critical limitations in existing approaches and may have broad applications in future gene-environment interaction studies.

Separating Scarring Effect and Selection of Early-Life Exposures With Genetic Data.

Causal life course research examining consequences of early-life exposures has largely relied on associations between early-life environments and later-life outcomes using exogenous environmental shocks. Nonetheless, even with (quasi-)randomized early-life exposures, these associations may reflect not only causation ("scarring") but also selection (i.e., which members are included in data assessing later life). Investigating this selection and its impacts on estimated effects of early-life conditions has, however, often been ignored because of a lack of pre-exposure data. This study proposes an approach for assessing and correcting selection, separately from scarring, using genetic measurements. Because genetic measurements are determined at the time of conception, any associations with early-life exposures should be interpreted as selection. Using data from the UK Biobank, we find that in utero exposure to a higher area-level infant mortality rate is associated with genetic predispositions correlated with better educational attainment and health. These findings point to the direction and magnitude of selection from this exposure. Corrections for this selection in examinations of effects of exposure on later educational attainment suggest underestimates of 26-74%; effects on other life course outcomes also vary across selection correction methods.

Estimating genetic nurture with summary statistics of multigenerational genome-wide association studies.

Marginal effect estimates in genome-wide association studies (GWAS) are mixtures of direct and indirect genetic effects. Existing methods to dissect these effects require family-based, individual-level genetic, and phenotypic data with large samples, which is difficult to obtain in practice. Here, we propose a statistical framework to estimate direct and indirect genetic effects using summary statistics from GWAS conducted on own and offspring phenotypes. Applied to birth weight, our method showed nearly identical results with those obtained using individual-level data. We also decomposed direct and indirect genetic effects of educational attainment (EA), which showed distinct patterns of genetic correlations with 45 complex traits. The known genetic correlations between EA and higher height, lower body mass index, less-active smoking behavior, and better health outcomes were mostly explained by the indirect genetic component of EA. In contrast, the consistently identified genetic correlation of autism spectrum disorder (ASD) with higher EA resides in the direct genetic component. A polygenic transmission disequilibrium test showed a significant overtransmission of the direct component of EA from healthy parents to ASD probands. Taken together, we demonstrate that traditional GWAS approaches, in conjunction with offspring phenotypic data collection in existing cohorts, could greatly benefit studies on genetic nurture and shed important light on the interpretation of genetic associations for human complex traits.

Transcriptome-wide transmission disequilibrium analysis identifies novel risk genes for autism spectrum disorder.

Recent advances in consortium-scale genome-wide association studies (GWAS) have highlighted the involvement of common genetic variants in autism spectrum disorder (ASD), but our understanding of their etiologic roles, especially the interplay with rare variants, is incomplete. In this work, we introduce an analytical framework to quantify the transmission disequilibrium of genetically regulated gene expression from parents to offspring. We applied this framework to conduct a transcriptome-wide association study (TWAS) on 7,805 ASD proband-parent trios, and replicated our findings using 35,740 independent samples. We identified 31 associations at the transcriptome-wide significance level. In particular, we identified POU3F2 (p = 2.1E-7), a transcription factor mainly expressed in developmental brain. Gene targets regulated by POU3F2 showed a 2.7-fold enrichment for known ASD genes (p = 2.0E-5) and a 2.7-fold enrichment for loss-of-function de novo mutations in ASD probands (p = 7.1E-5). These results provide a novel connection between rare and common variants, whereby ASD genes affected by very rare mutations are regulated by an unlinked transcription factor affected by common genetic variations.

Apolipoprotein E moderates the association between non-APOE polygenic risk score for Alzheimer's disease and aging on preclinical cognitive function.

INTRODUCTION: Variation in preclinical cognitive decline suggests additional genetic factors related to Alzheimer's disease (eg, a non-APOE polygenic risk score [PRS]) may interact with the APOE ε4 allele to influence cognitive decline. METHODS: We tested the PRS × APOE ε4 × age interaction on preclinical cognition using longitudinal data from the Wisconsin Registry for Alzheimer's Prevention. All analyses were fitted using a linear mixed-effects model and adjusted for within individual/family correlation among 1190 individuals. RESULTS: We found statistically significant PRS × APOE ε4 × age interactions on immediate learning (P = 0.038), delayed recall (P < 0.001), and Preclinical Alzheimer's Cognitive Composite 3 score (P = 0.026). PRS-related differences in overall and memory-related cognitive domains between people with and without APOE ε4 emerge around age 70, with a much stronger adverse PRS effect among APOE ε4 carriers. The findings were replicated in a population-based cohort. DISCUSSIONS: APOE ε4 can modify the association between PRS and cognition decline. HIGHLIGHTS: APOE ε4 can modify the association between polygenic risk scores (PRSs) and longitudinal cognition decline, with the modifying effects more pronounced when the PRS is constructed using a conservative P threshold (eg, P < 5e-8 ). The adverse genetic effect caused by the combined effect of the currently known genetic variants is more detrimental among APOE Îµ4 carriers around age 70. Individuals who are APOE Îµ4 carriers with high PRSs are the most vulnerable to the harmful effects caused by genetic burden.

Post-GWAS multiomic functional investigation of the TNIP1 locus in Alzheimer's disease highlights a potential role for GPX3.

INTRODUCTION: Recent genome-wide association studies (GWAS) have reported a genetic association with Alzheimer's disease (AD) at the TNIP1/GPX3 locus, but the mechanism is unclear. METHODS: We used cerebrospinal fluid (CSF) proteomics data to test (n = 137) and replicate (n = 446) the association of glutathione peroxidase 3 (GPX3) with CSF biomarkers (including amyloid and tau) and the GWAS-implicated variants (rs34294852 and rs871269). RESULTS: CSF GPX3 levels decreased with amyloid and tau positivity (analysis of variance P = 1.5 × 10-5) and higher CSF phosphorylated tau (p-tau) levels (P = 9.28 × 10-7). The rs34294852 minor allele was associated with decreased GPX3 (P = 0.041). The replication cohort found associations of GPX3 with amyloid and tau positivity (P = 2.56 × 10-6) and CSF p-tau levels (P = 4.38 × 10-9). DISCUSSION: These results suggest variants in the TNIP1 locus may affect the oxidative stress response in AD via altered GPX3 levels. HIGHLIGHTS: Cerebrospinal fluid (CSF) glutathione peroxidase 3 (GPX3) levels decreased with amyloid and tau positivity and higher CSF phosphorylated tau. The minor allele of rs34294852 was associated with lower CSF GPX3. levels when also controlling for amyloid and tau category. GPX3 transcript levels in the prefrontal cortex were lower in Alzheimer's disease than controls. rs34294852 is an expression quantitative trait locus for GPX3 in blood, neutrophils, and microglia.

Comparison of methods for estimating genetic correlation between complex traits using GWAS summary statistics.

Genetic correlation is the correlation of phenotypic effects by genetic variants across the genome on two phenotypes. It is an informative metric to quantify the overall genetic similarity between complex traits, which provides insights into their polygenic genetic architecture. Several methods have been proposed to estimate genetic correlation based on data collected from genome-wide association studies (GWAS). Due to the easy access of GWAS summary statistics and computational efficiency, methods only requiring GWAS summary statistics as input have become more popular than methods utilizing individual-level genotype data. Here, we present a benchmark study for different summary-statistics-based genetic correlation estimation methods through simulation and real data applications. We focus on two major technical challenges in estimating genetic correlation: marker dependency caused by linkage disequilibrium (LD) and sample overlap between different studies. To assess the performance of different methods in the presence of these two challenges, we first conducted comprehensive simulations with diverse LD patterns and sample overlaps. Then we applied these methods to real GWAS summary statistics for a wide spectrum of complex traits. Based on these experiments, we conclude that methods relying on accurate LD estimation are less robust in real data applications due to the imprecision of LD obtained from reference panels. Our findings offer guidance on how to choose appropriate methods for genetic correlation estimation in post-GWAS analysis.

Sex differences in the genetic architecture of cognitive resilience to Alzheimer's disease.

Approximately 30% of elderly adults are cognitively unimpaired at time of death despite the presence of Alzheimer's disease neuropathology at autopsy. Studying individuals who are resilient to the cognitive consequences of Alzheimer's disease neuropathology may uncover novel therapeutic targets to treat Alzheimer's disease. It is well established that there are sex differences in response to Alzheimer's disease pathology, and growing evidence suggests that genetic factors may contribute to these differences. Taken together, we sought to elucidate sex-specific genetic drivers of resilience. We extended our recent large scale genomic analysis of resilience in which we harmonized cognitive data across four cohorts of cognitive ageing, in vivo amyloid PET across two cohorts, and autopsy measures of amyloid neuritic plaque burden across two cohorts. These data were leveraged to build robust, continuous resilience phenotypes. With these phenotypes, we performed sex-stratified [n (males) = 2093, n (females) = 2931] and sex-interaction [n (both sexes) = 5024] genome-wide association studies (GWAS), gene and pathway-based tests, and genetic correlation analyses to clarify the variants, genes and molecular pathways that relate to resilience in a sex-specific manner. Estimated among cognitively normal individuals of both sexes, resilience was 20-25% heritable, and when estimated in either sex among cognitively normal individuals, resilience was 15-44% heritable. In our GWAS, we identified a female-specific locus on chromosome 10 [rs827389, ß (females) = 0.08, P (females) = 5.76 × 10-09, ß (males) = -0.01, P(males) = 0.70, ß (interaction) = 0.09, P (interaction) = 1.01 × 10-04] in which the minor allele was associated with higher resilience scores among females. This locus is located within chromatin loops that interact with promoters of genes involved in RNA processing, including GATA3. Finally, our genetic correlation analyses revealed shared genetic architecture between resilience phenotypes and other complex traits, including a female-specific association with frontotemporal dementia and male-specific associations with heart rate variability traits. We also observed opposing associations between sexes for multiple sclerosis, such that more resilient females had a lower genetic susceptibility to multiple sclerosis, and more resilient males had a higher genetic susceptibility to multiple sclerosis. Overall, we identified sex differences in the genetic architecture of resilience, identified a female-specific resilience locus and highlighted numerous sex-specific molecular pathways that may underly resilience to Alzheimer's disease pathology. This study illustrates the need to conduct sex-aware genomic analyses to identify novel targets that are unidentified in sex-agnostic models. Our findings support the theory that the most successful treatment for an individual with Alzheimer's disease may be personalized based on their biological sex and genetic context.

Identification of plasma metabolites associated with modifiable risk factors and endophenotypes reflecting Alzheimer's disease pathology.

Modifiable factors can influence the risk for Alzheimer's disease (AD) and serve as targets for intervention; however, the biological mechanisms linking these factors to AD are unknown. This study aims to identify plasma metabolites associated with modifiable factors for AD, including MIND diet, physical activity, smoking, and caffeine intake, and test their association with AD endophenotypes to identify their potential roles in pathophysiological mechanisms. The association between each of the 757 plasma metabolites and four modifiable factors was tested in the wisconsin registry for Alzheimer's prevention cohort of initially cognitively unimpaired, asymptomatic middle-aged adults. After Bonferroni correction, the significant plasma metabolites were tested for association with each of the AD endophenotypes, including twelve cerebrospinal fluid (CSF) biomarkers, reflecting key pathophysiologies for AD, and four cognitive composite scores. Finally, causal mediation analyses were conducted to evaluate possible mediation effects. Analyses were performed using linear mixed-effects regression. A total of 27, 3, 23, and 24 metabolites were associated with MIND diet, physical activity, smoking, and caffeine intake, respectively. Potential mediation effects include beta-cryptoxanthin in the association between MIND diet and preclinical Alzheimer cognitive composite score, hippurate between MIND diet and immediate learning, glutamate between physical activity and CSF neurofilament light, and beta-cryptoxanthin between smoking and immediate learning. Our study identified several plasma metabolites that are associated with modifiable factors. These metabolites can be employed as biomarkers for tracking these factors, and they provide a potential biological pathway of how modifiable factors influence the human body and AD risk.

The Big (Genetic) Sort? A Research Note on Migration Patterns and Their Genetic Imprint in the United Kingdom.

This research note reinvestigates Abdellaoui et al.'s (2019) findings that genetically selective migration may lead to persistent and accumulating socioeconomic and health inequalities between types (coal mining or non-coal mining) of places in the United Kingdom. Their migration measure classified migrants who moved to the same type of place (coal mining to coal mining or non-coal mining to non-coal mining) into "stay" categories, preventing them from distinguishing migrants from nonmigrants. We reinvestigate the question of genetically selective migration by examining migration patterns between places rather than place types and find genetic selectivity in whether people migrate and where. For example, we find evidence of positive selection: people with genetic variants correlated with better education moved from non-coal mining to coal mining places with our measure of migration. Such findings were obscured in earlier work that could not distinguish nonmigrants from migrants.

Understanding Internal Migration: A Research Note Providing an Assessment of Migration Selection With Genetic Data.

Migration is selective, resulting in inequalities between migrants and nonmigrants. However, investigating migration selection is empirically challenging because combined pre- and post-migration data are rarely available. We propose an alternative approach to assessing internal migration selection by integrating genetic data, enabling an investigation of migration selection with cross-sectional data collected post-migration. Using data from the UK Biobank, we utilized standard tools from statistical genetics to conduct a genome-wide association study (GWAS) for migration distance. We then calculated genetic correlations to compare GWAS results for migration with those for other characteristics. Given that individual genetics are determined at conception, these analyses allow a unique exploration of the association between pre-migration characteristics and migration. Results are generally consistent with the healthy migrant literature: genetics correlated with longer migration distance are associated with higher socioeconomic status and better health. We also extended the analysis to 53 traits and found novel correlations between migration and several physical health, mental health, personality, and sociodemographic traits.

Neuropathology-based APOE genetic risk score better quantifies Alzheimer's risk.

INTRODUCTION: Apolipoprotein E (APOE) ε4-carrier status or ε4 allele count are included in analyses to account for the APOE genetic effect on Alzheimer's disease (AD); however, this does not account for protective effects of APOE ε2 or heterogeneous effect of ε2, ε3, and ε4 haplotypes. METHODS: We leveraged results from an autopsy-confirmed AD study to generate a weighted risk score for APOE (APOE-npscore). We regressed cerebrospinal fluid (CSF) amyloid and tau biomarkers on APOE variables from the Wisconsin Registry for Alzheimer's Prevention (WRAP), Wisconsin Alzheimer's Disease Research Center (WADRC), and Alzheimer's Disease Neuroimaging Initiative (ADNI). RESULTS: The APOE-npscore explained more variance and provided a better model fit for all three CSF measures than APOE ε4-carrier status and ε4 allele count. These findings were replicated in ADNI and observed in subsets of cognitively unimpaired (CU) participants. DISCUSSION: The APOE-npscore reflects the genetic effect on neuropathology and provides an improved method to account for APOE in AD-related analyses.

Re-examining the relationship between education and adult mental health in the UK: A research note.

Previous studies using variation in education arising from compulsory schooling laws have found no causal effects of education on mental health in the UK. We re-examine the relationship between education and mental health in the UK by taking a different approach: sibling fixed-effects with controls for polygenic scores (summary measures of genetic predisposition) for educational attainment and adult depressive symptoms. We find that higher educational attainment is associated with better adult mental health, that sibling controls reduce these associations by ~40-70% but important associations remain and find evidence for non-monotonic effects. We also find suggestive evidence that education partially "rescues" genetic predictors of poor mental health.

Genome-wide association study reveals sex-specific genetic architecture of facial attractiveness.

Facial attractiveness is a complex human trait of great interest in both academia and industry. Literature on sociological and phenotypic factors associated with facial attractiveness is rich, but its genetic basis is poorly understood. In this paper, we conducted a genome-wide association study to discover genetic variants associated with facial attractiveness using 4,383 samples in the Wisconsin Longitudinal Study. We identified two genome-wide significant loci, highlighted a handful of candidate genes, and demonstrated enrichment for heritability in human tissues involved in reproduction and hormone synthesis. Additionally, facial attractiveness showed strong and negative genetic correlations with BMI in females and with blood lipids in males. Our analysis also suggested sex-specific selection pressure on variants associated with lower male attractiveness. These results revealed sex-specific genetic architecture of facial attractiveness and provided fundamental new insights into its genetic basis.

Multi-ethnic analysis shows genetic risk and environmental predictors interact to influence 25(OH)D concentration and optimal vitamin D intake.

25-Hydroxyvitamin D (25(OH)D) concentration is a complex trait with genetic and environmental predictors that may determine how much vitamin D exposure is required to reach optimal concentration. Interactions between continuous measures of a polygenic score (PGS) and vitamin D intake (PGS*intake) or available ultraviolet (UV) radiation (PGS*UV) were evaluated in individuals of African (n = 1,099) or European (n = 8,569) ancestries. Interaction terms and joint effects (main and interaction terms) were tested using one-degree of freedom (1-DF) and 2-DF models, respectively. Models controlled for age, sex, body mass index, cohort, and dietary intake/available UV. In addition, in participants achieving Institute of Medicine (IOM) vitamin D intake recommendations, 25(OH)D was evaluated by level PGS. The 2-DF PGS*intake, 1-DF PGS*UV, and 2-DF PGS*UV results were statistically significant in participants of European ancestry (p = 3.3 × 10-18 , p = 2.1 × 10-2 , and p = 2.4 × 10-19 , respectively), but not in those of African ancestry. In European-ancestry participants reaching IOM vitamin D intake guidelines, the percent of participants achieving adequate 25(OH)D ( >20 ng/ml) increased as genetic risk decreased (72% vs. 89% in highest vs. lowest risk; p = .018). Available UV radiation and vitamin D intake interact with genetics to influence 25(OH)D. Individuals with higher genetic risk may require more vitamin D exposure to maintain optimal 25(OH)D concentrations.

An evaluation of noncoding genome annotation tools through enrichment analysis of 15 genome-wide association studies.

Functionally annotating genetic variations is an essential yet challenging topic in human genetics research. As large consortia including ENCODE and Roadmap Epigenomics Project continue to generate high-throughput transcriptomic and epigenomic data, many computational frameworks have been developed to integrate these experimental data to predict functionality of genetic variations in both protein-coding and noncoding regions. Here, we compare a number of recently developed annotation frameworks for noncoding regions through enrichment analysis on genome-wide association studies (GWASs). We also compare several different strategies to quantify enrichment using GWAS summary statistics. Our analyses highlight the importance of jointly modeling context-specific annotations with genome-wide data in providing statistically powerful and biologically interpretable enrichment for complex disease associations. Our findings provide insights into when and how computational genome annotations may benefit future complex disease studies on the genome-wide scale.

Leveraging functional annotation to identify genes associated with complex diseases.

To increase statistical power to identify genes associated with complex traits, a number of transcriptome-wide association study (TWAS) methods have been proposed using gene expression as a mediating trait linking genetic variations and diseases. These methods first predict expression levels based on inferred expression quantitative trait loci (eQTLs) and then identify expression-mediated genetic effects on diseases by associating phenotypes with predicted expression levels. The success of these methods critically depends on the identification of eQTLs, which may not be functional in the corresponding tissue, due to linkage disequilibrium (LD) and the correlation of gene expression between tissues. Here, we introduce a new method called T-GEN (Transcriptome-mediated identification of disease-associated Genes with Epigenetic aNnotation) to identify disease-associated genes leveraging epigenetic information. Through prioritizing SNPs with tissue-specific epigenetic annotation, T-GEN can better identify SNPs that are both statistically predictive and biologically functional. We found that a significantly higher percentage (an increase of 18.7% to 47.2%) of eQTLs identified by T-GEN are inferred to be functional by ChromHMM and more are deleterious based on their Combined Annotation Dependent Depletion (CADD) scores. Applying T-GEN to 207 complex traits, we were able to identify more trait-associated genes (ranging from 7.7% to 102%) than those from existing methods. Among the identified genes associated with these traits, T-GEN can better identify genes with high (>0.99) pLI scores compared to other methods. When T-GEN was applied to late-onset Alzheimer's disease, we identified 96 genes located at 15 loci, including two novel loci not implicated in previous GWAS. We further replicated 50 genes in an independent GWAS, including one of the two novel loci.