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.

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 Hyalella species flock of Lake Titicaca (Crustacea: Amphipoda): perspectives and drawbacks of dna-based identification

Lake Titicaca, in the High Andes of Perú and Bolivia, harbours the world’s third most speciose ancient-lake amphipod radiation on record. A minimum of nineteen species of Hyalella derived from at least five independent colonization episodes concentrate in this high altitude water body, although the actual species number present has not yet been established and could be much higher. Herein, we take advantage of the description of three new species (H. krolli, H. gonzalezi, and H. hirsuta) and the re-description of other two (H. solida and H. nefrens) to assess the feasibility of adopting a dna-based identification approach to resolve the magnitude of this highly speciose amphipod assemblage. A Bayesian phylogenetic analysis of the evolutionary relationships among South American Hyalella cox1 haplotypes, including those of four out of the five species dealt with herein, shows a great disagreement between taxonomic units delimited under morphological and genetic data, hampering species identification exclusively based on cox1 dna barcode sequences.

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.

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.

Eduardo Moralejo et al. report a phylogenetic reconstruction tracing the origin and progression of a European outbreak of the almond scorch disease pathogen Xylella fastidiosa (Xf). Their data suggest Xf was introduced into Europe via grafting from infected Californian buds and was subsequently spread by the meadow spittlebug to multiple plant hosts.

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.

Morphological and molecular species boundaries in the Hyalella species flock of Lake Titicaca (Crustacea: Amphipoda)

The Hyalella species diversity in the high-altitude water bodies of the Andean Altiplano is addressed using mitochondrial cox1 sequences and implementing different molecular species delimitation criteria. We have recorded the presence of five major genetic lineages in the Altiplano, of which one seems to be exclusive to Lake Titicaca and nearby areas, whereas the rest occur also in other regions of South America. Eleven out of 36 South American entities diagnosed by molecular delimitation criteria in our study are likely endemic to the Titicaca and neighbouring water bodies. We have detected a remarkable disagreement between morphology and genetic data in the Titicacan Hyalella, with occurrence of several cases of the same morpho-species corresponding to several Molecular Operational Taxonomic Units (MOTUs), some even distantly related, and other instances where a particular MOTU is shared by a morphologically heterogeneous array of species, including species with body smooth and others with body heavily armoured. Species diversification and incongruence between morphological and molecular boundaries within this species assemblage may be associated to the sharp changes in hydrological conditions experienced by the water bodies of the Altiplano in the past, which included dramatic fluctuations in water level and salinity of Lake Titicaca. Such environmental shifts could have triggered rapid morphological changes and ecological differentiation within the Hyalella assemblage, followed by phenotypic convergence among the diverse lineages. Factors such as phenotypic plasticity, incomplete lineage sorting or admixture between divergent lineages might lie also at the root of the morphological-genetic incongruence described herein.

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.

Molecular systematics of Haploginglymus, a genus of subterranean amphipods endemic to the Iberian Peninsula (Amphipoda: Niphargidae)

The molecular systematics of the subterranean amphipod genus Haploginglymus is addressed through the phylogenetic analysis of three DNA gene fragments (nuclear ribosomal 28S and protein- coding Histone 3, plus mitochondrial Cytochrome c Oxidase subunit I). We take advantage of the description of a new species from southern Spain (Haploginglymus geos sp. nov.) to assess the singularity of this genus endemic to the Iberian Peninsula and the inclusion of the morphologically aberrant H. morenoi within Haploginglymus. Our results corroborate the monophyly of the family Niphargidae but shows Niphargus to be paraphyletic as it currently stands, with Haploginglymus appearing nested within it. A strongly supported sister-group relationship between niphargids and the (thalassoid) pseudoniphargids is recovered as well, but we propose the Niphargidae should continue to be considered as a primary limnic group for biogeographic purposes despite its presumed marine derivation. Our findings are in agreement with previous studies that suggest the family Niphargidae originated in the late Cretaceous in the NE Atlantic, from where it eventually expanded across continental Europe.

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.

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.

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.

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.

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.

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.