Pleiotropy in complex traits: challenges and strategies

Genetic Approaches to Understanding Psychiatric Disease

Human genetic studies have been the driving force in bringing to light the underlying biology of psychiatric conditions. As these studies fill in the gaps in our knowledge of the mechanisms at play, we will be better equipped to design therapies in rational and targeted ways, or repurpose existing therapies in previously unanticipated ways. This review is intended for those unfamiliar with psychiatric genetics as a field and provides a primer on different modes of genetic variation, the technologies currently used to probe them, and concepts that provide context for interpreting the gene-phenotype relationship. Like other subfields in human genetics, psychiatric genetics is moving from microarray technology to sequencing-based approaches as barriers of cost and expertise are removed, and the ramifications of this transition are discussed here. A summary is then given of recent genetic discoveries in a number of neuropsychiatric conditions, with particular emphasis on neurodevelopmental conditions. The general impact of genetics on drug development has been to underscore the extensive etiological heterogeneity in seemingly cohesive diagnostic categories. Consequently, the path forward is not in therapies hoping to reach large swaths of patients sharing a clinically defined diagnosis, but rather in targeting patients belonging to specific "biotypes" defined through a combination of objective, quantifiable data, including genotype.

Genetics: Implications for Prevention and Management of Coronary Artery Disease

An exciting new era has dawned for the prevention and management of coronary artery disease (CAD) utilizing genetic risk variants. The recent identification of over 60 susceptibility loci for CAD confirms not only the importance of established risk factors, but also the existence of many novel causal pathways that are expected to improve our understanding of the genetic basis of CAD and facilitate the development of new therapeutic agents over time. Concurrently, Mendelian randomization studies have provided intriguing insights on the causal relationship between CAD-related traits, and highlight the potential benefits of long-term modifications of risk factors. Last, genetic risk scores of CAD may serve not only as prognostic, but also as predictive markers, and carry the potential to considerably improve the delivery of established prevention strategies. This review will summarize the evolution and discovery of genetic risk variants for CAD and their current and future clinical applications.

Multiple phenotype association tests using summary statistics in genome-wide association studies

We study in this article jointly testing the associations of a genetic variant with correlated multiple phenotypes using the summary statistics of individual phenotype analysis from Genome-Wide Association Studies (GWASs). We estimated the between-phenotype correlation matrix using the summary statistics of individual phenotype GWAS analyses, and developed genetic association tests for multiple phenotypes by accounting for between-phenotype correlation without the need to access individual-level data. Since genetic variants often affect multiple phenotypes differently across the genome and the between-phenotype correlation can be arbitrary, we proposed robust and powerful multiple phenotype testing procedures by jointly testing a common mean and a variance component in linear mixed models for summary statistics. We computed the p-values of the proposed tests analytically. This computational advantage makes our methods practically appealing in large-scale GWASs. We performed simulation studies to show that the proposed tests maintained correct type I error rates, and to compare their powers in various settings with the existing methods. We applied the proposed tests to a GWAS Global Lipids Genetics Consortium summary statistics data set and identified additional genetic variants that were missed by the original single-trait analysis.

Genetic and environmental sources of covariation between early drinking and adult functioning

The vast majority of individuals initiate alcohol consumption for the first time in adolescence. Given the widespread nature of its use and evidence that adolescents may be especially vulnerable to its effects, there is concern about the long-term detrimental impact of adolescent drinking on adult functioning. While some researchers have suggested that genetic processes may confound the relationship, the mechanisms linking drinking and later adjustment remain unclear. The current study utilized a genetically informed sample and biometric modeling to examine the nature of the familial influences on this association and identify the potential for genetic confounding. The sample was drawn from the Minnesota Twin Family Study (MTFS), a longitudinal study consisting of 2,764 twins assessed in 2 cohorts at regular follow-ups from age 17 to age 29 (older cohort) or age 11 to age 29 (younger cohort). A broad range of adult measures was included assessing substance use, antisocial behavior, personality, socioeconomic status, and social functioning. A bivariate Cholesky decomposition was used to examine the common genetic and environmental influences on adolescent drinking and each of the measures of adult adjustment. The results revealed that genetic factors and nonshared environmental influences were generally most important in explaining the relationship between adolescent drinking and later functioning. While the presence of nonshared environmental influences on the association are not inconsistent with a causal impact of adolescent drinking, the findings suggest that many of the adjustment issues associated with adolescent alcohol consumption are best understood as genetically influenced vulnerabilities. (PsycINFO Database Record

Lactose Intolerance (LCT-13910C>T) Genotype Is Associated with Plasma 25-Hydroxyvitamin D Concentrations in Caucasians: A Mendelian Randomization Study

Background: The LCT-13910C>T gene variant is associated with lactose intolerance (LI) in different ethnic groups. Individuals with LI often limit or avoid dairy consumption, a major dietary source of vitamin D in North America, which may lead to inadequate vitamin D intake.Objective: The objective was to determine the prevalence of genotypes predictive of LI in different ethnic groups living in Canada and to determine whether the LCT genotype is associated with plasma 25(OH)D concentrations.Methods: Blood samples were drawn from a total of 1495 men and women aged 20-29 y from the Toronto Nutrigenomics and Health Study for genotyping and plasma 25(OH)D analysis. Intakes of dairy were assessed by using a 196-item food frequency questionnaire. The prevalence of LCT-13910C>T genotypes was compared by using χ² analysis. Using a Mendelian randomization approach, we examined the association between LCT genotypes and 25(OH)D concentrations.Results: Approximately 32% of Caucasians, 99% of East Asians, 74% of South Asians, and 59% of those with other or mixed ethnicities had the CC genotype associated with LI. Compared with those with the TT genotype, those with the CC genotype had a lower mean ± SE total dairy intake (2.15 ± 0.09 compared with 2.67 ± 0.12 servings/d, P = 0.003), a lower skim-milk intake (0.20 ± 0.03 compared with 0.46 ± 0.06 servings/d, P = 0.0004), and a lower plasma 25(OH)D concentration (63 ± 1.9 compared with 75.8 ± 2.4 nmol/L, P < 0.0001). The CT and CC genotypes were associated with a 50% and a 2-fold increased risk, respectively, of a suboptimal plasma 25(OH)D concentration (<75 nmol/L).Conclusions: In Caucasians, the CC genotype that predicts LI is associated with a lower plasma 25(OH)D concentration, which is attributable at least in part to a lower intake of dairy, particularly skim milk. Increased risk of suboptimal concentrations of vitamin D was also observed among those with the CT genotype, suggesting an intermediate effect of the heterozygous genotype.

Genetics of common complex diseases: a view from Iceland

In the past decade, large scale genotyping has led to discoveries of numerous sequence variants that confer increased risk of many common complex diseases. Interestingly, a substantial proportion of pioneering genetic work has originated from the small nation of Iceland and has been facilitated by an extensive genealogy database. We provide examples of relevant observations made so far in several major disease categories central to internal medicine practice. Some of these findings offer new mechanistic clues into the pathophysiology of common disorders and may suggest novel approaches in diagnosis and drug therapy. However, a number of unresolved issues remain that will be subject of future research, driven by recent advances in high-throughput sequencing of the genome. At the same time, we are ready to begin transforming the abundant existing genetic data into practical clinical knowledge with the aim of improving the delivery of medical care. The era of precision medicine has arrived.

Genetic evidence for role of integration of fast and slow neurotransmission in schizophrenia

The most recent genome-wide association studies (GWAS) of schizophrenia (SCZ) identified hundreds of risk variants potentially implicated in the disease. Further, novel statistical methodology designed for polygenic architecture revealed more potential risk variants. This can provide a link between individual genetic factors and the mechanistic underpinnings of SCZ. Intriguingly, a large number of genes coding for ionotropic and metabotropic receptors for various neurotransmitters-glutamate, γ-aminobutyric acid (GABA), dopamine, serotonin, acetylcholine and opioids-and numerous ion channels were associated with SCZ. Here, we review these findings from the standpoint of classical neurobiological knowledge of neuronal synaptic transmission and regulation of electrical excitability. We show that a substantial proportion of the identified genes are involved in intracellular cascades known to integrate 'slow' (G-protein-coupled receptors) and 'fast' (ionotropic receptors) neurotransmission converging on the protein DARPP-32. Inspection of the Human Brain Transcriptome Project database confirms that that these genes are indeed expressed in the brain, with the expression profile following specific developmental trajectories, underscoring their relevance to brain organization and function. These findings extend the existing pathophysiology hypothesis by suggesting a unifying role of dysregulation in neuronal excitability and synaptic integration in SCZ. This emergent model supports the concept of SCZ as an 'associative' disorder-a breakdown in the communication across different slow and fast neurotransmitter systems through intracellular signaling pathways-and may unify a number of currently competing hypotheses of SCZ pathophysiology.

Identification of genetic loci shared between schizophrenia and the Big Five personality traits

Schizophrenia is associated with differences in personality traits, and recent studies suggest that personality traits and schizophrenia share a genetic basis. Here we aimed to identify specific genetic loci shared between schizophrenia and the Big Five personality traits using a Bayesian statistical framework. Using summary statistics from genome-wide association studies (GWAS) on personality traits in the 23andMe cohort (n = 59,225) and schizophrenia in the Psychiatric Genomics Consortium cohort (n = 82,315), we evaluated overlap in common genetic variants. The Big Five personality traits neuroticism, extraversion, openness, agreeableness and conscientiousness were measured using a web implementation of the Big Five Inventory. Applying the conditional false discovery rate approach, we increased discovery of genetic loci and identified two loci shared between neuroticism and schizophrenia and six loci shared between openness and schizophrenia. The study provides new insights into the relationship between personality traits and schizophrenia by highlighting genetic loci involved in their common genetic etiology.

Phenome-wide scanning identifies multiple diseases and disease severity phenotypes associated with HLA variants

Although many phenotypes have been associated with variants in human leukocyte antigen (HLA) genes, the full phenotypic impact of HLA variants across all diseases is unknown. We imputed HLA genomic variation from two populations of 28,839 and 8431 European ancestry individuals and tested association of HLA variation with 1368 phenotypes. A total of 104 four-digit and 92 two-digit HLA allele phenotype associations were significant in both discovery and replication cohorts, the strongest being HLA-DQB1*03:02 and type 1 diabetes. Four previously unidentified associations were identified across the spectrum of disease with two- and four-digit HLA alleles and 10 with nonsynonymous variants. Some conditions associated with multiple HLA variants and stronger associations with more severe disease manifestations were identified. A comprehensive, publicly available catalog of clinical phenotypes associated with HLA variation is provided. Examining HLA variant disease associations in this large data set allows comprehensive definition of disease associations to drive further mechanistic insights.

Genetic architecture distinguishes systemic juvenile idiopathic arthritis from other forms of juvenile idiopathic arthritis: clinical and therapeutic implications

OBJECTIVES

Juvenile idiopathic arthritis (JIA) is a heterogeneous group of conditions unified by the presence of chronic childhood arthritis without an identifiable cause. Systemic JIA (sJIA) is a rare form of JIA characterised by systemic inflammation. sJIA is distinguished from other forms of JIA by unique clinical features and treatment responses that are similar to autoinflammatory diseases. However, approximately half of children with sJIA develop destructive, long-standing arthritis that appears similar to other forms of JIA. Using genomic approaches, we sought to gain novel insights into the pathophysiology of sJIA and its relationship with other forms of JIA.

METHODS

We performed a genome-wide association study of 770 children with sJIA collected in nine countries by the International Childhood Arthritis Genetics Consortium. Single nucleotide polymorphisms were tested for association with sJIA. Weighted genetic risk scores were used to compare the genetic architecture of sJIA with other JIA subtypes.

RESULTS

The major histocompatibility complex locus and a locus on chromosome 1 each showed association with sJIA exceeding the threshold for genome-wide significance, while 23 other novel loci were suggestive of association with sJIA. Using a combination of genetic and statistical approaches, we found no evidence of shared genetic architecture between sJIA and other common JIA subtypes.

CONCLUSIONS

The lack of shared genetic risk factors between sJIA and other JIA subtypes supports the hypothesis that sJIA is a unique disease process and argues for a different classification framework. Research to improve sJIA therapy should target its unique genetics and specific pathophysiological pathways.

An Analytic Solution to the Computation of Power and Sample Size for Genetic Association Studies under a Pleiotropic Mode of Inheritance

Our motivation here is to calculate the power of 3 statistical tests used when there are genetic traits that operate under a pleiotropic mode of inheritance and when qualitative phenotypes are defined by use of thresholds for the multiple quantitative phenotypes. Specifically, we formulate a multivariate function that provides the probability that an individual has a vector of specific quantitative trait values conditional on having a risk locus genotype, and we apply thresholds to define qualitative phenotypes (affected, unaffected) and compute penetrances and conditional genotype frequencies based on the multivariate function. We extend the analytic power and minimum-sample-size-necessary (MSSN) formulas for 2 categorical data-based tests (genotype, linear trend test [LTT]) of genetic association to the pleiotropic model. We further compare the MSSN of the genotype test and the LTT with that of a multivariate ANOVA (Pillai). We approximate the MSSN for statistics by linear models using a factorial design and ANOVA. With ANOVA decomposition, we determine which factors most significantly change the power/MSSN for all statistics. Finally, we determine which test statistics have the smallest MSSN. In this work, MSSN calculations are for 2 traits (bivariate distributions) only (for illustrative purposes). We note that the calculations may be extended to address any number of traits. Our key findings are that the genotype test usually has lower MSSN requirements than the LTT. More inclusive thresholds (top/bottom 25% vs. top/bottom 10%) have higher sample size requirements. The Pillai test has a much larger MSSN than both the genotype test and the LTT, as a result of sample selection. With these formulas, researchers can specify how many subjects they must collect to localize genes for pleiotropic phenotypes.

Efficient inference for genetic association studies with multiple outcomes

Combined inference for heterogeneous high-dimensional data is critical in modern biology, where clinical and various kinds of molecular data may be available from a single study. Classical genetic association studies regress a single clinical outcome on many genetic variants one by one, but there is an increasing demand for joint analysis of many molecular outcomes and genetic variants in order to unravel functional interactions. Unfortunately, most existing approaches to joint modeling are either too simplistic to be powerful or are impracticable for computational reasons. Inspired by Richardson and others (2010, Bayesian Statistics 9), we consider a sparse multivariate regression model that allows simultaneous selection of predictors and associated responses. As Markov chain Monte Carlo (MCMC) inference on such models can be prohibitively slow when the number of genetic variants exceeds a few thousand, we propose a variational inference approach which produces posterior information very close to that of MCMC inference, at a much reduced computational cost. Extensive numerical experiments show that our approach outperforms popular variable selection methods and tailored Bayesian procedures, dealing within hours with problems involving hundreds of thousands of genetic variants and tens to hundreds of clinical or molecular outcomes.

Multivariate simulation framework reveals performance of multi-trait GWAS methods

Burgeoning availability of genome-wide association study (GWAS) results and national biobank data has led to growing interest in performing multi-trait genetic analyses. Numerous multi-trait GWAS methods that exploit either summary statistics or individual-level data have been developed, but their relative performance is unclear. Here we develop a simulation framework to model the complex networks underlying multivariate genetic epidemiology, enabling the vast model space of genetic effects on multiple correlated traits to be explored systematically. We perform a comprehensive comparison of the leading multi-trait GWAS methods, finding: (1) method performance is highly sensitive to the specific combination of genetic effects and phenotypic correlations, (2) most of the current multivariate methods have remarkably similar statistical power, and (3) multivariate methods may offer a substantial increase in the discovery of genetic variants over the standard univariate approach. We believe our findings offer the clearest picture to date of the relative performance of multi-trait GWAS methods and act as a guide for method selection. We provide a web application and open-source software program implementing our simulation framework, for: (i) further benchmarking of multivariate GWAS methods, (ii) power calculations for multivariate genetic studies, and (iii) generating data for testing any multivariate method in genetic epidemiology.

Shared regulatory sites are abundant in the human genome and shed light on genome evolution and disease pleiotropy

Large-scale gene expression datasets are providing an increasing understanding of the location of cis-eQTLs in the human genome and their role in disease. However, little is currently known regarding the extent of regulatory site-sharing between genes. This is despite it having potentially wide-ranging implications, from the determination of the way in which genetic variants may shape multiple phenotypes to the understanding of the evolution of human gene order. By first identifying the location of non-redundant cis-eQTLs, we show that regulatory site-sharing is a relatively common phenomenon in the human genome, with over 10% of non-redundant regulatory variants linked to the expression of multiple nearby genes. We show that these shared, local regulatory sites are linked to high levels of chromatin looping between the regulatory sites and their associated genes. In addition, these co-regulated gene modules are found to be strongly conserved across mammalian species, suggesting that shared regulatory sites have played an important role in shaping human gene order. The association of these shared cis-eQTLs with multiple genes means they also appear to be unusually important in understanding the genetics of human phenotypes and pleiotropy, with shared regulatory sites more often linked to multiple human phenotypes than other regulatory variants. This study shows that regulatory site-sharing is likely an underappreciated aspect of gene regulation and has important implications for the understanding of various biological phenomena, including how the two and three dimensional structures of the genome have been shaped and the potential causes of disease pleiotropy outside coding regions.

Drosophila and genome-wide association studies: a review and resource for the functional dissection of human complex traits

Human genome-wide association studies (GWAS) have successfully identified thousands of susceptibility loci for common diseases with complex genetic etiologies. Although the susceptibility variants identified by GWAS usually have only modest effects on individual disease risk, they contribute to a substantial burden of trait variation in the overall population. GWAS also offer valuable clues to disease mechanisms that have long proven to be elusive. These insights could lead the way to breakthrough treatments; however, several challenges hinder progress, making innovative approaches to accelerate the follow-up of results from GWAS an urgent priority. Here, we discuss the largely untapped potential of the fruit fly, Drosophila melanogaster, for functional investigation of findings from human GWAS. We highlight selected examples where strong genomic conservation with humans along with the rapid and powerful genetic tools available for flies have already facilitated fine mapping of association signals, elucidated gene mechanisms, and revealed novel disease-relevant biology. We emphasize current research opportunities in this rapidly advancing field, and present bioinformatic analyses that systematically explore the applicability of Drosophila for interrogation of susceptibility signals implicated in more than 1000 human traits, based on all GWAS completed to date. Thus, our discussion is targeted at both human geneticists seeking innovative strategies for experimental validation of findings from GWAS, as well as the Drosophila research community, by whom ongoing investigations of the implicated genes will powerfully inform our understanding of human disease.

A systematic SNP selection approach to identify mechanisms underlying disease aetiology: linking height to post-menopausal breast and colorectal cancer risk

Data from GWAS suggest that SNPs associated with complex diseases or traits tend to co-segregate in regions of low recombination, harbouring functionally linked gene clusters. This phenomenon allows for selecting a limited number of SNPs from GWAS repositories for large-scale studies investigating shared mechanisms between diseases. For example, we were interested in shared mechanisms between adult-attained height and post-menopausal breast cancer (BC) and colorectal cancer (CRC) risk, because height is a risk factor for these cancers, though likely not a causal factor. Using SNPs from public GWAS repositories at p-values < 1 × 10⁻⁵ and a genomic sliding window of 1 mega base pair, we identified SNP clusters including at least one SNP associated with height and one SNP associated with either post-menopausal BC or CRC risk (or both). SNPs were annotated to genes using HapMap and GRAIL and analysed for significantly overrepresented pathways using ConsensuspathDB. Twelve clusters including 56 SNPs annotated to 26 genes were prioritised because these included at least one height- and one BC risk- or CRC risk-associated SNP annotated to the same gene. Annotated genes were involved in Indian hedgehog signalling (p-value = 7.78 × 10⁻⁷) and several cancer site-specific pathways. This systematic approach identified a limited number of clustered SNPs, which pinpoint potential shared mechanisms linking together the complex phenotypes height, post-menopausal BC and CRC.

Identifying Pleiotropic Genes in Genome-Wide Association Studies for Multivariate Phenotypes with Mixed Measurement Scales

We propose a multivariate genome-wide association test for mixed continuous, binary, and ordinal phenotypes. A latent response model is used to estimate the correlation between phenotypes with different measurement scales so that the empirical distribution of the Fisher’s combination statistic under the null hypothesis is estimated efficiently. The simulation study shows that our proposed correlation estimation methods have high levels of accuracy. More importantly, our approach conservatively estimates the variance of the test statistic so that the type I error rate is controlled. The simulation also shows that the proposed test maintains the power at the level very close to that of the ideal analysis based on known latent phenotypes while controlling the type I error. In contrast, conventional approaches–dichotomizing all observed phenotypes or treating them as continuous variables–could either reduce the power or employ a linear regression model unfit for the data. Furthermore, the statistical analysis on the database of the Study of Addiction: Genetics and Environment (SAGE) demonstrates that conducting a multivariate test on multiple phenotypes can increase the power of identifying markers that may not be, otherwise, chosen using marginal tests. The proposed method also offers a new approach to analyzing the Fagerström Test for Nicotine Dependence as multivariate phenotypes in genome-wide association studies.

Landscape of Pleiotropic Proteins Causing Human Disease: Structural and System Biology Insights

Pleiotropy is the phenomenon by which the same gene can result in multiple phenotypes. Pleiotropic proteins are emerging as important contributors to rare and common disorders. Nevertheless, little is known on the mechanisms underlying pleiotropy and the characteristic of pleiotropic proteins. We analyzed disease-causing proteins reported in UniProt and observed that 12% are pleiotropic (variants in the same protein cause more than one disease). Pleiotropic proteins were enriched in deleterious and rare variants, but not in common variants. Pleiotropic proteins were more likely to be involved in the pathogenesis of neoplasms, neurological, and circulatory diseases and congenital malformations, whereas non-pleiotropic proteins in endocrine and metabolic disorders. Pleiotropic proteins were more essential and had a higher number of interacting partners compared with non-pleiotropic proteins. Significantly more pleiotropic than non-pleiotropic proteins contained at least one intrinsically long disordered region (P < 0.001). Deleterious variants occurring in structurally disordered regions were more commonly found in pleiotropic, rather than non-pleiotropic proteins. In conclusion, pleiotropic proteins are an important contributor to human disease. They represent a biologically different class of proteins compared with non-pleiotropic proteins and a better understanding of their characteristics and genetic variants can greatly aid in the interpretation of genetic studies and drug design.

SCOPA and META-SCOPA: software for the analysis and aggregation of genome-wide association studies of multiple correlated phenotypes

Background

Genome-wide association studies (GWAS) of single nucleotide polymorphisms (SNPs) have been successful in identifying loci contributing genetic effects to a wide range of complex human diseases and quantitative traits. The traditional approach to GWAS analysis is to consider each phenotype separately, despite the fact that many diseases and quantitative traits are correlated with each other, and often measured in the same sample of individuals. Multivariate analyses of correlated phenotypes have been demonstrated, by simulation, to increase power to detect association with SNPs, and thus may enable improved detection of novel loci contributing to diseases and quantitative traits.

Results

We have developed the SCOPA software to enable GWAS analysis of multiple correlated phenotypes. The software implements “reverse regression” methodology, which treats the genotype of an individual at a SNP as the outcome and the phenotypes as predictors in a general linear model. SCOPA can be applied to quantitative traits and categorical phenotypes, and can accommodate imputed genotypes under a dosage model. The accompanying META-SCOPA software enables meta-analysis of association summary statistics from SCOPA across GWAS. Application of SCOPA to two GWAS of high-and low-density lipoprotein cholesterol, triglycerides and body mass index, and subsequent meta-analysis with META-SCOPA, highlighted stronger association signals than univariate phenotype analysis at established lipid and obesity loci. The META-SCOPA meta-analysis also revealed a novel signal of association at genome-wide significance for triglycerides mapping to GPC5 (lead SNP rs71427535, p = 1.1x10⁻⁸), which has not been reported in previous large-scale GWAS of lipid traits.

Conclusions

The SCOPA and META-SCOPA software enable discovery and dissection of multiple phenotype association signals through implementation of a powerful reverse regression approach.

Genome-wide association test of multiple continuous traits using imputed SNPs

More and more large cohort studies have conducted or are conducting genome-wide association studies (GWAS) to reveal the genetic components of many complex human diseases. These large cohort studies often collected a broad array of correlated phenotypes that reflect common physiological processes. By jointly analyzing these correlated traits, we can gain more power by aggregating multiple weak effects and shed light on the mechanisms underlying complex human diseases. The majority of existing multi-trait association test methods are based on jointly modeling the multivariate traits conditional on the genotype as covariate, and can readily accommodate the imputed SNPs by using their imputed dosage as a covariate. An alternative class of multi-trait association tests is based on the inverted regression, which models the distribution of genotypes conditional on the covariate and multivariate traits, and has been shown to have competitive performance. To our knowledge, all existing inverted regression approaches have implicitly used the "best-guess" genotypes, which is not efficient and known to lead to dramatic power loss, and there have not been any proposed methods of incorporating imputation uncertainty into inverted regressions. In this work, we propose a general and efficient framework that can account for the imputation uncertainty to further improve the association test power of inverted regression models for imputed SNPs. We demonstrate through extensive numerical studies that the proposed method has competitive performance. We further illustrate its usefulness by application to association test of diabetes-related glycemic traits in the Atherosclerosis Risk in Communities (ARIC) Study.

Dissection of Insertion-Deletion Variants within Differentially Expressed Genes Involved in Wood Formation in Populus

Short insertions and deletions (InDels) are one of the major genetic variants and are distributed widely across the genome; however, few investigations of InDels have been conducted in long-lived perennial plants. Here, we employed a combination of RNA-seq and population resequencing to identify InDels within differentially expressed (DE) genes underlying wood formation in a natural population of Populus tomentosa (435 individuals) and utilized InDel-based association mapping to detect the causal variants under additive, dominance, and epistasis underlying growth and wood properties. In the present paper, 5,482 InDels detected from 629 DE genes showed uneven distributions throughout all 19 chromosomes, and 95.9% of these loci were diallelic InDels. Seventy-four InDels (positive false discovery rate q ≤ 0.10) from 68 genes exhibited significant additive/dominant effects on 10 growth and wood-properties, with an average of 14.7% phenotypic variance explained. Potential pleiotropy was observed in one-third of the InDels (representing 24 genes). Seven genes exhibited significantly differential expression among the genotypic classes of associated InDels, indicating possible important roles for these InDels. Epistasis analysis showed that overlapping interacting genes formed unique interconnected networks for each trait, supporting the putative biochemical links that control quantitative traits. Therefore, the identification and utilization of InDels in trees will be recognized as an effective marker system for molecular marker-assisted breeding applications, and further facilitate our understanding of quantitative genomics.

Cross-Phenotype Association Analysis Using Summary Statistics from GWAS

For over a decade, genome-wide association studies (GWAS) have been a major tool for detecting genetic variants underlying complex traits. Recent studies have demonstrated that the same variant or gene can be associated with multiple traits, and such associations are termed cross-phenotype (CP) associations. CP association analysis can improve statistical power by searching for variants that contribute to multiple traits, which is often relevant to pleiotropy. In this chapter, we discuss existing statistical methods for analyzing association between a single marker and multivariate phenotypes, we introduce a general approach, CPASSOC, to detect the CP associations, and explain how to conduct the analysis in practice.

More than Meets the Eye: A Primer for "Timing of Locomotor Recovery from Anoxia Modulated by the white Gene in Drosophila melanogaster"

SummaryA single gene might have several functions within an organism, and so mutational loss of that gene has multiple effects across different physiological systems in the organism. Though the white gene in Drosophila melanogaster was identified originally for its effect on fly eye color, an article by Xiao and Robertson in the June 2016 issue of GENETICS describes a function for the white gene in the response of Drosophila to oxygen deprivation. This Primer article provides background information on the white gene, the phenomenon of pleiotropy, and the molecular and genetic approaches used in the study to demonstrate a new behavioral function for the white gene.

A comparison study of multivariate fixed models and Gene Association with Multiple Traits (GAMuT) for next-generation sequencing

In this paper, extensive simulations are performed to compare two statistical methods to analyze multiple correlated quantitative phenotypes: (1) approximate F-distributed tests of multivariate functional linear models (MFLM) and additive models of multivariate analysis of variance (MANOVA), and (2) Gene Association with Multiple Traits (GAMuT) for association testing of high-dimensional genotype data. It is shown that approximate F-distributed tests of MFLM and MANOVA have higher power and are more appropriate for major gene association analysis (i.e., scenarios in which some genetic variants have relatively large effects on the phenotypes); GAMuT has higher power and is more appropriate for analyzing polygenic effects (i.e., effects from a large number of genetic variants each of which contributes a small amount to the phenotypes). MFLM and MANOVA are very flexible and can be used to perform association analysis for (i) rare variants, (ii) common variants, and (iii) a combination of rare and common variants. Although GAMuT was designed to analyze rare variants, it can be applied to analyze a combination of rare and common variants and it performs well when (1) the number of genetic variants is large and (2) each variant contributes a small amount to the phenotypes (i.e., polygenes). MFLM and MANOVA are fixed effect models that perform well for major gene association analysis. GAMuT can be viewed as an extension of sequence kernel association tests (SKAT). Both GAMuT and SKAT are more appropriate for analyzing polygenic effects and they perform well not only in the rare variant case, but also in the case of a combination of rare and common variants. Data analyses of European cohorts and the Trinity Students Study are presented to compare the performance of the two methods.

Epigenetic alterations in blood mirror age-associated DNA methylation and gene expression changes in human liver

AIM

To study the impact of aging on DNA methylation and mRNA expression in human liver.

EXPERIMENTAL PROCEDURES

We analysed genome-wide DNA methylation and gene expression in human liver samples using Illumina 450K and HumanHT12 expression BeadChip arrays.

RESULTS

DNA methylation analysis of ∼455,000 CpG sites in human liver revealed that age was significantly associated with altered DNA methylation of 20,396 CpG sites. Comparison of liver methylation data with published methylation data in other tissues showed that vast majority of the age-associated significant CpG sites overlapped between liver and blood, whereas a smaller overlap was found between liver and pancreatic islets or adipose tissue, respectively. We identified 151 genes whose liver expression also correlated with age.

CONCLUSIONS

We identified age-associated DNA methylation and expression changes in human liver that are partly reflected by epigenetic alterations in blood.

Genetic variants linked to education predict longevity

Educational attainment is associated with many health outcomes, including longevity. It is also known to be substantially heritable. Here, we used data from three large genetic epidemiology cohort studies (Generation Scotland, n = ∼17,000; UK Biobank, n = ∼115,000; and the Estonian Biobank, n = ∼6,000) to test whether education-linked genetic variants can predict lifespan length. We did so by using cohort members' polygenic profile score for education to predict their parents' longevity. Across the three cohorts, meta-analysis showed that a 1 SD higher polygenic education score was associated with ∼2.7% lower mortality risk for both mothers (total ndeaths = 79,702) and ∼2.4% lower risk for fathers (total ndeaths = 97,630). On average, the parents of offspring in the upper third of the polygenic score distribution lived 0.55 y longer compared with those of offspring in the lower third. Overall, these results indicate that the genetic contributions to educational attainment are useful in the prediction of human longevity.

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

Biomechanisms of Comorbidity: Reviewing Integrative Analyses of Multi-omics Datasets and Electronic Health Records

OBJECTIVES

Disease comorbidity is a pervasive phenomenon impacting patients' health outcomes, disease management, and clinical decisions. This review presents past, current and future research directions leveraging both phenotypic and molecular information to uncover disease similarity underpinning the biology and etiology of disease comorbidity.

METHODS

We retrieved ~130 publications and retained 59, ranging from 2006 to 2015, that comprise a minimum number of five diseases and at least one type of biomolecule. We surveyed their methods, disease similarity metrics, and calculation of comorbidities in the electronic health records, if present.

RESULTS

Among the surveyed studies, 44% generated or validated disease similarity metrics in context of comorbidity, with 60% being published in the last two years. As inputs, 87% of studies utilized intragenic loci and proteins while 13% employed RNA (mRNA, LncRNA or miRNA). Network modeling was predominantly used (35%) followed by statistics (28%) to impute similarity between these biomolecules and diseases. Studies with large numbers of biomolecules and diseases used network models or naïve overlap of disease-molecule associations, while machine learning, statistics, and information retrieval were utilized in smaller and moderate sized studies. Multiscale computations comprising shared function, network topology, and phenotypes were performed exclusively on proteins.

CONCLUSION

This review highlighted the growing methods for identifying the molecular mechanisms underpinning comorbidities that leverage multiscale molecular information and patterns from electronic health records. The survey unveiled that intergenic polymorphisms have been overlooked for similarity imputation compared to their intragenic counterparts, offering new opportunities to bridge the mechanistic and similarity gaps of comorbidity.

Comparison of Heritability Estimation and Linkage Analysis for Multiple Traits Using Principal Component Analyses

A disease trait often can be characterized by multiple phenotypic measurements that can provide complementary information on disease etiology, physiology, or clinical manifestations. Given that multiple phenotypes may be correlated and reflect common underlying genetic mechanisms, the use of multivariate analysis of multiple traits may improve statistical power to detect genes and variants underlying complex traits. The literature, however, has been unclear as to the optimal approach for analyzing multiple correlated traits. In this study, heritability and linkage analysis was performed for six obstructive sleep apnea hypopnea syndrome (OSAHS) related phenotypes, as well as principal components of the phenotypes and principal components of the heritability (PCHs) using the data from Cleveland Family Study, which include both African and European American families. Our study demonstrates that principal components generally result in higher heritability and linkage evidence than individual traits. Furthermore, the PCHs can be transferred across populations, strongly suggesting that these PCHs reflect traits with common underlying genetic mechanisms for OSAHS across populations. Thus, PCHs can provide useful traits for using data on multiple phenotypes and for genetic studies of trans-ethnic populations.

Molecular Genetic Contributions to Social Deprivation and Household Income in UK Biobank

Individuals with lower socio-economic status (SES) are at increased risk of physical and mental illnesses and tend to die at an earlier age [1, 2, 3]. Explanations for the association between SES and health typically focus on factors that are environmental in origin [4]. However, common SNPs have been found collectively to explain around 18% of the phenotypic variance of an area-based social deprivation measure of SES [5]. Molecular genetic studies have also shown that common physical and psychiatric diseases are partly heritable [6]. It is possible that phenotypic associations between SES and health arise partly due to a shared genetic etiology. We conducted a genome-wide association study (GWAS) on social deprivation and on household income using 112,151 participants of UK Biobank. We find that common SNPs explain 21% of the variation in social deprivation and 11% of household income. Two independent loci attained genome-wide significance for household income, with the most significant SNP in each of these loci being rs187848990 on chromosome 2 and rs8100891 on chromosome 19. Genes in the regions of these SNPs have been associated with intellectual disabilities, schizophrenia, and synaptic plasticity. Extensive genetic correlations were found between both measures of SES and illnesses, anthropometric variables, psychiatric disorders, and cognitive ability. These findings suggest that some SNPs associated with SES are involved in the brain and central nervous system. The genetic associations with SES obviously do not reflect direct causal effects and are probably mediated via other partly heritable variables, including cognitive ability, personality, and health.

•

Common SNPs explain 21% of social deprivation and 11% of household income

•

Two loci attained genome-wide significance for household income

•

Genes in these loci have been linked to synaptic plasticity

•

Genetic correlations were found between both measures of SES and many other traits

Shared genetic aetiology between cognitive functions and physical and mental health in UK Biobank (N=112 151) and 24 GWAS consortia

Causes of the well-documented association between low levels of cognitive functioning and many adverse neuropsychiatric outcomes, poorer physical health and earlier death remain unknown. We used linkage disequilibrium regression and polygenic profile scoring to test for shared genetic aetiology between cognitive functions and neuropsychiatric disorders and physical health. Using information provided by many published genome-wide association study consortia, we created polygenic profile scores for 24 vascular-metabolic, neuropsychiatric, physiological-anthropometric and cognitive traits in the participants of UK Biobank, a very large population-based sample (N=112 151). Pleiotropy between cognitive and health traits was quantified by deriving genetic correlations using summary genome-wide association study statistics and to the method of linkage disequilibrium score regression. Substantial and significant genetic correlations were observed between cognitive test scores in the UK Biobank sample and many of the mental and physical health-related traits and disorders assessed here. In addition, highly significant associations were observed between the cognitive test scores in the UK Biobank sample and many polygenic profile scores, including coronary artery disease, stroke, Alzheimer's disease, schizophrenia, autism, major depressive disorder, body mass index, intracranial volume, infant head circumference and childhood cognitive ability. Where disease diagnosis was available for UK Biobank participants, we were able to show that these results were not confounded by those who had the relevant disease. These findings indicate that a substantial level of pleiotropy exists between cognitive abilities and many human mental and physical health disorders and traits and that it can be used to predict phenotypic variance across samples.

Molecular genetic approaches to understanding the comorbidity of psychiatric disorders

Epidemiologic studies demonstrating high rates of co-occurrence among psychiatric disorders at the population level have contributed to large literatures focused on identifying the causal mechanisms underlying the patterns of co-occurrence among these disorders. Such efforts have long represented a core focus of developmental psychopathologists and have more recently been supported by the Research Domain Criteria initiative developed by the NIMH, which provides a further framework for how the hypothesized mechanisms can be studied at different levels of analysis. The present overview focuses on molecular genetic approaches that are being used currently to study the etiology of psychiatric disorders, and how these approaches have been applied in efforts to understand the biological mechanisms that give rise to comorbid conditions. The present report begins with a review of molecular genetic approaches used to identify individual variants that confer risk for multiple disorders and the intervening biological mechanisms that contribute to their comorbidity. This is followed by a review of molecular genetic approaches that use genetic data in aggregate to examine these questions, and concludes with a discussion of how developmental psychopathologists are uniquely positioned to apply these methods in a way that will further our understanding of the causal factors that contribute to the development of comorbid conditions.

Pleiotropic Meta-Analyses of Longitudinal Studies Discover Novel Genetic Variants Associated with Age-Related Diseases

Age-related diseases may result from shared biological mechanisms in intrinsic processes of aging. Genetic effects on age-related diseases are often modulated by environmental factors due to their little contribution to fitness or are mediated through certain endophenotypes. Identification of genetic variants with pleiotropic effects on both common complex diseases and endophenotypes may reveal potential conflicting evolutionary pressures and deliver new insights into shared genetic contribution to healthspan and lifespan. Here, we performed pleiotropic meta-analyses of genetic variants using five NIH-funded datasets by integrating univariate summary statistics for age-related diseases and endophenotypes. We investigated three groups of traits: (1) endophenotypes such as blood glucose, blood pressure, lipids, hematocrit, and body mass index, (2) time-to-event outcomes such as the age-at-onset of diabetes mellitus (DM), cancer, cardiovascular diseases (CVDs) and neurodegenerative diseases (NDs), and (3) both combined. In addition to replicating previous findings, we identify seven novel genome-wide significant loci (< 5e-08), out of which five are low-frequency variants. Specifically, from Group 2, we find rs7632505 on 3q21.1 in SEMA5B, rs460976 on 21q22.3 (1 kb from TMPRSS2) and rs12420422 on 11q24.1 predominantly associated with a variety of CVDs, rs4905014 in ITPK1 associated with stroke and heart failure, rs7081476 on 10p12.1 in ANKRD26 associated with multiple diseases including DM, CVDs, and NDs. From Group 3, we find rs8082812 on 18p11.22 and rs1869717 on 4q31.3 associated with both endophenotypes and CVDs. Our follow-up analyses show that rs7632505, rs4905014, and rs8082812 have age-dependent effects on coronary heart disease or stroke. Functional annotation suggests that most of these SNPs are within regulatory regions or DNase clusters and in linkage disequilibrium with expression quantitative trait loci, implying their potential regulatory influence on the expression of nearby genes. Our mediation analyses suggest that the effects of some SNPs are mediated by specific endophenotypes. In conclusion, these findings indicate that loci with pleiotropic effects on age-related disorders tend to be enriched in genes involved in underlying mechanisms potentially related to nervous, cardiovascular and immune system functions, stress resistance, inflammation, ion channels and hematopoiesis, supporting the hypothesis of shared pathological role of infection, and inflammation in chronic age-related diseases.

Multivariate Analysis of Anthropometric Traits Using Summary Statistics of Genome-Wide Association Studies from GIANT Consortium

Meta-analysis of single trait for multiple cohorts has been used for increasing statistical power in genome-wide association studies (GWASs). Although hundreds of variants have been identified by GWAS, these variants only explain a small fraction of phenotypic variation. Cross-phenotype association analysis (CPASSOC) can further improve statistical power by searching for variants that contribute to multiple traits, which is often relevant to pleiotropy. In this study, we performed CPASSOC analysis on the summary statistics from the Genetic Investigation of ANthropometric Traits (GIANT) consortium using a novel method recently developed by our group. Sex-specific meta-analysis data for height, body mass index (BMI), and waist-to-hip ratio adjusted for BMI (WHRadjBMI) from discovery phase of the GIANT consortium study were combined using CPASSOC for each trait as well as 3 traits together. The conventional meta-analysis results from the discovery phase data of GIANT consortium studies were used to compare with that from CPASSOC analysis. The CPASSOC analysis was able to identify 17 loci associated with anthropometric traits that were missed by conventional meta-analysis. Among these loci, 16 have been reported in literature by including additional samples and 1 is novel. We also demonstrated that CPASSOC is able to detect pleiotropic effects when analyzing multiple traits.

A machine learning-based framework to identify type 2 diabetes through electronic health records

OBJECTIVE

To discover diverse genotype-phenotype associations affiliated with Type 2 Diabetes Mellitus (T2DM) via genome-wide association study (GWAS) and phenome-wide association study (PheWAS), more cases (T2DM subjects) and controls (subjects without T2DM) are required to be identified (e.g., via Electronic Health Records (EHR)). However, existing expert based identification algorithms often suffer in a low recall rate and could miss a large number of valuable samples under conservative filtering standards. The goal of this work is to develop a semi-automated framework based on machine learning as a pilot study to liberalize filtering criteria to improve recall rate with a keeping of low false positive rate.

MATERIALS AND METHODS

We propose a data informed framework for identifying subjects with and without T2DM from EHR via feature engineering and machine learning. We evaluate and contrast the identification performance of widely-used machine learning models within our framework, including k-Nearest-Neighbors, Naïve Bayes, Decision Tree, Random Forest, Support Vector Machine and Logistic Regression. Our framework was conducted on 300 patient samples (161 cases, 60 controls and 79 unconfirmed subjects), randomly selected from 23,281 diabetes related cohort retrieved from a regional distributed EHR repository ranging from 2012 to 2014.

RESULTS

We apply top-performing machine learning algorithms on the engineered features. We benchmark and contrast the accuracy, precision, AUC, sensitivity and specificity of classification models against the state-of-the-art expert algorithm for identification of T2DM subjects. Our results indicate that the framework achieved high identification performances (∼0.98 in average AUC), which are much higher than the state-of-the-art algorithm (0.71 in AUC).

DISCUSSION

Expert algorithm-based identification of T2DM subjects from EHR is often hampered by the high missing rates due to their conservative selection criteria. Our framework leverages machine learning and feature engineering to loosen such selection criteria to achieve a high identification rate of cases and controls.

CONCLUSIONS

Our proposed framework demonstrates a more accurate and efficient approach for identifying subjects with and without T2DM from EHR.

Challenges and novel approaches for investigating molecular mediation

Understanding mediation is useful for identifying intermediates lying between an exposure and an outcome which, when intervened upon, will block (some or all of) the causal pathway between the exposure and outcome. Mediation approaches used in conventional epidemiology have been adapted to understanding the role of molecular intermediates in situations of high-dimensional omics data with varying degrees of success. In particular, the limitations of observational epidemiological study including confounding, reverse causation and measurement error can afflict conventional mediation approaches and may lead to incorrect conclusions regarding causal effects. Solutions to analysing mediation which overcome these problems include the use of instrumental variable methods such as Mendelian randomization, which may be applied to evaluate causality in increasingly complex networks of omics data.

Cognitive ability and risk for substance misuse in men: genetic and environmental correlations in a longitudinal nation-wide family study

AIMS

To investigate the association in males between cognitive ability in late adolescence and subsequent substance misuse-related events, and to study the underlying genetic and environmental correlations.

DESIGN

A population-based longitudinal study with three different family-based designs. Cox proportional hazards models were conducted to investigate the association at the individual level. Bivariate quantitative genetic modelling in (1) full brothers and maternal half-brothers, (2) full brothers reared together and apart and (3) monozygotic and dizygotic twin brothers was used to estimate genetic and environmental correlations.

SETTING

Register-based study in Sweden.

PARTICIPANTS

The full sample included 1 402 333 Swedish men born 1958-91 and conscripted at mean age 18.2 [standard deviation (SD) = 0.5] years. A total of 1 361 066 men who had no substance misuse events before cognitive assessment at mandatory military conscription were included in the Cox regression models, with a follow-up time of up to 35.6 years.

MEASURES

Cognitive ability was assessed at conscription with the Swedish Enlistment Battery. Substance misuse events included alcohol- and drug-related court convictions, medical treatments and deaths, available from governmental registries.

FINDINGS

Lower cognitive ability in late adolescence predicted an increased risk for substance misuse events [hazard ratio (HR) for a 1-stanine unit decrease in cognitive ability: 1.29, 95% confidence interval (CI) = 1.29-1.30]. The association was somewhat attenuated within clusters of full brothers (HR = 1.21, 95% CI = 1.20-1.23). Quantitative genetic analyses indicated that the association was due primarily to genetic influences; the genetic correlations ranged between -0.39 (95% CI = -0.45, -0.34) and -0.52 (95% CI -0.55, -0.48) in the three different designs.

CONCLUSIONS

Shared genetic influences appear to underlie the association between low cognitive ability and subsequent risk for substance misuse events among Swedish men.

GNB5 mutation causes a novel neuropsychiatric disorder featuring attention deficit hyperactivity disorder, severely impaired language development and normal cognition

Background

Neuropsychiatric disorders are common forms of disability in humans. Despite recent progress in deciphering the genetics of these disorders, their phenotypic complexity continues to be a major challenge. Mendelian neuropsychiatric disorders are rare but their study has the potential to unravel novel mechanisms that are relevant to their complex counterparts.

Results

In an extended consanguineous family, we identified a novel neuropsychiatric phenotype characterized by severe speech impairment, variable expressivity of attention deficit hyperactivity disorder (ADHD), and motor delay. We identified the disease locus through linkage analysis on 15q21.2, and exome sequencing revealed a novel missense variant in GNB5. GNB5 encodes an atypical β subunit of the heterotrimeric GTP-binding proteins (Gβ5). Gβ5 is enriched in the central nervous system where it forms constitutive complexes with members of the regulator of G protein signaling family of proteins to modulate neurotransmitter signaling that affects a number of neurobehavioral outcomes. Here, we show that the S81L mutant form of Gβ5 has significantly impaired activity in terminating responses that are elicited by dopamine.

Conclusions

We demonstrate that these deficits originate from the impaired expression of the mutant Gβ5 protein, resulting in the decreased ability to stabilize regulator of G protein signaling complexes. Our data suggest that this novel neuropsychiatric phenotype is the human equivalent of Gnb5 deficiency in mice, which manifest motor deficits and hyperactivity, and highlight a critical role of Gβ5 in normal behavior as well as language and motor development in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1061-6) contains supplementary material, which is available to authorized users.

Improved methods for multi-trait fine mapping of pleiotropic risk loci

MOTIVATION

Genome-wide association studies (GWAS) have identified thousands of regions in the genome that contain genetic variants that increase risk for complex traits and diseases. However, the variants uncovered in GWAS are typically not biologically causal, but rather, correlated to the true causal variant through linkage disequilibrium (LD). To discern the true causal variant(s), a variety of statistical fine-mapping methods have been proposed to prioritize variants for functional validation.

RESULTS

In this work we introduce a new approach, fastPAINTOR, that leverages evidence across correlated traits, as well as functional annotation data, to improve fine-mapping accuracy at pleiotropic risk loci. To improve computational efficiency, we describe an new importance sampling scheme to perform model inference. First, we demonstrate in simulations that by leveraging functional annotation data, fastPAINTOR increases fine-mapping resolution relative to existing methods. Next, we show that jointly modeling pleiotropic risk regions improves fine-mapping resolution compared to standard single trait and pleiotropic fine mapping strategies. We report a reduction in the number of SNPs required for follow-up in order to capture 90% of the causal variants from 23 SNPs per locus using a single trait to 12 SNPs when fine-mapping two traits simultaneously. Finally, we analyze summary association data from a large-scale GWAS of lipids and show that these improvements are largely sustained in real data.

AVAILABILITY AND IMPLEMENTATION

The fastPAINTOR framework is implemented in the PAINTOR v3.0 package which is publicly available to the research community http://bogdan.bioinformatics.ucla.edu/software/paintor CONTACT: gkichaev@ucla.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types

Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10(-8) seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type-specific expression quantitative trait locus and enhancer-gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10(-5) in the three-cancer meta-analysis.

SIGNIFICANCE

We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052-67. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 932.

Strong genetic overlap between executive functions and intelligence

Executive functions (EFs) are cognitive processes that control, monitor, and coordinate more basic cognitive processes. EFs play instrumental roles in models of complex reasoning, learning, and decision making, and individual differences in EFs have been consistently linked with individual differences in intelligence. By middle childhood, genetic factors account for a moderate proportion of the variance in intelligence, and these effects increase in magnitude through adolescence. Genetic influences on EFs are very high, even in middle childhood, but the extent to which these genetic influences overlap with those on intelligence is unclear. We examined genetic and environmental overlap between EFs and intelligence in a racially and socioeconomically diverse sample of 811 twins ages 7 to 15 years (M = 10.91, SD = 1.74) from the Texas Twin Project. A general EF factor representing variance common to inhibition, switching, working memory, and updating domains accounted for substantial proportions of variance in intelligence, primarily via a genetic pathway. General EF continued to have a strong, genetically mediated association with intelligence even after controlling for processing speed. Residual variation in general intelligence was influenced only by shared and nonshared environmental factors, and there remained no genetic variance in general intelligence that was unique of EF. Genetic variance independent of EF did remain, however, in a more specific perceptual reasoning ability. These results provide evidence that genetic influences on general intelligence are highly overlapping with those on EF. (PsycINFO Database Record

Leveraging electronic health records to study pleiotropic effects on bipolar disorder and medical comorbidities

Patients with bipolar disorder (BD) have a high prevalence of comorbid medical illness. However, the mechanisms underlying these comorbidities with BD are not well known. Certain genetic variants may have pleiotropic effects, increasing the risk of BD and other medical illnesses simultaneously. In this study, we evaluated the association of BD-susceptibility genetic variants with various medical conditions that tend to co-exist with BD, using electronic health records (EHR) data linked to genome-wide single-nucleotide polymorphism (SNP) data. Data from 7316 Caucasian subjects were used to test the association of 19 EHR-derived phenotypes with 34 SNPs that were previously reported to be associated with BD. After Bonferroni multiple testing correction, P<7.7 × 10(-5) was considered statistically significant. The top association findings suggested that the BD risk alleles at SNP rs4765913 in CACNA1C gene and rs7042161 in SVEP1 may be associated with increased risk of 'cardiac dysrhythmias' (odds ratio (OR)=1.1, P=3.4 × 10(-3)) and 'essential hypertension' (OR=1.1, P=3.5 × 10(-3)), respectively. Although these associations are not statistically significant after multiple testing correction, both genes have been previously implicated with cardiovascular phenotypes. Moreover, we present additional evidence supporting these associations, particularly the association of the SVEP1 SNP with hypertension. This study shows the potential for EHR-based analyses of large cohorts to discover pleiotropic effects contributing to complex psychiatric traits and commonly co-occurring medical conditions.

Statistical Methods for Testing Genetic Pleiotropy

Genetic pleiotropy is when a single gene influences more than one trait. Detecting pleiotropy and understanding its causes can improve the biological understanding of a gene in multiple ways, yet current multivariate methods to evaluate pleiotropy test the null hypothesis that none of the traits are associated with a variant; departures from the null could be driven by just one associated trait. A formal test of pleiotropy should assume a null hypothesis that one or no traits are associated with a genetic variant. For the special case of two traits, one can construct this null hypothesis based on the intersection-union (IU) test, which rejects the null hypothesis only if the null hypotheses of no association for both traits are rejected. To allow for more than two traits, we developed a new likelihood-ratio test for pleiotropy. We then extended the testing framework to a sequential approach to test the null hypothesis that [Formula: see text] traits are associated, given that the null of k traits are associated was rejected. This provides a formal testing framework to determine the number of traits associated with a genetic variant, while accounting for correlations among the traits. By simulations, we illustrate the type I error rate and power of our new methods; describe how they are influenced by sample size, the number of traits, and the trait correlations; and apply the new methods to multivariate immune phenotypes in response to smallpox vaccination. Our new approach provides a quantitative assessment of pleiotropy, enhancing current analytic practice.

eMERGE Phenome-Wide Association Study (PheWAS) identifies clinical associations and pleiotropy for stop-gain variants

BACKGROUND

We explored premature stop-gain variants to test the hypothesis that variants, which are likely to have a consequence on protein structure and function, will reveal important insights with respect to the phenotypes associated with them. We performed a phenome-wide association study (PheWAS) exploring the association between a selected list of functional stop-gain genetic variants (variation resulting in truncated proteins or in nonsense-mediated decay) and an extensive group of diagnoses to identify novel associations and uncover potential pleiotropy.

RESULTS

In this study, we selected 25 stop-gain variants: 5 stop-gain variants with previously reported phenotypic associations, and a set of 20 putative stop-gain variants identified using dbSNP. For the PheWAS, we used data from the electronic MEdical Records and GEnomics (eMERGE) Network across 9 sites with a total of 41,057 unrelated patients. We divided all these samples into two datasets by equal proportion of eMERGE site, sex, race, and genotyping platform. We calculated single effect associations between these 25 stop-gain variants and ICD-9 defined case-control diagnoses. We also performed stratified analyses for samples of European and African ancestry. Associations were adjusted for sex, site, genotyping platform and the first three principal components to account for global ancestry. We identified previously known associations, such as variants in LPL associated with hyperglyceridemia indicating that our approach was robust. We also found a total of three significant associations with p < 0.01 in both datasets, with the most significant replicating result being LPL SNP rs328 and ICD-9 code 272.1 "Disorder of Lipoid metabolism" (pdiscovery = 2.59x10-6, preplicating = 2.7x10-4). The other two significant replicated associations identified by this study are: variant rs1137617 in KCNH2 gene associated with ICD-9 code category 244 "Acquired Hypothyroidism" (pdiscovery = 5.31x103, preplicating = 1.15x10-3) and variant rs12060879 in DPT gene associated with ICD-9 code category 996 "Complications peculiar to certain specified procedures" (pdiscovery = 8.65x103, preplicating = 4.16x10-3).

CONCLUSION

In conclusion, this PheWAS revealed novel associations of stop-gained variants with interesting phenotypes (ICD-9 codes) along with pleiotropic effects.

Network biology concepts in complex disease comorbidities

The co-occurrence of diseases can inform the underlying network biology of shared and multifunctional genes and pathways. In addition, comorbidities help to elucidate the effects of external exposures, such as diet, lifestyle and patient care. With worldwide health transaction data now often being collected electronically, disease co-occurrences are starting to be quantitatively characterized. Linking network dynamics to the real-life, non-ideal patient in whom diseases co-occur and interact provides a valuable basis for generating hypotheses on molecular disease mechanisms, and provides knowledge that can facilitate drug repurposing and the development of targeted therapeutic strategies.