Multivariate GWAS Research Articles

Genomic structural equation modeling reveals latent phenotypes in the human cortex with distinct genetic architecture

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs, although sensitivity analyses indicated that other structures were plausible. The multivariate GWASs of the GIBNs identified 74 genome-wide significant (GWS) loci (p < 5 × 10−8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed between attention deficit hyperactivity disorder (ADHD) and major depressive disorder (MDD). CT GIBNs displayed a negative genetic correlation with alcohol dependence. Even though we observed model instability in our application of genomic SEM to high-dimensional data, jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across neuroimaging phenotypes offers new insights into the genetics of cortical structure and relationships to psychopathology.

Digital phenotyping from wearables using AI characterizes psychiatric disorders and identifies genetic associations.

Psychiatric disorders are complex and influenced by both genetic and environmental factors. However, studying the full spectrum of these disorders is hindered by practical limitations on measuring human behavior. This highlights the need for novel technologies that can measure behavioral changes at an intermediate level between diagnosis and genotype. Wearable devices are a promising tool in precision medicine, since they can record physiological measurements over time in response to environmental stimuli and do so at low cost and minimal invasiveness. Here we analyzed wearable and genetic data from a cohort of the Adolescent Brain Cognitive Development study. We generated >250 wearable-derived features and used them as intermediate phenotypes in an interpretable AI modeling framework to assign risk scores and classify adolescents with psychiatric disorders. Our model identifies key physiological processes and leverages their temporal patterns to achieve a higher performance than has been previously possible. To investigate how these physiological processes relate to the underlying genetic architecture of psychiatric disorders, we also utilized these intermediate phenotypes in univariate and multivariate GWAS. We identified a total of 29 significant genetic loci and 52 psychiatric-associated genes, including ELFN1 and ADORA3. These results show that wearable-derived continuous features enable a more precise representation of psychiatric disorders and exhibit greater detection power compared to categorical diagnostic labels. In summary, we demonstrate how consumer wearable technology can facilitate dimensional approaches in precision psychiatry and uncover etiological linkages between behavior and genetics.

Maize green leaf area index dynamics: genetic basis of a new secondary trait for grain yield in optimal and drought conditions

Key messageGreen Leaf Area Index dynamics is a promising secondary trait for grain yield and drought tolerance. Multivariate GWAS is particularly well suited to identify the genetic determinants of the green leaf area index dynamics.Improvement of maize grain yield is impeded by important genotype-environment interactions, especially under drought conditions. The use of secondary traits, that are correlated with yield, more heritable and less prone to genotype-environment interactions, can increase breeding efficiency. Here, we studied the genetic basis of a new secondary trait: the green leaf area index (GLAI) dynamics over the maize life cycle. For this, we used an unmanned aerial vehicle to characterize the GLAI dynamics of a diverse panel in well-watered and water-deficient trials in two years. From the dynamics, we derived 24 traits (slopes, durations, areas under the curve), and showed that six of them were heritable traits representative of the panel diversity. To identify the genetic determinants of GLAI, we compared two genome-wide association approaches: a univariate (single-trait) method and a multivariate (multi-trait) method combining GLAI traits, grain yield, and precocity. The explicit modeling of correlation structure between secondary traits and grain yield in the multivariate mixed model led to 2.5 times more associations detected. A total of 475 quantitative trait loci (QTLs) were detected. The genetic architecture of GLAI traits appears less complex than that of yield with stronger-effect QTLs that are more stable between environments. We also showed that a subset of GLAI QTLs explains nearly one fifth of yield variability across a larger environmental network of 11 water-deficient trials. GLAI dynamics is a promising grain yield secondary trait in optimal and drought conditions, and the detected QTLs could help to increase breeding efficiency through a marker-assisted approach.

Multivariate GWAS of Alzheimer’s disease CSF biomarker profiles implies GRIN2D in synaptic functioning

BackgroundGenome-wide association studies (GWAS) of Alzheimer’s disease (AD) have identified several risk loci, but many remain unknown. Cerebrospinal fluid (CSF) biomarkers may aid in gene discovery and we previously demonstrated that six CSF biomarkers (β-amyloid, total/phosphorylated tau, NfL, YKL-40, and neurogranin) cluster into five principal components (PC), each representing statistically independent biological processes. Here, we aimed to (1) identify common genetic variants associated with these CSF profiles, (2) assess the role of associated variants in AD pathophysiology, and (3) explore potential sex differences.MethodsWe performed GWAS for each of the five biomarker PCs in two multi-center studies (EMIF-AD and ADNI). In total, 973 participants (n = 205 controls, n = 546 mild cognitive impairment, n = 222 AD) were analyzed for 7,433,949 common SNPs and 19,511 protein-coding genes. Structural equation models tested whether biomarker PCs mediate genetic risk effects on AD, and stratified and interaction models probed for sex-specific effects.ResultsFive loci showed genome-wide significant association with CSF profiles, two were novel (rs145791381 [inflammation] and GRIN2D [synaptic functioning]) and three were previously described (APOE, TMEM106B, and CHI3L1). Follow-up analyses of the two novel signals in independent datasets only supported the GRIN2D locus, which contains several functionally interesting candidate genes. Mediation tests indicated that variants in APOE are associated with AD status via processes related to amyloid and tau pathology, while markers in TMEM106B and CHI3L1 are associated with AD only via neuronal injury/inflammation. Additionally, seven loci showed sex-specific associations with AD biomarkers.ConclusionsThese results suggest that pathway and sex-specific analyses can improve our understanding of AD genetics and may contribute to precision medicine.

Meta-Analysis and Multivariate GWAS Analyses in 80,950 Individuals of African Ancestry Identify Novel Variants Associated with Blood Pressure Traits.

High blood pressure (HBP) has been implicated as a major risk factor for cardiovascular diseases in several populations, including individuals of African ancestry. Despite the elevated burden of HBP-induced cardiovascular diseases in Africa and other populations of African descent, limited genetic studies have been carried out to explore the genetic mechanism driving this phenomenon. We performed genome-wide association univariate and multivariate analyses of both systolic (SBP) and diastolic blood pressure (DBP) traits in 80,950 individuals of African ancestry. We used summary statistics data from six independent cohorts, including the African Partnership for Chronic Disease Research (APCDR), the UK Biobank, and the Million Veteran Program (MVP). FUMA was used to annotate, prioritize, visualize, and interpret our findings to gain a better understanding of the molecular mechanism(s) underlying the genetics of BP traits. Finally, we undertook a Bayesian fine-mapping analysis to identify potential causal variants. Our meta-analysis identified 10 independent variants associated with SBP and 9 with DBP traits. Whilst our multivariate GWAS method identified 21 independent signals, 18 of these SNPs have been previously identified. SBP was linked to gene sets involved in biological processes such as synapse assembly and cell-cell adhesion via plasma membrane adhesion. Of the 19 independent SNPs identified in the BP meta-analysis, only 11 variants had posterior probability (PP) of > 50%, including one novel variant: rs562545 (MOBP, PP = 77%). To facilitate further research and fine-mapping of high-risk loci/variants in highly susceptible groups for cardiovascular disease and other related traits, large-scale genomic datasets are needed. Our findings highlight the importance of including ancestrally diverse populations in large GWASs and the need for diversity in genetic research.

Genetic dissection of cassava brown streak disease in a genomic selection population.

Cassava brown streak disease (CBSD) is a major threat to food security in East and central Africa. Breeding for resistance against CBSD is the most economical and sustainable way of addressing this challenge. This study seeks to assess the (1) performance of CBSD incidence and severity; (2) identify genomic regions associated with CBSD traits and (3) candidate genes in the regions of interest, in the Cycle 2 population of the National Crops Resources Research Institute. A total of 302 diverse clones were screened, revealing that CBSD incidence across growing seasons was 44%. Severity scores for both foliar and root symptoms ranged from 1.28 to 1.99 and 1.75 to 2.28, respectively across seasons. Broad sense heritability ranged from low to high (0.15 - 0.96), while narrow sense heritability ranged from low to moderate (0.03 - 0.61). Five QTLs, explaining approximately 19% phenotypic variation were identified for CBSD severity at 3 months after planting on chromosomes 1, 13, and 18 in the univariate GWAS analysis. Multivariate GWAS analysis identified 17 QTLs that were consistent with the univariate analysis including additional QTLs on chromosome 6. Seventy-seven genes were identified in these regions with functions such as catalytic activity, ATP-dependent activity, binding, response to stimulus, translation regulator activity, transporter activity among others. These results suggest variation in virulence in the C2 population, largely due to genetics and annotated genes in these QTLs regions may play critical roles in virus initiation and replication, thus increasing susceptibility to CBSD.

Multivariate GWAS of psychiatric disorders and their cardinal symptoms reveal two dimensions of cross-cutting genetic liabilities

SUMMARYUnderstanding which biological pathways are specific versus general across diagnostic categories and levels of symptom severity is critical to improving nosology and treatment of psychopathology. Here, we combine transdiagnostic and dimensional approaches to genetic discovery for the first time, conducting a novel multivariate genome-wide association study of eight psychiatric symptoms and disorders broadly related to mood disturbance and psychosis. We identify two transdiagnostic genetic liabilities that distinguish between common forms of psychopathology versus rarer forms of serious mental illness. Biological annotation revealed divergent genetic architectures that differentially implicated prenatal neurodevelopment and neuronal function and regulation. These findings inform psychiatric nosology and biological models of psychopathology, as they suggest that the severity of mood and psychotic symptoms present in serious mental illness may reflect a difference in kind rather than merely in degree.

Multivariate GWAS of Structural Dental Anomalies and Dental Caries in a Multi-Ethnic Cohort.

Odontogenesis is a complex process, where disruption can result in dental anomalies and/or increase the risk of developing dental caries. Based on previous studies, certain dental anomalies tend to co-occur in patients, suggesting that these traits may share common genetic and etiological components. The main goal of this study was to implement a multivariate genome-wide association study approach to identify genetic variants shared between correlated structural dental anomalies and dental caries. Our cohort (N = 3,579) was derived from the Pittsburgh Orofacial Clefts Study, where multiple dental traits were assessed in both the unaffected relatives of orofacial cleft (OFC) cases (n = 2,187) and unaffected controls (n = 1,392). We identified four multivariate patterns of correlated traits in this data: tooth agenesis, impaction, and rotation (AIR); enamel hypoplasia, displacement, and rotation (HDR); displacement, rotation, and mamelon (DRM); and dental caries, tooth agenesis and enamel hypoplasia (CAH). We analyzed each of these four models using genome-wide multivariate tests of association. No genome-wide statistically significant results were found, but we identified multiple suggestive association signals (P < 10-5) near genes with known biological roles during tooth development, including ADAMTS9 and PRICKLE2 associated with AIR; GLIS3, WDR72, and ROR2 associated with HDR and DRM; ROBO2 associated with DRM; BMP7 associated with HDR; and ROBO1, SMAD2, and MSX2 associated with CAH. This is the first study to investigate genetic associations for multivariate patterns of correlated dental anomalies and dental caries. Further studies are needed to replicate these results in independent cohorts.

YOUR ILLNESS ON DRUGS: THE GENETICS OF SUBSTANCE USE INFORMS OUR RESPONSE TO SOMATIC DISEASES

Overall Abstract: SUDs are among the top contributors to all cause mortality. This session highlights cross-cutting genetic and genomic analyses that illuminate the relationships between substance use disorders (SUDs) and aspects of health and disease beyond the known relationships with other psychiatric disorders. This is a topic that has gotten very little attention in the past. It is timely, given multiple global health crises, as it addresses the role of SUDs in respiratory and cardiometabolic illness, COVID-19 and HIV/AIDS – each among the leading contributors to mortality. The first talk will showcase findings from a phenome-wide analysis (PheWAS) of alcohol and cannabis involvement, using medical center data. It demonstrates that genetic liability for substance use has a profound effect on multiple bodily systems. For example, genetic liability for cannabis use predicted respiratory illness while that for alcohol use predicted heart disease. Individuals with SUDs are at higher risk to develop severe cases of COVID-19. The second talk examines genetic correlations between COVID-19 severity phenotypes and 1014 heritable traits, and reveals that the largest genetic causality proportions were for alcohol use frequency, smoking and cannabis use. Third, we pivot from the current pandemic to the ongoing HIV epidemic to discuss the impact of cocaine abuse among HIV+ individuals. The findings showed cocaine use increased viral load, worsened progression, accelerated mortality, and increase the latent viral reservoir. Specific DNA methylation and gene expression differences that appear to be biological mediators of these relationships were uncovered. The last talk addresses the mediating role of behavioral disinhibition, which is genetically correlated with SUD, using a multivariate GWAS approach to creating a broad genetic measure of externalizing behavior (EXT). Pooling data from ∼1.5 million people, more than 500 genetic loci were discovered. A PheWAS showed that EXT-PRS was associated with 255 of the 1335 disease phenotypes tested (FDR < 0.05). Individuals with higher EXT-PRS showed poorer health in nearly every bodily system. To conclude, the panel will discuss the importance of considering the genetics of SUDs more broadly, to better understand both psychiatric and somatic disorders and the enormous health burden arising from them. Disclosure: Nothing to disclose.

Multivariate genome wide association and network analysis of subcortical imaging phenotypes in Alzheimer\u2019s disease

BackgroundGenome-wide association studies (GWAS) have identified many individual genes associated with brain imaging quantitative traits (QTs) in Alzheimer’s disease (AD). However single marker level association discovery may not be able to address the underlying biological interactions with disease mechanism.ResultsIn this paper, we used the MGAS (Multivariate Gene-based Association test by extended Simes procedure) tool to perform multivariate GWAS on eight AD-relevant subcortical imaging measures. We conducted multiple iPINBPA (integrative Protein-Interaction-Network-Based Pathway Analysis) network analyses on MGAS findings using protein-protein interaction (PPI) data, and identified five Consensus Modules (CMs) from the PPI network. Functional annotation and network analysis were performed on the identified CMs. The MGAS yielded significant hits within APOE, TOMM40 and APOC1 genes, which were known AD risk factors, as well as a few new genes such as LAMA1, XYLB, HSD17B7P2, and NPEPL1. The identified five CMs were enriched by biological processes related to disorders such as Alzheimer’s disease, Legionellosis, Pertussis, and Serotonergic synapse.ConclusionsThe statistical power of coupling MGAS with iPINBPA was higher than traditional GWAS method, and yielded new findings that were missed by GWAS. This study provides novel insights into the molecular mechanism of Alzheimer’s Disease and will be of value to novel gene discovery and functional genomic studies.

An expanded analysis framework for multivariate GWAS connects inflammatory biomarkers to functional variants and disease.

Multivariate methods are known to increase the statistical power to detect associations in the case of shared genetic basis between phenotypes. They have, however, lacked essential analytic tools to follow-up and understand the biology underlying these associations. We developed a novel computational workflow for multivariate GWAS follow-up analyses, including fine-mapping and identification of the subset of traits driving associations (driver traits). Many follow-up tools require univariate regression coefficients which are lacking from multivariate results. Our method overcomes this problem by using Canonical Correlation Analysis to turn each multivariate association into its optimal univariate Linear Combination Phenotype (LCP). This enables an LCP-GWAS, which in turn generates the statistics required for follow-up analyses. We implemented our method on 12 highly correlated inflammatory biomarkers in a Finnish population-based study. Altogether, we identified 11 associations, four of which (F5, ABO, C1orf140 and PDGFRB) were not detected by biomarker-specific analyses. Fine-mapping identified 19 signals within the 11 loci and driver trait analysis determined the traits contributing to the associations. A phenome-wide association study on the 19 representative variants from the signals in 176,899 individuals from the FinnGen study revealed 53 disease associations (p < 1 × 10-4). Several reported pQTLs in the 11 loci provided orthogonal evidence for the biologically relevant functions of the representative variants. Our novel multivariate analysis workflow provides a powerful addition to standard univariate GWAS analyses by enabling multivariate GWAS follow-up and thus promoting the advancement of powerful multivariate methods in genomics.

Genome-wide association studies of ionomic and agronomic traits in USDA mini core collection of rice and comparative analyses of different mapping methods

BackgroundRice is an important human staple food vulnerable to heavy metal contamination leading to serious concerns. High yield with low heavy metal contamination is a common but highly challenging goal for rice breeders worldwide due to lack of genetic knowledge and markers.ResultsTo identify candidate QTLs and develop molecular markers for rice yield and heavy metal content, a total of 191 accessions from the USDA Rice mini-core collection with over 3.2 million SNPs were employed to investigate the QTLs. Sixteen ionomic and thirteen agronomic traits were analyzed utilizing two univariate (GLM and MLM) and two multivariate (MLMM and FarmCPU) GWAS methods. 106, 47, and 97 QTLs were identified for ionomics flooded, ionomics unflooded, and agronomic traits, respectively, with the criterium of p-value < 1.53 × 10− 8, which was determined by the Bonferroni correction for p-value of 0.05. While 49 (~ 20%) of the 250 QTLs were coinciding with previously reported QTLs/genes, about 201 (~ 80%) were new. In addition, several new candidate genes involved in ionomic and agronomic traits control were identified by analyzing the DNA sequence, gene expression, and the homologs of the QTL regions. Our results further showed that each of the four GWAS methods can identify unique as well as common QTLs, suggesting that using multiple GWAS methods can complement each other in QTL identification, especially by combining univariate and multivariate methods.ConclusionsWhile 49 previously reported QTLs/genes were rediscovered, over 200 new QTLs for ionomic and agronomic traits were found in the rice genome. Moreover, multiple new candidate genes for agronomic and ionomic traits were identified. This research provides novel insights into the genetic basis of both ionomic and agronomic variations in rice, establishing the foundation for marker development in breeding and further investigation on reducing heavy-metal contamination and improving crop yields. Finally, the comparative analysis of the GWAS methods showed that each method has unique features and different methods can complement each other.

Combined multivariate factor analysis and GWAS for milk fatty acids trait in Comisana sheep breed.

Combined multivariate factor analysis and GWAS for milk fatty acids trait in Comisana sheep breed.

A Fast Method for Estimating Statistical Power of Multivariate GWAS in Real Case Scenarios: Examples from the Field of Imaging Genetics.

In GWAS of imaging phenotypes (e.g., by the ENIGMA and CHARGE consortia), the growing number of phenotypes considered presents a statistical challenge that other fields are not experiencing (e.g. psychiatry and the Psychiatric Genetics Consortium). However, the multivariate nature of MRI measurements may also be an advantage as many of the MRI phenotypes are correlated and multivariate methods could be considered. Here, we compared the statistical power of a multivariate GWAS versus the current univariate approach, which consists of multiple univariate analyses. To do so, we used results from twin models to estimate pertinent vectors of SNP effect sizes on brain imaging phenotypes, as well as the residual correlation matrices, necessary to estimate analytically the statistical power. We showed that for subcortical structure volumes and hippocampal subfields, a multivariate GWAS yields similar statistical power to the current univariate approach. Our analytical approach is as accurate but ~ 1000 times faster than simulations and we have released the code to facilitate the investigation of other scenarios, may they be outside the field of imaging genetics.

A Brief Critique of the TATES Procedure.

The Trait-based test that uses the Extended Simes procedure (TATES) was developed as a method for conducting multivariate GWAS for correlated phenotypes whose underlying genetic architecture is complex. In this paper, we provide a brief methodological critique of the TATES method using simulated examples and a mathematical proof. Our simulated examples using correlated phenotypes show that the Type I error rate is higher than expected, and that more TATES p values fall outside of the confidence interval relative to expectation. Thus the method may result in systematic inflation when used with correlated phenotypes. In a mathematical proof we further demonstrate that the distribution of TATES p values deviates from expectation in a manner indicative of inflation. Our findings indicate the need for caution when using TATES for multivariate GWAS of correlated phenotypes.

GW-SEM: A Statistical Package to Conduct Genome-Wide Structural Equation Modeling.

Improving the accuracy of phenotyping through the use of advanced psychometric tools will increase the power to find significant associations with genetic variants and expand the range of possible hypotheses that can be tested on a genome-wide scale. Multivariate methods, such as structural equation modeling (SEM), are valuable in the phenotypic analysis of psychiatric and substance use phenotypes, but these methods have not been integrated into standard genome-wide association analyses because fitting a SEM at each single nucleotide polymorphism (SNP) along the genome was hitherto considered to be too computationally demanding. By developing a method that can efficiently fit SEMs, it is possible to expand the set of models that can be tested. This is particularly necessary in psychiatric and behavioral genetics, where the statistical methods are often handicapped by phenotypes with large components of stochastic variance. Due to the enormous amount of data that genome-wide scans produce, the statistical methods used to analyze the data are relatively elementary and do not directly correspond with the rich theoretical development, and lack the potential to test more complex hypotheses about the measurement of, and interaction between, comorbid traits. In this paper, we present a method to test the association of a SNP with multiple phenotypes or a latent construct on a genome-wide basis using a diagonally weighted least squares (DWLS) estimator for four common SEMs: a one-factor model, a one-factor residuals model, a two-factor model, and a latent growth model. We demonstrate that the DWLS parameters and p-values strongly correspond with the more traditional full information maximum likelihood parameters and p-values. We also present the timing of simulations and power analyses and a comparison with and existing multivariate GWAS software package.