Multi-trait Genome-wide Association Study Research Articles

Multi-trait genome-wide association study of opioid addiction: OPRM1 and beyond

Opioid addiction (OA) is moderately heritable, yet only rs1799971, the A118G variant in OPRM1, has been identified as a genome-wide significant association with OA and independently replicated. We applied genomic structural equation modeling to conduct a GWAS of the new Genetics of Opioid Addiction Consortium (GENOA) data together with published studies (Psychiatric Genomics Consortium, Million Veteran Program, and Partners Health), comprising 23,367 cases and effective sample size of 88,114 individuals of European ancestry. Genetic correlations among the various OA phenotypes were uniformly high (rg > 0.9). We observed the strongest evidence to date for OPRM1: lead SNP rs9478500 (p = 2.56 × 10–9). Gene-based analyses identified novel genome-wide significant associations with PPP6C and FURIN. Variants within these loci appear to be pleiotropic for addiction and related traits.

Separation of the effects of two reduced height (Rht) genes and genomic background to select for less Fusarium head blight of short-strawed winter wheat (Triticum aestivum L.) varieties.

FHB resistance shared pleiotropic loci with plant height and anther retention. Genomic prediction allows to select for genomic background reducing FHB susceptibility in the presence of the dwarfing allele Rht-D1b. With the high interest for semi-dwarf cultivars in wheat, finding locally adapted resistance sources against Fusarium head blight (FHB) and FHB-neutral reduced height (Rht) genes is of utmost relevance. In this study, 401 genotypes of European origin without/with dwarfing alleles of Rht-D1 and/or Rht24 were analysed across five environments on FHB severity and the morphological traits such as plant height(PH), anther retention(AR), number of spikelets per ear, ear length and ear density. Data were analysed by combined correlation and path analyses, association mapping and coupling single- and multi-trait genome-wide association studies (ST-GWAS and MT-GWAS, respectively) and genomic prediction (GP). All FHB data were corrected for flowering date or heading stage. High genotypic correlation (rg = 0.74) and direct path effect (0.57) were detected between FHB severity and anther retention (AR). Moderate correlation (rg = -0.55) was found between FHB severity and plant height (PH) with a high indirect path via AR (-0.31). Indirect selection for FHB resistance should concentrate on AR and PH. ST-GWAS identified 25 quantitative trait loci (QTL) for FHB severity, PH and AR, while MT-GWAS detected six QTL across chromosomes 2A, 4D, 5A, 6B and 7B conveying pleiotropic effects on the traits. Rht-D1b was associated with high AR and FHB susceptibility. Our study identified a promising positively acting pleiotropic QTL on chromosome 7B which can be utilized to improve FHB resistance while reducing PH and AR. Rht-D1b genotypes having a high resistance genomic background exhibited lower FHB severity and AR. The use of GP for estimating the genomic background was more effective than selection of GWAS-detected markers. We demonstrated that GP has a great potential and should be exploited by selecting for semi-dwarf winter wheat genotypes with higher FHB resistance due to their genomic background.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce.

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

Multi-trait and cross-population genome-wide association studies across autoimmune and allergic diseases identify shared and distinct genetic component

ObjectivesAutoimmune and allergic diseases are outcomes of the dysregulation of the immune system. Our study aimed to elucidate differences or shared components in genetic backgrounds between autoimmune and allergic diseases.MethodsWe...

Association of Novel Loci With Keratoconus Susceptibility in a Multitrait Genome-Wide Association Study of the UK Biobank Database and Canadian Longitudinal Study on Aging

Keratoconus can be a debilitating corneal ectasia in which the cornea thins, bulges, and steepens into a conical shape. Early features of keratoconus include myopia and irregular astigmatism, which affect vision and can be treated with contact lenses, collagen cross-linking, or, in advanced cases, corneal transplant. Recent estimates of the prevalence of keratoconus based on results of Scheimpflug imaging in young adults are as high as 1.2%. However, obtaining very large keratoconus data sets for a genome-wide association study (GWAS) is problematic because few population studies include Scheimpflug imaging and because severe keratoconus is relatively rare. To identify novel keratoconus loci using corneal resistance factor (CRF) and central corneal thickness (CCT). This multitrait GWAS used European ancestry CRF data from UK Biobank (UKB) (n = 105 427) and the Canadian Longitudinal Study on Aging (CLSA) (n = 18 307) and European ancestry CCT data from the International Glaucoma Genetics Consortium (IGGC) (n = 17 803). The CRF and CCT variants in published keratoconus data sets (4669 cases and 116 547 controls) were compared. The data set from UKB was compiled March 24, 2020; data were released from the CLSA in July 2020; and IGGC data were available from May 1, 2018. Association of CRF and CCT variants with keratoconus risk. The GWAS included 4 cohorts: 105 427 UKB European ancestry (56 134 women [53.2%] and 49 293 men [46.7%]; mean [SD] age, 57 [8] years), 5029 UKB South Asian ancestry (2368 women [47.1%] and 2661 men [52.9%]; mean [SD] age, 54 [8] years), 902 UKB East Asian ancestry (622 women [68.9%] and 280 men [31.0%]; mean [SD] age, 53 [8] years), and 18 307 CLSA European ancestry (9260 women [50.6%] and 9047 men [49.4%]; mean [SD] age, 63 [10] years) participants. A total of 369 CRF and 233 CCT loci were identified, including 36 novel CRF loci and 114 novel CCT loci. Twenty-nine CRF loci and 24 CCT loci were associated with keratoconus. Polygenic risk scores (PRS) were constructed using CRF- and CCT-associated variants and published keratoconus variants. The PRS result showed that adding a CRF- or CCT-based PRS to the keratoconus PRS from previously published variants improved the prediction area under the receiver operating characteristic curve (from 0.705 to 0.756 for CRF and from 0.715 to 0.755 for CCT). These findings support the use of multitrait modeling of corneal parameters in a relatively large data set to identify new keratoconus risk loci and enhance polygenic risk score models.

Genetic Dissection of Heat Stress Tolerance in Faba Bean (Vicia faba L.) Using GWAS.

Heat waves are expected to become more frequent and intense, which will impact faba bean cultivation globally. Conventional breeding methods are effective but take considerable time to achieve breeding goals, and, therefore, the identification of molecular markers associated with key genes controlling heat tolerance can facilitate and accelerate efficient variety development. We phenotyped 134 accessions in six open field experiments during summer seasons at Terbol, Lebanon, at Hudeiba, Sudan, and at Central Ferry, WA, USA from 2015 to 2018. These accessions were genotyped using genotyping by sequencing (GBS), and 10,794 high quality single nucleotide polymorphisms (SNPs) were discovered. These accessions were clustered in one diverse large group, although several discrete groups may exist surrounding it. Fifteen lines belonging to different botanical groups were identified as tolerant to heat. SNPs associated with heat tolerance using single-trait (ST) and multi-trait (MT) genome-wide association studies (GWASs) showed 9 and 11 significant associations, respectively. Through the annotation of the discovered significant SNPs, we found that SNPs from transcription factor helix–loop–helix bHLH143-like S-adenosylmethionine carrier, putative pentatricopeptide repeat-containing protein At5g08310, protein NLP8-like, and photosystem II reaction center PSB28 proteins are associated with heat tolerance.

Genome-Wide Association Analysis Reveals Genetic Loci and Candidate Genes for Chest, Abdominal, and Waist Circumferences in Two Duroc Pig Populations

Chest circumference (CC), abdominal circumference (AC), and waist circumference (WC) are regarded as important indicators for improving economic traits because they can reflect the growth and physiological status in pigs. However, the genetic architecture of CC, AC, and WC is still elusive. Here, we performed single-trait and multi-trait genome-wide association studies (GWASs) for CC, AC, and WC in 2,206 American origin Duroc (AOD) and 2,082 Canadian origin Duroc (COD) pigs. As a result, one novel quantitative trait locus (QTL) on Sus scrofa chromosome (SSC) one was associated with CC and AC in COD pigs, which spans 6.92 Mb (from 170.06 to 176.98 Mb). Moreover, multi-trait GWAS identified 21 significant SNPs associated with the three conformation traits, indicating the multi-trait GWAS is a powerful statistical approach that uncovers pleiotropic locus. Finally, the three candidate genes (ITGA11, TLE3, and GALC) were selected that may play a role in the conformation traits. Further bioinformatics analysis indicated that the candidate genes for the three conformation traits mainly participated in sphingolipid metabolism and lysosome pathways. For all we know, this study was the first GWAS for WC in pigs. In general, our findings further reveal the genetic architecture of CC, AC, and WC, which may offer a useful reference for improving the conformation traits in pigs.

Comparison of Single-Trait and Multi-Trait Genome-Wide Association Models and Inclusion of Correlated Traits in the Dissection of the Genetic Architecture of a Complex Trait in a Breeding Program

Unknown genetic architecture makes it difficult to characterize the genetic basis of traits and associated molecular markers because of the complexity of small effect quantitative trait loci (QTLs), environmental effects, and difficulty in phenotyping. Seedling emergence of wheat (Triticum aestivum L.) from deep planting, has a poorly understood genetic architecture, is a vital factor affecting stand establishment and grain yield, and is historically correlated with coleoptile length. This study aimed to dissect the genetic architecture of seedling emergence while accounting for correlated traits using one multi-trait genome-wide association study (MT-GWAS) model and three single-trait GWAS (ST-GWAS) models. The ST-GWAS models included one single-locus model [mixed-linear model (MLM)] and two multi-locus models [fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK)]. We conducted GWAS using two populations. The first population consisted of 473 varieties from a diverse association mapping panel phenotyped from 2015 to 2019. The second population consisted of 279 breeding lines phenotyped in 2015 in Lind, WA, with 40,368 markers. We also compared the inclusion of coleoptile length and markers associated with reduced height as covariates in our ST-GWAS models. ST-GWAS found 107 significant markers across 19 chromosomes, while MT-GWAS found 82 significant markers across 14 chromosomes. The FarmCPU and BLINK models, including covariates, were able to identify many small effect markers while identifying large effect markers on chromosome 5A. By using multi-locus model breeding, programs can uncover the complex nature of traits to help identify candidate genes and the underlying architecture of a trait, such as seedling emergence.

Abstract MP16: Multi-trait Gwas Of Atherosclerosis Detects Novel Loci And Potential Therapeutic Targets

Atherosclerosis, which is the narrowing of the arterial walls via accumulation of cholesterol-rich arterial plaques, is the leading cause of vascular disease worldwide, including myocardial infarction and ischemic stroke. Although atherosclerosis affects arteries throughout the body, previous genome-wide association studies (GWAS) have been performed on specific atherosclerotic phenotypes such as coronary artery disease (CAD) and peripheral artery disease (PAD). There is substantial evidence to suggest that these more specific atherosclerosis phenotypes share a common genetic etiology. We performed a series of multi-trait GWAS using combinations of two atherosclerosis traits and seven atherosclerosis risk factor traits and detected 31 novel pleiotropic loci. We performed these multi-trait GWAS using the N-GWAMA multi-trait GWAS method and summary statistics for CAD (van der Harst et al. 2018), PAD (Klarin et al. 2019), body mass index (Pulit et al. 2019), type II diabetes (Vujkovic et al. 2020), smoking initiation (Wootton et al. 2020), and lipid traits (Klarin et al. 2019). We identified candidate causal genes for 14 of these loci through colocalization analysis with GTEx expression quantitative trait locus (eQTL) data. VDAC2 and PCSK6 are two candidate causal genes that our results and previous literature suggest are potential therapeutic targets. VDAC2 eQTLs in aorta and tibial artery colocalized with a multi-trait GWAS signal detected in the CAD PAD multi-trait GWAS. Previous work has shown that VDAC2 regulates apoptosis, and our results suggest increased VDAC2 expression in smooth muscle cells could increase smooth muscle cell accumulation in atherosclerotic plaques. A sQTL (splicing QTL) for PCSK6 in liver colocalized with a multi-trait GWAS signal between PAD and LDL. Further analysis of the sQTL signal suggested that the effect allele correlates with a more active isoform of PCSK6 , which could increase lipid fractions and risk of atherosclerosis. These results show that joint analysis of atherosclerotic disease traits and their risk factors allows for identification of unified biology that may offer the opportunity for therapeutic manipulation.

Genetic determinants of blood-cell traits influence susceptibility to childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Despite overlap between genetic risk loci for ALL and hematologic traits, the etiological relevance of dysregulated blood-cell homeostasis remains unclear. We investigated this question in a genome-wide association study (GWAS) of childhood ALL (2,666 affected individuals, 60,272 control individuals) and a multi-trait GWAS of nine blood-cell indices in the UK Biobank. We identified 3,000 blood-cell-trait-associated (p < 5.0× 10-8) variants, explaining 4.0% to 23.9% of trait variation and including 115 loci associated with blood-cell ratios (LMR, lymphocyte-to-monocyte ratio; NLR, neutrophil-to-lymphocyte ratio; PLR, platelet-to-lymphocyte ratio). ALL susceptibility was genetically correlated with lymphocyte counts (rg = 0.088, p = 4.0× 10-4) and PLR (rg = -0.072, p = 0.0017). In Mendelian randomization analyses, genetically predicted increase in lymphocyte counts was associated with increased ALL risk (odds ratio [OR] = 1.16, p = 0.031) and strengthened after accounting for other cell types (OR = 1.43, p = 8.8× 10-4). We observed positive associations with increasing LMR (OR = 1.22, p = 0.0017) and inverse effects for NLR (OR = 0.67, p = 3.1× 10-4) and PLR (OR = 0.80, p = 0.002). Our study shows that a genetically induced shift toward higher lymphocyte counts, overall and in relation to monocytes, neutrophils, and platelets, confers an increased susceptibility to childhood ALL.

Multitrait GWAS to connect disease variants and biological mechanisms.

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

Genetic analyses of gynecological disease identify genetic relationships between uterine fibroids and endometrial cancer, and a novel endometrial cancer genetic risk region at the WNT4 1p36.12 locus.

Endometriosis, polycystic ovary syndrome (PCOS) and uterine fibroids have been proposed as endometrial cancer risk factors; however, disentangling their relationships with endometrial cancer is complicated due to shared risk factors and comorbidities. Using genome-wide association study (GWAS) data, we explored the relationships between these non-cancerous gynecological diseases and endometrial cancer risk by assessing genetic correlation, causal relationships and shared risk loci. We found significant genetic correlation between endometrial cancer and PCOS, and uterine fibroids. Adjustment for genetically predicted body mass index (a risk factor for PCOS, uterine fibroids and endometrial cancer) substantially attenuated the genetic correlation between endometrial cancer and PCOS but did not affect the correlation with uterine fibroids. Mendelian randomization analyses suggested a causal relationship between only uterine fibroids and endometrial cancer. Gene-based analyses revealed risk regions shared between endometrial cancer and endometriosis, and uterine fibroids. Multi-trait GWAS analysis of endometrial cancer and the genetically correlated gynecological diseases identified a novel genome-wide significant endometrial cancer risk locus at 1p36.12, which replicated in an independent endometrial cancer dataset. Interrogation of functional genomic data at 1p36.12 revealed biologically relevant genes, including WNT4 which is necessary for the development of the female reproductive system. In summary, our study provides genetic evidence for a causal relationship between uterine fibroids and endometrial cancer. It further provides evidence that the comorbidity of endometrial cancer, PCOS and uterine fibroids may partly be due to shared genetic architecture. Notably, this shared architecture has revealed a novel genome-wide risk locus for endometrial cancer.

A conditional multi-trait sequence GWAS discovers pleiotropic candidate genes and variants for sheep wool, skin wrinkle and breech cover traits

BackgroundImputation to whole-genome sequence is now possible in large sheep populations. It is therefore of interest to use this data in genome-wide association studies (GWAS) to investigate putative causal variants and genes that underpin economically important traits. Merino wool is globally sought after for luxury fabrics, but some key wool quality attributes are unfavourably correlated with the characteristic skin wrinkle of Merinos. In turn, skin wrinkle is strongly linked to susceptibility to “fly strike” (Cutaneous myiasis), which is a major welfare issue. Here, we use whole-genome sequence data in a multi-trait GWAS to identify pleiotropic putative causal variants and genes associated with changes in key wool traits and skin wrinkle.ResultsA stepwise conditional multi-trait GWAS (CM-GWAS) identified putative causal variants and related genes from 178 independent quantitative trait loci (QTL) of 16 wool and skin wrinkle traits, measured on up to 7218 Merino sheep with 31 million imputed whole-genome sequence (WGS) genotypes. Novel candidate gene findings included the MAT1A gene that encodes an enzyme involved in the sulphur metabolism pathway critical to production of wool proteins, and the ESRP1 gene. We also discovered a significant wrinkle variant upstream of the HAS2 gene, which in dogs is associated with the exaggerated skin folds in the Shar-Pei breed.ConclusionsThe wool and skin wrinkle traits studied here appear to be highly polygenic with many putative candidate variants showing considerable pleiotropy. Our CM-GWAS identified many highly plausible candidate genes for wool traits as well as breech wrinkle and breech area wool cover.

A web-based survey on various symptoms of computer vision syndrome and the genetic understanding based on a multi-trait genome-wide association study

A variety of eye-related symptoms due to the overuse of digital devices is collectively referred to as computer vision syndrome (CVS). In this study, a web-based survey about mind and body functions, including eye strain, was conducted on 1998 Japanese volunteers. To investigate the biological mechanisms behind CVS, a multi-trait genome-wide association study (GWAS), a multivariate analysis on individual-level multivariate data, was performed based on the structural equation modeling methodology assuming a causal pathway for a genetic variant to influence each symptom via a single common latent variable. Twelve loci containing lead variants with a suggestive level of significance were identified. Two loci showed relatively strong signals and were associated with TRABD2B relative to the Wnt signaling pathway and SDK1 having neuronal adhesion and immune functions, respectively. By utilizing publicly available eQTL data, colocalization between GWAS and eQTL signals for four loci was detected, and a locus on 2p25.3 showed a strong colocalization (PPH4 > 0.9) on retinal MYT1L, known to play an important role in neuronal differentiation. This study suggested that the use of multivariate questionnaire data and multi-trait GWAS can lead to biologically reasonable findings and enhance our genetic understanding of complex relationships among symptoms related to CVS.

Genome-Wide Association Analysis of Growth Curve Parameters in Chinese Simmental Beef Cattle.

Simple SummaryComplex traits that require observations over multiple time points for the same individual are called longitudinal traits. Understanding the genetic architecture of beef cattle growth cannot be limited simply to a genome-wide association study (GWAS) for body weight at any specific ages, but should be extended to a more general purpose by considering the longitudinal weight–age using a growth curve approach. We compared three nonlinear models that described the body weight data of Chinese Simmental beef cattle. The parameters of the suitable model were treated as phenotypes of single-trait GWAS and multi-trait GWAS. We identified 87 significant single nucleotide polymorphisms (SNPs) associated with body weight. Many candidate genes associated with body weight were identified which may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.The objective of the present study was to perform a genome-wide association study (GWAS) for growth curve parameters using nonlinear models that fit original weight–age records. In this study, data from 808 Chinese Simmental beef cattle that were weighed at 0, 6, 12, and 18 months of age were used to fit the growth curve. The Gompertz model showed the highest coefficient of determination (R2 = 0.954). The parameters’ mature body weight (A), time-scale parameter (b), and maturity rate (K) were treated as phenotypes for single-trait GWAS and multi-trait GWAS. In total, 9, 49, and 7 significant SNPs associated with A, b, and K were identified by single-trait GWAS; 22 significant single nucleotide polymorphisms (SNPs) were identified by multi-trait GWAS. Among them, we observed several candidate genes, including PLIN3, KCNS3, TMCO1, PRKAG3, ANGPTL2, IGF-1, SHISA9, and STK3, which were previously reported to associate with growth and development. Further research for these candidate genes may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.

Multi-trait GWAS using imputed high-density genotypes from whole-genome sequencing identifies genes associated with body traits in Nile tilapia

BackgroundBody traits are generally controlled by several genes in vertebrates (i.e. polygenes), which in turn make them difficult to identify through association mapping. Increasing the power of association studies by combining approaches such as genotype imputation and multi-trait analysis improves the ability to detect quantitative trait loci associated with polygenic traits, such as body traits.ResultsA multi-trait genome-wide association study (mtGWAS) was performed to identify quantitative trait loci (QTL) and genes associated with body traits in Nile tilapia (Oreochromis niloticus) using genotypes imputed to whole-genome sequences (WGS). To increase the statistical power of mtGWAS for the detection of genetic associations, summary statistics from single-trait genome-wide association studies (stGWAS) for eight different body traits recorded in 1309 animals were used. The mtGWAS increased the statistical power from the original sample size from 13 to 44%, depending on the trait analyzed. The better resolution of the WGS data, combined with the increased power of the mtGWAS approach, allowed the detection of significant markers which were not previously found in the stGWAS. Some of the lead single nucleotide polymorphisms (SNPs) were found within important functional candidate genes previously associated with growth-related traits in other terrestrial species. For instance, we identified SNP within the α1,6-fucosyltransferase (FUT8), solute carrier family 4 member 2 (SLC4A2), A disintegrin and metalloproteinase with thrombospondin motifs 9 (ADAMTS9) and heart development protein with EGF like domains 1 (HEG1) genes, which have been associated with average daily gain in sheep, osteopetrosis in cattle, chest size in goats, and growth and meat quality in sheep, respectively.ConclusionsThe high-resolution mtGWAS presented here allowed the identification of significant SNPs, linked to strong functional candidate genes, associated with body traits in Nile tilapia. These results provide further insights about the genetic variants and genes underlying body trait variation in cichlid fish with high accuracy and strong statistical support.

Association of Genetic Variation With Cirrhosis: A Multi-Trait Genome-Wide Association and Gene–Environment Interaction Study

In contrast to most other common diseases, few genetic variants have been identified that impact risk of cirrhosis. We aimed to identify new genetic variants that predispose to cirrhosis, to test whether such variants, aggregated into a polygenic score, enable genomic risk stratification, and to test whether alcohol intake or body mass index interacts with polygenic predisposition. We conducted a multi-trait genome-wide association study combining cirrhosis and alanine aminotransferase levels performed in 5 discovery studies (UK Biobank, Vanderbilt BioVU, Atherosclerosis Risk in Communities study, and 2 case-control studies including 4829 individuals with cirrhosis and 72,705 controls and 362,539 individuals with alanine aminotransferase levels). Identified variants were replicated in 3 studies (Partners HealthCare Biobank, FinnGen, and Biobank Japan including 3554 individuals with cirrhosis and 343,826 controls). A polygenic score was tested in Partners HealthCare Biobank. Five previously reported and 7 newly identified genetic variants were associated with cirrhosis in both the discovery studies multi-trait genome-wide association study (P < 5× 10-8) and the replication studies (P < .05), including a missense variant in the APOE gene and a noncoding variant near EFN1A. These 12 variants were used to generate a polygenic score. Among Partners HealthCare Biobank individuals, high polygenic score-defined as the top quintile of the distribution-was associated with significantly increased risk of cirrhosis (odds ratio, 2.26; P < .001) and related comorbidities compared with the lowest quintile. Risk was even more pronounced among those with extreme polygenic risk (top 1% of the distribution, odds ratio, 3.16; P < .001). The impact of extreme polygenic risk was substantially more pronounced in those with elevated alcohol consumption or body mass index. Modeled as risk by age 75 years, probability of cirrhosis with extreme polygenic risk was 13.7%, 20.1%, and 48.2% among individuals with no or modest, moderate, and increased alcohol consumption, respectively (Pinteraction < .001). Similarly, probability among those with extreme polygenic risk was 6.5%, 10.3%, and 19.5% among individuals with normal weight, overweight, and obesity, respectively (Pinteraction < .001). Twelve independent genetic variants, 7 of which are newly identified in this study, conferred risk for cirrhosis. Aggregated into a polygenic score, these variants identified a subset of the population at substantially increased risk who are most susceptible to the hepatotoxic effects of excess alcohol consumption or obesity.

A Multiple-Trait Bayesian Variable Selection Regression Method for Integrating Phenotypic Causal Networks in Genome-Wide Association Studies

Bayesian regression methods that incorporate different mixture priors for marker effects are used in multi-trait genomic prediction. These methods can also be extended to genome-wide association studies (GWAS). In multiple-trait GWAS, incorporating the underlying causal structures among traits is essential for comprehensively understanding the relationship between genotypes and traits of interest. Therefore, we develop a GWAS methodology, SEM-Bayesian alphabet, which, by applying the structural equation model (SEM), can be used to incorporate causal structures into multi-trait Bayesian regression methods. SEM-Bayesian alphabet provides a more comprehensive understanding of the genotype-phenotype mapping than multi-trait GWAS by performing GWAS based on indirect, direct and overall marker effects. The superior performance of SEM-Bayesian alphabet was demonstrated by comparing its GWAS results with other similar multi-trait GWAS methods on real and simulated data. The software tool JWAS offers open-source routines to perform these analyses.

Joint association analysis method to dissect complex genetic architecture of multiple genetically related traits

Genome-wide association study (GWAS) has been a standard approach to discover the genetic determinants underlying complex traits. It is a major challenge in GWAS how to improve analysis power, uncover complex genetic correlation, and reveal gene-gene and gene-environment interactions through integrated analysis of multiple genetically related traits. To combat these challenges, we proposed a mixed linear model-based joint association analysis method for multiple traits, which include epistasis and gene-environment interaction in the mapping model and utilize within-trait variance and between-trait covariance simultaneously; A F-statistics based on Wilks statistics is used to test the significance of each SNP and paired interacted SNPs, each genetic effects of QTS are estimated and tested by the MCMC method based on a QTS full model. Simulations showed that the multi-trait GWAS method could provide increased power in detecting pleiotropic loci affecting more than one trait, and can unbiasedly estimate effects of QTS. To demonstrate the performance of the proposed method, we analyzed four blood lipid traits in Multi-Ethnic Study of Atherosclerosis (MESA) Cohort and two yield-related traits in a rice immortalized F2 dataset. A software package was developed for the proposed method.