Phylogeographic and Paleoclimatic Modelling Tools Improve Our Understanding of the Biogeographic History of Hierophis viridiflavus (Colubridae).

Phylogeographic and paleoclimatic modelling studies have been combined to infer the role of Pleistocene climatic oscillations as drivers of the genetic structure and distribution of Mediterranean taxa. For the European whip snake, Hierophis viridiflavus, previous studies based on paleoclimatic modelling have depicted a low reliability in the pattern of past climatic suitability across the central Mediterranean Basin, which barely fits the species' genetic structure. In this study, we combined phylogeographic and paleoclimatic modelling tools to improve our understanding of the biogeographic history of H. viridiflavus, particularly extending the sampling and phylogeographic inferences to previously under-sampled regions. Phylogeographic analyses recovered two major clades that diverged at the beginning of the Pleistocene and had diversified in different ways by the late Pleistocene: the east clade (composed of three subclades) and the west clade (with no further structure). Paleoclimatic models highlighted the temperate character of H. viridiflavus, indicating range contractions during both the last inter-glacial and last glacial maximum periods. Range expansions from southern-located climatic refugia likely occurred in the Bølling-Allerød and Middle Holocene periods, which are supported by signals of demographic growth in the west clade and South-East-North subclade. Overall, this work improves our understanding of the historical biogeography of H. viridiflavus, providing further insights into the evolutionary processes that occurred in the Mediterranean Basin hotspot.

Mitogenomics and hidden-trait models reveal the role of phoresy and host shifts in the diversification of parasitoid blister beetles (Coleoptera: Meloidae).

Changes in life history traits are often considered speciation triggers and can have dramatic effects on the evolutionary history of a lineage. Here, we examine the consequences of changes in two life history traits, host-type and phoresy, in the hypermetamorphic blister beetles, Meloidae. Subfamilies Nemognathinae and Meloinae exhibit a complex life cycle involving multiple metamorphoses and parasitoidism. Most genera and tribes are bee-parasitoids, and include phoretic or nonphoretic species, while two tribes feed on grasshopper eggs. These different life strategies are coupled with striking differences in species richness among clades. We generated a mitogenomic phylogeny for Nemognathinae and Meloinae, confirming the monophyly of these two clades, and used the dated phylogeny to explore the association between diversification rates and changes in host specificity and phoresy, using state-dependent speciation and extinction (SSE) models that include the effect of hidden traits. To account for the low taxon sampling, we implemented a phylogenetic-taxonomic approach based on birth-death simulations, and used a Bayesian framework to integrate parameter and phylogenetic uncertainty. Results show that the ancestral hypermetamorphic Meloidae was a nonphoretic bee-parasitoid, and that transitions towards a phoretic bee-parasitoid and grasshopper parasitoidism occurred multiple times. Nonphoretic bee-parasitoid lineages exhibit significantly higher relative extinction and lower diversification rates than phoretic bee-and grasshopper-parasitoids, but no significant differences were found between the latter two strategies. This suggests that Orthopteran host shifts and phoresy contributed jointly to the evolutionary success of the parasitoid meloidae. We also demonstrate that SSE models can be used to identify hidden traits coevolving with the focal trait in driving a lineage's diversification dynamics.

Metagenomics Analysis Reveals an Extraordinary Inner Bacterial Diversity in Anisakids (Nematoda: Anisakidae) L3 Larvae.

L3 larvae of anisakid nematodes are an important problem for the fisheries industry and pose a potential risk for human health by acting as infectious agents causing allergies and as potential vectors of pathogens and microrganisms. In spite of the close bacteria-nematode relationship very little is known of the anisakids microbiota. Fresh fish could be contaminated by bacteria vectored in the cuticle or in the intestine of anisakids when the L3 larvae migrate through the muscles. As a consequence, the bacterial inoculum will be spread, with potential effects on the quality of the fish, and possible clinical effects cannot be discarded. A total of 2,689,113 16S rRNA gene sequences from a total of 113 L3 individuals obtained from fish captured along the FAO 27 fishing area were studied. Bacteria were taxonomically characterized through 1803 representative operational taxonomic units (OTUs) sequences. Fourteen phyla, 31 classes, 52 orders, 129 families and 187 genera were unambiguously identified. We have found as part of microbiome an average of 123 OTUs per L3 individual. Diversity indices (Shannon and Simpson) indicate an extraordinary diversity of bacteria at an OTU level. There are clusters of anisakids individuals (samples) defined by the associated bacteria which, however, are not significantly related to fish hosts or anisakid taxa. This suggests that association or relationship among bacteria in anisakids, exists without the influence of fishes or nematodes. The lack of relationships with hosts of anisakids taxa has to be expressed by the association among bacterial OTUs or other taxonomical levels which range from OTUs to the phylum level. There are significant biological structural associations of microbiota in anisakid nematodes which manifest in clusters of bacteria ranging from phylum to genus level, which could also be an indicator of fish contamination or the geographic zone of fish capture. Actinobacteria, Aquificae, Firmicutes, and Proteobacteria are the phyla whose abundance value discriminate for defining such structures.

Limited long-distance dispersal success in a Western European fairy shrimp evidenced by nuclear and mitochondrial lineage structuring.

Anostraca are known by their ability for long-distance dispersal, but the existence in several species of deep, geographically structured mtDNA lineages suggests their populations are subjected to allopatric differentiation, isolation, and prevalence of local scale dispersion. Tanymastix stagnalis is one of the most widespread species of Anostraca and previous studies revealed an unclear geographical pattern of mtDNA genetic diversity. Here, we analyze populations from the Iberian and Italian Peninsulas, Central Europe, and Scandinavia, with the aim to characterize the patterns of genetic diversity in a spatio-temporal framework using mtDNA and nuclear markers to test gene flow among close populations. For these aims we built a time-calibrated phylogeny and carried out Bayesian phylogeographic analyses using a continuous diffusion model. Our results indicated that T. stagnalis presents a deeply structured genetic diversity, including 7 ancient lineages, some of them even predating the Pleistocene. The Iberian Peninsula harbors high diversity of lineages, with strong isolation and recent absence of gene flow between populations. Dispersal at local scale seems to be the prevailing dispersal mode of T. stagnalis, which exhibits a pattern of isolation-by-distance in the Iberian Peninsula. We remark the vulnerability of most of these lineages, given the limited known geographic distribution of some of them, and the high risk of losing important evolutionary potential for the species.

The classification and phylogenetic status of Jekelius (Reitterius) punctulatus (Jekel, 1866) and Jekelius (Jekelius) brullei (Jekel, 1866) (Coleoptera: Geotrupidae) using molecular data.

The phylogenetic placement of Jekelius brullei (Jekel, 1866) and J. punctulatus (Jekel, 1866) (Coleoptera: Geotrupidae) was assessed using mitochondrial and nuclear molecular data to discern contrasting nomenclatural views provided by López-Colón (1996) and the Catalogue of Palaearctic Coleoptera (Löbl et al. 2006). Our results support both the monophyletic and classification status of the genus Jekelius López-Colón, 1989; and the splitting of the genus into the subgenera Jekelius López-Colón, 1989 and Reitterius López-Colón, 1996. The basal phylogenetic placement of these two species also suggests an oriental origin for Jekelius within the western Palaearctic region. Finally, we include a potential distributional map of Jekelius (Reitterius) punctulatus (Jekel, 1866) based on an exhaustive search of occurrence data.

Proteomic profiling and characterization of differential allergens in the nematodes Anisakis simplex sensu stricto and A. pegreffii.

The parasite species complex Anisakis simplex sensu lato (Anisakis simplex sensu stricto; (A. simplex s.s.), A. pegreffii, A. simplex C) is the main cause of severe anisakiasis (allergy) worldwide and is now an important health matter. In this study, the relationship of this Anisakis species complex and their allergenic capacities is assessed by studying the differences between the two most frequent species (A. simplex s.s., A. pegreffii) and their hybrid haplotype by studying active L3 larvae parasiting Merluccius merluccius. They were compared by 2D gel electrophoresis and parallel Western blot (2DE gels were hybridized with pools of sera from Anisakis allergenic patients). Unambiguous spot differences were detected and protein assignation was made by MALDI-TOF/TOF analysis or de novo sequencing. Seventy-five gel spots were detected and the corresponding proteins were identified. Differentially expressed proteins for A. simplex s.s., A. pegreffii, and their hybrid are described and results are statistically supported. Twenty-eight different allergenic proteins are classified according to different families belonging to different biological functions. These proteins are described for the first time as antigenic and potentially new allergens in Anisakis. Comparative proteomic analyses of allergenic capacities are useful for diagnosis, epidemiological surveys, and clinical research. All MS data have been deposited in the ProteomeXchange with identifier PXD000662 (http://proteomecentral.proteomexchange.org/dataset/PXD000662).