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.

Anthropization level of Lascaux Cave microbiome shown by regional-scale comparisons of pristine and anthropized caves.

Limestone areas across the world develop karstic caves, which are populated by a wide range of macro- and microorganisms. Many of these caves display Paleolithic art or outstanding speleothems, and in the last century they have been subjected to anthropization due to touristic management and intense human frequentation. Despite their cultural importance and associated conservation issues, the impact of anthropization on cave biodiversity is not known. Here, we show that anthropization is associated with specific cave biota modifications. We compared diversity in four pristine caves, four anthropized show caves, and the iconic Lascaux Cave with even stronger anthropization. The predominant microbial higher taxa were the same in all caves, but the most anthropized cave (Lascaux) was unique as it differed from the eight others by a higher proportion of Bacteroidetes bacteria and the absence of Euryarchaeota and Woesearchaeota archaea. Anthropization resulted in lower diversity and altered community structure for bacteria and archaea on cave walls, especially in Lascaux, but with a more limited effect on microeukaryotes and arthropods. Our findings fill a key gap in our understanding of the response of karstic communities to anthropization, by revealing that tourism-related anthropization impacts on the prokaryotic microbiome rather than on eukaryotic residents, and that it shapes cave biota irrespective of cave natural features.

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.

Nomenclature for the Nameless: A Proposal for an Integrative Molecular Taxonomy of Cryptic Diversity Exemplified by Planktonic Foraminifera.

Investigations of biodiversity, biogeography, and ecological processes rely on the identification of "species" as biologically significant, natural units of evolution. In this context, morphotaxonomy only provides an adequate level of resolution if reproductive isolation matches morphological divergence. In many groups of organisms, morphologically defined species often disguise considerable genetic diversity, which may be indicative of the existence of cryptic species. The diversity hidden by morphological species can be disentangled through genetic surveys, which also provide access to data on the ecological distribution of genetically circumscribed units. These units can be identified by unique DNA sequence motifs and allow studies of evolutionary and ecological processes at different levels of divergence. However, the nomenclature of genetically circumscribed units within morphological species is not regulated and lacks stability. This represents a major obstacle to efforts to synthesize and communicate data on genetic diversity for multiple stakeholders. We have been confronted with such an obstacle in our work on planktonic foraminifera, where the stakeholder community is particularly diverse, involving geochemists, paleoceanographers, paleontologists, and biologists, and the lack of stable nomenclature beyond the level of formal morphospecies prevents effective transfer of knowledge. To circumvent this problem, we have designed a stable, reproducible, and flexible nomenclature system for genetically circumscribed units, analogous to the principles of a formal nomenclature system. Our system is based on the definition of unique DNA sequence motifs collocated within an individual, their typification (in analogy with holotypes), utilization of their hierarchical phylogenetic structure to define levels of divergence below that of the morphospecies, and a set of nomenclature rules assuring stability. The resulting molecular operational taxonomic units remain outside the domain of current nomenclature codes, but are linked to formal morphospecies as regulated by the codes. Subsequently, we show how this system can be applied to classify genetically defined units using the SSU rDNA marker in planktonic foraminifera and we highlight its potential use for other groups of organisms where similarly high levels of connectivity between molecular and formal taxonomies can be achieved.

Mitochondrial genomes reveal the extinct Hippidion as an outgroup to all living equids.

Hippidions were equids with very distinctive anatomical features. They lived in South America 2.5 million years ago (Ma) until their extinction approximately 10 000 years ago. The evolutionary origin of the three known Hippidion morphospecies is still disputed. Based on palaeontological data, Hippidion could have diverged from the lineage leading to modern equids before 10 Ma. In contrast, a much later divergence date, with Hippidion nesting within modern equids, was indicated by partial ancient mitochondrial DNA sequences. Here, we characterized eight Hippidion complete mitochondrial genomes at 3.4-386.3-fold coverage using target-enrichment capture and next-generation sequencing. Our dataset reveals that the two morphospecies sequenced (H. saldiasi and H. principale) formed a monophyletic clade, basal to extant and extinct Equus lineages. This contrasts with previous genetic analyses and supports Hippidion as a distinct genus, in agreement with palaeontological models. We date the Hippidion split from Equus at 5.6-6.5 Ma, suggesting an early divergence in North America prior to the colonization of South America, after the formation of the Panamanian Isthmus 3.5 Ma and the Great American Biotic Interchange.

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.

The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton.

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

A collaborative backbone resource for comparative studies of subterranean evolution: The World Asellidae database.

Transition to novel environments, such as groundwater colonization by surface organisms, provides an excellent research ground to study phenotypic evolution. However, interspecific comparative studies on evolution to groundwater life are few because of the challenge in assembling large ecological and molecular resources for species-rich taxa comprised of surface and subterranean species. Here, we make available to the scientific community an operational set of working tools and resources for the Asellidae, a family of freshwater isopods containing hundreds of surface and subterranean species. First, we release the World Asellidae database (WAD) and its web application, a sustainable and FAIR solution to producing and sharing data and biological material. WAD provides access to thousands of species occurrences, specimens, DNA extracts and DNA sequences with rich metadata ensuring full scientific traceability. Second, we perform a large-scale dated phylogenetic reconstruction of Asellidae to support phylogenetic comparative analyses. Of 424 terminal branches, we identify 34 pairs of surface and subterranean species representing independent replicates of the transition from surface water to groundwater. Third, we exemplify the usefulness of WAD for documenting phenotypic shifts associated with colonization of subterranean habitats. We provide the first phylogenetically controlled evidence that body size of males decreases relative to that of females upon groundwater colonization, suggesting competition for rare receptive females selects for smaller, more agile males in groundwater. By making these tools and resources widely accessible, we open up new opportunities for exploring how phenotypic traits evolve in response to changes in selective pressures and trade-offs during groundwater colonization.

Microbiome analysis in Lascaux Cave in relation to black stain alterations of rock surfaces and collembola.

Anthropization of Palaeolithic caves open for tourism may favour collembola invasion and result in the formation of black stains attributed to pigmented fungi. However, ecological processes underpinning black stain formation are not fully understood. Here, we tested the hypotheses that black stains from the Apse room of Lascaux Cave display a specific microbiota enriched in pigmented fungi, and that collembola thriving on the stains have the potential to consume and disseminate these black fungi. Metabarcoding showed that the microbiota of black stains and neighbouring unstained parts strongly differed, with in black stains a higher prevalence of Ochroconis and other pigmented fungi and the strong regression of Pseudomonas bacteria (whose isolates inhibited in vitro the growth of pigmented fungi). Isotopic analyses indicated that Folsomia candida collembola thriving on stains could feed on black stain in situ and assimilate the pigmented fungi they were fed with in vitro. They could carry these fungi and disseminate them when tested with complex black stains from Lascaux. This shows that black stain formation is linked to the development of pigmented fungi, which coincides with the elimination of antagonistic pseudomonads, and points towards a key role of F. candida collembola in the dynamics of pigmented fungi.

Renewal of planktonic foraminifera diversity after the Cretaceous Paleogene mass extinction by benthic colonizers.

The biotic crisis following the end-Cretaceous asteroid impact resulted in a dramatic renewal of pelagic biodiversity. Considering the severe and immediate effect of the asteroid impact on the pelagic environment, it is remarkable that some of the most affected pelagic groups, like the planktonic foraminifera, survived at all. Here we queried a surface ocean metabarcoding dataset to show that calcareous benthic foraminifera of the clade Globothalamea are able to disperse actively in the plankton, and we show using molecular clock phylogeny that the modern planktonic clades originated from different benthic ancestors that colonized the plankton after the end-Cretaceous crisis. We conclude that the diversity of planktonic foraminifera has been the result of a constant leakage of benthic foraminifera diversity into the plankton, continuously refueling the planktonic niche, and challenge the classical interpretation of the fossil record that suggests that Mesozoic planktonic foraminifera gave rise to the modern communities.

Plant resistance to mechanical stress: evidence of an avoidance-tolerance trade-off.

External mechanical forces resulting from the pressure exerted by wind or water movement are a major stress factor for plants and may cause regular disturbances in many ecosystems. A plant's ability to resist these forces relies either on minimizing the forces encountered by the plant (avoidance strategy), or on maximizing its resistance to breakage (tolerance strategy). We investigated plant resistance strategies using aquatic vegetation as a model, and examined whether avoidance and tolerance are negatively correlated. We tested the avoidance-tolerance correlation across 28 species using a phylogenetically corrected analysis, after construction of a molecular phylogeny for the species considered. Different species demonstrated contrasting avoidance and tolerance and we demonstrated a significant negative relationship between the two strategies, which suggests an avoidance-tolerance trade-off. Negative relationships may result from costs that each strategy incurs or from constraints imposed by physical laws on plant tissues. The existence of such a trade-off has important ecological and evolutionary consequences. It would lead to constraints on the evolution and variation of both strategies, possibly limiting their evolution and may constrain many morphological, anatomical and architectural traits that underlie avoidance and tolerance.

Scent evolution in Chinese roses.

The phenolic methyl ether 3,5-dimethoxytoluene (DMT) is a major scent compound of many modern rose varieties, and its fragrance participates in the characteristic "tea scent" that gave their name to Tea and Hybrid Tea roses. Among wild roses, phenolic methyl ether (PME) biosynthesis is restricted to Chinese rose species, but the progenitors of modern roses included both European and Chinese species (e.g., Rosa chinensis cv Old Blush), so this trait was transmitted to their hybrid progeny. The last steps of the biosynthetic pathways leading to DMT involve two methylation reactions catalyzed by the highly similar orcinol O-methyltransferases (OOMT) 1 and 2. OOMT1 and OOMT2 enzymes exhibit different substrate specificities that are consistent with their operating sequentially in DMT biosynthesis. Here, we show that these different substrate specificities are mostly due to a single amino acid polymorphism in the phenolic substrate binding site of OOMTs. An analysis of the OOMT gene family in 18 species representing the diversity of the genus Rosa indicated that only Chinese roses possess both the OOMT2 and the OOMT1 genes. In addition, we provide evidence that the Chinese-rose-specific OOMT1 genes most probably evolved from an OOMT2-like gene that has homologues in the genomes of all extant roses. We propose that the emergence of the OOMT1 gene may have been a critical step in the evolution of scent production in Chinese roses.

Bedrock radioactivity influences the rate and spectrum of mutation.

All organisms on Earth are exposed to low doses of natural radioactivity but some habitats are more radioactive than others. Yet, documenting the influence of natural radioactivity on the evolution of biodiversity is challenging. Here, we addressed whether organisms living in naturally more radioactive habitats accumulate more mutations across generations using 14 species of waterlice living in subterranean habitats with contrasted levels of radioactivity. We found that the mitochondrial and nuclear mutation rates across a waterlouse species' genome increased on average by 60% and 30%, respectively, when radioactivity increased by a factor of three. We also found a positive correlation between the level of radioactivity and the probability of G to T (and complementary C to A) mutations, a hallmark of oxidative stress. We conclude that even low doses of natural bedrock radioactivity influence the mutation rate possibly through the accumulation of oxidative damage, in particular in the mitochondrial genome.

Freeze tolerance evolution among anurans: Frequency and timing of appearance.

Despite numerous mechanistic studies on physiological responses supporting freeze tolerance in anurans, few have addressed the evolutionary significance of this trait. We thus investigated the phylogenetic relationships among anuran species whose freeze tolerance has been assessed and in combination with new data on freezing tolerance of two closely related species of the European brown frogs (Rana temporaria and Rana dalmatina). The species we studied exhibited short survival times in frozen state (around 8 h for both species). Phylogenetic analysis suggests that freeze tolerance evolved at least two times among Ranidae and one or two times among Hylidae and never in Bufonidae. Furthermore, in order to assess the timing of divergence of this character we used a relaxed molecular clock created, and found that the most recent separation between a freeze tolerant species and a freeze intolerant species dates from 15.9 +/- 7.6 Myr (Rana arvalis and R. temporaria). The comparison between these two species thus represents the best current model to understand freeze tolerance evolution. Addressing the evolution of this trait with such large-scale approaches will not only improve our understanding of cold hardiness strategies, but might also create a framework guiding future comparative studies.

DNA from extinct giant lemurs links archaeolemurids to extant indriids.

BACKGROUND: Although today 15% of living primates are endemic to Madagascar, their diversity was even greater in the recent past since dozens of extinct species have been recovered from Holocene excavation sites. Among them were the so-called "giant lemurs" some of which weighed up to 160 kg. Although extensively studied, the phylogenetic relationships between extinct and extant lemurs are still difficult to decipher, mainly due to morphological specializations that reflect ecology more than phylogeny, resulting in rampant homoplasy. RESULTS: Ancient DNA recovered from subfossils recently supported a sister relationship between giant "sloth" lemurs and extant indriids and helped to revise the phylogenetic position of Megaladapis edwardsi among lemuriformes, but several taxa - such as the Archaeolemuridae - still await analysis. We therefore used ancient DNA technology to address the phylogenetic status of the two archaeolemurid genera (Archaeolemur and Hadropithecus). Despite poor DNA preservation conditions in subtropical environments, we managed to recover 94- to 539-bp sequences for two mitochondrial genes among 5 subfossil samples. CONCLUSION: This new sequence information provides evidence for the proximity of Archaeolemur and Hadropithecus to extant indriids, in agreement with earlier assessments of their taxonomic status (Primates, Indrioidea) and in contrast to recent suggestions of a closer relationship to the Lemuridae made on the basis of analyses of dental developmental and postcranial characters. These data provide new insights into the evolution of the locomotor apparatus among lemurids and indriids.

Horizontal gene transfer and phylogenetics.

The initial analysis of complete genomes has suggested that horizontal gene transfer events are very frequent between microorganisms. This could potentially render the inference, and even the concept itself, of the organismal phylogeny impossible. However, a coherent phylogenetic pattern has recently emerged from an analysis of about a hundred genes, the so-called 'core', strongly suggesting that it is possible to infer the phylogeny of prokaryotes. Also, estimation of the frequency of horizontal gene transfers at the genome level in a phylogenetic context seems to indicate that it is rather low, although of significant biological impact. Nevertheless, it should be emphasized that the history of microorganisms cannot be properly represented by the phylogeny of the core, which represents only a tiny fraction of the genome. This history, even if horizontal gene transfers are rare, should be represented by a network surrounding the core phylogeny.

PFR²: a curated database of planktonic foraminifera 18S ribosomal DNA as a resource for studies of plankton ecology, biogeography and evolution.

Planktonic foraminifera (Rhizaria) are ubiquitous marine pelagic protists producing calcareous shells with conspicuous morphology. They play an important role in the marine carbon cycle, and their exceptional fossil record serves as the basis for biochronostratigraphy and past climate reconstructions. A major worldwide sampling effort over the last two decades has resulted in the establishment of multiple large collections of cryopreserved individual planktonic foraminifera samples. Thousands of 18S rDNA partial sequences have been generated, representing all major known morphological taxa across their worldwide oceanic range. This comprehensive data coverage provides an opportunity to assess patterns of molecular ecology and evolution in a holistic way for an entire group of planktonic protists. We combined all available published and unpublished genetic data to build PFR(2), the Planktonic foraminifera Ribosomal Reference database. The first version of the database includes 3322 reference 18S rDNA sequences belonging to 32 of the 47 known morphospecies of extant planktonic foraminifera, collected from 460 oceanic stations. All sequences have been rigorously taxonomically curated using a six-rank annotation system fully resolved to the morphological species level and linked to a series of metadata. The PFR(2) website, available at http://pfr2.sb-roscoff.fr, allows downloading the entire database or specific sections, as well as the identification of new planktonic foraminiferal sequences. Its novel, fully documented curation process integrates advances in morphological and molecular taxonomy. It allows for an increase in its taxonomic resolution and assures that integrity is maintained by including a complete contingency tracking of annotations and assuring that the annotations remain internally consistent.

SSU rDNA divergence in planktonic foraminifera: molecular taxonomy and biogeographic implications.

The use of planktonic foraminifera in paleoceanography requires taxonomic consistency and precise assessment of the species biogeography. Yet, ribosomal small subunit (SSUr) DNA analyses have revealed that most of the modern morpho-species of planktonic foraminifera are composed of a complex of several distinct genetic types that may correspond to cryptic or pseudo-cryptic species. These genetic types are usually delimitated using partial sequences located at the 3'end of the SSUrDNA, but typically based on empirical delimitation. Here, we first use patristic genetic distances calculated within and among genetic types of the most common morpho-species to show that intra-type and inter-type genetic distances within morpho-species may significantly overlap, suggesting that genetic types have been sometimes inconsistently defined. We further apply two quantitative and independent methods, ABGD (Automatic Barcode Gap Detection) and GMYC (General Mixed Yule Coalescent) to a dataset of published and newly obtained partial SSU rDNA for a more objective assessment of the species status of these genetic types. Results of these complementary approaches are highly congruent and lead to a molecular taxonomy that ranks 49 genetic types of planktonic foraminifera as genuine (pseudo)cryptic species. Our results advocate for a standardized sequencing procedure allowing homogenous delimitations of (pseudo)cryptic species. On the ground of this revised taxonomic framework, we finally provide an integrative taxonomy synthesizing geographic, ecological and morphological differentiations that can occur among the genuine (pseudo)cryptic species. Due to molecular, environmental or morphological data scarcities, many aspects of our proposed integrative taxonomy are not yet fully resolved. On the other hand, our study opens up the potential for a correct interpretation of environmental sequence datasets.

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.