Univariate Genome-wide Association Studies Research Articles

Identification of genetic basis of brain imaging by group sparse multi-task learning leveraging summary statistics

Brain imaging genetics is an evolving neuroscience topic aiming to identify genetic variations related to neuroimaging measurements of interest. Traditional linear regression methods have shown success, but their reliance on individual-level imaging and genetic data limits their applicability. Herein, we proposed S-GsMTLR, a group sparse multi-task linear regression method designed to harness summary statistics from genome-wide association studies (GWAS) of neuroimaging quantitative traits. S-GsMTLR directly employs GWAS summary statistics, bypassing the requirement for raw imaging genetic data, and applies multivariate multi-task sparse learning to these univariate GWAS results. It amalgamates the strengths of conventional sparse learning methods, including sophisticated modeling techniques and efficient feature selection. Additionally, we implemented a rapid optimization strategy to alleviate computational burdens by identifying genetic variants associated with phenotypes of interest across the entire chromosome. We first evaluated S-GsMTLR using summary statistics derived from the Alzheimer's Disease Neuroimaging Initiative. The results were remarkably encouraging, demonstrating its comparability to conventional methods in modeling and identification of risk loci. Furthermore, our method was evaluated with two additional GWAS summary statistics datasets: One focused on white matter microstructures and the other on whole brain imaging phenotypes, where the original individual-level data was unavailable. The results not only highlighted S-GsMTLR's ability to pinpoint significant loci but also revealed intriguing structures within genetic variations and loci that went unnoticed by GWAS. These findings suggest that S-GsMTLR is a promising multivariate sparse learning method in brain imaging genetics. It eliminates the need for original individual-level imaging and genetic data while demonstrating commendable modeling and feature selection capabilities.

Multi-trait GWAS for diverse ancestries: mapping the knowledge gap

BackgroundApproximately 95% of samples analyzed in univariate genome-wide association studies (GWAS) are of European ancestry. This bias toward European ancestry populations in association screening also exists for other analyses and methods that are often developed and tested on European ancestry only. However, existing data in non-European populations, which are often of modest sample size, could benefit from innovative approaches as recently illustrated in the context of polygenic risk scores.MethodsHere, we extend and assess the potential limitations and gains of our multi-trait GWAS pipeline, JASS (Joint Analysis of Summary Statistics), for the analysis of non-European ancestries. To this end, we conducted the joint GWAS of 19 hematological traits and glycemic traits across five ancestries (European (EUR), admixed American (AMR), African (AFR), East Asian (EAS), and South-East Asian (SAS)).ResultsWe detected 367 new genome-wide significant associations in non-European populations (15 in Admixed American (AMR), 72 in African (AFR) and 280 in East Asian (EAS)). New associations detected represent 5%, 17% and 13% of associations in the AFR, AMR and EAS populations, respectively. Overall, multi-trait testing increases the replication of European associated loci in non-European ancestry by 15%. Pleiotropic effects were highly similar at significant loci across ancestries (e.g. the mean correlation between multi-trait genetic effects of EUR and EAS ancestries was 0.88). For hematological traits, strong discrepancies in multi-trait genetic effects are tied to known evolutionary divergences: the ARKC1 loci, which is adaptive to overcome p.vivax induced malaria.ConclusionsMulti-trait GWAS can be a valuable tool to narrow the genetic knowledge gap between European and non-European populations.

Multivariate canonical correlation analysis identifies additional genetic variants for chronic kidney disease

Chronic kidney diseases (CKD) have genetic associations with kidney function. Univariate genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN), two complementary kidney function markers. However, it is unknown whether additional SNPs for kidney function can be identified by multivariate statistical analysis. To address this, we applied canonical correlation analysis (CCA), a multivariate method, to two individual-level CKD genotype datasets, and metaCCA to two published GWAS summary statistics datasets. We identified SNPs previously associated with kidney function by published univariate GWASs with high replication rates, validating the metaCCA method. We then extended discovery and identified previously unreported lead SNPs for both kidney function markers, jointly. These showed expression quantitative trait loci (eQTL) colocalisation with genes having significant differential expression between CKD and healthy individuals. Several of these identified lead missense SNPs were predicted to have a functional impact, including in SLC14A2. We also identified previously unreported lead SNPs that showed significant correlation with both kidney function markers, jointly, in the European ancestry CKDGen, National Unified Renal Translational Research Enterprise (NURTuRE)-CKD and Salford Kidney Study (SKS) datasets. Of these, rs3094060 colocalised with FLOT1 gene expression and was significantly more common in CKD cases in both NURTURE-CKD and SKS, than in the general population. Overall, by using multivariate analysis by CCA, we identified additional SNPs and genes for both kidney function and CKD, that can be prioritised for further CKD analyses.

Principal and independent genomic components of brain structure and function.

The highly polygenic and pleiotropic nature of behavioural traits, psychiatric disorders and structural and functional brain phenotypes complicate mechanistic interpretation of related genome-wide association study (GWAS) signals, thereby obscuring underlying causal biological processes. We propose genomic principal and independent component analysis (PCA, ICA) to decompose a large set of univariate GWAS statistics of multimodal brain traits into more interpretable latent genomic components. Here we introduce and evaluate this novel methods various analytic parameters and reproducibility across independent samples. Two UK Biobank GWAS summary statistic releases of 2240 imaging-derived phenotypes (IDPs) were retrieved. Genome-wide beta-values and their corresponding standard-error scaled z-values were decomposed using genomic PCA/ICA. We evaluated variance explained at multiple dimensions up to 200. We tested the inter-sample reproducibility of output of dimensions 5, 10, 25 and 50. Reproducibility statistics of the respective univariate GWAS served as benchmarks. Reproducibility of 10-dimensional PCs and ICs showed the best trade-off between model complexity and robustness and variance explained (PCs: |rz - max| = 0.33, |rraw - max| = 0.30; ICs: |rz - max| = 0.23, |rraw - max| = 0.19). Genomic PC and IC reproducibility improved substantially relative to mean univariate GWAS reproducibility up to dimension 10. Genomic components clustered along neuroimaging modalities. Our results indicate that genomic PCA and ICA decompose genetic effects on IDPs from GWAS statistics with high reproducibility by taking advantage of the inherent pleiotropic patterns. These findings encourage further applications of genomic PCA and ICA as fully data-driven methods to effectively reduce the dimensionality, enhance the signal to noise ratio and improve interpretability of high-dimensional multitrait genome-wide analyses.

Guidelines for Evaluating the Comparability of Down-Sampled GWAS Summary Statistics

Proprietary genetic datasets are valuable for boosting the statistical power of genome-wide association studies (GWASs), but their use can restrict investigators from publicly sharing the resulting summary statistics. Although researchers can resort to sharing down-sampled versions that exclude restricted data, down-sampling reduces power and might change the genetic etiology of the phenotype being studied. These problems are further complicated when using multivariate GWAS methods, such as genomic structural equation modeling (Genomic SEM), that model genetic correlations across multiple traits. Here, we propose a systematic approach to assess the comparability of GWAS summary statistics that include versus exclude restricted data. Illustrating this approach with a multivariate GWAS of an externalizing factor, we assessed the impact of down-sampling on (1) the strength of the genetic signal in univariate GWASs, (2) the factor loadings and model fit in multivariate Genomic SEM, (3) the strength of the genetic signal at the factor level, (4) insights from gene-property analyses, (5) the pattern of genetic correlations with other traits, and (6) polygenic score analyses in independent samples. For the externalizing GWAS, although down-sampling resulted in a loss of genetic signal and fewer genome-wide significant loci; the factor loadings and model fit, gene-property analyses, genetic correlations, and polygenic score analyses were found robust. Given the importance of data sharing for the advancement of open science, we recommend that investigators who generate and share down-sampled summary statistics report these analyses as accompanying documentation to support other researchers’ use of the summary statistics.

A multivariate genome-wide association study of psycho-cardiometabolic multimorbidity.

Coronary artery disease (CAD), type 2 diabetes (T2D) and depression are among the leading causes of chronic morbidity and mortality worldwide. Epidemiological studies indicate a substantial degree of multimorbidity, which may be explained by shared genetic influences. However, research exploring the presence of pleiotropic variants and genes common to CAD, T2D and depression is lacking. The present study aimed to identify genetic variants with effects on cross-trait liability to psycho-cardiometabolic diseases. We used genomic structural equation modelling to perform a multivariate genome-wide association study of multimorbidity (Neffective = 562,507), using summary statistics from univariate genome-wide association studies for CAD, T2D and major depression. CAD was moderately genetically correlated with T2D (rg = 0.39, P = 2e-34) and weakly correlated with depression (rg = 0.13, P = 3e-6). Depression was weakly correlated with T2D (rg = 0.15, P = 4e-15). The latent multimorbidity factor explained the largest proportion of variance in T2D (45%), followed by CAD (35%) and depression (5%). We identified 11 independent SNPs associated with multimorbidity and 18 putative multimorbidity-associated genes. We observed enrichment in immune and inflammatory pathways. A greater polygenic risk score for multimorbidity in the UK Biobank (N = 306,734) was associated with the co-occurrence of CAD, T2D and depression (OR per standard deviation = 1.91, 95% CI = 1.74-2.10, relative to the healthy group), validating this latent multimorbidity factor. Mendelian randomization analyses suggested potentially causal effects of BMI, body fat percentage, LDL cholesterol, total cholesterol, fasting insulin, income, insomnia, and childhood maltreatment. These findings advance our understanding of multimorbidity suggesting common genetic pathways.

The genetic architecture of the corpus callosum and its genetic overlap with common neuropsychiatric diseases

BackgroundThe corpus callosum (CC) is the main structure transferring information between the cerebral hemispheres. Although previous large-scale genome-wide association study (GWAS) has illustrated the genetic architecture of white matter integrity of CC, CC volume is less stressed. MethodsUsing MRI data from 33,861 individuals in UK Biobank, we conducted univariate and multivariate GWAS for CC fractional anisotropy (FA) and volume with PLINK 2.0 and MOSTest. All discovered SNPs in the multivariate framework were functionally annotated in FUMA v1.3.8. In the meanwhile, a series of gene property analyses was conducted simultaneously. In addition, we estimated genetic relationship between CC metrics and other neuropsychiatric traits and diseases. ResultsWe identified a total of 36 and 82 significant genomic loci for CC FA and volume (P < 5 × 10−8). And 53 and 27 genes were respectively mapped by four mapping strategies. For CC volume, gene-set analysis revealed pathways mainly relating to cell migration; cell-type analysis found the top enrichment in neuroglia while for CC FA in GABAergic neurons. Furthermore, we found a lot of genetic overlap and shared loci between CC FA and volume and common neuropsychiatric diseases. DiscussionCollectively, this study helps to better understand the genetic architecture of whole CC and CC subregions. However, the way to divide CC FA and volume in our study restricts the interpretations of our results. Future work will be needed to pay attention to the genetic structure of white matter volume, and an appropriate division of CC may help to better understand CC structure.

Univariate and multivariate genome-wide association studies for hematological traits in Murciano-Granadina goats.

Hematological traits are important indicators of health status, and they are frequently used as criteria for clinical diagnosis. In humans, the genomic architecture of blood traits has been investigated in depth and thousands of associations with genetic variants have been found. In contrast, the association between marker genotypes and the variation of hematological traits has not been investigated in goats yet. Herewith, we have recorded 12 hematological parameters in 882 Murciano-Granadina goats that were also genotyped with the Goat SNP50 BeadChip (Illumina). Performance of a univariate genome-wide association study (GWAS) made it possible to detect one genomic region on goat chromosome (CHI) 21 (19.2-19.5 Mb) associated, at the genome-wide level of significance, with 4 red blood cell traits. The three markers displaying the highest significances were rs268272996 (CHI21: 19225290 bp), rs268273004 (CHI21: 19565629 bp) and rs268239059 (CHI13: 9615190 bp). Consistently, a multivariate GWAS indicated that the rs268273004 marker on chromosome 21 is associated with seven blood cell traits. Interestingly, this marker maps close to the FA Complementation Group I (FANCI) gene (CHI21: 20021947-20077025 bp), which is functionally related to Fanconi anemia, a syndrome characterized by bone marrow failure, aplastic anemia, and congenital disorders. We have also uncovered additional chromosome-wide significant associations between genetic markers and erythrocyte and leukocyte traits in the univariate GWAS. These findings evidence that the phenotypic variation of hematological traits in goats is regulated, at least to some extent, by polygenic determinants distributed in multiple chromosomes.

The genetic architecture of human amygdala volumes and their overlap with common brain disorders

The amygdala is a crucial interconnecting structure in the brain that performs several regulatory functions, yet its genetic architectures and involvement in brain disorders remain largely unknown. We carried out the first multivariate genome-wide association study (GWAS) of amygdala subfield volumes in 27,866 UK Biobank individuals. The whole amygdala was segmented into nine nuclei groups using Bayesian amygdala segmentation. The post-GWAS analysis allowed us to identify causal genetic variants in phenotypes at the SNP, locus, and gene levels, as well as genetic overlap with brain health-related traits. We further generalized our GWAS in Adolescent Brain Cognitive Development (ABCD) cohort. The multivariate GWAS identified 98 independent significant variants within 32 genomic loci associated (P < 5 × 10−8) with amygdala volume and its nine nuclei. The univariate GWAS identified significant hits for eight of the ten volumes, tagging 14 independent genomic loci. Overall, 13 of the 14 loci identified in the univariate GWAS were replicated in the multivariate GWAS. The generalization in ABCD cohort supported the GWAS results with the 12q23.2 (RNA gene RP11-210L7.1) being discovered. All of these imaging phenotypes are heritable, with heritability ranging from 15% to 27%. Gene-based analyses revealed pathways relating to cell differentiation/development and ion transporter/homeostasis, with the astrocytes found to be significantly enriched. Pleiotropy analyses revealed shared variants with neurological and psychiatric disorders under the conjFDR threshold of 0.05. These findings advance our understanding of the complex genetic architectures of amygdala and their relevance in neurological and psychiatric disorders.

Multivariate GWAS analysis reveals loci associated with liver functions in continental African populations.

Liver disease is any condition that causes liver damage and inflammation and may likely affect the function of the liver. Vital biochemical screening tools that can be used to evaluate the health of the liver and help diagnose, prevent, monitor, and control the development of liver disease are known as liver function tests (LFT). LFTs are performed to estimate the level of liver biomarkers in the blood. Several factors are associated with differences in concentration levels of LFTs in individuals, such as genetic and environmental factors. The aim of our study was to identify genetic loci associated with liver biomarker levels with a shared genetic basis in continental Africans, using a multivariate genome-wide association study (GWAS) approach. We used two distinct African populations, the Ugandan Genome Resource (UGR = 6,407) and South African Zulu cohort (SZC = 2,598). The six LFTs used in our analysis were: aspartate transaminase (AST), alanine transaminase (ALT), alkaline phosphatase (ALP), gamma-glutamyl transferase (GGT), total bilirubin, and albumin. A multivariate GWAS of LFTs was conducted using the exact linear mixed model (mvLMM) approach implemented in GEMMA and the resulting P-values were presented in Manhattan and quantile-quantile (QQ) plots. First, we attempted to replicate the findings of the UGR cohort in SZC. Secondly, given that the genetic architecture of UGR is different from that of SZC, we further undertook similar analysis in the SZC and discussed the results separately. A total of 59 SNPs reached genome-wide significance (P = 5x10-8) in the UGR cohort and with 13 SNPs successfully replicated in SZC. These included a novel lead SNP near the RHPN1 locus (lead SNP rs374279268, P-value = 4.79x10-9, Effect Allele Frequency (EAF) = 0.989) and a lead SNP at the RGS11 locus (lead SNP rs148110594, P-value = 2.34x10-8, EAF = 0.928). 17 SNPs were significant in the SZC, while all the SNPs fall within a signal on chromosome 2, rs1976391 mapped to UGT1A was identified as the lead SNP within this region. Using multivariate GWAS method improves the power to detect novel genotype-phenotype associations for liver functions not found with the standard univariate GWAS in the same dataset.

Predicting Gene Regulatory Interactions Using Natural Genetic Variation.

Genome-wide association studies (GWAS) are a powerful tool to elucidate the genotype-phenotype map. Although GWAS are usually used to assess simple univariate associations between genetic markers and traits of interest, it is also possible to infer the underlying genetic architecture and to predict gene regulatory interactions. In this chapter, we describe the latest methods and tools to perform GWAS by calculating permutation-based significance thresholds. For this purpose, we first provide guidelines on univariate GWAS analyses that are extended in the second part of this chapter to more complex models that enable the inference of gene regulatory networks and how these networks vary.

The link between liver fat and cardiometabolic diseases is highlighted by genome-wide association study of MRI-derived measures of body composition

Obesity and associated morbidities, metabolic associated fatty liver disease (MAFLD) included, constitute some of the largest public health threats worldwide. Body composition and related risk factors are known to be heritable and identification of their genetic determinants may aid in the development of better prevention and treatment strategies. Recently, large-scale whole-body MRI data has become available, providing more specific measures of body composition than anthropometrics such as body mass index. Here, we aimed to elucidate the genetic architecture of body composition, by conducting genome-wide association studies (GWAS) of these MRI-derived measures. We ran both univariate and multivariate GWAS on fourteen MRI-derived measurements of adipose and muscle tissue distribution, derived from scans from 33,588 White European UK Biobank participants (mean age of 64.5 years, 51.4% female). Through multivariate analysis, we discovered 100 loci with distributed effects across the body composition measures and 241 significant genes primarily involved in immune system functioning. Liver fat stood out, with a highly discoverable and oligogenic architecture and the strongest genetic associations. Comparison with 21 common cardiometabolic traits revealed both shared and specific genetic influences, with higher mean heritability for the MRI measures (h2 = .25 vs. .13, p = 1.8x10−7). We found substantial genetic correlations between the body composition measures and a range of cardiometabolic diseases, with the strongest correlation between liver fat and type 2 diabetes (rg = .49, p = 2.7x10−22). These findings show that MRI-derived body composition measures complement conventional body anthropometrics and other biomarkers of cardiometabolic health, highlighting the central role of liver fat, and improving our knowledge of the genetic architecture of body composition and related diseases.

Genome-wide Association Meta-analysis of Childhood and Adolescent Internalizing Symptoms.

To investigate the genetic architecture of internalizing symptoms in childhood and adolescence. In 22 cohorts, multiple univariate genome-wide association studies (GWASs) were performed using repeated assessments of internalizing symptoms, in a total of 64,561 children and adolescents between 3 and 18 years of age. Results were aggregated in meta-analyses that accounted for sample overlap, first using all available data, and then using subsets of measurements grouped by rater, age, and instrument. The meta-analysis of overall internalizing symptoms (INToverall) detected no genome-wide significant hits and showed low single nucleotide polymorphism (SNP) heritability (1.66%, 95% CI= 0.84-2.48%, neffective= 132,260). Stratified analyses indicated rater-based heterogeneity in genetic effects, with self-reported internalizing symptoms showing the highest heritability (5.63%, 95% CI= 3.08%-8.18%). The contribution of additive genetic effects on internalizing symptoms appeared to be stable over age, with overlapping estimates of SNP heritability from early childhood to adolescence. Genetic correlations were observed with adult anxiety, depression, and the well-being spectrum (|rg| > 0.70), as well as with insomnia, loneliness, attention-deficit/hyperactivity disorder, autism, and childhood aggression (range |rg|= 0.42-0.60), whereas there were no robust associations with schizophrenia, bipolar disorder, obsessive-compulsive disorder, or anorexia nervosa. Genetic correlations indicate that childhood and adolescent internalizing symptoms share substantial genetic vulnerabilities with adult internalizing disorders and other childhood psychiatric traits, which could partially explain both the persistence of internalizing symptoms over time and the high comorbidity among childhood psychiatric traits. Reducing phenotypic heterogeneity in childhood samples will be key in paving the way to future GWAS success.

A pharmacogenomic assessment of psychiatric adverse drug reactions to levetiracetam.

ObjectiveLevetiracetam (LEV) is an effective antiseizure medicine, but 10%–20% of people treated with LEV report psychiatric side‐effects, and up to 1% may have psychotic episodes. Pharmacogenomic predictors of these adverse drug reactions (ADRs) have yet to be identified. We sought to determine the contribution of both common and rare genetic variation to psychiatric and behavioral ADRs associated with LEV.MethodsThis case‐control study compared cases of LEV‐associated behavioral disorder (n = 149) or psychotic reaction (n = 37) to LEV‐exposed people with no history of psychiatric ADRs (n = 920). All samples were of European ancestry. We performed genome‐wide association study (GWAS) analysis comparing those with LEV ADRs to controls. We estimated the polygenic risk scores (PRS) for schizophrenia and compared cases with LEV‐associated psychotic reaction to controls. Rare variant burden analysis was performed using exome sequence data of cases with psychotic reactions (n = 18) and controls (n = 122).ResultsUnivariate GWAS found no significant associations with either LEV‐associated behavioural disorder or LEV‐psychotic reaction. PRS analysis showed that cases of LEV‐associated psychotic reaction had an increased PRS for schizophrenia relative to contr ols (p = .0097, estimate = .4886). The rare‐variant analysis found no evidence of an increased burden of rare genetic variants in people who had experienced LEV‐associated psychotic reaction relative to controls.SignificanceThe polygenic burden for schizophrenia is a risk factor for LEV‐associated psychotic reaction. To assess the clinical utility of PRS as a predictor, it should be tested in an independent and ideally prospective cohort. Larger sample sizes are required for the identification of significant univariate common genetic signals or rare genetic signals associated with psychiatric LEV ADRs.

Multivariate genome-wide association study models to improve prediction of Crohn’s disease risk and identification of potential novel variants

BackgroundCrohn's disease (CD) is a type of inflammatory bowel disease (IBD) that affects the gastrointestinal tract with diverse symptoms. At present, genome-wide association studies (GWAS) has discovered more than 140 genetic loci associated with CD from several datasets. Using the usual univariate GWAS methods, researchers have discovered common variants with small effects. Univariate methods assume independence among the variants that miss subtle combinatorial signals. Multivariate approaches have improved risk prediction and have complemented univariate methods for elucidating the etiology of complex traits and potential novel associations. However, the current multivariate models for CD have been assessed for three datasets (published from 2006 to 2008) under unrelated methodological settings showing a broad performance spectrum. Notably, these multivariate studies do not analyze potential novel variants. Here, we aimed to perform a robust multivariate analysis of a CD dataset different from the one commonly used, and we used the information yielded by the models to identify whether the generated models could provide additional information about the potential novel variants of CD. MethodsTherefore, we compared different multivariate methods and models, LASSO (least absolute shrinkage and selection operator), XGBoost, random forest (RF), Bootstrap stage-wise model selection (BSWiMS), and LDpred, using a strict random subsampling approach to predict the CD risk using a recent GWAS dataset, United Kingdom IBD IBD Genetics Consortium (UKIBDGC), made available in 2017, that had not been used for CD prediction studies.In addition, we assessed the effect of common strategies by increasing and decreasing the number of single-nucleotide polymorphism (SNP) markers (using genotype imputation and linkage disequilibrium (LD)–clumping). ResultsWe found that the LDpred model without any imputation was the best model among all the tested models for predicting the CD risk (area under the receiver operating characteristic curve (AUROC) = 0.667 ± 0.024) in this dataset. We validated the best models using a second dataset (National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) IBD Genetics Consortium, which was previously used in CD prediction studies) in which LDpred was also the best method with a similar performance (AUROC = 0.634 ± 0.009). Based on the importance of the variants yielded by the multivariate models, we identified an unnoticed region within chromosome 6, tagged by SNP rs4945943; this region was close to the gene MARCKS, which appeared to contribute to CD risk. ConclusionsThis research is the first multivariate prediction analysis applied to the UKIBDGC dataset. Our robust multivariate setting analysis enabled us to identify a potential variant that contributed to the CD risk. Multivariate methods are valuable tools for identifying genes that contribute to disease risk.

Shared genetic influences on resting-state functional networks of the brain.

The amplitude of activation in brain resting state networks (RSNs), measured with resting‐state functional magnetic resonance imaging, is heritable and genetically correlated across RSNs, indicating pleiotropy. Recent univariate genome‐wide association studies (GWASs) explored the genetic underpinnings of individual variation in RSN activity. Yet univariate genomic analyses do not describe the pleiotropic nature of RSNs. In this study, we used a novel multivariate method called genomic structural equation modeling to model latent factors that capture the shared genomic influence on RSNs and to identify single nucleotide polymorphisms (SNPs) and genes driving this pleiotropy. Using summary statistics from GWAS of 21 RSNs reported in UK Biobank (N = 31,688), the genomic latent factor analysis was first conducted in a discovery sample (N = 21,081), and then tested in an independent sample from the same cohort (N = 10,607). In the discovery sample, we show that the genetic organization of RSNs can be best explained by two distinct but correlated genetic factors that divide multimodal association networks and sensory networks. Eleven of the 17 factor loadings were replicated in the independent sample. With the multivariate GWAS, we found and replicated nine independent SNPs associated with the joint architecture of RSNs. Further, by combining the discovery and replication samples, we discovered additional SNP and gene associations with the two factors of RSN amplitude. We conclude that modeling the genetic effects on brain function in a multivariate way is a powerful approach to learn more about the biological mechanisms involved in brain function.

Pleiotropy of Alzheimer’s Disease and Educational Attainment: Insights from the Summary Statistics

Epidemiological studies report beneficial associations of higher educational attainment (EDU) with Alzheimer’s disease (AD). Prior genome-wide association studies (GWAS) also reported variants associated with AD and EDU separately. The analysis of pleiotropic predisposition to these phenotypes may shed light on EDU-related protection against AD. We examined pleiotropic predisposition to AD and EDU using Fisher’s method and omnibus test applied to summary statistics for single nucleotide polymorphisms (SNPs) associated with AD and EDU in large-scale univariate GWAS at suggestive-effect (5×10-8

Pleiotropy of Alzheimer’s Disease and Educational Attainment: Insights from the Summary Statistics

Epidemiological studies report beneficial associations of higher educational attainment (EDU) with Alzheimer’s disease (AD). Prior genome-wide association studies (GWAS) also reported variants associated with AD and EDU separately. The analysis of pleiotropic predisposition to these phenotypes may shed light on EDU-related protection against AD. We examined pleiotropic predisposition to AD and EDU using Fisher’s method and omnibus test applied to summary statistics for single nucleotide polymorphisms (SNPs) associated with AD and EDU in large-scale univariate GWAS at suggestive-effect (5×10-8

Pleiotropic predisposition to Alzheimer's disease and educational attainment: insights from the summary statistics analysis.

Epidemiological studies report beneficial associations of higher educational attainment (EDU) with Alzheimer’s disease (AD). Prior genome-wide association studies (GWAS) also reported variants associated with AD and EDU separately. The analysis of pleiotropic associations with these phenotypes may shed light on EDU-related protection against AD. We performed pleiotropic meta-analyses using Fisher’s method and omnibus test applied to summary statistics for single nucleotide polymorphisms (SNPs) associated with AD and EDU in large-scale univariate GWAS at suggestive-effect (5 × 10−8 < p < 0.1) and genome-wide (p ≤ 5 × 10−8) significance levels. We report 53 SNPs that attained p ≤ 5 × 10−8 at least in one of the pleiotropic meta-analyses and were reported in the univariate GWAS at 5 × 10−8 < p < 0.1. Of them, there were 46 pleiotropic SNPs according to Fisher’s method. Additionally, Fisher’s method identified 25 of 206 SNPs with pleiotropic effects, which attained p ≤ 5 × 10−8 in the univariate GWAS. We showed that a large fraction of the pleiotropic associations was affected by a counterintuitive phenomenon of antagonistic genetic heterogeneity, which explains the increase, rather than decrease, of the significance of the pleiotropic associations in the omnibus test. Functional enrichment analysis showed that apart from cancers, gene set harboring the non-pleiotropic SNPs was characterized by late-onset AD and neurodevelopmental disorders. The pleiotropic gene set was characterized by a broad spectrum of progressive neurological and neuromuscular diseases and immune-mediated conditions, including progressive motor neuropathy, multiple sclerosis, Parkinson’s disease, and severe AD. Our results suggest that disentangling genes harboring variants with and without pleiotropic associations with AD and EDU is promising for dissecting heterogeneity in biological mechanisms of AD.Supplementary InformationThe online version contains supplementary material available at 10.1007/s11357-021-00484-1.

Unraveling the Genetic Architecture of Two Complex, Stomata-Related Drought-Responsive Traits by High-Throughput Physiological Phenotyping and GWAS in Cowpea (Vigna. Unguiculata L. Walp).

Drought is one of the most devasting and frequent abiotic stresses in agriculture. While many morphological, biochemical and physiological indicators are being used to quantify plant drought responses, stomatal control, and hence the transpiration and photosynthesis regulation through it, is of particular importance in marking the plant capacity of balancing stress response and yield. Due to the difficulties in simultaneous, large-scale measurement of stomatal traits such as sensitivity and speed of stomatal closure under progressive soil drought, forward genetic mapping of these important behaviors has long been unavailable. The recent emerging phenomic technologies offer solutions to identify the water relations of whole plant and assay the stomatal regulation in a dynamic process at the population level. Here, we report high-throughput physiological phenotyping of water relations of 106 cowpea accessions under progressive drought stress, which, in combination of genome-wide association study (GWAS), enables genetic mapping of the complex, stomata-related drought responsive traits “critical soil water content” (θcri) and “slope of transpiration rate declining” (KTr). The 106 accessions showed large variations in θcri and KTr, indicating that they had broad spectrum of stomatal control in response to soil water deficit, which may confer them different levels of drought tolerance. Univariate GWAS identified six and fourteen significant SNPs associated with θcri and KTr, respectively. The detected SNPs distributed in nine chromosomes and accounted for 8.7–21% of the phenotypic variation, suggesting that both stomatal sensitivity to soil drought and the speed of stomatal closure to completion were controlled by multiple genes with moderate effects. Multivariate GWAS detected ten more significant SNPs in addition to confirming eight of the twenty SNPs as detected by univariate GWAS. Integrated, a final set of 30 significant SNPs associated with stomatal closure were reported. Taken together, our work, by combining phenomics and genetics, enables forward genetic mapping of the genetic architecture of stomatal traits related to drought tolerance, which not only provides a basis for molecular breeding of drought resistant cultivars of cowpea, but offers a new methodology to explore the genetic determinants of water budgeting in crops under stressful conditions in the phenomics era.