Amplification of the 1731 LTR retrotransposon in Drosophila melanogaster cultured cells: Origin of neocopies and impact on the genome

Local assembly of long reads enables phylogenomics of transposable elements in a polyploid cell line

Animal cell lines often undergo extreme genome restructuring events, including polyploidy and segmental aneuploidy that can impede de novo whole-genome assembly (WGA). In some species like Drosophila, cell lines also exhibit massive proliferation of transposable elements (TEs). To better understand the role of transposition during animal cell culture, we sequenced the genome of the tetraploid Drosophila S2R+ cell line using long-read and linked-read technologies. WGAs for S2R+ were highly fragmented and generated variable estimates of TE content across sequencing and assembly technologies. We therefore developed a novel WGA-independent bioinformatics method called TELR that identifies, locally assembles, and estimates allele frequency of TEs from long-read sequence data (https://github.com/bergmanlab/telr). Application of TELR to a ∼130x PacBio dataset for S2R+ revealed many haplotype-specific TE insertions that arose by transposition after initial cell line establishment and subsequent tetraploidization. Local assemblies from TELR also allowed phylogenetic analysis of paralogous TEs, which revealed that proliferation of TE families in vitro can be driven by single or multiple source lineages. Our work provides a model for the analysis of TEs in complex heterozygous or polyploid genomes that are recalcitrant to WGA and yields new insights into the mechanisms of genome evolution in animal cell culture.

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 melanogaster retrotransposon and inverted repeat-derived endogenous siRNAs are differentially processed in distinct cellular locations

BACKGROUND

Endogenous small interfering (esi)RNAs repress mRNA levels and retrotransposon mobility in Drosophila somatic cells by poorly understood mechanisms. 21 nucleotide esiRNAs are primarily generated from retrotransposons and two inverted repeat (hairpin) loci in Drosophila culture cells in a Dicer2 dependent manner. Additionally, proteins involved in 3' end processing, such as Symplekin, CPSF73 and CPSR100, have been recently implicated in the esiRNA pathway.

RESULTS

Here we present evidence of overlap between two essential RNA metabolic pathways: esiRNA biogenesis and mRNA 3' end processing. We have identified a nucleus-specific interaction between the essential esiRNA cleavage enzyme Dicer2 (Dcr2) and Symplekin, a component of the core cleavage complex (CCC) required for 3' end processing of all eukaryotic mRNAs. This interaction is mediated by the N-terminal 271 amino acids of Symplekin; CCC factors CPSF73 and CPSF100 do not contact Dcr2. While Dcr2 binds the CCC, Dcr2 knockdown does not affect mRNA 3' end formation. RNAi-depletion of CCC components Symplekin and CPSF73 causes perturbations in esiRNA abundance that correlate with fluctuations in retrotransposon and hairpin esiRNA precursor levels. We also discovered that esiRNAs generated from retrotransposons and hairpins have distinct physical characteristics including a higher predominance of 22 nucleotide hairpin-derived esiRNAs and differences in 3' and 5' base preference. Additionally, retrotransposon precursors and derived esiRNAs are highly enriched in the nucleus while hairpins and hairpin derived esiRNAs are predominantly cytoplasmic similar to canonical mRNAs. RNAi-depletion of either CPSF73 or Symplekin results in nuclear retention of both hairpin and retrotransposon precursors suggesting that polyadenylation indirectly affects cellular localization of Dcr2 substrates.

CONCLUSIONS

Together, these observations support a novel mechanism in which differences in localization of esiRNA precursors impacts esiRNA biogenesis. Hairpin-derived esiRNAs are generated in the cytoplasm independent of Dcr2-Symplekin interactions, while retrotransposons are processed in the nucleus.

Recurrent emergence of structural variants of LTR retrotransposon CsRn1 evolving novel expression strategy and their selective expansion in a carcinogenic liver fluke, Clonorchis sinensis

Autonomous retrotransposons, in which replication and transcription are coupled, encode the essential gag and pol genes as a fusion or separate overlapping form(s) that are expressed in single transcripts regulated by a common upstream promoter. The element-specific expression strategies have driven development of relevant translational recoding mechanisms including ribosomal frameshifting to satisfy the protein stoichiometry critical for the assembly of infectious virus-like particles. Retrotransposons with different recoding strategies exhibit a mosaic distribution pattern across the diverse families of reverse transcribing elements, even though their respective distributions are substantially skewed towards certain family groups. However, only a few investigations to date have focused on the emergence of retrotransposons evolving novel expression strategy and causal genetic drivers of the structural variants. In this study, the bulk of genomic and transcribed sequences of a Ty3/gypsy-like CsRn1 retrotransposon in Clonorchis sinensis were analyzed for the comprehensive examination of its expression strategy. Our results demonstrated that structural variants with single open reading frame (ORF) have recurrently emerged from precedential CsRn1 copies encoding overlapping gag-pol ORFs by a single-nucleotide insertion in an upstream region of gag stop codon. In the parasite genome, some of the newly evolved variants appeared to undergo proliferative burst as active master lineages together with their ancestral copies. The genetic event was similarly observed in Opisthorchis viverrini, the closest neighbor of C. sinensis, whereas the resulting structural variants might have failed to overcome purifying selection and comprised minor remnant copies in the Opisthorchis genome.

Antisense Transcription of Retrotransposons in Drosophila: An Origin of Endogenous Small Interfering RNA Precursors

Movement of transposons causes insertions, deletions, and chromosomal rearrangements potentially leading to premature lethality in Drosophila melanogaster. To repress these elements and combat genomic instability, eukaryotes have evolved several small RNA-mediated defense mechanisms. Specifically, in Drosophila somatic cells, endogenous small interfering (esi)RNAs suppress retrotransposon mobility. EsiRNAs are produced by Dicer-2 processing of double-stranded RNA precursors, yet the origins of these precursors are unknown. We show that most transposon families are transcribed in both the sense (S) and antisense (AS) direction in Dmel-2 cells. LTR retrotransposons Dm297, mdg1, and blood, and non-LTR retrotransposons juan and jockey transcripts, are generated from intraelement transcription start sites with canonical RNA polymerase II promoters. We also determined that retrotransposon antisense transcripts are less polyadenylated than sense. RNA-seq and small RNA-seq revealed that Dicer-2 RNA interference (RNAi) depletion causes a decrease in the number of esiRNAs mapping to retrotransposons and an increase in expression of both S and AS retrotransposon transcripts. These data support a model in which double-stranded RNA precursors are derived from convergent transcription and processed by Dicer-2 into esiRNAs that silence both sense and antisense retrotransposon transcripts. Reduction of sense retrotransposon transcripts potentially lowers element-specific protein levels to prevent transposition. This mechanism preserves genomic integrity and is especially important for Drosophila fitness because mobile genetic elements are highly active.

What makes transposable elements move in the Drosophila genome?

Transposable elements (TEs), by their capacity of moving and inducing mutations in the genome, are considered important drivers of species evolution. The successful invasions of TEs in genomes, despite their mutational properties, are an apparent paradox. TEs' transposition is usually strongly regulated to low value, but in some cases these elements can also show high transposition rates, which has been associated sometimes to changes in environmental conditions. It is evident that factors susceptible to induce transpositions in natural populations contribute to TE perpetuation. Different factors were proposed as causative agents of TE mobilization in a wide range of organisms: biotic and abiotic stresses, inter- and intraspecific crosses and populational factors. However, there is no clear evidence of the factors capable of inducing TE mobilization in Drosophila, and data on laboratory stocks show contradictory results. The aim of this review is to have an update critical revision about mechanisms promoting transposition of TEs in Drosophila, and to provide to the readers a global vision of the dynamics of these genomic elements in the Drosophila genome.

Blanks, a nuclear siRNA/dsRNA-binding complex component, is required for Drosophila spermiogenesis

Small RNAs and a diverse array of protein partners control gene expression in eukaryotes through a variety of mechanisms. By combining siRNA affinity chromatography and mass spectrometry, we have identified the double-stranded RNA-binding domain protein Blanks to be an siRNA- and dsRNA-binding protein from Drosophila S2 cells. We find that Blanks is a nuclear factor that contributes to the efficiency of RNAi. Biochemical fractionation of a Blanks-containing complex shows that the Blanks complex is unlike previously described RNA-induced silencing complexes and associates with the DEAD-box helicase RM62, a protein previously implicated in RNA silencing. In flies, Blanks is highly expressed in testes tissues and is necessary for postmeiotic spermiogenesis, but loss of Blanks is not accompanied by detectable transposon derepression. Instead, genes related to innate immunity pathways are up-regulated in blanks mutant testes. These results reveal Blanks to be a unique component of a nuclear siRNA/dsRNA-binding complex that contributes to essential RNA silencing-related pathways in the male germ line.

The transcriptional diversity of 25 Drosophila cell lines

Drosophila melanogaster cell lines are important resources for cell biologists. Here, we catalog the expression of exons, genes, and unannotated transcriptional signals for 25 lines. Unannotated transcription is substantial (typically 19% of euchromatic signal). Conservatively, we identify 1405 novel transcribed regions; 684 of these appear to be new exons of neighboring, often distant, genes. Sixty-four percent of genes are expressed detectably in at least one line, but only 21% are detected in all lines. Each cell line expresses, on average, 5885 genes, including a common set of 3109. Expression levels vary over several orders of magnitude. Major signaling pathways are well represented: most differentiation pathways are "off" and survival/growth pathways "on." Roughly 50% of the genes expressed by each line are not part of the common set, and these show considerable individuality. Thirty-one percent are expressed at a higher level in at least one cell line than in any single developmental stage, suggesting that each line is enriched for genes characteristic of small sets of cells. Most remarkable is that imaginal disc-derived lines can generally be assigned, on the basis of expression, to small territories within developing discs. These mappings reveal unexpected stability of even fine-grained spatial determination. No two cell lines show identical transcription factor expression. We conclude that each line has retained features of an individual founder cell superimposed on a common "cell line" gene expression pattern.

Genomic distribution of retrotransposons 297, 1731, copia, mdg1 and roo in the Drosophila melanogaster species subgroup

The intragenomic distribution of five retrotransposon families (297, 1731, copia, mdg1 and roo) in the species of the melanogaster complex was studied by comparing results of the Southern blotting technique in males and females with those of in situ hybridization. The degree of structural polymorphism of each family in the different species was also investigated by restriction enzyme analysis. It was found that genomic distribution is a trait that depends on the family and species. The distribution of roo is mainly euchromatic in the four species and 1731 is heterochromatic, but the distribution of families 297, copia and mdg1 is markedly different in the melanogaster and simulans clades. These families were mainly euchromatic in D. melanogaster but heterochromatic in its sibling species. In the simulans clade most copia and mdg1 elements are located on chromosome Y. Differences in genomic distribution are unrelated with structural conservation. The relation of intragenomic distribution to phylogeny, transpositional activity and the role of the host genome are discussed.

MERE1, a low-copy-number copia-type retroelement in Medicago truncatula active during tissue culture

We have identified an active Medicago truncatula copia-like retroelement called Medicago RetroElement1-1 (MERE1-1) as an insertion in the symbiotic NSP2 gene. MERE1-1 belongs to a low-copy-number family in the sequenced Medicago genome. These copies are highly related, but only three of them have a complete coding region and polymorphism exists between the long terminal repeats of these different copies. This retroelement family is present in all M. truncatula ecotypes tested but also in other legume species like Lotus japonicus. It is active only during tissue culture in both R108 and Jemalong Medicago accessions and inserts preferentially in genes.

Epigenomic consequences of immortalized plant cell suspension culture

Plant cells grown in culture exhibit genetic and epigenetic instability. Using a combination of chromatin immunoprecipitation and DNA methylation profiling on tiling microarrays, we have mapped the location and abundance of histone and DNA modifications in a continuously proliferating, dedifferentiated cell suspension culture of Arabidopsis. We have found that euchromatin becomes hypermethylated in culture and that a small percentage of the hypermethylated genes become associated with heterochromatic marks. In contrast, the heterochromatin undergoes dramatic and very precise DNA hypomethylation with transcriptional activation of specific transposable elements (TEs) in culture. High throughput sequencing of small interfering RNA (siRNA) revealed that TEs activated in culture have increased levels of 21-nucleotide (nt) siRNA, sometimes at the expense of the 24-nt siRNA class. In contrast, TEs that remain silent, which match the predominant 24-nt siRNA class, do not change significantly in their siRNA profiles. These results implicate RNA interference and chromatin modification in epigenetic restructuring of the genome following the activation of TEs in immortalized cell culture.

Changes in chromatin structure correlate with transcriptional activity of nucleolar rDNA in polytene chromosomes

Ribosomal DNA genes (rDNA) are found in tandem arrays of hundreds of repeated genes, but only a fraction of these genes are actively transcribed. The regulatory mechanism controlling the transition between active and inactive rDNA in higher eukaryotes is vital for cell survival. Here, we show that the nucleolus from Drosophila salivary gland cells contains two levels of chromatin organization reflecting differences in transcriptional activity: Decondensed chromatin is highly occupied with TATA-box-binding protein (TBP), phosphorylated H3S10, and acetylated H3K14, suggesting that rDNA in decondensed nucleolar areas is actively transcribed. Condensed chromatin lacks TBP, phosphorylated H3S10, or acetylated H3K14 and is enriched in the rDNA retrotransposons R1 and R2. The data show that R1 and R2 retrotransposons are not actively transcribed in salivary glands and may lead to the epigenetic silencing of flanking rDNA genes and that the silencing mechanisms of these sequences might be partially independent of heterochromatin formation by methylation of histone H3 at lysine 9 and binding of heterochromatin protein 1.

Endogenous RNA interference provides a somatic defense against Drosophila transposons

BACKGROUND

Because of the mutagenic consequences of mobile genetic elements, elaborate defenses have evolved to restrict their activity. A major system that controls the activity of transposable elements (TEs) in flies and vertebrates is mediated by Piwi-interacting RNAs (piRNAs), which are approximately 24-30 nucleotide RNAs that are bound by Piwi-class effectors. The piRNA system is thought to provide primarily a germline defense against TE activity.

RESULTS

Here, we describe a second system that represses Drosophila TEs by using endogenous small interfering RNAs (siRNAs), which are 21 nucleotide, 3'-end-modified RNAs that are dependent on Dicer-2 and Argonaute-2. In contrast to piRNAs, we find that the TE-siRNA system is active in somatic tissues, and particularly so in various immortalized cell lines. Analysis of the patterns and properties of TE-derived small RNAs reveals further distinctions between TE regions and genomic loci that are converted into piRNAs and siRNAs, respectively. Finally, functional tests show that many transposon transcripts accumulate to higher levels in cells and animal tissues that are deficient for Dicer-2 or Argonaute-2.

CONCLUSIONS

Drosophila utilizes two small-RNA systems to restrict transposon activity in the germline (mostly via piRNAs) and in the soma (mostly via siRNAs).

Fluorescence in situ hybridization analysis of hobo, mdg1 and Dm412 transposable elements reveals genomic instability following the Drosophila melanogaster genome sequencing

The genome of Drosophila melanogaster strain y cn bw sp has been sequenced and the transposable elements insertion sites have been determined. We hybridized fluorescence-labeled probes directed to the hobo transposon, Dm412 and mdg1 retrotransposons to polytene chromosomes and compared the observed sites to those published in the annotated genome sequence. We observed an almost twofold increase in the number of hobo hybridization sites (46 found as compared to 24 annotated sites). There was no evidence that the hobo transposition rate is slowing over the 10-year period. The patterns of Dm412 and mdg1 sites have changed less dramatically since the time of genome sequencing. Three novel Dm412 hybridization sites were detected while 4 out of 30 annotated sites were missing. Only one additional mdg1 site was found, while 1 out of 29 annotated sites has been lost.