Search

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

The Phasmida genus Candovia comprises nine traditionally recognized species, all endemic to Australia. In this study, Candovia diversity is explored through molecular species-delimitation analyses using the COIFol gene fragment and phylogenetic inferences leveraging seven additional mitochondrial and nuclear loci. Molecular results were integrated with morphological observations, leading us to confirm the already described species and to the delineation of several new taxa and of the new genus Paracandovia. New Candovia species from various parts of Queensland and New South Wales are described and illustrated (C. alata sp. nov., C. byfieldensis sp. nov., C. dalgleishae sp. nov., C. eungellensis sp. nov., C. karasi sp. nov., C. koensi sp. nov. andC. wollumbinensis sp. nov.). New combinations are proposed and species removed from synonymy with the erection of the new genus Paracandovia (P. cercata stat. rev., comb. nov., P. longipes stat. rev., comb. nov., P. pallida comb. nov., P. peridromes comb. nov., P. tenera stat. rev., comb. nov.). Phylogenetic analyses suggest that the egg capitulum may have independently evolved multiple times throughout the evolutionary history of these insects. Furthermore, two newly described species represent the first taxa with fully developed wings in this previously considered apterous clade.

The West-Palaearctic Colobopsis ant populations have long been considered a single species (Colobopsis truncata). We studied the diversity of this species by employing a multidisciplinary approach and combining data from our surveys, museum and private collections, and citizen science platforms. As a result, we have revealed the existence of a second species, which we describe as Colobopsis imitans sp. nov., distributed allopatrically from Co. truncata and living in the Maghreb, Sicily and southern Iberia. While the pigmentation of Co. truncata is reminiscent of Dolichoderus quadripunctatus, that of Co. imitans is similar to Crematogaster scutellaris, with which Co. imitans lives in close spatial association, and whose foraging trails it habitually follows, similar to Camponotus lateralis and other ant-mimicking ants. The isolation between Co. imitans and Co. truncata seems to have occurred relatively recently because of significant, yet not extreme, morphometric differentiation, and to mtDNA polyphyly. Both Co. imitans and Co. truncata appear to employ mimicry of an unpalatable or aggressive ant species as an important defensive strategy; this ‘choice’ of a different model species is motivated by biogeographic reasons and appears to act as a critical evolutionary driver of their diversification.

The salticid spider Habrocestum graecum Dalmas, 1920, until now only known from Greece, is for the first time recorded in Italy. Observations on ecology and behavior are also reported and pictures of its habitus and genitalia are provided. Furthermore, the first DNA barcode sequence for H. graecum is produced and made publicly available. The species has been observed in Puglia, in South-Eastern Italy, and a trans-Ionian dispersal pattern is most likely the cause of its presence both in Greece and Southern Italy, as reported for other taxa with similar distribution in different animal groups.

The myth of a “typical” mitochondrial genome (mtDNA) is a rock-hard belief in the field of genetics, at least for the animal kingdom [...]

Stingless bees often associate with termites, while association with ants is uncommon due to the high costs related to the aggressiveness of ants. Few combative genera of stingless bees can associate with a larger number of hosts, including ants. Here, we describe for the first time the association between a stingless bee (Partamona testacea) and the aggressive predator bullet ant (Paraponera clavata). In the study area, the colonies of P. testacea we spotted were all associated with P. clavata or Atta sp. ants. COI mitochondrial gene sequences of bees associated with bullet ants and leafcutter ants did not show any divergence, thus revealing no evidence that the two nesting strategies represent a case of cryptic speciation and specialisation on specific hosts. Bees are not unarmed respect to the bullet ants; when ants attempted to penetrate in the colony entrance, they were dragged inside the nest and covered by a resin-like substance. Behavioural experiments focused on ants in arenas and focused on bees at their nest entrance proved that the ants are significantly less aggressive toward associated bees and that guard bees are less alarmed when associated ants are presented. We verified by Gas Chromatography and Mass Spectrometry that P. testacea maintains its species-specific cuticular signature in the association with different ant species and that ants and bees possess typical colony signatures. The differential behavioural responses expressed toward associated colony members by both species are likely based on learning these heterospecific cuticular signatures as it occurs in ant parabiotic associations.

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans), the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans), the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively), even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases, moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

The term drywood termites was introduced for those species nesting and feeding in wood with low or no moisture. The term is used as a synonym for the family Kalotermitidae, although not all kalotermitids exhibit the same substrate prefer- ence. There are, in fact, many species living in rotten or damp wood logs (such as Neotermes). Kalotermitidae are placed among the lower ter- mites, a paraphyletic group that includes early- branching termite lineages and is characterized by the presence of symbiotic flagellates in their gut.

Although DNA transposons often generated internal deleted derivatives such as miniature inverted-repeat transposable elements, short internally deleted elements (SIDEs) derived from nonlong terminal-repeat retrotransposons are rare. Here, we found a novel SIDE, named Persaeus, that originated from the chicken repeat 1 (CR1) retrotransposon Zenon and it has been found widespread in Lepidoptera insects. Our findings suggested that Persaeus and the partner Zenon have experienced a transposition burst in their host genomes and the copy number of Persaeus and Zenon in assayed genomes are significantly correlated. Accordingly, the activity though age analysis indicated that the replication wave of Persaeus coincided with that of Zenon. Phylogenetic analyses suggested that Persaeus may have evolved at least four times independently, and that it has been vertically transferred into its host genomes. Together, our results provide new insights into the evolution dynamics of SIDEs and its partner non-LTRs.

Crustaceans of the order Notostraca (Branchiopoda) are distributed worldwide and are known for the remarkable morphological stasis between their extant and Permian fossil species. Moreover, these crustaceans show relevant ecological traits and a wide range of reproductive strategies. However, genomic studies on notostracans are fairly limited. Here, we present the genome sequences of two notostracan taxa, Lepidurus arcticus and Lepidurus apus lubbocki. Taking advantage of the small genome sizes (~0.11pg) of these taxa, genomes were sequenced for one individual per species with one run on the Illumina HiSeq X platform. We finally assembled 73.2Mbp (L.arcticus) and 90.3Mbp (L.apus lubbocki) long genomes. Assemblies cover up to 84% of the estimated genome size, with a gene completeness >97% for both genomes. In total, 13%–16% of the assembled genomes consist of repeats, and based on read mapping, L.apus lubbocki shows a significantly lower transposable element content than L.arcticus. The analysis of 2,376 orthologous genes indicates an ~7% divergence between the two Lepidurus taxa, with a nucleotide substitution rate significantly lower than that of Daphnia taxa. Ka/Ks analysis suggests purifying selection in both branchiopod lineages, raising the question of whether the low substitution rate of Lepidurus is correlated with morphological conservation or is linked to specific biological traits. Our analysis demonstrates that, in these organisms, it is possible to obtain high-quality draft genomes from single individuals with a relatively low sequencing effort. This result makes Lepidurus and Notostraca interesting models for genomic studies at taxonomic, ecological and evolutionary levels.

Inferring the selective forces that orthologous genes underwent across different lineages can help us understand the evolutionary processes that have shaped their extant diversity and the phenotypes they underlie. The most widespread metric to estimate the selection regimes of coding genes—across sites and phylogenies—is the ratio of nonsynonymous to synonymous substitutions (dN/dS, also known as ω). Nowadays, modern sequencing technologies and the large amount of already available sequence data allow the retrieval of thousands of orthologous genes across large numbers of species. Nonetheless, the tools available to explore selection regimes are not designed to automatically process all genes, and their practical usage is often restricted to the single‐copy ones which are found across all species considered (i.e., ubiquitous genes). This approach limits the scale of the analysis to a fraction of single‐copy genes, which can be as low as an order of magnitude in respect to those which are not consistently found in all species considered (i.e., nonubiquitous genes). Here, we present a workflow named BASE that—leveraging the CodeML framework—eases the inference and interpretation of gene selection regimes in the context of comparative genomics. Although a number of bioinformatics tools have already been developed to facilitate this kind of analyses, BASE is the first to be specifically designed to allow the integration of nonubiquitous genes in a straightforward and reproducible manner. The workflow—along with all relevant documentation—is available at github.com/for‐giobbe/BASE., Comparative genomics analyses for selection are often restricted to the subset of genes which is shared by all the species considered (i.e., ubiquitous genes) while those which are not found across some of them (i.e., nonubiquitous genes) are often overlooked, due to the lack of automated approaches to include them. Yet, disregarding such a large portion of genes may potentially conceal important evolutionary processes. For this reason, we developed a novel workflow—named BASE—which allows the integration of nonubiquitous genes in a straightforward and reproducible manner.

Phasmatodea species diversity lies almost entirely within its suborder Euphasmatodea, which exhibits a pantropical distribution and is considered to derive from a recent and rapid evolutionary radiation. To shed light on Euphasmatodea origins and diversification, we assembled the mitogenomes of 17 species from transcriptomic sequencing data and analysed them along with 22 already available Phasmatodea mitogenomes and 33 mitogenomes representing most of the Polyneoptera lineages. Maximum Likelihood and Bayesian Inference approaches retrieved consistent topologies, both showing the widespread conflict between phylogenetic approaches and traditional systematics. We performed a divergence time analysis leveraging ten fossil specimens representative of most polyneopteran lineages: the time tree obtained supports an older radiation of the clade with respect to previous hypotheses. Euphasmatodea diversification is inferred to have started ~ 187 million years ago, suggesting that the Triassic-Jurassic mass extinction and the breakup of Pangea could have contributed to the process. We also investigated Euphasmatodea mitogenomes patterns of dN, dS and dN/dS ratio throughout our time-tree, trying to characterize the selective regime which may have shaped the clade evolution.

The concept that complex ancestral traits can never be recovered after their loss is still widely accepted, despite phylogenetic and molecular approaches suggest instances where phenotypes may have been lost throughout the evolutionary history of a clade and subsequently reverted back in derived lineages. One of the first and most notable examples of such a process is wing evolution in phasmids; this polyneopteran order of insects, which comprises stick and leaf insects, has played a central role in initiating a long-standing debate on the topic. In this study, a novel and comprehensive time tree including over 300 Phasmatodea species is used as a framework for investigating wing evolutionary patterns in the clade. Despite accounting for several possible biases and sources of uncertainty, macroevolutionary analyses consistently revealed multiple reversals to winged states taking place after their loss, and reversibility is coupled with higher species diversification rates. Our findings support a loss of or reduction in wings that occurred in the lineage leading to the extant phasmid most recent common ancestor, and brachyptery is inferred to be an unstable state unless co-opted for nonaerodynamic adaptations. We also explored how different assumptions of wing reversals probability could impact their inference: we found that until reversals are assumed to be over 30 times more unlikely than losses, they are consistently inferred despite uncertainty in tree and model parameters. Our findings demonstrate that wing evolution is a reversible and dynamic process in phasmids and contribute to our understanding of complex trait evolution. [Dollo’s law; Phasmatodea; phylogenetic comparative methods; polyneoptera; reversals; wing.]

The concept that complex ancestral traits can never be re-acquired after their loss has grown popular since its initial formulation and its often referred to as Dollos law. Nonetheless, several macroevolutionary evidences - along with molecular ones - suggest instances where complex phenotypes could have been lost throughout a clade evolutionary history and subsequently reverted to their former state in derived lineages. One of the first and most notable rejection of Dollos law is represented by wing evolution in phasmids: this polyneopteran order of insects - which comprises stick and leaf insects - has played a central role in initiating a long-standing debate on the topic. In this study, a new and comprehensive molecular phylogeny of over 300 Phasmatodea species is used as a framework for investigating wings evolutionary patterns in the clade, taking into consideration several sources of uncertainty and all the methodological recommendations which have been proposed to test Dollos law rejection. Macroevolutionary analyses support a dynamic and reversible evolution of wings, with multiple transitions to ancestral states taking place after their loss. Our findings suggest that neither wings or flight have acted as drivers of Phasmatodea species diversification and that brachyptery is an unstable state, when not co-opted for non-aerodynamic adaptations. We also explored the impact on our results of different assumptions relative to the probability of reversals and losses: we found that until reversals are assumed over 30 times more unlikely than losses, they are consistently inferred despite uncertainty in tree and model parameters. Our findings demonstrate that wings evolution can be a reversible and dynamic process in phasmids and contribute to shape our understanding of complex phenotypes evolution.

The stick-insect genus Pijnackeria includes four diploid bisexual and two polyploid (3n, 4n) parthenogenetic species. Earlier analyses of the tetraploid parthenogen P. hispanica using mitochondrial markers allowed tracing its maternal ancestry to Pijnackeria originis, while no maternal nuclear contribution was found, thus suggesting an androgenetic and hybrid origin. The recently described Pijnackeria recondita—showing, among other features, a specific antennal structure linking it to the tetraploid parthenogen—prompted us to check whether the new species could be theunknown paternal ancestorof P. hispanica. In this work, we use karyology and molecular analysis of the mitochondrial gene cytochrome c oxidase subunit 2 (cox2) and the nuclear gene elongation factor 1 subunit α (ef1-α) to investigate the origin of such a complex tetraploid hybrid parthenogen. The molecular analysis supported P. recondita as being a paternal ancestor of P. hispanica, but also suggested that two more fathering species have to be taken into account: P. barbarae and the unknown paternal ancestor of the triploid hybrid P. masettii. Therefore, P. hispanica is apparently a polyphyletic chimeric androgen, which we propose to indicate as an androgenetic complex. Our data also revealed that P. hispanica is between 1.96 Myr and 3.31 Myr old, making it the oldest parthenogenetic taxon discovered among insects.

Transposable elements are widely distributed within genomes where they may significantly impact their evolution and cell functions. Short interspersed elements (SINEs) are non-autonomous, fast-evolving elements, but some of them carry a highly conserved domain (HCD), whose sequence remained substantially unchanged throughout the metazoan evolution. SINEs carrying the HCD called V are absent in amniote genomes, but V-like sequences were found within the miniature inverted-repeat transposable element (MITE) MER6 in Homo sapiens. In the present work, the genomic distribution and evolution of MER6 are investigated, in order to reconstruct the origin of human V domain and to envisage its possible functional role. The analysis of 85 tetrapod genomes revealed that MER6 and its variant MER6A are found in primates, while only the MER6A variant was found in bats and eulipotyphlans. These MITEs appeared no longer active, in line with literature data on mammalian DNA transposons. Moreover, they appeared to have originated from a Mariner element found in turtles and from a V-SINE from bony fishes. MER6 insertions were found within genes and conserved in mRNAs: in line with previous hypothesis on functional role of HCDs, the MER6 V domain may be important for cell function also in mammals.

Background Transposable elements (TEs) are common and often present with high copy numbers in cellular genomes. Unlike in cellular organisms, TEs were previously thought to be either rare or absent in viruses. Almost all reported TEs display only one or two copies per viral genome. In addition, the discovery of pandoraviruses with genomes up to 2.5-Mb emphasizes the need for biologists to rethink the fundamental nature of the relationship between viruses and cellular life. Results Herein, we performed the first comprehensive analysis of miniature inverted-repeat transposable elements (MITEs) in the 5170 viral genomes for which sequences are currently available. Four hundred and fifty one copies of ten miniature inverted-repeat transposable elements (MITEs) were found and each MITE had reached relatively large copy numbers (some up to 90) in viruses. Eight MITEs belonging to two DNA superfamilies (hobo/Activator/Tam3 and Chapaev–Mirage–CACTA) were for the first time identified in viruses, further expanding the organismal range of these two superfamilies. TEs may play important roles in shaping the evolution of pandoravirus genomes, which were here found to be very rich in MITEs. We also show that putative autonomous partners of seven MITEs are present in the genomes of viral hosts, suggesting that viruses may borrow the transpositional machinery of their cellular hosts’ autonomous elements to spread MITEs and colonize their own genomes. The presence of seven similar MITEs in viral hosts, suggesting horizontal transfers (HTs) as the major mechanism for MITEs propagation. Conclusions Our discovery highlights that TEs contribute to shape genome evolution of pandoraviruses. We concluded that as for cellular organisms, TEs are part of the pandoraviruses’ diverse mobilome. Electronic supplementary material The online version of this article (10.1186/s13100-018-0125-4) contains supplementary material, which is available to authorized users.

Background The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes. Results Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values. Conclusions Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution. Electronic supplementary material The online version of this article (10.1186/s40851-019-0131-5) contains supplementary material, which is available to authorized users.

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

Terminal repeat retrotransposons in miniature (TRIMs) are small non-autonomous LTR retrotransposons consisting of two terminal direct repeats surrounding a short internal domain. The detection and characterization of these elements has been mainly limited to plants. Here we present the frst fnding of a TRIM element in bivalves, and among the frst known in the kingdom Animalia. Class Bivalvia has high ecological and commercial importance in marine ecosystems and aquaculture, and, in recent years, an increasing number of genomic studies has addressed to these organisms. We have identifed biv-TRIM in several bivalve species: Donax trunculus, Ruditapes decussatus, R. philippinarum, Venerupis corrugata, Polititapes rhomboides, Venus verrucosa, Dosinia exoleta, Glycymeris glycymeris, Cerastoderma edule, Magallana gigas, Mytilus galloprovincialis. biv-TRIM has several characteristics typical for this group of elements, exhibiting diferent variations. In addition to canonically structured elements, solo-TDRs and tandem repeats were detected. The presence of this element in the genome of each species is

Termites are eusocial insects that live in colonies characterized by cooperative behavior, where most individuals forego their own reproduction and help to raise the offspring of a few nest members. This may impact on the colony abilities to adapt and exploit the environment. In subterranean termites of the genus Reticulitermes, a new colony is settled by a single, heterosexual pair of winged individuals (primary reproductives; kings and queens); as colony maturate and/or upon founders death, neotenics (secondary reproductive; SRs) develop and contribute to offspring production. This leads to inbreeding as SRs are genetically related. In some species, including the Italian R. lucifugus, primary queens produce secondary queens by parthenogenesis (P); then, secondary queens will mate with the primary king extending the genetic contribution of the primary queen through time (Asexual Queen Succession, AQS), eventually helping to avoid inbreeding. As a consequence of the increasing genetic contribution of the primary queen, AQS colonies show a female-biased sex ratio of winged reproducers. Here we present population genetics and sex ratio analyses of R. lucifugus, and compare results with data on two other species: R. flavipes and R. grassei. Genetic data indicated the presence of AQS in R. lucifugus and its absence in the other two species, accordingly to winged sex ratio observations. Data confirmed that P occurs through a mechanism of terminal fusion. Moreover, winged reproductives resulted all produced by amphigony, confirming that P is only used for secondary queen production. Finally, morphometric and biomass analyses performed on R. lucifugus winged reproductives suggested a larger investment in female sex in AQS colonies. The AQS system has been observed to be scattered among termites, occurring through different mechanisms and leading to different outcomes. This variability still need to be explained in the light of social behavior evolution.

Transposable elements (TEs) are self-replicating, mobile DNA sequences which constitute a significant fraction of eukaryotic genomes. They are generally considered selfish DNA, as their replication and random insertion may have deleterious effects on genome functionalities, although some beneficial effects and evolutionary potential have been recognized. Short interspersed elements (SINEs) are non-autonomous TEs with a modular structure: a small RNA-related head, a body, and a long interspersed element-related tail. Despite their high turnover rate and de novo emergence, the body may retain highly conserved domains (HCDs) shared among divergent SINE families: in metazoans, at least nine HCD-SINEs have been recognized. Data mining on public molecular databases allowed the retrieval of 16 new HCD-SINE families from cnidarian, molluscs, arthropods, and vertebrates. Tracking the ancestry of HCDs on the metazoan phylogeny revealed that some of them date back to the Radiata-Bilateria split. Moreover, phylogenetic and age versus divergence analyses of the most ancient HCDs suggested that long-term vertical inheritance is the rule, with few horizontal transfer events. We suggest that the evolutionary conservation of HCDs may be linked to their potential to serve as recombination hotspots. This indirectly affects host genomes by maintaining active and diverse SINE lineages, whose insertions may impact (either positively or negatively) on the evolution of the genome.

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.

Biological invasions are among key factors of ecological changes, and social insects appear as highly successful invasive animals. Subterranean termites of the holarctic genus Reticulitermes are present in Europe with six native and one invasive (the nearctic R. flavipes) species. The species R. urbis shows a disjunct distribution in the Western Balkans, Eastern Italy and Southern France. Previous molecular and population genetics data suggested that the taxon originated from the Balkans, and that Italian and French populations are invasive, but it is still unknown how many introduction events occurred and from which Balkan source populations. To address these questions, a population genetics analysis was performed on a larger sampling than previous studies, using mitochondrial cytochrome oxidase II and 6 microsatellite markers on 47 colonies collected across the whole distribution area. Mitochondrial analysis confirmed the presence of two major lineages where colonies from Balkans, Italy, and France intermingle. Similarly, microsatellite loci analysis indicated the presence of two genetic clusters, though not corresponding to the two mitochondrial clades, each including colonies from the three sampled areas and with individuals showing mixed cluster membership. Overall, French and Italian populations showed indications of bottleneck (reduced genetic diversity and change of allele frequencies) and do not appear genetically differentiated from the Balkan population. Results presented here support a history of multiple introductions in Italy and France, in a scenario consistent with continuous exchanges between native and invasive areas, as expected along human trades routes.

At least 12 genes (FH, HIF2A, MAX, NF1, RET, SDHA, SDHB, SDHC, SDHD, SDHAF2, TMEM127,andVHL) have been implicated in inherited predisposition to phaeochromocytoma (PCC), paraganglioma (PGL), or head and neck paraganglioma (HNPGL) and a germline mutation may be detected in more than 30% of cases. Knowledge of somatic mutations contributing to PCC/PGL/HNPGL pathogenesis has received less attention though mutations inHRAS, HIF2A, NF1, RET,andVHLhave been reported. To further elucidate the role of somatic mutation in PCC/PGL/HNPGL tumourigenesis, we employed a next generation sequencing strategy to analyse “mutation hotspots” in 50 human cancer genes. Mutations were identified forHRAS(c.37G>C; p.G13R and c.182A>G; p.Q61R) in 7.1% (6/85); forBRAF(c.1799T>A; p.V600E) in 1.2% (1/85) of tumours; and forTP53(c.1010G>A; p.R337H) in 2.35% (2/85) of cases. Twenty-one tumours harboured mutations in inherited PCC/PGL/HNPGL genes and noHRAS, BRAF, orTP53mutations occurred in this group. Combining our data with previous reports ofHRASmutations in PCC/PGL we find that the mean frequency ofHRAS/BRAFmutations in sporadic PCC/PGL is 8.9% (24/269) and in PCC/PGL with an inherited gene mutation 0% (0/148) suggesting thatHRAS/BRAFmutations and inherited PCC/PGL genes mutations might be mutually exclusive. We report the first evidence forBRAFmutations in the pathogenesis of PCC/PGL/HNPGL.

Biological invasions are a key factor in the ecological changes, and social insects are among the most successful invasive animals. Phylogeography and population genetics can provide detailed information on the introduction routes and invasion biology. Reticulitermes urbis is a subterranean termite native of Balkan Peninsula and introduced in Southern France and Eastern Italy. Introductions in these countries probably resulted from anthropic activities and invasive populations currently live in both forests and urban areas. To identify the source population and to infer the minimum number of introductions, we analyzed mitochondrial (COII) and nuclear (6 microsatellites) loci on 7 native and 39 introduced colonies (15 from France; 24 from Italy). Mitochondrial analysis confirmed the presence of two major lineages in which native and introduced populations cluster together, irrespective of the sampling area (native or invasive range). Microsatellite loci analysis identified two genetic clusters, each including individuals from both native and introduced samples. Moreover the analysis of molecular variance evidenced very low genetic differentiation between the three considered areas (native, French or Italian range). Finally, data showed that introduced populations are less variable in both mitochondrial and nuclear markers, suggesting that introductions might have induced a loss of genetic diversity. On the whole, results supported multiple introductions into the two invasive ranges of France and Italy. This is consistent with links between native and invasive areas protracted in time, as expected in cases of human trades routes. Historical data of human settlement and kingdoms might provide explanations for R. urbis modern distribution.

The only native European dry-wood termites belong to the genus Kalotermes. Until recently, these termites were thought to belong to a single species, Kalotermes flavicollis. However, recent studies described two new species: K. italicus (Italy) and K. phoenicae (Cyprus, Middle-East). Moreover, molecular phylogenetic studies revealed divergent K. flavicollis sublineages in Corse-Sardinia and Southern France. As such lineages have been shown to frequently merge in same colonies, it is important to better understand the evolution of this genus in Europe. To widen the analyses on Kalotermes phylogeography, 22 colonies collected in 11 European localities were analyzed together with those obtained from previous studies, sequencing 912 bp of the mitochondrial genome (COI/tRNA-Leu/COII) in two individuals/colony. Phylogenetic analyses indicate that (i) K. italicus spreads from Central Italy to South France, (ii) the Corse- Sardinian lineage can also be found in the Italian peninsula, and (iii) the Southern France lineage extends to the Iberian peninsula. A time-calibrated phylogeny, computed on two mitochondrial markers (COI/tRNA- Leu/COII; 16S) on representative samples, showed that (i) extant European Kalotermes started differentiating ~32 Myr ago, with the separation of K. phoenicae; (ii) K. italicus and K. flavicollis split ~15 Myr ago; (iii) K. flavicollis sublineages diverged during and after the Messinian salinity crisis, ~6 Myr ago. Another European termite genus, Reticulitermes, comprises in the same area seven species/subspecies and all taxa showed marked intraspecific differentiation. Interestingly, recent analyses indicated that Reticulitermes differentiated in half time with respect to Kalotermes. The different nesting ecology and/or habitat constraints may explain this significant difference.

Short interspersed elements (SINEs) are non-autonomous retrotransposons. Although they usually show fast evolutionary rates, in some instances highly conserved domains (HCDs) have been observed in elements with otherwise divergent sequences and from distantly related species. Here, we document the life history of two HCD-SINE families in the elephant shark Callorhinchus milii, one specific to the holocephalan lineage (CmiSINEs) and another one (SacSINE1-CM) with homologous elements in sharks and the coelacanth (SacSINE1s, LmeSINE1s). The analyses of their relationships indicated that these elements share the same 3′-tail, which would have allowed both elements to rise to high copy number by exploiting the C. milii L2-2_CM long interspersed element (LINE) enzymes. Molecular clock analysis on SINE activity in C. milii genome evidenced two replication bursts occurring right after two major events in the holocephalan evolution: the end-Permian mass extinction and the radiation of modern Holocephali. Accordingly, the same analysis on the coelacanth homologous elements, LmeSINE1, identified a replication wave close to the split age of the two extant Latimeria species. The genomic distribution of the studied SINEs pointed out contrasting results: some elements were preferentially sorted out from gene regions, but accumulated in flanking regions, while others appear more conserved within genes. Moreover, data from the C. milii transcriptome suggest that these SINEs could be involved in miRNA biogenesis and may be targets for miRNA-based regulation.

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.

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.