Using Whole-Genome Sequence Data to Predict Quantitative Trait Phenotypes in Drosophila melanogaster

Predicting organismal phenotypes from genotype data is important for plant and animal breeding, medicine, and evolutionary biology. Genomic-based phenotype prediction has been applied for single-nucleotide polymorphism (SNP) genotyping platforms, but not using complete genome sequences. Here, we report genomic prediction for starvation stress resistance and startle response in Drosophila melanogaster, using ∼2.5 million SNPs determined by sequencing the Drosophila Genetic Reference Panel population of inbred lines. We constructed a genomic relationship matrix from the SNP data and used it in a genomic best linear unbiased prediction (GBLUP) model. We assessed predictive ability as the correlation between predicted genetic values and observed phenotypes by cross-validation, and found a predictive ability of 0.239±0.008 (0.230±0.012) for starvation resistance (startle response). The predictive ability of BayesB, a Bayesian method with internal SNP selection, was not greater than GBLUP. Selection of the 5% SNPs with either the highest absolute effect or variance explained did not improve predictive ability. Predictive ability decreased only when fewer than 150,000 SNPs were used to construct the genomic relationship matrix. We hypothesize that predictive power in this population stems from the SNP–based modeling of the subtle relationship structure caused by long-range linkage disequilibrium and not from population structure or SNPs in linkage disequilibrium with causal variants. We discuss the implications of these results for genomic prediction in other organisms

Epistasis dominates the genetic architecture of Drosophila quantitative traits

Epistasis-nonlinear genetic interactions between polymorphic loci-is the genetic basis of canalization and speciation, and epistatic interactions can be used to infer genetic networks affecting quantitative traits. However, the role that epistasis plays in the genetic architecture of quantitative traits is controversial. Here, we compared the genetic architecture of three Drosophila life history traits in the sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and a large outbred, advanced intercross population derived from 40 DGRP lines (Flyland). We assessed allele frequency changes between pools of individuals at the extremes of the distribution for each trait in the Flyland population by deep DNA sequencing. The genetic architecture of all traits was highly polygenic in both analyses. Surprisingly, none of the SNPs associated with the traits in Flyland replicated in the DGRP and vice versa. However, the majority of these SNPs participated in at least one epistatic interaction in the DGRP. Despite apparent additive effects at largely distinct loci in the two populations, the epistatic interactions perturbed common, biologically plausible, and highly connected genetic networks. Our analysis underscores the importance of epistasis as a principal factor that determines variation for quantitative traits and provides a means to uncover genetic networks affecting these traits. Knowledge of epistatic networks will contribute to our understanding of the genetic basis of evolutionarily and clinically important traits and enhance predictive ability at an individualized level in medicine and agricultur

Association between Sex-Biased Gene Expression and Mutations with Sex-Specific Phenotypic Consequences in Drosophila

Genome-wide mRNA transcription profiles reveal widespread molecular sexual dimorphism or “sex-biased” gene expression, yet the relationship between molecular and phenotypic sexual dimorphism remains unclear. A major unresolved question is whether sex-biased genes typically perform male- and female-specific functions (whether these genes have sex-biased phenotypic or fitness consequences) or have similar functional importance for both sexes. To elucidate the relationship between sex-biased transcription and sex-biased fitness consequences, we analyzed a large data set of lethal, visible, and sterile mutations that have been mapped to the Drosophila melanogaster genome. The data permitted us to classify genes according to their sex-specific mutational effects and to infer the relationship between sex-biased transcription level and sex-specific fitness consequences. We find that mutations in female-biased genes are (on average) more deleterious to females than to males and that mutations in male-biased genes tend to be more deleterious to males than to females. Nevertheless, mutations in most sex-biased genes have similar phenotypic consequences for both sexes, which suggests that sex-biased transcription is not necessarily associated with functional genetic differentiation between males and females. These results have interesting implications for the evolution of sexual dimorphism and sex-specific adaptation

Detection of regulator genes and eQTLs in gene networks

Genetic differences between individuals associated to quantitative phenotypic traits, including disease states, are usually found in non-coding genomic regions. These genetic variants are often also associated to differences in expression levels of nearby genes (they are "expression quantitative trait loci" or eQTLs for short) and presumably play a gene regulatory role, affecting the status of molecular networks of interacting genes, proteins and metabolites. Computational systems biology approaches to reconstruct causal gene networks from large-scale omics data have therefore become essential to understand the structure of networks controlled by eQTLs together with other regulatory genes, and to generate detailed hypotheses about the molecular mechanisms that lead from genotype to phenotype. Here we review the main analytical methods and softwares to identify eQTLs and their associated genes, to reconstruct co-expression networks and modules, to reconstruct causal Bayesian gene and module networks, and to validate predicted networks in silico.Comment: minor revision with typos corrected; review article; 24 pages, 2 figure

Phenotypic Plasticity of the Drosophila Transcriptome

Phenotypic plasticity is the ability of a single genotype to produce different phenotypes in response to changing environments. We assessed variation in genome-wide gene expression and four fitness-related phenotypes of an outbred Drosophila melanogaster population under 20 different physiological, social, nutritional, chemical, and physical environments; and we compared the phenotypically plastic transcripts to genetically variable transcripts in a single environment. The environmentally sensitive transcriptome consists of two transcript categories, which comprise ∼15% of expressed transcripts. Class I transcripts are genetically variable and associated with detoxification, metabolism, proteolysis, heat shock proteins, and transcriptional regulation. Class II transcripts have low genetic variance and show sexually dimorphic expression enriched for reproductive functions. Clustering analysis of Class I transcripts reveals a fragmented modular organization and distinct environmentally responsive transcriptional signatures for the four fitness-related traits. Our analysis suggests that a restricted environmentally responsive segment of the transcriptome preserves the balance between phenotypic plasticity and environmental canalization

Enteric Microbiome Metabolites Correlate with Response to Simvastatin Treatment

Although statins are widely prescribed medications, there remains considerable variability in therapeutic response. Genetics can explain only part of this variability. Metabolomics is a global biochemical approach that provides powerful tools for mapping pathways implicated in disease and in response to treatment. Metabolomics captures net interactions between genome, microbiome and the environment. In this study, we used a targeted GC-MS metabolomics platform to measure a panel of metabolites within cholesterol synthesis, dietary sterol absorption, and bile acid formation to determine metabolite signatures that may predict variation in statin LDL-C lowering efficacy. Measurements were performed in two subsets of the total study population in the Cholesterol and Pharmacogenetics (CAP) study: Full Range of Response (FR), and Good and Poor Responders (GPR) were 100 individuals randomly selected from across the entire range of LDL-C responses in CAP. GPR were 48 individuals, 24 each from the top and bottom 10% of the LDL-C response distribution matched for body mass index, race, and gender. We identified three secondary, bacterial-derived bile acids that contribute to predicting the magnitude of statin-induced LDL-C lowering in good responders. Bile acids and statins share transporters in the liver and intestine; we observed that increased plasma concentration of simvastatin positively correlates with higher levels of several secondary bile acids. Genetic analysis of these subjects identified associations between levels of seven bile acids and a single nucleotide polymorphism (SNP), rs4149056, in the gene encoding the organic anion transporter SLCO1B1. These findings, along with recently published results that the gut microbiome plays an important role in cardiovascular disease, indicate that interactions between genome, gut microbiome and environmental influences should be considered in the study and management of cardiovascular disease. Metabolic profiles could provide valuable information about treatment outcomes and could contribute to a more personalized approach to therapy

Expression QTL Modules as Functional Components Underlying Higher-Order Phenotypes

Systems genetics studies often involve the mapping of numerous regulatory relations between genetic loci and expression traits. These regulatory relations form a bipartite network consisting of genetic loci and expression phenotypes. Modular network organizations may arise from the pleiotropic and polygenic regulation of gene expression. Here we analyzed the expression QTL (eQTL) networks derived from expression genetic data of yeast and mouse liver and found 65 and 98 modules respectively. Computer simulation result showed that such modules rarely occurred in randomized networks with the same number of nodes and edges and same degree distribution. We also found significant within-module functional coherence. The analysis of genetic overlaps and the evidences from biomedical literature have linked some eQTL modules to physiological phenotypes. Functional coherence within the eQTL modules and genetic overlaps between the modules and physiological phenotypes suggests that eQTL modules may act as functional units underlying the higher-order phenotypes

Overexpression of Myocilin in the Drosophila Eye Activates the Unfolded Protein Response: Implications for Glaucoma

Glaucoma is the world's second leading cause of bilateral blindness with progressive loss of vision due to retinal ganglion cell death. Myocilin has been associated with congenital glaucoma and 2-4% of primary open angle glaucoma (POAG) cases, but the pathogenic mechanisms remain largely unknown. Among several hypotheses, activation of the unfolded protein response (UPR) has emerged as a possible disease mechanism.We used a transgenic Drosophila model to analyze whole-genome transcriptional profiles in flies that express human wild-type or mutant MYOC in their eyes. The transgenic flies display ocular fluid discharge, reflecting ocular hypertension, and a progressive decline in their behavioral responses to light. Transcriptional analysis shows that genes associated with the UPR, ubiquitination, and proteolysis, as well as metabolism of reactive oxygen species and photoreceptor activity undergo altered transcriptional regulation. Following up on the results from these transcriptional analyses, we used immunoblots to demonstrate the formation of MYOC aggregates and showed that the formation of such aggregates leads to induction of the UPR, as evident from activation of the fluorescent UPR marker, xbp1-EGFP. CONCLUSIONS / SIGNIFICANCE: Our results show that aggregation of MYOC in the endoplasmic reticulum activates the UPR, an evolutionarily conserved stress pathway that culminates in apoptosis. We infer from the Drosophila model that MYOC-associated ocular hypertension in the human eye may result from aggregation of MYOC and induction of the UPR in trabecular meshwork cells. This process could occur at a late age with wild-type MYOC, but might be accelerated by MYOC mutants to account for juvenile onset glaucoma

Making Informed Choices about Microarray Data Analysis

This article describes the typical stages in the analysis of microarray data for non-specialist researchers in systems biology and medicine. Particular attention is paid to significant data analysis issues that are commonly encountered among practitioners, some of which need wider airing. The issues addressed include experimental design, quality assessment, normalization, and summarization of multiple-probe data. This article is based on the ISMB 2008 tutorial on microarray data analysis. An expanded version of the material in this article and the slides from the tutorial can be found at http://www.people.vcu.edu/~mreimers/OGMDA/index.html