piRNA clusters as a main source of small RNAs in the animal germline

PIWI subfamily Argonaute proteins and small RNAs bound to them (PIWI interacting RNA, piRNA) control mobilization of transposable elements (TE) in the animal germline. piRNAs are generated by distinct genomic regions termed piRNA clusters. piRNA clusters are often extensive loci enriched in damaged fragments of TEs. New TE integration into piRNA clusters causes production of TE-specific piRNAs and repression of cognate sequences. piRNAs are thought to be generated from long single-stranded precursors encoded by piRNA clusters. Special chromatin structures might be essential to distinguish these genomic loci as a source for piRNAs. In this review, we present recent findings on the structural organization of piRNA clusters and piRNA biogenesis in Drosophila and other organisms, which are important for understanding a key epigenetic mechanism that provides defense against TE expansion.     

piRNA dynamics in divergent zebrafish strains reveal long-lasting maternal influence on zygotic piRNA profiles

Transposable elements (TEs) are mobile genetic elements that can have many deleterious effects on the fitness of their host. The germline-specific PIWI pathway guards the genome against TEs, deriving its specificity from sequence complementarity between PIWI-bound small RNAs (piRNAs) and the TEs. The piRNAs are derived from so-called piRNA clusters. Recent studies have demonstrated that the piRNA repertoire can be adjusted to accommodate recent TE invasions by capturing invading TEs in piRNA loci. Thus far, no information concerning piRNA divergence is available from vertebrates. We present piRNA analyses of two relatively divergent zebrafish strains. We find that significant differences in the piRNA populations have accumulated, most notably among active class I TEs. This divergence can be split into differences in piRNA abundance per element and differences in sense/antisense polarity ratios. In crosses between animals of the different strains, many of these differences are resolved in the progeny. However, some differences remain, often leaning to the maternally contributed piRNA population. These differences can be detected at least two generations later. Our data illustrate, for the first time, the fluidity of piRNA populations in vertebrates and how the established diversity is transmitted to future generations.     

RECURRENT AND RECENT SELECTIVE SWEEPS IN THE piRNA PATHWAY

Uncontrolled transposable element (TE) insertions and excisions can cause chromosome breaks and mutations with dramatic deleterious effects. The PIWI interacting RNA (piRNA) pathway functions as an adaptive TE silencing system during germline development. Several essential piRNA pathway proteins appear to be rapidly evolving, suggesting that TEs and the silencing machinery may be engaged in a classical “evolutionary arms race.” Using a variety of molecular evolutionary and population genetic approaches, we find that the piRNA pathway genes rhino, krimper, and aubergine show patterns suggestive of extensive recurrent positive selection across Drosophila species. We speculate that selection on these proteins reflects crucial roles in silencing unfamiliar elements during vertical and horizontal transmission of TEs into naïve populations and species, respectively.     

Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.     

Illumina TruSeq Synthetic Long-Reads Empower De Novo Assembly and Resolve Complex,Highly-Repetitive Transposable Elements

High-throughput DNA sequencing technologies have revolutionized genomic analysis, including the de novo assembly of whole genomes. Nevertheless, assembly of complex genomes remains challenging, in part due to the presence of dispersed repeats which introduce ambiguity during genome reconstruction. Transposable elements (TEs) can be particularly problematic, especially for TE families exhibiting high sequence identity, high copy number, or complex genomic arrangements. While TEs strongly affect genome function and evolution, most current de novo assembly approaches cannot resolve long, identical, and abundant families of TEs. Here, we applied a novel Illumina technology called TruSeq synthetic long-reads, which are generated through highly-parallel library preparation and local assembly of short read data and which achieve lengths of 1.5–18.5 Kbp with an extremely low error rate (0.03% per base). To test the utility of this technology, we sequenced and assembled the genome of the model organism Drosophila melanogaster (reference genome strain y; cn, bw, sp) achieving an N50 contig size of 69.7 Kbp and covering 96.9% of the euchromatic chromosome arms of the current reference genome. TruSeq synthetic long-read technology enables placement of individual TE copies in their proper genomic locations as well as accurate reconstruction of TE sequences. We entirely recovered and accurately placed 4,229 (77.8%) of the 5,434 annotated transposable elements with perfect identity to the current reference genome. As TEs are ubiquitous features of genomes of many species, TruSeq synthetic long-reads, and likely other methods that generate long-reads, offer a powerful approach to improve de novo assemblies of whole genomes.     

Genomic insights into the adaptation and evolution of the nautilus,an ancient but evolving “living fossil”

The nautilus, commonly known as a “living fossil,” is endangered and may be at risk of extinction. The lack of genomic information hinders a thorough understanding of its biology and evolution, which can shed light on the conservation of this endangered species. Here, we report the first high-quality chromosome-level genome assembly of Nautilus pompilius. The assembled genome size comprised 785.15 Mb. Comparative genomic analyses indicated that transposable elements (TEs) and large-scale genome reorganizations may have driven lineage-specific evolution in the cephalopods. Remarkably, evolving conserved genes and recent TE insertion activities were identified in N. pompilius, and we speculate that these findings reflect the strong adaptability and long-term survival of the nautilus. We also identified gene families that are potentially responsible for specific adaptation and evolution events. Our study provides unprecedented insights into the specialized biology and evolution of N. pompilius, and the results serve as an important resource for future conservation genomics of the nautilus and closely related species.     

TERAD: Extraction of transposable element composition from RADseq data

Transposable elements (TEs) – selfish DNA sequences that can move within the genome – comprise a large proportion of the genomes of many organisms. Although low‐coverage whole‐genome sequencing can be used to survey TE composition, it is noneconomical for species with large quantities of DNA. Here, we utilize restriction‐site associated DNA sequencing (RADSeq) as an alternative method to survey TE composition. First, we demonstrate in silico that double digest restriction‐site associated DNA sequencing (ddRADseq) markers contain the same TE compositions as whole genome assemblies across arthropods. Next, we show empirically using eight Synalpheus snapping shrimp species with large genomes that TE compositions from ddRADseq and low‐coverage whole‐genome sequencing are comparable within and across species. Finally, we develop a new bioinformatic pipeline, TERAD, to extract TE compositions from RADseq data. Our study expands the utility of RADseq to study the repeatome, making comparative studies of genome structure for species with large genomes more tractable and affordable.     

Patterns of piRNA Regulation in Drosophila Revealed through Transposable Element Clade Inference

Transposable elements (TEs) are self-replicating “genetic parasites” ubiquitous to eukaryotic genomes. In addition to conflict between TEs and their host genomes, TEs of the same family are in competition with each other. They compete for the same genomic niches while experiencing the same regime of copy-number selection. This suggests that competition among TEs may favor the emergence of new variants that can outcompete their ancestral forms. To investigate the sequence evolution of TEs, we developed a method to infer clades: collections of TEs that share SNP variants and represent distinct TE family lineages. We applied this method to a panel of 85 Drosophila melanogaster genomes and found that the genetic variation of several TE families shows significant population structure that arises from the population-specific expansions of single clades. We used population genetic theory to classify these clades into younger versus older clades and found that younger clades are associated with a greater abundance of sense and antisense piRNAs per copy than older ones. Further, we find that the abundance of younger, but not older clades, is positively correlated with antisense piRNA production, suggesting a general pattern where hosts preferentially produce antisense piRNAs from recently active TE variants. Together these findings suggest a pattern whereby new TE variants arise by mutation and then increase in copy number, followed by the host producing antisense piRNAs that may be used to silence these emerging variants.     

piRNAs Are Associated with Diverse Transgenerational Effects on Gene and Transposon Expression in a Hybrid Dysgenic Syndrome of D. virilis

Sexual reproduction allows transposable elements (TEs) to proliferate, leading to rapid divergence between populations and species. A significant outcome of divergence in the TE landscape is evident in hybrid dysgenic syndromes, a strong form of genomic incompatibility that can arise when (TE) family abundance differs between two parents. When TEs inherited from the father are absent in the mother''s genome, TEs can become activated in the progeny, causing germline damage and sterility. Studies in Drosophila indicate that dysgenesis can occur when TEs inherited paternally are not matched with a pool of corresponding TE silencing PIWI-interacting RNAs (piRNAs) provisioned by the female germline. Using the D. virilis syndrome of hybrid dysgenesis as a model, we characterize the effects that divergence in TE profile between parents has on offspring. Overall, we show that divergence in the TE landscape is associated with persisting differences in germline TE expression when comparing genetically identical females of reciprocal crosses and these differences are transmitted to the next generation. Moreover, chronic and persisting TE expression coincides with increased levels of genic piRNAs associated with reduced gene expression. Combined with these effects, we further demonstrate that gene expression is idiosyncratically influenced by differences in the genic piRNA profile of the parents that arise though polymorphic TE insertions. Overall, these results support a model in which early germline events in dysgenesis establish a chronic, stable state of both TE and gene expression in the germline that is maintained through adulthood and transmitted to the next generation. This work demonstrates that divergence in the TE profile is associated with diverse piRNA-mediated transgenerational effects on gene expression within populations.     

DeviaTE: Assembly‐free analysis and visualization of mobile genetic element composition

Transposable elements (TEs) are selfish DNA sequences that multiply within host genomes. They are present in most species investigated so far at varying degrees of abundance and sequence diversity. The TE composition may not only vary between but also within species and could have important biological implications. Variation in prevalence among populations may for example indicate a recent TE invasion, whereas sequence variation could indicate the presence of hyperactive or inactive forms. Gaining unbiased estimates of TE composition is thus vital for understanding the evolutionary dynamics of transposons. To this end, we developed DeviaTE, a tool to analyse and visualize TE abundance using Illumina or Sanger sequencing reads. Our tool requires sequencing reads of one or more samples (tissue, individual or population) and consensus sequences of TEs. It generates a table and a visual representation of TE composition. This allows for an intuitive assessment of coverage, sequence divergence, segregating SNPs and indels, as well as the presence of internal and terminal deletions. By contrasting the coverage between TEs and single copy genes, DeviaTE derives unbiased estimates of TE abundance. We show that naive approaches, which do not consider regions spanned by internal deletions, may substantially underestimate TE abundance. Using published data we demonstrate that DeviaTE can be used to study the TE composition within samples, identify clinal variation in TEs, compare TE diversity among species, and monitor TE invasions. Finally we present careful validations with publicly available and simulated data. DeviaTE is implemented in Python and distributed under the GPLv3 ( https://github.com/W-L/deviaTE ).