The Importance of Baseline Health Surveillance Efforts in Freshwater Fish Conservation Using the Threatened Iberian Leuciscids as an Example.

Freshwater fish species are experiencing the highest decline among vertebrates in this century. Although a great effort has been made to identify and tackle threats to the conservation of this taxa, several knowledge gaps still exist particularly for noncommercial endangered species, including considerations regarding fish health status. These species face deteriorating environmental conditions in their natural habitats that may lead to stress and increased risk for infectious disease outbreaks. Establishing health surveillance is crucial to identify and predict physiologic disruption in fish populations. Additionally, information retrieved may be used to direct targeted efforts to contribute to improving the conservation status of these species. We used threatened Iberian leuciscids as a case study to discuss the current knowledge regarding their health surveillance and to suggest recommendations for the establishment of practical health assessments that can benefit conservation plans for these species and be implemented in threatened or endangered freshwater fish species plans globally.

Three mitochondrial lineages and no Atlantic-Mediterranean barrier for the bogue Boops boops across its widespread distribution.

Marine species exhibiting wide distributional ranges are frequently subdivided into discrete genetic units over limited spatial scales. This is often due to specific life-history traits or oceanographic barriers that prevent gene flow. Fine-scale sampling studies revealed distinct phylogeographic patterns in the northeastern Atlantic and the Mediterranean, ranging from panmixia to noticeable population genetic structure. Here, we used mitochondrial sequence data to analyse connectivity in the bogue Boops boops throughout most of its widespread distribution. Our results identified the existence of three clades, one comprising specimens from the Azores and eastern Atlantic/Mediterranean, another with individuals from the Canary Islands, Madeira and Cape Verde archipelagos, and the third with samples from Mauritania only. One of the branches of the northern subtropical gyre (Azores Current) that drifts towards the Gulf of Cádiz promotes a closer connection between the Azores, southern Portugal and the Mediterranean B. boops populations. The Almería-Oran Front, widely recognised as an oceanographic barrier for many organisms to cross the Atlantic-Mediterranean divide, does not seem to affect the dispersal of this benthopelagic species. The southward movement of the Cape Verde Frontal Zone during the winter, combined with the relatively short duration of the pelagic larval stage of B. boops, may be potential factors for preventing the connectivity between the Atlantic oceanic archipelagos and Mauritania shaping the genetic signature of this species.

Aeromonas spp. Prevalence, Virulence, and Antimicrobial Resistance in an Ex Situ Program for Threatened Freshwater Fish-A Pilot Study with Protective Measures.

Ex situ breeding programs are important conservation tools for endangered freshwater fish. However, developing husbandry techniques that decrease the likelihood of disease, antimicrobial resistance, and virulence determinants acquisition during this process is challenging. In this pilot study, we conducted a captivity experiment with Portuguese nase (Iberochondrostoma lusitanicum), a critically endangered leuciscid species, to investigate the influence of simple protective measures (i.e., material disinfection protocols and animal handling with gloves) on the dynamics of a potential pathogenic genus, Aeromonas, as well as its virulence profiles and antimicrobial resistance signatures. Our findings show that antimicrobial resistance in Aeromonas spp. collected from I. lusitanicum significantly increased during the extent of the assay (5 weeks), with all isolates collected at the end of the study classified as multidrug-resistant. Additionally, humans handling fishes without protective measures were colonized by Aeromonas spp. The use of protective measures suggested a decreasing trend in Aeromonas spp. prevalence in I. lusitanicum, while bacterial isolates displayed significantly lower virulence index values when virulence phenotypical expression was tested at 22 °C. Despite this study representing an initial trial, which needs support from further research, protective measures tested are considered a simple tool to be applied in ex situ breeding programs for aquatic animals worldwide. Furthermore, current results raise concern regarding antimicrobial resistance amplification and zoonotic transmission of Aeromonas spp. in aquatic ex situ programs.

Sympatric threatened Iberian leuciscids exhibit differences in Aeromonas diversity and skin lesions' prevalence.

Assessments regarding health aspects of Iberian leuciscids are limited. There is currently an information gap regarding effects of infectious diseases on these populations and their role as a possible conservation threat. Moreover, differences in susceptibility to particular agents, such as Aeromonas spp., by different species/populations is not clear. To understand potential differences in Aeromonas diversity and load, as well as in the prevalence and proportion of skin lesions, in fishes exposed to similar environmental conditions, an observational study was implemented. Using a set of 12 individuals belonging to two sympatric Iberian leuciscid species (Squalius pyrenaicus and Iberochondrostoma lusitanicum), the skin lesion score in each individual was analyzed. Furthermore, a bacterial collection of Aeromonas spp. isolated from each individual was created and isolates' load was quantified by plate counting, identified at species level using a multiplex-PCR assay and virulence profiles established using classical phenotypic methods. The similarity relationships of the isolates were evaluated using a RAPD analysis. The skin lesion score was significantly higher in S. pyrenaicus, while the Aeromonas spp. load did not differ between species. When analyzing Aeromonas species diversity between fishes, different patterns were observed. A predominance of A. hydrophila was detected in S. pyrenaicus individuals, while I. lusitanicum individuals displayed a more diverse structure. Similarly, the virulence index of isolates from S. pyrenaicus was higher, mostly due to the isolated Aeromonas species. Genomic typing clustered the isolates mainly by fish species and skin lesion score. Specific Aeromonas clusters were associated with higher virulence indexes. Current results suggest potential differences in susceptibility to Aeromonas spp. at the fish species/individual level, and constitute important knowledge for proper wildlife management through the signalization of at-risk fish populations and hierarchization of conservation measures.

Climatic Alterations Influence Bacterial Growth, Biofilm Production and Antimicrobial Resistance Profiles in Aeromonas spp.

Climate change is expected to create environmental disruptions that will impact a wide array of biota. Projections for freshwater ecosystems include severe alterations with gradients across geographical areas. Life traits in bacteria are modulated by environmental parameters, but there is still uncertainty regarding bacterial responses to changes caused by climatic alterations. In this study, we used a river water microcosm model to evaluate how Aeromonas spp., an important pathogenic and zoonotic genus ubiquitary in aquatic ecosystems, responds to environmental variations of temperature and pH as expected by future projections. Namely, we evaluated bacterial growth, biofilm production and antimicrobial resistance profiles of Aeromonas species in pure and mixed cultures. Biofilm production was significantly influenced by temperature and culture, while temperature and pH affected bacterial growth. Reversion of antimicrobial susceptibility status occurred in the majority of strains and tested antimicrobial compounds, with several combinations of temperature and pH contributing to this effect. Current results highlight the consequences that bacterial genus such as Aeromonas will experience with climatic alterations, specifically how their proliferation and virulence and phenotypic resistance expression will be modulated. Such information is fundamental to predict and prevent future outbreaks and deleterious effects that these bacterial species might have in human and animal populations.

Molecular Epidemiology, Virulence Traits and Antimicrobial Resistance Signatures of Aeromonas spp. in the Critically Endangered Iberochondrostoma lusitanicum Follow Geographical and Seasonal Patterns.

Despite the fact that freshwater fish populations are experiencing severe declines worldwide, our knowledge on the interaction between endangered populations and pathogenic agents remains scarce. In this study, we investigated the prevalence and structure of Aeromonas communities isolated from the critically endangered Iberochondrostoma lusitanicum, a model species for threatened Iberian leuciscids, as well as health parameters in this species. Additionally, we evaluated the virulence profiles, antimicrobial resistance signatures and genomic relationships of the Aeromonas isolates. Lesion prevalence, extension and body condition were deeply affected by location and seasonality, with poorer performances in the dry season. Aeromonas composition shifted among seasons and was also different across river streams. The pathogenic potential of the isolates significantly increased during the dry season. Additionally, isolates displaying clinically relevant antimicrobial resistance phenotypes (carbapenem and fluroquinolone resistance) were detected. As it inhabits intermittent rivers, often reduced to disconnected pools during the summer, the dry season is a critical period for I. lusitanicum, with lower general health status and a higher potential of infection by Aeromonas spp. Habitat quality seems a determining factor on the sustainable development of this fish species. Also, these individuals act as reservoirs of important antimicrobial resistant bacteria with potential implications for public health.

Genetic hypervariability of a Northeastern Atlantic venomous rockfish.

BACKGROUND: Understanding the interplay between climate and current and historical factors shaping genetic diversity is pivotal to infer changes in marine species range and communities' composition. A phylogeographical break between the Atlantic and the Mediterranean has been documented for several marine organisms, translating into limited dispersal between the two basins. METHODS: In this study, we screened the intraspecific diversity of 150 individuals of the Madeira rockfish (Scorpaena maderensis) across its distributional range (seven sampling locations in the Atlantic and Mediterranean basins) using the mitochondrial control region and the nuclear S7 first intron. RESULTS: The present work is the most comprehensive study done for this species, yielding no genetic structure across sampled locations and no detectable Atlantic-Mediterranean break in connectivity. Our results reveal deep and hyper-diverse bush-like genealogies with large numbers of singletons and very few shared haplotypes. The genetic hyper-diversity found for the Madeira rockfish is relatively uncommon in rocky coastal species, whose dispersal capability is limited by local oceanographic patterns. The effect of climate warming on the distribution of the species is discussed.

Genetic population structure of the Blackspot seabream (Pagellus bogaraveo): contribution of mtDNA control region to fisheries management.

Marine fisheries management models have traditionally considered biological parameters and geopolitical boundaries. The result is the existence of fisheries management units that do not match genetic populations. However, this panorama is changing with the contribution of genetic and genomic data. Pagellus bogaraveo is a commercially important sparid in the northeast Atlantic, with three stock components being considered by ICES: the Celtic Sea and Bay of Biscay, Atlantic Iberian waters and the Azores. The northern stock collapsed (1975-1985) and is essential to characterize the genetic makeup of the species, particularly in the Iberian Peninsula, where it is managed as a single stock. The mitochondrial control region was used to screen the intraspecific diversity and population structure of individuals from six locations across the species range. The genetic diversity found is similar among sites, and there is differentiation between the Azores and the remaining locations.

Against all odds: a tale of marine range expansion with maintenance of extremely high genetic diversity.

The displacement of species from equatorial latitudes to temperate locations following the increase in sea surface temperatures is among the significant reported consequences of climate change. Shifts in the distributional ranges of species result in fish communities tropicalisation, i.e., high latitude colonisations by typically low latitude distribution species. These movements create new interactions between species and new trophic assemblages. The Senegal seabream, Diplodus bellottii, may be used as a model to understand the population genetics of these invasions. In the last decades, this species has undergone an outstanding range expansion from its African area of origin to the Atlantic coast of the Iberian Peninsula, where now occurs abundantly. Mitochondrial and nuclear markers revealed a striking high haplotypic nucleotide and genetic diversity values, along with significant population differentiation throughout the present-day geographical range of the Senegal seabream. These results are not consistent with the central-marginal hypothesis, nor with the expectations of a leptokurtic distribution of individuals, as D. bellottii seems to be able to retain exceptional levels of diversity in marginal and recently colonised areas. We discuss possible causes for hyperdiversity and lack of geographical structure and subsequent implications for fisheries.

Time matters: genetic composition and evaluation of effective population size in temperate coastal fish species.

BACKGROUND: Extensive knowledge on the genetic characterization of marine organisms has been assembled, mainly concerning the spatial distribution and structuring of populations. Temporal monitoring assesses not only the stability in genetic composition but also its trajectory over time, providing critical information for the accurate forecast of changes in genetic diversity of marine populations, particularly important for both fisheries and endangered species management. We assessed fluctuations in genetic composition among different sampling periods in the western Portuguese shore in three fish species. METHODS: White seabream Diplodus sargus, sand smelt Atherina presbyter and shanny Lipophrys pholis were chosen, because of their genetic patterns in distinct ecological environments, insight into historical and contemporary factors influencing population effective size (N e ), and degree of commercial exploitation. Samples were obtained near Lisbon between 2003 and 2014 and screened for genetic variation with mitochondrial and nuclear markers. Analyses included genealogies, genetic diversities, temporal structures and contemporary N e . RESULTS: For mtDNA no temporal structure was detected, while for nDNA significant differences were recorded between some sampling periods for the shanny and the sand smelt. Haplotype networks revealed deep genealogies, with various levels of diversification. The shanny revealed a smaller N e /generation when compared to the other species, which, in turn, revealed no evidence of genetic drift for most study periods. These results highlight the fact that temporal variations in genetic pool composition should be considered when evaluating the population structure of fish species with long distance dispersal, which are more vulnerable to recruitment fluctuations.

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe.

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (FST ~0.1), over a short coastline (<60 km) in the North Sea-Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post-glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.

Interbreeding between local and translocated populations of a cleaner fish in an experimental mesocosm predicts risk of disrupted local adaptation.

Translocation of organisms within or outside its native range carries the risk of modifying the community of the recipient ecosystems and induces gene flow between locally adapted populations or closely related species. In this study, we evaluated the genetic consequences of large-scale translocation of cleaner wrasses that has become a common practice within the salmon aquaculture industry in northern Europe to combat sea lice infestation. A major concern with this practice is the potential for hybridization of escaped organisms with the local, recipient wrasse population, and thus potentially introduce exogenous alleles and breaking down coadapted gene complexes in local populations. We investigated the potential threat for such genetic introgressions in a large seminatural mesocosm basin. The experimental setting represented a simulated translocation of corkwing wrasse (Symphodus melops) that occurs on a large scale in the Norwegian salmon industry. Parentage assignment analysis of mesocosm's offspring revealed 30% (195 out of 651 offspring) interbreeding between the two populations, despite their being genetically (F ST = 0.094, p < 0.05) and phenotypically differentiated. Moreover, our results suggest that reproductive fitness of the translocated western population doubled that of the local southern population. Our results confirm that human translocations may overcome the impediments imposed by natural habitat discontinuities and urge for immediate action to manage the genetic resources of these small benthic wrasses.

A continuous genome assembly of the corkwing wrasse (Symphodus melops).

The wrasses (Labridae) are one of the most successful and species-rich families of the Perciformes order of teleost fish. Its members display great morphological diversity, and occupy distinct trophic levels in coastal waters and coral reefs. The cleaning behaviour displayed by some wrasses, such as corkwing wrasse (Symphodus melops), is of particular interest for the salmon aquaculture industry to combat and control sea lice infestation as an alternative to chemicals and pharmaceuticals. There are still few genome assemblies available within this fish family for comparative and functional studies, despite the rapid increase in genome resources generated during the past years. Here, we present a highly continuous genome assembly of the corkwing wrasse using PacBio SMRT sequencing (x28.8) followed by error correction with paired-end Illumina data (x132.9). The present genome assembly consists of 5040 contigs (N50 = 461,652 bp) and a total size of 614 Mbp, of which 8.5% of the genome sequence encode known repeated elements. The genome assembly covers 94.21% of highly conserved genes across ray-finned fish species. We find evidence for increased copy numbers specific for corkwing wrasse possibly highlighting diversification and adaptive processes in gene families including N-linked glycosylation (ST8SIA6) and stress response kinases (HIPK1). By comparative analyses, we discover that de novo repeats, often not properly investigated during genome annotation, encode hundreds of immune-related genes. This new genomic resource, together with the ballan wrasse (Labrus bergylta), will allow for in-depth comparative genomics as well as population genetic analyses for the understudied wrasses.

Asymmetrical dispersal and putative isolation-by-distance of an intertidal blenniid across the Atlantic-Mediterranean divide.

Transition zones are of high evolutionary interest because unique patterns of spatial variation are often retained. Here, we investigated the phylogeography of the peacock blenny, Salaria pavo, a small marine intertidal fish that inhabits rocky habitats of the Mediterranean and the adjacent Atlantic Ocean. We screened 170 individuals using mitochondrial and nuclear sequence data from eight locations. Four models of genetic structure were tested: panmixia, isolation-by-distance, secondary contact and phylogeographic break. Results indicated clear asymmetric migration from the Mediterranean to the Atlantic but only marginally supported the isolation-by-distance model. Additionally, the species displays an imprint of demographic expansion compatible with the last glacial maximum. Although the existence of a refugium in the Mediterranean cannot be discarded, the ancestral lineage most likely originated in the Atlantic, where most of the genetic diversity occurs.

Historical gene flow constraints in a northeastern Atlantic fish: phylogeography of the ballan wrasse Labrus bergylta across its distribution range.

The distribution and demographic patterns of marine organisms in the north Atlantic were largely shaped by climatic changes during the Pleistocene, when recurrent glacial maxima forced them to move south or to survive in northern peri-glacial refugia. These patterns were also influenced by biological and ecological factors intrinsic to each species, namely their dispersion ability. The ballan wrasse (Labrus bergylta), the largest labrid fish along Europe's continental margins, is a target for fisheries and aquaculture industry. The phylogeographic pattern, population structure, potential glacial refugia and recolonization routes for this species were assessed across its full distribution range, using mitochondrial and nuclear markers. The existence of a marked population structure can reflect both recolonization from three distinct glacial refugia and current and past oceanographic circulation patterns. Although isolated in present times, shared haplotypes between continental and Azores populations and historical exchange of migrants in both directions point to a common origin of L. bergylta. This situation is likely to be maintained and/or accentuated by current circulation patterns in the north Atlantic, and may lead to incipient speciation in the already distinct Azorean population. Future monitoring of this species is crucial to evaluate how this species is coping with current environmental changes.

Broad-scale sampling of primary freshwater fish populations reveals the role of intrinsic traits, inter-basin connectivity, drainage area and latitude on shaping contemporary patterns of genetic diversity.

Background. Worldwide predictions suggest that up to 75% of the freshwater fish species occurring in rivers with reduced discharge could be extinct by 2070 due to the combined effect of climate change and water abstraction. The Mediterranean region is considered to be a hotspot of freshwater fish diversity but also one of the regions where the effects of climate change will be more severe. Iberian cyprinids are currently highly endangered, with over 68% of the species raising some level of conservation concern. Methods. During the FISHATLAS project, the Portuguese hydrographical network was extensively covered (all the 34 river basins and 47 sub-basins) in order to contribute with valuable data on the genetic diversity distribution patterns of native cyprinid species. A total of 188 populations belonging to 16 cyprinid species of Squalius, Luciobarbus, Achondrostoma, Iberochondrostoma, Anaecypris and Pseudochondrostoma were characterized, for a total of 3,678 cytochrome b gene sequences. Results. When the genetic diversity of these populations was mapped, it highlighted differences among populations from the same species and between species with identical distribution areas. Factors shaping the contemporary patterns of genetic diversity were explored and the results revealed the role of latitude, inter-basin connectivity, migratory behaviour, species maximum size, species range and other species intrinsic traits in determining the genetic diversity of sampled populations. Contrastingly, drainage area and hydrological regime (permanent vs. temporary) seem to have no significant effect on genetic diversity. Species intrinsic traits, maximum size attained, inter-basin connectivity and latitude explained over 30% of the haplotype diversity variance and, generally, the levels of diversity were significantly higher for smaller sized species, from connected and southerly river basins. Discussion. Targeting multiple co-distributed species of primary freshwater fish allowed us to assess the relative role of historical versus contemporary factors affecting genetic diversity. Since different patterns were detected for species with identical distribution areas we postulate that contemporary determinants of genetic diversity (species' intrinsic traits and landscape features) must have played a more significant role than historical factors. Implications for conservation in a context of climate change and highly disturbed habitats are detailed, namely the need to focus management and conservation actions on intraspecific genetic data and to frequently conduct combined genetic and demographic surveys.

Not so sluggish: the success of the Felimare picta complex (Gastropoda, Nudibranchia) crossing Atlantic biogeographic barriers.

The molecular phylogeny of the Atlanto-Mediterranean species of the genus Felimare, particularly those attributed to the species F. picta, was inferred using two mitochondrial markers (16S and COI). A recent revision of the Chromodorididae clarified the taxonomic relationships at the family level redefining the genus Felimare. However, conflicting taxonomic classifications have been proposed for a restrict group of taxa with overlapping morphological characteristics and geographical distributions designated here as the Felimare picta complex. Three major groups were identified: one Mediterranean and amphi-Atlantic group; a western Atlantic group and a tropical eastern Atlantic group. F. picta forms a paraphyletic group since some subspecies are more closely related with taxa traditionaly classified as independent species (e.g. F. zebra) than with other subspecies with allopatric distributions (e.g. F. picta picta and F. picta tema). Usually, nudibranchs have adhesive demersal eggs, short planktonic larval phases and low mobility as adults unless rafting on floating materials occurs. Surprisingly however, the phylogeny of the F. picta complex suggests that they successfully cross main Atlantic biogeographic barriers including the mid-Atlantic barrier. This ability to cross different biogeographic barriers may be related to F. picta's distinct life history and ecological traits. Compared to other Chromodorididae F. picta has larger eggs and planktotrophic larvae which could be related to a longer planktonic phase.

Absence of consistent genetic differentiation among several morphs of Actinia (Actiniaria: Actiniidae) occurring in the Portuguese coast.

Actinia equina, the beadlet sea anemone, is a very labile species, displaying variable colour patterns, broad habitat choice and diverse modes of reproduction. Historically, studies using genetic markers such as allozymes and differences in habitat choice lead several authors to propose that different colour morphs could represent different species. One of the species defined was A. fragacea. In this paper, the relationships between brown, red and green colour morphs of A. equina and A. fragacea were studied, using two DNA fragments (one mitochondrial and one nuclear). Individuals were sampled from three different areas in Portugal separated by a maximum distance of 500 km. This is the first study applying direct sequencing of selected gene fragments to approach the validity of Actinia morphs as different genetic entities. The results show that, at least in the Portuguese coast, these colour morphs do not correspond to the two valid species recognized in the literature. The existence of cryptic species is discussed. 

Metapopulations in temporary streams - the role of drought-flood cycles in promoting high genetic diversity in a critically endangered freshwater fish and its consequences for the future.

Genetic factors have direct and indirect impacts in the viability of endangered species. Assessing their genetic diversity levels and population structure is thus fundamental for conservation and management. In this paper we use mitochondrial and nuclear markers to address phylogeographic and demographic data on the critically endangered Anaecypris hispanica, using a broad sampling set which covered its known distribution area in the Iberian Peninsula. Our results showed that the populations of A. hispanica are strongly differentiated (high and significant ФST and FST values, corroborated by the results from AMOVA and SAMOVA) and genetically diversified. We suggest that the restricted gene flow between populations may have been potentiated by ecological, hydrological and anthropogenic causes. Bayesian skyline plots revealed a signal for expansion for all populations (tMRCA between 68kya and 1.33Mya) and a genetic diversity latitudinal gradient was detected between the populations from the Upper (more diversified) and the Lower (less diversified) Guadiana river basin. We postulate a Pleistocenic westwards colonization route for A. hispanica in the Guadiana river basin, which is in agreement with the tempo and mode of paleoevolution of this drainage. The colonization of River Guadalquivir around 60kya with migrants from the Upper Guadiana, most likely by stream capture, is also suggested. This study highlights the view that critically endangered species facing range retreats (about 47% of its known populations have disappeared in the last 15years) are not necessarily small and genetically depleted. However, the extinction risk is not negligible since A. hispanica faces the combined effect of several deterministic and stochastic negative factors and, moreover, recolonization events after localized extinctions are very unlikely to occur due to the strong isolation of populations and to the patchily ecologically-conditioned distribution of fish. The inferred species distribution models highlight the significant contribution of temperature seasonality and isothermality to A. hispanica occurrence in Guadiana environments and emphasize the importance of stable climatic conditions for the preservation of this species. Given the strong population structure, high percentage of private haplotypes and virtual absence of inter-basin gene flow we suggest that each A. hispanica population should be considered as an independent Operational Conservation Unit and that ex-situ and in-situ actions should be conducted in parallel to allow for the long-term survival of the species and the preservation of the genetic integrity of its populations.

Efficient isolation of polymorphic microsatellites from high-throughput sequence data based on number of repeats.

Transcriptome data are a good resource to develop microsatellites due to their potential in targeting candidate genes. However, developing microsatellites can be a time-consuming enterprise due to the numerous primer pairs to be tested. Therefore, the use of methodologies that make it efficient to identify polymorphic microsatellites is desirable. Here we used a 62,038 contigs transcriptome assembly, obtained from pyrosequencing a peacock blenny (Salaria pavo) multi-tissue cDNA library, to mine for microsatellites and in silico evaluation of their polymorphism. A total of 4190 microsatellites were identified in 3670 unique unigenes, and from these microsatellites, in silico polymorphism was detected in 733. We selected microsatellites based either on their in silico polymorphism and annotation results or based only on their number of repeats. Using these two approaches, 28 microsatellites were successfully amplified in twenty-six individuals, and all but 2 were found to be polymorphic, being the first genetic markers for this species. Our results showed that the strategy of selection based on number of repeats is more efficient in obtaining polymorphic microsatellites than the strategy of in silico polymorphism (allelic richness was 8.2±3.85 and 4.56±2.45 respectively). This study demonstrates that combining the knowledge of number of repeats with other predictors of variability, for example in silico microsatellite polymorphism, improves the rates of polymorphism, yielding microsatellites with higher allelic richness, and decreases the number of monomorphic microsatellites obtained.