Drosophila Females Undergo Genome Expansion after Interspecific Hybridization

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.

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

Different structural variants of roo retrotransposon are active in Drosophila melanogaster

Retrotransposon roo is one of the most active elements in Drosophila melanogaster. The level of nucleotide diversity between copies of roo is very low but structural variation in the 5'-UTR is considerable. Transposition of roo at high frequency (around 5 × 10^-2 per generation) has been shown previously in the set of mutation accumulation lines named Oviedo. Here we isolated thirteen individual insertions by inverse PCR and sequenced the 5' end of the elements (between 1663 and 2039 nt) including the LTR, the 5'-UTR and a fragment of 661 nucleotides from the ORF, to study whether the new transposed copies come from a unique variant (the master copy model) or different elements are able to move (the transposon model). The elements in the Oviedo lines presented the same structural variants as the reference genome. Different structural variants were active, a behaviour compatible with the "transposon model" in which the copies localized in multiple sites in the genome are able to transpose. At the level of sequence, the copies of roo in our lines are highly similar to the elements in the reference genome. The phylogenetic tree shows a shallow diversification with unsupported nodes denoting that all the elements currently active are very young. This observation together with the great polymorphism in insertion sites implies a rapid turnover of the elements.

Drosophila Interspecific Hybridization Causes A Deregulation of the piRNA Pathway Genes

Almost all eukaryotes have transposable elements (TEs) against which they have developed defense mechanisms. In the Drosophila germline, the main transposable element (TE) regulation pathway is mediated by specific Piwi-interacting small RNAs (piRNAs). Nonetheless, for unknown reasons, TEs sometimes escape cellular control during interspecific hybridization processes. Because the piRNA pathway genes are involved in piRNA biogenesis and TE control, we sequenced and characterized nine key genes from this pathway in Drosophila buzzatii and Drosophila koepferae species and studied their expression pattern in ovaries of both species and their F1 hybrids. We found that gene structure is, in general, maintained between both species and that two genes—armitage and aubergine—are under positive selection. Three genes—krimper, methyltransferase 2, and zucchini—displayed higher expression values in hybrids than both parental species, while others had RNA levels similar to the parental species with the highest expression. This suggests that the overexpression of some piRNA pathway genes can be a primary response to hybrid stress. Therefore, these results reinforce the hypothesis that TE deregulation may be due to the protein incompatibility caused by the rapid evolution of these genes, leading to a TE silencing failure, rather than to an underexpression of piRNA pathway genes.

Transposable Element Misregulation Is Linked to the Divergence between Parental piRNA Pathways in Drosophila Hybrids

Interspecific hybridization is a genomic stress condition that leads to the activation of transposable elements (TEs) in both animals and plants. In hybrids between Drosophila buzzatii and Drosophila koepferae, mobilization of at least 28 TEs has been described. However, the molecular mechanisms underlying this TE release remain poorly understood. To give insight on the causes of this TE activation, we performed a TE transcriptomic analysis in ovaries (notorious for playing a major role in TE silencing) of parental species and their F1 and backcrossed (BC) hybrids. We find that 15.2% and 10.6% of the expressed TEs are deregulated in F1 and BC1 ovaries, respectively, with a bias toward overexpression in both cases. Although differences between parental piRNA (Piwi-interacting RNA) populations explain only partially these results, we demonstrate that piRNA pathway proteins have divergent sequences and are differentially expressed between parental species. Thus, a functional divergence of the piRNA pathway between parental species, together with some differences between their piRNA pools, might be at the origin of hybrid instabilities and ultimately cause TE misregulation in ovaries. These analyses were complemented with the study of F1 testes, where TEs tend to be less expressed than in D. buzzatii. This can be explained by an increase in piRNA production, which probably acts as a defence mechanism against TE instability in the male germline. Hence, we describe a differential impact of interspecific hybridization in testes and ovaries, which reveals that TE expression and regulation are sex-biased.

From Heterochromatin to Long Noncoding RNAs in Drosophila: Expanding the Arena of Gene Function and Regulation

Recent years have witnessed a remarkable interest in exploring the significance of pervasive noncoding transcripts in diverse eukaryotes. Classical cytogenetic studies using the Drosophila model system unraveled the perplexing attributes and "functions" of the "gene"-poor heterochromatin. Recent molecular studies in the fly model are likewise revealing the very diverse and significant roles played by long noncoding RNAs (lncRNAs) in development, gene regulation, chromatin organization, cell and nuclear architecture, etc. There has been a rapid increase in the number of identified lncRNAs, although a much larger number still remains unknown. The diversity of modes of actions and functions of the limited number of Drosophila lncRNAs, which have been examined, already reflects the profound roles of such RNAs in generating and sustaining the biological complexities of eukaryotes. Several of the known Drosophila lncRNAs originate as independent sense or antisense transcripts from promoter or intergenic, intronic, or 5'/3'-UTR regions, while many of them are independent genes that produce only lncRNAs or coding as well as noncoding RNAs. The different lncRNAs affect chromatin organization (local or large-scale pan-chromosomal), transcription, RNA processing/stability, or translation either directly through interaction with their target DNA sequences or indirectly by acting as intermediary molecules for specific regulatory proteins or may act as decoys/sinks, or storage sites for specific proteins or groups of proteins, or may provide a structural framework for the assembly of substructures in nucleus/cytoplasm. It is interesting that many of the "functions" alluded to heterochromatin in earlier cytogenetic studies appear to find correlates with the known subtle as well as far-reaching actions of the different small and long noncoding RNAs. Further studies exploiting the very rich and powerful genetic and molecular resources available for the Drosophila model are expected to unravel the mystery underlying the long reach of ncRNAs.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Dissecting the Satellite DNA Landscape in Three Cactophilic Drosophila Sequenced Genomes

Eukaryote genomes are replete with repetitive DNAs. This class includes tandemly repeated satellite DNAs (satDNA) which are among the most abundant, fast evolving (yet poorly studied) genomic components. Here, we used high-throughput sequencing data from three cactophilic Drosophila species, D. buzzatii, D. seriema, and D. mojavensis, to access and study their whole satDNA landscape. In total, the RepeatExplorer software identified five satDNAs, three previously described (pBuM, DBC-150 and CDSTR198) and two novel ones (CDSTR138 and CDSTR130). Only pBuM is shared among all three species. The satDNA repeat length falls within only two classes, between 130 and 200 bp or between 340 and 390 bp. FISH on metaphase and polytene chromosomes revealed the presence of satDNA arrays in at least one of the following genomic compartments: centromeric, telomeric, subtelomeric, or dispersed along euchromatin. The chromosomal distribution ranges from a single chromosome to almost all chromosomes of the complement. Fiber-FISH and sequence analysis of contigs revealed interspersion between pBuM and CDSTR130 in the microchromosomes of D. mojavensis Phylogenetic analyses showed that the pBuM satDNA underwent concerted evolution at both interspecific and intraspecific levels. Based on RNA-seq data, we found transcription activity for pBuM (in D. mojavensis) and CDSTR198 (in D. buzzatii) in all five analyzed developmental stages, most notably in pupae and adult males. Our data revealed that cactophilic Drosophila present the lowest amount of satDNAs (1.9-2.9%) within the Drosophila genus reported so far. We discuss how our findings on the satDNA location, abundance, organization, and transcription activity may be related to functional aspects.