Theoretical and empirical studies have shown differential management of transposable elements in organisms with different reproductive strategies. To investigate this issue, we analysed the R2 retroelement structure and variability in parthenogenetic and bisexual populations of Bacillus rossius stick insects, as well as insertions inheritance in the offspring of parthenogenetic isolates and of crosses. The B. rossius genome hosts a functional (R2Br(fun) ) and a degenerate (R2Br(deg) ) element, their presence correlating with neither reproductive strategies nor population distribution. The median-joining network method indicated that R2Br(fun) duplicates through a multiple source model, while R2Br(deg) is apparently still duplicating via a master gene model. Offspring analyses showed that unisexual and bisexual offspring have a similar number of R2Br-occupied sites. Multiple or recent shifts from gonochoric to parthenogenetic reproduction may explain the observed data. Moreover, insertion frequency spectra show that higher-frequency insertions in unisexual offspring significantly outnumber those in bisexual offspring. This suggests that unisexual offspring eliminate insertions with lower efficiency. A comparison with simulated insertion frequencies shows that inherited insertions in unisexual and bisexual offspring are significantly different from the expectation. On the whole, different mechanisms of R2 elimination in unisexual vs bisexual offspring and a complex interplay between recombination effectiveness, natural selection and time can explain the observed data.

Transposable elements (TEs) are selfish genetic elements whose self-replication is contrasted by the host genome. In this context, host reproductive strategies are predicted to impact on both TEs load and activity. The presence and insertion distribution of the non-LTR retrotransposon R2 was here studied in populations of the strictly bisexual Bacillus grandii maretimi and of the obligatory parthenogenetic Bacillus atticus atticus. Furthermore, data were also obtained from the offspring of selected B. a. atticus females. At the population level, the gonochoric B. g. maretimi showed a significantly higher R2 load than the obligatory parthenogenetic B. a. atticus. The comparison with bisexual and unisexual Bacillus rossius populations showed that their values were higher than those recorded for B. a. atticus and similar, or even higher, than those of B. g. maretimi. Consistently, an R2 load reduction is scored in B. a. atticus offspring even if with a great variance. On the whole, data here produced indicate that in the obligatory unisexual B. a. atticus R2 is active and that mechanisms of molecular turnover are effective. Furthermore, progeny analyses show that, at variance of the facultative parthenogenetic B. rossius, the R2 activity is held at a lower rate. Modeling parental-offspring inheritance, suggests that in B. a. atticus recombination plays a major role in eliminating insertions rather than selection, as previously suggested for unisexual B. rossius progeny, even if in both cases a high variance is observed. In addition to this, mechanisms of R2 silencing or chances of clonal selection cannot be ruled out.

The repetitive DNA content of the stick insect species Bacillus rossius (facultative parthenogenetic), Bacillus grandii (gonochoric), and Bacillus atticus (obligate parthenogenetic) was analyzed through the survey of random genomic libraries roughly corresponding to 0.006% of the genome. By repeat masking, 19 families of transposable elements were identified (two LTR and six non-LTR retrotransposons; 11 DNA transposons). Moreover, a de novo analysis revealed, among the three libraries, the first MITE family observed in polyneopteran genomes. On the whole, transposable element abundance represented 23.3% of the genome in B. rossius, 22.9% in B. atticus, and 18% in B. grandii. Tandem repeat content in the three libraries is much lower: 1.32%, 0.64%, and 1.86% in B. rossius, B. grandii, and B. atticus, respectively. Microsatellites are the most abundant in all species. Minisatellites were only found in B. rossius and B. atticus, and five monomers belonging to the Bag320 satellite family were detected in B. atticus. Assuming the survey provides adequate representation of the relative genome, the obligate parthenogenetic species (B. atticus), compared with the other two species analyzed, does not show a lower transposable element content, as expected from some theoretical and empirical studies.

Horizontal transfer (HT) is an event in which the genetic material is transferred from one species to another, even if distantly related, and it has been demonstrated as a possible essential part of the lifecycle of transposable elements (TEs). However, previous studies on the non-LTR R2 retrotransposon, a metazoan-wide distributed element, indicated its vertical transmission since the Radiata-Bilateria split. Here we present the first possible instances of R2 HT in stick insects of the genus Bacillus (Phasmida). Six R2 elements were characterized in the strictly bisexual subspecies B. grandii grandii, B. grandii benazzii and B. grandii maretimi and in the obligatory parthenogenetic taxon B. atticus. These elements were compared with those previously retrieved in the facultative parthenogenetic species B. rossius. Phylogenetic inconsistencies between element and host taxa, and age versus divergence analyses agree and support at least two HT events. These HT events can be explained by taking into consideration the complex Bacillus reproductive biology, which includes also hybridogenesis, gynogenesis and androgenesis. Through these non-canonical reproductive modes, R2 elements may have been transferred between Bacillus genomes. Our data suggest, therefore, a possible role of hybridization for TEs survival and the consequent reshaping of involved genomes.

No abstract