Homology-Free Detection of Transposable Elements Unveils Their Dynamics in Three Ecologically Distinct Rhodnius Species.

Transposable elements (TEs) are widely distributed repetitive sequences in the genomes across the tree of life, and represent an important source of genetic variability. Their distribution among genomes is specific to each lineage. A phenomenon associated with this feature is the sudden expansion of one or several TE families, called bursts of transposition. We previously proposed that bursts of the Mariner family (DNA transposons) contributed to the speciation of Rhodnius prolixus Stål, 1859. This hypothesis motivated us to study two additional species of the R. prolixus complex: Rhodnius montenegrensis da Rosa et al., 2012 and Rhodnius marabaensis Souza et al., 2016, together with a new, de novo annotation of the R. prolixus repeatome using unassembled short reads. Our analysis reveals that the total amount of TEs present in Rhodnius genomes (19% to 23.5%) is three to four times higher than that expected based on the original quantifications performed for the original genome description of R. prolixus. We confirm here that the repeatome of the three species is dominated by Class II elements of the superfamily Tc1-Mariner, as well as members of the LINE order (Class I). In addition to R. prolixus, we also identified a recent burst of transposition of the Mariner family in R. montenegrensis and R. marabaensis, suggesting that this phenomenon may not be exclusive to R. prolixus. Rather, we hypothesize that whilst the expansion of Mariner elements may have contributed to the diversification of the R. prolixus-R. robustus species complex, the distinct ecological characteristics of these new species did not drive the general evolutionary trajectories of these TEs.

Transposition burst of mariner-like elements in the sequenced genome of Rhodnius prolixus.

Transposable elements (TEs) are widespread in insect's genomes. However, there are wide differences in the proportion of the total DNA content occupied by these repetitive sequences in different species. We have analyzed the TEs present in R. prolixus (vector of the Chagas disease) and showed that 3.0% of this genome is occupied by Class II TEs, belonging mainly to the Tc1-mariner superfamily (1.65%) and MITEs (1.84%). Interestingly, most of this genomic content is due to the expansion of two subfamilies belonging to: irritans himar, a well characterized subfamily of mariners, and prolixus1, one of the two novel subfamilies here described. The high amount of sequences in these subfamilies suggests that bursts of transposition occurred during the life cycle of this family. In an attempt to characterize these elements, we performed an in silico analysis of the sequences corresponding to the DDD/E domain of the transposase gene. We performed an evolutionary analysis including network and Bayesian coalescent-based methods in order to infer the dynamics of the amplification, as well as to estimate the time of the bursts identified in these subfamilies. Given our data, we hypothesized that the TE expansions occurred around the time of speciation of R. prolixus around 1.4 mya. This suggestion lays on the "Transposon Model" of TE evolution, in which the members of a TE population that are replicative active are present at multiple loci in the genome, but their replicative potential varies, and of the "Life Cycle Model" that states that when present-day TEs have been involved in amplification bursts, they share an ancestral copy that dates back to this initial amplification.