Search

Joint analysis of multiple correlated phenotypes for genome-wide association studies (GWAS) can identify and interpret pleiotropic loci which are essential to understand pleiotropy in diseases and complex traits. Meanwhile, constructing a network based on associations between phenotypes and genotypes provides a new insight to analyze multiple phenotypes, which can explore whether phenotypes and genotypes might be related to each other at a higher level of cellular and organismal organization. In this paper, we first develop a bipartite signed network by linking phenotypes and genotypes into a Genotype and Phenotype Network (GPN). The GPN can be constructed by a mixture of quantitative and qualitative phenotypes and is applicable to binary phenotypes with extremely unbalanced case-control ratios in large-scale biobank datasets. We then apply a powerful community detection method to partition phenotypes into disjoint network modules based on GPN. Finally, we jointly test the association between multiple phenotypes in a network module and a single nucleotide polymorphism (SNP). Simulations and analyses of 72 complex traits in the UK Biobank show that multiple phenotype association tests based on network modules detected by GPN are much more powerful than those without considering network modules. The newly proposed GPN provides a new insight to investigate the genetic architecture among different types of phenotypes. Multiple phenotypes association studies based on GPN are improved by incorporating the genetic information into the phenotype clustering. Notably, it might broaden the understanding of genetic architecture that exists between diagnoses, genes, and pleiotropy. Author summary: Biological pleiotropy refers to a SNP or gene that has a direct biological influence on more than one phenotypic trait, which can offer significant insights in understanding the complex genotype-phenotype relationships. Network analyses provide an integrative approach to characterize complex genomic associations by linking phenotypes and genotypes into a Genotype and Phenotype Network (GPN). Jointly analyzing multiple phenotypes and incorporating the genetic information into the phenotype clustering may increase the statistical power to discover the cross-phenotype association and pleiotropy. We evaluate our proposed multiple phenotype association tests based on network modules detected by GPN for 72 EHR-derived phenotypes in the diseases of the musculoskeletal system and connective tissue in the UK Biobank. From the post-GWAS analyses, we observe that the test based on GPN can identify more significantly enriched biological pathways than that without considering the network modules. Meanwhile, some of the uniquely identified SNPs by the test based on GPN are also colocalized in the eQTL study of the gene expression in the Muscle Skeletal tissue. [ABSTRACT FROM AUTHOR]

Transcriptional cis-regulatory modules, e.g., enhancers, control the time and location of metazoan gene expression. While changes in enhancers can provide a powerful force for evolution, there is also significant deep conservation of enhancers for developmentally important genes, with function and sequence characteristics maintained over hundreds of millions of years of divergence. Not well understood, however, is how the overall regulatory composition of a locus evolves, with important outstanding questions such as how many enhancers are conserved vs. novel, and to what extent are the locations of conserved enhancers within a locus maintained? We begin here to address these questions with a comparison of the respective single-minded (sim) loci in the two dipteran species Drosophila melanogaster (fruit fly) and Aedes aegypti (mosquito). sim encodes a highly conserved transcription factor that mediates development of the arthropod embryonic ventral midline. We identify two enhancers in the A. aegypti sim locus and demonstrate that they function equivalently in both transgenic flies and transgenic mosquitoes. One A. aegypti enhancer is highly similar to known Drosophila counterparts in its activity, location, and autoregulatory capability. The other differs from any known Drosophila sim enhancers with a novel location, failure to autoregulate, and regulation of expression in a unique subset of midline cells. Our results suggest that the conserved pattern of sim expression in the two species is the result of both conserved and novel regulatory sequences. Further examination of this locus will help to illuminate how the overall regulatory landscape of a conserved developmental gene evolves. Author summary: The expression patterns and roles of genes, especially those involved in core developmental processes, are often conserved over vast evolutionary distances. Paradoxically, the DNA sequences surrounding these genes, which contain the cis-regulatory sequences (enhancers) that regulate gene expression, tend to be highly diverged. The manner and extent to which enhancers are functionally conserved, and how the overall organization of regulatory sequences within a locus is preserved or restructured, is not well understood. In this paper, we investigate these questions by identifying enhancers controlling expression of a master nervous system regulatory gene named sim in the mosquito Aedes aegypti, and comparing their functions and locations to those in the well-characterized sim locus of the fruit fly Drosophila melanogaster. Our results suggest that the two species generate identical patterns of sim expression through a mix of conserved and novel regulatory sequences. Continued exploration of the sim locus in these two species will help to build a comprehensive picture of how a regulatory locus for a master developmental regulator has evolved. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWAS) have identified many genetic loci associated with complex traits and diseases in the past 20 years. Multiple heritable covariates may be added into GWAS regression models to estimate direct effects of genetic variants on a focal trait, or to improve the power by accounting for environmental effects and other sources of trait variations. When one or more covariates are causally affected by both genetic variants and hidden confounders, adjusting for them in GWAS will produce biased estimation of SNP effects, known as collider bias. Several approaches have been developed to correct collider bias through estimating the bias by Mendelian randomization (MR). However, these methods work for only one covariate, some of which utilize MR methods with relatively strong assumptions, both of which may not hold in practice. In this paper, we extend the bias-correction approaches in two aspects: first we derive an analytical expression for the collider bias in the presence of multiple covariates, then we propose estimating the bias using a robust multivariable MR (MVMR) method based on constrained maximum likelihood (called MVMR-cML), allowing the presence of invalid instrumental variables (IVs) and correlated pleiotropy. We also established the estimation consistency and asymptotic normality of the new bias-corrected estimator. We conducted simulations to show that all methods mitigated collider bias under various scenarios. In real data analyses, we applied the methods to two GWAS examples, the first a GWAS of waist-hip ratio with adjustment for only one covariate, body-mass index (BMI), and the second a GWAS of BMI adjusting metabolomic principle components as multiple covariates, illustrating the effectiveness of bias correction. Author summary: Genome-wide association studies (GWAS) are powerful in identifying genetic variants influencing complex traits and diseases. However, adjusting for heritable covariates in GWAS may introduce collider bias when both genetic variants and confounders may causally influence these covariates. In this study, for the first time we derived the analytical form of the bias term in GWAS with multiple covariates, enabling bias estimation and correction using any MVMR method. On the other hand, many existing MVMR methods may not be robust to invalid IVs and are designed for independent samples. Since GWAS data of multiple traits are needed, overlapping samples become inevitable. Hence, while investigating the performance of many MVMR methods, we mainly adopt MVMR-cML, a novel MVMR approach robust to invalid IVs and sample overlap. Our simulations underscore that most MVMR methods effectively reduce collider bias across various scenarios. Furthermore, by accounting for correlations among GWAS statistics, as well as the linkage disequilibrium (LD) between the target SNP and IVs, we establish the consistency and asymptotic normality of the bias-corrected estimator based on MVMR-cML. The application of our bias-correction approach to two published GWAS data examples illustrates its utility and efficacy. [ABSTRACT FROM AUTHOR]