In the deep Pacific and Indian oceans, Alvinellid worms diversified about 100 million years ago to colonize a variety of hydrothermal vent environments. It has been suggested that the last common ancestor of this family was a thermophilic species. However, the evolutionary history of these worms is complex, with putative gene flow among ancestors. In this study, we investigated the early evolution of the family in relation to the diversification of thermal niches. We demonstrated that phylogenetic histories of alvinellid species greatly vary along chromosomes, possibly due to allele introgression between nascent ecotypes. We performed sequence reconstructions of ancestral cytosolic MDH and Cu/Zn SOD under the two most frequent phylogenies encountered along the genome. In silico simulations of folding stability for these two enzymes were highly correlated with their biophysical characterization (micro-calorimetry and differential scanning fluorimetry), and were used to predict the thermostability of thousands of alternatively reconstructed proteins. We further generalized the reconstruction at the proteome scale, taking the amino-acid usage bias as a proxy for folding stability in a new phylogenetic model accounting for amino-acid variations over time. Both approaches agree that the last common ancestors of Alvinellidae gained highly stable proteins, comparable to proteins of present-day thermotolerant species. Our computational protocol allowed us to validate this scenario under multiple phylogenetic hypotheses and different sequence reconstruction models. Considered together, Alvinellidae shed light on how metazoan species diversify in order to colonize different thermal niches, combining adaptive mutation and selection on genetic variants.

How the interplay of biotic and abiotic factors shapes current genetic diversity at the community level remains an open question, particularly in the deep sea. Comparative phylogeography of multiple species can reveal the influence of past climatic events, geographic barriers, and species life history traits on spatial patterns of genetic structure across lineages. To shed light on the factors that shape community-level genetic variation and to improve our understanding of deep-sea biogeographic patterns, we conducted a comparative population genomics study on seven hydrothermal vent species co-distributed in the Back-Arc Basins of the Southwest Pacific region. Using ddRAD-seq, we compared the range-wide distribution of genomic diversity across species and discovered a shared phylogeographic break. Demogenetic inference revealed shared histories of lineage divergence and a secondary contact. Low levels of asymmetric gene flow probably occurred in most species between the Woodlark and North Fiji basins, but the exact location of contact zones varied from species to species. For two species, we found individuals from the two lineages co-occurring in sympatry in Woodlark Basin. Although species exhibit congruent patterns of spatial structure (Eastern vs. Western sites), they also show variation in the degree of divergence among lineages across the suture zone. Our results also show heterogeneous gene flow across the genome, indicating possible partial reproductive isolation between lineages and early speciation. Our comparative study highlights the pivotal role of historical and contemporary factors, underscoring the need for a comprehensive approach-especially in addressing knowledge gaps on the life history traits of deep-sea species.

The Alvinellidae are a family of worms that are endemic to deep-sea hydrothermal vents in the Pacific and Indian Oceans. These annelid worms, a sister group to the Ampharetidae, occupy a wide range of thermal habitats. The family includes the most thermotolerant marine animals described to date such as the Pompeii worm Alvinella pompejana, and other species living at much lower temperatures such as Paralvinella grasslei or Paralvinella pandorae. The phylogeny of this family has not been studied extensively. It is, however, a complex case where molecular phylogenies have given conflicting results, especially concerning the monophyletic or polyphyletic nature of the genus Paralvinella.We carried out a comprehensive study of the phylogeny of this family using the best molecular data currently available from RNAseq datasets, leading to several hundred orthologous transcripts for 11 of the 14 species currently described or in description. The results obtained by the two most popular phylogenetic inference models (using either gene concatenation with maximum likelihood, or a coalescent-based model from gene trees) were compared using a series of ampharetid and terebellid outgroups.Our study shows that the global phylogenetic signal favors the hypothesis of paraphyly for the Paralvinella genus, with P. pandorae being sister species of the other Alvinellidae.However, a high number of gene trees also supports the hypothesis of alternative trees in which the monophyly of the Paralvinella genus, as initially proposed by Desbruyères and Laubier, is valid with the species P. pandorae and Paralvinella unidentata being sister species. According to molecular dating, the radiation of the Alvinellidae was rapid and took place in a short period of time between 80 and 110 million years ago. This is reflected at the genomic scale by gene trees equally separated between different phylogenetic hypothesis, showing high rates of incomplete lineage sorting between the first lineages of the Alvinellidae and probable gene transfers. Although different genomic regions seem to have different phylogenetic stories in the early step of the alvinellid radiation, our study supports the view that the two P. pandorae species can be grouped into a separate genus (possibly Nautalvinella) and that the Miralvinella subgenus, defined by Desbruyères and Laubier, is not valid anymore.

Microsatellite data (12 loci) of 14 populations of Laminaria digitata covering the entire species range

Understanding drivers of biodiversity patterns is essential to evaluate the potential impact of deep-sea mining on ecosystems resilience. While the South West Pacific forms an independent biogeographic province for hydrothermal vent fauna, different degrees of connectivity among basins were previously reported for a variety of species depending on their ability to disperse. In this study, we compared phylogeographic patterns of several vent gastropods across South West Pacific back-arc basins and the newly-discovered La Scala site on the Woodlark Ridge by analysing their genetic divergence using a barcoding approach. We focused on six genera of vent gastropods widely distributed in the region: Lepetodrilus , Symmetromphalus , Lamellomphalus , Shinkailepas , Desbruyeresia and Provanna . A wide-range sampling was conducted at different vent fields across the Futuna Volcanic Arc, the Manus, Woodlark, North Fiji, and Lau Basins, during the CHUBACARC cruise in 2019. The Cox1 -based genetic structure of geographic populations was examined for each taxon to delineate putative cryptic species and assess potential barriers or contact zones between basins. Results showed contrasted phylogeographic patterns among species, even between closely related species. While some species are widely distributed across basins (i.e. Shinkailepas tollmanni , Desbruyeresia melanioides and Lamellomphalus ) without evidence of strong barriers to gene flow, others are restricted to one (i.e. Shinkailepas tufari complex of cryptic species, Desbruyeresia cancellata and D . costata ). Other species showed intermediate patterns of isolation with different lineages separating the Manus Basin from the Lau/North Fiji Basins (i.e. Lepetodrilus schrolli , Provanna and Symmetromphalus spp.). Individuals from the Woodlark Basin were either endemic to this area (though possibly representing intermediate OTUs between the Manus Basin and the other eastern basins populations) or, coming into contact from these basins, highlighting the stepping-stone role of the Woodlark Basin in the dispersal of the South West Pacific vent fauna. Results are discussed according to the dispersal ability of species and the geological history of the South West Pacific.

Deep hydrothermal vents are highly fragmented and unstable habitats at all temporal and spatial scales. Such environmental dynamics likely play a non-negligible role in speciation. Little is, however, known about the evolutionary processes that drive population-level differentiation and vent species isolation and, more specifically, how geography and habitat specialisation interplay in the species history of divergence. In this study, the species range and divergence of Alviniconcha snails that occupy active Western Pacific vent fields was assessed by using sequence variation data of the mitochondrial Cox1 gene, RNAseq, and ddRAD-seq. Combining morphological description and sequence datasets of the three species across five basins, we confirmed that A. kojimai, A. boucheti, and A. strummeri, while partially overlapping over their range, display high levels of divergence in the three genomic compartments analysed that usually encompass values retrieved for reproductively isolated species with divergences rang from 9% to 12.5% (mtDNA) and from 2% to 3.1% (nuDNA). Moreover, the three species can be distinguished on the basis of their external morphology by observing the distribution of bristles and the shape of the columella. According to this sampling, A. boucheti and A. kojimai form an east-to-west species abundance gradient, whereas A. strummeri is restricted to the Futuna Arc/Lau and North Fiji Basins. Surprisingly, population models with both gene flow and population size heterogeneities among genomes indicated that these three species are still able to exchange genes due to secondary contacts at some localities after a long period of isolation.

The stalked barnacle Pollicipes pollicipes inhabits rocky shores from the Atlantic coasts of Brittany (France) to Senegal. Because of the culinary traditions of southern Europe, stalked barnacles represent an important target species for local fisheries on the Iberian Peninsula. To manage this fishery sustainably, it is therefore important to assess the dynamics of local populations over the Iberian coast, and how they are interconnected at a wider scale using finely tuned genetic markers. In this work, a new enriched library of GT microsatellites for P. pollicipes was prepared and sequenced using Ion Torrent™ Next Gen-Sequencing Technology. 1,423 adults and juveniles were sampled in 15 localities of three geographic regions: southern Portugal, Galicia and Asturias (both in northern Spain). Twenty polymorphic loci arranged in five multiplex PCRs were then tested and validated as new molecular tools to address the spatial and temporal genetic patterns of P. pollicipes . Our results revealed high genetic diversity among adults. However, juveniles were genetically more structured than their adult counterparts, which alternatively displayed much more connectivity among the three studied regions. The lack of spatial genetic heterogeneity in adults may be due to the overlapping of several generations of settlers coming from different geographic origins, which mainly depends on the orientation of residual currents along the coast during reproduction. The genetic differentiation of juveniles may indeed be congruent with Iberian Peninsula hydrodynamics, which can produce chaotic genetic patchiness (CGP) at small temporal scales due to sweepstake reproductive success, collective dispersal and/or self-recruitment. Remarkably, most of the genetic heterogeneity of juveniles found in this work was located in Galicia, which could represent an admixture between distinct metapopulations or an old refuge for the most northern populations. To conclude, high genetic variation in P. pollicipes can lead to the false impression of population panmixia at the Iberian scale by masking more restricted and current-driven larval exchanges between regions. This possibility should be taken into consideration for further specific management and conservation plans for the species over the Iberian Peninsula.

Reconstructing ancestral characters on a phylogeny is an arduous task because the observed states at the tips of the tree correspond to a single realization of the underlying evolutionary process. Recently, it was proposed that ancestral traits can be indirectly estimated with the help of molecular data, based on the fact that life history traits influence substitution rates. Here we challenge these new approaches in the Cetartiodactyla, a clade of large mammals which, according to paleontology, derive from small ancestors. Analysing transcriptome data in 41 species, of which 22 were newly sequenced, we provide a dated phylogeny of the Cetartiodactyla and report a significant effect of body mass on the overall substitution rate, the synonymous vs. non-synonymous substitution rate and the dynamics of GC-content. Our molecular comparative analysis points toward relatively small Cetartiodactyla ancestors, in agreement with the fossil record, even though our data set almost exclusively consists of large species. This analysis demonstrates the potential of phylogenomic methods for ancestral trait reconstruction and gives credit to recent suggestions that the ancestor to placental mammals was a relatively large and long-lived animal.

In the North-east (NE) Atlantic, most intertidal fucoids and warm-temperate kelps show unique low-latitude gene pools matching long-term climatic refugia. For cold-temperate kelps data are scarcer despite their unique cultural, ecological and economic significance. Here we test whether the amphi-Atlantic range ofLaminaria digitatais derived from past glacial survival (and vicariance) in both NE and North-west (NW) Atlantic refugia (as suggested by niche modelling), or post-glacial (re)colonization (as suggested by low mtDNA divergence). We screened 14 populations from across the species range for 12 microsatellite loci to identify and map major gene pools and refugia. We assessed if NW Atlantic survival was supported by unique endemic variation, and if genetic diversity and structure were, as predicted from larger hindcasted glacial ranges, higher in the NE Atlantic. Microsatellite data subdividedL. digitatainto three main genetic groups matching Brittany, northern Europe and the NW Atlantic, with finer-scale sub-structuring within European clusters. The relatively diverse NE Atlantic lineages probably survived the Last Glacial Maximum along unglaciated periglacial shorelines of the Armorican and Celtic Seas (Brittany cluster) and Ireland (northern European cluster), and remain well differentiated despite their relative proximity. The unique Brittany gene pool, at the contemporary European rear edge, is projected to disappear in the near future under high greenhouse gas emission scenarios. Low allelic diversity and low endemism in the NW Atlantic are consistent with recent post-glacial colonization from Europe, challenging the long-standing hypothesis ofin situglacial survival. Confusion withHedophyllum nigripesmay have led to underestimation of regional diversity ofL. digitata, but also to overestimation of its presence along putative trans-Atlantic migration routes. Partial incongruence between modelling and genetic-based biogeographic inferences highlights the benefits of comparing both approaches to understand how shifting climatic conditions affect marine species distributions and explain large-scale patterns of spatial genetic structure.

Temperature is one of the most important range-limiting factors for many seaweeds. Driven by the recent climatic changes, rapid northward shifts of species’ distribution ranges can potentially modify the phylogeographic signature of Last Glacial Maximum. We explored this question in detail in the cold-tolerant kelp species Saccharina latissima, using microsatellites and double digest restriction site-associated DNA sequencing ( ddRAD-seq) derived single nucleotide polymorphisms (SNPs) to analyze the genetic diversity and structure in 11 sites spanning the entire European Atlantic latitudinal range of this species. In addition, we checked for statistical correlation between genetic marker allele frequencies and three environmental proxies (sea surface temperature, salinity, and water turbidity). Our findings revealed that genetic diversity was significantly higher for the northernmost locality (Spitsbergen) compared to the southern ones (Northern Iberia), which we discuss in light of the current state of knowledge on phylogeography of S. latissima and the potential influence of the recent climatic changes on the population structure of this species. Seven SNPs and 12 microsatellite alleles were found to be significantly associated with at least one of the three environmental variables. We speculate on the putative adaptive functions of the genes associated with the outlier markers and the importance of these markers for successful conservation and aquaculture strategies for S. latissima in this age of rapid global change.

Biological invasions are occurring at increasing rates since the onset of the twentieth century. While ports and marinas have been identified as a major point-of-entry for the spread of marine non-indigenous species (NIS), their relationships with wild habitats however needs further scrutiny. We had the rare opportunity to monitor the real-time colonization dynamics of a newly-built marina by the notorious invasive kelp Undaria pinnatifida in the Bay-of-Morlaix, Brittany (France). Field surveys (> 20,000 individuals geo-localized) were combined with genetic analyses (10 microsatellite loci, N = 890 individuals) over 3 years (i.e., 6 generations in the study area). Regarding the colonization dynamics, a dramatic snow-ball effect was documented over time, with local density reaching up to 50 individuals per m after 2 years. Assignment tests showed that the primary colonizers came from neighboring populations established in natural rocky reefs. A shift towards a self-sustaining population was however observed the following year, with 44% of self-assignment. These processes are best explained by (i) life history traits, notably rapid growth and selfing, and (ii) natural dispersal within the marina combined with human-mediated dispersal—through leisure boating- over longer distances. Spill-over effects have been previously documented, and here also reported, to explain the expansion of U. pinnatifida from marinas to the wild. We showed that the on-going ocean sprawl also offers a perfect arena for spill-back events (i.e., spread from natural habitats to artificial structures), highlighting the need for careful surveillance of newly built infrastructures.

Ports and farms are well-known primary introduction hot spots for marine non-indigenous species (NIS). The extent to which these anthropogenic habitats are sustainable sources of propagules and influence the evolution of NIS in natural habitats was examined in the edible seaweed Undaria pinnatifida, native to Asia and introduced to Europe in the 1970s. Following its deliberate introduction 40 years ago along the French coast of the English Channel, this kelp is now found in three contrasting habitat types: farms, marinas and natural rocky reefs. In the light of the continuous spread of this NIS, it is imperative to better understand the processes behind its sustainable establishment in the wild. In addition, developing effective management plans to curtail the spread of U. pinnatifida requires determining how the three types of populations interact with one another. In addition to an analysis using microsatellite markers, we developed, for the first time in a kelp, a ddRAD-sequencing technique to genotype 738 individuals sampled in 11 rocky reefs, 12 marinas, and two farms located along ca. 1,000 km of coastline. As expected, the RAD-seq panel showed more power than the microsatellite panel for identifying fine-grained patterns. However, both panels demonstrated habitat-specific properties of the study populations. In particular, farms displayed very low genetic diversity and no inbreeding conversely to populations in marinas and natural rocky reefs. In addition, strong, but chaotic regional genetic structure, was revealed, consistent with human-mediated dispersal (e.g., leisure boating). We also uncovered a tight relationship between populations in rocky reefs and those in nearby marinas, but not with nearby farms, suggesting spillover from marinas into the wild. At last, a temporal survey spanning 20 generations showed that wild populations are now self-sustaining, albeit there was no evidence for local adaptation to any of the three habitats. These findings highlight that limiting the spread of U. pinnatifida requires efficient management policies that also target marinas.

Morphological changes-such as dark spots, twisted stipes and deformed blades-have been observed in wild and cultivated Saccharina latissima. The putative cause for the disease symptoms is the filamentous endophytic brown alga Laminarionema elsbetiae, which is known to invade stipes and fronds of its hosts. Little is known about this interaction and its occurrence in the field, although former studies indicated high endophyte prevalence in kelp populations. Previous epidemiological studies on kelp endophytes were mainly based on the examination of microscopic sections, followed by time-consuming isolation and cultivation steps in order to identify the endophyte and a reliable method to quantify endophyte infections was missing. As a novel approach, we established and validated a qPCR assay for relative quantification of the endophyte L. elsbetiae within its host S. latissima, which allows to examine both, the prevalence of endophytic algae and the severity of infections. The assay was shown to be highly specific and suitable to reliably detect small amounts of endophyte DNA in the host. Using this method, we detected very high endophyte prevalence in the investigated kelp populations, up to 100% in young S. latissima sporophytes in Brittany during spring. Furthermore, our results suggest that Saccharina sporophytes are infected early in their life and that seasonality and environmental factors have a significant impact on infection rates. In the future, this approach could also be applied to study other host-endophyte pairs using specific primers.

BACKGROUND:Contamination is a well-known but often neglected problem in molecular biology. Here, we investigated the prevalence of cross-contamination among 446 samples from 116 distinct species of animals, which were processed in the same laboratory and subjected to subcontracted transcriptome sequencing.RESULTS:Using cytochrome oxidase 1 as a barcode, we identified a minimum of 782 events of between-species contamination, with approximately 80% of our samples being affected. An analysis of laboratory metadata revealed a strong effect of the sequencing center: nearly all the detected events of between-species contamination involved species that were sent the same day to the same company. We introduce new methods to address the amount of within-species, between-individual contamination, and to correct for this problem when calling genotypes from base read counts.CONCLUSIONS:We report evidence for pervasive within-species contamination in this data set, and show that classical population genomic statistics, such as synonymous diversity, the ratio of non-synonymous to synonymous diversity, inbreeding coefficient FIT, and Tajima's D, are sensitive to this problem to various extents. Control analyses suggest that our published results are probably robust to the problem of contamination. Recommendations on how to prevent or avoid contamination in large-scale population genomics/molecular ecology are provided based on this analysis.

The fixation probability of a recessive beneficial mutation is increased on the X or Z chromosome, relative to autosomes, because recessive alleles carried by X or Z are exposed to selection in the heterogametic sex. This leads to an increased dN/dS ratio on sex chromosomes relative to autosomes, a pattern called the "fast-X" or "fast-Z" effect. Besides positive selection, the strength of genetic drift and the efficacy of purifying selection, which affect the rate of molecular evolution, might differ between sex chromosomes and autosomes. Disentangling the complex effects of these distinct forces requires the genome-wide analysis of polymorphism, divergence and gene expression data in a variety of taxa. Here we study the influence of hemizygosity of the Z chromosome in Maniola jurtina and Pyronia tithonus, two species of butterflies (Lepidoptera, Nymphalidae, Satyrinae). Using transcriptome data, we compare the strength of positive and negative selection between Z and autosomes accounting for sex-specific gene expression. We show that M. jurtina and P. tithonus do not experience a faster, but rather a slightly slower evolutionary rate on the Z than on autosomes. Our analysis failed to detect a significant difference in adaptive evolutionary rate between Z and autosomes, but comparison of male-biased, unbiased and female-biased Z-linked genes revealed an increased efficacy of purifying selection against recessive deleterious mutations in female-biased Z-linked genes. This probably contributes to the lack of fast-Z evolution of satyrines. We suggest that the effect of hemizygosity on the fate of recessive deleterious mutations should be taken into account when interpreting patterns of molecular evolution in sex chromosomes vs. autosomes.

Comparative population genetics in asexual vs. sexual species offers the opportunity to investigate the impact of asexuality on genome evolution. Here, we analyse coding sequence polymorphism and divergence patterns in the fascinating Lineus ribbon worms, a group of marine, carnivorous nemerteans with unusual regeneration abilities, and in which asexual reproduction by fissiparity is documented. The population genomics of the fissiparous L. pseudolacteus is characterized by an extremely high level of heterozygosity and unexpectedly elevated π N /π S ratio, in apparent agreement with theoretical expectations under clonal evolution. Analysis of among‐species allele sharing and read‐count distribution, however, reveals that L. pseudolacteus is a triploid hybrid between Atlantic populations of L. sanguineus and L. lacteus . We model and quantify the relative impact of hybridity, polyploidy and asexuality on molecular variation patterns in L. pseudolacteus and conclude that (i) the peculiarities of L. pseudolacteus population genomics result in the first place from hybridization and (ii) the accumulation of new mutations through the Meselson effect is more than compensated by processes of heterozygosity erosion, such as gene conversion or gene copy loss. This study illustrates the complexity of the evolutionary processes associated with asexuality and identifies L. pseudolacteus as a promising model to study the first steps of polyploid genome evolution in an asexual context.

Termites are eusocial insects related to cockroaches that feed on lignocellulose. These insects are key species in ecosystems since they recycle a large amount of nutrients but also are pests, exerting major economic impacts. Knowledge on the molecular pathways underlying reproduction, caste differentiation or lignocellulose digestion would largely benefit from additional transcriptomic data. This study focused on transcriptomes of secondary reproductive females (nymphoid neotenics). Thirteen transcriptomes were used: 10 of Reticulitermes flavipes and R. grassei sequenced from a previous study, and two transcriptomes of R. lucifugus sequenced for the present study. After transcriptome assembly and read mapping, we examined interspecific variations of genes expressed by termites or gut microorganisms. A total of 18,323 orthologous gene clusters were detected. Functional annotation and taxonomic assignment were performed on a total of 41,287 predicted contigs in the three termite species. Between the termite species studied, functional categories of genes were comparable. Gene ontology (GO) terms analysis allowed the discovery of 9 cellulases and a total of 79 contigs potentially involved in 11 enzymatic activities used in wood metabolism. Altogether, results of this study illustrate the strong potential for the use of comparative interspecific transcriptomes, representing a complete resource for future studies including differentially expressed genes between castes or SNP analysis for population genetics.

Because mating systems affect population genetics and ecology, they are expected to impact the molecular evolution of species. Self-fertilizing species experience reduced effective population size, recombination rates, and heterozygosity, which in turn should decrease the efficacy of natural selection, both adaptive and purifying, and the strength of meiotic drive processes such as GC-biased gene conversion. The empirical evidence is only partly congruent with these predictions, depending on the analyzed species, some, but not all, of the expected effects have been observed. One possible reason is that self-fertilization is an evolutionary dead-end, so that most current selfers recently evolved self-fertilization, and their genome has not yet been strongly impacted by selfing. Here, we investigate the molecular evolution of two groups of freshwater snails in which mating systems have likely been stable for several millions of years. Analyzing coding sequence polymorphism, divergence, and expression levels, we report a strongly reduced genetic diversity, decreased efficacy of purifying selection, slower rate of adaptive evolution, and weakened codon usage bias/GC-biased gene conversion in the selfer Galba compared with the outcrosser Physa, in full agreement with theoretical expectations. Our results demonstrate that self-fertilization, when effective in the long run, is a major driver of population genomic and molecular evolutionary processes. Despite the genomic effects of selfing, Galba truncatula seems to escape the demographic consequences of the genetic load. We suggest that the particular ecology of the species may buffer the negative consequences of selfing, shedding new light on the dead-end hypothesis.

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life-history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome-wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA-seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes – a marker of the effective population size – in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.