Genetical toxicogenomics in Drosophila identifies master-modulatory loci that are regulated by developmental exposure to lead

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper quickly results in cell and tissue damage that can range in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes also respond to other non-essential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the response to heavy metal stress. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource (DSPR) using a combination of differential expression analysis and expression quantitative trait locus (eQTL) mapping. Differential expression analysis revealed clear patterns of tissue-specific expression, primarily driven by a more pronounced gene expression response in gut tissue. eQTL mapping of gene expression under control and copper conditions as well as for the change in gene expression following copper exposure (copper response eQTL) revealed hundreds of genes with tissue-specific local cis-eQTL and many distant trans-eQTL. eQTL associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited genotype by environment effects on gene expression under copper stress, illuminating several tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight many candidate genes for future functional characterization.

10 Years of Toxicogenomics section in Frontiers in Genetics: Past discoveries and Future Perspectives

The Frontiers Media family has over 200 journals, which are each headed by usually one Field Chief Editor, and several thousand specialty sections, which are each headed by one or more Specialty Chief Editors. The year 2021 was the 10th anniversary of the founding of the Frontiers in Genetics journal and the Frontiers in Toxicogenomics specialty section of this journal. In 2021, we also announce one of the newest of the Frontiers journals-Frontiers in Toxicology which is part of the Frontiers Media family of journals but independent of Frontiers in Genetics. Dr. Ruden is the founding, and currently sole, Specialty Chief Editor of Frontiers in Toxicogenomics and one of 9 Specialty Chief Editors of Frontiers in Toxicology. As of 2021, Frontiers in Toxicogenomics has published over 138 articles and has over 370 Editors including 90 Associate Editors and 280 Review Editors. The Frontiers in Genetics impact factor was initially approximately 2.5 when it was first listed in PubMed in 2015 and has risen steadily to its current value of 4.8, which is typical for the majority of the over 200 Frontiers journals that have established impact factors. In this overview of the first decade of Frontiers in Toxicogenomics, we discuss the top 5 articles with the highest Scopus citations, which were all written in the first few years of the journal. The article with the highest number of citations, with 353 Scopus over 600 Google Scholar citations, and the highest average number of citations (67) that steadily increased from 10 citations in 2013 to 119 citations in 2021, was written in 2012 by Dr. Ruden's laboratory and titled, "Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift." The five most influential authors who published in the journal in the past 10 years based on Scopus citations of a particular paper are Dr. Ruden's laboratory, with 353 Scopus citations for the SnpSift paper mentioned above; Drs. Brock Christensen and Carmen J. Marsit, with 86 Scopus citations for their review, "Epigenomics in environmental health"; Dr. Michael Aschner and colleagues, with 61 Scopus citations for their paper "Genetic factors and manganese-induced neurotoxicity"; and Dr. Sandra C. dos Santos and colleagues, with 59 Scopus citations for their paper, "Yeast toxicogenomics: genome-wide responses to chemical stresses with impact in environmental health, pharmacology, and biotechnology." While the top 5 papers were published in the early years of the journal, we will also discuss a more recent article published in 2018 on a comparison of RNA-seq and microarray methods by Dr. Michael Liguori's laboratory, "Comparison of RNA-Seq and Microarray Gene Expression Platforms for the Toxicogenomic Evaluation of Liver From Short-Term Rat Toxicity Studies," that far exceeds the number of downloads and views of all the other articles published in the first 10 years of the journal and will likely be a top cited paper in the second decade highlights of this journal. Finally, we discuss where the Frontiers in Toxicogenomics specialty journal and the Frontiers in Toxicology journal will go to advance the field of toxicogenomics, and more generally, toxicology, in the future.

Drug resistance in NSCLC is associated with tumor micro-environment

BACKGROUND

Tumor cell resistance to chemotherapy is the most critical factor that influences the prognosis of cancer patients. It is generally believed that drug resistance is caused by genetic alterations in tumor cells; however, the relationship between drug resistance and the tumor microenvironment (TME) has not been adequately studied. Herein, we successfully identified drug resistance and sensitivity clusters using single-cell transcriptome sequencing data from GSE149383 and established a proportional hazards model to find genes that affected prognosis. The results showed that marker genes between resistant and sensitive clusters were significantly associated with the TME; additionally, the model showed good reliability. Furthermore, we used bulk RNA-seq data to analyze the expression of CD24 and CYP1B1, which revealed little difference in the levels of the two genes in normal and tumor tissues but a significant difference in their expression between drug-resistant and -sensitive cells. In conclusion, our study demonstrated a link between drug resistance and the TME, and we found that CD24 and CYP1B1 may be key regulators of drug resistance development in tumor cells via altering the TME.

Pre-adult exposure to three heavy metals leads to changes in the head transcriptome of adult flies

We examined the effect of developmental exposure to three heavy metals - cadmium, copper, and lead - on gene expression in adult head tissue in the model organism Drosophila melanogaster . All metals affected development time and/or gene expression level. While variation in the response to each metal was apparent, two differentially-expressed genes were upregulated in response to all three metal treatments, and 11 genes were downregulated in two of the three treatments. Our work reveals that developmental metal exposure has the potential to have long-lasting, metal-specific effects on gene expression in adults, even after the metal stress has been removed.

Metal-metal interaction and metal toxicity: a comparison between mammalian and D. melanogaster

1. Non-essential heavy metals such as mercury (Hg), arsenic (As), cadmium (Cd), and aluminium (Al) are useless to organisms and have shown extensive toxic effects. Previous studies show that two main molecular mechanisms of metal toxicity are oxidative stress and metal-metal interaction which can disrupt metal homeostasis.2. In this paper, we mainly illustrate metal toxicity and metal-metal interaction through examples in mammalians and D. melanogaster (fruit fly).3. We describe the interference of metal homeostasis by metal-metal interactions in three aspects including replacement, cellular transporter competition, and disruption of the regulation mechanism, and elaborate the mechanisms of metal toxicity to better deal with the challenges of heavy metal pollution and related health problems.

Characterizing the genetic basis of copper toxicity in Drosophila reveals a complex pattern of allelic, regulatory, and behavioral variation

A range of heavy metals are required for normal cell function and homeostasis. However, the anthropogenic release of metal compounds into soil and water sources presents a pervasive health threat. Copper is one of many heavy metals that negatively impacts diverse organisms at a global scale. Using a combination of quantitative trait locus (QTL) mapping and RNA sequencing in the Drosophila Synthetic Population Resource, we demonstrate that resistance to the toxic effects of ingested copper in D. melanogaster is genetically complex and influenced by allelic and expression variation at multiple loci. QTL mapping identified several QTL that account for a substantial fraction of heritability. Additionally, we find that copper resistance is impacted by variation in behavioral avoidance of copper and may be subject to life-stage specific regulation. Gene expression analysis further demonstrated that resistant and sensitive strains are characterized by unique expression patterns. Several of the candidate genes identified via QTL mapping and RNAseq have known copper-specific functions (e.g., Ccs, Sod3, CG11825), and others are involved in the regulation of other heavy metals (e.g., Catsup, whd). We validated several of these candidate genes with RNAi suggesting they contribute to variation in adult copper resistance. Our study illuminates the interconnected roles that allelic and expression variation, organism life stage, and behavior play in copper resistance, allowing a deeper understanding of the diverse mechanisms through which metal pollution can negatively impact organisms.

Monitoring the Heavy Metal Lead Inside Living Drosophila with a FRET-Based Biosensor

The harmful impact of the heavy metal lead on human health has been known for years. However, materials that contain lead remain in the environment. Measuring the blood lead level (BLL) is the only way to officially evaluate the degree of exposure to lead. The so-called "safe value" of the BLL seems to unreliably represent the secure threshold for children. In general, lead's underlying toxicological mechanism remains unclear and needs to be elucidated. Therefore, we developed a novel genetically encoded fluorescence resonance energy transfer (FRET)-based lead biosensor, Met-lead, and applied it to transgenic Drosophila to perform further investigations. We combined Met-lead with the UAS-GAL4 system to the sensor protein specifically expressed within certain regions of fly brains. Using a suitable imaging platform, including a fast epifluorescent or confocal laser-scanning/two-photon microscope with high resolution, we recorded the changes in lead content inside fly brains ex vivo and in vivo and at different life stages. The blood-brain barrier was found to play an important role in the protection of neurons in the brain against damage due to the heavy metal lead, either through food or microinjection into the abdomen. Met-lead has the potential to be a powerful tool for the sensing of lead within living organisms by employing either a fast epi-FRET microscope or high-resolution brain imaging.

Oxidative stress caused by lead (Pb) induces iron deficiency in Drosophila melanogaster

Toxic elements exposure disturbs the homeostasis of essential elements in organisms, but the mechanism remains elusive. In this study, we demonstrated that Drosophila melanogaster exposed to Lead (Pb, a pervasive environmental threat to human health) exhibited various health defects, including retarded development, decreased survival rate, impaired mobility and reduced egg production. These phenotypes could be significantly modulated by either intervention of dietary iron levels or altering expression of genes involved in iron metabolism. Further study revealed that Pb exposure leads to systemic iron deficiency. Strikingly, reactive oxygen species (ROS) clearance significantly increased iron uptake by restoring the expression of iron metabolism genes in the midgut and subsequently attenuated Pb toxicity. This study highlights the role of ROS in Pb induced iron dyshomeostasis and provides unique insights into understanding the mechanism of Pb toxicity and suggests ideal ways to attenuate Pb toxicity by iron supplementation therapy or ROS clearance.

Chronic lead (Pb) exposure results in diminished hemocyte count and increased susceptibility to bacterial infection in Drosophila melanogaster

Heavy metal Pb is a common toxic pollutant present in our environment adversely affecting health of the living organisms. Recent studies suggest positive correlation between heavy metal exposure and immune dysfunction and present work utilizes Drosophila to address this issue in relation to Pb exposure. In-vivo Pb toxicity was established by dietary intake where essential parameters like development and life span were found to be hampered and augmented upon metallothionein B (mtnB) downregulation hinting towards potential role of mtnB in Pb detoxification. Further response of Drosophila to B. subtilis bacterial infection was monitored by carrying out oral infections. Pb fed flies showed increased susceptibility to infection as compared to their controls. Since Drosophila hemocytes play dual role as immune cells, we checked for the total hemocyte count and found significant decrease in hemocyte numbers in Pb fed larvae. Both crystal cells and plasmatocytes, the two major hemocytes in third instar larval hemolymph were reduced. However we did not find any visible morphological changes in Giemsa stained hemocytes. Crystal cells are crucial for synthesis and release of phenoloxidase (PO), an enzyme required for melanin clot synthesis and deposition. PO activity assessed from total hemolymph protein isolates was found to be substantially decreased in Pb raised animals. Results were also confirmed by spot test and native gel activity assay of PO. Overall our results suggest immunotoxic effect of Pb through decrease in hemocyte count including crystal cell which in turn leads to decreased PO activity and increased susceptibility to B. subtilis.

Lead Modulates trans- and cis-Expression Quantitative Trait Loci (eQTLs) in Drosophila melanogaster Heads

Lead exposure has long been one of the most important topics in global public health because it is a potent developmental neurotoxin. Here, an eQTL analysis, which is the genome-wide association analysis of genetic variants with gene expression, was performed. In this analysis, the male heads of 79 Drosophila melanogaster inbred lines from Drosophila Synthetic Population Resource (DSPR) were treated with or without developmental exposure, from hatching to adults, to 250 μM lead acetate [Pb(C₂H₃O₂)₂]. The goal was to identify genomic intervals that influence the gene-expression response to lead. After detecting 1798 cis-eQTLs and performing an initial trans-eQTL analysis, we focused our analysis on lead-sensitive "trans-eQTL hotspots," defined as genomic regions that are associated with a cluster of genes in a lead-dependent manner. We noticed that the genes associated with one of the 14 detected trans-eQTL hotspots, Chr 2L: 6,250,000 could be roughly divided into two groups based on their differential expression profile patterns and different categories of function. This trans-eQTL hotspot validates one identified in a previous study using different recombinant inbred lines. The expression of all the associated genes in the trans-eQTL hotspot was visualized with hierarchical clustering analysis. Besides the overall expression profile patterns, the heatmap displayed the segregation of differential parental genetic contributions. This suggested that trans-regulatory regions with different genetic contributions from the parental lines have significantly different expression changes after lead exposure. We believe this study confirms our earlier study, and provides important insights to unravel the genetic variation in lead susceptibility in Drosophila model.

Experimental Protocol for Using Drosophila As an Invertebrate Model System for Toxicity Testing in the Laboratory

Emergent properties and external factors (population-level and ecosystem-level interactions, in particular) play important roles in mediating ecologically-important endpoints, though they are rarely considered in toxicological studies. D. melanogaster is emerging as a toxicology model for the behavioral, neurological, and genetic impacts of toxicants, to name a few. More importantly, species in the genus Drosophila can be utilized as a model system for an integrative framework approach to incorporate emergent properties and answer ecologically-relevant questions in toxicology research. The aim of this paper is to provide a protocol for exposing species in the genus Drosophila to pollutants to be used as a model system for a range of phenotypic outputs and ecologically-relevant questions. More specifically, this protocol can be used to 1) link multiple biological levels of organization and understand the impact of toxicants on both individual- and population-level fitness; 2) test the impact of toxicants at different stages of developmental exposure; 3) test multigenerational and evolutionary implications of pollutants; and 4) test multiple contaminants and stressors simultaneously.

Zinc Detoxification: A Functional Genomics and Transcriptomics Analysis in Drosophila melanogaster Cultured Cells

Cells require some metals, such as zinc and manganese, but excess levels of these metals can be toxic. As a result, cells have evolved complex mechanisms for maintaining metal homeostasis and surviving metal intoxication. Here, we present the results of a large-scale functional genomic screen in Drosophila cultured cells for modifiers of zinc chloride toxicity, together with transcriptomics data for wild-type or genetically zinc-sensitized cells challenged with mild zinc chloride supplementation. Altogether, we identified 47 genes for which knockdown conferred sensitivity or resistance to toxic zinc or manganese chloride treatment, and >1800 putative zinc-responsive genes. Analysis of the 'omics data points to the relevance of ion transporters, glutathione (GSH)-related factors, and conserved disease-associated genes in zinc detoxification. Specific genes identified in the zinc screen include orthologs of human disease-associated genes CTNS, PTPRN (also known as IA-2), and ATP13A2 (also known as PARK9). We show that knockdown of red dog mine (rdog; CG11897), a candidate zinc detoxification gene encoding an ABCC-type transporter family protein related to yeast cadmium factor (YCF1), confers sensitivity to zinc intoxication in cultured cells, and that rdog is transcriptionally upregulated in response to zinc stress. As there are many links between the biology of zinc and other metals and human health, the 'omics data sets presented here provide a resource that will allow researchers to explore metal biology in the context of diverse health-relevant processes.

Identification of Splicing Quantitative Trait Loci (sQTL) in Drosophila melanogaster with Developmental Lead (Pb2+) Exposure

Lead (Pb) poisoning has been a major public health issue globally and the recent Flint water crisis has drawn nation-wide attention to its effects. To better understand how lead plays a role as a neurotoxin, we utilized the Drosophila melanogaster model to study the genetic effects of lead exposure during development and identified lead-responsive genes. In our previous studies, we have successfully identified hundreds of lead-responsive expression QTLs (eQTLs) by using RNA-seq analysis on heads collected from the Drosophila Synthetic Population Resource. Cis-eQTLs, also known as allele-specific expression (ASE) polymorphisms, are generally single-nucleotide polymorphisms in the promoter regions of genes that affect expression of the gene, such as by inhibiting the binding of transcription factors. Trans-eQTLs are genes that regulate mRNA levels for many genes, and are generally thought to be SNPs in trans-acting transcription or translation factors. In this study, we focused our attention on alternative splicing events that are affected by lead exposure. Splicing QTLs (sQTLs), which can be caused by SNPs that alter splicing or alternative splicing (AS), such as by changing the sequence-specific binding affinity of splicing factors to the pre-mRNA. We applied two methods in search for sQTLs by using RNA-seq data from control and lead-exposed w¹¹¹⁸Drosophila heads. First, we used the fraction of reads in a gene that falls in each exon as the phenotype. Second, we directly compared the transcript counts among the various splicing isoforms as the phenotype. Among the 1,236 potential Pb-responsive sQTLs (p < 0.0001, FDR < 0.39), mostly cis-sQTLs, one of the most distinct genes is Dscam1 (Down Syndrome Cell Adhesion Molecule), which has over 30,000 potential alternative splicing isoforms. We have also identified a candidate Pb-responsive trans-sQTL hotspot that appears to regulate 129 genes that are enriched in the “cation channel” gene ontology category, suggesting a model in which alternative splicing of these channels might lead to an increase in the elimination of Pb²⁺ from the neurons encoding these channels. To our knowledge, this is the first paper that uses sQTL analyses to understand the neurotoxicology of an environmental toxin in any organism, and the first reported discovery of a candidate trans-sQTL hotspot.

Accumulation, elimination, sequestration, and genetic variation of lead (Pb2+) loads within and between generations of Drosophila melanogaster

We examined accumulation, sequestration, elimination, and genetic variation for lead (Pb) loads within and between generations of Drosophila melanogaster. Flies were reared in control or leaded medium at various doses and tested for their Pb loads at different stages of development (larvae, eclosion, newly-eclosed adults, and mature adults). Pb loads were tested using Inductively Coupled Plasma Mass Spectrometry (ICP-MS). We found that D. melanogaster readily accumulated Pb throughout their lifespan and the levels of accumulation increased with Pb exposure in the medium. Wandering third-instar larvae accumulated more Pb than mature adults; this phenomenon may be due to elimination of Pb in the pupal cases during eclosion and/or depuration in adults post-eclosion. The accumulated Pb in mature adults was not transferred to F1 mature adult offspring. Using a set of recombinant inbred strains, we identified a quantitative trait locus for adult Pb loads and found that genetic variation accounted for 34% of the variance in Pb load. We concluded that D. melanogaster is a useful model organism for evaluating changes in Pb loads during development, as well as between generations. Furthermore, we found that genetic factors can influence Pb loads; this provides an essential foundation for evaluating phenotypic variation induced by the toxic effects of Pb.

Armour GE, Thakkar E, Mackay TFC, Anholt RRH. A Drosophila model for toxicogenomics: Genetic variation in susceptibility to heavy metal exposure

The genetic factors that give rise to variation in susceptibility to environmental toxins remain largely unexplored. Studies on genetic variation in susceptibility to environmental toxins are challenging in human populations, due to the variety of clinical symptoms and difficulty in determining which symptoms causally result from toxic exposure; uncontrolled environments, often with exposure to multiple toxicants; and difficulty in relating phenotypic effect size to toxic dose, especially when symptoms become manifest with a substantial time lag. Drosophila melanogaster is a powerful model that enables genome-wide studies for the identification of allelic variants that contribute to variation in susceptibility to environmental toxins, since the genetic background, environmental rearing conditions and toxic exposure can be precisely controlled. Here, we used extreme QTL mapping in an outbred population derived from the D. melanogaster Genetic Reference Panel to identify alleles associated with resistance to lead and/or cadmium, two ubiquitous environmental toxins that present serious health risks. We identified single nucleotide polymorphisms (SNPs) associated with variation in resistance to both heavy metals as well as SNPs associated with resistance specific to each of them. The effects of these SNPs were largely sex-specific. We applied mutational and RNAi analyses to 33 candidate genes and functionally validated 28 of them. We constructed networks of candidate genes as blueprints for orthologous networks of human genes. The latter not only provided functional contexts for known human targets of heavy metal toxicity, but also implicated novel candidate susceptibility genes. These studies validate Drosophila as a translational toxicogenomics gene discovery system.

Although physiological effects of environmental toxins are well documented, we know little about the genetic factors that determine individual variation in susceptibility to toxins. Such information is difficult to obtain in human populations due to heterogeneity in genetic background and environmental exposure, and the diversity of symptoms and time lag with which they appear after toxic exposure. Here, we show that the fruit fly, Drosophila, can serve as a powerful genetic model system to elucidate the genetic underpinnings that contribute to individual variation in resistance to toxicity, using lead and cadmium exposure as an experimental paradigm. We identified genes that harbor genetic variants that contribute to individual variation in resistance to heavy metal exposure. Furthermore, we constructed genetic networks on which we could superimpose human counterparts of Drosophila genes. We were able to place human genes previously implicated in heavy metal toxicity in biological context and identify novel targets for heavy metal toxicity. Thus, we demonstrate that based on evolutionary conservation of fundamental biological processes, we can use Drosophila as a powerful translational model for toxicogenomics studies.

Drosophila melanogaster Models of Metal-Related Human Diseases and Metal Toxicity

Iron, copper and zinc are transition metals essential for life because they are required in a multitude of biological processes. Organisms have evolved to acquire metals from nutrition and to maintain adequate levels of each metal to avoid damaging effects associated with its deficiency, excess or misplacement. Interestingly, the main components of metal homeostatic pathways are conserved, with many orthologues of the human metal-related genes having been identified and characterized in Drosophila melanogaster. Drosophila has gained appreciation as a useful model for studying human diseases, including those caused by mutations in pathways controlling cellular metal homeostasis. Flies have many advantages in the laboratory, such as a short life cycle, easy handling and inexpensive maintenance. Furthermore, they can be raised in a large number. In addition, flies are greatly appreciated because they offer a considerable number of genetic tools to address some of the unresolved questions concerning disease pathology, which in turn could contribute to our understanding of the metal metabolism and homeostasis. This review recapitulates the metabolism of the principal transition metals, namely iron, zinc and copper, in Drosophila and the utility of this organism as an experimental model to explore the role of metal dyshomeostasis in different human diseases. Finally, a summary of the contribution of Drosophila as a model for testing metal toxicity is provided.

Asymmetrical positive assortative mating induced by developmental lead (Pb2+) exposure in a model system, Drosophila melanogaster

Anthropogenic pollutants have the potential to disrupt reproductive strategies. Little is known about how lead (Pb²⁺) exposure disrupts individual-level responses in reproductive behaviors, which are important for fitness. Drosophila melanogaster was used as a model system to determine the effects of: 1) developmental lead exposure on pre-mating reproductive behaviors (i.e., mate preference), and 2) lead exposure and mating preferences on fitness in the F₀ parental generation and F₁ un-exposed offspring. Wild-type strains of D. melanogaster were reared from egg stage to adulthood in control or leaded medium (250 μM PbAc) and tested for differences in: mate preference, male song performance, sex pheromone expression, fecundity, mortality, and body weight. F₀ leaded females preferentially mated with leaded males (i.e., asymmetrical positive assortative mating) in 2-choice tests. This positive assortative mating was mediated by the females (and not the males) and was dependent upon context and developmental exposure to Pb. Neither the courtship song nor the sex pheromone profile expressed by control and leaded males mediated the positive assortative mating in leaded females. Leaded females did not incur a fitness cost in terms of reduced fecundity, increased mortality, or decreased body weight by mating with leaded males. These results suggest that sublethal exposure to lead during development can alter mate preferences in adults, but not fitness measures once lead exposure has been removed. We suggest that changes in mate preference may induce fitness costs, as well as long-term population and multi-generational implications, if pollution is persistent in the environment.

The Genetic Basis for Variation in Sensitivity to Lead Toxicity in Drosophila melanogaster

BACKGROUND

Lead toxicity presents a worldwide health problem, especially due to its adverse effects on cognitive development in children. However, identifying genes that give rise to individual variation in susceptibility to lead toxicity is challenging in human populations.

OBJECTIVES

Our goal was to use Drosophila melanogaster to identify evolutionarily conserved candidate genes associated with individual variation in susceptibility to lead exposure.

METHODS

To identify candidate genes associated with variation in susceptibility to lead toxicity, we measured effects of lead exposure on development time, viability and adult activity in the Drosophila melanogaster Genetic Reference Panel (DGRP) and performed genome-wide association analyses to identify candidate genes. We used mutants to assess functional causality of candidate genes and constructed a genetic network associated with variation in sensitivity to lead exposure, on which we could superimpose human orthologs.

RESULTS

We found substantial heritabilities for all three traits and identified candidate genes associated with variation in susceptibility to lead exposure for each phenotype. The genetic architectures that determine variation in sensitivity to lead exposure are highly polygenic. Gene ontology and network analyses showed enrichment of genes associated with early development and function of the nervous system.

CONCLUSIONS

Drosophila melanogaster presents an advantageous model to study the genetic underpinnings of variation in susceptibility to lead toxicity. Evolutionary conservation of cellular pathways that respond to toxic exposure allows predictions regarding orthologous genes and pathways across phyla. Thus, studies in the D. melanogaster model system can identify candidate susceptibility genes to guide subsequent studies in human populations.

CITATION

Zhou S, Morozova TV, Hussain YN, Luoma SE, McCoy L, Yamamoto A, Mackay TF, Anholt RR. 2016. The genetic basis for variation in sensitivity to lead toxicity in Drosophila melanogaster. Environ Health Perspect 124:1062-1070; http://dx.doi.org/10.1289/ehp.1510513.

Gene Expression Evolves under a House-of-Cards Model of Stabilizing Selection

Divergence in gene regulation is hypothesized to underlie much of phenotypic evolution, but the role of natural selection in shaping the molecular phenotype of gene expression continues to be debated. To resolve the mode of gene expression, evolution requires accessible theoretical predictions for the effect of selection over long timescales. Evolutionary quantitative genetic models of phenotypic evolution can provide such predictions, yet those predictions depend on the underlying hypotheses about the distributions of mutational and selective effects that are notoriously difficult to disentangle. Here, we draw on diverse genomic data sets including expression profiles of natural genetic variation and mutation accumulation lines, empirical estimates of genomic mutation rates, and inferences of genetic architecture to differentiate contrasting hypotheses for the roles of stabilizing selection and mutation in shaping natural expression variation. Our analysis suggests that gene expression evolves in a domain of phenotype space well fit by the House-of-Cards (HC) model. Although the strength of selection inferred is sensitive to the number of loci controlling gene expression, the model is not. The consistency of these results across evolutionary time from budding yeast through fruit fly implies that this model is general and that mutational effects on gene expression are relatively large. Empirical estimates of the genetic architecture of gene expression traits imply that selection provides modest constraints on gene expression levels for most genes, but that the potential for regulatory evolution is high. Our prediction using data from laboratory environments should encourage the collection of additional data sets allowing for more nuanced parameterizations of HC models for gene expression.

Epigenetics as an answer to Darwin's "special difficulty," Part 2: natural selection of metastable epialleles in honeybee castes

In a recent perspective in this journal, Herb (2014) discussed how epigenetics is a possible mechanism to circumvent Charles Darwin's "special difficulty" in using natural selection to explain the existence of the sterile-fertile dimorphism in eusocial insects. Darwin's classic book "On the Origin of Species by Means of Natural Selection" explains how natural selection of the fittest individuals in a population can allow a species to adapt to a novel or changing environment. However, in bees and other eusocial insects, such as ants and termites, there exist two or more castes of genetically similar females, from fertile queens to multiple sub-castes of sterile workers, with vastly different phenotypes, lifespans, and behaviors. This necessitates the selection of groups (or kin) rather than individuals in the evolution of honeybee hives, but group and kin selection theories of evolution are controversial and mechanistically uncertain. Also, group selection would seem to be prohibitively inefficient because the effective population size of a colony is reduced from thousands to a single breeding queen. In this follow-up perspective, we elaborate on possible mechanisms for how a combination of both epigenetics, specifically, the selection of metastable epialleles, and genetics, the selection of mutations generated by the selected metastable epialleles, allows for a combined means for selection amongst the fertile members of a species to increase colony fitness. This "intra-caste evolution" hypothesis is a variation of the epigenetic directed genetic error hypothesis, which proposes that selected metastable epialleles increase genetic variability by directing mutations specifically to the epialleles. Natural selection of random metastable epialleles followed by a second round of natural selection of random mutations generated by the metastable epialleles would allow a way around the small effective population size of eusocial insects.

Xenotransplantation models to study the effects of toxicants on human fetal tissues

Many diseases that manifest throughout the lifetime are influenced by factors affecting fetal development. Fetal exposure to xenobiotics, in particular, may influence the development of adult diseases. Established animal models provide systems for characterizing both developmental biology and developmental toxicology. However, animal model systems do not allow researchers to assess the mechanistic effects of toxicants on developing human tissue. Human fetal tissue xenotransplantation models have recently been implemented to provide human-relevant mechanistic data on the many tissue-level functions that may be affected by fetal exposure to toxicants. This review describes the development of human fetal tissue xenotransplant models for testis, prostate, lung, liver, and adipose tissue, aimed at studying the effects of xenobiotics on tissue development, including implications for testicular dysgenesis, prostate disease, lung disease, and metabolic syndrome. The mechanistic data obtained from these models can complement data from epidemiology, traditional animal models, and in vitro studies to quantify the risks of toxicant exposures during human development.

Whole adult organism transcriptional profiling of acute metal exposures in male zebrafish

Background

A convergence of technological breakthroughs in the past decade has facilitated the development of rapid screening tools for biomarkers of toxicant exposure and effect. Platforms using the whole adult organism to evaluate the genome-wide response to toxicants are especially attractive. Recent work demonstrates the feasibility of this approach in vertebrates using the experimentally robust zebrafish model. In the present study, we evaluated gene expression changes in whole adult male zebrafish following an acute 24 hr high dose exposure to three metals with known human health risks. Male adult zebrafish were exposed to nickel chloride, cobalt chloride or sodium dichromate concentrations corresponding to their respective 96 hr LC₂₀, LC₄₀ and LC₆₀. Histopathology was performed on a subset of metal-exposed zebrafish to phenotypically anchor transcriptional changes associated with each metal.

Results

Comparative analysis identified subsets of differentially expressed transcripts both overlapping and unique to each metal. Application of gene ontology (GO) and transcription factor (TF) enrichment algorithms revealed a number of key biological processes perturbed by metal poisonings and the master transcriptional regulators mediating gene expression changes. Metal poisoning differentially activated biological processes associated with ribosome biogenesis, proteosomal degradation, and p53 signaling cascades, while repressing oxygen-generating pathways associated with amino acid and lipid metabolism. Despite appreciable effects on gene regulation, nickel poisoning did not induce any morphological alterations in male zebrafish organs and tissues. Histopathological effects of cobalt remained confined to the olfactory system, while chromium targeted the gills, pharynx, and intestinal mucosa. A number of enriched transcription factors mediated the observed gene response to metal poisoning, including known targets such as p53, HIF1α, and the myc oncogene, and novel regulatory factors such as XBP1, GATA6 and HNF3β.

Conclusions

This work uses an experimentally innovative approach to capture global responses to metal poisoning and provides mechanistic insights into metal toxicity.

Systems-level approaches reveal conservation of trans-regulated genes in the rat and genetic determinants of blood pressure in humans

AIMS

Human genome-wide association studies (GWAS) of hypertension identified only few susceptibility loci with large effect that were replicated across populations. The vast majority of genes detected by GWAS has small effect and the regulatory mechanisms through which these genetic variants cause disease remain mostly unclear. Here, we used comparative genomics between human and an established rat model of hypertension to explore the transcriptional mechanisms mediating the effect of genes identified in 15 hypertension GWAS.

METHODS AND RESULTS

Time series analysis of radiotelemetric blood pressure (BP) was performed to assess 11 parameters of BP variation in recombinant inbred strains derived from the spontaneously hypertensive rat. BP data were integrated with ∼27 000 expression quantative trait loci (eQTLs) mapped across seven tissues, detecting >8000 significant associations between eQTL genes and BP variation in the rat. We then compiled a large catalogue of human genes from GWAS of hypertension and identified a subset of 2292 rat-human orthologous genes. Expression levels for 795 (34%) of these genes correlated with BP variation across rat tissues: 51 genes were cis-regulated, whereas 459 were trans-regulated and enriched for 'calcium signalling pathway' (P = 9.6 × 10(-6)) and 'ion channel' genes (P = 3.5 × 10(-7)), which are important determinants of hypertension. We identified 158 clusters of trans-eQTLs, annotated the underlying 'master regulator' genes and found significant over-representation in the human hypertension gene set (enrichment P = 5 × 10(-4)).

CONCLUSION

We showed extensive conservation of trans-regulated genes and their master regulators between rat and human hypertension. These findings reveal that small-effect genes associated with hypertension by human GWAS are likely to exert their action through coordinate regulation of pathogenic pathways.

Developmental lead effects on behavior and brain gene expression in male and female BALB/cAnNTac mice

Lead (Pb) was one of the first poisons identified, and the developing nervous system is particularly vulnerable to its toxic effects. Relatively low, subclinical doses, of Pb that produce no overt signs of encephalopathy can affect cognitive, emotional, and motor functions. In the present study, the effects of developmental Pb-exposure on behavioral performance and gene expression in BALB/cAnNTac mice were evaluated. Pups were exposed to Pb from gestational-day (gd) 8 to postnatal-day (pnd) 21 and later evaluated in exploratory behavior, rotarod, Morris water maze, and resident-intruder assays as adults. Pb-exposure caused significant alterations in exploratory behavior and water maze performance during the probe trial, but rotarod performance was not affected. Pb-exposed males displayed violent behavior towards their cage mates, but not to a stranger in the resident-intruder assay. Gene expression analysis at pnd21 by microarray and qRT-PCR was performed to provide a molecular link to the behavior changes that were observed. Pb strongly up-regulated gene expression within the signaling pathways of mitogen activated protein kinases (MAPKs), extra-cellular matrix (ECM) receptor, focal adhesion, and vascular endothelial growth-factor (VEGF), but Pb down-regulated gene expression within the pathways for glycan structures-biosynthesis 1, purine metabolism, and N-glycan biosynthesis. Pb increased transcription of genes for major histocompatibility (MHC) proteins, the chemokine Ccl28, chemokine receptors, IL-7, IL7R, and proteases. The qRT-PCR analysis indicated an increase of gene expression in the whole brain for caspase 1 and NOS2. Analysis of IL-1β, caspase 1, NOS2, Trail, IL-18 and IL-33 gene expression of brain regions indicated that Pb perturbed the inter-regional expression pattern of pro-inflammatory genes. Brain region protein concentrations for IL-10, an anti-inflammatory cytokine, showed a significant decrease only within the cortex region. Results indicate that Pb differentially affects the behavior of male and female mice in that females did less exploration and the males were selectively more aggressive. Gene expression data pointed to evidence of neuroinflammation in the brain of both female and male mice. Pb had more of an effect in the males on expression of vomeronasal receptor genes associated with odor detection and social behavior.

Drosophila melanogaster as a model for lead neurotoxicology and toxicogenomics research

Drosophila melanogaster is an excellent model animal for studying the neurotoxicology of lead. It has been known since ancient Roman times that long-term exposure to low levels of lead results in behavioral abnormalities, such as what is now known as attention deficit hyperactivity disorder (ADHD). Because lead alters mechanisms that underlie developmental neuronal plasticity, chronic exposure of children, even at blood lead levels below the current CDC community action level (10 μg/dl), can result in reduced cognitive ability, increased likelihood of delinquency, behaviors associated with ADHD, changes in activity level, altered sensory function, delayed onset of sexual maturity in girls, and changes in immune function. In order to better understand how lead affects neuronal plasticity, we will describe recent findings from a Drosophila behavioral genetics laboratory, a Drosophila neurophysiology laboratory, and a Drosophila quantitative genetics laboratory who have joined forces to study the effects of lead on the Drosophila nervous system. Studying the effects of lead on Drosophila nervous system development will give us a better understanding of the mechanisms of Pb neurotoxicity in the developing human nervous system.

Population differences in transcript-regulator expression quantitative trait loci

Gene expression quantitative trait loci (eQTL) are useful for identifying single nucleotide polymorphisms (SNPs) associated with diseases. At times, a genetic variant may be associated with a master regulator involved in the manifestation of a disease. The downstream target genes of the master regulator are typically co-expressed and share biological function. Therefore, it is practical to screen for eQTLs by identifying SNPs associated with the targets of a transcript-regulator (TR). We used a multivariate regression with the gene expression of known targets of TRs and SNPs to identify TReQTLs in European (CEU) and African (YRI) HapMap populations. A nominal p-value of <1×10(-6) revealed 234 SNPs in CEU and 154 in YRI as TReQTLs. These represent 36 independent (tag) SNPs in CEU and 39 in YRI affecting the downstream targets of 25 and 36 TRs respectively. At a false discovery rate (FDR) = 45%, one cis-acting tag SNP (within 1 kb of a gene) in each population was identified as a TReQTL. In CEU, the SNP (rs16858621) in Pcnxl2 was found to be associated with the genes regulated by CREM whereas in YRI, the SNP (rs16909324) was linked to the targets of miRNA hsa-miR-125a. To infer the pathways that regulate expression, we ranked TReQTLs by connectivity within the structure of biological process subtrees. One TReQTL SNP (rs3790904) in CEU maps to Lphn2 and is associated (nominal p-value = 8.1×10(-7)) with the targets of the X-linked breast cancer suppressor Foxp3. The structure of the biological process subtree and a gene interaction network of the TReQTL revealed that tumor necrosis factor, NF-kappaB and variants in G-protein coupled receptors signaling may play a central role as communicators in Foxp3 functional regulation. The potential pleiotropic effect of the Foxp3 TReQTLs was gleaned from integrating mRNA-Seq data and SNP-set enrichment into the analysis.

Mapping gene-environment interactions at regulatory polymorphisms: insights into mechanisms of phenotypic variation

Genetic effects on gene regulation make a substantial contribution to phenotypic diversity, yet their mechanisms remain elusive. Here, we discuss the potential insights to be gained from mapping gene-environment interactions at regulatory polymorphisms (i.e., genetic variation that affects gene expression under specific environmental conditions). We highlight a novel statistical method to identify specific patterns of gene-environment interaction at these regulatory polymorphisms. Reviewing its application to a study that mapped gene expression in the presence and absence of glucocorticoids, we discuss the mechanistic insights that this approach provides.

Peripheral blood signatures of lead exposure

BACKGROUND

Current evidence indicates that even low-level lead (Pb) exposure can have detrimental effects, especially in children. We tested the hypothesis that Pb exposure alters gene expression patterns in peripheral blood cells and that these changes reflect dose-specific alterations in the activity of particular pathways.

METHODOLOGY/PRINCIPAL FINDING

Using Affymetrix Mouse Genome 430 2.0 arrays, we examined gene expression changes in the peripheral blood of female Balb/c mice following exposure to per os lead acetate trihydrate or plain drinking water for two weeks and after a two-week recovery period. Data sets were RMA-normalized and dose-specific signatures were generated using established methods of supervised classification and binary regression. Pathway activity was analyzed using the ScoreSignatures module from GenePattern.

CONCLUSIONS/SIGNIFICANCE

The low-level Pb signature was 93% sensitive and 100% specific in classifying samples a leave-one-out crossvalidation. The high-level Pb signature demonstrated 100% sensitivity and specificity in the leave-one-out crossvalidation. These two signatures exhibited dose-specificity in their ability to predict Pb exposure and had little overlap in terms of constituent genes. The signatures also seemed to reflect current levels of Pb exposure rather than past exposure. Finally, the two doses showed differential activation of cellular pathways. Low-level Pb exposure increased activity of the interferon-gamma pathway, whereas high-level Pb exposure increased activity of the E2F1 pathway.

Identification of methylmercury tolerance gene candidates in Drosophila

Methylmercury (MeHg) is a ubiquitous environmental contaminant that preferentially targets the developing nervous system. Variable outcomes of prenatal MeHg exposure within a population point to a genetic component that regulates MeHg toxicity. We therefore sought to identify fundamental MeHg tolerance genes using the Drosophila model for genetic and molecular dissection of a MeHg tolerance trait. We observe autosomal dominance in a MeHg tolerance trait (development on MeHg food) in both wild-derived and laboratory-selected MeHg-tolerant strains of flies. We performed whole-genome transcript profiling of larval brains of tolerant (laboratory selected) and nontolerant (control) strains in the presence and absence of MeHg stress. Pairwise transcriptome comparisons of four conditions (+/-selection and +/-MeHg) identified a "down-down-up" expression signature, whereby MeHg alone and selection alone resulted in a greater number of downregulated transcripts, and the combination of selection + MeHg resulted in a greater number of upregulated transcripts. Functional annotation cluster analyses showed enrichment for monooxygenases/oxidoreductases, which include cytochrome P450 (CYP) family members. Among the 10 CYPs upregulated with selection + MeHg in tolerant strains, CYP6g1, previously identified as the dichlorodiphenyl trichloroethane resistance allele in flies, was the most highly expressed and responsive to MeHg. Among all the genes, Turandot A (TotA), an immune pathway-regulated humoral response gene, showed the greatest upregulation with selection + MeHg. Neural-specific transgenic overexpression of TotA enhanced MeHg tolerance during pupal development. Identification of TotA and CYP genes as MeHg tolerance genes is an inroad to investigating the conserved function of immune signaling and phase I metabolism pathways in MeHg toxicity and tolerance in higher organisms.

Pb2+: an endocrine disruptor in Drosophila?

Environmental exposure to Pb(2+) affects hormone-mediated responses in vertebrates. To help establish the fruit fly, Drosophila melanogaster, as a model system for studying such disruption, we describe effects of Pb(2+) on hormonally regulated traits. These include duration of development, longevity, females' willingness to mate, fecundity and adult locomotor activity. Developmental Pb(2+) exposure has been shown to affect gene expression in a specific region of the Drosophila genome (approximately 122 genes) involved in lead-induced changes in adult locomotion and to affect regulation of intracellular calcium levels associated with neuronal activity at identified synapses in the larval neuromuscular junction. We suggest ways in which Drosophila could become a new model system for the study of endocrine disruptors at genetic, neural and behavioral levels of analysis, particularly by use of genomic methods. This will facilitate efforts to distinguish between behavioral effects of Pb(2+) caused by direct action on neural mechanisms versus effects of Pb(+2) on behavior mediated through endocrine disruption.