Lara Konecny-Dupré

Timetree of Aselloidea Reveals Species Diversification Dynamics in Groundwater

A key challenge for biologists is to document and explain global patterns of diversification in a wide range of environments. Here, we explore patterns of continental-scale diversification in a groundwater species-rich clade, the superfamily Aselloidea (Pancrustacea: Isopoda). Our analyses supported a constant diversification rate during most of the course of Aselloidea evolution, until 4-15 Ma when diversification rates started to decrease. This constant accumulation of lineages challenges the view that groundwater species diversification in temperate regions might have been primarily driven by major changes in physical environment leading to the extinction of surface populations and subsequent synchronous isolation of multiple groundwater populations. Rather than acting synchronously over broad geographic regions, factors causing extinction of surface populations and subsequent reproductive isolation of groundwater populations may act in a local and asynchronous manner, thereby resulting in a constant speciation rate over time. Our phylogeny also revealed several cases of parapatric distributions among closely related surface-water and groundwater species suggesting that species diversification could also arise from a process of disruptive selection along the surface-subterranean environmental gradient. Our results call for re-evaluating the spatial scale and timing of factors causing diversification events in groundwater.

Less effective selection leads to larger genomes

The evolutionary origin of the striking genome size variations found in eukaryotes remains enigmatic. The effective size of populations, by controlling selection efficacy, is expected to be a key parameter underlying genome size evolution. However, this hypothesis has proved difficult to investigate using empirical data sets. Here, we tested this hypothesis using 22 de novo transcriptomes and low-coverage genomes of asellid isopods, which represent 11 independent habitat shifts from surface water to resource-poor groundwater. We show that these habitat shifts are associated with higher transcriptome-wide [Formula: see text] After ruling out the role of positive selection and pseudogenization, we show that these transcriptome-wide [Formula: see text] increases are the consequence of a reduction in selection efficacy imposed by the smaller effective population size of subterranean species. This reduction is paralleled by an important increase in genome size (25% increase on average), an increase also confirmed in subterranean decapods and mollusks. We also control for an adaptive impact of genome size on life history traits but find no correlation between body size, or growth rate, and genome size. We show instead that the independent increases in genome size measured in subterranean isopods are the direct consequence of increasing invasion rates by repeat elements, which are less efficiently purged out by purifying selection. Contrary to selection efficacy, polymorphism is not correlated to genome size. We propose that recent demographic fluctuations and the difficulty of observing polymorphism variation in polymorphism-poor species can obfuscate the link between effective population size and genome size when polymorphism data are used alone.

Bayesian phylogeographic inferences reveal contrasting colonization dynamics among European groundwater isopods.

The potentially important role of northern microrefugia during postglacial dispersal is challenging the view of southern Europe as a refuge and source area of European biota. In groundwaters, large geographic ranges of presumably good dispersers are increasingly suspected to consist of assemblages of cryptic species with narrow ranges. Moreover, a large species range, even when confirmed by molecular evidence, tells us little about the spatiotemporal dynamics of dispersal. Here, we used phylogenetic inferences, species delineation methods and Bayesian phylogeographic diffusion models to test for the likelihood of postglacial colonization from distant refugia among five morphospecies of Proasellus (Isopoda, Asellidae). All morphospecies except one were monophyletic, but they comprised a total of 15–17 cryptic species. Three cryptic species retained ranges that spanned a distance >650 km, similar to that of the nominal morphospecies. Bayesian diffusion models based on mitochondrial markers revealed considerable spatiotemporal heterogeneity in dispersal rates, suggesting that short‐time dispersal windows were instrumental in shaping species ranges. Only one species was found to experience a recent, presumably postglacial, range expansion. The Jura and Alpine foothills probably played a major role in maintaining diversity within Proasellus in northern regions by acting both as diversification hotspots and Pleistocene refugia. Gaining insight into the spatiotemporal heterogeneity of dispersal rates revealed contrasting colonization dynamics among species that were not consistent with a global postglacial colonization of Europe from distant refugia.

Do cryptic species matter in macroecology? Sequencing European groundwater crustaceans yields smaller ranges but does not challenge biodiversity determinants

Ecologists increasingly rely on molecular delimitation methods (MMs) to identify species boundaries, thereby potentially increasing the number of putative species because of the presence of morphologically cryptic species. It has been argued that cryptic species could challenge our understanding of what determine large-scale biodiversity patterns which have traditionally been documented from morphology alone. Here, we used morphology and three MMs to derive four different sets of putative species among the European groundwater crustaceans. Then, we used regression models to compare the relative importance of spatial heterogeneity, productivity and historical climates, in shaping species richness and range size patterns across sets of putative species. We tested three predictions. First, MMs would yield many more putative species than morphology because groundwater is a constraining environment allowing little morphological changes. Second, for species richness, MMs would increase the importance of spatial heterogeneity because cryptic species are more likely along physical barriers separating ecologically similar regions than along resource gradients promoting ecologically-based divergent selection. Third, for range size, MMs would increase the importance of historical climates because of reduced and asymmetrical fragmentation of large morphological species ranges at northern latitudes. MMs yielded twice more putative species than morphology and decreased by 10-fold the average species range size. Yet, MMs strengthened the mid-latitude ridge of high species richness and the Rapoport effect of increasing range size at higher latitudes. Species richness predictors did not vary between morphology and MMs but the latter increased the proportion of variance in range size explained by historical climates. These findings demonstrate that our knowledge of groundwater biodiversity determinants is robust to overlooked cryptic species because the latter are homogeneously distributed along environmental gradients. Yet, our findings call for incorporating multiple species delimitation methods into the analysis of large-scale biodiversity patterns across a range of taxa and ecosystems.

No Evidence That Nitrogen Limitation Influences the Elemental Composition of Isopod Transcriptomes and Proteomes

The field of stoichiogenomics aims at understanding the influence of nutrient limitations on the elemental composition of the genome, transcriptome, and proteome. The 20 amino acids and the 4 nt differ in the number of nutrients they contain, such as nitrogen (N). Thus, N limitation shall theoretically select for changes in the composition of proteins or RNAs through preferential use of N-poor amino acids or nucleotides, which will decrease the N-budget of an organism. While these N-saving mechanisms have been evidenced in microorganisms, they remain controversial in multicellular eukaryotes. In this study, we used 13 surface and subterranean isopod species pairs that face strongly contrasted N limitations, either in terms of quantity or quality. We combined in situ nutrient quantification and transcriptome sequencing to test if N limitation selected for N-savings through changes in the expression and composition of the transcriptome and proteome. No evidence of N-savings was found in the total N-budget of transcriptomes or proteomes or in the average protein N-cost. Nevertheless, subterranean species evolving in N-depleted habitats displayed lower N-usage at their third codon positions. To test if this convergent compositional change was driven by natural selection, we developed a method to detect the strand-asymmetric signature that stoichiogenomic selection should leave in the substitution pattern. No such signature was evidenced, indicating that the observed stoichiogenomic-like patterns were attributable to nonadaptive processes. The absence of stoichiogenomic signal despite strong N limitation within a powerful phylogenetic framework casts doubt on the existence of stoichiogenomic mechanisms in metazoans.

Microsatellite Development and First Population Size Estimates for the Groundwater Isopod Proasellus walteri

Effective population size (N e) is one of the most important parameters in, ecology, evolutionary and conservation biology; however, few studies of N e in surface freshwater organisms have been published to date. Even fewer studies have been carried out in groundwater organisms, although their evolution has long been considered to be particularly constrained by small N e. In this study, we estimated the contemporary effective population size of the obligate groundwater isopod: Proaselluswalteri (Chappuis, 1948). To this end, a genomic library was enriched for microsatellite motifs and sequenced using 454 GS-FLX technology. A total of 54,593 reads were assembled in 10,346 contigs or singlets, of which 245 contained candidate microsatellite sequences with suitable priming sites. Ninety-six loci were tested for amplification, polymorphism and multiplexing properties, of which seven were finally selected for N e estimation. Linkage disequilibrium and approximate Bayesian computation methods revealed that N e in this small interstitial groundwater isopod could reach large sizes (> 585 individuals). Our results suggest that environmental conditions in groundwater, while often referred to as extreme, are not necessarily associated with small N e.

Life History Traits Impact the Nuclear Rate of Substitution but Not the Mitochondrial Rate in Isopods

The rate of molecular evolution varies widely among species. Life history traits (LHTs) have been proposed as a major driver of these variations. However, the relative contribution of each trait is poorly understood. Here, we test the influence of metabolic rate (MR), longevity, and generation time (GT) on the nuclear and mitochondrial synonymous substitution rates using a group of isopod species that have made multiple independent transitions to subterranean environments. Subterranean species have repeatedly evolved a lower MR, a longer lifespan and a longer GT. We assembled the nuclear transcriptomes and the mitochondrial genomes of 13 pairs of closely related isopods, each pair composed of one surface and one subterranean species. We found that subterranean species have a lower rate of nuclear synonymous substitution than surface species whereas the mitochondrial rate remained unchanged. We propose that this decoupling between nuclear and mitochondrial rates comes from different DNA replication processes in these two compartments. In isopods, the nuclear rate is probably tightly controlled by GT alone. In contrast, mitochondrial genomes appear to replicate and mutate at a rate independent of LHTs. These results are incongruent with previous studies, which were mostly devoted to vertebrates. We suggest that this incongruence can be explained by developmental differences between animal clades, with a quiescent period during female gametogenesis in mammals and birds which imposes a nuclear and mitochondrial rate coupling, as opposed to the continuous gametogenesis observed in most arthropods.

Integrating multiple species criteria and species hypotheses in subterranean biology

Decades of debates around the species problem have resulted in the emergence of a unified species concept with multiple criteria to delimit species taxa. Many biologists now agree to consider species as separately evolving segments of metapopulation lineages (i.e. the species concept), and to consider species taxa (i.e. the elementary units used in biodiversity science) as scientific hypotheses of separately evolving entities. In this framework, sets of species hypotheses are generated using different criteria (i.e. morphological distinguishability, genetic isolation) that mirror the properties expressed by species at different times and sequential orders during the extended and heterogeneous process of speciation. This conceptual and methodological advance in taxonomy has several implications for biodiversity science. First, species taxa represent a heterogeneous set of hypotheses whose properties are contingent on the heterogeneous, continuous and extended nature of speciation. Second, species databases need to integrate information on the diverse properties of species by attributing specimens to multiple species hypotheses generated using different delimitation criteria. Third, biodiversity science at large can provide novel insights into biodiversity processes by incorporating multiple species hypotheses into the analysis of biodiversity patterns. Here, we show how these implications ‡,‡ § | ¶ # ‡ ‡ # ¶ ‡