Faecal Microbiota Divergence in Allopatric Populations of Podarcis lilfordi and P. pityusensis, Two Lizard Species Endemic to the Balearic Islands.

Gut microbial communities provide essential functions to their hosts and are known to influence both their ecology and evolution. However, our knowledge of these complex associations is still very limited in reptiles. Here we report the 16S rRNA gene faecal microbiota profiles of two lizard species endemic to the Balearic archipelago (Podarcis lilfordi and P. pityusensis), encompassing their allopatric range of distribution through a noninvasive sampling, as an alternative to previous studies that implied killing specimens of these IUCN endangered and near-threatened species, respectively. Both lizard species showed a faecal microbiome composition consistent with their omnivorous trophic ecology, with a high representation of cellulolytic bacteria taxa. We also identified species-specific core microbiota signatures and retrieved lizard species, islet ascription, and seasonality as the main factors in explaining bacterial community composition. The different Balearic Podarcis populations are characterised by harbouring a high proportion of unique bacterial taxa, thus reinforcing their view as unique and divergent evolutionary entities.

The age and diversification of metacrangonyctid subterranean amphipod crustaceans revisited.

The Metacrangonyctidae are a small family of amphipod crustaceans of marine origin found only in subterranean continental waters. They display a broad but punctuated distribution between the Caribbean and the Arabian Peninsula, with major disjunctions either due to vicariance by plate tectonics or to occurrence of recent episodes of long-distance transoceanic dispersal. We re-examine the phylogeny of the family and the time frame for its diversification using mitochondrial genomes in the light of two key taxa recently discovered, from Oman (Arabian Peninsula) and the Rif area of Morocco, respectively. We also use a novel fossil calibration scheme of the mitogenome phylogeny. Results of previous analyses based on palaeogeographic calibrations are not contradicted by the new approach, with vicariance by plate tectonics remaining as the main explanatory factor for the amphi-Atlantic distribution displayed by this ancient group of subterranean amphipods.

Angiostrongylus cantonensis in North African hedgehogs as vertebrate hosts, Mallorca, Spain, October 2018.

In October 2018, two Atelerix algirus hedgehogs were admitted to the Wildlife Rehabilitation Hospital in Mallorca (Balearic Islands, Spain) with signs of acute neurological disease. Necropsy detected immature, fully developed nematodes in the subarachnoid space of both hedgehogs, including a gravid female worm. DNA-based molecular tools confirmed the nematode as Angiostrongylus cantonensis, an important aetiological agent of eosinophilic meningitis in humans. So far this zoonotic parasite in has not been reported in western European wildlife.

Species delimitation and mitogenome phylogenetics in the subterranean genus Pseudoniphargus (Crustacea: Amphipoda).

The amphi-Atlantic distributions exhibited by many thalassoid stygobiont (obligate subterranean) crustaceans have been explained by fragmentation by plate tectonics of ancestral shallow water marine populations. The amphipod stygobiont genus Pseudoniphargus is distributed across the Mediterranean region but also in the North Atlantic archipelagos of Bermuda, Azores, Madeira and the Canaries. We used species delimitation methods and mitogenome phylogenetic analyses to clarify the species diversity and evolutionary relationships within the genus and timing their diversification. Analyses included samples from the Iberian Peninsula, northern Morocco, the Balearic, Canarian, Azores and Madeira archipelagoes plus Bermuda. In most instances, morphological and molecular-based species delimitation analyses yielded consistent results. Notwithstanding, in a few cases either incipient speciation with no involvement of detectable morphological divergence or species crypticism were the most plausible explanations for the disagreement found between morphological and molecular species delimitations. Phylogenetic analyses based on a robust calibrated mitochondrial tree suggested that Pseudoniphargus lineages have a younger age than for other thalassoid amphipods displaying a disjunct distribution embracing both sides of the Atlantic Ocean. A major split within the family was estimated to occur at the Paleocene, when a lineage from Northern Iberian Peninsula diverged from the rest of pseudoniphargids. Species diversification in the peri-Mediterranean area was deduced to occur in early Miocene to Tortonian times, while in the Atlantic islands it started in the Pliocene. Our results show that the current distribution pattern of Pseudoniphargus resulted from a complex admix of relatively ancient vicariance events and several episodes of long- distance dispersal.

Niche modelling and landscape genetics of the yellow-legged hornet (Vespa velutina): An integrative approach for evaluating central-marginal population dynamics in Europe.

Genetic diversity is an important biological trait for a successful invasion. During the expansion across a new territory, an invasive species may face unprecedented ecological conditions that will determine its demography and genetic diversity. The first record of the yellow-legged hornet (Vespa velutina) in Europe dates back to 2004 in France, from where it has successfully spread through a large territory in the continent, including Italy, Spain and Portugal. Integrative approaches offer a powerful strategy to detect and understand patterns of genetic variation in central and marginal populations. Here, we have analysed the relationship between genetic diversity parameters inferred from 15 V. velutina nuclear DNA microsatellite loci, and geographical and environmental drivers, such as the distance to the introduction focus, environmental suitability and distance to native and invasive niche centroids. Our results revealed a central-marginal dynamic, where allelic richness decreased towards the edge of the expansion range. The low environmental suitability of the territories invaded by marginal populations could prevent a diverse population from establishing and reducing the genetic diversity in populations at the expansion edge. Moreover, Markov chain Monte Carlo analysis showed both geographical and environmental distances were influencing population genetic differentiation. This study highlights the importance of combining genetic analysis with geographical and environmental drivers to understand genetic trends of invasive species to new environment.

DNA metabarcoding the diet of Podarcis lizards endemic to the Balearic Islands.

Dietary studies are essential to unravel the functioning of ecosystems and ultimately to understand biodiversity. This task, which at first may seem simple, becomes especially complex in those cases of omnivorous species with highly variable diets. In this regard, the emergence of next-generation DNA sequencing methodologies represents a powerful tool to address the problem. Here we implement a high-throughput metabarcoding strategy based on the analysis of four molecular markers aimed at sequencing both mitochondrial (animal prey) and chloroplast (diet plants) genome fragments from fecal samples of two lizard species endemic to the Balearic Archipelago (Podarcis lilfordi and P. pityusensis) obtained through non-invasive methods. The results allowed for the characterization of their diets with a high degree of taxonomic detail and have contributed a large number of new trophic records. The reported diets are based mainly on the consumption of arthropods, mollusks and plants from a diversity of taxonomic orders, as well as carrion and marine subsidies. Our analyses also reveal inter- and intra-specific differences both in terms of seasonality and geographical distribution of the sampled lizard populations. These molecular findings provide new insights into the trophic interactions of these threatened endemic lizards in their unique and isolated ecosystems.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea).

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson-Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Phylogenomics of the Hyalella amphipod species-flock of the Andean Altiplano.

Species diversification in ancient lakes has enabled essential insights into evolutionary theory as they embody an evolutionary microcosm compared to continental terrestrial habitats. We have studied the high-altitude amphipods of the Andes Altiplano using mitogenomic, nuclear ribosomal and single-copy nuclear gene sequences obtained from 36 Hyalella genomic libraries, focusing on species of the Lake Titicaca and other water bodies of the Altiplano northern plateau. Results show that early Miocene South American lineages have recently (late Pliocene or early Pleistocene) diversified in the Andes with a striking morphological convergence among lineages. This pattern is consistent with the ecological opportunities (access to unoccupied resources, initial relaxed selection on ecologically-significant traits and low competition) offered by the lacustrine habitats established after the Andean uplift.

Space and Vine Cultivar Interact to Determine the Arbuscular Mycorrhizal Fungal Community Composition.

The interest in the use of microbes as biofertilizers is increasing in recent years as the demands for sustainable cropping systems become more pressing. Although very widely used as biofertilizers, arbuscular mycorrhizal (AM) fungal associations with specific crops have received little attention and knowledge is limited, especially in the case of vineyards. In this study, the AM fungal community associated with soil and roots of a vineyard on Mallorca Island, Spain was characterized by DNA sequencing to resolve the relative importance of grape variety on their diversity and composition. Overall, soil contained a wider AM fungal diversity than plant roots, and this was found at both taxonomic and phylogenetic levels. The major effect on community composition was associated with sample type, either root or soil material, with a significant effect for the variety of the grape. This effect interacted with the spatial distribution of the plants. Such an interaction revealed a hierarchical effect of abiotic and biotic factors in shaping the composition of AM fungal communities. Our results have direct implications for the understanding of plant-fungal assemblages and the potential functional differences across plants in vineyard cropping.

Phylogenetic inference enables reconstruction of a long-overlooked outbreak of almond leaf scorch disease (Xylella fastidiosa) in Europe.

The recent introductions of the bacterium Xylella fastidiosa (Xf) into Europe are linked to the international plant trade. However, both how and when these entries occurred remains poorly understood. Here, we show how almond scorch leaf disease, which affects ~79% of almond trees in Majorca (Spain) and was previously attributed to fungal pathogens, was in fact triggered by the introduction of Xf around 1993 and subsequently spread to grapevines (Pierce's disease). We reconstructed the progression of almond leaf scorch disease by using broad phylogenetic evidence supported by epidemiological data. Bayesian phylogenetic inference predicted that both Xf subspecies found in Majorca, fastidiosa ST1 (95% highest posterior density, HPD: 1990-1997) and multiplex ST81 (95% HPD: 1991-1998), shared their most recent common ancestors with Californian Xf populations associated with almonds and grapevines. Consistent with this chronology, Xf-DNA infections were identified in tree rings dating to 1998. Our findings uncover a previously unknown scenario in Europe and reveal how Pierce's disease reached the continent.

DNA barcoding insect-host plant associations.

Short-sequence fragments ('DNA barcodes') used widely for plant identification and inventorying remain to be applied to complex biological problems. Host-herbivore interactions are fundamental to coevolutionary relationships of a large proportion of species on the Earth, but their study is frequently hampered by limited or unreliable host records. Here we demonstrate that DNA barcodes can greatly improve this situation as they (i) provide a secure identification of host plant species and (ii) establish the authenticity of the trophic association. Host plants of leaf beetles (subfamily Chrysomelinae) from Australia were identified using the chloroplast trnL(UAA) intron as barcodes amplified from beetle DNA extracts. Sequence similarity and phylogenetic analyses provided precise identifications of each host species at tribal, generic and specific levels, depending on the available database coverage in various plant lineages. The 76 species of Chrysomelinae included-more than 10 per cent of the known Australian fauna-feed on 13 plant families, with preference for Australian radiations of Myrtaceae (eucalypts) and Fabaceae (acacias). Phylogenetic analysis of beetles shows general conservation of host association but with rare host shifts between distant plant lineages, including a few cases where barcodes supported two phylogenetically distant host plants. The study demonstrates that plant barcoding is already feasible with the current publicly available data. By sequencing plant barcodes directly from DNA extractions made from herbivorous beetles, strong physical evidence for the host association is provided. Thus, molecular identification using short DNA fragments brings together the detection of species and the analysis of their interactions.

Updated checklist of Balearic leaf beetles (Coleoptera: Chrysomelidae).

The first updated checklist of Balearic leaf beetles (Chrysomelidae) since 1960 is presented here, evincing the presence of 118 species. This estimation is clearly lower than the 141 species reported in the only list available to date (Jolivet, 1953), and the dissimilarity is even more pronounced if we take into account that 22 new species have been added during this period. The possible explanations for these differences are discussed. The main island in the archipelago holds most of the species (Mallorca, 113 spp.), followed by Menorca (71 spp.), Eivissa (39 spp.) and Formentera (19 spp.). Thus, the Gymnesian islands (Mallorca and Menorca) are more species-rich than the Pityusic ones (Eivissa and Formentera). The number of species per island is significantly correlated with their respective areas not only for the Balearic but also for the much larger western Mediterranean islands of Corsica and Sardinia, and these abundances are not related with their nearness to the closest mainland. Among the different subfamilies and tribes, the Balearic flea-beetles (Alticinae) are clearly more prevalent whereas on the contrary, the Clytrini are less represented in comparison with the nearest mainland (Iberian Peninsula). The presented checklist includes four endemic species, Cryptocephalus majoricensis (Mallorca, Menorca and Formentera), C. tramuntanae (Mallorca), Cyrtonus majoricensis (Mallorca) and Timarcha balearica (Mallorca and Menorca). Furthermore, two adventitious species, Monoxia obesula and Epitrix hirtipennis of North American origin, have been reported for the first time in the Balearic Islands, in agreement with previous findings in other Mediterranean countries.

Phylogenetic evidence that both ancient vicariance and dispersal have contributed to the biogeographic patterns of anchialine cave shrimps.

Cave shrimps from the genera Typhlatya, Stygiocaris and Typhlopatsa (Atyidae) are restricted to specialised coastal subterranean habitats or nearby freshwaters and have a highly disconnected distribution (Eastern Pacific, Caribbean, Atlantic, Mediterranean, Madagascar, Australia). The combination of a wide distribution and a limited dispersal potential suggests a large-scale process has generated this geographic pattern. Tectonic plates that fragment ancestral ranges (vicariance) has often been assumed to cause this process, with the biota as passive passengers on continental blocks. The ancestors of these cave shrimps are believed to have inhabited the ancient Tethys Sea, with three particular geological events hypothesised to have led to their isolation and divergence; (1) the opening of the Atlantic Ocean, (2) the breakup of Gondwana, and (3) the closure of the Tethys Seaway. We test the relative contribution of vicariance and dispersal in the evolutionary history of this group using mitochondrial genomes to reconstruct phylogenetic and biogeographic scenarios with fossil-based calibrations. Given that the Australia/Madagascar shrimp divergence postdates the Gondwanan breakup, our results suggest both vicariance (the Atlantic opening) and dispersal. The Tethys closure appears not to have been influential, however we hypothesise that changing marine currents had an important early influence on their biogeography.

The complete mitochondrial genome of the cave shrimp Typhlatya miravetensis (Decapoda, Atyidae) and its systematic position.

The complete mitochondrial genome of Typhlatya miravetensis from one of its only three known localities (Ullal de la Rabla de Miravet, Castellón, Spain) is presented here. The mitogenome is 15,865 bp in length and includes the standard set of two rRNAs, two non-coding regions plus 13 protein-coding genes. The later have been used to perform a phylogenetic analysis together with other Caridea representatives with mitogenome data in GenBank, inferring a close relationship with the Hawaiian volcano shirmp (Halocaridina rubra) within the family Atyidae.

Mitochondrial genome rearrangements at low taxonomic levels: three distinct mitogenome gene orders in the genus Pseudoniphargus (Crustacea: Amphipoda).

A comparison of mitochondrial genomes of three species of the amphipod Pseudoniphargus revealed the occurrence of a surprisingly high level of gene rearrangement involving protein-coding genes that is a rare phenomenon at low taxonomic levels. The three Pseudoniphargus mitogenomes also display a unique gene arrangement with respect to either the presumed Pancrustacean order or those known for other amphipods. Relative long non-coding sequences appear adjacent to the putative breakage points involved in gene rearrangements of protein coding genes. Other details of the newly obtained mitochondrial genomes - e.g., gene content, nucleotide composition and codon usage - are similar to those found in the mitogenomes of other amphipod species studied. They all contain the typical mitochondrial genome set consisting of 13 protein-coding genes, 22 tRNAs, and two rRNAS, as well as a large control region. The secondary structures and characteristics of tRNA and ribosomal mitochondrial genes of these three species are also discussed.

The mitogenome of the amphipod Hyalella lucifugax (Crustacea) and its phylogenetic placement.

The mitogenome of Hyalella lucifugax from Lake Titicaca, obtained using Illumina NGS technology, is described. The mitogenome attains 14,994 bp (although the control region could not be completed) and comprises the standard set of 2 rRNAs, 13 protein-coding genes, and 22 tRNA genes, plus two non-coding regions. A phylogenetic analysis based on the protein-coding mitochondrial genes from representatives from all amphipod genera with available sequences in GenBank recovers the monophyly of H. lucifugax with the superfamily Talitroidea.

High-throughput biodiversity analysis: Rapid assessment of species richness and ecological interactions of Chrysomelidae (Coleoptera) in the tropics.

Biodiversity assessment has been the focus of intense debate and conceptual and methodological advances in recent years. The cultural, academic and aesthetic impulses to recognise and catalogue the diversity in our surroundings, in this case of living objects, is furthermore propelled by the urgency of understanding that we may be responsible for a dramatic reduction of biodiversity, comparable in magnitude to geological mass extinctions. One of the most important advances in this attempt to characterise biodiversity has been incorporating DNA-based characters and molecular taxonomy tools to achieve faster and more efficient species delimitation and identification, even in hyperdiverse tropical biomes. In this assay we advocate for a broad understanding of Biodiversity as the inventory of species in a given environment, but also the diversity of their interactions, with both aspects being attainable using molecular markers and phylogenetic approaches. We exemplify the suitability and utility of this framework for large-scale biodiversity assessment with the results of our ongoing projects trying to characterise the communities of leaf beetles and their host plants in several tropical setups. Moreover, we propose that approaches similar to ours, establishing the inventories of two ecologically inter-related and species-rich groups of organisms, such as insect herbivores and their angiosperm host-plants, can serve as the foundational stone to anchor a comprehensive assessment of diversity, also in tropical environments, by subsequent addition of trophic levels.

New contributions to the molecular systematics and the evolution of host-plant associations in the genus Chrysolina (Coleoptera, Chrysomelidae, Chrysomelinae).

The taxonomic circumscription of the large and diverse leaf beetle genus Chrysolina Motschulsky is not clear, and its discrimination from the closely related genus Oreina Chevrolat has classically been controversial. In addition, the subgeneric arrangement of the species is unstable, and proposals segregating Chrysolina species into new genera have been recently suggested. In this context, the availability of a phylogenetic framework would provide the basis for a stable taxonomic system, but the existing phylogenies are based on few taxa and have low resolution. In the present study we perform a phylogenetic analysis based on mitochondrial (cox1 and rrnL) and nuclear (H3) DNA sequences from a sample of fifty-two Chrysolina species representing almost half of the subgeneric diversity of the group (thirty out of sixty-five subgenera) and most of the morphological, ecological and karyological variation in the genus. In addition, five Oreina species from two subgenera have also been analysed. The resulting phylogeny is used to evaluate some of the most relevant taxonomic hypotheses for Chrysolina, and also to reconstruct its ancestral host plant associations in a Bayesian framework. Our findings support the paraphyly of Chrysolina as currently defined due to the inclusion of Oreina, the monophyly of the Chrysolina (plus Oreina) species including the divergent Chrysolina (Polysticta) vigintimaculata (Clark, 1864), and enable inferences of deep-level evolutionary relationships among the studied subgenera. The plant family Lamiaceae is inferred as the ancestral host of the study group, whose evolution is characterized by continuous host-shifting among pre-existing host plant families. Some Chrysolina clades include mixtures of species with different levels of diet breadth, indicating that niche width has varied through time.