Pleiotropy in complex traits: challenges and strategies

Polygenic risk score identifies associations between sleep duration and diseases determined from an electronic medical record biobank

STUDY OBJECTIVES

We aimed to detect cross-sectional phenotype and polygenic risk score (PRS) associations between sleep duration and prevalent diseases using the Partners Biobank, a hospital-based cohort study linking electronic medical records (EMR) with genetic information.

METHODS

Disease prevalence was determined from EMR, and sleep duration was self-reported. A PRS for sleep duration was derived using 78 previously associated SNPs from genome-wide association studies (GWAS) for self-reported sleep duration. We tested for associations between (1) self-reported sleep duration and 22 prevalent diseases (n = 30 251), (2) the PRS and self-reported sleep duration (n = 6903), and (3) the PRS and the 22 prevalent diseases (n = 16 033). For observed PRS-disease associations, we tested causality using two-sample Mendelian randomization (MR).

RESULTS

In the age-, sex-, and race-adjusted model, U-shaped associations were observed for sleep duration and asthma, depression, hypertension, insomnia, obesity, obstructive sleep apnea, and type 2 diabetes, where both short and long sleepers had higher odds for these diseases than normal sleepers (p < 2.27 × 10-3). Next, we confirmed associations between the PRS and longer sleep duration (0.65 ± 0.19 SD minutes per effect allele; p = 7.32 × 10-04). The PRS collectively explained 1.4% of the phenotypic variance in sleep duration. After adjusting for age, sex, genotyping array, and principal components of ancestry, we observed that the PRS was also associated with congestive heart failure (CHF; p = 0.015), obesity (p = 0.019), hypertension (p = 0.039), restless legs syndrome (RLS; p = 0.041), and insomnia (p = 0.049). Associations were maintained following additional adjustment for obesity status, except for hypertension and insomnia. For all diseases, except RLS, carrying a higher genetic burden of the 78 sleep duration-increasing alleles (i.e. higher sleep duration PRS) associated with lower odds for prevalent disease. In MR, we estimated causal associations between genetically defined longer sleep duration with decreased risk of CHF (inverse variance weighted [IVW] OR per minute of sleep [95% CI] = 0.978 [0.961-0.996]; p = 0.019) and hypertension (IVW OR [95% CI] = 0.993 [0.986-1.000]; p = 0.049), and increased risk of RLS (IVW OR [95% CI] = 1.018 [1.000-1.036]; p = 0.045).

CONCLUSIONS

By validating the PRS for sleep duration and identifying cross-phenotype associations, we lay the groundwork for future investigations on the intersection between sleep, genetics, clinical measures, and diseases using large EMR datasets.

Comprehensive identification of pleiotropic loci for body fat distribution using the NHGRI-EBI Catalog of published genome-wide association studies

We conducted a hypothesis-free cross-trait analysis for waist-to-hip ratio adjusted for body mass index (WHR_adjBMI ) loci derived through genome-wide association studies (GWAS). Summary statistics from published GWAS were used to capture all WHR_adjBMI single-nucleotide polymorphisms (SNPs), and their proxy SNPs were identified. These SNPs were used to extract cross-trait associations between WHR_adjBMI SNPs and other traits through the NHGRI-EBI GWAS Catalog. Pathway analysis was conducted for pleiotropic WHR_adjBMI SNPs. We found 160 WHR_adjBMI SNPs and 3675 proxy SNPs. Cross-trait analysis identified 239 associations, of which 100 were for obesity traits. The remaining 139 associations were filtered down to 101 unique linkage disequilibrium block associations, which were grouped into 13 categories: lipids, red blood cell traits, white blood cell counts, inflammatory markers and autoimmune diseases, type 2 diabetes-related traits, adiponectin, cancers, blood pressure, height, neuropsychiatric disorders, electrocardiography changes, urea measurement, and others. The highest number of cross-trait associations were found for triglycerides (n = 10), high-density lipoprotein cholesterol (n = 9), and reticulocyte counts (n = 8). Pathway analysis for WHR_adjBMI pleiotropic SNPs found immune function pathways as the top canonical pathways. Results from our original methodology indicate a novel genetic association between WHR_adjBMI and reticulocyte counts and highlight the pleiotropy between abdominal obesity, immune pathways, and other traits.

GWAS and Beyond: Using Omics Approaches to Interpret SNP Associations

PURPOSE OF REVIEW

Neurodegenerative diseases, neuropsychiatric disorders, and related traits have highly complex etiologies but are also highly heritable and identifying the causal genes and biological pathways underlying these traits may advance the development of treatments and preventive strategies. While many genome-wide association studies (GWAS) have successfully identified variants contributing to polygenic neurodegenerative and neuropsychiatric phenotypes including Alzheimer's disease (AD), schizophrenia (SCZ), and bipolar disorder (BPD) amongst others, interpreting the biological roles of significantly-associated variants in the genetic architecture of these traits remains a significant challenge. Here we review several 'omics' approaches which attempt to bridge the gap from associated genetic variants to phenotype by helping define the functional roles of GWAS loci in the development of neuropsychiatric disorders and traits.

RECENT FINDINGS

Several common 'omics' approaches have been applied to examine neuropsychiatric traits, such as nearest-gene mapping, trans-ethnic fine mapping, annotation enrichment analysis, transcriptomic analysis, and pathway analysis, and each of these approaches has strengths and limitations in providing insight into biological mechanisms. One popular emerging method is the examination of tissue-specific genetically-regulated gene expression (GReX), which aggregates the genetic variants' effects at the gene-level. Furthermore, proteomic, metabolomic, and microbiomic studies and phenome-wide association studies will further enhance our understanding of neuropsychiatric traits.

SUMMARY

GWAS has been applied to neuropsychiatric traits for a decade, but our understanding about the biological function of identified variants remains limited. Today, technological advancements have created analytical approaches for integrating transcriptomics, metabolomics, proteomics, pharmacology and toxicology as tools for understanding the functional roles of genetics variants. These data, as well as the broader clinical information provided by electronic health records, can provide additional insight and complement genomic analyses.

A Geometric Perspective on the Power of Principal Component Association Tests in Multiple Phenotype Studies

Joint analysis of multiple phenotypes can increase statistical power in genetic association studies. Principal component analysis, as a popular dimension reduction method, especially when the number of phenotypes is high-dimensional, has been proposed to analyze multiple correlated phenotypes. It has been empirically observed that the first PC, which summarizes the largest amount of variance, can be less powerful than higher order PCs and other commonly used methods in detecting genetic association signals. In this paper, we investigate the properties of PCA-based multiple phenotype analysis from a geometric perspective by introducing a novel concept called principal angle. A particular PC is powerful if its principal angle is 0° and is powerless if its principal angle is 90°. Without prior knowledge about the true principal angle, each PC can be powerless. We propose linear, non-linear and data-adaptive omnibus tests by combining PCs. We demonstrate that the Wald test is a special quadratic PC-based test. We show that the omnibus PC test is robust and powerful in a wide range of scenarios. We study the properties of the proposed methods using power analysis and eigen-analysis. The subtle differences and close connections between these combined PC methods are illustrated graphically in terms of their rejection boundaries. Our proposed tests have convex acceptance regions and hence are admissible. The p-values for the proposed tests can be efficiently calculated analytically and the proposed tests have been implemented in a publicly available R package MPAT. We conduct simulation studies in both low and high dimensional settings with various signal vectors and correlation structures. We apply the proposed tests to the joint analysis of metabolic syndrome related phenotypes with data sets collected from four international consortia to demonstrate the effectiveness of the proposed combined PC testing procedures.

Exploring Impact of Rare Variation in Systemic Lupus Erythematosus by a Genome Wide Imputation Approach

The importance of low frequency and rare variation in complex disease genetics is difficult to estimate in patient populations. Genome-wide association studies are therefore, underpowered to detect rare variation. We have used a combined approach of genome-wide-based imputation with a highly stringent sequence kernel association (SKAT) test and a case-control burden test. We identified 98 candidate genes containing rare variation that in aggregate show association with SLE many of which have recognized immunological function, but also function and expression related to relevant tissues such as the joints, skin, blood or central nervous system. In addition we also find that there is a significant enrichment of genes annotated for disease-causing mutations in the OMIM database, suggesting that in complex diseases such as SLE, such mutations may be involved in subtle or combined phenotypes or could accelerate specific organ abnormalities found in the disease. We here provide an important resource of candidate genes for SLE.

Multivariate Methods for Meta-Analysis of Genetic Association Studies

Multivariate meta-analysis of genetic association studies and genome-wide association studies has received a remarkable attention as it improves the precision of the analysis. Here, we review, summarize and present in a unified framework methods for multivariate meta-analysis of genetic association studies and genome-wide association studies. Starting with the statistical methods used for robust analysis and genetic model selection, we present in brief univariate methods for meta-analysis and we then scrutinize multivariate methodologies. Multivariate models of meta-analysis for a single gene-disease association studies, including models for haplotype association studies, multiple linked polymorphisms and multiple outcomes are discussed. The popular Mendelian randomization approach and special cases of meta-analysis addressing issues such as the assumption of the mode of inheritance, deviation from Hardy-Weinberg Equilibrium and gene-environment interactions are also presented. All available methods are enriched with practical applications and methodologies that could be developed in the future are discussed. Links for all available software implementing multivariate meta-analysis methods are also provided.

Assessing pleiotropy and mediation in genetic loci associated with chronic obstructive pulmonary disease

Genetic association studies have increasingly recognized variant effects on multiple phenotypes. Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease with environmental and genetic causes. Multiple genetic variants have been associated with COPD, many of which show significant associations to additional phenotypes. However, it is unknown if these associations represent biological pleiotropy or if they exist through correlation of related phenotypes ("mediated pleiotropy"). Using 6,670 subjects from the COPDGene study, we describe the association of known COPD susceptibility loci with other COPD-related phenotypes and distinguish if these act directly on the phenotypes (i.e., biological pleiotropy) or if the association is due to correlation (i.e., mediated pleiotropy). We identified additional associated phenotypes for 13 of 25 known COPD loci. Tests for pleiotropy between genotype and associated outcomes were significant for all loci. In cases of significant pleiotropy, we performed mediation analysis to test if SNPs had a direct association to phenotype. Most loci showed a mediated effect through the hypothesized causal pathway. However, many loci also had direct associations, suggesting causal explanations (i.e., emphysema leading to reduced lung function) are incomplete. Our results highlight the high degree of pleiotropy in complex disease-associated loci and provide novel insights into the mechanisms underlying COPD.

Evolutionary highways to persistent bacterial infection

Persistent infections require bacteria to evolve from their naïve colonization state by optimizing fitness in the host via simultaneous adaptation of multiple traits, which can obscure evolutionary trends and complicate infection management. Accordingly, here we screen 8 infection-relevant phenotypes of 443 longitudinal Pseudomonas aeruginosa isolates from 39 young cystic fibrosis patients over 10 years. Using statistical modeling, we map evolutionary trajectories and identify trait correlations accounting for patient-specific influences. By integrating previous genetic analyses of 474 isolates, we provide a window into early adaptation to the host, finding: (1) a 2–3 year timeline of rapid adaptation after colonization, (2) variant “naïve” and “adapted” states reflecting discordance between phenotypic and genetic adaptation, (3) adaptive trajectories leading to persistent infection via three distinct evolutionary modes, and (4) new associations between phenotypes and pathoadaptive mutations. Ultimately, we effectively deconvolute complex trait adaptation, offering a framework for evolutionary studies and precision medicine in clinical microbiology.

The pathogen Pseudomonas aeruginosa undergoes complex trait adaptation within cystic fibrosis patients. Here, Bartell, Sommer, and colleagues use statistical modeling of longitudinal isolates to characterize the joint genetic and phenotypic evolutionary trajectories of P. aeruginosa within hosts.

Investigation of multi-trait associations using pathway-based analysis of GWAS summary statistics

Background

Genome-wide association studies (GWAS) have been successful in identifying disease-associated genetic variants. Recently, an increasing number of GWAS summary statistics have been made available to the research community, providing extensive repositories for studies of human complex diseases. In particular, cross-trait associations at the genetic level can be beneficial from large-scale GWAS summary statistics by using genetic variants that are associated with multiple traits. However, direct assessment of cross-trait associations using susceptibility loci has been challenging due to the complex genetic architectures in most diseases, calling for advantageous methods that could integrate functional interpretation and imply biological mechanisms.

Results

We developed an analytical framework for systematic integration of cross-trait associations. It incorporates two different approaches to detect enriched pathways and requires only summary statistics. We demonstrated the framework using 25 traits belonging to four phenotype groups. Our results revealed an average of 54 significantly associated pathways (ranged between 18 and 175) per trait. We further proved that pathway-based analysis provided increased power to estimate cross-trait associations compared to gene-level analysis. Based on Fisher’s Exact Test (FET), we identified a total of 24 (53) pairs of trait-trait association at adjusted p_FET < 1 × 10^− 3 (p_FET < 0.01) among the 25 traits. Our trait-trait association network revealed not only many relationships among the traits within the same group but also novel relationships among traits from different groups, which warrants further investigation in future.

Conclusions

Our study revealed that risk variants for 25 different traits aggregated in particular biological pathways and that these pathways were frequently shared among traits. Our results confirmed known mechanisms and also suggested several novel insights into the etiology of multi-traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5373-7) contains supplementary material, which is available to authorized users.

"Mendelian Randomization" Approach in Economic Assessment of Health Conditions

The increased prevalence of non-communicable chronic diseases (NCDs) is reflected in the rising economic burden of health conditions. Observational studies conducted in health economics research are detecting associations of NCDs or related risk factors with economic measures like health insurance, economic inequalities, accessibility of jobs, education, annual income, health expenditure, etc. The inferences of such relationships do not prove causation and are limited to associations which are many times influenced by confounding factors and reverse causation. Mendelian randomization (MR) approach is a useful method for exploring causal relations between modifiable risk factors and measures of health economics. The application of MR in economic assessment of health conditions has been started and is producing fruitful results.

Multi-SKAT: General framework to test for rare-variant association with multiple phenotypes

In genetic association analysis, a joint test of multiple distinct phenotypes can increase power to identify sets of trait-associated variants within genes or regions of interest. Existing multiphenotype tests for rare variants make specific assumptions about the patterns of association with underlying causal variants, and the violation of these assumptions can reduce power to detect association. Here, we develop a general framework for testing pleiotropic effects of rare variants on multiple continuous phenotypes using multivariate kernel regression (Multi-SKAT). Multi-SKAT models affect sizes of variants on the phenotypes through a kernel matrix and perform a variance component test of association. We show that many existing tests are equivalent to specific choices of kernel matrices with the Multi-SKAT framework. To increase power of detecting association across tests with different kernel matrices, we developed a fast and accurate approximation of the significance of the minimum observed P value across tests. To account for related individuals, our framework uses random effects for the kinship matrix. Using simulated data and amino acid and exome-array data from the METabolic Syndrome In Men (METSIM) study, we show that Multi-SKAT can improve power over single-phenotype SKAT-O test and existing multiple-phenotype tests, while maintaining Type I error rate.

Precision medicine review: rare driver mutations and their biophysical classification

How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.

Crop genome-wide association study: a harvest of biological relevance

With the advent of rapid genotyping and next-generation sequencing technologies, genome-wide association study (GWAS) has become a routine strategy for decoding genotype-phenotype associations in many species. More than 1000 such studies over the last decade have revealed substantial genotype-phenotype associations in crops and provided unparalleled opportunities to probe functional genomics. Beyond the many 'hits' obtained, this review summarizes recent efforts to increase our understanding of the genetic architecture of complex traits by focusing on non-main effects including epistasis, pleiotropy, and phenotypic plasticity. We also discuss how these achievements and the remaining gaps in our knowledge will guide future studies. Synthetic association is highlighted as leading to false causality, which is prevalent but largely underestimated. Furthermore, validation evidence is appealing for future GWAS, especially in the context of emerging genome-editing technologies.

Identification of pleiotropic genetic variants affecting osteoporosis risk in a Korean elderly cohort

Pleiotropy has important implications for understanding the genetic basis and risk assessment of osteoporosis. Our aim was to identify pleiotropic genetic variants associated with the development of osteoporosis and predict osteoporosis risk by leveraging pleiotropic variants. We evaluated the effects of 21 conventional risk factors and 185 single-nucleotide polymorphisms (SNPs) in 63 inflammation- and metabolism-related genes on osteoporosis risk in a community-based Korean cohort study of 1025 participants, the Hallym Aging Study. Ten nongenetic factors, including sex (female) and hematocrit level, and 12 SNPs across ten genes showed evidence of association with incident osteoporosis in 270 initially osteoporosis-free subjects who completed a 6-year follow up. Three gene variants, rs1801282 (PPARG-Pro12Ala, hazard ratio (HR) = 3.26, P = 0.008), rs1408282 (near EPHA7, HR = 1.87, P = 0.002), and rs2076212 (PNPLA3-Gly115Cys, HR = 2.24, P = 0.024), were associated with significant differences in survival among the three genotype groups (P_diff = 0.042, 0.003, and 0.048, respectively). Individuals in the highest polygenic risk score tertile were 27.9 fold more likely to develop osteoporosis than those in the lowest tertile (P = 0.004). The PPARG gene in particular was a hub pleiotropic gene in the epistasis network. Our findings highlight pleiotropic modulations of metabolism- and inflammation-related genes in the development of osteoporosis and demonstrate the contribution of pleiotropic genetic variants in prediction of osteoporosis risk.

Detecting potential pleiotropy across cardiovascular and neurological diseases using univariate, bivariate, and multivariate methods on 43,870 individuals from the eMERGE network

The link between cardiovascular diseases and neurological disorders has been widely observed in the aging population. Disease prevention and treatment rely on understanding the potential genetic nexus of multiple diseases in these categories. In this study, we were interested in detecting pleiotropy, or the phenomenon in which a genetic variant influences more than one phenotype. Marker-phenotype association approaches can be grouped into univariate, bivariate, and multivariate categories based on the number of phenotypes considered at one time. Here we applied one statistical method per category followed by an eQTL colocalization analysis to identify potential pleiotropic variants that contribute to the link between cardiovascular and neurological diseases. We performed our analyses on ~530,000 common SNPs coupled with 65 electronic health record (EHR)-based phenotypes in 43,870 unrelated European adults from the Electronic Medical Records and Genomics (eMERGE) network. There were 31 variants identified by all three methods that showed significant associations across late onset cardiac- and neurologic- diseases. We further investigated functional implications of gene expression on the detected "lead SNPs" via colocalization analysis, providing a deeper understanding of the discovered associations. In summary, we present the framework and landscape for detecting potential pleiotropy using univariate, bivariate, multivariate, and colocalization methods. Further exploration of these potentially pleiotropic genetic variants will work toward understanding disease causing mechanisms across cardiovascular and neurological diseases and may assist in considering disease prevention as well as drug repositioning in future research.

Genetic architecture underlying the lignin biosynthesis pathway involves noncoding RNAs and transcription factors for growth and wood properties in Populus

Lignin provides structural support in perennial woody plants and is a complex phenolic polymer derived from phenylpropanoid pathway. Lignin biosynthesis is regulated by coordinated networks involving transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). However, the genetic networks underlying the lignin biosynthesis pathway for tree growth and wood properties remain unknown. Here, we used association genetics (additive, dominant and epistasis) and expression quantitative trait nucleotide (eQTN) mapping to decipher the genetic networks for tree growth and wood properties in 435 unrelated individuals of Populus tomentosa. We detected 124 significant associations (P ≤ 6.89E-05) for 10 growth and wood property traits using 30 265 single nucleotide polymorphisms from 203 lignin biosynthetic genes, 81 TF genes, 36 miRNA genes and 71 lncRNA loci, implying their common roles in wood formation. Epistasis analysis uncovered 745 significant pairwise interactions, which helped to construct proposed genetic networks of lignin biosynthesis pathway and found that these regulators might affect phenotypes by linking two lignin biosynthetic genes. eQTNs were used to interpret how causal genes contributed to phenotypes. Lastly, we investigated the possible functions of the genes encoding 4-coumarate: CoA ligase and cinnamate-4-hydroxylase in wood traits using epistasis, eQTN mapping and enzymatic activity assays. Our study provides new insights into the lignin biosynthesis pathway in poplar and enables the novel genetic factors as biomarkers for facilitating genetic improvement of trees.

Biological clock function is linked to proactive and reactive personality types

Background

Many physiological processes in our body are controlled by the biological clock and show circadian rhythmicity. It is generally accepted that a robust rhythm is a prerequisite for optimal functioning and that a lack of rhythmicity can contribute to the pathogenesis of various diseases. Here, we tested in a heterogeneous laboratory zebrafish population whether and how variation in the rhythmicity of the biological clock is associated with the coping styles of individual animals, as assessed in a behavioural assay to reliably measure this along a continuum between proactive and reactive extremes.

Results

Using RNA sequencing on brain samples, we demonstrated a prominent difference in the expression level of genes involved in the biological clock between proactive and reactive individuals. Subsequently, we tested whether this correlation between gene expression and coping style was due to a consistent change in the level of clock gene expression or to a phase shift or to altered amplitude of the circadian rhythm of gene expression. Our data show a remarkable individual variation in amplitude of the clock gene expression rhythms, which was also reflected in the fluctuating concentrations of melatonin and cortisol, and locomotor activity. This variation in rhythmicity showed a strong correlation with the coping style of the individual, ranging from robust rhythms with large amplitudes in proactive fish to a complete absence of rhythmicity in reactive fish. The rhythmicity of the proactive fish decreased when challenged with constant light conditions whereas the rhythmicity of reactive individuals was not altered.

Conclusion

These results shed new light on the role of the biological clock by demonstrating that large variation in circadian rhythmicity of individuals may occur within populations. The observed correlation between coping style and circadian rhythmicity suggests that the level of rhythmicity forms an integral part of proactive or reactive coping styles.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0618-0) contains supplementary material, which is available to authorized users.

PCA-Based Multiple-Trait GWAS Analysis: A Powerful Model for Exploring Pleiotropy

Principal component analysis (PCA) is a potential approach that can be applied in multiple-trait genome-wide association studies (GWAS) to explore pleiotropy, as well as increase the power of quantitative trait loci (QTL) detection. In this study, the relationship of test single nucleotide polymorphisms (SNPs) was determined between single-trait GWAS and PCA-based GWAS. We found that the estimated pleiotropic quantitative trait nucleotides (QTNs) β * ^ were in most cases larger than the single-trait model estimations ( β 1 ^ and β 2 ^ ). Analysis using the simulated data showed that PCA-based multiple-trait GWAS has improved statistical power for detecting QTL compared to single-trait GWAS. For the minor allele frequency (MAF), when the MAF of QTNs was greater than 0.2, the PCA-based model had a significant advantage in detecting the pleiotropic QTNs, but when its MAF was reduced from 0.2 to 0, the advantage began to disappear. In addition, as the linkage disequilibrium (LD) of the pleiotropic QTNs decreased, its detection ability declined in the co-localization effect model. Furthermore, on the real data of 1141 Simmental cattle, we applied the PCA model to the multiple-trait GWAS analysis and identified a QTL that was consistent with a candidate gene, MCHR2, which was associated with presoma muscle development in cattle. In summary, PCA-based multiple-trait GWAS is an efficient model for exploring pleiotropic QTNs in quantitative traits.

Chromatin interactions and expression quantitative trait loci reveal genetic drivers of multimorbidities

Clinical studies of non-communicable diseases identify multimorbidities that suggest a common set of predisposing factors. Despite the fact that humans have ~24,000 genes, we do not understand the genetic pathways that contribute to the development of multimorbid non-communicable disease. Here we create a multimorbidity atlas of traits based on pleiotropy of spatially regulated genes. Using chromatin interaction and expression Quantitative Trait Loci (eQTL) data, we analyse 20,782 variants (p < 5 × 10-6) associated with 1351 phenotypes to identify 16,248 putative spatial eQTL-eGene pairs that are involved in 76,013 short- and long-range regulatory interactions (FDR < 0.05) in different human tissues. Convex biclustering of spatial eGenes that are shared among phenotypes identifies complex interrelationships between nominally different phenotype-associated SNPs. Our approach enables the simultaneous elucidation of variant interactions with target genes that are drivers of multimorbidity, and those that contribute to unique phenotype associated characteristics.

Polygenic Score for Alzheimer Disease and cognition: The mediating role of personality

Alzheimer's disease (AD) polygenic risk score (PGS) is associated with lower cognitive functioning even among older individuals without dementia. We tested the hypothesis that personality traits mediate the association between AD genetic risk and cognitive functioning. Participants (N > 7,000, aged 50-99 years old) from the Health and Retirement Study were genotyped and completed personality and cognition tests at baseline. Cognition was assessed again four years later. Bootstrap analysis revealed that a higher AD polygenic risk score was associated with lower cognitive scores at baseline through higher neuroticism, lower conscientiousness, and lower levels of the industriousness facet of conscientiousness. In addition, a higher polygenic score for AD was associated with decline in cognition over four years in part through higher neuroticism and lower conscientiousness. The findings support the hypothesis that the genetic vulnerability for AD contributes to cognitive functioning in part through its association with personality traits.

Diagnostics for Pleiotropy in Mendelian Randomization Studies: Global and Individual Tests for Direct Effects

Diagnosing pleiotropy is critical for assessing the validity of Mendelian randomization (MR) analyses. The popular MR-Egger method evaluates whether there is evidence of bias-generating pleiotropy among a set of candidate genetic instrumental variables. In this article, we propose a statistical method-global and individual tests for direct effects (GLIDE)-for systematically evaluating pleiotropy among the set of genetic variants (e.g., single nucleotide polymorphisms (SNPs)) used for MR. As a global test, simulation experiments suggest that GLIDE is nearly uniformly more powerful than the MR-Egger method. As a sensitivity analysis, GLIDE is capable of detecting outliers in individual variant-level pleiotropy, in order to obtain a refined set of genetic instrumental variables. We used GLIDE to analyze both body mass index and height for associations with colorectal cancer risk in data from the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry (multiple studies). Among the body mass index-associated SNPs and the height-associated SNPs, several individual variants showed evidence of pleiotropy. Removal of these potentially pleiotropic SNPs resulted in attenuation of respective estimates of the causal effects. In summary, the proposed GLIDE method is useful for sensitivity analyses and improves the validity of MR.

Revisiting the evolutionary origins of obesity: lazy versus peppy-thrifty genotype hypothesis

The recent global obesity epidemic is attributed to major societal and environmental changes, such as excessive energy intake and sedentary lifestyle. However, exposure to 'obesogenic' environments does not necessarily result in obesity at the individual level, as 40-75% of body mass index variation in population is attributed to genetic differences. The thrifty genotype theory posits that genetic variants promoting efficient food sequestering and optimal deposition of fat during periods of food abundance were evolutionarily advantageous for the early hunter-gatherer and were positively selected. However, the thrifty genotype is likely too simplistic and fails to provide a justification for the complex distribution of obesity predisposing gene variants and for the broad range of body mass index observed in diverse ethnic groups. This review proposes that gene pleiotropy may better account for the variability in the distribution of obesity susceptibility alleles across modern populations. We outline the lazy-thrifty versus peppy-thrifty genotype hypothesis and detail the body of evidence in the literature in support of this novel concept. Future population genetics and mathematical modelling studies that account for pleiotropy may further improve our understanding of the evolutionary origins of the current obesity epidemic.

Genetic and Environmental Overlaps Among Sasang Constitution Types: A Multivariate Twin Study

According to the Sasang theory, humans can be categorized into one of the four Sasang constitution (SC) types. The four SC types are Tae-Yang (TY), Tae-Eum (TE), So-Yang (SY), and So-Eum (SE), which are determined mainly on the basis of anthropometric characteristics, personality, and the balance of the physiological functions of the major organ systems. There is a growing recognition in the complementary and alternative medicine area that SC types have the potential to be a useful scientific tool for the prevention, diagnosis, and treatment of diseases (Cooper, Evidence Based Complementary and Alternative Medicine, Vol. 6 (Suppl. 1), 2009, pp. 1–3). The main purposes of the present study are to estimate genetic and environmental influences on SC types, and to explore genetic and environmental correlations that affect phenotypic associations among the SC types. In total, 1,742 (365 monozygotic male, 173 dizygotic male, 675 monozygotic female, 271 dizygotic female, and 258 opposite-sex dizygotic) twins (mean age = 19.1 ± 3.1 year) completed a Sasang questionnaire. Univariate and multivariate model-fitting analyses were performed. Total (additive and non-additive) genetic influences were 71% for males and 81% for females in TE, 70% for males and 71% for females in SE, and 47% for both sexes in SY. Non-additive genetic effects were substantial, and shared environmental influences were negligible in most SC types. Multivariate model-fitting analysis revealed that non-additive genetic and individual-specific environmental correlations between TE and SE were -0.92 (95% CI [-0.89, -0.93]) and -0.62 (95% CI [-0.57, -0.68]), respectively. The corresponding estimates were -0.55 (95% CI [-0.48, -0.61]) and -0.44 (95% CI [-0.37, -0.51]) between TE and SY and 0.19 (95% CI [0.09, 0.29]) and -0.40 (95% CI [-0.32, -0.47]) between SE and SY. These results suggest that the phenotypic associations among SC types may be mediated by pleiotropic mechanism of genes.

Phenome-wide association studies across large population cohorts support drug target validation

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

Genetic pleiotropy between mood disorders, metabolic, and endocrine traits in a multigenerational pedigree

Bipolar disorder (BD) is a mental disorder characterized by alternating periods of depression and mania. Individuals with BD have higher levels of early mortality than the general population, and a substantial proportion of this is due to increased risk for comorbid diseases. To identify the molecular events that underlie BD and related medical comorbidities, we generated imputed whole-genome sequence data using a population-specific reference panel for an extended multigenerational Old Order Amish pedigree (n = 394), segregating BD and related disorders. First, we investigated all putative disease-causing variants at known Mendelian disease loci present in this pedigree. Second, we performed genomic profiling using polygenic risk scores (PRS) to establish each individual's risk for several complex diseases. We identified a set of Mendelian variants that co-occur in individuals with BD more frequently than their unaffected family members, including the R3527Q mutation in APOB associated with hypercholesterolemia. Using PRS, we demonstrated that BD individuals from this pedigree were enriched for the same common risk alleles for BD as the general population (β = 0.416, p = 6 × 10^-4). Furthermore, we find evidence for a common genetic etiology between BD risk and polygenic risk for clinical autoimmune thyroid disease (p = 1 × 10^-4), diabetes (p = 1 × 10^-3), and lipid traits such as triglyceride levels (p = 3 × 10^-4) in the pedigree. We identify genomic regions that contribute to the differences between BD individuals and unaffected family members by calculating local genetic risk for independent LD blocks. Our findings provide evidence for the extensive genetic pleiotropy that can drive epidemiological findings of comorbidities between diseases and other complex traits.

A polygenic p factor for major psychiatric disorders

It has recently been proposed that a single dimension, called the p factor, can capture a person's liability to mental disorder. Relevant to the p hypothesis, recent genetic research has found surprisingly high genetic correlations between pairs of psychiatric disorders. Here, for the first time, we compare genetic correlations from different methods and examine their support for a genetic p factor. We tested the hypothesis of a genetic p factor by applying principal component analysis to matrices of genetic correlations between major psychiatric disorders estimated by three methods-family study, genome-wide complex trait analysis, and linkage-disequilibrium score regression-and on a matrix of polygenic score correlations constructed for each individual in a UK-representative sample of 7 026 unrelated individuals. All disorders loaded positively on a first unrotated principal component, which accounted for 57, 43, 35, and 22% of the variance respectively for the four methods. Our results showed that all four methods provided strong support for a genetic p factor that represents the pinnacle of the hierarchical genetic architecture of psychopathology.

Heritability informed power optimization (HIPO) leads to enhanced detection of genetic associations across multiple traits

Genome-wide association studies have shown that pleiotropy is a common phenomenon that can potentially be exploited for enhanced detection of susceptibility loci. We propose heritability informed power optimization (HIPO) for conducting powerful pleiotropic analysis using summary-level association statistics. We find optimal linear combinations of association coefficients across traits that are expected to maximize non-centrality parameter for the underlying test statistics, taking into account estimates of heritability, sample size variations and overlaps across the traits. Simulation studies show that the proposed method has correct type I error, robust to population stratification and leads to desired genome-wide enrichment of association signals. Application of the proposed method to publicly available data for three groups of genetically related traits, lipids (N = 188,577), psychiatric diseases (Ncase = 33,332, Ncontrol = 27,888) and social science traits (N ranging between 161,460 to 298,420 across individual traits) increased the number of genome-wide significant loci by 12%, 200% and 50%, respectively, compared to those found by analysis of individual traits. Evidence of replication is present for many of these loci in subsequent larger studies for individual traits. HIPO can potentially be extended to high-dimensional phenotypes as a way of dimension reduction to maximize power for subsequent genetic association testing.

Prediction of Schizophrenia Diagnosis by Integration of Genetically Correlated Conditions and Traits

Schizophrenia is genetically heterogeneous and comorbid with many conditions. In this study, we explored polygenic scores (PGSs) from genetically related conditions and traits to predict schizophrenia diagnosis using both logistic regression and deep neural network (DNN) models. We used the combined Molecular Genetics of Schizophrenia and Swedish Schizophrenia Case Control Study (MGS + SSCCS) data for training and testing the models, and used the Clinical Antipsychotic Trials for Intervention Effectiveness (CATIE) data as independent validation. We screened 28 conditions and traits comorbid with schizophrenia to identify traits as potential predictors and used LASSO regression to select predictors for model construction. We investigated how PGS calculation influenced model performance. We found that the inclusion of comorbid traits improved model performance and PGSs calculated from two traits were more generalizable in independent validation. With a DNN model using 19 PGS predictors, we accomplished a prediction accuracy of 0.813 and an AUC of 0.905 in the MGS + SSCCS data. When this model was validated with the CATIE data, it achieved an accuracy of 0.721 and AUC of 0.747. Our results indicate that PGSs alone may not be sufficient to predict schizophrenia accurately and the inclusion of behavioral and clinical data may be necessary for more accurate prediction model.

The use of electronic health records for psychiatric phenotyping and genomics

The widespread adoption of electronic health record (EHRs) in healthcare systems has created a vast and continuously growing resource of clinical data and provides new opportunities for population-based research. In particular, the linking of EHRs to biospecimens and genomic data in biobanks may help address what has become a rate-limiting study for genetic research: the need for large sample sizes. The principal roadblock to capitalizing on these resources is the need to establish the validity of phenotypes extracted from the EHR. For psychiatric genetic research, this represents a particular challenge given that diagnosis is based on patient reports and clinician observations that may not be well-captured in billing codes or narrative records. This review addresses the opportunities and pitfalls in EHR-based phenotyping with a focus on their application to psychiatric genetic research. A growing number of studies have demonstrated that diagnostic algorithms with high positive predictive value can be derived from EHRs, especially when structured data are supplemented by text mining approaches. Such algorithms enable semi-automated phenotyping for large-scale case-control studies. In addition, the scale and scope of EHR databases have been used successfully to identify phenotypic subgroups and derive algorithms for longitudinal risk prediction. EHR-based genomics are particularly well-suited to rapid look-up replication of putative risk genes, studies of pleiotropy (phenomewide association studies or PheWAS), investigations of genetic networks and overlap across the phenome, and pharmacogenomic research. EHR phenotyping has been relatively under-utilized in psychiatric genomic research but may become a key component of efforts to advance precision psychiatry.

How can genetics help understand the relationship between cognitive dysfunction and schizophrenia?

Despite the consistent finding that cognitive dysfunction is a core characteristic of schizophrenia (SCZ), little is known about the underlying pathophysiology. Recent progress in human genetics, driven by large genome-wide association studies (GWAS), has provided new data about the genetic architecture of complex human traits, including cognition and SCZ. Novel analytical tools have provided unprecedented opportunities to leverage the large amount of information from GWAS. Here we review the latest findings related to genetic architecture and risk genes of SCZ and cognitive functions, and recent findings of overlapping genetic factors. The recent GWAS of SCZ implicate over 100 risk gene loci, each with a small effect. A similar genetic architecture seems to be present in cognitive domains, suggesting that these phenotypes are highly polygenic. Further, GWAS have revealed more than 20 gene loci associated with cognitive traits, including intelligence, general cognition (g-factor), reaction time and verbal-numerical reasoning. Several gene loci have been implicated in educational attainment, a proxy measure of cognitive function. Recently, overlapping gene loci were found between education and SCZ, and between SCZ and cognitive traits, suggesting common genetic risk between SCZ and cognitive dysfunction. Mathematical modeling of GWAS of cognition and SCZ indicate that only a fraction of the heritability is identified. The evidence suggests a polygenic architecture for SCZ and cognitive functions, and a large degree of shared genetic risk. This indicates novel molecular genetic mechanisms and strengthens the notion that SCZ is more likely a part of the normal distribution and not a separate entity.

Adjustment for covariates using summary statistics of genome-wide association studies

Linear regression is a standard approach to identify genetic variants associated with continuous traits in genome-wide association studies (GWAS). In a standard epidemiology study, linear regression is often performed with adjustment for covariates to estimate the independent effect of a predictor variable or to improve statistical power by reducing residual variability. However, it is problematic to adjust for heritable covariates in genetic association analysis. Here, we propose a new method that utilizes summary statistics of the covariate from additional samples for reducing the residual variability and hence improves statistical power. Our simulation study showed that the proposed methodology can maintain a good control of Type I error and can achieve much higher power than a simple linear regression. The method is illustrated by an application to the GWAS results from the Genetic Investigation of Anthropometric Traits consortium.

Methods and results from the genome-wide association group at GAW20

BACKGROUND

This paper summarizes the contributions from the Genome-wide Association Study group (GWAS group) of the GAW20. The GWAS group contributions focused on topics such as association tests, phenotype imputation, and application of empirical kinships. The goals of the GWAS group contributions were varied. A real or a simulated data set based on the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study was employed by different methods. Different outcomes and covariates were considered, and quality control procedures varied throughout the contributions.

RESULTS

The consideration of heritability and family structure played a major role in some contributions. The inclusion of family information and adaptive weights based on data were found to improve power in genome-wide association studies. It was proven that gene-level approaches are more powerful than single-marker analysis. Other contributions focused on the comparison between pedigree-based kinship and empirical kinship matrices, and investigated similar results in heritability estimation, association mapping, and genomic prediction. A new approach for linkage mapping of triglyceride levels was able to identify a novel linkage signal.

CONCLUSIONS

This summary paper reports on promising statistical approaches and findings of the members of the GWAS group applied on real and simulated data which encompass the current topics of epigenetic and pharmacogenomics.

Exploring the Genetic Correlation Between Growth and Immunity Based on Summary Statistics of Genome-Wide Association Studies

The relationship between growth and immune phenotypes has been presented in the context of physiology and energy allocation theory, but has rarely been explained genetically in humans. As more summary statistics of genome-wide association studies (GWAS) become available, it is increasingly possible to explore the genetic relationship between traits at the level of genome-wide summary statistics. In this study, publicly available summary statistics of growth and immune related traits were used to evaluate the genetic correlation coefficients between immune and growth traits, as well as the cause and effect relationship between them. In addition, pleiotropic variants and KEGG pathways were identified. As a result, we found negative correlations between birthweight and immune cell count phenotypes, a positive correlation between childhood head circumference and eosinophil counts (EO), and positive or negative correlations between childhood body mass index and immune phenotypes. Statistically significant negative effects of immune cell count phenotypes on human height, and a slight but significant negative influence of human height on allergic disease were also observed. A total of 98 genomic regions were identified as containing variants potentially related to both immunity and growth. Some variants, such as rs3184504 located in SH2B3, rs13107325 in SLC39A8, and rs1260326 located in GCKR, which have been identified to be pleiotropic SNPs among other traits, were found to also be related to growth and immune traits in this study. Meanwhile, the most frequent overlapping KEGG pathways between growth and immune phenotypes were autoimmune related pathways. Pleiotropic pathways such as the adipocytokine signaling pathway and JAK-STAT signaling pathway were also identified to be significant. The results of this study indicate the complex genetic relationship between growth and immune phenotypes, and reveal the genetic background of their correlation in the context of pleiotropy.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

Genetic variants in microRNA biogenesis genes as novel indicators for secondary growth in Populus

MicroRNAs (miRNAs) function as key regulators of complex traits, but how genetic alterations in miRNA biogenesis genes (miRBGs) affect quantitative variation has not been elucidated. We conducted transcript analyses and association genetics to investigate how miRBGs, miRNA genes (MIRNAs) and their respective targets contribute to secondary growth in a natural population of 435 Populus tomentosa individuals. This analysis identified 29 843 common single-nucleotide polymorphisms (SNPs; frequency > 0.10) within 682 genes (80 miRBGs, 152 MIRNAs, and 457 miRNA targets). Single-SNP association analysis found SNPs in 234 candidate genes exhibited significant additive/dominant effects on phenotypes. Among these, specific candidates that associated with the same traits produced 791 miRBG-MIRNA-target combinations, suggesting possible genetic miRBG-MIRNA and MIRNA-target interactions, providing an important clue for the regulatory mechanisms of miRBGs. Multi-SNP association found 4672 epistatic pairs involving 578 genes that showed significant associations with traits and identified 106 miRBG-MIRNA-target combinations. Two multi-hierarchical networks were constructed based on correlations of miRBG-miRNA and miRNA-target expression to further probe the mechanisms of trait diversity underlying changes in miRBGs. Our study opens avenues for the investigation of miRNA function in perennial plants and underscored miRBGs as potentially modulating quantitative variation in traits.

Inferring Causal Relationships Between Risk Factors and Outcomes from Genome-Wide Association Study Data

An observational correlation between a suspected risk factor and an outcome does not necessarily imply that interventions on levels of the risk factor will have a causal impact on the outcome (correlation is not causation). If genetic variants associated with the risk factor are also associated with the outcome, then this increases the plausibility that the risk factor is a causal determinant of the outcome. However, if the genetic variants in the analysis do not have a specific biological link to the risk factor, then causal claims can be spurious. We review the Mendelian randomization paradigm for making causal inferences using genetic variants. We consider monogenic analysis, in which genetic variants are taken from a single gene region, and polygenic analysis, which includes variants from multiple regions. We focus on answering two questions: When can Mendelian randomization be used to make reliable causal inferences, and when can it be used to make relevant causal inferences?

The Emerging Immunogenetic Architecture of Schizophrenia

Schizophrenia is a severe psychiatric disorder of complex etiology. Immune processes have long been proposed to contribute to the development of schizophrenia, and accumulating evidence supports immune involvement in at least a subset of cases. In recent years, large-scale genetic studies have provided new insights into the role of the immune system in this disease. Here, we provide an overview of the immunogenetic architecture of schizophrenia based on findings from genome-wide association studies (GWAS). First, we review individual immune loci identified in secondary analyses of GWAS, which implicate over 30 genes expressed in both immune and brain cells. The function of the proteins encoded by these immune candidates highlight the role of the complement system, along with regulation of apoptosis in both immune and neuronal cells. Next, we review hypothesis-free pathway analyses which have so far been inconclusive with respect to identifying immune pathways involved in schizophrenia. Finally, we explore the genetic overlap between schizophrenia and immune-mediated diseases. Although there have been some inconsistencies across studies, genome-wide pleiotropy has been reported between schizophrenia and Crohn's disease, multiple sclerosis, rheumatoid arthritis, systemic lupus erythematosus, type 1 diabetes, and ulcerative colitis. Overall, there are multiple lines of evidence supporting the role of immune genes in schizophrenia. Current evidence suggests that specific immune pathways are involved-likely those with dual functions in the central nervous system. Future studies focused on further elucidating the relevant pathways hold the potential to identify novel biomarkers and therapeutic targets for schizophrenia.

An Atlas of Genetic Variation Linking Pathogen-Induced Cellular Traits to Human Disease

Pathogens have been a strong driving force for natural selection. Therefore, understanding how human genetic differences impact infection-related cellular traits can mechanistically link genetic variation to disease susceptibility. Here we report the Hi-HOST Phenome Project (H2P2): a catalog of cellular genome-wide association studies (GWAS) comprising 79 infection-related phenotypes in response to 8 pathogens in 528 lymphoblastoid cell lines. Seventeen loci surpass genome-wide significance for infection-associated phenotypes ranging from pathogen replication to cytokine production. We combined H2P2 with clinical association data from patients to identify a SNP near CXCL10 as a risk factor for inflammatory bowel disease. A SNP in the transcriptional repressor ZBTB20 demonstrated pleiotropy, likely through suppression of multiple target genes, and was associated with viral hepatitis. These data are available on a web portal to facilitate interpreting human genome variation through the lens of cell biology and should serve as a rich resource for the research community.

Genome-wide identification of directed gene networks using large-scale population genomics data

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10-10), among which transcription factors were overrepresented (Fisher's P = 3.3 × 10-7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser.

LLR: a latent low-rank approach to colocalizing genetic risk variants in multiple GWAS

Motivation

Genome-wide association studies (GWAS), which genotype millions of single nucleotide polymorphisms (SNPs) in thousands of individuals, are widely used to identify the risk SNPs underlying complex human phenotypes (quantitative traits or diseases). Most conventional statistical methods in GWAS only investigate one phenotype at a time. However, an increasing number of reports suggest the ubiquity of pleiotropy, i.e. many complex phenotypes sharing common genetic bases. This motivated us to leverage pleiotropy to develop new statistical approaches to joint analysis of multiple GWAS.

Results

In this study, we propose a latent low-rank (LLR) approach to colocalizing genetic risk variants using summary statistics. In the presence of pleiotropy, there exist risk loci that affect multiple phenotypes. To leverage pleiotropy, we introduce a low-rank structure to modulate the probabilities of the latent association statuses between loci and phenotypes. Regarding the computational efficiency of LLR, a novel expectation-maximization-path (EM-path) algorithm has been developed to greatly reduce the computational cost and facilitate model selection and inference. We demonstrate the advantages of LLR over competing approaches through simulation studies and joint analysis of 18 GWAS datasets.

Availability and implementation

The LLR software is available on https://sites.google.com/site/liujin810822.

Contact

macyang@ust.hk.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Trans-Ethnic Polygenic Analysis Supports Genetic Overlaps of Lumbar Disc Degeneration With Height, Body Mass Index, and Bone Mineral Density

Lumbar disc degeneration (LDD) is age-related break-down in the fibrocartilaginous joints between lumbar vertebrae. It is a major cause of low back pain and is conventionally assessed by magnetic resonance imaging (MRI). Like most other complex traits, LDD is likely polygenic and influenced by both genetic and environmental factors. However, genome-wide association studies (GWASs) of LDD have uncovered few susceptibility loci due to the limited sample size. Previous epidemiology studies of LDD also reported multiple heritable risk factors, including height, body mass index (BMI), bone mineral density (BMD), lipid levels, etc. Genetics can help elucidate causality between traits and suggest loci with pleiotropic effects. One such approach is polygenic score (PGS) which summarizes the effect of multiple variants by the summation of alleles weighted by estimated effects from GWAS. To investigate genetic overlaps of LDD and related heritable risk factors, we calculated the PGS of height, BMI, BMD and lipid levels in a Chinese population-based cohort with spine MRI examination and a Japanese case-control cohort of lumbar disc herniation (LDH) requiring surgery. Because most large-scale GWASs were done in European populations, PGS of corresponding traits were created using weights from European GWASs. We calibrated their prediction performance in independent Chinese samples, then tested associations with MRI-derived LDD scores and LDH affection status. The PGS of height, BMI, BMD and lipid levels were strongly associated with respective phenotypes in Chinese, but phenotype variances explained were lower than in Europeans which would reduce the power to detect genetic overlaps. Despite of this, the PGS of BMI and lumbar spine BMD were significantly associated with LDD scores; and the PGS of height was associated with the increased the liability of LDH. Furthermore, linkage disequilibrium score regression suggested that, osteoarthritis, another degenerative disorder that shares common features with LDD, also showed genetic correlations with height, BMI and BMD. The findings suggest a common key contribution of biomechanical stress to the pathogenesis of LDD and will direct the future search for pleiotropic genes.

Education, Smoking, and Cohort Change: Forwarding a Multidimensional Theory of the Environmental Moderation of Genetic Effects

We introduce a genetic correlation by environment interaction model [(rG)xE] which allows for social environmental moderation of the genetic relationship between two traits. To empirically demonstrate the significance of the (rG)xE perspective, we use genome wide information from respondents of the Health and Retirement Study (HRS; n = 8,181; birth years 1920-1959) and the National Longitudinal Study of Adolescent to Adult Health (Add Health; n = 4,347; birth years 1974-1983) to examine whether the genetic correlation (rG) between education and smoking has increased over historical time. Genetic correlation estimates (rG_HRS = -0.357; rG_{Add Health} = -0.729) support this hypothesis. Using polygenic scores for educational attainment, we show that this is not due to latent indicators of intellectual capacity, and we highlight the importance of education itself as an explanation of the increasing genetic correlation. Analyses based on contextual variation the milieus of the Add Health respondents corroborate key elements of the birth cohort analyses. We argue that the increasing overlap with respect to genes associated with educational attainment and smoking is a fundamentally social process involving complex process of selection based on observable behaviors that may be linked to genotype.

Pleiotropy-robust Mendelian randomization

Background

The potential of Mendelian randomization studies is rapidly expanding due to: (i) the growing power of genome-wide association study (GWAS) meta-analyses to detect genetic variants associated with several exposures; and (ii) the increasing availability of these genetic variants in large-scale surveys. However, without a proper biological understanding of the pleiotropic working of genetic variants, a fundamental assumption of Mendelian randomization (the exclusion restriction) can always be contested.

Methods

We build upon and synthesize recent advances in the literature on instrumental variables (IVs) estimation that test and relax the exclusion restriction. Our pleiotropy-robust Mendelian randomization (PRMR) method first estimates the degree of pleiotropy, and in turn corrects for it. If (i) a subsample exists for which the genetic variants do not affect the exposure; (ii) the selection into this subsample is not a joint consequence of the IV and the outcome; (iii) pleiotropic effects are homogeneous, PRMR obtains unbiased estimates of causal effects.

Results

Simulations show that existing MR methods produce biased estimators for realistic forms of pleiotropy. Under the aforementioned assumptions, PRMR produces unbiased estimators. We illustrate the practical use of PRMR by estimating the causal effect of: (i) tobacco exposure on body mass index (BMI); (ii) prostate cancer on self-reported health; and (iii) educational attainment on BMI in the UK Biobank data.

Conclusions

PRMR allows for instrumental variables that violate the exclusion restriction due to pleiotropy, and it corrects for pleiotropy in the estimation of the causal effect. If the degree of pleiotropy is unknown, PRMR can still be used as a sensitivity analysis.

High-sensitivity C-reactive protein, low-grade systemic inflammation and type 2 diabetes mellitus: A two-sample Mendelian randomization study

BACKGROUND AND AIMS

The role of inflammation in type 2 diabetes mellitus (T2D) remains unclear. We investigated the associations of high sensitivity C-reactive protein (hsCRP) concentration with T2D and glycemic traits using two-sample Mendelian Randomization.

METHODS AND RESULTS

We used publically available summary-statistics data from genome-wide association studies on T2D (DIAGRAM: 12 171 cases; 56 862 controls) and glycemic traits (MAGIC: 46 186 participants without diabetes mellitus). We combined the effects of the genetic instrumental variables through inverse-variance weighting (IVW), and MR-Egger regression and weighted-median estimation as sensitivity analyses which take into account potential violations (e.g., directional pleiotropy) of the assumptions of instrumental variable analyses. Analyses were conducted using 15 known hsCRP genetic instruments among which 6 instruments are hsCRP specific and not involved in inflammatory processes beyond hsCRP concentration regulation. Though we found no association between the combined effect of the genetic instrumental variables for hsCRP and T2D with IVW (odds ratio per 1 ln [hsCRP in mg/L]: 1.15; 95% confidence interval: 0.93, 1.42), we found associations for T2D with MR-Egger regression and weighted-median estimation (odds ratio with 95% confidence interval per 1 ln [hsCRP in mg/L], MR-Egger regression: 1.29; 1.08, 1.49; weighted-median estimator: 1.21; 1.02, 1.39). We found no association with T2D for the combination of hsCRP-specific genetic instruments nor did we found associations with glycemic traits in any of the analyses.

CONCLUSION

Evidence was provided for a potential causal association between hsCRP and T2D, but only after considering directional pleiotropy. However, hsCRP was not causally associated with glycemic traits.

Evaluating the potential role of pleiotropy in Mendelian randomization studies

Pleiotropy, the phenomenon of a single genetic variant influencing multiple traits, is likely widespread in the human genome. If pleiotropy arises because the single nucleotide polymorphism (SNP) influences one trait, which in turn influences another ('vertical pleiotropy'), then Mendelian randomization (MR) can be used to estimate the causal influence between the traits. Of prime focus among the many limitations to MR is the unprovable assumption that apparent pleiotropic associations are mediated by the exposure (i.e. reflect vertical pleiotropy), and do not arise due to SNPs influencing the two traits through independent pathways ('horizontal pleiotropy'). The burgeoning treasure trove of genetic associations yielded through genome wide association studies makes for a tantalizing prospect of phenome-wide causal inference. Recent years have seen substantial attention devoted to the problem of horizontal pleiotropy, and in this review we outline how newly developed methods can be used together to improve the reliability of MR.

Testing an optimally weighted combination of common and/or rare variants with multiple traits

Recently, joint analysis of multiple traits has become popular because it can increase statistical power to identify genetic variants associated with complex diseases. In addition, there is increasing evidence indicating that pleiotropy is a widespread phenomenon in complex diseases. Currently, most of existing methods test the association between multiple traits and a single genetic variant. However, these methods by analyzing one variant at a time may not be ideal for rare variant association studies because of the allelic heterogeneity as well as the extreme rarity of rare variants. In this article, we developed a statistical method by testing an optimally weighted combination of variants with multiple traits (TOWmuT) to test the association between multiple traits and a weighted combination of variants (rare and/or common) in a genomic region. TOWmuT is robust to the directions of effects of causal variants and is applicable to different types of traits. Using extensive simulation studies, we compared the performance of TOWmuT with the following five existing methods: gene association with multiple traits (GAMuT), multiple sequence kernel association test (MSKAT), adaptive weighting reverse regression (AWRR), single-TOW, and MANOVA. Our results showed that, in all of the simulation scenarios, TOWmuT has correct type I error rates and is consistently more powerful than the other five tests. We also illustrated the usefulness of TOWmuT by analyzing a whole-genome genotyping data from a lung function study.