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

Systems genomics of salinity stress response in rice

Populations can adapt to stressful environments through changes in gene expression. However, the role of gene regulation in mediating stress response and adaptation remains largely unexplored. Here, we use an integrative field dataset obtained from 780 plants of Oryza sativa ssp. indica (rice) grown in a field experiment under normal or moderate salt stress conditions to examine selection and evolution of gene expression variation under salinity stress conditions. We find that salinity stress induces increased selective pressure on gene expression. Further, we show that trans-eQTLs rather than cis-eQTLs are primarily associated with rice's gene expression under salinity stress, potentially via a few master-regulators. Importantly, and contrary to the expectations, we find that cis-trans reinforcement is more common than cis-trans compensation which may be reflective of rice diversification subsequent to domestication. We further identify genetic fixation as the likely mechanism underlying this compensation/reinforcement. Additionally, we show that cis- and trans-eQTLs are under different selection regimes, giving us insights into the evolutionary dynamics of gene expression variation. By examining genomic, transcriptomic, and phenotypic variation across a rice population, we gain insights into the molecular and genetic landscape underlying adaptive salinity stress responses, which is relevant for other crops and other stresses.

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.

Molecular mechanisms underpinning quantitative resistance to Phytophthora sojae in Glycine max using a systems genomics approach

Expression of quantitative disease resistance in many host-pathogen systems is controlled by genes at multiple loci, each contributing a small effect to the overall response. We used a systems genomics approach to study the molecular underpinnings of quantitative disease resistance in the soybean-Phytophthora sojae pathosystem, incorporating expression quantitative trait loci (eQTL) mapping and gene co-expression network analysis to identify the genes putatively regulating transcriptional changes in response to inoculation. These findings were compared to previously mapped phenotypic (phQTL) to identify the molecular mechanisms contributing to the expression of this resistance. A subset of 93 recombinant inbred lines (RILs) from a Conrad × Sloan population were inoculated with P. sojae isolate 1.S.1.1 using the tray-test method; RNA was extracted, sequenced, and the normalized read counts were genetically mapped from tissue collected at the inoculation site 24 h after inoculation from both mock and inoculated samples. In total, more than 100,000 eQTLs were mapped. There was a switch from predominantly cis-eQTLs in the mock treatment to an almost entirely nonoverlapping set of predominantly trans-eQTLs in the inoculated treatment, where greater than 100-fold more eQTLs were mapped relative to mock, indicating vast transcriptional reprogramming due to P. sojae infection occurred. The eQTLs were organized into 36 hotspots, with the four largest hotspots from the inoculated treatment corresponding to more than 70% of the eQTLs, each enriched for genes within plant-pathogen interaction pathways. Genetic regulation of trans-eQTLs in response to the pathogen was predicted to occur through transcription factors and signaling molecules involved in plant-pathogen interactions, plant hormone signal transduction, and MAPK pathways. Network analysis identified three co-expression modules that were correlated with susceptibility to P. sojae and associated with three eQTL hotspots. Among the eQTLs co-localized with phQTLs, two cis-eQTLs with putative functions in the regulation of root architecture or jasmonic acid, as well as the putative master regulators of an eQTL hotspot nearby a phQTL, represent candidates potentially underpinning the molecular control of these phQTLs for resistance.

Viruses in Laboratory Drosophila and Their Impact on Host Gene Expression

Drosophila melanogaster has one of the best characterized antiviral immune responses among invertebrates. However, relatively few easily transmitted natural virus isolates are available, and so many Drosophila experiments have been performed using artificial infection routes and artificial host-virus combinations. These may not reflect natural infections, especially for subtle phenotypes such as gene expression. Here, to explore the laboratory virus community and to better understand how natural virus infections induce changes in gene expression, we have analysed seven publicly available D. melanogaster transcriptomic sequencing datasets that were originally sequenced for projects unrelated to virus infection. We have found ten known viruses-including five that have not been experimentally isolated-but no previously unknown viruses. Our analysis of host gene expression revealed that numerous genes were differentially expressed in flies that were naturally infected with a virus. For example, flies infected with nora virus showed patterns of gene expression consistent with intestinal vacuolization and possible host repair via the upd3 JAK/STAT pathway. We also found marked sex differences in virus-induced differential gene expression. Our results show that natural virus infection in laboratory Drosophila does indeed induce detectable changes in gene expression, suggesting that this may form an important background condition for experimental studies in the laboratory.

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.

cis-regulatory variation modulates susceptibility to enteric infection in the Drosophila genetic reference panel

BACKGROUND

Resistance to enteric pathogens is a complex trait at the crossroads of multiple biological processes. We have previously shown in the Drosophila Genetic Reference Panel (DGRP) that resistance to infection is highly heritable, but our understanding of how the effects of genetic variants affect different molecular mechanisms to determine gut immunocompetence is still limited.

RESULTS

To address this, we perform a systems genetics analysis of the gut transcriptomes from 38 DGRP lines that were orally infected with Pseudomonas entomophila. We identify a large number of condition-specific, expression quantitative trait loci (local-eQTLs) with infection-specific ones located in regions enriched for FOX transcription factor motifs. By assessing the allelic imbalance in the transcriptomes of 19 F1 hybrid lines from a large round robin design, we independently attribute a robust cis-regulatory effect to only 10% of these detected local-eQTLs. However, additional analyses indicate that many local-eQTLs may act in trans instead. Comparison of the transcriptomes of DGRP lines that were either susceptible or resistant to Pseudomonas entomophila infection reveals nutcracker as the only differentially expressed gene. Interestingly, we find that nutcracker is linked to infection-specific eQTLs that correlate with its expression level and to enteric infection susceptibility. Further regulatory analysis reveals one particular eQTL that significantly decreases the binding affinity for the repressor Broad, driving differential allele-specific nutcracker expression.

CONCLUSIONS

Our collective findings point to a large number of infection-specific cis- and trans-acting eQTLs in the DGRP, including one common non-coding variant that lowers enteric infection susceptibility.

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.

A 24 h Age Difference Causes Twice as Much Gene Expression Divergence as 100 Generations of Adaptation to a Novel Environment

Gene expression profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. Here, we show that a 24 h age difference of Drosophilasimulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. Our results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework.