Viral nervous necrosis resistance in gilthead sea bream (Sparus aurata) at the larval stage: heritability and accuracy of genomic prediction with different training and testing settings.

BACKGROUND: The gilthead sea bream (Sparus aurata) has long been considered resistant to viral nervous necrosis (VNN), until recently, when significant mortalities caused by a reassortant nervous necrosis virus (NNV) strain were reported. Selective breeding to enhance resistance against NNV might be a preventive action. In this study, 972 sea bream larvae were subjected to a NNV challenge test and the symptomatology was recorded. All the experimental fish and their parents were genotyped using a genome-wide single nucleotide polymorphism (SNP) array consisting of over 26,000 markers. RESULTS: Estimates of pedigree-based and genomic heritabilities of VNN symptomatology were consistent with each other (0.21, highest posterior density interval at 95% (HPD95%): 0.1-0.4; 0.19, HPD95%: 0.1-0.3, respectively). The genome-wide association study suggested one genomic region, i.e., in linkage group (LG) 23 that might be involved in sea bream VNN resistance, although it was far from the genome-wide significance threshold. The accuracies (r) of the predicted estimated breeding values (EBV) provided by three Bayesian genomic regression models (Bayes B, Bayes C, and Ridge Regression) were consistent and on average were equal to 0.90 when assessed in a set of cross-validation (CV) procedures. When genomic relationships between training and testing sets were minimized, accuracy decreased greatly (r = 0.53 for a validation based on genomic clustering, r = 0.12 for a validation based on a leave-one-family-out approach focused on the parents of the challenged fish). Classification of the phenotype using the genomic predictions of the phenotype or using the genomic predictions of the pedigree-based, all data included, EBV as classifiers was moderately accurate (area under the ROC curve 0.60 and 0.66, respectively). CONCLUSIONS: The estimate of the heritability for VNN symptomatology indicates that it is feasible to implement selective breeding programs for increased resistance to VNN of sea bream larvae/juveniles. Exploiting genomic information offers the opportunity of developing prediction tools for VNN resistance, and genomic models can be trained on EBV using all data or phenotypes, with minimal differences in classification performance of the trait phenotype. In a long-term view, the weakening of the genomic ties between animals in the training and test sets leads to decreased genomic prediction accuracies, thus periodical update of the reference population with new data is mandatory.

Resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax L.): heritability and relationships with body weight, cortisol concentration, and antibody titer.

BACKGROUND: Susceptibility of European sea bass (Dicentrarchus labrax L.) to viral nervous necrosis (VNN) is well-known. Interest towards selective breeding as a tool to enhance genetic resistance in this species has increased sharply due to the major threat represented by VNN for farmed sea bass and limitations concerning specific therapeutical measures. A sea bass experimental population (N = 650) was challenged with nervous necrosis virus (NNV) to investigate genetic variation in VNN mortality. In addition, relationships of this trait with serum cortisol concentration after stress exposure, antibody titer against NNV antigens, and body weight at a fixed age were studied to identify potential indicator traits of VNN resistance. RESULTS: The estimate of heritability for VNN mortality was moderate and ranged from 0.15 (HPD95%, 95% highest posterior density interval: 0.02, 0.31) to 0.23 (HPD95%: 0.06, 0.47). Heritability estimates for cortisol concentration, antibody titer, and body weight were 0.19 (HPD95%: 0.07, 0.34), 0.36 (HPD95%: 0.16, 0.59) and 0.57 (HPD95%: 0.33, 0.84), respectively. Phenotypic relationships between traits were trivial and not statistically significant, except for the estimated correlation between antibody titer and body weight (0.24). Genetic correlations of mortality with body weight or antibody titer (- 0.39) exhibited a 0.89 probability of being negative. A negligible genetic correlation between mortality and cortisol concentration was detected. Antibody titer was estimated to be positively correlated with body weight (0.49). CONCLUSIONS: Antibody titer against NNV offers the opportunity to use indirect selection to enhance resistance, while the use of cortisol concentration as an indicator trait in breeding programs for VNN resistance is questionable. The estimate of heritability for VNN mortality indicates the feasibility of selective breeding to enhance resistance to NNV and raises attention to the development of genomic prediction tools to simplify testing procedures for selection candidates.

Depuration processes affect the Vibrio community in the microbiota of the Manila clam, Ruditapes philippinarum.

As filter-feeders, bivalve molluscs accumulate Vibrio into edible tissues. Consequently, an accurate assessment of depuration procedures and the characterization of the persistent Vibrio community in depurated shellfish represent a key issue to guarantee food safety in shellfish products. The present study investigated changes in the natural Vibrio community composition of the Ruditapes philippinarum microbiota with specific focus on human pathogenic species. For this purpose, the study proposed a MLSA-NGS approach (rRNA 16S, recA and pyrH) for the detection and identification of Vibrio species. Clam microbiota were analysed before and after depuration procedures performed in four depuration plants, using culture-dependent and independent approaches. Microbiological counts and NGS data revealed differences in terms of both contamination load and Vibrio community between depuration plants. The novel MLSA-NGS approach allowed for a clear definition of the Vibrio species specific to each depuration plant. Specifically, depurated clam microbiota showed presence of human pathogenic species. Ozone treatments and the density of clams in the depuration tank probably influenced the level of contamination and the Vibrio community composition. The composition of Vibrio community specific to each plant should be carefully evaluated during the risk assessment to guarantee a food-safe shellfish-product for the consumer.

Molecular insights into post-mating immune response in a fish with internal fertilization.

The tight connection between immunity and reproduction has been studied for decades. However, basic knowledge at the molecular level of the effect of mating on immune function is still lacking in many taxa. Determining whether and how the immune system is engaged after mating is a crucial step in understanding post-mating mechanisms of reproduction and sexual selection. Here, we study the transcriptional changes in immunity-related genes caused by the ejaculate in the female reproductive tract using a model species for sexual selection studies, the guppy Poecilia reticulata. To study changes triggered by the ejaculate only, rather than caused by mating, we used artificial inseminations to transfer ejaculate into females. We then compared gene expression in the reproductive tract (gonoduct and ovary) of females artificially inseminated either with ejaculate or with a control solution, after 1 hr and after 6 hr. Overall, contact with ejaculate caused short-term changes in the expression of immune-related genes in the female reproductive tract, with a complex pattern of up- and down-regulation of immune-related pathways, but with clear indication of a marked down-regulation of the immune system shortly after ejaculate contact. This suggests a link between immune function and processes occurring between mating and fertilization in this species.

Gene expression profile analysis of Manila clam (Ruditapes philippinarum) hemocytes after a Vibrio alginolyticus challenge using an immune-enriched oligo-microarray.

BACKGROUND: The Manila clam (Ruditapes philippinarum) is a cultured bivalve with worldwide commercial importance, and diseases cause high economic losses. For this reason, interest in the immune genes in this species has recently increased. The present work describes the construction of the first R. philippinarum microarray containing immune-related hemocyte sequences and its application to study the gene transcription profiles of hemocytes from clams infected with V. alginolyticus through a time course. RESULTS: The complete set of sequences from R. philippinarum available in the public databases and the hemocyte sequences enriched in immune transcripts were assembled successfully. A total of 12,156 annotated sequences were used to construct the 8 × 15 k oligo-microarray. The microarray experiments yielded a total of 579 differentially expressed transcripts. Using the gene expression results, the associated Gene Ontology terms and the enrichment analysis, we found different response mechanisms throughout the experiment. Genes related to signaling, transcription and apoptosis, such as IL-17D, NF-κB or calmodulin, were typically expressed as early as 3 hours post-challenge (hpc), while characteristic immune genes, such as PGRPs, FREPs and defense proteins appeared later at 8 hpc. This immune-triggering response could have affected a high number of processes that seemed to be activated 24 hpc to overcome the Vibrio challenge, including the expression of many cytoskeleton molecules, which is indicative of the active movement of hemocytes. In fact functional studies showed an increment in apoptosis, necrosis or cell migration after the infection. Finally, 72 hpc, activity returned to normal levels, and more than 50% of the genes were downregulated in a negative feedback of all of the previously active processes. CONCLUSIONS: Using a new version of the R. philippinarum oligo-microarray, a putative timing for the response against a Vibrio infection was established. The key point to overcome the challenge seemed to be 8 hours after the challenge, when we detected immune functions that could lead to the destruction of the pathogen and the activation of a wide variety of processes related to homeostasis and defense. These results highlight the importance of a fast response in bivalves and the effectiveness of their innate immune system.

Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius).

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome-wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large- and fine-scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F(CT) = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F(CT) range 0.275-0.705) and fine-scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.

Deep transcriptome sequencing of Pecten maximus hemocytes: a genomic resource for bivalve immunology.

Pecten maximus, the king scallop, is a bivalve species with important commercial value for both fisheries and aquaculture, traditionally consumed in several European countries. Major problems in larval rearing, however, still limit hatchery-based seed production. High mortalities during early larval stages, likely related to bacterial pathogens, represent the most relevant bottleneck. To address this issue, understanding host defense mechanisms against microbes is extremely important. In this study next-generation RNA-sequencing was carried on scallop hemocytes. To enrich for immune-related transcripts, cDNA libraries from hemocytes challenged in vivo with inactivated-Vibrio anguillarum and in vitro with pathogen-associated molecular patterns, as well as unchallenged controls, were sequenced yielding 216,444,674 sequence reads. De novo assembly of the scallop hemocyte transcriptome consisted of 73,732 contigs (31% annotated). A total of 934 contigs encoded proteins with a known immune function, grouped into several functional categories. Particular attention was reserved to Toll-like receptors (TLRs), a family of pattern recognition receptors (PRRs) involved in non-self recognition. Through mining the scallop hemocyte transcriptome, at least four TLRs could be identified. The organization of canonical TLR domains demonstrated that single cysteine cluster and multiple cysteine cluster TLRs co-exist in this species. In addition, preliminary data concerning their mRNA level following bacterial challenge suggested that different members of this family could exhibit opposite responses to pathogenic stimuli. Finally, a global analysis of differential expression comparing gene-expression levels in in vitro and in vivo stimulated hemocytes against controls provided evidence on a large set of transcripts involved in the great scallop immune response.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea).

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam.

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

The mitochondrial genome of the ascalaphid owlfly Libelloides macaronius and comparative evolutionary mitochondriomics of neuropterid insects.

BACKGROUND: The insect order Neuroptera encompasses more than 5,700 described species. To date, only three neuropteran mitochondrial genomes have been fully and one partly sequenced. Current knowledge on neuropteran mitochondrial genomes is limited, and new data are strongly required. In the present work, the mitochondrial genome of the ascalaphid owlfly Libelloides macaronius is described and compared with the known neuropterid mitochondrial genomes: Megaloptera, Neuroptera and Raphidioptera. These analyses are further extended to other endopterygotan orders. RESULTS: The mitochondrial genome of L. macaronius is a circular molecule 15,890 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. The gene order of this newly sequenced genome is unique among Neuroptera and differs from the ancestral type of insects in the translocation of trnC. The L. macaronius genome shows the lowest A+T content (74.50%) among known neuropterid genomes. Protein-coding genes possess the typical mitochondrial start codons, except for cox1, which has an unusual ACG. Comparisons among endopterygotan mitochondrial genomes showed that A+T content and AT/GC-skews exhibit a broad range of variation among 84 analyzed taxa. Comparative analyses showed that neuropterid mitochondrial protein-coding genes experienced complex evolutionary histories, involving features ranging from codon usage to rate of substitution, that make them potential markers for population genetics/phylogenetics studies at different taxonomic ranks. The 22 tRNAs show variable substitution patterns in Neuropterida, with higher sequence conservation in genes located on the α strand. Inferred secondary structures for neuropterid rrnS and rrnL genes largely agree with those known for other insects. For the first time, a model is provided for domain I of an insect rrnL. The control region in Neuropterida, as in other insects, is fast-evolving genomic region, characterized by AT-rich motifs. CONCLUSIONS: The new genome shares many features with known neuropteran genomes but differs in its low A+T content. Comparative analysis of neuropterid mitochondrial genes showed that they experienced distinct evolutionary patterns. Both tRNA families and ribosomal RNAs show composite substitution pathways. The neuropterid mitochondrial genome is characterized by a complex evolutionary history.

Combining Culture-Dependent and Culture-Independent Methods: New Methodology Insight on the Vibrio Community of Ruditapes philippinarum.

Vibrios represent a natural contaminant of seafood products. V. alginolyticus, V. cholerae, V. parahaemolyticus and V. vulnificus are the most hazardous species to human health. Given the worldwide consumption of mollusc products, reliable detection of Vibrio species is recommended to prevent human vibriosis. In this study, culture-dependent and -independent methods were compared and integrated to implement knowledge of the Manila clam Vibrio community composition. Here, 16S and recA-pyrH metabarcoding were applied to compare the microbial communities of homogenate clam samples (culture-independent method) and their culture-derived samples plated on three different media (culture-dependent method). In addition, a subset of plated clam samples was investigated using shotgun metagenomics. Homogenate metabarcoding characterized the most abundant taxa (16S) and Vibrio species (recA-pyrH). Culture-dependent metabarcoding detected the cultivable taxa, including rare species. Moreover, marine agar medium was found to be a useful substrate for the recovery of several Vibrio species, including the main human pathogenic ones. The culture-dependent shotgun metagenomics detected all the main human pathogenic Vibrio species and a higher number of vibrios with respect to the recA-pyrH metabarcoding. The study revealed that integration of culture-dependent and culture-independent methods might be a valid approach for the characterization of Vibrio biodiversity.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish.

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

fshr: a fish sex-determining locus shows variable incomplete penetrance across flathead grey mullet populations.

Whole-genome sequencing data were produced from a single flathead grey mullet female and assembled into a draft genome sequence, whereas publicly available sequence data were used to obtain a male draft sequence. Two pools, each consisting of 60 unrelated individuals, respectively, of male and female fish were analyzed using Pool-Sequencing. Mapping and analysis of Pool-Seq data against the draft genome(s) revealed >30 loci potentially associated with sex, the most promising locus of which, encoding the follicle-stimulating hormone receptor (fshr) and harboring two missense variants, was genotyped on 245 fish from four Mediterranean populations. Genotype data showed that fshr represents a previously unknown sex-determining locus, although the incomplete association pattern between fshr genotype and sex-phenotype, the variability of such pattern across different populations, and the presence of other candidate loci reveal that a greater complexity underlies sex determination in the flathead grey mullet.

Long-lasting effects of chronic exposure to chemical pollution on the hologenome of the Manila clam.

Chronic exposure to pollutants affects natural populations, creating specific molecular and biochemical signatures. In the present study, we tested the hypothesis that chronic exposure to pollutants might have substantial effects on the Manila clam hologenome long after removal from contaminated sites. To reach this goal, a highly integrative approach was implemented, combining transcriptome, genetic and microbiota analyses with the evaluation of biochemical and histological profiles of the edible Manila clam Ruditapes philippinarum, as it was transplanted for 6 months from the polluted area of Porto Marghera (PM) to the clean area of Chioggia (Venice lagoon, Italy). One month post-transplantation, PM clams showed several modifications to its resident microbiota, including an overrepresentation of the opportunistic pathogen Arcobacter spp. This may be related to the upregulation of several immune genes in the PM clams, potentially representing a host response to the increased abundance of deleterious bacteria. Six months after transplantation, PM clams demonstrated a lower ability to respond to environmental/physiological stressors related to the summer season, and the hepatopancreas-associated microbiota still showed different compositions among PM and CH clams. This study confirms that different stressors have predictable effects in clams at different biological levels and demonstrates that chronic exposure to pollutants leads to long-lasting effects on the animal hologenome. In addition, no genetic differentiation between samples from the two areas was detected, confirming that PM and CH clams belong to a single population. Overall, the obtained responses were largely reversible and potentially related to phenotypic plasticity rather than genetic adaptation. The results here presented will be functional for the assessment of the environmental risk imposed by chemicals on an economically important bivalve species.

Population structure, demographic history, and selective processes: contrasting evidences from mitochondrial and nuclear markers in the European spiny lobster Palinurus elephas (Fabricius, 1787).

The European spiny lobster Palinurus elephas (Fabricius, 1787) is an ecologically and economically important species inhabiting a wide geographic range that extends from the North-east Atlantic and Azores to the Eastern Mediterranean. We investigated the population structure and evolutionary history of this species by both mitochondrial and microsatellite markers. Ten population samples covering a large part of the species distribution range (three samples from the Atlantic Ocean and seven from the Mediterranean Sea) were analyzed for a portion of the mitochondrial control region and seven polymorphic microsatellite loci. Both markers rejected the hypothesis of panmixia identifying two differentiated gene pools. The control region clearly distinguished the Atlantic and Mediterranean populations in two genetically separated groups. Microsatellites, also revealed two groups roughly associated to the Atlantic-Mediterranean separation, however, the Azores sample did not conform to this geographic scheme. Discrepancy between mitochondrial and nuclear markers emerged also when reconstructing the history of the species. Neutrality tests of the mitochondrial sequences indicated a departure from mutation-drift equilibrium that, combined to the mismatch analysis, pointed toward a sudden population expansion in both Atlantic and Mediterranean gene pools. Unexpectedly, microsatellites did not identify any signal of population expansion neither in the Atlantic pool nor in the Mediterranean one.

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.

Genetic and spatial characterization of the red fox (Vulpes vulpes) population in the area stretching between the Eastern and Dinaric Alps and its relationship with rabies and canine distemper dynamics.

Information on the population dynamics of a reservoir species have been increasingly adopted to understand and eventually predict the dispersal patterns of infectious diseases throughout an area. Although potentially relevant, to date there are no studies which have investigated the genetic structure of the red fox population in relation to infectious disease dynamics. Therefore, we genetically and spatially characterised the red fox population in the area stretching between the Eastern and Dinaric Alps, which has been affected by both distemper and rabies at different time intervals. Red foxes collected from north-eastern Italy, Austria, Slovenia and Croatia between 2006-2012, were studied using a set of 21 microsatellite markers. We confirmed a weak genetic differentiation within the fox population using Bayesian clustering analyses, and we were able to differentiate the fox population into geographically segregated groups. Our finding might be due to the presence of geographical barriers that have likely influenced the distribution of the fox population, limiting in turn gene flow and spread of infectious diseases. Focusing on the Italian red fox population, we observed interesting variations in the prevalence of both diseases among distinct fox clusters, with the previously identified Italy 1 and Italy 2 rabies as well as distemper viruses preferentially affecting different sub-groups identified in the study. Knowledge of the regional-scale population structure can improve understanding of the epidemiology and spread of diseases. Our study paves the way for an integrated approach for disease control coupling pathogen, host and environmental data to inform targeted control programs in the future.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The population genomics of yellowfin tuna (Thunnus albacares) at global geographic scale challenges current stock delineation.

Yellowfin tuna, Thunnus albacares, is one of the most important seafood commodities in the world. Despite its great biological and economic importance, conflicting evidence arises from classical genetic and tagging studies concerning the yellowfin tuna population structure at local and global oceanic scales. Access to more powerful and cost effective genetic tools would represent the first step towards resolving the population structure of yellowfin tuna across its distribution range. Using a panel of 939 neutral Single Nucleotide Polymorphisms (SNPs), and the most comprehensive data set of yellowfin samples available so far, we found genetic differentiation among the Atlantic, Indian and Pacific oceans. The genetic stock structure analysis carried out with 33 outlier SNPs, putatively under selection, identified discrete populations within the Pacific Ocean and, for the first time, also within the Atlantic Ocean. Stock assessment approaches that consider genetic differences at neutral and adaptive genomic loci should be routinely implemented to check the status of the yellowfin tuna, prevent illegal trade, and develop more sustainable management measures.

Genomic analysis of Sparus aurata reveals the evolutionary dynamics of sex-biased genes in a sequential hermaphrodite fish.

Sexual dimorphism is a fascinating subject in evolutionary biology and mostly results from sex-biased expression of genes, which have been shown to evolve faster in gonochoristic species. We report here genome and sex-specific transcriptome sequencing of Sparus aurata, a sequential hermaphrodite fish. Evolutionary comparative analysis reveals that sex-biased genes in S. aurata are similar in number and function, but evolved following strikingly divergent patterns compared with gonochoristic species, showing overall slower rates because of stronger functional constraints. Fast evolution is observed only for highly ovary-biased genes due to female-specific patterns of selection that are related to the peculiar reproduction mode of S. aurata, first maturing as male, then as female. To our knowledge, these findings represent the first genome-wide analysis on sex-biased loci in a hermaphrodite vertebrate species, demonstrating how having two sexes in the same individual profoundly affects the fate of a large set of evolutionarily relevant genes.