Search

The deep ocean is the largest biome on Earth and yet it is among the least studied environments of our planet. Life at great depths requires several specific adaptations; however, their molecular mechanisms remain understudied. We examined patterns of positive selection in 416 genes from four brittle star (Ophiuroidea) families displaying replicated events of deep-sea colonization (288 individuals from 216 species). We found consistent signatures of molecular convergence in functions related to protein biogenesis, including protein folding and translation. Five genes were recurrently positively selected, including chaperonin-containing TCP-1 subunit α (CCTα), which is essential for protein folding. Molecular convergence was detected at the functional and gene levels but not at the amino-acid level. Pressure-adapted proteins are expected to display higher stability to counteract the effects of denaturation. We thus examined in silico local protein stability of CCTα across the ophiuroid tree of life (967 individuals from 725 species) in a phylogenetically corrected context and found that deep-sea-adapted proteins display higher stability within and next to the substrate-binding region, which was confirmed by in silico global protein stability analyses. This suggests that CCTα displays not only structural but also functional adaptations to deep-water conditions. The CCT complex is involved in the folding of ∼10% of newly synthesized proteins and has previously been categorized as a “cold-shock” protein in numerous eukaryotes. We thus propose that adaptation mechanisms to cold and deep-sea environments may be linked and highlight that efficient protein biogenesis, including protein folding and translation, is a key metabolic deep-sea adaptation.

Reticulated evolution -i.e. secondary introgression/admixture between sister taxa- is increasingly recognized as a key evolutionary process that may play a role in structuring infra-specific genetic variation, and possibly promoting adaptation. Mytilus spp. is an ideal system to assess its importance, because these marine mussels form semi-isolated species that remain reproductively compatible over large time-scales. It includes three taxa that hybridize in the Northern Hemisphere (M. edulis, M. galloprovincialis and M. trossulus) and two taxa of uncertain ancestry in the Southern Hemisphere (M. platensis: South America and the Kerguelen Islands; and M. planulatus: Australasia). The Kerguelen mussels are of particular interest to investigate the potential role of admixture in enhancing micro-geographic structure, as they inhabit a small and isolated island in the Southern Ocean characterized by a highly heterogeneous environment, and genomic reticulation between Northern and Southern lineages has been suspected. Here, we extended a previous analysis by using targeted-sequencing data (51,878 SNPs) across the three Northern species and the Kerguelen population, coupled with a panel of 33 SNPs genotyped on 695 mussels across 35 sites in the Kerguelen Islands. The panel was enriched with ancestry-informative SNPs, i.e. SNPs that were more differentiated than the genomic average between Northern lineages, to evaluate whether reticulated evolution contributed to micro-geographic structure. We first showed that the Kerguelen population belongs to a divergent Southern lineage, most related to M. edulis mussels, that experienced secondary admixture with non-indigenous Northern species. We then demonstrated that the Kerguelen mussels were significantly differentiated over small spatial distance, and that this local genetic structure was associated with environmental variation and mostly revealed by ancestry-informative markers. Although local adaptation can explain the association with the environment we believe it more likely that environment variables better describe population connectivity than geographic distance. Our study highlights genetic connectivity of populations is more easily revealed by non-equilibrium secondary introgression clines at a subset of loci, while association with the environment should not be hastily advocated to support adaptation from admixture variation.

Why sex chromosomes turn over and remain undifferentiated in some taxa, whereas they degenerate in others, is still an area of ongoing research. The recurrent occurrence of homologous and homomorphic sex chromosomes in distantly related taxa suggests their independent evolution or continued recombination since their first emergence. Fishes display a great diversity of sex-determining systems. Here, we focus on sex chromosome evolution in haplochromines, the most species-rich lineage of cichlid fishes. We investigate sex-specific signatures in the Pseudocrenilabrus philander species complex, which belongs to a haplochromine genus found in many river systems and ichthyogeographic regions in northern, eastern, central, and southern Africa. Using whole-genome sequencing and population genetic, phylogenetic, and read-coverage analyses, we show that one population of P. philander has an XX–XY sex-determining system on LG7 with a large region of suppressed recombination. However, in a second bottlenecked population, we did not find any sign of a sex chromosome. Interestingly, LG7 also carries an XX–XY system in the phylogenetically more derived Lake Malawi haplochromine cichlids. Although the genomic regions determining sex are the same in Lake Malawi cichlids and P. philander, we did not find evidence for shared ancestry, suggesting that LG7 evolved as sex chromosome at least twice in haplochromine cichlids. Hence, our work provides further evidence for the labile nature of sex determination in fishes and supports the hypothesis that the same genomic regions can repeatedly and rapidly be recruited as sex chromosomes in more distantly related lineages.

Understanding the drivers and dynamics of diversification is a central topic in evolutionary biology. Here, we investigated the dynamics of diversification in the cichlid fish Astatotilapia burtoni that diverged along a lake-stream environmental gradient. Whole-genome and morphometric analyses revealed that divergent selection was essential at the early stages of diversification, but that periods in allopatry were likely involved towards the completion of speciation. While morphological differentiation was continuous, genomic differentiation was not, as shown by two clearly separated categories of genomic differentiation. Reproductive isolation increased along a continuum of genomic divergence, with a “grey zone” of speciation at ∼0.1% net nucleotide divergence. The quantification of the extent of (non-)parallelism in nine lake-stream population pairs from four cichlid species by means of multivariate analyses revealed one parallel axis of genomic and morphological differentiation among seven lake-stream systems. Finally, we found that parallelism was higher when ancestral lake populations were more similar.

Description, African cichlid fishes reveal the dynamics of speciation and the predictability of evolution., Understanding the dynamics of speciation is a central topic in evolutionary biology. Here, we investigated how morphological and genomic differentiation accumulated along the speciation continuum in the African cichlid fish Astatotilapia burtoni. While morphological differentiation was continuously distributed across different lake-stream population pairs, we found that there were two categories with respect to genomic differentiation, suggesting a “gray zone” of speciation at ~0.1% net nucleotide divergence. Genomic differentiation was increased in the presence of divergent selection and drift compared to drift alone. The quantification of phenotypic and genetic parallelism in four cichlid species occurring along a lake-stream environmental contrast revealed parallel and antiparallel components in rapid adaptive divergence, and morphological convergence in species replicates inhabiting the same environments. Furthermore, we show that the extent of parallelism was higher when ancestral populations were more similar. Our study highlights the complementary roles of divergent selection and drift on speciation and parallel evolution.

International audience; The conspicuous Mediterranean brittle star Ophioderma longicauda (Bruzelius, 1805) has been discovered to represent a cryptic species complex, consisting of six nuclear clusters with contrasting reproductive modes (broadcast spawners and brooders). Here, O. longicauda is re-described. It is distinguished by a dark reddish-brown colouration both dorsally and on the ventral disc, and multiple tumid dorsal arm plates. One eastern Mediterranean brooding cluster is described as O. zibrowii sp. nov., characterized by a dark olive-green colour both dorsally and on the ventral disc, and single dorsal arm plates. Another brooder is described from Tunisia as O. hybrida sp. nov., with a highly variable morphology that reflects its origin by hybridization of O. longicauda and a brooder (possibly O. zibrowii sp. nov.), leaving the third brooding cluster as morphologically indistinguishable at this point and possibly conspecific with one of the others. The West-African O. guineense Greef, 1882 is resurrected as a valid species, differing morphologically from O. longicauda by predominantly single dorsal arm plates and light green or creamy white ventral side. Also from West Africa, O. africana sp. nov. is described, characterized by a dark brown colour, dorsally and ventrally, and single dorsal arm plates.

The deep ocean is the largest biome on Earth and yet it is among the least studied environments of our planet. Life at great depths requires several specific adaptations, however their molecular mechanisms remain understudied. We examined patterns of positive selection in 416 genes from four brittle star (Ophiuroidea) families displaying replicated events of deep-sea colonization (288 individuals from 216 species). We found consistent signatures of molecular convergence in functions related to protein biogenesis, including protein folding and translation. Five genes were recurrently positively selected, including CCTα (Chaperonin Containing TCP-1 subunit α), which is essential for protein folding. Molecular convergence was detected at the functional and gene levels but not at the amino-acid level. Pressure-adapted proteins are expected to display higher stability to counteract the effects of denaturation. We thus examinedin silicolocal protein stability of CCTα across the ophiuroid tree of life (967 individuals from 725 species) in a phylogenetically-corrected context and found that deep sea-adapted proteins display higher stability within and next to the substrate-binding region, which was confirmed byin silicoglobal protein stability analyses. This suggests that CCTα not only displays structural but also functional adaptations to deep water conditions. The CCT complex is involved in the folding of ∼10% of newly synthesized proteins and has previously been categorized as ‘cold-shock’ protein in numerous eukaryotes. We thus propose that adaptation mechanisms to cold and deep-sea environments may be linked and highlight that efficient protein biogenesis, including protein folding and translation, are key metabolic deep-sea adaptations.

Understanding how reproductive barriers evolve and which barriers contribute to speciation requires the examination of organismal lineages that are still in the process of diversification and the study of the full range of reproductive barriers acting at different life stages. Lake and river ecotypes of the East African cichlid fish Astatotilapia burtoni show habitat-specific adaptations, despite different levels of genetic differentiation, and thus represent an ideal model to study the evolution of reproductive barriers. To evaluate the degree of reproductive isolation between genetically divergent lake and river populations, we performed a mesocosm mating experiment in a semi-natural setting at Lake Tanganyika. We assessed reproductive isolation in the presence of male-male competition by analyzing survival and growth rates of introduced adults and their reproductive success from genetic parentage of surviving offspring. The genetically divergent river population showed reduced fitness in terms of survival, growth rate, and mating success in a lake-like environment. Hybrid offspring between different populations showed intermediate survival consistent with extrinsic postzygotic reproductive barriers. Our results suggest that both prezygotic (immigrant inviability) and postzygotic reproductive barriers contribute to divergence, and highlight the value of assessing multiple reproductive barriers acting at different stages and in natural contexts to understand speciation mechanisms.

International audience; Closely related species are key models to investigate mechanisms leading to reproductive isolation and early stages of diversification, also at the genomic level. The brittle star cryptic species complex Ophioderma longicauda encompasses the sympatric broadcast-spawning species C3 and the internal brooding species C5. Here, we used de novo transcriptome sequencing and assembly in two closely related species displaying contrasting reproductive modes to compare their genetic diversity and to investigate the role of natural selection in reproductive isolation. We reconstructed 20 146 and 22 123 genes for C3 and C5, respectively, and characterized a set of 12 229 orthologs. Genetic diversity was 1.5-2 times higher in C3 compared to C5, confirming that species with low parental investment display higher levels of genetic diversity. Forty-eight genes were the targets of positive diversifying selection during the evolution of the two species. Notably, two genes (NHE and TetraKCNG) are sperm-specific ion channels involved in sperm motility. Ancestral sequence reconstructions show that natural selection targeted the two genes in the brooding species. This may result from an adaptation to the novel environmental conditions surrounding sperm in the brooding species, either directly affecting sperm or via an increase in male/female conflict. This phenomenon could have promoted prezygotic reproductive isolation between C3 and C5. Finally, the sperm receptors to egg chemoattractants differed between C3 and C5 in the ligand-binding region. We propose that mechanisms of species-specific gamete recognition in brittle stars occur during sperm chemotaxis (sperm attraction towards the eggs), contrary to other marine invertebrates where prezygotic barriers to interspecific hybridization typically occur before sperm-egg fusion.

Cryptic species are numerous in the marine environment. The brittle star Ophioderma longicauda is composed of six mitochondrial lineages, encompassing brooders, which form a monophyletic group, and broadcasters, from which the brooders are derived. To clarify the species limits within O. longicauda, we compared the reproductive status of the sympatric lineages L1 and L3 (defined after sequencing a portion of the mitochondrial gene COI) during the month of May in Greece. In addition, we genotyped a nuclear marker, intron i51. Each L3 female was brooding, whereas all L1 specimens displayed full gonads, suggesting temporal pre-zygotic isolation between brooders and broadcasters. Statistical differences were found among lineages in morphology and bathymetric distribution. Finally, the intron i51 was polymorphic in L1 (60 individuals), but monomorphic in L3 (109 individuals), confirming the absence of gene flow between brooders and broadcasters. In conclusion, the broadcasting lineage L1 and the brooding lineage L3 are different biological species., Les espèces cryptiques sont nombreuses dans l’environnement marin. L’ophiure Ophioderma longicauda est composée de six lignées mitochondriales, comprenant des incubants, qui forment un groupe monophylétique, et des individus à larves, desquels les incubants sont dérivés. Afin de clarifier les limites d’espèces chez O. longicauda, nous avons comparé l’état de reproduction des lignées sympatriques L1 et L3 (définies après séquençage d’une partie du gène mitochondrial COI) durant le mois de mai en Grèce. De plus, nous avons génotypé un marqueur nucléaire, l’intron i51. Chaque femelle L3 était incubante tandis que tous les individus L1 avaient les gonades pleines, suggérant un isolement pré-zygotique temporel entre les incubants et les individus à larves. Des différences statistiques de morphologie et de distribution bathymétrique ont été trouvées entre lignées. Finalement, l’intron i51 était polymorphe chez L1 (60 individus), mais monomorphe chez L3 (109 individus), confirmant l’absence de flux de gènes entre les incubants et dispersants. En conclusion, la lignée à larves L1 et la lignée incubante L3 sont des espèces biologiques différentes.

Ribosomal DNA is commonly used as a marker for protist phylogeny and taxonomy because of its ubiquity and its expected species specificity thanks to the mechanism of concerted evolution. However, numerous studies reported the occurrence of intragenomic (intra-individual) polymorphism in various protists and particularly in Foraminifera. To infer to what extent the SSU rDNA intragenomic variability occurs in Foraminifera, we studied 16 foraminiferal species belonging to single-chambered monothalamids and multi-chambered Globothalamea, with one to six individuals per species. We performed single-cell DNA extractions and PCRs of a 600 bp fragment of SSU rDNA, and sequenced 9 to 23 clones per individual for a total of 818 sequences. We found intragenomic variability in almost all species, even after excluding singleton mutations. Intra-individual sequence divergence ranged from 0 to 5.15% and was higher than 1% in 11 species. Variability was usually located at the end of stem-loop structures and included compensatory single nucleotide polymorphisms and expansion segments polymorphisms. However, the polymorphisms did not change the secondary structure of the rRNA. Our results suggest a non-concerted evolution of rRNA genes in Foraminifera. The origin of this variability and its implications for species identification in environmental DNA studies are discussed.

International audience; Genetic diversity is the amount of variation observed between DNA sequences from distinct individuals of a given species. This pivotal concept of population genetics has implications for species health, domest-ication, management and conservation. Levels of genetic diversity seem to vary greatly in natural populations and species, but the determinants of this variation, and particularly the relative influences of species biology and ecology versus population history, are still largely mysterious 1,2. Here we show that the diversity of a species is predictable , and is determined in the first place by its ecological strategy. We investigated the genome-wide diversity of 76 non-model animal species by sequencing the transcriptome of two to ten individuals in each species. The distribution of genetic diversity between species revealed no detectable influence of geographic range or invasive status but was accurately predicted by key species traits related to parental investment: long-lived or low-fecundity species with brooding ability were genetically less diverse than short-lived or highly fecund ones. Our analysis demonstrates the influence of long-term life-history strategies on species response to short-term environmental perturbations, a result with immediate implications for conservation policies. Since the early studies of evolutionary genetics, there has been no understanding of how and why genetic diversity levels vary between species. This old puzzle, considered four decades ago as 'the central problem in population genetics'

Protists are key players in microbial communities, yet our understanding of their role in ecosystem functioning is seriously impeded by difficulties in identification of protistan species and their quantification. Current microscopy-based methods used for determining the abundance of protists are tedious and often show a low taxonomic resolution. Recent development of next-generation sequencing technologies offered a very powerful tool for studying the richness of protistan communities. Still, the relationship between abundance of species and number of sequences remains subjected to various technical and biological biases. Here, we test the impact of some of these biological biases on sequence abundance of SSU rRNA gene in foraminifera. First, we quantified the rDNA copy number and rRNA expression level of three species of foraminifera by qPCR. Then, we prepared five mock communities with these species, two in equal proportions and three with one species ten times more abundant. The libraries of rDNA and cDNA of the mock communities were constructed, Sanger sequenced and the sequence abundance was calculated. The initial species proportions were compared to the raw sequence proportions as well as to the sequence abundance normalized by rDNA copy number and rRNA expression level per species. Our results showed that without normalization, all sequence data differed significantly from the initial proportions. After normalization, the congruence between the number of sequences and number of specimens was much better. We conclude that without normalization, species abundance determination based on sequence data was not possible because of the effect of biological biases. Nevertheless, by taking into account the variation of rDNA copy number and rRNA expression level we were able to infer species abundance, suggesting that our approach can be successful in controlled conditions.