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

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.

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.

Novel role of prostate cancer risk variant rs7247241 on PPP1R14A isoform transition through allelic TF binding and CpG methylation

Although previous studies identified numerous single nucleotide polymorphisms (SNPs) and their target genes predisposed to prostate cancer (PrCa) risks, SNP-related splicing associations are rarely reported. In this study, we applied distance-based sQTL analysis (sQTLseekeR) using RNA-seq and SNP genotype data from benign prostate tissue (n = 467), and identified significant associations in 3344 SNP-transcript pairs (P ≤ 0.05) at PrCa risk loci. We characterized a common SNP (rs7247241) and its target gene (PPP1R14A) located in chr19q13, a sQTL with risk allele T associated with upregulation of long isoform (P = 9.99E-7). We confirmed the associations in both TCGA (P = 2.42E-24) and GTEX prostate cohorts (P = 9.08E-78). To functionally characterize this SNP, we performed chromatin Immunoprecipitation qPCR and confirmed stronger CTCF and PLAGL2 binding in rs7247241 C than T allele. We found that CTCF binding enrichment was negatively associated with methylation level at the SNP site in human cell lines (r = -0.58). Bisulfite sequencing showed consistent association of rs7247241-T allele with nearby sequence CpG hypermethylation in prostate cell lines and tissues. Importantly, the methylation level at CpG sites nearest to the CTCF binding and first exon splice-in (ψ) of PPP1R14A were significantly associated with aggressive phenotype in the TCGA PrCa cohort. Meanwhile, the long isoform of the gene also promoted cell proliferation. Taken together, with the most updated gene annotations, we reported a set of sQTL associated with multiple traits related to human prostate diseases, and revealed a unique role of PrCa risk SNP rs7247241 on PPP1R14A isoform transition.

ExonSkipDB: functional annotation of exon skipping event in human

Exon skipping (ES) is reported to be the most common alternative splicing event due to loss of functional domains/sites or shifting of the open reading frame (ORF), leading to a variety of human diseases and considered therapeutic targets. To date, systematic and intensive annotations of ES events based on the skipped exon units in cancer and normal tissues are not available. Here, we built ExonSkipDB, the ES annotation database available at https://ccsm.uth.edu/ExonSkipDB/, aiming to provide a resource and reference for functional annotation of ES events in multiple cancer and tissues to identify therapeutically targetable genes in individual exon units. We collected 14 272 genes that have 90 616 and 89 845 ES events across 33 cancer types and 31 normal tissues from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx). For the ES events, we performed multiple functional annotations. These include ORF assignment of exon skipped transcript, studies of lost protein functional features due to ES events, and studies of exon skipping events associated with mutations and methylations based on multi-omics evidence. ExonSkipDB will be a unique resource for cancer and drug research communities to identify therapeutically targetable exon skipping events.

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.

Genome-Wide Identification of Splicing Quantitative Trait Loci (sQTLs) in Diverse Ecotypes of Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs contributes to transcriptome diversity and enables plants to generate different protein isoforms from a single gene and/or fine-tune gene expression during different development stages and environmental changes. Although AS is pervasive, the genetic basis for differential isoform usage in plants is still emerging. In this study, we performed genome-wide analysis in 666 geographically distributed diverse ecotypes of Arabidopsis thaliana to identify genomic regions [splicing quantitative trait loci (sQTLs)] that may regulate differential AS. These ecotypes belong to different microclimatic conditions and are part of the relict and non-relict populations. Although sQTLs were spread across the genome, we observed enrichment for trans-sQTL (trans-sQTLs hotspots) on chromosome one. Furthermore, we identified several sQTL (911) that co-localized with trait-linked single nucleotide polymorphisms (SNP) identified in the Arabidopsis genome-wide association studies (AraGWAS). Many sQTLs were enriched among circadian clock, flowering, and stress-responsive genes, suggesting a role for differential isoform usage in regulating these important processes in diverse ecotypes of Arabidopsis. In conclusion, the current study provides a deep insight into SNPs affecting isoform ratios/genes and facilitates a better mechanistic understanding of trait-associated SNPs in GWAS studies. To the best of our knowledge, this is the first report of sQTL analysis in a large set of Arabidopsis ecotypes and can be used as a reference to perform sQTL analysis in the Brassicaceae family. Since whole genome and transcriptome datasets are available for these diverse ecotypes, it could serve as a powerful resource for the biological interpretation of trait-associated loci, splice isoform ratios, and their phenotypic consequences to help produce more resilient and high yield crop varieties.

Environmental influences on RNA processing: Biochemical, molecular and genetic regulators of cellular response

RNA processing has emerged as a key mechanistic step in the regulation of the cellular response to environmental perturbation. Recent work has uncovered extensive remodeling of transcriptome composition upon environmental perturbation and linked the impacts of this molecular plasticity to health and disease outcomes. These isoform changes and their underlying mechanisms are varied-involving alternative sites of transcription initiation, alternative splicing, and alternative cleavage at the 3' end of the mRNA. The mechanisms and consequences of differential RNA processing have been characterized across a range of common environmental insults, including chemical stimuli, immune stimuli, heat stress, and cancer pathogenesis. In each case, there are perturbation-specific contributions of local (cis) regulatory elements or global (trans) factors and downstream consequences. Overall, it is clear that choices in isoform usage involve a balance between the usage of specific genetic elements (i.e., splice sites, polyadenylation sites) and the timing at which certain decisions are made (i.e., transcription elongation rate). Fine-tuned cellular responses to environmental perturbation are often dependent on the genetic makeup of the cell. Genetic analyses of interindividual variation in splicing have identified genetic effects on splicing that contribute to variation in complex traits. Finally, the increase in the number of tissue types and environmental conditions analyzed for RNA processing is paralleled by the need to develop appropriate analytical tools. The combination of large datasets, novel methods and conditions explored promises to provide a much greater understanding of the role of RNA processing response in human phenotypic variation. This article is categorized under: RNA Processing > RNA Editing and Modification RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

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.