The development of terrestrial ecosystems emerging after glacier retreat.

The global retreat of glaciers is dramatically altering mountain and high-latitude landscapes, with new ecosystems developing from apparently barren substrates1-4. The study of these emerging ecosystems is critical to understanding how climate change interacts with microhabitat and biotic communities and determines the future of ice-free terrains1,5. Here, using a comprehensive characterization of ecosystems (soil properties, microclimate, productivity and biodiversity by environmental DNA metabarcoding6) across 46 proglacial landscapes worldwide, we found that all the environmental properties change with time since glaciers retreated, and that temperature modulates the accumulation of soil nutrients. The richness of bacteria, fungi, plants and animals increases with time since deglaciation, but their temporal patterns differ. Microorganisms colonized most rapidly in the first decades after glacier retreat, whereas most macroorganisms took longer. Increased habitat suitability, growing complexity of biotic interactions and temporal colonization all contribute to the increase in biodiversity over time. These processes also modify community composition for all the groups of organisms. Plant communities show positive links with all other biodiversity components and have a key role in ecosystem development. These unifying patterns provide new insights into the early dynamics of deglaciated terrains and highlight the need for integrated surveillance of their multiple environmental properties5.

Local climate modulates the development of soil nematode communities after glacier retreat.

The worldwide retreat of glaciers is causing a faster than ever increase in ice-free areas that are leading to the emergence of new ecosystems. Understanding the dynamics of these environments is critical to predicting the consequences of climate change on mountains and at high latitudes. Climatic differences between regions of the world could modulate the emergence of biodiversity and functionality after glacier retreat, yet global tests of this hypothesis are lacking. Nematodes are the most abundant soil animals, with keystone roles in ecosystem functioning, but the lack of global-scale studies limits our understanding of how the taxonomic and functional diversity of nematodes changes during the colonization of proglacial landscapes. We used environmental DNA metabarcoding to characterize nematode communities of 48 glacier forelands from five continents. We assessed how different facets of biodiversity change with the age of deglaciated terrains and tested the hypothesis that colonization patterns are different across forelands with different climatic conditions. Nematodes colonized ice-free areas almost immediately. Both taxonomic and functional richness quickly increased over time, but the increase in nematode diversity was modulated by climate, so that colonization started earlier in forelands with mild summer temperatures. Colder forelands initially hosted poor communities, but the colonization rate then accelerated, eventually leveling biodiversity differences between climatic regimes in the long term. Immediately after glacier retreat, communities were dominated by colonizer taxa with short generation time and r-ecological strategy but community composition shifted through time, with increased frequency of more persister taxa with K-ecological strategy. These changes mostly occurred through the addition of new traits instead of their replacement during succession. The effects of local climate on nematode colonization led to heterogeneous but predictable patterns around the world that likely affect soil communities and overall ecosystem development.

Ice-inhabiting species of Bdelloidea Rotifera reveal a pre-Quaternary ancestry in the Arctic cryosphere.

Historical climate data indicate that the Earth has passed through multiple geological periods with much warmer-than-present climates, including epochs of the Miocene (23-5.3 mya BP) with temperatures 3-4°C above present, and more recent interglacial stages of the Quaternary, for example, Marine Isotope Stage 11c (approx. 425-395 ka BP) and Middle Holocene thermal maximum (7.5-4.2 ka BP), during which continental glaciers may have melted entirely. Such warm periods would have severe consequences for ice-obligate fauna in terms of their distribution, biodiversity and population structure. To determine the impacts of these climatic events in the Nordic cryosphere, we surveyed ice habitats throughout mainland Norway and Svalbard ranging from maritime glaciers to continental ice patches (i.e. non-flowing, inland ice subjected to deep freezing overwinter), finding particularly widespread populations of ice-inhabiting bdelloid rotifers. Combined mitochondrial and nuclear DNA sequencing identified approx. 16 undescribed, species-level rotifer lineages that revealed an ancestry predating the Quaternary (> 2.58 mya). These rotifers also displayed robust freeze/thaw tolerance in laboratory experiments. Collectively, these data suggest that extensive ice refugia, comparable with stable ice patches across the contemporary Norwegian landscape, persisted in the cryosphere over geological time, and may have facilitated the long-term survival of ice-obligate Metazoa before and throughout the Quaternary.

Body-size shifts in aquatic and terrestrial urban communities.

Body size is intrinsically linked to metabolic rate and life-history traits, and is a crucial determinant of food webs and community dynamics1,2. The increased temperatures associated with the urban-heat-island effect result in increased metabolic costs and are expected to drive shifts to smaller body sizes 3 . Urban environments are, however, also characterized by substantial habitat fragmentation 4 , which favours mobile species. Here, using a replicated, spatially nested sampling design across ten animal taxonomic groups, we show that urban communities generally consist of smaller species. In addition, although we show urban warming for three habitat types and associated reduced community-weighted mean body sizes for four taxa, three taxa display a shift to larger species along the urbanization gradients. Our results show that the general trend towards smaller-sized species is overruled by filtering for larger species when there is positive covariation between size and dispersal, a process that can mitigate the low connectivity of ecological resources in urban settings 5 . We thus demonstrate that the urban-heat-island effect and urban habitat fragmentation are associated with contrasting community-level shifts in body size that critically depend on the association between body size and dispersal. Because body size determines the structure and dynamics of ecological networks 1 , such shifts may affect urban ecosystem function.

Back to the roots, desiccation and radiation resistances are ancestral characters in bdelloid rotifers.

BACKGROUND: Bdelloid rotifers are micro-invertebrates distributed worldwide, from temperate latitudes to the most extreme areas of the planet like Antarctica or the Atacama Desert. They have colonized any habitat where liquid water is temporarily available, including terrestrial environments such as soils, mosses, and lichens, tolerating desiccation and other types of stress such as high doses of ionizing radiation (IR). It was hypothesized that bdelloid desiccation and radiation resistance may be attributed to their potential ability to repair DNA double-strand breaks (DSBs). Here, these properties are investigated and compared among nine bdelloid species collected from both mild and harsh habitats, addressing the correlation between the ability of bdelloid rotifers to survive desiccation and their capacity to repair massive DNA breakage in a phylogenetically explicit context. Our research includes both specimens isolated from habitats that experience frequent desiccation (at least 1 time per generation), and individuals sampled from habitats that rarely or never experienced desiccation. RESULTS: Our analysis reveals that DNA repair prevails in somatic cells of both desiccation-tolerant and desiccation-sensitive bdelloid species after exposure to X-ray radiation. Species belonging to both categories are able to withstand high doses of ionizing radiation, up to 1000 Gy, without experiencing any negative effects on their survival. However, the fertility of two desiccation-sensitive species, Rotaria macrura and Rotaria rotatoria, was more severely impacted by low doses of radiation than that of desiccation-resistant species. Surprisingly, the radioresistance of desiccation-resistant species is not related to features of their original habitat. Indeed, bdelloids isolated from Atacama Desert or Antarctica were not characterized by a higher radioresistance than species found in more temperate environments. CONCLUSIONS: Tolerance to desiccation and radiation are supported as ancestral features of bdelloid rotifers, with a group of species of the genus Rotaria having lost this trait after colonizing permanent water habitats. Together, our results provide a comprehensive overview of the evolution of desiccation and radiation resistance among bdelloid rotifers.

Metabarcoding data reveal vertical multitaxa variation in topsoil communities during the colonization of deglaciated forelands.

Ice-free areas are expanding worldwide due to dramatic glacier shrinkage and are undergoing rapid colonization by multiple lifeforms, thus representing key environments to study ecosystem development. It has been proposed that the colonization dynamics of deglaciated terrains is different between surface and deep soils but that the heterogeneity between communities inhabiting surface and deep soils decreases through time. Nevertheless, tests of this hypothesis remain scarce, and it is unclear whether patterns are consistent among different taxonomic groups. Here, we used environmental DNA metabarcoding to test whether community diversity and composition of six groups (Eukaryota, Bacteria, Mycota, Collembola, Insecta, and Oligochaeta) differ between the surface (0-5 cm) and deeper (7.5-20 cm) soil at different stages of development and across five Alpine glaciers. Taxonomic diversity increased with time since glacier retreat and with soil evolution. The pattern was consistent across groups and soil depths. For Eukaryota and Mycota, alpha-diversity was highest at the surface. Time since glacier retreat explained more variation of community composition than depth. Beta-diversity between surface and deep layers decreased with time since glacier retreat, supporting the hypothesis that the first 20 cm of soil tends to homogenize through time. Several molecular operational taxonomic units of bacteria and fungi were significant indicators of specific depths and/or soil development stages, confirming the strong functional variation of microbial communities through time and depth. The complexity of community patterns highlights the importance of integrating information from multiple taxonomic groups to unravel community variation in response to ongoing global changes.

Every fifth published metagenome is not available to science.

Have you ever sought to use metagenomic DNA sequences reported in scientific publications? Were you successful? Here, we reveal that metagenomes from no fewer than 20% of the papers found in our literature search, published between 2016 and 2019, were not deposited in a repository or were simply inaccessible. The proportion of inaccessible data within the literature has been increasing year-on-year. Noncompliance with Open Data is best predicted by the scientific discipline of the journal. The number of citations, journal type (e.g., Open Access or subscription journals), and publisher are not good predictors of data accessibility. However, many publications in high-impact factor journals do display a higher likelihood of accessible metagenomic data sets. Twenty-first century science demands compliance with the ethical standard of data sharing of metagenomes and DNA sequence data more broadly. Data accessibility must become one of the routine and mandatory components of manuscript submissions-a requirement that should be applicable across the increasing number of disciplines using metagenomics. Compliance must be ensured and reinforced by funders, publishers, editors, reviewers, and, ultimately, the authors.

Zooplankton as a Transitional Host for Escherichia coli in Freshwater.

This study shows that Escherichia coli can be temporarily enriched in zooplankton under natural conditions and that these bacteria can belong to different phylogroups and sequence types (STs), including environmental, clinical, and animal isolates. We isolated 10 E. coli strains and sequenced the genomes of two of them. Phylogenetically, the two isolates were closer to strains isolated from poultry meat than to freshwater E. coli, albeit their genomes were smaller than those of the poultry isolates. After isolation and fluorescent protein tagging of strains ED1 and ED157, we show that Daphnia sp. can take up these strains and release them alive again, thus becoming a temporary host for E. coli. In a chemostat experiment, we show that this association does not prolong bacterial long-term survival, but at low abundances it also does not significantly reduce bacterial numbers. We demonstrate that E. coli does not belong to the core microbiota of Daphnia, suffers from competition by the natural Daphnia microbiota, but can profit from its carapax to survive in water. All in all, this study suggests that the association of E. coli with Daphnia is only temporary, but the cells are viable therein, and this might allow encounters with other bacteria for genetic exchange and potential genomic adaptation to the freshwater environment. IMPORTANCE The contamination of freshwater with feces-derived bacteria is a major concern regarding drinking water acquisition and recreational activities. Ecological interactions promoting their persistence are still very scarcely studied. This study, which analyses the survival of E. coli in the presence of zooplankton, is thus of ecological and water safety relevance.

Freshwater zooplankton microbiome composition is highly flexible and strongly influenced by the environment.

The association with microbes in plants and animals is known to be beneficial for host's survival and fitness, but the generality of the effect of the microbiome is still debated. For some animals, similarities in microbiome composition reflect taxonomic relatedness of the hosts, a pattern termed phylosymbiosis. The mechanisms behind the pattern could be due to co-evolution and/or to correlated ecological constraints. General conclusions are hampered by the fact that available knowledge is highly dominated by microbiomes from model species. Here, we addressed the issue of the generality of phylosymbiosis by analysing the species-specificity of microbiomes across different species of freshwater zooplankton, including rotifers, cladocerans, and copepods, coupling field surveys and experimental manipulations. We found that no signal of phylosymbiosis was present, and that the proportion of "core" microbial taxa, stable and consistent within each species, was very low. Changes in food and temperature under laboratory experimental settings revealed that the microbiome of freshwater zooplankton is highly flexible and can be influenced by the external environment. Thus, the role of co-evolution, strict association, and interaction with microbes within the holobiont concept highlighted for vertebrates, corals, sponges, and other animals does not seem to be supported for all animals, at least not for freshwater zooplankton. Zooplankton floats in the environment where both food and bacteria that can provide help in digesting such food are available. In addition, there is probably redundancy for beneficial bacterial functions in the environment, not allowing a strict host-microbiome association to originate and persist.

Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots.

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

Comparative genomics of bdelloid rotifers: Insights from desiccating and nondesiccating species.

Bdelloid rotifers are a class of microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasised desiccation as the cause of several unusual genomic features in this group. Here, we present high-quality whole-genome sequences of 3 bdelloid species: Rotaria macrura and R. magnacalcarata, which are both desiccation intolerant, and Adineta ricciae, which is desiccation tolerant. In combination with the published assembly of A. vaga, which is also desiccation tolerant, we apply a comparative genomics approach to evaluate the potential effects of desiccation tolerance and asexuality on genome evolution in bdelloids. We find that ancestral tetraploidy is conserved among all 4 bdelloid species, but homologous divergence in obligately aquatic Rotaria genomes is unexpectedly low. This finding is contrary to current models regarding the role of desiccation in shaping bdelloid genomes. In addition, we find that homologous regions in A. ricciae are largely collinear and do not form palindromic repeats as observed in the published A. vaga assembly. Consequently, several features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred nonmetazoan genes in both desiccating and nondesiccating bdelloid species and show that this unusual feature is not shared by other animal phyla, even those with desiccation-tolerant representatives. These comparisons call into question the proposed role of desiccation in mediating horizontal genetic transfer.

Effects of Ailanthus altissima Invasion and Removal on High-Biodiversity Mediterranean Grasslands.

Ailanthus altissima is one of the worst invasive plants in Europe in several habitat types, including high-biodiversity grasslands. The aim of this work was to evaluate the impact of the invasive A. altissima on high-biodiversity grassland vegetation and the effects of its removal on the recovery of native plant communities. The study area was within the Alta Murgia National Park (SE Italy). Seventeen vegetation quadrats were sampled in invaded grasslands and nine quadrats were sampled in nearby uninvaded areas. A. altissima was removed from six quadrats, which were sampled for two years after plant removal. Cluster analysis and non-metric multidimensional scaling ordination were used to identify and visualize the general vegetation pattern. Generalised Linear Models with different error structures were used to analyse the effects of A. altissima on native grasslands and vegetation recovery after removal. Results showed that the invasion of A. altissima changed drastically the community composition, reduced plant richness and diversity. Invaded stands had a greater presence of ruderal and widely distributed taxa, as opposed to a lesser presence of endemic and Mediterranean ones. The differences in the community composition between invaded and uninvaded quadrats became clearly detectable when A. altissima plants exceeded a threshold of 1 m of height and 50% of coverage. After A. altissima removal, the recovery of the grassland community was not completely achieved after two years.

Urbanization drives cross-taxon declines in abundance and diversity at multiple spatial scales.

The increasing urbanization process is hypothesized to drastically alter (semi-)natural environments with a concomitant major decline in species abundance and diversity. Yet, studies on this effect of urbanization, and the spatial scale at which it acts, are at present inconclusive due to the large heterogeneity in taxonomic groups and spatial scales at which this relationship has been investigated among studies. Comprehensive studies analysing this relationship across multiple animal groups and at multiple spatial scales are rare, hampering the assessment of how biodiversity generally responds to urbanization. We studied aquatic (cladocerans), limno-terrestrial (bdelloid rotifers) and terrestrial (butterflies, ground beetles, ground- and web spiders, macro-moths, orthopterans and snails) invertebrate groups using a hierarchical spatial design, wherein three local-scale (200 m × 200 m) urbanization levels were repeatedly sampled across three landscape-scale (3 km × 3 km) urbanization levels. We tested for local and landscape urbanization effects on abundance and species richness of each group, whereby total richness was partitioned into the average richness of local communities and the richness due to variation among local communities. Abundances of the terrestrial active dispersers declined in response to local urbanization, with reductions up to 85% for butterflies, while passive dispersers did not show any clear trend. Species richness also declined with increasing levels of urbanization, but responses were highly heterogeneous among the different groups with respect to the richness component and the spatial scale at which urbanization impacts richness. Depending on the group, species richness declined due to biotic homogenization and/or local species loss. This resulted in an overall decrease in total richness across groups in urban areas. These results provide strong support to the general negative impact of urbanization on abundance and species richness within habitat patches and highlight the importance of considering multiple spatial scales and taxa to assess the impacts of urbanization on biodiversity.

Antibiotic disturbance affects aquatic microbial community composition and food web interactions but not community resilience.

Notwithstanding the fundamental role that environmental microbes play for ecosystem functioning, data on how microbes react to disturbances are still scarce, and most factors that confer stability to microbial communities are unknown. In this context, antibiotic discharge into the environment is considered a worldwide threat for ecosystems with potential risks to human health. We therefore tested resilience of microbial communities challenged by the presence of an antibiotic. In a continuous culture experiment, we compared the abundance, composition and diversity of microbial communities undisturbed or disturbed by the constant addiction of tetracycline in low (10 µg/L) or intermediate (100 µg/L) concentration (press disturbance). Further, the bacterial communities in the three treatments had to face the sudden pulse disturbance of adding an allochthonous bacterium (Escherichia coli). Tetracycline, even at low concentrations, affected microbial communities by changing their phylogenetic composition and causing cell aggregation. This, however, did not coincide with a reduced microbial diversity, but was mainly caused by a shift in dominance of specific bacterial families. Moreover, the less disturbed community (10 µg/L tetracycline) was sometimes more similar to the control and sometimes more similar to heavily disturbed community (100 µg/L tetracycline). All in all, we could not see a pattern where the communities disturbed with antibiotics were less resilient to a second disturbance introducing E. coli, but they seemed to be able to buffer the input of the allochthonous strain in a similar manner as the control.

Species and hybrids in the genus Diaphanosoma Fischer, 1850 (Crustacea: Branchiopoda: Cladocera).

Cladocerans are well-studied planktonic crustaceans, especially those of the genus Daphnia in which interesting evolutionary questions have been addressed on speciation processes. The aim of the present study is to demonstrate that other genera of cladocerans show similar levels of cryptic diversity, intraspecific gene flow, and thus become useful model systems for comparison. In order to do so, we chose the genus Diaphanosoma, widespread in tropical and temperate areas. We started with a survey of species diversity in the genus Diaphanosoma in Asia using a morphological approach, then obtained sequences from a mitochondrial and a nuclear marker from multiple individuals of different species, performed tests on DNA taxonomy and molecular phylogenies, and assessed the role of hybridization in explaining the cases of mitonuclear discordance. The results are that cryptic diversity occurs in Diaphanosoma, and mitonuclear discordance was found in about 6% of the sequenced animals. Past hybridization is supported as the most likely explanation for the discordance: no evidence was found of first generation hybrids with heterozygous sequences. Our analysis on patterns of genetic diversity in Diaphanosoma supports similarities and differences with what is known in Daphnia.

Lack of host specificity of copepod crustaceans associated with mushroom corals in the Red Sea.

The radiation of symbiotic copepods (Crustacea: Copepoda) living in association with stony corals (Cnidaria: Scleractinia) is considered host-specific and linked to the phylogenetic diversification of their hosts. However, symbiotic copepods are poorly investigated, occurrence records are mostly anecdotal, and no explicit analysis exists regarding their relationship with the hosts. Here, we analysed the occurrence of symbiotic copepods on different co-occurring and phylogenetically closely related scleractinian corals. We used an innovative approach of DNA extraction from single microscopic specimens that preserves the shape of the organisms for integrative morphological studies. The rationale of the study involved: (i) sampling of mushroom corals (Fungiidae) belonging to 13 species and eight genera on different reefs along the Saudi coastline in the Red Sea, (ii) extraction of all the associated copepods, (iii) morphological screening and identification of copepod species, (iv) use of DNA taxonomy on mitochondrial and nuclear markers to determine species boundaries for morphologically unknown copepod species, (v) reconstruction of phylogenies to understand their evolutionary relationships, and (vi) analysis of the ecological drivers of the occurrence, diversity and host specificity of the copepods. The seven species of coral-associated copepods, all new to science, did not show any statistically significant evidence of host-specificity or other pattern of ecological association. We thus suggest that, contrary to most assumptions and previous anecdotal evidence on this coral-copepod host-symbiont system, the association between copepods and their host corals is rather labile, not strict, and not phylogenetically constrained, changing our perception on evolutionary patterns and processes in symbiotic copepods.

Environmental filtering and phylogenetic clustering correlate with the distribution patterns of cryptic protist species.

The community composition of any group of organisms should theoretically be determined by a combination of assembly processes including resource partitioning, competition, environmental filtering, and phylogenetic legacy. Environmental DNA studies have revealed a huge diversity of protists in all environments, raising questions about the ecological significance of such diversity and the degree to which they obey to the same rules as macroscopic organisms. The fast-growing cultivable protist species on which hypotheses are usually experimentally tested represent only a minority of the protist diversity. Addressing these questions for the lesser known majority can only be inferred through observational studies. We conducted an environmental DNA survey of the genus Nebela, a group of closely related testate (shelled) amoeba species, in different habitats within Sphagnum-dominated peatlands. Identification based on the mitochondrial cytochrome c oxidase 1 gene, allowed species-level resolution as well as phylogenetic reconstruction. Community composition varied strongly across habitats and associated environmental gradients. Species showed little overlap in their realized niche, suggesting resource partitioning, and a strong influence of environmental filtering driving community composition. Furthermore, phylogenetic clustering was observed in the most nitrogen-poor samples, supporting phylogenetic inheritance of adaptations in the group of N. guttata. This study showed that the studied free-living unicellular eukaryotes follow to community assembly rules similar to those known to determine plant and animal communities; the same may be true for much of the huge functional and taxonomic diversity of protists.

Genetic spatial structure of an anchialine cave annelid indicates connectivity within - but not between - islands of the Great Bahama Bank.

Land-locked anchialine blue holes are karstic sinkholes and caves with tidally influenced, vertically stratified water bodies that harbor endemic fauna exhibiting variable troglomorphic features. These habitats represent island-like systems, which can serve to elucidate evolutionary and biogeographic processes at local scales. We investigated whether the 'continuous spelean corridor' hypothesis may elucidate the biogeographical distributions of the stygobitic annelid Pelagomacellicephala iliffei (Polynoidae) collected from the Great Bahama and Caicos Banks of the Bahamas Archipelago. Phylogenetic reconstructions were performed using Bayesian Inference on individual and combined datasets of three molecular markers (16S rDNA, COI, 18S rDNA) and species delimitation employed three widely accepted methods in DNA taxonomy, namely GMYC, bPTP, and ABGD. Mantel tests were used to test the effect of geography on genetic structure. Using these analyses, we recovered five independently evolving entities of the focal species across four islands of the Great Bahama Bank including Cat, Eleuthera, Exumas, and Long. Genetic data yielded strong correlations between islands and phylogenetic entities, signifying independent evolutionary histories within anchialine caves across the platform. The island of Eleuthera showed intra-island gene flow and dispersal capabilities between blue holes separated by 115km, providing evidence of a crevicular spelean corridor within the island. However, no evidence of inter-island dispersal is present in the analyzed system. Consistent with previous biogeographic studies of cave crustaceans, the major barriers shaping the cave biota of the Bahamas Archipelago appears to be the deep trenches and channels separating the Bahamian banks.

Integrative Taxonomy Recognizes Evolutionary Units Despite Widespread Mitonuclear Discordance: Evidence from a Rotifer Cryptic Species Complex.

Mitonuclear discordance across taxa is increasingly recognized as posing a major challenge to species delimitation based on DNA sequence data. Integrative taxonomy has been proposed as a promising framework to help address this problem. However, we still lack compelling empirical evidence scrutinizing the efficacy of integrative taxonomy in relation to, for instance, complex introgression scenarios involving many species. Here, we report remarkably widespread mitonuclear discordance between about 15 mitochondrial and 4 nuclear Brachionus calyciflorus groups identified using different species delimitation approaches. Using coalescent-, Bayesian admixture-, and allele sharing-based methods with DNA sequence or microsatellite data, we provide strong evidence in support of hybridization as a driver of the observed discordance. We then describe our combined molecular, morphological, and ecological approaches to resolving phylogenetic conflict and inferring species boundaries. Species delimitations based on the ITS1 and 28S nuclear DNA markers proved a more reliable predictor of morphological variation than delimitations using the mitochondrial COI gene. A short-term competition experiment further revealed systematic differences in the competitive ability between two of the nuclear-delimited species under six different growth conditions, independent of COI delimitations; hybrids were also observed. In light of these findings, we discuss the failure of the COI marker to estimate morphological stasis and morphological plasticity in the B. calyciflorus complex. By using B. calyciflorus as a representative case, we demonstrate the potential of integrative taxonomy to guide species delimitation in the presence of mitonuclear phylogenetic conflicts.

Estimating the magnitude of morphoscapes: how to measure the morphological component of biodiversity in relation to habitats using geometric morphometrics.

Ecological indicators are currently developed to account for the different facets of loss of biological diversity due to direct or indirect effects of human activities. Most ecological indicators include species richness as a metric. Others, such as functional traits and phylogenetic diversity, account for differences in species, even when species richness is the same. Here, we describe and apply a different indicator, called morphoscape dimension, accounting for morphological variability across habitats in a geographical region. We use the case of ground beetles (Coleoptera: Carabidae) in four different habitats in the Po Plain in Northern Italy to exemplify how to quantify the magnitude of the morphological space (i.e. the dimension of the morphoscape) occupied by the species in each habitat using geometric morphometrics. To this aim, we employed a variety of metrics of morphological disparity related to univariate size, and more complex multivariate shape and form. Our 'proof of concept' suggests that metrics assessing size and form might largely tend to simply mirror the information provided by species richness, whereas shape morphoscape disparity may be able to account for non-trivial differences in species traits amongst habitats. This is indicated by the woodland morphoscape being on average bigger than that of crops, the most species-rich habitat, despite having almost 20% less species. We conclude suggesting that the analysis of morphoscape dimension has the potential to become a new additional and complimentary tool in the hands of conservation biologists and ecologists to explore and quantify habitat complexity and inform decisions on management and conservation based on a wide set of ecological indicators.