Selection-driven adaptation to the extreme Antarctic environment in the Emperor penguin.

The eco-evolutionary history of penguins is characterised by shifting from temperate to cold environments. Breeding in Antarctica, the Emperor penguin appears as an extreme outcome of this process, with unique features related to insulation, heat production and energy management. However, whether this species actually diverged from a less cold-adapted ancestor, more ecologically similar to its sister species, the King penguin, is still an open question. As the Antarctic colonisation likely resulted in vast changes in selective pressure experienced by the Emperor penguin, the relative quantification of the genomic signatures of selection, unique to each sister species, could answer this question. Applying phylogeny-based selection tests on 7651 orthologous genes, we identified a more pervasive selection shift in the Emperor penguin than in the King penguin, supporting the hypothesis that its extreme cold adaptation is a derived state. Furthermore, among candidate genes under selection, four (TRPM8, LEPR, CRB1, and SFI1) were identified before in other cold-adapted homeotherms, like the woolly Mammoth, while other 161 genes can be assigned to biological functions relevant to cold adaptation identified in previous studies. Location and structural effects of TRPM8 substitutions in Emperor and King penguin lineages support their functional role with putative divergent effects on thermal adaptation. We conclude that extreme cold adaptation in the Emperor penguin largely involved unique genetic options which, however, affect metabolic and physiological traits common to other cold-adapted homeotherms.

Low impact of different SNP panels from two building-loci pipelines on RAD-Seq population genomic metrics: case study on five diverse aquatic species.

BACKGROUND: The irruption of Next-generation sequencing (NGS) and restriction site-associated DNA sequencing (RAD-seq) in the last decade has led to the identification of thousands of molecular markers and their genotyping for refined genomic screening. This approach has been especially useful for non-model organisms with limited genomic resources. Many building-loci pipelines have been developed to obtain robust single nucleotide polymorphism (SNPs) genotyping datasets using a de novo RAD-seq approach, i.e. without reference genomes. Here, the performances of two building-loci pipelines, STACKS 2 and Meyer's 2b-RAD v2.1 pipeline, were compared using a diverse set of aquatic species representing different genomic and/or population structure scenarios. Two bivalve species (Manila clam and common edible cockle) and three fish species (brown trout, silver catfish and small-spotted catshark) were studied. Four SNP panels were evaluated in each species to test both different building-loci pipelines and criteria for SNP selection. Furthermore, for Manila clam and brown trout, a reference genome approach was used as control. RESULTS: Despite different outcomes were observed between pipelines and species with the diverse SNP calling and filtering steps tested, no remarkable differences were found on genetic diversity and differentiation within species with the SNP panels obtained with a de novo approach. The main differences were found in brown trout between the de novo and reference genome approaches. Genotyped vs missing data mismatches were the main genotyping difference detected between the two building-loci pipelines or between the de novo and reference genome comparisons. CONCLUSIONS: Tested building-loci pipelines for selection of SNP panels seem to have low influence on population genetics inference across the diverse case-study scenarios here studied. However, preliminary trials with different bioinformatic pipelines are suggested to evaluate their influence on population parameters according with the specific goals of each study.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea.

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species' evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Contrasting life-history traits of two sympatric smooth-hound species: implication for vulnerability.

In this study, life-history traits (maximum and average size, size at maturity and fecundity) of two congeneric smooth-hounds, Mustelus mustelus and Mustelus punctulatus, which share a geographical distribution and experience a similar fishing exploitation, were estimated and compared between species. The results indicated a lower maximum and average size, a lower size at maturity and a higher fecundity in M. punctulatus compared with those in M. mustelus. Considering that these two species co-occur in the same areas and are caught by the same fishing gears, the results indicate a higher vulnerability to exploitation of M. mustelus compared with that of M. punctulatus.

Draft genome assembly and transcriptome data of the icefish Chionodraco myersi reveal the key role of mitochondria for a life without hemoglobin at subzero temperatures.

Antarctic fish belonging to Notothenioidei represent an extraordinary example of radiation in the cold. In addition to the absence of hemoglobin, icefish show a number of other striking peculiarities including large-diameter blood vessels, high vascular densities, mitochondria-rich muscle cells, and unusual mitochondrial architecture. In order to investigate the bases of icefish adaptation to the extreme Southern Ocean conditions we sequenced the complete genome of the icefish Chionodraco myersi. Comparative analyses of the icefish genome with those of other teleost species, including two additional white-blooded and five red-blooded notothenioids, provided a new perspective on the evolutionary loss of globin genes. Muscle transcriptome comparative analyses against red-blooded notothenioids as well as temperate fish revealed the peculiar regulation of genes involved in mitochondrial function in icefish. Gene duplication and promoter sequence divergence were identified as genome-wide patterns that likely contributed to the broad transcriptional program underlying the unique features of icefish mitochondria.

Connectivity Among Populations of the Top Shell Gibbula divaricata in the Adriatic Sea.

Genetic connectivity studies are essential to understand species diversity and genetic structure and to assess the role of potential factors affecting connectivity, thus enabling sound management and conservation strategies. Here, we analyzed the patterns of genetic variability in the marine snail Gibbula divaricata from five coastal locations in the central-south Adriatic Sea (central Mediterranean) and one in the adjacent northern Ionian Sea, using 21 described polymorphic microsatellite loci. Observed and expected heterozygosity varied from 0.582 to 0.635 and 0.684 to 0.780, respectively. AMOVA analyses showed that 97% of genetic variation was observed within populations. Nevertheless, significant, although small, genetic differentiation was found among nearly all of the pairwise F ST comparisons. Over a general pattern of panmixia, three groups of populations were identified: eastern Adriatic populations, western Adriatic populations, and a third group represented by the single northern Ionian Sea population. Nonetheless, migration and gene flow were significant between these groups. Gibbula divaricata is thought to have a limited dispersal capacity related to its lecithotrophic trochophore larval stage. Our results indicated high levels of self-recruitment and gene flow that is mainly driven through coastline dispersion, with populations separated by the lack of suitable habitats or deep waters. This stepping-stone mode of dispersion together with the high levels of self-recruitment could lead to higher levels of population structuring and differentiation along the Adriatic Sea. Large effective population sizes and episodic events of long-distance dispersal might be responsible for the weak differentiation observed in the analyzed populations. In summary, the circulation system operating in this region creates natural barriers for dispersion that, together with life-history traits and habitat requirements, certainly affect connectivity in G. divaricata. However, this scenario of potential differentiation seems to be overridden by sporadic events of long-distance dispersal across barriers and large effective population sizes.

Along-shelf connectivity and circumpolar gene flow in Antarctic silverfish (Pleuragramma antarctica).

The Antarctic silverfish (Pleuragramma antarctica) is a critically important forage species with a circumpolar distribution and is unique among other notothenioid species for its wholly pelagic life cycle. Previous studies have provided mixed evidence of population structure over regional and circumpolar scales. The aim of the present study was to test the recent population hypothesis for Antarctic silverfish, which emphasizes the interplay between life history and hydrography in shaping connectivity. A total of 1067 individuals were collected over 25 years from different locations on a circumpolar scale. Samples were genotyped at fifteen microsatellites to assess population differentiation and genetic structuring using clustering methods, F-statistics, and hierarchical analysis of variance. A lack of differentiation was found between locations connected by the Antarctic Slope Front Current (ASF), indicative of high levels of gene flow. However, gene flow was significantly reduced at the South Orkney Islands and the western Antarctic Peninsula where the ASF is absent. This pattern of gene flow emphasized the relevance of large-scale circulation as a mechanism for circumpolar connectivity. Chaotic genetic patchiness characterized population structure over time, with varying patterns of differentiation observed between years, accompanied by heterogeneous standard length distributions. The present study supports a more nuanced version of the genetic panmixia hypothesis that reflects physical-biological interactions over the life history.

Mediterranean Bioconstructions Along the Italian Coast.

Marine bioconstructions are biodiversity-rich, three-dimensional biogenic structures, regulating key ecological functions of benthic ecosystems worldwide. Tropical coral reefs are outstanding for their beauty, diversity and complexity, but analogous types of bioconstructions are also present in temperate seas. The main bioconstructions in the Mediterranean Sea are represented by coralligenous formations, vermetid reefs, deep-sea cold-water corals, Lithophyllum byssoides trottoirs, coral banks formed by the shallow-water corals Cladocora caespitosa or Astroides calycularis, and sabellariid or serpulid worm reefs. Bioconstructions change the morphological and chemicophysical features of primary substrates and create new habitats for a large variety of organisms, playing pivotal roles in ecosystem functioning. In spite of their importance, Mediterranean bioconstructions have not received the same attention that tropical coral reefs have, and the knowledge of their biology, ecology and distribution is still fragmentary. All existing data about the spatial distribution of Italian bioconstructions have been collected, together with information about their growth patterns, dynamics and connectivity. The degradation of these habitats as a consequence of anthropogenic pressures (pollution, organic enrichment, fishery, coastal development, direct physical disturbance), climate change and the spread of invasive species was also investigated. The study of bioconstructions requires a holistic approach leading to a better understanding of their ecology and the application of more insightful management and conservation measures at basin scale, within ecologically coherent units based on connectivity: the cells of ecosystem functioning.

Population genomics meet Lagrangian simulations: Oceanographic patterns and long larval duration ensure connectivity among Paracentrotus lividus populations in the Adriatic and Ionian seas.

Connectivity between populations influences both their dynamics and the genetic structuring of species. In this study, we explored connectivity patterns of a marine species with long-distance dispersal, the edible common sea urchin Paracentrotus lividus, focusing mainly on the Adriatic-Ionian basins (Central Mediterranean). We applied a multidisciplinary approach integrating population genomics, based on 1,122 single nucleotide polymorphisms (SNPs) obtained from 2b-RAD in 275 samples, with Lagrangian simulations performed with a biophysical model of larval dispersal. We detected genetic homogeneity among eight population samples collected in the focal Adriatic-Ionian area, whereas weak but significant differentiation was found with respect to two samples from the Western Mediterranean (France and Tunisia). This result was not affected by the few putative outlier loci identified in our dataset. Lagrangian simulations found a significant potential for larval exchange among the eight Adriatic-Ionian locations, supporting the hypothesis of connectivity of P. lividus populations in this area. A peculiar pattern emerged from the comparison of our results with those obtained from published P. lividus cytochrome b (cytb) sequences, the latter revealing genetic differentiation in the same geographic area despite a smaller sample size and a lower power to detect differences. The comparison with studies conducted using nuclear markers on other species with similar pelagic larval durations in the same Adriatic-Ionian locations indicates species-specific differences in genetic connectivity patterns and warns against generalizing single-species results to the entire community of rocky shore habitats.

Population genomics of an endemic Mediterranean fish: differentiation by fine scale dispersal and adaptation.

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. For species with large effective population sizes a low number of markers may fail to identify population structure. A solution of this shortcoming can be high-throughput sequencing that allows genotyping thousands of markers on a genome-wide approach while facilitating the detection of genetic structuring shaped by selection. We used Genotyping-by-Sequencing (GBS) on 176 individuals of the endemic East Atlantic peacock wrasse (Symphodus tinca), from 6 locations in the Adriatic and Ionian seas. We obtained a total of 4,155 polymorphic SNPs and we observed two strong barriers to gene flow. The first one differentiated Tremiti Islands, in the northwest, from all the other locations while the second one separated east and south-west localities. Outlier SNPs potentially under positive selection and neutral SNPs both showed similar patterns of structuring, although finer scale differentiation was unveiled with outlier loci. Our results reflect the complexity of population genetic structure and demonstrate that both habitat fragmentation and positive selection are on play. This complexity should be considered in biodiversity assessments of different taxa, including non-model yet ecologically relevant organisms.

Unravelling population genetic structure with mitochondrial DNA in a notional panmictic coastal crab species: sample size makes the difference.

BACKGROUND: The extent of genetic structure of a species is determined by the amount of current gene flow and the impact of historical and demographic factors. Most marine invertebrates have planktonic larvae and consequently wide potential dispersal, so that genetic uniformity should be common. However, phylogeographic investigations reveal that panmixia is rare in the marine realm. Phylogeographic patterns commonly coincide with geographic transitions acting as barriers to gene flow. In the Mediterranean Sea and adjoining areas, the best known barriers are the Atlantic-Mediterranean transition, the Siculo-Tunisian Strait and the boundary between Aegean and Black seas. Here, we perform the so far broadest phylogeographic analysis of the crab Pachygrapsus marmoratus, common across the north-eastern Atlantic Ocean, Mediterranean and Black seas. Previous studies revealed no or weak genetic structuring at meso-geographic scale based on mtDNA, while genetic heterogeneity at local scale was recorded with microsatellites, even if without clear geographic patterns. Continuing the search for phylogeographic signal, we here enlarge the mtDNA dataset including 51 populations and covering most of the species' distribution range. RESULTS: This enlarged dataset provides new evidence of three genetically separable groups, corresponding to the Portuguese Atlantic Ocean, Mediterranean Sea plus Canary Islands, and Black Sea. Surprisingly, hierarchical AMOVA and Principal Coordinates Analysis agree that our Canary Islands population is closer to western Mediterranean populations than to mainland Portugal and Azores populations. Within the Mediterranean Sea, we record genetic homogeneity, suggesting that population connectivity is unaffected by the transition between the western and eastern Mediterranean. The Mediterranean metapopulation seems to have experienced a relatively recent expansion around 100,000 years ago. CONCLUSIONS: Our results suggest that the phylogeographic pattern of P. marmoratus is shaped by the geological history of Mediterranean and adjacent seas, restricted current gene flow among different marginal seas, and incomplete lineage sorting. However, they also caution from exclusively testing well-known biogeographic barriers, thereby neglecting other possible phylogeographic patterns. Mostly, this study provides evidence that a geographically exhaustive dataset is necessary to detect shallow phylogeographic structure within widespread marine species with larval dispersal, questioning all studies where species have been categorized as panmictic based on numerically and geographically limited datasets.

Transcriptomic profiling of male European eel (Anguilla anguilla) livers at sexual maturity.

The European eel Anguilla anguilla has a complex life cycle that includes freshwater, seawater and morphologically distinct stages as well as two extreme long distance migrations. Eels do not feed as they migrate across the Atlantic to the Sargasso Sea but nevertheless reach sexual maturity before spawning. It is not yet clear how existing energy stores are used to reach the appropriate developmental state for reproduction. Since the liver is involved in energy metabolism, protein biosynthesis and endocrine regulation it is expected to play a key role in the regulation of reproductive development. We therefore used microarrays to identify genes that may be involved in this process. Using this approach, we identified 231 genes that were expressed at higher and 111 genes that were expressed at lower levels in sexually mature compared with immature males. The up-regulated set includes genes involved in lipid metabolism, fatty acid synthesis and transport, mitochondrial function, steroid transport and bile acid metabolism. Several genes with putative enzyme functions were also expressed at higher levels at sexual maturity while genes involved in immune system processes and protein biosynthesis tended to be down-regulated at this stage. By using a high-throughput approach, we have identified a subset of genes that may be linked with the mobilization of energy stores for sexual maturation and migration. These results contribute to an improved understanding of eel reproductive biology and provide insight into the role of the liver in other teleosts with a long distance spawning migrations.

Multiple paternity and hybridization in two smooth-hound sharks.

Multiple paternity appears to be a common trait of elasmobranch mating systems, with its occurrence likely driven by convenience, due to females seeking to minimize the stress of male harassment. Here we use molecular markers to analyse the frequency of multiple paternity in two related viviparous sharks, Mustelus mustelus and Mustelus punctulatus. We first applied molecular methods to assign pregnant females, embryos and additional reference adults (N = 792) to one of the two species. Paternity analysis was performed using a total of 9 polymorphic microsatellites on 19 females and 204 embryos of M. mustelus, and on 13 females and 303 embryos of M. punctulatus. Multiple paternity occurs in both species, with 47% of M. mustelus and 54% of M. punctulatus litters sired by at least two fathers. Female fecundity is not influenced by multiple mating and in 56% of polyandrous litters paternity is skewed, with one male siring most of the pups. Genetic analyses also revealed hybridization between the two species, with a M. punctulatus female bearing pups sired by a M. mustelus male. The frequency of polyandrous litters in these species is consistent with aspects of their reproductive biology, such as synchronous ovulation and possible occurrence of breeding aggregations.

New molecular tools for the identification of 2 endangered smooth-hound sharks, Mustelus mustelus and Mustelus punctulatus.

The smooth-hounds represent a significant proportion of the elasmobranch catch in the Adriatic basin of the Mediterranean Sea, where the common (Mustelus mustelus) and blackspotted (Mustelus punctulatus) smooth-hounds co-occur. The 2 species share several morphological and morphometric characters that lead to frequent misidentification. In order to provide information useful for their species identification, we performed a morphological identification of several Mustelus specimens to select individuals unambiguously attributed to 1 of the 2 species, and assayed these with 3 new molecular tests. First, we developed and validated a mitochondrial DNA assay based on species-specific amplification of the cytochrome c oxidase subunit 1 (COI). Second, a fragment analysis of 15 microsatellites cross-amplified from several triakid species was performed to identify diagnostic loci. Finally, a length difference was identified in the internal transcribed spacer 2 (ITS2) region and a diagnostic test based on its amplification was established. All the samples classified morphologically as M. mustelus and M. punctulatus showed a species-specific profile using all the 3 molecular tests. In addition, cross-amplification of microsatellites allowed identification of 9 highly polymorphic loci that will be useful for the study of the mating system and population differentiation of the 2 species.

Changes in the gene expression profiles of the brains of male European eels (Anguilla anguilla) during sexual maturation.

BACKGROUND: The vertebrate brain plays a critical role in the regulation of sexual maturation and reproduction by integrating environmental information with developmental and endocrine status. The European eel Anguilla anguilla is an important species in which to better understand the neuroendocrine factors that control reproduction because it is an endangered species, has a complex life cycle that includes two extreme long distance migrations with both freshwater and seawater stages and because it occupies a key position within the teleost phylogeny. At present, mature eels have never been caught in the wild and little is known about most aspects of reproduction in A. anguilla. The goal of this study was to identify genes that may be involved in sexual maturation in experimentally matured eels. For this, we used microarrays to compare the gene expression profiles of sexually mature to immature males. RESULTS: Using a false discovery rate of 0.05, a total of 1,497 differentially expressed genes were identified. Of this set, 991 were expressed at higher levels in brains (forebrain and midbrain) of mature males while 506 were expressed at lower levels relative to brains of immature males. The set of up-regulated genes includes genes involved in neuroendocrine processes, cell-cell signaling, neurogenesis and development. Interestingly, while genes involved in immune system function were down-regulated in the brains of mature males, changes in the expression levels of several receptors and channels were observed suggesting that some rewiring is occurring in the brain at sexual maturity. CONCLUSIONS: This study shows that the brains of eels undergo major changes at the molecular level at sexual maturity that may include re-organization at the cellular level. Here, we have defined a set of genes that help to understand the molecular mechanisms controlling reproduction in eels. Some of these genes have previously described functions while many others have roles that have yet to be characterized in a reproductive context. Since most of the genes examined here have orthologs in other vertebrates, the results of this study will contribute to the body of knowledge concerning reproduction in vertebrates as well as to an improved understanding of eel biology.

First evidence of inbreeding, relatedness and chaotic genetic patchiness in the holoplanktonic jellyfish Pelagia noctiluca (Scyphozoa, Cnidaria).

Genetic drift and non-random mating seldom influence species with large breeding populations and high dispersal potential, characterized by unstructured gene pool and panmixia at a scale lower than the minimum dispersal range of individuals. In the present study, a set of nine microsatellite markers was developed and used to investigate the spatio-temporal genetic patterns of the holoplanktonic jellyfish Pelagia noctiluca (Scyphozoa) in the Southern Tyrrhenian Sea. Homozygote excess was detected at eight loci, and individuals exhibited intra-population relatedness higher than expected by chance in at least three samples. This result was supported by the presence of siblings in at least 5 out 8 samples, 4 of which contained full-sib in addition to half-sib dyads. Having tested and ruled out alternative explanations as null alleles, our results suggest the influence of reproductive and behavioural features in shaping the genetic structure of P. noctiluca, as outcomes of population genetics analyses pointed out. Indeed, the genetic differentiation among populations was globally small but highlighted: a) a spatial genetic patchiness uncorrelated with distance between sampling locations, and b) a significant genetic heterogeneity between samples collected in the same locations in different years. Therefore, despite its extreme dispersal potential, P. noctiluca does not maintain a single homogenous population, but rather these jellyfish appear to have intra-bloom localized recruitment and/or individual cohesiveness, whereby siblings more likely swarm together as a single group and remain close after spawning events. These findings provide the first evidence of family structures and consequent genetic patchiness in a species with highly dispersive potential throughout its whole life cycle, contributing to understanding the patterns of dispersal and connectivity in marine environments.

Genome evolution in the cold: Antarctic icefish muscle transcriptome reveals selective duplications increasing mitochondrial function.

Antarctic notothenioids radiated over millions of years in subzero waters, evolving peculiar features, such as antifreeze glycoproteins and absence of heat shock response. Icefish, family Channichthyidae, also lack oxygen-binding proteins and display extreme modifications, including high mitochondrial densities in aerobic tissues. A genomic expansion accompanying the evolution of these fish was reported, but paucity of genomic information limits the understanding of notothenioid cold adaptation. We reconstructed and annotated the first skeletal muscle transcriptome of the icefish Chionodraco hamatus providing a new resource for icefish genomics (http://compgen.bio.unipd.it/chamatusbase/, last accessed December 12, 2012). We exploited deep sequencing of this energy-dependent tissue to test the hypothesis of selective duplication of genes involved in mitochondrial function. We developed a bioinformatic approach to univocally assign C. hamatus transcripts to orthology groups extracted from phylogenetic trees of five model species. Chionodraco hamatus duplicates were recorded for each orthology group allowing the identification of duplicated genes specific to the icefish lineage. Significantly more duplicates were found in the icefish when transcriptome data were compared with whole-genome data of model species. Indeed, duplicated genes were significantly enriched in proteins with mitochondrial localization, involved in mitochondrial function and biogenesis. In cold conditions and without oxygen-carrying proteins, energy production is challenging. The combination of high mitochondrial densities and the maintenance of duplicated genes involved in mitochondrial biogenesis and aerobic respiration might confer a selective advantage by improving oxygen diffusion and energy supply to aerobic tissues. Our results provide new insights into the genomic basis of icefish cold adaptation.

Body size correlates with fertilization success but not gonad size in grass goby territorial males.

In fish species with alternative male mating tactics, sperm competition typically occurs when small males that are unsuccessful in direct contests steal fertilization opportunities from large dominant males. In the grass goby Zosterisessor ophiocephalus, large territorial males defend and court females from nest sites, while small sneaker males obtain matings by sneaking into nests. Parentage assignment of 688 eggs from 8 different nests sampled in the 2003-2004 breeding season revealed a high level of sperm competition. Fertilization success of territorial males was very high but in all nests sneakers also contributed to the progeny. In territorial males, fertilization success correlated positively with male body size. Gonadal investment was explored in a sample of 126 grass gobies collected during the period 1995-1996 in the same area (61 territorial males and 65 sneakers). Correlation between body weight and testis weight was positive and significant for sneaker males, while correlation was virtually equal to zero in territorial males. That body size in territorial males is correlated with fertilization success but not gonad size suggests that males allocate much more energy into growth and relatively little into sperm production once the needed size to become territorial is attained. The increased paternity of larger territorial males might be due to a more effective defense of the nest in comparison with smaller territorial males.

Sequencing and characterization of striped venus transcriptome expand resources for clam fishery genetics.

BACKGROUND: The striped venus Chamelea gallina clam fishery is among the oldest and the largest in the Mediterranean Sea, particularly in the inshore waters of northern Adriatic Sea. The high fishing pressure has lead to a strong stock abundance decline, enhanced by several irregular mortality events. The nearly complete lack of molecular characterization limits the available genetic resources for C. gallina. We achieved the first transcriptome of this species with the aim of identifying an informative set of expressed genes, potential markers to assess genetic structure of natural populations and molecular resources for pathogenic contamination detection. METHODOLOGY/PRINCIPAL FINDINGS: The 454-pyrosequencing of a normalized cDNA library of a pool C. gallina adult individuals yielded 298,494 raw reads. Different steps of reads assembly and filtering produced 36,422 contigs of high quality, one half of which (18,196) were annotated by similarity. A total of 111 microsatellites and 20,377 putative SNPs were identified. A panel of 13 polymorphic transcript-linked microsatellites was developed and their variability assessed in 12 individuals. Remarkably, a scan to search for contamination sequences of infectious origin indicated the presence of several Vibrionales species reported to be among the most frequent clam pathogen's species. Results reported in this study were included in a dedicated database available at http://compgen.bio.unipd.it/chameleabase. CONCLUSIONS/SIGNIFICANCE: This study represents the first attempt to sequence and de novo annotate the transcriptome of the clam C. gallina. The availability of this transcriptome opens new perspectives in the study of biochemical and physiological role of gene products and their responses to large and small-scale environmental stress in C. gallina, with high throughput experiments such as custom microarray or targeted re-sequencing. Molecular markers, such as the already optimized EST-linked microsatellites and the discovered SNPs will be useful to estimate effects of demographic processes and to detect minute levels of population structuring.

Surviving in a toxic world: transcriptomics and gene expression profiling in response to environmental pollution in the critically endangered European eel.

BACKGROUND: Genomic and transcriptomic approaches have the potential for unveiling the genome-wide response to environmental perturbations. The abundance of the catadromous European eel (Anguilla anguilla) stock has been declining since the 1980s probably due to a combination of anthropogenic and climatic factors. In this paper, we explore the transcriptomic dynamics between individuals from high (river Tiber, Italy) and low pollution (lake Bolsena, Italy) environments, which were measured for 36 PCBs, several organochlorine pesticides and brominated flame retardants and nine metals. RESULTS: To this end, we first (i) updated the European eel transcriptome using deep sequencing data with a total of 640,040 reads assembled into 44,896 contigs (Eeelbase release 2.0), and (ii) developed a transcriptomic platform for global gene expression profiling in the critically endangered European eel of about 15,000 annotated contigs, which was applied to detect differentially expressed genes between polluted sites. Several detoxification genes related to metabolism of pollutants were upregulated in the highly polluted site, including genes that take part in phase I of the xenobiotic metabolism (CYP3A), phase II (glutathione-S-transferase) and oxidative stress (glutathione peroxidase). In addition, key genes in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated at the Tiber site relative to the Bolsena site. CONCLUSIONS: Together with the induced high expression of detoxification genes, the suggested lowered expression of genes supposedly involved in metabolism suggests that pollution may also be associated with decreased respiratory and energy production.