A snapshot of histone modifications within transposable elements in Drosophila wild type strains

PAD4 controls tumor immunity via restraining the MHC class II machinery in macrophages

Tumor-associated macrophages (TAMs) shape tumor immunity and therapeutic efficacy. However, it is poorly understood whether and how post-translational modifications (PTMs) intrinsically affect the phenotype and function of TAMs. Here, we reveal that peptidylarginine deiminase 4 (PAD4) exhibits the highest expression among common PTM enzymes in TAMs and negatively correlates with the clinical response to immune checkpoint blockade. Genetic and pharmacological inhibition of PAD4 in macrophages prevents tumor progression in tumor-bearing mouse models, accompanied by an increase in macrophage major histocompatibility complex (MHC) class II expression and T cell effector function. Mechanistically, PAD4 citrullinates STAT1 at arginine 121, thereby promoting the interaction between STAT1 and protein inhibitor of activated STAT1 (PIAS1), and the loss of PAD4 abolishes this interaction, ablating the inhibitory role of PIAS1 in the expression of MHC class II machinery in macrophages and enhancing T cell activation. Thus, the PAD4-STAT1-PIAS1 axis is an immune restriction mechanism in macrophages and may serve as a cancer immunotherapy target.

ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

A Quantitative, Genome-Wide Analysis in Drosophila Reveals Transposable Elements' Influence on Gene Expression Is Species-Specific

Transposable elements (TEs) are parasite DNA sequences that are able to move and multiply along the chromosomes of all genomes. They can be controlled by the host through the targeting of silencing epigenetic marks, which may affect the chromatin structure of neighboring sequences, including genes. In this study, we used transcriptomic and epigenomic high-throughput data produced from ovarian samples of several Drosophila melanogaster and Drosophila simulans wild-type strains, in order to finely quantify the influence of TE insertions on gene RNA levels and histone marks (H3K9me3 and H3K4me3). Our results reveal a stronger epigenetic effect of TEs on ortholog genes in D. simulans compared with D. melanogaster. At the same time, we uncover a larger contribution of TEs to gene H3K9me3 variance within genomes in D. melanogaster, which is evidenced by a stronger correlation of TE numbers around genes with the levels of this chromatin mark in D. melanogaster. Overall, this work contributes to the understanding of species-specific influence of TEs within genomes. It provides a new light on the considerable natural variability provided by TEs, which may be associated with contrasted adaptive and evolutionary potentials.

Transposons contribute to the functional diversification of the head, gut, and ovary transcriptomes across Drosophila natural strains

Transcriptomes are dynamic, with cells, tissues, and body parts expressing particular sets of transcripts. Transposable elements (TEs) are a known source of transcriptome diversity; however, studies often focus on a particular type of chimeric transcript, analyze single body parts or cell types, or are based on incomplete TE annotations from a single reference genome. In this work, we have implemented a method based on de novo transcriptome assembly that minimizes the potential sources of errors while identifying a comprehensive set of gene-TE chimeras. We applied this method to the head, gut, and ovary dissected from five Drosophila melanogaster natural strains, with individual reference genomes available. We found that ∼19% of body part-specific transcripts are gene-TE chimeras. Overall, chimeric transcripts contribute a mean of 43% to the total gene expression, and they provide protein domains for DNA binding, catalytic activity, and DNA polymerase activity. Our comprehensive data set is a rich resource for follow-up analysis. Moreover, because TEs are present in virtually all species sequenced to date, their role in spatially restricted transcript expression is likely not exclusive to the species analyzed in this work.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

The population genetic structure approach adds new insights into the evolution of plant LTR retrotransposon lineages

Long terminal repeat retrotransposons (LTR-RTs) in plant genomes differ in abundance, structure and genomic distribution, reflecting the large number of evolutionary lineages. Elements within lineages can be considered populations, in which each element is an individual in its genomic environment. In this way, it would be reasonable to apply microevolutionary analyses to understand transposable element (TE) evolution, such as those used to study the genetic structure of natural populations. Here, we applied a Bayesian method to infer genetic structure of populations together with classical phylogenetic and dating tools to analyze LTR-RT evolution using the monocot Setaria italica as a model species. In contrast to a phylogeny, the Bayesian clusterization method identifies populations by assigning individuals to one or more clusters according to the most probabilistic scenario of admixture, based on genetic diversity patterns. In this work, each LTR-RT insertion was considered to be one individual and each LTR-RT lineage was considered to be a single species. Nine evolutionary lineages of LTR-RTs were identified in the S. italica genome that had different genetic structures with variable numbers of clusters and levels of admixture. Comprehensive analysis of the phylogenetic, clusterization and time of insertion data allowed us to hypothesize that admixed elements represent sequences that harbor ancestral polymorphic sequence signatures. In conclusion, application of microevolutionary concepts in genome evolution studies is suitable as a complementary approach to phylogenetic analyses to address the evolutionary history and functional features of TEs.

TEtools facilitates big data expression analysis of transposable elements and reveals an antagonism between their activity and that of piRNA genes

Over recent decades, substantial efforts have been made to understand the interactions between host genomes and transposable elements (TEs). The impact of TEs on the regulation of host genes is well known, with TEs acting as platforms of regulatory sequences. Nevertheless, due to their repetitive nature it is considerably hard to integrate TE analysis into genome-wide studies. Here, we developed a specific tool for the analysis of TE expression: TEtools. This tool takes into account the TE sequence diversity of the genome, it can be applied to unannotated or unassembled genomes and is freely available under the GPL3 (https://github.com/l-modolo/TEtools). TEtools performs the mapping of RNA-seq data obtained from classical mRNAs or small RNAs onto a list of TE sequences and performs differential expression analyses with statistical relevance. Using this tool, we analyzed TE expression from five Drosophila wild-type strains. Our data show for the first time that the activity of TEs is strictly linked to the activity of the genes implicated in the piwi-interacting RNA biogenesis and therefore fits an arms race scenario between TE sequences and host control genes.

Stress affects the epigenetic marks added by natural transposable element insertions in Drosophila melanogaster

Transposable elements are emerging as an important source of cis-acting regulatory sequences and epigenetic marks that could influence gene expression. However, few studies have dissected the role of specific transposable element insertions on epigenetic gene regulation. Bari-Jheh is a natural transposon that mediates resistance to oxidative stress by adding cis-regulatory sequences that affect expression of nearby genes. In this work, we integrated publicly available ChIP-seq and piRNA data with chromatin immunoprecipitation experiments to get a more comprehensive picture of Bari-Jheh molecular effects. We showed that Bari-Jheh was enriched for H3K9me3 in nonstress conditions, and for H3K9me3, H3K4me3 and H3K27me3 in oxidative stress conditions, which is consistent with expression changes in adjacent genes. We further showed that under oxidative stress conditions, H3K4me3 and H3K9me3 spread to the promoter region of Jheh1 gene. Finally, another insertion of the Bari1 family was associated with increased H3K27me3 in oxidative stress conditions suggesting that Bari1 histone marks are copy-specific. We concluded that besides adding cis-regulatory sequences, Bari-Jheh influences gene expression by affecting the local chromatin state.

Identification of misexpressed genetic elements in hybrids between Drosophila-related species

Crosses between close species can lead to genomic disorders, often considered to be the cause of hybrid incompatibility, one of the initial steps in the speciation process. How these incompatibilities are established and what are their causes remain unclear. To understand the initiation of hybrid incompatibility, we performed reciprocal crosses between two species of Drosophila (D. mojavensis and D. arizonae) that diverged less than 1 Mya. We performed a genome-wide transcriptomic analysis on ovaries from parental lines and on hybrids from reciprocal crosses. Using an innovative procedure of co-assembling transcriptomes, we show that parental lines differ in the expression of their genes and transposable elements. Reciprocal hybrids presented specific gene categories and few transposable element families misexpressed relative to the parental lines. Because TEs are mainly silenced by piwi-interacting RNAs (piRNAs), we hypothesize that in hybrids the deregulation of specific TE families is due to the absence of such small RNAs. Small RNA sequencing confirmed our hypothesis and we therefore propose that TEs can indeed be major players of genome differentiation and be implicated in the first steps of genomic incompatibilities through small RNA regulation.

Parallel evolution of male germline epigenetic poising and somatic development in animals

Changes in gene regulation frequently underlie changes in morphology during evolution, and differences in chromatin state have been linked with changes in anatomical structure and gene expression across evolutionary time. Here we assess the relationship between evolution of chromatin state in germ cells and evolution of gene regulatory programs governing somatic development. We examined the poised (H3K4me3/H3K27me3 bivalent) epigenetic state in male germ cells from five mammalian and one avian species. We find that core genes poised in germ cells from multiple amniote species are ancient regulators of morphogenesis that sit at the top of transcriptional hierarchies controlling somatic tissue development, whereas genes that gain poising in germ cells from individual species act downstream of core poised genes during development in a species-specific fashion. We propose that critical regulators of animal development gained an epigenetically privileged state in germ cells, manifested in amniotes by H3K4me3/H3K27me3 poising, early in metazoan evolution.

Transposable element islands facilitate adaptation to novel environments in an invasive species

Adaptation requires genetic variation, but founder populations are generally genetically depleted. Here we sequence two populations of an inbred ant that diverge in phenotype to determine how variability is generated. Cardiocondyla obscurior has the smallest of the sequenced ant genomes and its structure suggests a fundamental role of transposable elements (TEs) in adaptive evolution. Accumulations of TEs (TE islands) comprising 7.18% of the genome evolve faster than other regions with regard to single-nucleotide variants, gene/exon duplications and deletions and gene homology. A non-random distribution of gene families, larvae/adult specific gene expression and signs of differential methylation in TE islands indicate intragenomic differences in regulation, evolutionary rates and coalescent effective population size. Our study reveals a tripartite interplay between TEs, life history and adaptation in an invasive species.

Genetic variation is key to species evolution. Here the authors sequence two phenotypically distinct populations of the ant Cardiocondyla obscurior, and find accumulations of transposable elements correlating with genetic variation that may have a role in differentiation, adaptation and speciation.

Maternally deposited germline piRNAs silence the tirant retrotransposon in somatic cells

Transposable elements (TEs), whose propagation can result in severe damage to the host genome, are silenced in the animal gonad by Piwi-interacting RNAs (piRNAs). piRNAs produced in the ovaries are deposited in the embryonic germline and initiate TE repression in the germline progeny. Whether the maternally transmitted piRNAs play a role in the silencing of somatic TEs is however unknown. Here we show that maternally transmitted piRNAs from the tirant retrotransposon in Drosophila are required for the somatic silencing of the TE and correlate with an increase in histone H3K9 trimethylation an active tirant copy.