When does a parasite become a disease? eDNA unravels complex host-pathogen dynamics across environmental stress gradients in wild salmonid populations.

Infectious diseases stem from disrupted interactions among hosts, parasites, and the environment. Both abiotic and biotic factors can influence infection outcomes by shaping the abundance of a parasite's infective stages, as well as the host's ability to fight infection. However, disentangling these mechanisms within natural ecosystems remains challenging. Here, combining environmental DNA analysis and niche modelling at a regional scale, we uncovered the biotic and abiotic drivers of an infectious disease of salmonid fish, triggered by the parasite Tetracapsuloides bryosalmonae. We found that the occurrence and abundance of the parasite in the water-i.e., the propagule pressure- were mainly correlated to the abundances of its two primary hosts, the bryozoan Fredericella sultana and the fish Salmo trutta, but poorly to local abiotic environmental stressors. In contrast, the occurrence and abundance of parasites within fish hosts-i.e., proxies for disease emergence-were closely linked to environmental stressors (water temperature, agricultural activities, dams), and to a lesser extent to parasite propagule pressure. These results suggest that pathogen distribution alone cannot predict the risk of disease in wildlife, and that local anthropogenic stressors may play a pivotal role in disease emergence among wild host populations, likely by modulating the hosts' immune response. Our study sheds light on the intricate interplay between biotic and abiotic factors in shaping pathogen distribution and raises concerns about the effects of global change on pathogen emergence.

[Brucellosis in Alpine ibex: 10 years of research and expert assessments]. / La brucellose du bouquetin des Alpes - Un exemple de dix années de recherche et d'expertise.

Brucellosis due to Brucella melitensis affects domestic and wild ruminants, as well as other mammals, including humans. Despite France being officially free of bovine brucellosis since 2005, two human cases of Brucella melitensis infection in the French Alps in 2012 led to the discovery of one infected cattle herd and of one infected population of wild Alpine ibex (Capra ibex). In this review, we present the results of 10 years of research on the epidemiology of brucellosis in this population of Alpine ibex. We also discuss the insights brought by research and expert assessments on the efficacy of disease management strategies used to mitigate brucellosis in the French Alps.

Title: La brucellose du bouquetin des Alpes - Un exemple de dix années de recherche et d'expertise. Abstract: La brucellose à Brucella melitensis touche les ruminants domestiques et sauvages, ainsi que d'autres mammifères, dont les humains. Bien que la France soit officiellement indemne depuis 2005, deux cas humains reportés en Haute-Savoie en 2012 ont conduit à la découverte de l'infection dans un élevage bovin et chez les bouquetins des Alpes (Capra ibex) du massif du Bargy. Nous présentons dans cette synthèse les principales découvertes de ces dix dernières années sur le système brucellose-bouquetins. Nous discuterons également de l'apport de la recherche et de l'expertise sur l'évaluation de l'efficacité des mesures de gestion sanitaire mises en place dans le massif du Bargy pour lutter contre la brucellose.

Quantifying heritability and estimating evolutionary potential in the wild when individuals that share genes also share environments.

Accurate heritability estimates for fitness-related traits are required to predict an organism's ability to respond to global change. Heritability estimates are theoretically expected to be inflated if, due to limited dispersal, individuals that share genes are also likely to share similar environments. However, if relatives occupy similar environments due, at least partly, to genetic variation for habitat selection, then accounting for environmental similarity in quantitative genetic models may result in diminished heritability estimates in wild populations. This potential issue has been pointed out in the literature, but has not been evaluated by empirical studies. Here, we investigate whether environmental similarity among individuals can be partly explained by genetic variation for habitat selection, and how this link potentially blurs estimates for heritability in fitness-related traits. Using intensive GPS monitoring, we quantified home-range habitat composition for 293 roe deer inhabiting a heterogeneous landscape to assess environmental similarity. To investigate if environmental similarity might harbour genetic variation, we combined genome-wide data in a quantitative genetic framework to evaluate genetic variation for home-range habitat composition, which is partly the result of habitat selection at settlement. Finally, we explored how environmental similarity affects heritability estimates for behaviours related to the risk avoidance-resource acquisition trade-off (i.e. being in open habitat and distance to roads) and proxies of individual performance (i.e. body mass and hind foot length). We found substantial heritability for home-range habitat composition, with estimates ranging from 0.40 (proportion of meadows) to 0.85 (proportion of refuge habitat). Accounting for similarity in habitat composition between relatives decreased the heritability estimates for both behavioural and morphological traits (reduction ranging from 55% to 100% and from 22% to 41% respectively). As a consequence, only half of these heritability estimates remained significantly different from zero. Our results show that similar genotypes occupy similar environments, which could lead to heritable variation being incorrectly attributed to environmental effects. To accurately distinguish the sources of phenotypic variation and predict the ability of organisms to respond to global change, it is necessary to develop quantitative genetic studies investigating the mechanisms underpinning environmental similarity among relatives.

Pathogen-mediated selection favours the maintenance of innate immunity gene polymorphism in a widespread wild ungulate.

Toll-like receptors (TLR) play a central role in recognition and host frontline defence against a wide range of pathogens. A number of recent studies have shown that TLR genes (Tlrs) often exhibit large polymorphism in natural populations. Yet, there is little knowledge on how this polymorphism is maintained and how it influences disease susceptibility in the wild. In previous work, we showed that some Tlrs exhibit similarly high levels of genetic diversity as genes of the Major Histocompatibility Complex (MHC), and signatures of contemporary balancing selection in roe deer (Capreolus capreolus), the most abundant cervid species in Europe. Here, we investigated the evolutionary mechanisms by which pathogen-mediated selection could shape this innate immunity genetic diversity by examining the relationships between Tlr (Tlr2, Tlr4 and Tlr5) genotypes (heterozygosity status and presence of specific alleles) and infections with Toxoplasma and Chlamydia, two widespread intracellular pathogens known to cause reproductive failure in ungulates. We showed that Toxoplasma and Chlamydia exposures vary significantly across years and landscape features with few co-infection events detected and that the two pathogens exert antagonistic selection on Tlr2 polymorphism. By contrast, we found limited support for Tlr heterozygote advantage. Our study confirmed the importance of looking beyond Mhc genes in wildlife immunogenetic studies. It also emphasized the necessity to consider multiple pathogen challenges and their spatiotemporal variation to improve our understanding of vertebrate defence evolution against pathogens.

Maternal effects shape offspring physiological condition but do not senesce in a wild mammal.

In vertebrates, offspring survival often decreases with increasing maternal age. While many studies have reported a decline in fitness-related traits of offspring with increasing maternal age, the study of senescence in maternal effect through age-specific changes in offspring physiological condition is still at its infancy. We assessed the influence of maternal age and body mass on offspring physiological condition in two populations of roe deer (Capreolus capreolus) subjected to markedly different environmental conditions. We measured seven markers to index body condition and characterize the immune profile in 86 fawns which became recently independent of their known-aged mothers. We did not find striking effects of maternal age on offspring physiological condition measured at 8 months of age. This absence of evidence for senescence in maternal effects is likely due to the strong viability selection observed in the very first months of life in this species. Offspring physiological condition was, on the other hand, positively influenced by maternal body mass. Between-population differences in environmental conditions experienced by fawns also influenced their average body condition and immune phenotype. Fawns facing food limitation displayed lower values in some markers of body condition (body mass and haemoglobin levels) than those living in good quality habitat. They also allocated preferentially to humoral immunity, contrary to those living in good conditions, which allocated more to cellular response. These results shed a new light on the eco-physiological pathways mediating the relationship between mother's mass and offspring condition.

Genetic consequences of social structure in the golden-crowned sifaka.

Many species are structured in social groups (SGs) where individuals exhibit complex mating strategies. Yet, most population genetic studies ignore SGs either treating them as small random-mating units or focusing on a higher hierarchical level (the population). Empirical studies acknowledging SGs have found an overall excess of heterozygotes within SGs and usually invoke inbreeding avoidance strategies to explain this finding. However, there is a lack of null models against which ecological theories can be tested and inbreeding avoidance quantified. Here, we investigate inbreeding (deviation from random mating) in an endangered forest-dwelling pair-living lemur species (Propithecus tattersalli). In particular, we measure the inbreeding coefficient (FIS) in empirical data at different scales: SGs, sampling sites and forest patches. We observe high excess of heterozygotes within SGs. The magnitude of this excess is highly dependent on the sampling scheme: while offspring are characterised by a high excess of heterozygotes (FIS < 0), the reproductive pair does not show dramatic departures from Hardy-Weinberg expectations. Moreover, the heterozygosity excess disappears at larger geographic scales (sites and forests). We use a modelling framework that incorporates details of the sifaka mating system but does not include active inbreeding avoidance mechanisms. The simulated data show that, although apparent "random mating" or even inbreeding may occur at the "population" level, outbreeding is maintained within SGs. Altogether our results suggest that social structure leads to high levels of outbreeding without the need for active inbreeding avoidance mechanisms. Thus, demonstrating and measuring the existence of active inbreeding avoidance mechanisms may be more difficult than usually assumed.

Demography and adaptation promoting evolutionary transitions in a mammalian genus that diversified during the Pleistocene.

Species that evolved in temperate regions during the Pleistocene experienced periods of extreme climatic transitions. Consequent population fragmentation and dynamics had the potential to generate small, isolated populations where the influence of genetic drift would be expected to be strong. We use comparative genomics to assess the evolutionary influence of historical demographics and natural selection through a series of transitions associated with the formation of the genus Capreolus, speciation within this genus during the Quaternary and during divergence among European roe deer (C. capreolus) populations. Our analyses were facilitated by the generation of a new high-coverage reference genome for the Siberian roe deer (C. pygargus). We find progressive reductions in effective population size (Ne ), despite very large census sizes in modern C. capreolus populations and show that low Ne has impacted the C. capreolus genome, reducing diversity and increasing linkage disequilibrium. Even so, we find evidence for natural selection shared among C. capreolus populations, including a historically documented founder population that has been through a severe bottleneck. During each phylogenetic transition there is evidence for selection (from dN/dS and nucleotide diversity tests), including at loci associated with diapause (delayed embryonic development), a phenotype restricted to this genus among the even-toed ungulates. Together these data allow us to assess expectations for the origin and diversification of a mammalian genus during a period of extreme environmental change.

Genetic epidemiology of the Alpine ibex reservoir of persistent and virulent brucellosis outbreak.

While it is now broadly accepted that inter-individual variation in the outcomes of host-pathogen interactions is at least partially genetically controlled, host immunogenetic characteristics are rarely investigated in wildlife epidemiological studies. Furthermore, most immunogenetic studies in the wild focused solely on the major histocompatibility complex (MHC) diversity despite it accounts for only a fraction of the genetic variation in pathogen resistance. Here, we investigated immunogenetic diversity of the Alpine ibex (Capra ibex) population of the Bargy massif, reservoir of a virulent outbreak of brucellosis. We analysed the polymorphism and associations with disease resistance of the MHC Class II Drb gene and several non-MHC genes (Toll-like receptor genes, Slc11A1) involved in the innate immune response to Brucella in domestic ungulates. We found a very low neutral genetic diversity and a unique MHC Drb haplotype in this population founded few decades ago from a small number of individuals. By contrast, other immunity-related genes have maintained polymorphism and some showed significant associations with the brucellosis infection status hence suggesting a predominant role of pathogen-mediated selection in their recent evolutionary trajectory. Our results highlight the need to monitor immunogenetic variation in wildlife epidemiological studies and to look beyond the MHC.

Pedigree-free quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer.

Assessing the evolutionary potential of animal populations in the wild is crucial to understanding how they may respond to selection mediated by rapid environmental change (e.g. habitat loss and fragmentation). A growing number of studies have investigated the adaptive role of behaviour, but assessments of its genetic basis in a natural setting remain scarce. We combined intensive biologging technology with genome-wide data and a pedigree-free quantitative genetic approach to quantify repeatability, heritability and evolvability for a suite of behaviours related to the risk avoidance-resource acquisition trade-off in a wild roe deer (Capreolus capreolus) population inhabiting a heterogeneous, human-dominated landscape. These traits, linked to the stress response, movement and space-use behaviour, were all moderately to highly repeatable. Furthermore, the repeatable among-individual component of variation in these traits was partly due to additive genetic variance, with heritability estimates ranging from 0.21 ± 0.08 to 0.70 ± 0.11 and evolvability ranging from 1.1% to 4.3%. Changes in the trait mean can therefore occur under hypothetical directional selection over just a few generations. To the best of our knowledge, this is the first empirical demonstration of additive genetic variation in space-use behaviour in a free-ranging population based on genomic relatedness data. We conclude that wild animal populations may have the potential to adjust their spatial behaviour to human-driven environmental modifications through microevolutionary change.

Brucella melitensis Rev.1 vaccination generates a higher shedding risk of the vaccine strain in Alpine ibex (Capra ibex) compared to the domestic goat (Capra hircus).

Epidemiological investigations implemented in wild and domestic ruminants evidenced a reservoir for Brucella in Capra ibex in the French Alps. Vaccination was considered as a possible way to control Brucella infection in this wildlife population. Twelve ibexes and twelve goats were allocated into four groups housed separately, each including six males or six non-pregnant females. Four to five animals were vaccinated and one or two animals were contact animals. Half of the animals were necropsied 45 days post-vaccination (pv), and the remaining ones at 90 days pv. Additional samples were collected 20 and 68 days pv to explore bacterial distribution in organs and humoral immunity. Neither clinical signs nor Brucella-specific lesions were observed and all vaccinated animals seroconverted. Brucella distribution and antibody profiles were highly contrasted between both species. Proportion of infected samples was significantly higher in ibex compared to goats and decreased between 45 and 90 days pv. Two male ibex presented urogenital excretion at 20 or 45 days pv. The bacterial load was higher 45 days in ibexes compared to goats, whereas it remained moderate to low 90 days pv in both species with large variability between animals. In this experiment, differences between species remained the main source of variation, with low impact of other individual factors. To conclude, multiplicative and shedding capacity of Rev.1 was much higher in ibex compared to goats within 90 days. These results provide initial information on the potential use in natura of a commercial vaccine.

RAD-sequencing for estimating genomic relatedness matrix-based heritability in the wild: A case study in roe deer.

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree-free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long-term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD-sequencing for estimating heritability in a free-ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step-by-step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM-based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix-based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low-frequency variants, affect the relatedness structure of the GRM, generating lower h2 estimates. Our work illustrates the huge potential of RAD-sequencing for estimating GRM-based heritability in virtually any natural population.

Climate change and human colonization triggered habitat loss and fragmentation in Madagascar.

The relative effect of past climate fluctuations and anthropogenic activities on current biome distribution is subject to increasing attention, notably in biodiversity hot spots. In Madagascar, where humans arrived in the last ~4 to 5,000 years, the exact causes of the demise of large vertebrates that cohabited with humans are yet unclear. The prevailing narrative holds that Madagascar was covered with forest before human arrival and that the expansion of grasslands was the result of human-driven deforestation. However, recent studies have shown that vegetation and fauna structure substantially fluctuated during the Holocene. Here, we study the Holocene history of habitat fragmentation in the north of Madagascar using a population genetics approach. To do so, we infer the demographic history of two northern Madagascar neighbouring, congeneric and critically endangered forest dwelling lemur species-Propithecus tattersalli and Propithecus perrieri-using population genetic analyses. Our results highlight the necessity to consider population structure and changes in connectivity in demographic history inferences. We show that both species underwent demographic fluctuations which most likely occurred after the mid-Holocene transition. While mid-Holocene climate change probably triggered major demographic changes in the two lemur species range and connectivity, human settlements that expanded over the last four millennia in northern Madagascar likely played a role in the loss and fragmentation of the forest cover.

Immunogenetic heterogeneity in a widespread ungulate: the European roe deer (Capreolus capreolus).

Understanding how immune genetic variation is shaped by selective and neutral processes in wild populations is of prime importance in both evolutionary biology and epidemiology. The European roe deer (Capreolus capreolus) has considerably expanded its distribution range these last decades, notably by colonizing agricultural landscapes. This range shift is likely to have led to bottlenecks and increased roe deer exposure to a new range of pathogens that until recently predominantly infected humans and domestic fauna. We therefore investigated the historical and contemporary forces that have shaped variability in a panel of genes involved in innate and acquired immunity in roe deer, including Mhc-Drb and genes encoding cytokines or toll-like receptors (TLRs). Together, our results suggest that genetic drift is the main contemporary evolutionary force shaping immunogenetic variation within populations. However, in contrast to the classical view, we found that some innate immune genes involved in micropathogen recognition (e.g. Tlrs) continue to evolve dynamically in roe deer in response to pathogen-mediated positive selection. Most studied Tlrs (Tlr2, Tlr4 and Tlr5) had similarly high levels of amino acid diversity in the three studied populations including one recently established in southwestern France that showed a clear signature of genetic bottleneck. Tlr2 implicated in the recognition of Gram-positive bacteria in domestic ungulates, showed strong evidence of balancing selection. The high immunogenetic variation revealed here implies that roe deer are able to cope with a wide spectrum of pathogens and to respond rapidly to emerging infectious diseases.

Reduced microsatellite heterozygosity does not affect natal dispersal in three contrasting roe deer populations.

Although theoretical studies have predicted a link between individual multilocus heterozygosity and dispersal, few empirical studies have investigated the effect of individual heterozygosity on dispersal propensity or distance. We investigated this link using measures of heterozygosity at 12 putatively neutral microsatellite markers and natal dispersal behaviour in three contrasting populations of European roe deer (Capreolus capreolus), a species displaying pre-saturation condition-dependent natal dispersal. We found no effect of individual heterozygosity on either dispersal propensity or dispersal distance. Average heterozygosity was similar across the three studied populations, but dispersal propensity and distance differed markedly among them. In Aurignac, dispersal propensity and distance were positively related to individual body mass, whereas there was no detectable effect of body mass on dispersal behaviour in Chizé and Trois Fontaines. We suggest that we should expect both dispersal propensity and distance to be greater when heterozygosity is lower only in those species where dispersal behaviour is driven by density-dependent competition for resources.

A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range.

In tropical regions, most primary ecosystems have been replaced by mosaic landscapes in which species must cope with a large shift in the distribution of their habitat and associated food resources. Primates are particularly vulnerable to habitat modifications. Most species persist in small fragments surrounded by complex human-mediated matrices whose structure and connectivity may strongly influence their dispersal and feeding behavior. Behavioral plasticity appears to be a crucial parameter governing the ability of organisms to exploit the resources offered by new matrix habitats and thus to persist in fragmented habitats. In this study, we were interested in the dietary plasticity of the golden-crowned sifaka (Propithecus tattersalli), an endangered species of lemur, found only in the Daraina region in north-eastern Madagascar. We used a DNA-based approach combining the barcoding concept and Illumina next-generation sequencing to (i) describe the species diet across its entire range and (ii) evaluate the influence of landscape heterogeneity on diet diversity and composition. Faeces from 96 individuals were sampled across the entire species range and their contents were analyzed using the trnL metabarcoding approach. In parallel, we built a large DNA reference database based on a checklist of the plant species of the Daraina region. Our results suggest that golden-crowned sifakas exhibit remarkable dietary diversity with at least 130 plant species belonging to 80 genera and 49 different families. We highlighted an influence of both habitat type and openness on diet composition suggesting a high flexibility of foraging strategies. Moreover, we observed the presence of numerous cultivated and naturalized plants in the faeces of groups living in forest edge areas. Overall, our findings support our initial expectation that P. tattersalli is able to cope with the current level of alteration of the landscape and confirm our previous results on the distribution and the dispersal ability of this species.

Genetic data suggest a natural prehuman origin of open habitats in northern Madagascar and question the deforestation narrative in this region.

The impact of climate change and anthropogenic deforestation on biodiversity is of growing concern worldwide. Disentangling how past anthropogenic and natural factors contributed to current biome distribution is thus a crucial issue to understand their complex interactions on wider time scales and to improve predictions and conservation strategies. This is particularly important in biodiversity hotspots, such as Madagascar, dominated by large open habitats whose origins are increasingly debated. Although a dominant narrative argues that Madagascar was originally entirely covered by woodlands, which were destroyed by humans, a number of recent studies have suggested that past climatic fluctuations played a major role in shaping current biome distributions well before humans arrived. Here, we address the question of the origin of open habitats in the Daraina region in northern Madagascar, using a multiproxy approach combining population genetics modeling and remote-sensing analyses. We show that (i) contrary to most regions of Madagascar, the forest cover in Daraina remained remarkably stable over the past 60 y, and (ii) the golden-crowned sifaka (Propithecus tattersalli), a forest-dwelling lemur, underwent a strong population contraction before the arrival of the first humans, hence excluding an anthropogenic cause. Prehuman Holocene droughts may have led to a significant increase of grasslands and a reduction in the species' habitat. This contradicts the prevailing narrative that land cover changes are necessarily anthropogenic in Madagascar but does not preclude the later role played by humans in other regions in which recent lemur bottlenecks have been observed.

Landscape genetics of an endangered lemur (Propithecus tattersalli) within its entire fragmented range.

Habitat fragmentation may strongly reduce individuals' dispersal among resource patches and hence influence population distribution and persistence. We studied the impact of landscape heterogeneity on the dispersal of the golden-crowned sifaka (Propithecus tattersalli), an endangered social lemur species living in a restricted and highly fragmented landscape. We combined spatial analysis and population genetics methods to describe population units and identify the environmental factors which best predict the rates and patterns of genetic differentiation within and between populations. We used non-invasive methods to genotype 230 individuals at 13 microsatellites in all the main forest fragments of its entire distribution area. Our analyses suggest that the Manankolana River and geographical distance are the primary structuring factors, while a national road crossing the region does not seem to impede gene flow. Altogether, our results are in agreement with a limited influence of forest habitat connectivity on gene flow patterns (except for North of the species' range), suggesting that dispersal is still possible today among most forest patches for this species. Within forest patches, we find that dispersal is mainly among neighbouring social groups, hence confirming previous behavioural observations.

Spatial variation in density and total size estimates in fragmented primate populations: the golden-crowned sifaka (Propithecus tattersalli).

The golden-crowned sifaka (Propithecus tattersalli) is an endangered lemur species found only in the Daraina region, a very restricted area in north-eastern Madagascar. Its forest habitat is highly fragmented and expected to suffer from significant changes in the near future. The species is poorly known and only one census study, carried out in 2000, has ever been published. It is thus crucial to update the conservation status of the golden-crowned sifaka before major anthropogenic environmental changes take place. Using the line-transect approach, we estimated the species density in the main forest fragments located in both the peripheral and central parts of the distribution range, including both protected and unprotected areas. In parallel, we tried to determine whether an edge effect could be detected by comparing densities at different distances from the forest edges. We found important variation of sifaka densities among forest fragments. The total species abundance is thus difficult to determine, but we estimated that it is likely to be over 18,000, two to three times higher than previously thought. However, our data also suggested that most P. tattersalli live in forests located in the central part of the distribution range and that the estimated densities in the central part were high (>80 individuals/km(2)). Two forest fragments, found to host a large part of the total population, are currently outside the managed area and their incorporation to the managed area is strongly recommended. Lastly, as expected for a folivorous and not heavily hunted species, our results are consistent with the hypothesis that this species does not experience a clear edge effect, at least during the first half of the dry season. This could be due to a high resiliency to habitat fragmentation or to the fact that fragmentation has been going on for some time.

Molecular phylogeny in mytilids supports the wooden steps to deep-sea vents hypothesis.

Molecular data were used to study the diversity of mytilids associated with sunken-woods sampled in the Solomon Islands and discuss the 'wooden steps to deep-sea vent' hypothesis proposed by Distel et al. First, COI data used in a barcoding approach confirm the presence of four distinct species. Analyses of the 18S rDNA and COI dataset then confirmed that these sunken-wood mytilids belonged to a monophyletic group including all species from deep-sea reducing environments. Finally, we analyzed the relationships within this monophyletic group that include the Bathymodiolinae using a COI dataset and a combined analysis of mitochondrial COI and ND4 genes and nuclear rDNA 18S and 28S. Our study supported the 'wooden steps to deep-sea vent' hypothesis: one of the sunken-wood species had a basal position within the Bathymodiolionae, and all described vent and seep mussels included in our analyses were derived taxa within Bathymodiolinae.