Genes of the Unfolded Protein Response Pathway Harbor Risk Alleles for Primary Open Angle Glaucoma

The statistical power of genome-wide association (GWA) studies to detect risk alleles for human diseases is limited by the unfavorable ratio of SNPs to study subjects. This multiple testing problem can be surmounted with very large population sizes when common alleles of large effects give rise to disease status. However, GWA approaches fall short when many rare alleles may give rise to a common disease, or when the number of subjects that can be recruited is limited. Here, we demonstrate that this multiple testing problem can be overcome by a comparative genomics approach in which an initial genome-wide screen in a genetically amenable model organism is used to identify human orthologues that may harbor risk alleles for adult-onset primary open angle glaucoma (POAG). Glaucoma is a major cause of blindness, which affects over 60 million people worldwide. Several genes have been associated with juvenile onset glaucoma, but genetic factors that predispose to adult onset primary open angle glaucoma (POAG) remain largely unknown. Previous genome-wide analysis in a Drosophila ocular hypertension model identified transcripts with altered regulation and showed induction of the unfolded protein response (UPR) upon overexpression of transgenic human glaucoma-associated myocilin (MYOC). We selected 16 orthologous genes with 62 polymorphic markers and identified in two independent human populations two genes of the UPR that harbor POAG risk alleles, BIRC6 and PDIA5. Thus, effectiveness of the UPR in response to accumulation of misfolded or aggregated proteins may contribute to the pathogenesis of POAG and provide targets for early therapeutic intervention

Genetic basis of transcriptome diversity in Drosophila melanogaster

Understanding how DNA sequence variation is translated into variation for complex phenotypes has remained elusive but is essential for predicting adaptive evolution, for selecting agriculturally important animals and crops, and for personalized medicine. Gene expression may provide a link between variation in DNA sequence and organismal phenotypes, and its abundance can be measured efficiently and accurately. Here we quantified genomewide variation in gene expression in the sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP), increasing the annotated Drosophila transcriptome by 11%, including thousands of novel transcribed regions (NTRs). We found that 42%of the Drosophila transcriptome is genetically variable in males and females, including the NTRs, and is organized into modules of genetically correlated transcripts. We found that NTRs often were negatively correlated with the expression of protein-coding genes, which we exploited to annotate NTRs functionally. We identified regulatory variants for the mean and variance of gene expression, which have largely independent genetic control. Expression quantitative trait loci (eQTLs) for the mean, but not for the variance, of gene expression were concentrated near genes. Notably, the variance eQTLs often interacted epistatically with local variants in these genes to regulate gene expression. This comprehensive characterization of population-scale diversity of transcriptomes and its genetic basis in the DGRP is critically important for a systems understanding of quantitative trait variation

TSPO interacts with VDAC1 and triggers a ROS-mediated inhibition of mitochondrial quality control

The 18-kDa TSPO (translocator protein) localizes on the outer mitochondrial membrane (OMM) and participates in cholesterol transport. Here, we report that TSPO inhibits mitochondrial autophagy downstream of the PINK1-PARK2 pathway, preventing essential ubiquitination of proteins. TSPO abolishes mitochondrial relocation of SQSTM1/p62 (sequestosome 1), and consequently that of the autophagic marker LC3 (microtubule-associated protein 1 light chain 3), thus leading to an accumulation of dysfunctional mitochondria, altering the appearance of the network. Independent of cholesterol regulation, the modulation of mitophagy by TSPO is instead dependent on VDAC1 (voltage-dependent anion channel 1), to which TSPO binds, reducing mitochondrial coupling and promoting an overproduction of reactive oxygen species (ROS) that counteracts PARK2-mediated ubiquitination of proteins. These data identify TSPO as a novel element in the regulation of mitochondrial quality control by autophagy, and demonstrate the importance for cell homeostasis of its expression ratio with VDAC1

Obp56h Modulates Mating Behavior in Drosophila melanogaster

Social interactions in insects are driven by conspecific chemical signals that are detected via olfactory and gustatory neurons. Odorant binding proteins (Obps) transport volatile odorants to chemosensory receptors, but their effects on behaviors remain poorly characterized. Here, we report that RNAi knockdown of Obp56h gene expression in Drosophila melanogaster enhances mating behavior by reducing courtship latency. The change in mating behavior that results from inhibition of Obp56h expression is accompanied by significant alterations in cuticular hydrocarbon (CHC) composition, including reduction in 5-tricosene (5-T), an inhibitory sex pheromone produced by males that increases copulation latency during courtship. Whole genome RNA sequencing confirms that expression of Obp56h is virtually abolished in Drosophila heads. Inhibition of Obp56h expression also affects expression of other chemoreception genes, including upregulation of lush in both sexes and Obp83ef in females, and reduction in expression of Obp19b and Or19b in males. In addition, several genes associated with lipid metabolism, which underlies the production of cuticular hydrocarbons, show altered transcript abundances. Our data show that modulation of mating behavior through reduction of Obp56h is accompanied by altered cuticular hydrocarbon profiles and implicate 5-T as a possible ligand for Obp56h

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

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

Epigenetic Regulation of Learning and Memory by Drosophila EHMT/G9a

Behavioral phenotyping and genome-wide profiling of the histone modifier EHMT in Drosophila reveals a mechanism through which an epigenetic writer may control cognition

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

Genome-Wide Association Analysis of Oxidative Stress Resistance in Drosophila melanogaster

Background: Aerobic organisms are susceptible to damage by reactive oxygen species. Oxidative stress resistance is a quantitative trait with population variation attributable to the interplay between genetic and environmental factors. Drosophila melanogaster provides an ideal system to study the genetics of variation for resistance to oxidative stress. Methods and Findings: We used 167 wild-derived inbred lines of the Drosophila Genetic Reference Panel for a genomewide association study of acute oxidative stress resistance to two oxidizing agents, paraquat and menadione sodium bisulfite. We found significant genetic variation for both stressors. Single nucleotide polymorphisms (SNPs) associated with variation in oxidative stress resistance were often sex-specific and agent-dependent, with a small subset common for both sexes or treatments. Associated SNPs had moderately large effects, with an inverse relationship between effect size and allele frequency. Linear models with up to 12 SNPs explained 67–79 % and 56–66 % of the phenotypic variance for resistance to paraquat and menadione sodium bisulfite, respectively. Many genes implicated were novel with no known role in oxidative stress resistance. Bioinformatics analyses revealed a cellular network comprising DNA metabolism and neuronal development, consistent with targets of oxidative stress-inducing agents. We confirmed associations of seven candidate genes associated with natural variation in oxidative stress resistance through mutational analysis. Conclusions: We identified novel candidate genes associated with variation in resistance to oxidative stress that hav