Application of the ecosystem service concept at a small-scale: The cases of coralligenous habitats in the North-western Mediterranean Sea.

The understanding of ecosystem services is essential to support sustainable use and preservation of ecosystems. Coralligenous habitats, main contributors of the Mediterranean marine biodiversity, are yet understudied in term of services provided. This study presents an original small-scale approach to investigate the services provided by coralligenous habitats of a French study area consisting of two marine sites (Marseille and Port-Cros sites) in order to cover two contrasted anthropogenic pressure despite the small-scale. Our results are based on the opinions of 43 experts who ranked 15 services in terms of existence and level of importance for human well-being: supporting ecological functions were considered the most important, then provisioning and cultural services. Regulating services were considered uncertain due to a lack of knowledge. The small-scale approach highlighted a need for a referential frame to determine the existence of services (e.g. geographical and temporal scales, benefits and beneficiaries levels).

Positive selection on sperm ion channels in a brooding brittle star: consequence of life-history traits evolution.

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.

Morphological and genetic analyses reveal a cryptic species complex in the echinoid Echinocardium cordatum and rule out a stabilizing selection explanation.

Preliminary analyses revealed the presence of at least five mitochondrial clades within the widespread sea urchin Echinocardium cordatum (Spatangoida). In this study, we analyzed the genetic (two mitochondrial and two nuclear sequence loci) and morphological characteristics (20 indices) from worldwide samples of this taxon to establish the species limits, morphological diversity and differentiation. Co-occurring spatangoid species were also analyzed with mitochondrial DNA. The nuclear sequences confirm that mitochondrial lineages correspond to true genetic entities and reveal that two clades (named A and B1) hybridize in their sympatry area, although a more closely related pair of clades (B1 and B2), whose distributions widely overlap, does not display hybridization. The morphology of all E. cordatum clade pairs was significantly differentiated, but no morphological diagnostic character was evidenced. By contrast, other spatangoid species pairs that diverged more recently than the E. cordatum clades display clear diagnostic characters. Morphological diversity thus appears responsible for the absence of diagnostic characters, ruling out stabilizing selection, a classical explanation for cryptic species. Alternative classical explanations are (i) environmental plasticity or (ii) a high diversity of genes determining morphology, maintained by varying environmental conditions. We suggest a new hypothesis that the observed morphological diversity is selectively neutral and reflects high effective population sizes in the E. cordatum complex. It is supported by the higher abundance of this taxon compared with other taxa, a trend for the genetic and morphological diversity to be correlated in Europe, and the higher genetic and morphological diversities found in clades of E cordatum (except B1) than in other spatangoid samples in Europe. However, the Pacific clades do not confirm these trends.

Influence of the larval phase on connectivity: strong differences in the genetic structure of brooders and broadcasters in the Ophioderma longicauda species complex.

Closely related species with divergent life history traits are excellent models to infer the role of such traits in genetic diversity and connectivity. Ophioderma longicauda is a brittle star species complex composed of different genetic clusters, including brooders and broadcasters. These species diverged very recently and some of them are sympatric and ecologically syntopic, making them particularly suitable to study the consequences of their trait differences. At the scale of the geographic distribution of the broadcasters (Mediterranean Sea and northeastern Atlantic), we sequenced the mitochondrial marker COI and genotyped an intron (i51) for 788 individuals. In addition, we sequenced 10 nuclear loci newly developed from transcriptome sequences, for six sympatric populations of brooders and broadcasters from Greece. At the large scale, we found a high genetic structure within the brooders (COI: 0.07 < F(ST) < 0.65) and no polymorphism at the nuclear locus i51. In contrast, the broadcasters displayed lower genetic structure (0 < F(ST) < 0.14) and were polymorphic at locus i51. At the regional scale, the multilocus analysis confirmed the contrasting genetic structure between species, with no structure in the broadcasters (global F(ST) < 0.001) and strong structure in the brooders (global F(ST) = 0.49), and revealed a higher genetic diversity in broadcasters. Our study showed that the lecithotrophic larval stage allows on average a 50-fold increase in migration rates, a 280-fold increase in effective size and a threefold to fourfold increase in genetic diversity. Our work, investigating complementary genetic markers on sympatric and syntopic taxa, highlights the strong impact of the larval phase on connectivity and genetic diversity.

Comparative population genomics in animals uncovers the determinants of genetic diversity.

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, domestication, 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. 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.

PCR survey of 50 introns in animals: cross-amplification of homologous EPIC loci in eight non-bilaterian, protostome and deuterostome phyla.

Exon Primed Intron Crossing (EPIC) markers provide molecular tools that are susceptible to be variable within species while remaining amplifiable by PCR using potentially universal primers. In this study we tested the possibility of obtaining PCR products from 50 EPIC markers on 23 species belonging to seven different phyla (Porifera, Cnidaria, Arthropoda, Nematoda, Mollusca, Annelida, Echinodermata) using 70 new primer pairs. A previous study had identified and tested those loci in a dozen species, including another phylum, Urochordata (Chenuil et al., 2010). Results were contrasted among species. The best results were achieved with the oyster (Mollusca) where 28 loci provided amplicons susceptible to contain an intron according to their size. This was however not the case with the other mollusk Crepidula fornicata, which seems to have undergone a reduction in intron number or intron size. In the Porifera, 13 loci appeared susceptible to contain an intron, a surprisingly high number for this phylum considering its phylogenetic distance with genomic data used to design the primers. For two cnidarian species, numerous loci (24) were obtained. Ecdysozoan phyla (arthropods and nematodes) proved less successful than others as expected considering reports of their rapid rate of genome evolution and the worst results were obtained for several arthropods. Some general patterns among phyla arose, and we discuss how the results of this EPIC survey may give new insights into genome evolution of the study species. This work confirms that this set of EPIC loci provides an easy-to-use toolbox to identify genetic markers potentially useful for population genetics, phylogeography or phylogenetic studies for a large panel of metazoan species. We then argue that obtaining diploid sequence genotypes for these loci became simple and affordable owing to Next-Generation Sequencing development. Species surveyed in this study belong to several genera (Acanthaster, Alvinocaris, Aplysina, Aurelia, Crepidula, Eunicella, Hediste, Hemimysis, Litoditis, Lophelia, Mesopodopsis, Mya, Ophiocten, Ophioderma, Ostrea, Pelagia, Platynereis, Rhizostoma, Rimicaris), two of them, belonging to the family Vesicomydae and Eunicidae, could not be determined at the genus level.

How to infer reliable diploid genotypes from NGS or traditional sequence data: from basic probability to experimental optimization.

The use of diploid sequence markers is still challenging despite the good quality of the information they provide. There is a common problem to all sequencing approaches [traditional cloning and sequencing of PCR amplicons as well as next-generation sequencing (NGS)]: when no variation is found within the sequences from a given individual, homozygozity can never be asserted with certainty. As a consequence, sequence data from diploid markers are mostly analysed at the population (not the individual level) particularly in animal studies. This study aims at contributing to solve this. Using the Bayes theorem and the binomial law, useful results are derived, among which: (i) the number of sequence reads per individual (or sequencing depth) which is required to ensure, at a given probability threshold, that some heterozygotes are not considered erroneously as homozygotes, as a function of the observed heterozygozity (H(o) ) of the locus in the population; (ii) a way of estimating H(o) from low coverage NGS data; (iii) a way of testing the null hypothesis that a genetic marker corresponds to a single and diploid locus, in the absence of data from controlled crosses; (iv) strategies for characterizing sequence genotypes in populations minimizing the average number of sequence reads per individual; (v) a rationale to decide which are the variations that one needs to consider along the sequence, as a function of the sequencing depth affordable, the level of polymorphism desired and the risk of sequencing error. For traditional sequencing technology, optimal strategies appear surprisingly different from the usual empirical ones. The average number of sequence reads required to obtain 99% of fully determined genotypes never exceeds six, this value corresponding to the worst situation when H(o) equals 0.6. This threshold value of H(o) is strikingly stable when the tolerated proportion of nonfully resolved genotypes varies in a reasonable range. These results do not rely on the Hardy-Weinberg equilibrium assumption or on diallelism of nucleotidic sites.

Did vicariance and adaptation drive cryptic speciation and evolution of brooding in Ophioderma longicauda (Echinodermata: Ophiuroidea), a common Atlanto-Mediterranean ophiuroid?

Over the last decade, cryptic speciation has been discovered in an increasing number of taxa. Species complexes are useful models for the understanding of speciation processes. Motivated by the discovery of brooding specimens in the common Atlanto-Mediterranean broadcast spawning brittle star, Ophioderma longicauda, a recent study revealed the occurrence of divergent mitochondrial lineages. We analysed 218 specimens from 23 locations spread over the geographic range of the species with partial Cytochrome c Oxidase subunit I (COI) sequences. A subset of this sample was also surveyed with the internal transcribed spacer of the ribosomal DNA cluster (nuclear ITS-1). Our study revealed six highly divergent mitochondrial lineages, and the ITS-1 data confirmed that they most likely represent a species complex. Geographic ranges, abundances and genetic structures are contrasted among the putative cryptic species. Lineages in which brooding specimens have been found form a monophyletic group and are restricted to the Eastern Mediterranean basin, an oligotrophic zone. A phylogeny-trait association analysis revealed a phylogenetic signal for low 'chlorophyll a' values (our proxy for oligotrophy). An ecological shift related to the hyper oligotrophy of the Eastern Mediterranean region is therefore likely to have played a role in the evolution of brooding. This study revealed that a complex mixture of vicariance, population expansion, adaptive divergence and possibly high local diversification rates resulting from brooding has shaped the evolution of this species complex. The dating analysis showed that these events probably occurred in the Pleistocene epoch.

Phylogeography of the red coral (Corallium rubrum): inferences on the evolutionary history of a temperate gorgonian.

The red coral Corallium rubrum (Cnidaria, Octocorallia) is an exploited, long-lived sessile species from the Mediterranean Sea and the adjacent coastline in the Atlantic Ocean. Surveys of genetic variation using microsatellites have shown that populations of C. rubrum are characterized by strong differentiation at the local scale but a study of the phylogeography of this species was still lacking. Here, we used seven polymorphic microsatellite loci, together with sequence data from an intron of the elongation factor 1 (EF1) gene, to investigate the genetic structure of C. rubrum across its geographical range in the western Mediterranean Sea and in the Adriatic Sea. The EF1 sequences were also used to analyse the consequences of demographic fluctuations linked with past environmental change. Clustering analysis with microsatellite loci highlighted three to seven genetic groups with the distinction of North African and Adriatic populations; this distinction appeared significant with AMOVA and differentiation tests. Microsatellite and EF1 data extended the isolation by distance pattern previously observed for this species at the western Mediterranean scale. EF1 sequences confirmed the genetic differentiation observed between most samples with microsatellites. A statistical parsimony network of EF1 haplotypes provided no evidence of high sequence divergence among regions, suggesting no long-term isolation. Selective neutrality tests on microsatellites and EF1 were not significant but should be interpreted with caution in the case of EF1 because of the low sample sizes for this locus. Our results suggest that recent Quaternary environmental fluctuations had a limited impact on the genetic structure of C. rubrum.

Defining reproductively isolated units in a cryptic and syntopic species complex using mitochondrial and nuclear markers: the brooding brittle star, Amphipholis squamata (Ophiuroidea).

At a time when biodiversity is threatened, we are still discovering new species, and particularly in the marine realm. Delimiting species boundaries is the first step to get a precise idea of diversity. For sympatric species which are morphologically undistinguishable, using a combination of independent molecular markers is a necessary step to define separate species. Amphipholis squamata, a cosmopolitan brittle star, includes several very divergent mitochondrial lineages. These lineages appear totally intermixed in the field and studies on morphology and colour polymorphism failed to find any diagnostic character. Therefore, these mitochondrial lineages may be totally interbreeding presently. To test this hypothesis, we characterized the genetic structure of the complex in the French Mediterranean coast using sequences of mitochondrial DNA (16S) and for the first time, several nuclear DNA markers (introns and microsatellites). The data revealed six phylogenetic lineages corresponding to at least four biological species. These sibling species seem to live in syntopy. However, they seem to display contrasted levels of genetic diversity, suggesting they have distinct demographic histories and/or life-history traits. Genetic differentiation and isolation-by-distance within the French Mediterranean coasts are revealed in three lineages, as expected for a species without a free larval phase. Finally, although recombinant nuclear genotypes are common within mitochondrial lineages, the data set displays a total lack of heterozygotes, suggesting a very high selfing rate, a feature likely to have favoured the formation of the species complex.

Autosomal differences between males and females in hybrid zones: a first report from Barbus barbus and Barbus meridionalis (Cyprinidae).

Narrow hybrid zones are generally subjected to the action of two forces: dispersal, which tends to homogenise the hybridising taxa, and selection against hybrids, which, in contrast, produces steep clines of introgression for diagnostic markers. Although differences between sexes in dispersal abilities or in susceptibility to hybrid counterselection are common in hybrid zones, autosomal genetic differences between males and females have never been reported to our knowledge. Barbus barbus and Barbus meridionalis (Cyprinidae) form a hybrid zone along the Lergue river. By carrying out a genetic analysis of males and females in six samples from two central stations of the hybrid zone using codominant markers (six allozymes and four microsatellite loci), we revealed significant multilocus and monolocus differences between the sexes. This could reflect a genetic difference among sexes within a same cohort, caused either by a survival (or fertility) differential among sexes or by a sex-specific pattern of dispersal. Alternatively, this may be due to genetic differentiation between cohorts, since male and female barbel exhibit different maturation, growth and survival patterns leading to different age distributions among sexes, and particularly among reproducers.

A test of the hypothesis of an autopolyploid vs. allopolyploid origin for a tetraploid lineage: application to the genus barbus (Cyprinidae)

A new method is described for determination of the origin of polyploid lineages. It tests the hypothesis that a tetraploid lineage originated via autopolyploidization vs. allopolyploidization. The method is based on the hypothesis that, in the case of autopolyploidy, any genetic marker in the first tetraploid ancestor is represented by two copies (one for each homoeologous chromosome of the haploid complement), whereas in allopolyploidy some markers absent from one of the hybridizing species will display one copy at most. The model requires knowledge of the phylogeny (topology and branch lengths) of a sample of species descending from the same tetraploidization event, together with the number of homoeologous copies present in each species for a set of neutral markers. The likelihood of a given proportion of the markers being present in both homoeologous chromosome pairs of the ancestral tetraploid is expressed as a function of the deletion rate of a marker. In the case of an autopolyploid origin, this proportion equals one. A likelihood-ratio test was carried out to test this hypothesis. The method was used to examine five microsatellite loci in eight species of Barbus (sensu lato). Assuming the validity of the hypotheses on phylogenetic relationships and evolutionary rates, the test rejects the possibility that European tetraploid barbs originated through autopolyploidy. This is the first test that can reject autopolyploidy, and it would appear particularly useful for phylogenetic studies in taxa where hybridization is known and where, consequently, undetected reticulate evolution may impair phylogenetic reconstruction.

Evolution of the large-subunit ribosomal RNA binding site for protein L23/25.

The region of the large-subunit rRNA encompassing the D7 divergent domain is organized within eukaryotes in a patchwork of short conservative secondary-structure features interspersed with more rapidly evolving sequences. It contains the attachment site of protein L25 (E. coli L23), which binds rRNA in the first stages of ribosome assembly, suggesting a crucial importance of this region in ribosome elaboration and functioning. A better understanding of its roles requires a good knowledge of its mode of structural variation during the course of evolution. With this aim, we sequenced the D7 region for 24 new invertebrate species belonging to annelids, molluscs, arthropods, and eight other deep-branching invertebrate phyla. Their comparison allowed us to propose refinements in previous eukaryotic folding models. A detailed analysis of the pattern of variation at each position both within the D7 region and along the L23/25 sequence by reference to previous heterologous binding experiments gives new insight into the rRNA-protein contacts. We identified in the D7 region and L23/25, respectively, six and five positions presenting a pattern of variation compatible with experimental results, three of which show coincident variations which support their possible involvement in the rRNA-L23/25 binding.

Does polyploidy lead to fewer and shorter microsatellites in Barbus (Teleostei:Cyprinidae)?

Screening of a hybrid Barbus barbus-B. meridionalis genome was performed for CA, GA, TAT, TCT, TAG, TGT, TATT, TACT, ATCT motifs, and simultaneously on another fish species, tilapia S. melanotheron. Sequences of positive clones were obtained for Barbus and revealed that repetitive structure significantly depends on the motif: most TAT and TATT repeats contain small numbers of repeats, and these repeats are highly heterogeneous, whereas other motifs (we mainly obtained CA and GATA repeats) form longer and much more homogeneous arrays. Polymorphism data from five loci in two different species of barbel show that perfectly repetitive loci are much more variable than imperfect loci (TAT and TATT). We compared the frequency of positive clones for different repeat motifs between barbel and tilapia. For dinucleotide repeats (CA and GA), the comparison was extended to additional fish species, trout and sea bass, which were screened in nearly identical conditions for these motifs. The most salient feature of these comparisons reveals that arrays of dinucleotide motifs are significantly under-represented and shorter in Barbus than in other fish species. We propose an explanation that can account for most features of microsatellites characterizing the genome of barbel. A bias toward deletion affecting slipped-strand mispairing events would lead to shortening and loss of microsatellite loci. Such a bias would represent an efficient way of eliminating useless DNA from polyploidized species with an excessive amount of DNA.

Molecular phylogenetic study of a myrmecophyte symbiosis: did Leonardoxa/ ant associations diversify via cospeciation?

The Leonardoxa africana (Leguminosae: Caesalpinioideae) complex is a group of four closely related taxa (L1 to L4) exhibiting various grades of specificity and specialization in mutualistic associations with ants. Each of the two most specialized species, Leonardoxa taxon 3 (L3) and L. africana sensu stricto (L4), interacts with a specific species of formicine ant, respectively Aphomomyrmex after and Petalomyrmex phylax, which nests in specialized swollen twigs. These two monotypic genera are the sole African members of the tribe Myrmelachistini, and their occurrence in closely related plants suggested thehypothesis that the two associations L4/Petalomyrmex and L3/Aphomomyrmex are derived by cospeciation from an ancestral association. Phylogenies based on DNA sequences were reconstructed for the ants and compared with phylogenies available for the plants in order to test for this hypothesis of cospeciation. The resulting topologies suggest either that the association with myrmelachistine ants arose several times or that a plant species (L2) and an ant population split off from an ancestral association. Furthermore, dates of speciation events appear to differ between ants and corresponding plants. An estimate of at least 4 million years was obtained for the separation of Aphomomyrmex and Petalomyrmex, whereas biological, biogeographic, and molecular-genetic data suggest a much more recent divergence for the plants. Thus, we reject the hypothesis of cospeciation and conclude that Aphomomyrmex and Petalomyrmex independently colonized different taxa of Leonardoxa. This striking instance of parallel evolution supports the notion that specific ant-plant associations originated by ecological fitting of preadapted partners. We discuss alternative evolutionary scenarios that are consistent with molecular data.