Exposure to the pesticide tefluthrin causes developmental neurotoxicity in zebrafish.

BACKGROUND: The insecticide tefluthrin is widely used in agriculture, resulting in widespread pollution. Tefluthrin is a type I pyrethroid characterized by its high persistence in the environment. Understanding the mechanisms of toxicity of tefluthrin will improve its risk assessment. OBJECTIVES: We aimed to decipher the molecular modes of action of tefluthrin. METHODS: Phenotypic developmental toxicity was assessed by exposing zebrafish embryos and larvae to increasing concentrations of tefluthrin. Tg(mnx:mGFP) line was used to assess neurotoxicity. Multi-omics approaches including transcriptomics and lipidomics were applied to analyze RNA and lipid contents, respectively. Finally, an in-silico ligand-protein docking computational method was used to study a possible interaction between tefluthrin and a protein target. RESULTS: Tefluthrin exposure caused severe morphological malformations in zebrafish larvae, including motor neuron abnormalities. The differentially expressed genes were associated with neurotoxicity and metabolic disruption. Lipidomics analysis revealed a disruption in fatty acid, phospholipid, and lysophospholipid recycling. Protein docking modeling suggested that the LPCAT3 enzyme, which recycles lysophospholipids in the Land's cycle, directly interacts with tefluthrin. CONCLUSIONS: Tefluthrin exposure causes morphological and neuronal malformations in zebrafish larvae at nanomolar concentrations. Multi-omics results revealed a potential molecular initiating event i.e., inhibition of LPCAT3, and key events i.e., an altered lysophospholipid to phospholipid ratio, leading to the adverse outcomes of neurotoxicity and metabolic disruption.

Metabolomic analysis of Drosophila melanogaster larvae lacking pyruvate kinase.

Pyruvate kinase (Pyk) is a rate-limiting enzyme that catalyzes the final metabolic reaction in glycolysis. The importance of this enzyme, however, extends far beyond ATP production, as Pyk is also known to regulate tissue growth, cell proliferation, and development. Studies of this enzyme in Drosophila melanogaster are complicated by the fact that the fly genome encodes 6 Pyk paralogs whose functions remain poorly defined. To address this issue, we used sequence distance and phylogenetic approaches to demonstrate that the gene Pyk encodes the enzyme most similar to the mammalian Pyk orthologs, while the other 5 Drosophila Pyk paralogs have significantly diverged from the canonical enzyme. Consistent with this observation, metabolomic studies of 2 different Pyk mutant strains revealed that larvae lacking Pyk exhibit a severe block in glycolysis, with a buildup of glycolytic intermediates upstream of pyruvate. However, our analysis also unexpectedly reveals that pyruvate levels are unchanged in Pyk mutants, indicating that larval metabolism maintains pyruvate pool size despite severe metabolic limitations. Consistent with our metabolomic findings, a complementary RNA-seq analysis revealed that genes involved in lipid metabolism and protease activity are elevated in Pyk mutants, again indicating that loss of this glycolytic enzyme induces compensatory changes in other aspects of metabolism. Overall, our study provides both insight into how Drosophila larval metabolism adapts to disruption of glycolytic metabolism as well as immediate clinical relevance, considering that Pyk deficiency is the most common congenital enzymatic defect in humans.

Metabolomic analysis of Drosophila melanogaster larvae lacking Pyruvate kinase.

Pyruvate kinase (Pyk) is a rate-limiting enzyme that catalyzes the final metabolic reaction in glycolysis. The importance of this enzyme, however, extends far beyond ATP production, as Pyk is also known to regulate tissue growth, cell proliferation, and development. Studies of this enzyme in Drosophila melanogaster , however, are complicated by the fact that the fly genome encodes six Pyk paralogs whose functions remain poorly defined. To address this issue, we used sequence distance and phylogenetic approaches to demonstrate that the gene Pyk encodes the enzyme most similar to the mammalian Pyk orthologs, while the other five Drosophila Pyk paralogs have significantly diverged from the canonical enzyme. Consistent with this observation, metabolomic studies of two different Pyk mutant backgrounds revealed that larvae lacking Pyk exhibit a severe block in glycolysis, with a buildup of glycolytic intermediates upstream of pyruvate. However, our analysis also unexpectedly reveals that steady state pyruvate levels are unchanged in Pyk mutants, indicating that larval metabolism maintains pyruvate pool size despite severe metabolic limitations. Consistent with our metabolomic findings, a complementary RNA-seq analysis revealed that genes involved in lipid metabolism and peptidase activity are elevated in Pyk mutants, again indicating that loss of this glycolytic enzyme induces compensatory changes in other aspects of metabolism. Overall, our study provides both insight into how Drosophila larval metabolism adapts to disruption of glycolytic metabolism as well as immediate clinical relevance, considering that Pyk deficiency is the most common congenital enzymatic defect in humans.

CAGEE: Computational Analysis of Gene Expression Evolution.

Despite the increasing abundance of whole transcriptome data, few methods are available to analyze global gene expression across phylogenies. Here, we present a new software package (Computational Analysis of Gene Expression Evolution [CAGEE]) for inferring patterns of increases and decreases in gene expression across a phylogenetic tree, as well as the rate at which these changes occur. In contrast to previous methods that treat each gene independently, CAGEE can calculate genome-wide rates of gene expression, along with ancestral states for each gene. The statistical approach developed here makes it possible to infer lineage-specific shifts in rates of evolution across the genome, in addition to possible differences in rates among multiple tissues sampled from the same species. We demonstrate the accuracy and robustness of our method on simulated data and apply it to a data set of ovule gene expression collected from multiple self-compatible and self-incompatible species in the genus Solanum to test hypotheses about the evolutionary forces acting during mating system shifts. These comparisons allow us to highlight the power of CAGEE, demonstrating its utility for use in any empirical system and for the analysis of most morphological traits. Our software is available at https://github.com/hahnlab/CAGEE/.

Genome-Wide Profiling of Protein-DNA Interactions with Chromatin Endogenous Cleavage and High-Throughput Sequencing (ChEC-Seq ).

Interactions between regulatory proteins and specific genomic regions are critical for all chromatin-based processes, including transcription, DNA replication, and DNA repair. Genome-wide mapping of such interactions is most commonly performed with chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), but a number of orthogonal methods employing targeted enzymatic activity have also been introduced. We previously described a genome-wide implementation of chromatin endogenous cleavage (ChEC-Seq), wherein a protein of interest is fused to micrococcal nuclease (MNase) to enable targeted, calcium-dependent genomic cleavage. Here, we describe the ChEC-Seq protocol for use in budding yeast though it can be used in other organisms in conjunction with appropriate methods for introduction of an MNase fusion protein.

Architectural Mediator subunits are differentially essential for global transcription in Saccharomyces cerevisiae.

Mediator is a modular coactivator complex involved in the transcription of the majority of RNA polymerase II-regulated genes. However, the degrees to which individual core subunits of Mediator contribute to its activity have been unclear. Here, we investigate the contribution of two essential architectural subunits of Mediator to transcription in Saccharomyces cerevisiae. We show that acute depletion of the main complex scaffold Med14 or the head module nucleator Med17 is lethal and results in global transcriptional downregulation, though Med17 removal has a markedly greater negative effect. Consistent with this, Med17 depletion impairs preinitiation complex (PIC) assembly to a greater extent than Med14 removal. Co-depletion of Med14 and Med17 reduced transcription and TFIIB promoter occupancy similarly to Med17 ablation alone, indicating that the contributions of Med14 and Med17 to Mediator function are not additive. We propose that, while the structural integrity of complete Mediator and the head module are both important for PIC assembly and transcription, the head module plays a greater role in this process and is thus the key functional module of Mediator in this regard.

Mediator Is Essential for Small Nuclear and Nucleolar RNA Transcription in Yeast.

Eukaryotic RNA polymerase II (RNAPII) transcribes mRNA genes and non-protein-coding RNA (ncRNA) genes, including those encoding small nuclear and nucleolar RNAs (sn/snoRNAs). In metazoans, RNAPII transcription of sn/snoRNAs is facilitated by a number of specialized complexes, but no such complexes have been discovered in yeast. It has been proposed that yeast sn/snoRNA and mRNA expression relies on a set of common factors, but the extent to which regulators of mRNA genes function at yeast sn/snoRNA genes is unclear. Here, we investigated a potential role for the Mediator complex, essential for mRNA gene transcription, in sn/snoRNA gene transcription. We found that Mediator maps to sn/snoRNA gene regulatory regions and that rapid depletion of the essential structural subunit Med14 strongly reduces RNAPII and TFIIB occupancy as well as nascent transcription of sn/snoRNA genes. Deletion of Med3 and Med15, subunits of the activator-interacting Mediator tail module, does not affect Mediator recruitment to or RNAPII and TFIIB occupancy of sn/snoRNA genes. Our analyses suggest that Mediator promotes PIC formation and transcription at sn/snoRNA genes, expanding the role of this critical regulator beyond its known functions in mRNA gene transcription and demonstrating further mechanistic similarity between the transcription of mRNA and sn/snoRNA genes.

Analysis of Epstein-Barr virus-regulated host gene expression changes through primary B-cell outgrowth reveals delayed kinetics of latent membrane protein 1-mediated NF-κB activation.

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that dramatically reorganizes host gene expression to immortalize primary B cells. In this study, we analyzed EBV-regulated host gene expression changes following primary B-cell infection, both during initial proliferation and through transformation into lymphoblastoid cell lines (LCLs). While most EBV-regulated mRNAs were changed during the transition from resting, uninfected B cells through initial B-cell proliferation, a substantial number of mRNAs changed uniquely from early proliferation through LCL outgrowth. We identified constitutively and dynamically EBV-regulated biological processes, protein classes, and targets of specific transcription factors. Early after infection, genes associated with proliferation, stress responses, and the p53 pathway were highly enriched. However, the transition from early to long-term outgrowth was characterized by genes involved in the inhibition of apoptosis, the actin cytoskeleton, and NF-κB activity. It was previously thought that the major viral protein responsible for NF-κB activation, latent membrane protein 1 (LMP1), is expressed within 2 days after infection. Our data indicate that while this is true, LCL-level LMP1 expression and NF-κB activity are not evident until 3 weeks after primary B-cell infection. Furthermore, heterologous NF-κB activation during the first week after infection increased the transformation efficiency, while early NF-κB inhibition had no effect on transformation. Rather, inhibition of NF-κB was not toxic to EBV-infected cells until LMP1 levels and NF-κB activity were high. These data collectively highlight the dynamic nature of EBV-regulated host gene expression and support the notion that early EBV-infected proliferating B cells have a fundamentally distinct growth and survival phenotype from that of LCLs.

An ATM/Chk2-mediated DNA damage-responsive signaling pathway suppresses Epstein-Barr virus transformation of primary human B cells.

Epstein-Barr virus (EBV), an oncogenic herpesvirus that causes human malignancies, infects and immortalizes primary human B cells in vitro into indefinitely proliferating lymphoblastoid cell lines, which represent a model for EBV-induced tumorigenesis. The immortalization efficiency is very low, suggesting that an innate tumor suppressor mechanism is operative. We identify the DNA damage response (DDR) as a major component of the underlying tumor suppressor mechanism. EBV-induced DDR activation was not due to lytic viral replication, nor did the DDR marks colocalize with latent episomes. Rather, a transient period of EBV-induced hyperproliferation correlated with DDR activation. Inhibition of the DDR kinases ATM and Chk2 markedly increased transformation efficiency of primary B cells. Further, the viral latent oncoprotein EBNA3C was required to attenuate the EBV-induced DDR. We propose that heightened oncogenic activity in early cell divisions activates a growth-suppressive DDR that is attenuated by viral latency products to induce cell immortalization.