European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation

From population connectivity to the art of striping Russian dolls: the lessons from Pocillopora corals

Here, we examined the genetic variability in the coral genus Pocillopora, in particular within the Primary Species Hypothesis PSH09, identified by Gélin, Postaire, Fauvelot and Magalon (2017) using species delimitation methods [also named Pocillopora eydouxi/meandrina complex sensu, Schmidt-Roach, Miller, Lundgren, & Andreakis (2014)] and which was found to split into three secondary species hypotheses (SSH09a, SSH09b, and SSH09c) according to assignment tests using multi-locus genotypes (13 microsatellites). From a large sampling (2,507 colonies) achieved in three marine provinces [Western Indian Ocean (WIO), Tropical Southwestern Pacific (TSP), and Southeast Polynesia (SEP)], genetic structuring analysis conducted with two clustering analyses (structure and DAPC) using 13 microsatellites revealed that SSH09a was restricted to the WIO while SSH09b and SSH09c were almost exclusively in the TSP and SEP. More surprisingly, each SSH split into two to three genetically differentiated clusters, found in sympatry at the reef scale, leading to a pattern of nested hierarchical levels (PSH > SSH > cluster), each level hiding highly differentiated genetic groups. Thus, rather than structured populations within a single species, these three SSHs, and even the eight clusters, likely represent distinct genetic lineages engaged in a speciation process or real species. The issue is now to understand which hierarchical level (SSH, cluster, or even below) corresponds to the species one. Several hypotheses are discussed on the processes leading to this pattern of mixed clusters in sympatry, evoking formation of reproductive barriers, either by allopatric speciation or habitat selection.

Recombination in the eggs and sperm in a simultaneously hermaphroditic vertebrate

When there is no recombination (achiasmy) in one sex, it is in the heterogametic one. This observation is so consistent that it constitutes one of the few patterns in biology that may be regarded as a 'rule' and Haldane (Haldane 1922 J. Genet. 12, 101-109. (doi:10.1007/BF02983075)) proposed that it might be driven by selection against recombination in the sex chromosomes. Yet differences in recombination rates between the sexes (heterochiasmy) have also been reported in hermaphroditic species that lack sex chromosomes. In plants-the vast majority of which are hermaphroditic-selection at the haploid stage has been proposed to drive heterochiasmy. Yet few data are available for hermaphroditic animals, and barely any for hermaphroditic vertebrates. Here, we leverage reciprocal crosses between two black hamlets (Hypoplectrus nigricans, Serranidae), simultaneously hermaphroditic reef fishes from the wider Caribbean, to generate high-density egg- and sperm-specific linkage maps for each parent. We find globally higher recombination rates in the eggs, with dramatically pronounced heterochiasmy at the chromosome peripheries. We suggest that this pattern may be due to female meiotic drive, and that this process may be an important source of heterochiasmy in animals. We also identify a large non-recombining region that may play a role in speciation and local adaptation in Hypoplectrus.

Fantastic Beasts and How To Sequence Them: Ecological Genomics for Obscure Model Organisms

The application of genomic approaches to 'obscure model organisms' (OMOs), meaning species with no prior genomic resources, enables increasingly sophisticated studies of the genomic basis of evolution, acclimatization, and adaptation in real ecological contexts. I consider here ecological questions that can be addressed using OMOs, and indicate optimal sequencing and data-handling solutions for each case. With this I hope to promote the diversity of OMO-based projects that would capitalize on the peculiarities of the natural history of OMOs and could feasibly be completed within the scope of a single PhD thesis.

Genomics of adaptive divergence with chromosome-scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome-scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low-crossover centres relative to the high-crossover peripheries of chromosomes ("Chromosome Centre-Biased Differentiation", CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome-scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome-scale heterogeneity in crossover rate in evolutionary genomics.

The rises and falls of opsin genes in 59 ray-finned fish genomes and their implications for environmental adaptation

We studied the evolution of opsin genes in 59 ray-finned fish genomes. We identified the opsin genes and adjacent genes (syntenies) in each genome. Then we inferred the changes in gene copy number (N), syntenies, and tuning sites along each phylogenetic branch during evolution. The Exorh (rod opsin) gene has been retained in 56 genomes. Rh1, the intronless rod opsin gene, first emerged in ancestral Actinopterygii, and N increased to 2 by the teleost-specific whole genome duplication, but then decreased to 1 in the ancestor of Neoteleostei fishes. For cone opsin genes, the rhodopsin-like (Rh2) and long-wave-sensitive (LWS) genes showed great variation in N among species, ranging from 0 to 5 and from 0 to 4, respectively. The two short-wave-sensitive genes, SWS1 and SWS2, were lost in 23 and 6 species, respectively. The syntenies involving LWS, SWS2 and Rh2 underwent complex changes, while the evolution of the other opsin gene syntenies was much simpler. Evolutionary adaptation in tuning sites under different living environments was discussed. Our study provides a detailed view of opsin gene gains and losses, synteny changes and tuning site changes during ray-finned fish evolution.

Innate Cell-Mediated Cytotoxic Activity of European Sea Bass Leucocytes Against Nodavirus-Infected Cells: A Functional and RNA-seq Study

Nervous necrosis virus (NNV) causes high mortalities in several marine species. We aimed to evaluate the innate cell-mediated cytotoxic (CMC) activity of head-kidney leucocytes (HKLs) isolated from naïve European sea bass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata), a very susceptible and resistant fish species to NNV, respectively, against fish cell lines infected with NNV. Seabream HKLs showed significantly increased innate CMC activity against NNV-infected cells, compared to those uninfected, while sea bass HKLs failed to do so. Thus, we performed a RNA-seq study to identify genes related to the CMC activity of sea bass leucocytes. Thus, we found that sea bass HKLs incubated with DLB-1 cells alone (CMC_DLB1) or with NNV-infected DLB-1 cells (CMC_DLB1-NNV) showed very similar transcriptomic profiles and the GO analysis revealed that most of the up-regulated genes were related to immunity. Strikingly, when the CMC samples with and without NNV were compared, GO analysis revealed that most of the up-regulated genes in CMC_DLB1-NNV samples were related to metabolism and very few to immunity. This is also in agreement with the functional data. These data point to the escape of CMC activity by NNV infection as an important factor involved in the high susceptibility to nodavirus infections of European sea bass.

Genomic differentiation and patterns of gene flow between two long-tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions ("genomic differentiation islands") spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene-flow patterns for autosomes using RAD-seq data between two closely related species of long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.

Nutrigenomics and immune function in fish: new insights from omics technologies

The interplay between nutrition and immune system is well recognised, however the true integration of research between nutrition, animal energy status and immune function is still far from clear. In fish nutrition, especially for species maintained in aquaculture, formulated feeds are significantly different from the natural diet with recent changes in nutrient sources, especially with protein and oil sources now being predominated by terrestrial derived ingredients. Additionally, many feeds are now incorporated to health management and termed functional feeds, which are believed to improve fish health, reduce disease outbreaks and/or improve post-infection recovery. Using new omics technologies, including transcriptomics (microarray and RNA-seq) and proteomics, the impacts of nutrition on the immune system is becoming clearer. By using molecular pathway enrichment analysis, modules of genes can indicate how both local (intestinal) and systemic immune function are being altered. Although great progress has been made to define the changes in host immune function, understanding the interplay between fish nutrition, intestinal microbiome and immune system is only just beginning to emerge.

Small ocean temperature increases elicit stage-dependent changes in DNA methylation and gene expression in a fish, the European sea bass

In natural fish populations, temperature increases can result in shifts in important phenotypic traits. DNA methylation is an epigenetic mechanism mediating phenotypic changes. However, whether temperature increases of the magnitude predicted by the latest global warming models can affect DNA methylation is unknown. Here, we exposed European sea bass to moderate temperature increases in different periods within the first two months of age. We show that increases of even 2 °C in larvae significantly changed global DNA methylation and the expression of ecologically-relevant genes related to DNA methylation, stress response, muscle and organ formation, while 4 °C had no effect on juveniles. Furthermore, DNA methylation changes were more marked in larvae previously acclimated to a different temperature. The expression of most genes was also affected by temperature in the larvae but not in juveniles. In conclusion, this work constitutes the first study of DNA methylation in fish showing that temperature increases of the magnitude predicted by the latest global warming models result in stage-dependent alterations in global DNA methylation and gene expression levels. This study, therefore, provides insights on the possible consequences of climate change in fish mediated by genome-wide epigenetic modifications.

Unexpected mosaic distribution of two hybridizing sibling lineages in the teleplanically dispersing snail Stramonita haemastoma suggests unusual postglacial redistribution or cryptic invasion

Molecular approaches have proven efficient to identify cryptic lineages within single taxonomic entities. Sometimes these cryptic lineages maybe previously unreported or unknown invasive taxa. The genetic structure of the marine gastropod Stramonita haemastoma has been examined in the Western Mediterranean and North‐Eastern Atlantic populations with mtDNA COI sequences and three newly developed microsatellite markers. We identified two cryptic lineages, differentially fixed for alternative mtDNA COI haplogroups and significantly differentiated at microsatellite loci. The mosaic distribution of the two lineages is unusual for a warm‐temperate marine invertebrate with a teleplanic larval stage. The Atlantic lineage was unexpectedly observed as a patch enclosed in the north of the Western Mediterranean Sea between eastern Spain and the French Riviera, and the Mediterranean lineage was found in Macronesian Islands. Although cyto‐nuclear disequilibrium is globally maintained, asymmetric introgression occurs in the Spanish region where the two lineages co‐occur in a hybrid zone. A first interpretation of our results is mito‐nuclear discordance in a stable postglacial hybrid zone. Under this hypothesis, though, the location of genetic discontinuities would be unusual among planktonic dispersers. An alternative interpretation is that the Atlantic lineage, also found in Senegal and Venezuela, has been introduced by human activities in the Mediterranean area and is introgressing Mediterranean genes during its propagation, as theoretically expected. This second hypothesis would add an additional example to the growing list of cryptic marine invasions revealed by molecular studies.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

Genomewide patterns of variation in genetic diversity are shared among populations, species and higher-order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e ) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage-specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e . Utilizing population-level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population-scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST , PBS, d_xy ) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

Gene-centromere mapping in meiotic gynogenetic European seabass

Background

Fully isogenic lines in fish can be developed using “mitotic” gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families.

Results

From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly.

Conclusion

Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3826-z) contains supplementary material, which is available to authorized users.

Transcriptomic differentiation underlying marine-to-freshwater transitions in the South American silversides Odontesthes argentinensis and O. bonariensis (Atheriniformes)

Salinity gradients are critical habitat determinants for freshwater organisms. Silverside fishes in the genus Odontesthes have recently and repeatedly transitioned from marine to freshwater habitats, overcoming a strong ecological barrier. Genomic and transcriptomic changes involved in this kind of transition are only known for a few model species. We present new data and analyses of gene expression and microbiome composition in the gills of two closely related silverside species, marine O. argentinensis and freshwater O. bonariensis and find more than three thousand transcripts differentially expressed, with osmoregulatory/ion transport genes and immune genes showing very different expression patterns across species. Interspecific differences also involve more than one thousand transcripts with nonsynonymous SNPs in the coding sequences, most of which were not differentially expressed. In addition to characterizing gill transcriptomes from wild‐caught marine and freshwater fishes, we test experimentally the response to salinity increases by O. bonariensis collected from freshwater habitats. Patterns of expression in gill transcriptomes of O. bonariensis exposed to high salinity do not resemble O. argentinensis mRNA expression, suggesting lack of plasticity for adaptation to marine conditions in this species. The diversity of functions associated with both the differentially expressed set of transcripts and those with sequence divergence plus marked microbiome differences suggest that multiple abiotic and biotic factors in marine and freshwater habitats are driving transcriptomic differences between these species.

Identification of conserved genes triggering puberty in European sea bass males (Dicentrarchus labrax) by microarray expression profiling

BACKGROUND

Spermatogenesis is a complex process characterized by the activation and/or repression of a number of genes in a spatio-temporal manner. Pubertal development in males starts with the onset of the first spermatogenesis and implies the division of primary spermatogonia and their subsequent entry into meiosis. This study is aimed at the characterization of genes involved in the onset of puberty in European sea bass, and constitutes the first transcriptomic approach focused on meiosis in this species.

RESULTS

European sea bass testes collected at the onset of puberty (first successful reproduction) were grouped in stage I (resting stage), and stage II (proliferative stage). Transition from stage I to stage II was marked by an increase of 11ketotestosterone (11KT), the main fish androgen, whereas the transcriptomic study resulted in 315 genes differentially expressed between the two stages. The onset of puberty induced 1) an up-regulation of genes involved in cell proliferation, cell cycle and meiosis progression, 2) changes in genes related with reproduction and growth, and 3) a down-regulation of genes included in the retinoic acid (RA) signalling pathway. The analysis of GO-terms and biological pathways showed that cell cycle, cell division, cellular metabolic processes, and reproduction were affected, consistent with the early events that occur during the onset of puberty. Furthermore, changes in the expression of three RA nuclear receptors point at the importance of the RA-signalling pathway during this period, in agreement with its role in meiosis.

CONCLUSION

The results contribute to boost our knowledge of the early molecular and endocrine events that trigger pubertal development and the onset of spermatogenesis in fish. These include an increase in 11KT plasma levels and changes in the expression of several genes involved in cell proliferation, cell cycle progression, meiosis or RA-signalling pathway. Moreover, the results can be applied to study meiosis in this economically important fish species for Mediterranean countries, and may help to develop tools for its sustainable aquaculture.

Evolution of Th2 responses: characterization of IL-4/13 in sea bass (Dicentrarchus labrax L.) and studies of expression and biological activity

Th2 immunity is a primary host defence against metazoan pathogens and two of the important cytokines involved in this immune response in mammals are IL-4 and IL-13. Recently the origin and evolution of Th2 immune responses have been investigated in fish where a molecule with relatedness to both IL-4 and IL-13 is present, termed IL-4/13. Different IL-4/13 paralogues (IL-4/13 A and IL-4/13B) exist in teleost fish. In this paper, we have focused on the IL-4/13 isoforms found in the European sea bass (Dicentrarchus labrax L.). Two tandem duplicated but divergent IL-4/13 A isoforms and one IL-4/13B are present, a unique situation compared to other teleosts. These genes were studied in terms of their in vitro and in vivo transcript levels after different treatments and their biological activities after production of the recombinant isoforms. The results show that the presence of these three paralogues is associated with different activities, both in terms of their expression profiles and the ability of the proteins to modulate the expression of immune genes in head kidney leukocytes. It is clear that the initiation and control of type-2 responses in seabass is complex due to the presence of multiple IL-4/13 isoforms with overlapping but distinct activities.

Viral encephalopathy and retinopathy in aquaculture: a review

Viral encephalopathy and retinopathy (VER), otherwise known as viral nervous necrosis (VNN), is a major devastating threat for aquatic animals. Betanodaviruses have been isolated in at least 70 aquatic animal species in marine and in freshwater environments throughout the world, with the notable exception of South America. In this review, the main features of betanodavirus, including its diversity, its distribution and its transmission modes in fish, are firstly presented. Then, the existing diagnosis and detection methods, as well as the different control procedures of this disease, are reviewed. Finally, the potential of selective breeding, including both conventional and genomic selection, as an opportunity to obtain resistant commercial populations, is examined.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.

Extensive recent secondary contacts between four European white oak species

Historical trajectories of tree species during the late Quaternary have been well reconstructed through genetic and palaeobotanical studies. However, many congeneric tree species are interfertile, and the timing and contribution of introgression to species divergence during their evolutionary history remains largely unknown. We quantified past and current gene flow events between four morphologically divergent oak species (Quercus petraea, Q. robur, Q. pyrenaica, Q. pubescens), by two independent inference methods: diffusion approximation to the joint frequency spectrum (∂a∂i) and approximate Bayesian computation (ABC). For each pair of species, alternative scenarios of speciation allowing gene flow over different timescales were evaluated. Analyses of 3524 single nucleotide polymorphisms (SNPs) randomly distributed in the genome, showed that these species evolved in complete isolation for most of their history, but recently came into secondary contact, probably facilitated by the most recent period of postglacial warming. We demonstrated that: there was sufficient genetic differentiation before secondary contact for the accumulation of barriers to gene flow; and current European white oak genomes are a mosaic of genes that have crossed species boundaries and genes impermeable to gene flow.

Metapopulation patterns of additive and nonadditive genetic variance in the sea bass (Dicentrarchus labrax)

Describing and explaining the geographic within‐species variation in phenotypes (“phenogeography”) in the sea over a species distribution range is central to our understanding of a variety of eco‐evolutionary topics. However, phenogeographic studies that have a large potential to investigate adaptive variation are overcome by phylogeographic studies, still mainly focusing on neutral markers. How genotypic and phenotypic data could covary over large geographic scales remains poorly understood in marine species. We crossed 75 noninbred sires (five origins) and 26 dams (two origins; each side of a hybrid zone) in a factorial diallel cross in order to investigate geographic variation for early survival and sex ratio in the metapopulation of the European sea bass (Dicentrarchus labrax), a highly prized marine fish species. Full‐sib families (N = 1,950) were produced and reared in a common environment. Parentage assignment of 7,200 individuals was performed with seven microsatellite markers. Generalized linear models showed significant additive effects for both traits and pleiotropy between traits. A significant nonadditive genetic effect was detected. Different expression of traits and distinct relative performances were found for reciprocal crosses involving populations located on each side of the main hybrid zone located at the Almeria‐Oran front, illustrating asymmetric reproductive isolation. The poor fitness performance observed for the Western Mediterranean population of sea bass is discussed as it represents the main source of seed hatchery production, but also because it potentially illustrates nonadaptive introgression and maladaptation.

Immunoglobulin T from sea bass (Dicentrarchus labrax L.): molecular characterization, tissue localization and expression after nodavirus infection

Background

Immunoglobulins (Igs) are fundamental components of the adaptive immune system of vertebrates, with the IgT/IgZ isotype specific of Teleosts. In this paper we describe the identification of an IgT heavy chain from the European sea bass (Dicentrarchus labrax L.), its molecular characterization and tissue mRNA localization by in situ hybridization.

Results

Sea bass IgT consists of 552 aa (Accession Number KM410929) and it contains a putative 19 amino acids long signal peptide and one potential N-glycosylation site. The C-region consists of four C_H domains; each contains the cysteine and tryptophan residues required for their correct folding. Based on the recent sequencing of sea bass genome, we have identified five different genomic contigs bearing exons unequivocally pertaining to IgT (C_H2, C_H3 and C_H4), but none corresponded to a complete IgH locus as IgT sequences were found in the highly fragmented assembled genomic regions which could not be assigned to any major scaffold. The 3D structure of sea bass IgT has been modelled using the crystal structure of a mouse Ig gamma as a template, thus showing that the amino acid sequence is suitable for the expected topology referred to an immunoglobulin-like architecture. The basal expression of sea bass IgT and IgM in different organs has been analysed: gut and gills, important mucosal organs, showed high IgT transcripts levels and this was the first indication of the possible involvement of sea bass IgT in mucosal immune responses. Moreover, sea bass IgT expression increased in gills and spleen after infection with nodavirus, highlighting the importance of IgT in sea bass immune responses. In situ hybridization confirmed the presence of IgT transcripts in the gut and it revealed a differential expression along the intestinal tract, with a major expression in the posterior intestine, suggesting the hindgut as a site for the recruitment of IgT⁺ cells in this species. IgT transcripts were also found in gill filaments and parallel lamellae and, for the first time, we identified scattered IgT positive cells in the liver, with a strong signal in the hepatic parenchyma.

Conclusions

In conclusion, we performed a full molecular characterization of IgT in sea bass that points out its possible involvement in mucosal immune responses of this species.

Linkage Mapping and Comparative Genomics of Red Drum (Sciaenops ocellatus) Using Next-Generation Sequencing

Developments in next-generation sequencing allow genotyping of thousands of genetic markers across hundreds of individuals in a cost-effective manner. Because of this, it is now possible to rapidly produce dense genetic linkage maps for nonmodel species. Here, we report a dense genetic linkage map for red drum, a marine fish species of considerable economic importance in the southeastern United States and elsewhere. We used a prior microsatellite-based linkage map as a framework and incorporated 1794 haplotyped contigs derived from high-throughput, reduced representation DNA sequencing to produce a linkage map containing 1794 haplotyped restriction-site associated DNA (RAD) contigs, 437 anonymous microsatellites, and 44 expressed sequence-tag-linked microsatellites (EST-SSRs). A total of 274 candidate genes, identified from transcripts from a preliminary hydrocarbon exposure study, were localized to specific chromosomes, using a shared synteny approach. The linkage map will be a useful resource for red drum commercial and restoration aquaculture, and for better understanding and managing populations of red drum in the wild.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Construction of high-resolution recombination maps in Asian seabass

Background

A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.

Results

A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.

Conclusions

These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3462-z) contains supplementary material, which is available to authorized users.

Shedding Light on the Grey Zone of Speciation along a Continuum of Genomic Divergence

Speciation results from the progressive accumulation of mutations that decrease the probability of mating between parental populations or reduce the fitness of hybrids—the so-called species barriers. The speciation genomic literature, however, is mainly a collection of case studies, each with its own approach and specificities, such that a global view of the gradual process of evolution from one to two species is currently lacking. Of primary importance is the prevalence of gene flow between diverging entities, which is central in most species concepts and has been widely discussed in recent years. Here, we explore the continuum of speciation thanks to a comparative analysis of genomic data from 61 pairs of populations/species of animals with variable levels of divergence. Gene flow between diverging gene pools is assessed under an approximate Bayesian computation (ABC) framework. We show that the intermediate "grey zone" of speciation, in which taxonomy is often controversial, spans from 0.5% to 2% of net synonymous divergence, irrespective of species life history traits or ecology. Thanks to appropriate modeling of among-locus variation in genetic drift and introgression rate, we clarify the status of the majority of ambiguous cases and uncover a number of cryptic species. Our analysis also reveals the high incidence in animals of semi-isolated species (when some but not all loci are affected by barriers to gene flow) and highlights the intrinsic difficulty, both statistical and conceptual, of delineating species in the grey zone of speciation.

Isolated populations accumulate genetic differences across their genomes as they diverge, whereas gene flow between populations counteracts divergence and tends to restore genetic homogeneity. Speciation proceeds by the accumulation at specific loci of mutations that reduce the fitness of hybrids, therefore preventing gene flow—the so-called species barriers. Importantly, species barriers are expected to act locally within the genome, leading to the prediction of a mosaic pattern of genetic differentiation between populations at intermediate levels of divergence—the genic view of speciation. At the same time, linked selection also contributes to speed up differentiation in low-recombining and gene-dense regions. We used a modelling approach that accounts for both sources of genomic heterogeneity and explored a wide continuum of genomic divergence made by 61 pairs of species/populations in animals. Our analysis provides a unifying picture of the relationship between molecular divergence and ability to exchange genes. We show that the "grey zone" of speciation—the intermediate state in which species definition is controversial—spans from 0.5% to 2% of molecular divergence, with these thresholds being independent of species life history traits and ecology. Semi-isolated species, between which alleles can be exchanged at some but not all loci, are numerous, with the earliest species barriers being detected at divergences as low as 0.075%. These results have important implications regarding taxonomy, conservation biology, and the management of biodiversity.

An integrated genomic approach for the study of mandibular prognathism in the European seabass (Dicentrarchus labrax)

Skeletal anomalies in farmed fish are a relevant issue affecting animal welfare and health and causing significant economic losses. Here, a high-density genetic map of European seabass for QTL mapping of jaw deformity was constructed and a genome-wide association study (GWAS) was carried out on a total of 298 juveniles, 148 of which belonged to four full-sib families. Out of 298 fish, 107 were affected by mandibular prognathism (MP). Three significant QTLs and two candidate SNPs associated with MP were identified. The two GWAS candidate markers were located on ChrX and Chr17, both in close proximity with the peaks of the two most significant QTLs. Notably, the SNP marker on Chr17 was positioned within the Sobp gene coding region, which plays a pivotal role in craniofacial development. The analysis of differentially expressed genes in jaw-deformed animals highlighted the “nervous system development” as a crucial pathway in MP. In particular, Zic2, a key gene for craniofacial morphogenesis in model species, was significantly down-regulated in MP-affected animals. Gene expression data revealed also a significant down-regulation of Sobp in deformed larvae. Our analyses, integrating transcriptomic and GWA methods, provide evidence for putative mechanisms underlying seabass jaw deformity.

Molecular characterization and expression of Na+/K+-ATPase α1 isoforms in the European sea bass Dicentrarchus labrax osmoregulatory tissues following salinity transfer

The Na⁺/K⁺-ATPase (NKA) is considered as the main pump involved in active ion transport. In the European sea bass, Dicentrarchus labrax, we found two genes encoding for the alpha 1 subunit isoforms (NKA α1a and NKA α1b). NKA α1a and NKA α1b isoform amino acid (aa) sequences were compared through phylogeny and regarding key functional motifs between salmonids and other acanthomorph species. Analysis of aa sequences of both isoforms revealed a high degree of conservation across teleosts. The expression pattern of both nka α1a and nka α1b was measured in the gill, kidney and posterior intestine of fish in seawater (SW) and transferred to fresh water (FW) at different exposure times. Nka α1a was more expressed than nka α1b whatever the condition and the tissue analyzed. After long-term salinity acclimation (2.5 years) either in FW or SW, transcript levels of nka α1a were higher in the kidney followed by the posterior intestine and the gill. Compared to SW conditions, expression of nka α1a in FW was significantly increased or decreased, respectively, in gill and posterior intestine. In contrast, branchial nka α1b was significantly decreased in FW-acclimated fish. Short-term FW acclimation seems to rapidly increase nka α1a transcript levels in the kidney unlike in gill tissues where different gene expression levels are detected only after long-term acclimation.

Genome-wide nucleotide diversity of hatchery-reared Atlantic and Mediterranean strains of brown trout Salmo trutta compared to wild Mediterranean populations

A genome-wide assessment of diversity is provided for wild Mediterranean brown trout Salmo trutta populations from headwater tributaries of the Orb River and from Atlantic and Mediterranean hatchery-reared strains that have been used for stocking. Double-digest restriction-site-associated DNA sequencing (dd-RADseq) was performed and the efficiency of de novo and reference-mapping approaches to obtain individual genotypes was compared. Large numbers of single nucleotide polymorphism (SNP) markers with similar genome-wide distributions were discovered using both approaches (196 639 v. 121 016 SNPs, respectively), with c. 80% of the loci detected de novo being also found with reference mapping, using the Atlantic salmon Salmo salar genome as a reference. Lower mapping density but larger nucleotide diversity (π) was generally observed near extremities of linkage groups, consistent with regions of residual tetrasomic inheritance observed in salmonids. Genome-wide diversity estimates revealed reduced polymorphism in hatchery strains (π = 0·0040 and π = 0·0029 in Atlantic and Mediterranean strains, respectively) compared to wild populations (π = 0·0049), a pattern that was congruent with allelic richness estimated from microsatellite markers. Finally, pronounced heterozygote deficiency was found in hatchery strains (Atlantic F_IS = 0·18; Mediterranean F_IS = 0·42), indicating that stocking practices may affect the genetic diversity in wild populations. These new genomic resources will provide important tools to define better conservation strategies in S. trutta.

Having a BLAST: Searchable transcriptome resources for the gilthead sea bream and the European sea bass

The gilthead sea bream (Sparus aurata) and the European sea bass (Dicentrarchus labrax) are the most important aquaculture species in the Mediterranean Sea and since the last decade it has been seen an exponential increase in their available molecular resources. In order to improve accessibility to transcriptome resources, Expressed Sequence Tags (ESTs), mRNA sequences and raw read sequences were assembled and deposited in BLAST queryable databases. The publicly available sea bream and sea bass sequences (6.4 and 247.5 million) generated 45,094 and 68,117 assembled sequences, with, respectively, arithmetic mean size of 998 and 2125bp and N50 of 1302 and 2966bp. The assemblies will be regularly updated and new analytical tools added to the web server at http://sea.ccmar.ualg.pt.

Repeatability of cortisol stress response in the European sea bass (Dicentrarchus labrax) and transcription differences between individuals with divergent responses

Understanding the stress responses of organisms is of importance in the performance and welfare of farmed animals, including fish. Especially fish in aquaculture commonly face stressors, and better knowledge of their responses may assist in proper husbandry and selection of breeding stocks. European sea bass (Dicentrarchus labrax), a species with high cortisol concentrations, is of major importance in this respect. The main objectives of the present study were to assess the repeatability and consistency of cortisol stress response and to identify differences in liver transcription profiles of European sea bass individuals, showing a consistent low (LR) or high (HR) cortisol response. The progeny of six full sib families was used, and sampled for plasma cortisol after an acute stress challenge once per month, for four consecutive months. Results suggest that cortisol responsiveness was a repeatable trait with LR and HR fish showing low or high resting, free and post-stress cortisol concentrations respectively. Finally, the liver transcription profiles of LR and HR fish showed some important differences, indicating differential hepatic regulation between these divergent phenotypes. These transcription differences were related to various metabolic and immunological processes, with 169 transcripts being transcribed exclusively in LR fish and 161 exclusively in HR fish.

Adaptive evolution and segregating load contribute to the genomic landscape of divergence in two tree species connected by episodic gene flow

Speciation often involves repeated episodes of genetic contact between divergent populations before reproductive isolation (RI) is complete. Whole-genome sequencing (WGS) holds great promise for unravelling the genomic bases of speciation. We have studied two ecologically divergent, hybridizing species of the 'model tree' genus Populus (poplars, aspens, cottonwoods), Populus alba and P. tremula, using >8.6 million single nucleotide polymorphisms (SNPs) from WGS of population pools. We used the genomic data to (i) scan these species' genomes for regions of elevated and reduced divergence, (ii) assess key aspects of their joint demographic history based on genomewide site frequency spectra (SFS) and (iii) infer the potential roles of adaptive and deleterious coding mutations in shaping the genomic landscape of divergence. We identified numerous small, unevenly distributed genome regions without fixed polymorphisms despite high overall genomic differentiation. The joint SFS was best explained by ancient and repeated gene flow and allowed pinpointing candidate interspecific migrant tracts. The direction of selection (DoS) differed between genes in putative migrant tracts and the remainder of the genome, thus indicating the potential roles of adaptive divergence and segregating deleterious mutations on the evolution and breakdown of RI. Genes affected by positive selection during divergence were enriched for several functionally interesting groups, including well-known candidate 'speciation genes' involved in plant innate immunity. Our results suggest that adaptive divergence affects RI in these hybridizing species mainly through intrinsic and demographic processes. Integrating genomic with molecular data holds great promise for revealing the effects of particular genetic pathways on speciation.

Gene Expression Profiling Reveals Functional Specialization along the Intestinal Tract of a Carnivorous Teleostean Fish (Dicentrarchus labrax)

High-quality sequencing reads from the intestine of European sea bass were assembled, annotated by similarity against protein reference databases and combined with nucleotide sequences from public and private databases. After redundancy filtering, 24,906 non-redundant annotated sequences encoding 15,367 different gene descriptions were obtained. These annotated sequences were used to design a custom, high-density oligo-microarray (8 × 15 K) for the transcriptomic profiling of anterior (AI), middle (MI), and posterior (PI) intestinal segments. Similar molecular signatures were found for AI and MI segments, which were combined in a single group (AI-MI) whereas the PI outstood separately, with more than 1900 differentially expressed genes with a fold-change cutoff of 2. Functional analysis revealed that molecular and cellular functions related to feed digestion and nutrient absorption and transport were over-represented in AI-MI segments. By contrast, the initiation and establishment of immune defense mechanisms became especially relevant in PI, although the microarray expression profiling validated by qPCR indicated that these functional changes are gradual from anterior to posterior intestinal segments. This functional divergence occurred in association with spatial transcriptional changes in nutrient transporters and the mucosal chemosensing system via G protein-coupled receptors. These findings contribute to identify key indicators of gut functions and to compare different fish feeding strategies and immune defense mechanisms acquired along the evolution of teleosts.

Effects of Sodium Butyrate Treatment on Histone Modifications and the Expression of Genes Related to Epigenetic Regulatory Mechanisms and Immune Response in European Sea Bass (Dicentrarchus Labrax) Fed a Plant-Based Diet

Bacteria that inhabit the epithelium of the animals' digestive tract provide the essential biochemical pathways for fermenting otherwise indigestible dietary fibers, leading to the production of short-chain fatty acids (SCFAs). Of the major SCFAs, butyrate has received particular attention due to its numerous positive effects on the health of the intestinal tract and peripheral tissues. The mechanisms of action of this four-carbon chain organic acid are different; many of these are related to its potent regulatory effect on gene expression since butyrate is a histone deacetylase inhibitor that play a predominant role in the epigenetic regulation of gene expression and cell function. In the present work, we investigated in the European sea bass (Dicentrarchus labrax) the effects of butyrate used as a feed additive on fish epigenetics as well as its regulatory role in mucosal protection and immune homeostasis through impact on gene expression. Seven target genes related to inflammatory response and reinforcement of the epithelial defense barrier [tnfα (tumor necrosis factor alpha) il1β, (interleukin 1beta), il-6, il-8, il-10, and muc2 (mucin 2)] and five target genes related to epigenetic modifications [dicer1(double-stranded RNA-specific endoribonuclease), ehmt2 (euchromatic histone-lysine-N-methyltransferase 2), pcgf2 (polycomb group ring finger 2), hdac11 (histone deacetylase-11), and jarid2a (jumonji)] were analyzed in fish intestine and liver. We also investigated the effect of dietary butyrate supplementation on histone acetylation, by performing an immunoblotting analysis on liver core histone extracts. Results of the eight-week-long feeding trial showed no significant differences in weight gain or SGR (specific growth rate) of sea bass that received 0.2% sodium butyrate supplementation in the diet in comparison to control fish that received a diet without Na-butyrate. Dietary butyrate led to a twofold increase in the acetylation level of histone H4 at lysine 8, but showed no effect on the histone H3 at Lys9. Moreover, two different isoforms of histone H3 that might correspond to the H3.1 and H3.2 isoforms previously found in terrestrial animals were separated on the immunoblots. The expression of four (il1 β, il8, irf1, and tnfα) out of seven analyzed genes related to mucosal protection and inflammatory response was significantly different between the two analyzed tissues but only il10 showed differences in expression due to the interaction between tissue and butyrate treatment. In addition, butyrate caused significant changes in vivo in the expression of genes related to epigenetic regulatory mechanisms such as hdac11, ehmt2, and dicer1. Statistical analysis by two-way ANOVA for these genes showed not only significant differences due to the butyrate treatment, but also due to the interaction between tissue and treatment.

Inferring the demographic history underlying parallel genomic divergence among pairs of parasitic and nonparasitic lamprey ecotypes

Understanding the evolutionary mechanisms generating parallel genomic divergence patterns among replicate ecotype pairs remains an important challenge in speciation research. We investigated the genomic divergence between the anadromous parasitic river lamprey (Lampetra fluviatilis) and the freshwater-resident nonparasitic brook lamprey (Lampetra planeri) in nine population pairs displaying variable levels of geographic connectivity. We genotyped 338 individuals with RAD sequencing and inferred the demographic divergence history of each population pair using a diffusion approximation method. Divergence patterns in geographically connected population pairs were better explained by introgression after secondary contact, whereas disconnected population pairs have retained a signal of ancient migration. In all ecotype pairs, models accounting for differential introgression among loci outperformed homogeneous migration models. Generating neutral predictions from the inferred divergence scenarios to detect highly differentiated markers identified greater proportions of outliers in disconnected population pairs than in connected pairs. However, increased similarity in the most divergent genomic regions was found among connected ecotype pairs, indicating that gene flow was instrumental in generating parallelism at the molecular level. These results suggest that heterogeneous genomic differentiation and parallelism among replicate ecotype pairs have partly emerged through restricted introgression in genomic islands.

Genetic variability, local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South-East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South-East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. FST and principal components analysis revealed striking divergence between South-East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South-East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South-East Asia during mid-Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the 'genomic islands' scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid-Pleistocene.

Genomic heterogeneity of historical gene flow between two species of newts inferred from transcriptome data

The role of gene flow in species formation is a major unresolved issue in speciation biology. Progress in this area requires information on the long-term patterns of gene flow between diverging species. Here, we used thousands of single-nucleotide polymorphisms derived from transcriptome resequencing and a method modeling the joint frequency spectrum of these polymorphisms to reconstruct patterns of historical gene flow between two Lissotriton newts: L. vulgaris (Lv) and L. montandoni (Lm). We tested several models of divergence including complete isolation and various scenarios of historical gene flow. The model of secondary contact received the highest support. According to this model, the species split from their common ancestor ca. 5.5 million years (MY) ago, evolved in isolation for ca. 2 MY, and have been exchanging genes for the last 3.5 MY Demographic changes have been inferred in both species, with the current effective population size of ca. 0.7 million in Lv and 0.2 million in Lm. The postdivergence gene flow resulted in two-directional introgression which affected the genomes of both species, but was more pronounced from Lv to Lm. Interestingly, we found evidence for genomic heterogeneity of interspecific gene flow. This study demonstrates the complexity of long-term gene flow between distinct but incompletely reproductively isolated taxa which divergence was initiated millions of years ago.

The complete mitochondrial genome sequence of white perch Morone americana (Perciformes, Moronidae)

The complete mitochondrial genome of the white perch Morone americana is described in this study. The total length of M. americana mitogenome is 17,966 bp, consisting of 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNA genes and a noncoding control region. As with several other species in Moronidae, the ND6 gene in M. americana is found within the control region rather than at the canonical position between the ND5 and CYTB genes. Phylogenetic analysis based on CYTB gene places M. americana within Moronidae and confirms its close relationship with yellow perch (M. mississippiensis). This study provides a basic resource for further research on genetic structure and demographic history of M. americana.

Parallel genetic divergence among coastal-marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome-wide data set of RAD-derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal-marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine-scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.

Drospirenone intake alters plasmatic steroid levels and cyp17a1 expression in gonads of juvenile sea bass

Drospirenone (DRO) is one of the most widely used progestins in contraceptive treatments and hormone replacement therapies. The pharmacokinetics and potential toxicological effects of DRO were investigated in juvenile sea bass (Dicentrarchus labrax) exposed through the diet (0.01-10 μg DRO/g) for up to 31 days. DRO was detected in the blood (4-27 ng/mL) of fish exposed to the highest concentration, with no significant bioaccumulation over time and no alteration of hepatic metabolizing enzymes, namely, CYP1A and CYP3A-catalysed activities and UDP-glucuronyltransferase (UGT). Pregnenolone (P5), progesterone (P4), 17α-hydroxyprogesterone (17P4), 17α-hydroxypregnenolone (17P5), androstenedione (AD) and testosterone (T) were determined in plasma and gene expression of cyp17a1, cyp19a1a and cyp11β analysed by qRT-PCR in gonads. The significant increase in plasmatic levels of 17P5, 17P4 and AD detected after 31 days exposure to 10 ng DRO/g together with the increased expression of cyp17a1 in females evidence the ability of DRO to alter steroid synthesis at low intake concentrations (7 ng DRO/day). However, the potential consequences of this steroid shift for female reproduction remain to be investigated.

Conserved Anti-Müllerian Hormone: Anti-Müllerian Hormone Type-2 Receptor Specific Interaction and Intracellular Signaling in Teleosts

In higher vertebrates, anti-Müllerian hormone (AMH) is required for Müllerian duct regression in fetal males. AMH is also produced during postnatal life in both sexes regulating steroidogenesis and early stages of folliculogenesis. Teleosts lack Müllerian ducts, but Amh has been identified in several species including European sea bass. However, information on Amh type-2 receptor (Amhr2), the specific receptor for Amh binding, is restricted to a couple of fish species. Here, we report on cloning sea bass amhr2, the production of a recombinant sea bass Amh, and the functional analysis of this ligand-receptor couple. Phylogenetic analysis revealed that sea bass amhr2 segregates with Amhr2 from other vertebrates. This piscine receptor is capable of activating Smad proteins. Antibodies raised against sea bass Amh were used to study native and recombinant Amh, revealing proteins in the range of 66-70 kDa corresponding to the full length Amh. Once proteolytically treated, recombinant sea bass Amh generates a 12 kDa C-terminal mature protein, suggesting that contrary to what has been described for other fish Amh proteins, this protein is processed in a similar way as mammalian AMH. The mature sea bass Amh is a biologically active protein able to bind sea bass Amhr2 and, surprisingly, also human AMHR2. In prepubertal sea bass testes, Amh was detected by immunohistochemistry mostly in Sertoli cells surrounding early germ-cell generations. During spermatogenesis, a weaker staining signal could be observed in Sertoli cells surrounding spermatocytes.