Fish gut and skin microbiota dysbiosis induced by exposure to commercial sunscreen formulations.

UV filters (organic or mineral) present in sunscreen products are emerging contaminants of coastal aquatic environments. There is an urgent need to understand marine organisms responses to these compounds. In this study, we investigated the effect of exposure to dilutions of commercial sunscreen formulations on bacterial communities of mullet (Chelon sp.). The gut and skin mucus microbial communities were characterized using a metabarcoding approach targeting the 16S rRNA gene. Our results revealed that mullets had its own bacterial communities that differ from their surrounding habitats and specific to tissue. The dilutions of commercial sunscreens modified the relative abundance of Actinobacteroita, Bacteriodota and Proteobacteria for both gut and skin microbiota. They also allowed to bacteria affiliated to Mycobacterium, Nocardia and Tenacibaculum genera, known to house pathogenic species, to colonize the epithelium which may have implications for fish host health.

LAMP-based molecular sexing in a gonochoric marine bivalve (Macoma balthica rubra) with divergent sex-specific mitochondrial genomes.

Taking advantage of the unique system of doubly uniparental inheritance (DUI) of mitochondria, we developed a reliable molecular method to sex individuals of the marine bivalve Macoma balthica rubra. In species with DUI (~100 known bivalves), both sexes transmit their mitochondria: males bear both a male- and female-type mitogenome, while females bear only the female type. Male and female mitotypes are sufficiently divergent to reliably PCR-amplify them specifically. Loop-mediated isothermal amplification (LAMP) is a precise, economical and portable alternative to PCR for molecular sexing and we demonstrate its application in this context. We used 154 individuals sampled along the Atlantic coast of France and sexed microscopically by gonad examination to test for the congruence among gamete type, PCR sexing and LAMP sexing. We show an exact match among the sexing results from these three methods using the male and female mt-cox1 genes. DUI can be disrupted in inter-specific hybrids, causing unexpected distribution of mitogenomes, such as homoplasmic males or heteroplasmic females. To our knowledge, DUI disruption at the intra-specific scale has never been tested. We applied our sexing protocol to control for unexpected heteroplasmy caused by hybridization between divergent genetic lineages and found no evidence of disruption in the mode of mitochondrial inheritance in M. balthica rubra. We propose LAMP as a useful tool to accelerate eco-evolutionary studies of DUI. It offers the opportunity to investigate the potential role of, previously unaccounted-for, sex-specific patterns such as sexual selection or sex-specific dispersal bias in the evolution of free-spawning benthic species.

Growth in two deep-sea associates: the octocoral Pseudogorgia bellona and the euryalid snake star Asteroschema ajax.

The deep-sea octocoral Pseudochrysogorgia bellona was recently described from specimens sampled on the Chesterfield Plateau, off New Caledonia. It is morphologically and genetically similar to the con-familial Metallogorgia melanotrichos, which is known to closely associate with a species of brittle star, Ophiocreas oedipus. These latter two species have never been observed separately and are thought to grow synchronously. The morphological similarity between M. melanotrichos and P. bellona makes the latter another possible host for ophiuroids. However, no brittle star was associated with P. bellona specimens from the type collection. In 2017, 130 P. bellona colonies were sampled near the type locality, and 98% were associated with Asteroschema ajax, a species closely related to O. oedipus. Mitochondrial DNA analysis confirmed the morphological identifications of both P. bellona and A. ajax. Uni- and multivariate statistical analyses were used to characterize the morphological space of both species to test if larger ophiuroids are associated with larger corals. Two variables were measured to estimate the size of the coral (total height and diameter of the skeletal axis at its base) and 9 variables were used to characterize the brittle star (disc and arm morphology). Morphological variables representing the size for both species were significantly correlated (Spearman rank correlation coefficient: 50%, p < 0.001), suggesting that larger ophiuroids indeed associate with larger corals. This is one of the rare studies that allowed comparison of growth in associated deep-sea invertebrates.

The longest mitochondrial protein in metazoans is encoded by the male-transmitted mitogenome of the bivalve Scrobicularia plana.

Cytochrome c oxidase subunit II (COX2) is one of the three mitochondrially encoded proteins of the complex IV of the respiratory chain that catalyses the reduction of oxygen to water. The cox2 gene spans about 690 base pairs in most animal species and produces a protein composed of approximately 230 amino acids. We discovered an extreme departure from this pattern in the male-transmitted mitogenome of the bivalve Scrobicularia plana with doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA), which possesses an important in-frame insertion of approximately 4.8 kb in its cox2 gene. This feature-an enlarged male cox2 gene-is found in many species with DUI; the COX2 protein can be up to 420 amino acids long. Through RT-PCRs, immunoassays and comparative genetics, the evolution and functionality of this insertion in S. plana were characterized. The in-frame insertion is conserved among individuals from different populations and bears the signature of purifying selection seemingly indicating maintenance of functionality. Its transcription and translation were confirmed: this gene produces a polypeptide of 1892 amino acids, making it the largest metazoan COX2 protein known to date. We hypothesize that these extreme modifications in the COX2 protein affect the metabolism of mitochondria containing the male-transmitted mtDNA in Scrobicularia plana.

Complex mitogenomic rearrangements within the Pectinidae (Mollusca: Bivalvia).

BACKGROUND: Scallops (Bivalvia: Pectinidae) present extraordinary variance in both mitochondrial genome size, structure and content, even when compared to the extreme diversity documented within Mollusca and Bivalvia. In pectinids, mitogenome rearrangements involve protein coding and rRNA genes along with tRNAs, and different genome organization patterns can be observed even at the level of Tribes. Existing pectinid phylogenies fail to resolve some relationships in the family, Chlamydinae being an especially problematic group. RESULTS: In our study, we sequenced, annotated and characterized the mitochondrial genome of a member of Chlamydinae, Mimachlamys varia-a species of commercial interest and an effective bioindicator-revealing yet another novel gene arrangement in the Pectinidae. The phylogeny based on all mitochondrial protein coding and rRNA genes suggests the paraphyly of the Mimachlamys genus, further commending the taxonomic revision of the classification within the Chlamydinae subfamily. At the scale of the Pectinidae, we found that 15 sequence blocks are involved in mitogenome rearrangements, which behave as separate units. CONCLUSIONS: Our study reveals incongruities between phylogenies based on mitochondrial protein-coding versus rRNA genes within the Pectinidae, suggesting that locus sampling affects phylogenetic inference at the scale of the family. We also conclude that the available taxon sampling does not allow for understanding of the mechanisms responsible for the high variability of mitogenome architecture observed in the Pectinidae, and that unraveling these processes will require denser taxon sampling.

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds.

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Molecular genetic diversity of seaweeds morphologically related to Ulva rigida at three sites along the French Atlantic coast.

Foliose species of the genus Ulva are notoriously difficult to identify due to their variable morphological characteristics and high phenotypic plasticity. We reassessed the taxonomic status of several distromatic foliose Ulva spp., morphologically related to Ulva rigida, using DNA barcoding with the chloroplastic tufA and rbcL (for a subset of taxa) genes for 339 selected attached Ulva specimens collected from three intertidal rocky sites. Two of the collection sites were in Brittany and one site was in Vendée, along the Atlantic coast of France. Molecular analyses included several museum specimens and the holotype of Ulva armoricana Dion, Reviers & Coat. We identified five different tufA haplotypes using a combination of phylogenetic analysis, with the support of several recently sequenced holotypes and lectotypes, and a species delimitation method based on hierarchical clustering. Four haplotypes were supported by validly named species: Ulva australis Areschoug, Ulva fenestrata Postels & Ruprecht, Ulva lacinulata (Kützing) Wittrock and U. rigida C. Agardh. The later was additionally investigated using rbcL. The fifth haplotype represented exact sequence matches to an unnamed species from European Atlantic coasts. Our results support: (1) the synonymy of both U. rigida sensu Bliding non C. Agardh and U. armoricana with U. lacinulata. This finding is based on current genetic analysis of tufA from the U. armoricana holotype and recent molecular characterization of the lectotype of U. laetevirens, which is synonymous to U. australis, (2) the presence of U. australis as a misidentified introduced species in Brittany, and (3) the presence of U. fenestrata and U. rigida in southern Brittany. The taxonomic history of each species is discussed, highlighting issues within distromatic foliose taxa of the genus Ulva and the need to genetically characterize all its available type specimens.

Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria.

In animals, strictly maternal inheritance (SMI) of mitochondria is the rule, but one exception (doubly uniparental inheritance or DUI), marked by the transmission of sex-specific mitogenomes, has been reported in bivalves. Associated with DUI is a frequent modification of the mitochondrial cox2 gene, as well as additional sex-specific mitochondrial genes not involved in oxidative phosphorylation. With the exception of freshwater mussels (for 3 families of the order Unionida), these DUI-associated features have only been shown in few species [within Mytilidae (order Mytilida) and Veneridae (order Venerida)] because of the few complete sex-specific mitogenomes published for these orders. Here, we present the complete sex-specific mtDNAs of two recently-discovered DUI species in two families of the order Venerida, Scrobicularia plana (Semelidae) and Limecola balthica (Tellinidae). These species display the largest differences in genome size between sex-specific mitotypes in DUI species (>10 kb), as well as the highest mtDNA divergences (sometimes reaching >50%). An important in-frame insertion (>3.5 kb) in the male cox2 gene is partly responsible for the differences in genome size. The S. plana cox2 gene is the largest reported so far in the Kingdom Animalia. The mitogenomes may be carrying sex-specific genes, indicating that general mitochondrial features are shared among DUI species.

Data, time and money: evaluating the best compromise for inferring molecular phylogenies of non-model animal taxa.

For over a decade now, High Throughput sequencing (HTS) approaches have revolutionized phylogenetics, both in terms of data production and methodology. While transcriptomes and (reduced) genomes are increasingly used, generating and analyzing HTS datasets remain expensive, time consuming and complex for most non-model taxa. Indeed, a literature survey revealed that 74% of the molecular phylogenetics trees published in 2018 are based on data obtained through Sanger sequencing. In this context, our goal was to identify the strategy that would represent the best compromise among costs, time and robustness of the resulting tree. We sequenced and assembled 32 transcriptomes of the marine mollusk family Turridae, considered as a typical non-model animal taxon. From these data, we extracted the loci most commonly used in gastropod phylogenies (cox1, 12S, 16S, 28S, h3 and 18S), full mitogenomes, and a reduced nuclear transcriptome representation. With each dataset, we reconstructed phylogenies and compared their robustness and accuracy. We discuss the impact of missing data and the use of statistical tests, tree metrics, and supertree and supermatrix methods to further improve phylogenetic data acquisition pipelines. We evaluated the overall costs (time and money) in order to identify the best compromise for phylogenetic data sampling in non-model animal taxa. Although sequencing full mitogenomes seems to constitute the best compromise both in terms of costs and node support, they are known to induce biases in phylogenetic reconstructions. Rather, we recommend to systematically include loci commonly used for phylogenetics and taxonomy (i.e. DNA barcodes, rRNA genes, full mitogenomes, etc.) among the other loci when designing baits for capture.

Advances on the phylogenetic placement of the enigmatic octocoral Dendrobrachia Brook 1889.

The monogeneric family Dendrobrachiidae has been a taxonomic curiosity since its original description in 1889. Using one nuclear (18S) and two mitochondrial (mtMutS and cox1) genes, the phylogenetic placement of Dendrobrachiidae within the Octocorallia was investigated based on recently-collected specimens and museum collections. In particular, the relationship between Dendrobrachia and its suspected close allies from the Chrysogorgiidae and Ifalukellidae was examined. Phylogenetic inferences based on nuclear 18S sequences were inconclusive. On the basis of mitochondrial cox1 and mtMutS, the hypothesis that Dendrobrachia is closely related to the Chrysogorgiidae or the Ifalukellidae could be rejected with confidence. Dendrobrachia appeared distantly related to all octocoral groups considered here. However, the phylogenetic placement of Dendrobrachia could not be further clearly resolved.

Ocean warming combined with lower omega-3 nutritional availability impairs the cardio-respiratory function of a marine fish.

Highly unsaturated fatty acids of the omega-3 series (HUFA) are major constituents of cell membranes, yet are poorly synthesised de novo by consumers. Their production, mainly supported by aquatic microalgae, has been decreasing with global change. The consequences of such reductions may be profound for ectotherm consumers, as temperature tightly regulates the HUFA content in cell membranes, maintaining their functionality. Integrating individual, tissue and molecular approaches, we examined the consequences of the combined effects of temperature and HUFA depletion on the key cardio-respiratory functions of the golden grey mullet, an ectotherm grazer of high ecological importance. For 4 months, fish were exposed to two contrasting HUFA diets [4.8% eicosapentaenoic acid (EPA)+docosahexaenoic acid (DHA) on dry matter (DM) versus 0.2% EPA+DHA on DM] at 12 and 20°C. Ventricular force development coupled with gene expression profiles measured on cardiac muscle suggest that combining HUFA depletion with warmer temperatures leads to: (1) a proliferation of sarcolemmal and sarcoplasmic reticulum Ca2+ channels and (2) a higher force-generating ability by increasing extracellular Ca2+ influx via sarcolemmal channels when the heart has to sustain excessive effort due to stress and/or exercise. At the individual scale, these responses were associated with a greater aerobic scope, maximum metabolic rate and net cost of locomotion, suggesting the higher energy cost of this strategy. This impaired cardiac performance could have wider consequences for other physiological performance such as growth, reproduction or migration, all of which greatly depend on heart function.

Evidence for a duplicated mitochondrial region in Audubon's shearwater based on MinION sequencing.

Mitochondrial genetic markers have been extensively used to study the phylogenetics and phylogeography of many birds, including seabirds of the order Procellariiformes. Evidence suggests that part of the mitochondrial genome of Procellariiformes, especially albatrosses, is duplicated, but no DNA fragment covering the entire duplication has been sequenced. We sequenced the complete mitochondrial genome of a non-albatross species of Procellariiformes, Puffinus lherminieri (Audubon's shearwater) using the long-read MinION (ONT) technology. Two mitogenomes were assembled from the same individual, differing by 52 SNPs and in length. The shorter was 19 kb long while the longer was 21 kb, due to the presence of two identical copies of nad6, three tRNA, and two dissimilar copies of the control region (CR). Contrary to albatrosses, cob was not duplicated. We further detected a complex repeated region of undetermined length between the CR and 12S. Long-read sequencing suggests heteroplasmy and a novel arrangement within the duplicated region, indicating a complex evolution of the mitogenome in Procellariiformes.

Nomenclatural note on the homonymy between the octocoral genus Dendrobrachia Brook 1889 and the entoproct phylum Dendrobrachia Xian-Guan, Bergström, Xio-Ya and Jie 2006.

Dendrobrachia was proposed by Xian-Guan et al. (2006) as the name of a new phylum of entoproct-like animals from the Lower Cambrian (Yunnan Province, China). This name was however preoccupied by the octocoral Dendrobrachia Brook 1889, which position within the subclass Octocorallia is currently uncertain (Opresko and Bayer 1991; Berntson et al. 1999, 2001). Based on the principles of homonymy and priority (Articles 52 and 23, respectively; ICZN 1999), a substitute name for the entoproct Dendrobrachia is necessary, and I propose Dendrobrachion (based on the Greek brachion, for arm).

Short-Term and Long-Term Biological Effects of Chronic Chemical Contamination on Natural Populations of a Marine Bivalve.

Understanding the effects of chronic chemical contamination on natural populations of marine organisms is complex due to the combined effects of different types of pollutants and environmental parameters that can modulate the physiological responses to stress. Here, we present the effects of a chronic contamination in a marine bivalve by combining multiple approaches that provide information on individual and population health. We sampled variegated scallops (Mimachlamys varia) at sites characterized by different contaminants and contamination levels to study the short and long-term (intergenerational) responses of this species to physiological stress. We used biomarkers (SOD, MDA, GST, laccase, citrate synthase and phosphatases) as indicators of oxidative stress, immune system alteration, mitochondrial respiration and general metabolism, and measured population genetic diversity at each site. In parallel, concentration of 14 trace metals and 45 organic contaminants (PAHs, PCBs, pesticides) in tissues were measured. Scallops were collected outside and during their reproductive season to investigate temporal variability in contaminant and biomarker levels. Our analyses revealed that the levels of two biomarkers (Laccase-type phenoloxidase and malondialdehyde) were significantly correlated with Cd concentration. Additionally, we observed significant seasonal differences for four of the five biomarkers, which is likely due to the scallop reproductive status at time of sampling. As a source of concern, a location that was identified as a reference site on the basis of inorganic contaminant levels presented the same level of some persistent organic pollutants (DDT and its metabolites) than more impacted sites. Finally, potential long-term effects of heavy metal contamination were observed for variegated scallops as genetic diversity was depressed in the most polluted sites.

Species are hypotheses: avoid connectivity assessments based on pillars of sand.

Connectivity among populations determines the dynamics and evolution of populations, and its assessment is essential in ecology in general and in conservation biology in particular. The robust basis of any ecological study is the accurate delimitation of evolutionary units, such as populations, metapopulations and species. Yet a disconnect still persists between the work of taxonomists describing species as working hypotheses and the use of species delimitation by molecular ecologists interested in describing patterns of gene flow. This problem is particularly acute in the marine environment where the inventory of biodiversity is relatively delayed, while for the past two decades, molecular studies have shown a high prevalence of cryptic species. In this study, we illustrate, based on marine case studies, how the failure to recognize boundaries of evolutionary-relevant unit leads to heavily biased estimates of connectivity. We review the conceptual framework within which species delimitation can be formalized as falsifiable hypotheses and show how connectivity studies can feed integrative taxonomic work and vice versa. Finally, we suggest strategies for spatial, temporal and phylogenetic sampling to reduce the probability of inadequately delimiting evolutionary units when engaging in connectivity studies.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities.

BACKGROUND: Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform. RESULTS: The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations. CONCLUSIONS: Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.

Applicability of RAD-tag genotyping for interfamilial comparisons: empirical data from two cetaceans.

Restriction-site-associated DNA tag (RAD-tag) sequencing has become a popular approach to generate thousands of SNPs used to address diverse questions in population genomics. Comparatively, the suitability of RAD-tag genotyping to address evolutionary questions across divergent species has been the subject of only a few recent studies. Here, we evaluate the applicability of this approach to conduct genome-wide scans for polymorphisms across two cetacean species belonging to distinct families: the short-beaked common dolphin (Delphinus delphis; n = 5 individuals) and the harbour porpoise (Phocoena phocoena; n = 1 individual). Additionally, we explore the effects of varying two parameters in the Stacks analysis pipeline on the number of loci and level of divergence obtained. We observed a 34% drop in the total number of loci that were present in all individuals when analysing individuals from the distinct families compared with analyses restricted to intraspecific comparisons (i.e. within D. delphis). Despite relatively stringent quality filters, 3595 polymorphic loci were retrieved from our interfamilial comparison. Cetaceans have undergone rapid diversification, and the estimated divergence time between the two families is relatively recent (14-19 Ma). Thus, our results showed that, for this level of divergence, a large number of orthologous loci can still be genotyped using this approach, which is on par with two recent in silico studies. Our findings constitute one of the first empirical investigations using RAD-tag sequencing at this level of divergence and highlights the great potential of this approach in comparative studies and to address evolutionary questions.

marmap: A package for importing, plotting and analyzing bathymetric and topographic data in R.

In this communication we introduce marmap, a package designed for downloading, plotting and manipulating bathymetric and topographic data in R. marmap can query the ETOPO1 bathymetry and topography database hosted by the NOAA, use simple latitude-longitude-depth data in ascii format, and take advantage of the advanced plotting tools available in R to build publication-quality bathymetric maps. Functions to query data (bathymetry, sampling information…) are available interactively by clicking on marmap maps. Bathymetric and topographic data can also be used to calculate projected surface areas within specified depth/altitude intervals, and constrain the calculation of realistic shortest path distances. Such information can be used in molecular ecology, for example, to evaluate genetic isolation by distance in a spatially-explicit framework.

Deep-sea origin and in-situ diversification of chrysogorgiid octocorals.

The diversity, ubiquity and prevalence in deep waters of the octocoral family Chrysogorgiidae Verrill, 1883 make it noteworthy as a model system to study radiation and diversification in the deep sea. Here we provide the first comprehensive phylogenetic analysis of the Chrysogorgiidae, and compare phylogeny and depth distribution. Phylogenetic relationships among 10 of 14 currently-described Chrysogorgiidae genera were inferred based on mitochondrial (mtMutS, cox1) and nuclear (18S) markers. Bathymetric distribution was estimated from multiple sources, including museum records, a literature review, and our own sampling records (985 stations, 2345 specimens). Genetic analyses suggest that the Chrysogorgiidae as currently described is a polyphyletic family. Shallow-water genera, and two of eight deep-water genera, appear more closely related to other octocoral families than to the remainder of the monophyletic, deep-water chrysogorgiid genera. Monophyletic chrysogorgiids are composed of strictly (Iridogorgia Verrill, 1883, Metallogorgia Versluys, 1902, Radicipes Stearns, 1883, Pseudochrysogorgia Pante & France, 2010) and predominantly (Chrysogorgia Duchassaing & Michelotti, 1864) deep-sea genera that diversified in situ. This group is sister to gold corals (Primnoidae Milne Edwards, 1857) and deep-sea bamboo corals (Keratoisidinae Gray, 1870), whose diversity also peaks in the deep sea. Nine species of Chrysogorgia that were described from depths shallower than 200 m, and mtMutS haplotypes sequenced from specimens sampled as shallow as 101 m, suggest a shallow-water emergence of some Chrysogorgia species.

SNP detection from de novo transcriptome sequencing in the bivalve Macoma balthica: marker development for evolutionary studies.

Hybrid zones are noteworthy systems for the study of environmental adaptation to fast-changing environments, as they constitute reservoirs of polymorphism and are key to the maintenance of biodiversity. They can move in relation to climate fluctuations, as temperature can affect both selection and migration, or remain trapped by environmental and physical barriers. There is therefore a very strong incentive to study the dynamics of hybrid zones subjected to climate variations. The infaunal bivalve Macoma balthica emerges as a noteworthy model species, as divergent lineages hybridize, and its native NE Atlantic range is currently contracting to the North. To investigate the dynamics and functioning of hybrid zones in M. balthica, we developed new molecular markers by sequencing the collective transcriptome of 30 individuals. Ten individuals were pooled for each of the three populations sampled at the margins of two hybrid zones. A single 454 run generated 277 Mb from which 17K SNPs were detected. SNP density averaged 1 polymorphic site every 14 to 19 bases, for mitochondrial and nuclear loci, respectively. An [Formula: see text] scan detected high genetic divergence among several hundred SNPs, some of them involved in energetic metabolism, cellular respiration and physiological stress. The high population differentiation, recorded for nuclear-encoded ATP synthase and NADH dehydrogenase as well as most mitochondrial loci, suggests cytonuclear genetic incompatibilities. Results from this study will help pave the way to a high-resolution study of hybrid zone dynamics in M. balthica, and the relative importance of endogenous and exogenous barriers to gene flow in this system.