The evolutionary history of Drosophila buzzatii. XXXIV. The distribution of the retrotransposon Osvaldo in original and colonizing populations

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.

Changes of Osvaldo expression patterns in germline of male hybrids between the species Drosophila buzzatii and Drosophila koepferae

Hybridization between different genomes is a source of genomic instability, sometimes associated with transposable element (TE) mobilization. Previous work showed that hybridization between the species Drosophila buzzatii and Drosophila koepferae induced mobilization of different (TEs), the Osvaldo retrotransposon being the most unstable. However, we ignore the mechanisms involved in this transposition release in interspecific hybrids. In order to disentangle the mechanisms involved in this process, we performed Osvaldo expression studies in somatic and germinal tissues from hybrids and parental species. There was a trend towards increased Osvaldo expression in the somatic tissues of hybrid females and males, which was always significant in males compared to the parental species D. buzzatii but, not in females compared to maternal species D. koepferae. There were massive changes of Osvaldo expression in the testes, which varied depending on the hybrid generation and family. Moreover, Osvaldo hybridization signals, restricted to the apical and primary spermatocyte regions in parents, occupied broader region in the hybrids. In ovaries, there were no significant differences in Osvaldo expression rates between hybrids and the maternal species D. koepferae. The transcript location was restricted to ovarian nurse cells in both parents and hybrids, undetectable in some hybrids. This research highlights first, the existence of putative complex deregulation mechanisms different between sexes and cell types and second, disruption of Osvaldo activity particularly evident in testes from sterile hybrid males. Deeper studies of the total transcriptome in hybrids and parental species are necessary to gain a better knowledge of the TE deregulation pathways in hybrids.

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.

Insertional variability of four transposable elements and population structure of the midge Chironomus riparius (Diptera)

The dipteran Chironomus riparius is found across the entire Palearctic region; its larvae are among the most abundant macroinvertebrates inhabiting inland waterbodies. Chironomid larvae have been extensively used in ecotoxicological and cytogenetic research, but relatively little is known on the population structure of this species. Transposable elements (TEs) are DNA sequences that are capable of autonomous replication; the number and genomic location of TE insertions varies across individuals; this variability is increasingly being used in population studies. Several TEs had been characterized in Chironomids; this enabled the analysis of insertional variability of four different TEs in six natural populations of C. riparius from Italy, Bulgaria and Russia using a PCR-based method, transposon insertion display (TID). The method allows to obtain dominant markers, similar to AFLP. In all populations, TE insertions showed high individual polymorphism, while median copy numbers of the same TEs did not vary between populations. Analysis of molecular variance (AMOVA) detected significant differentiation between populations for three of the TEs; although no correlation between genetic and geographic distances was found, the corresponding population structures were found to be significantly correlated and indicate a degree of isolation by distance. TEs belonging to different classes have different mechanisms of replication, resulting in different transposition rates of mobilization; the finding of mostly concordant population structuring for three of the TEs indicates that population dynamics contributed significantly in shaping the detected insertional polymorphism.

Osvaldo and Isis retrotransposons as markers of the Drosophila buzzatii colonisation in Australia

Background

Transposable elements (TEs) constitute an important source of genetic variability owing to their jumping and regulatory properties, and are considered to drive species evolution. Several factors that are able to induce TE transposition in genomes have been documented (for example environmental stress and inter- and intra-specific crosses) but in many instances the reasons for TE mobilisation have yet to be elucidated. Colonising populations constitute an ideal model for studying TE behaviour and distribution as they are exposed to different environmental and new demographic conditions. In this study, the distribution of two TEs, Osvaldo and Isis, was examined in two colonising populations of D. buzzatii from Australia. Comparing Osvaldo copy numbers between Australian and Old World (reported in previous studies) colonisations provides a valuable tool for elucidating the colonisation process and the effect of new conditions encountered by colonisers on TEs.

Results

The chromosomal distributions of Osvaldo and Isis retrotransposons in two colonising populations of D. buzzatii from Australia revealed sites of high insertion frequency (>10%) and low frequency sites. Comparisons between Osvaldo insertion profiles in colonising populations from the Old World and Australia demonstrate a tendency towards a higher number of highly occupied sites with higher insertion frequency in the Old World than in Australian populations. Tests concerning selection against deleterious TE insertions indicate that Isis is more controlled by purifying selection than Osvaldo. The distribution of both elements on chromosomal arms follows a Poisson distribution and there are non-significant positive correlations between highly occupied sites and chromosomal inversions.

Conclusions

The occupancy profile of Osvaldo and Isis retrotransposons is characterised by the existence of high and low insertion frequency sites in the populations. These results demonstrate that Australian D. buzzatii populations were subjected to a founder effect during the colonisation process. Moreover, there are more sites with high insertion frequency in the Old World colonisation than in the Australian colonisation, indicating a probable stronger bottleneck effect in Australia. The results suggest that selection does not seem to play a major role, compared to demography, in the distribution of transposable elements in the Australian populations.

Distribution of the transposable elements bilbo and gypsy in original and colonizing populations of Drosophila subobscura

Background

Transposable elements (TEs) constitute a substantial amount of all eukaryotic genomes. They induce an important proportion of deleterious mutations by insertion into genes or gene regulatory regions. However, their mutational capabilities are not always adverse but can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution and activity in the genomes of populations under different environmental and demographic regimes, is important to understand their role in species evolution. In this work we study the chromosomal distribution of two TEs, gypsy and bilbo, in original and colonizing populations of Drosophila subobscura to reveal the putative effect of colonization on their insertion profile.

Results

Chromosomal frequency distribution of two TEs in one original and three colonizing populations of D. subobscura, is different. Whereas the original population shows a low insertion frequency in most TE sites, colonizing populations have a mixture of high (frequency ≥ 10%) and low insertion sites for both TEs. Most highly occupied sites are coincident among colonizing populations and some of them are correlated to chromosomal arrangements. Comparisons of TE copy number between the X chromosome and autosomes show that gypsy occupancy seems to be controlled by negative selection, but bilbo one does not.

Conclusion

These results are in accordance that TEs in Drosophila subobscura colonizing populations are submitted to a founder effect followed by genetic drift as a consequence of colonization. This would explain the high insertion frequencies of bilbo and gypsy in coincident sites of colonizing populations. High occupancy sites would represent insertion events prior to colonization. Sites of low frequency would be insertions that occurred after colonization and/or copies from the original population whose frequency is decreasing in colonizing populations. This work is a pioneer attempt to explain the chromosomal distribution of TEs in a colonizing species with high inversion polymorphism to reveal the putative effect of arrangements in TE insertion profiles. In general no associations between arrangements and TE have been found, except in a few cases where the association is very strong. Alternatively, founder drift effects, seem to play a leading role in TE genome distribution in colonizing populations.

Topi, an IS630/Tc1/mariner-type transposable element in the African malaria mosquito, Anopheles gambiae

IS630/Tc1/mariner elements are diverse and widespread within insects. The African malaria mosquito, Anopheles gambiae, contains over 30 families of IS630/Tc1/mariner elements although few have been studied in any detail. To examine the history of Topi elements in An. gambiae populations, Topi elements (n=73) were sampled from five distinct populations of An. gambiae from eastern and western Africa and evaluated with respect to copy number, nucleotide diversity and insertion site-occupancy frequency. Topi 1 and 2 elements were abundant (10-34 per diploid genome) and highly diverse (pi=0.051). Elements from mosquitoes collected in Nigeria were Topi 2 elements and those from mosquitoes collected in Mozambique were Topi 1 elements. Of the 49 Topi transposase open reading frames sequenced none were found to be identical. Intact elements with complete transposase open reading frames were common, although based on insertion site-occupancy frequency data it appeared that genetic drift was the major force acting on these IS630/Tc1/mariner-type elements. Topi 3 elements were not recovered from any of the populations sampled in this study and appear to be rare elements in An. gambiae, possibly due to a recent introduction.

The evolutionary history of Drosophila buzzatii. XXXVI. Molecular structural analysis of Osvaldo retrotransposon insertions in colonizing populations unveils drift effects in founder events

Previous work on transposable element distribution in colonizing populations of Drosophila buzzatii revealed a high frequency of occupancy in several chromosomal sites. Two explanatory hypotheses were advanced: the founder hypothesis, by which founder genetic drift was responsible, and the unstable hypothesis that assigns this unusual distribution to bursts of transposition toward some chromosomal sites. Here, we study the molecular structure of three euchromatic Osvaldo clones isolated from sites occupied at high (A4 and B9) and low frequency (B4) in colonizing populations, to test these hypotheses. Large insertions, duplications, and indels in the Osvaldo coding region and LTR were detected in the A4 clone and a truncated Osvaldo with many substitutions was found in the B9 clone. These altered sequences indicate that the two copies of this retroelement are precolonization insertions. Interestingly, the LTR of the A4 clone and the reverse transcriptase region of B9 show identical sequences in all colonizing populations indicating, most probably, that they are identical by descent. Moreover, Osvaldo is inserted at the same nucleotide site in all colonizing populations. On the other hand an almost identical LTR sequence, except by 1 base deletion, was found in the B4 clone compared to the canonical active Osvaldo element. These results suggest that Osvaldo copies in highly occupied sites are, most probably, identical by descent and strongly favor the founder hypothesis. On the other hand, low-insertion-frequency sites could represent recent transposition events. This work emphasizes the importance of molecular population studies to disentangle the effects of genetic drift and transposition in colonization.

hAT element population genetics in Anopheles gambiae s.l. in Mozambique

Herves is a functional Class II transposable element in Anopheles gambiae belonging to the hAT superfamily of elements. Class II transposable elements are used as gene vectors in this species and are also being considered as genetic drive agents for spreading desirable genes through natural populations as part of an effort to control malaria transmission. In this study, Herves was investigated in populations of Anopheles gambiae s.s., Anopheles arabiensis and Anopheles merus in Mozambique over a period of 2 years. The copy number of Herves within these three species was approximately 5 copies per diploid genome and did not differ among species or between years. Based on the insertion-site occupancy-frequency distribution and existing models of transposable element dynamics, Herves appears to be transpositionally active currently or, at least recently, in all species tested. Ninety-five percent of the individuals within the populations of the three species tested contained intact elements with complete Herves transposase genes and this is consistent with the idea that these elements are currently active.

Molecular characterization and genomic distribution of Isis: a new retrotransposon of Drosophila buzzatii

A new transposable element, Isis, is identified as a LTR retrotransposon in Drosophila buzzatii. DNA sequence analysis shows that Isis contains three long ORFs similar to gag, pol and env genes of retroviruses. The ORF1 exhibits sequence homology to matrix, capsid and nucleocapsid gag proteins and ORF2 encodes a putative protease (PR), a reverse transcriptase (RT), an Rnase H (RH) and an integrase (IN) region. The analysis of a putative env product, encoded by the env ORF3, shows a degenerated protein containing several stop codons. The molecular study of the putative proteins coded by this new element shows striking similarities to both Ulysses and Osvaldo elements, two LTR retrotransposons, present in D. virilis and D. buzzatii, respectively. Comparisons of the predicted Isis RT to several known retrotransposons show strong phylogenetic relationships to gypsy-like elements, particulary to Ulysses retrotransposon. Studies of Isis chromosomal distribution show a strong hybridization signal in centromeric and pericentromeric regions, and a scattered distribution along all chromosomal arms. The existence of insertional polymorphisms between different strains and high molecular weight bands by Southern blot suggests the existence of full-sized copies that have been active recently. The presence of euchromatic insertion sites coincident between Isis and Osvaldo could indicate preferential insertion sites of Osvaldo element into Isis sequence or vice versa. Moreover, the presence of Isis in different species of the buzzatii complex indicates the ancient origin of this element.

Insertional polymorphism of a non-LTR mobile element (NLRCth1) in European populations of Chironomus riparius (Diptera, Chironomidae) as detected by transposon insertion display

The midge Chironomus riparius is distributed all over the Palearctic region and is well characterized both at the morphological and cytogenetic levels. Here we describe a population study based on the insertional polymorphism of the retroposon NLRCth1, by means of a S-SAP (sequence-specific amplification polymorphism) derived technique (transposon insertion display; TID). While a previous study of allozyme polymorphism in Russian samples showed little variability, all the amplicons we identified are polymorphic. Genetic distances between 6 natural populations were calculated according to Nei and did not show a positive correlation with geographic distances. The genetic diversity detected among individuals of a given population was one order of magnitude higher than that among populations. However, the value of phi(ST) was significant (p < 0.001) and indicates that natural populations are more genetically differentiated than random samples of individuals.

Is genome size influenced by colonization of new environments in dipteran species?

Genome size differences are usually attributed to the amplification and deletion of various repeated DNA sequences, including transposable elements (TEs). Because environmental changes may promote modifications in the amount of these repeated sequences, it has been postulated that when a species colonizes new environments this could be followed by an increase in its genome size. We tested this hypothesis by estimating the genome size of geographically distinct populations of Drosophila ananassae, Drosophila malerkotliana, Drosophila melanogaster, Drosophila simulans, Drosophila subobscura, and Zaprionus indianus, all of which have known colonization capacities. There was no strong statistical differences between continents for most species. However, we found that populations of D. melanogaster from east Africa have smaller genomes than more recent populations. For species in which colonization is a recent event, the differences between genome sizes do not thus seem to be related to colonization history. These findings suggest either that genome size is seldom modified in a significant way during colonization or that it takes time for genome size of invading species to change significantly.

Insertion polymorphism of retrotransposable elements in populations of the insular, endemic species Drosophila madeirensis

The insertion site numbers of the retrotransposable elements (TE) 412, gypsy and bilbo were determined in individuals of five distinct natural populations of the endemic species Drosophila madeirensis from the island of Madeira. The TE distributions were compared to those of the paleartic, widespread and phylogenetically closely related species, D. subobscura. In situ hybridization and Southern blots showed that in D. madeirensis the number of insertion sites ranged between 10 and 15, three and six, and 35 and 42 for elements 412, gypsy and bilbo, respectively. The corresponding values for D. subobscura were similar. Two of these elements, 412 and gypsy, had very few insertions in the heterochromatin, unlike bilbo, which displayed a high heterochromatic insertion number. The Southern band polymorphism was very high, leading to within-population variation of 97.2%, whatever the population and the TE concerned. Using the polymorphic TE insertion sites as markers to analyse population structure by AMOVA, adapted for RAPD (Randomly Amplified Polymorphic DNA) data, we found small but significant genetic differences between the populations on Madeira. This slight differentiation, coupled with similar copy numbers for each TE between populations, suggests that the D. madeirensis species consists of a single, only slightly subdivided population. These data also show that insular populations and endemic species of Drosophila can have as many copies of TEs as more widespread species.

The role of the use of different host plants in the maintenance of the inversion polymorphism in the cactophilic Drosophila buzzatii

Inversion polymorphisms often have been associated with fitness variation. Cactophilic Drosophila buzzatii has been used widely for the study of the maintenance of chromosomal variation. The purpose of this paper is to address the relative importance of variable selection regimes associated with the use of three different host cacti and antagonistic pleiotropy in the maintenance of chromosomal variation. Using homokaryotypic stocks derived from several lines homozygous for four second-chromosome arrangements, we show that inversions significantly affect first-instar larva to adult viability (VT), developmental time (DT) and adult thorax length (TL). We also show that the effects of inversions on DT and VT are dependent on the cactus rearing media. The effects of polymorphic gene arrangements on life-history traits suggest the existence of trade-offs between early and late fitness components. The dosage of arrangement 2st, the ancestral gene order, was negatively correlated with DT and TL, whereas flies carrying the derived arrangements 2j and 2jq7 had longer DTs and larger TLs. Arrangements 2st and 2jq7 increased viability, at least in one of the cactus media tested. Our results suggest that environmental heterogeneity, as represented by the use of different cactus hosts and the trade-off between DT and TL, may be involved in the maintenance of the polymorphism. In addition, our data suggest that the chromosomal phylogeny may be decoupled from the evolution of the genes affecting life-history traits linked to the inversion system.