In silico mining and characterization of simple sequence repeats from gilthead sea bream (Sparus aurata) expressed sequence tags (EST-SSRs); PCR amplification, polymorphism evaluation and multiplexing and cross-species assays

The genetic population structure and temporal genetic stability of gilthead sea bream Sparus aurata populations in the Aegean and Ionian Seas, using microsatellite DNA markers

We examined 662 gilthead sea bream Sparus aurata from wild samples of the species in the Aegean and Ionian Seas, using 20 EST-linked microsatellite markers, in three multiplex panels, as well as seven anonymous loci. Most of the markers were revealed to be highly polymorphic. We found low genetic differentiation between the sampling stations/areas with total F_ST 0.002 (P < 0.05). Based on comparison of five temporal samples, our results indicate genetic data consistency over time for all tested samples, pointing to stable populations, despite reported repeated escape events. Our results confirm the genetic population structure previously observed in these specific areas, using by far more markers than in previous studies in both coding and non-coding DNA loci. The limited genetic structure and the temporal genetic stability indicate neither major genetic differentiation of local populations by geographic isolation nor influence from anthropogenic factors. These results provide a baseline for future reference in any management programme of both wild and farmed population of S. aurata as well as of other aquaculture species with a potential introgression among farmed and wild populations.

Absence of spatial genetic structure in common dentex (Dentex dentex Linnaeus, 1758) in the Mediterranean Sea as evidenced by nuclear and mitochondrial molecular markers

The common dentex, Dentex dentex, is a fish species which inhabits marine environments in the Mediterranean and Northeast Atlantic regions. This is an important species from an ecological, economic and conservation perspective, however critical information on its population genetic structure is lacking. Most samples were obtained from the Mediterranean Sea (17 sites) with an emphasis around Corsica (5 sites), plus one Atlantic Ocean site. This provided an opportunity to examine genetic structuring at local and broader scales to provide science based data for the management of fishing stocks in the region. Two mitochondrial regions were examined (D-loop and COI) along with eight microsatellite loci. The COI data was combined with publicly available sequences and demonstrated past misidentification of common dentex. All markers indicated the absence of population genetic structure from the Bay of Biscay to the eastern Mediterranean Sea. Bayesian approaches, as well as the statistical tests performed on the allelic frequencies from microsatellite loci, indicated low differentiation between samples; there was only a slight (p = 0.05) indication of isolation by distance. Common dentex is a marine fish species with a unique panmictic population in the Mediterranean and likely in the Atlantic Ocean as well.

Transcriptome-Derived Tetranucleotide Microsatellites and Their Associated Genes from the Giant Panda (Ailuropoda melanoleuca)

Recently, an increasing number of microsatellites or simple sequence repeats (SSRs) have been found and characterized from transcriptomes. Such SSRs can be employed as putative functional markers to easily tag corresponding genes, which play an important role in biomedical studies and genetic analysis. However, the transcriptome-derived SSRs for giant panda (Ailuropoda melanoleuca) are not yet available. In this work, we identified and characterized 20 tetranucleotide microsatellite loci from a transcript database generated from the blood of giant panda. Furthermore, we assigned their predicted transcriptome locations: 16 loci were assigned to untranslated regions (UTRs) and 4 loci were assigned to coding regions (CDSs). Gene identities of 14 transcripts contained corresponding microsatellites were determined, which provide useful information to study the potential contribution of SSRs to gene regulation in giant panda. The polymorphic information content (PIC) values ranged from 0.293 to 0.789 with an average of 0.603 for the 16 UTRs-derived SSRs. Interestingly, 4 CDS-derived microsatellites developed in our study were also polymorphic, and the instability of these 4 CDS-derived SSRs was further validated by re-genotyping and sequencing. The genes containing these 4 CDS-derived SSRs were embedded with various types of repeat motifs. The interaction of all the length-changing SSRs might provide a way against coding region frameshift caused by microsatellite instability. We hope these newly gene-associated biomarkers will pave the way for genetic and biomedical studies for giant panda in the future. In sum, this set of transcriptome-derived markers complements the genetic resources available for giant panda.

Microsatellites in Pursuit of Microbial Genome Evolution

Microsatellites or short sequence repeats are widespread genetic markers which are hypermutable 1-6 bp long short nucleotide motifs. Significantly, their applications in genetics are extensive due to their ceaseless mutational degree, widespread length variations and hypermutability skills. These features make them useful in determining the driving forces of evolution by using powerful molecular techniques. Consequently, revealing important questions, for example, what is the significance of these abundant sequences in DNA, what are their roles in genomic evolution? The answers of these important questions are hidden in the ways these short motifs contributed in altering the microbial genomes since the origin of life. Even though their size ranges from 1 -to- 6 bases, these repeats are becoming one of the most popular genetic probes in determining their associations and phylogenetic relationships in closely related genomes. Currently, they have been widely used in molecular genetics, biotechnology and evolutionary biology. However, due to limited knowledge; there is a significant gap in research and lack of information concerning hypermutational mechanisms. These mechanisms play a key role in microsatellite loci point mutations and phase variations. This review will extend the understandings of impacts and contributions of microsatellite in genomic evolution and their universal applications in microbiology.

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.

Development of the first standardised panel of two new microsatellite multiplex PCRs for gilthead seabream (Sparus aurata L.)

The high number of multiplex PCRs developed for gilthead seabream (Sparus aurata L.) from many different microsatellite markers does not allow comparison among populations. This highlights the need for developing a reproducible panel of markers, which can be used with safety and reliability by all users. In this study, the first standardised panel of two new microsatellite multiplex PCRs was developed for this species. Primers of 138 specific microsatellites from the genetic linkage map were redesigned and evaluated according to their genetic variability, allele size range and genotyping reliability. A protocol to identify and classify genotyping errors or potential errors was proposed to assess the reliability of each marker. Two new multiplex PCRs from the best assessed markers were designed with 11 markers in each, named SMsa1 and SMsa2 (SuperMultiplex Sparus aurata). Three broodstocks (59, 47 and 98 breeders) from different Spanish companies, and a sample of 80 offspring from each one, were analysed to validate the usefulness of these multiplexes in the parental assignation. It was possible to assign each offspring to a single parent pair (100% success) using the exclusion method with SMsa1 and/or SMsa2. In each genotyped a reference sample (Ref-sa) was used, and its DNA is available on request similar to the kits of bin set to genotype by genemapper (v.3.7) software (kit-SMsa1 and kit-SMsa2). This will be a robust and effective tool for pedigree analysis or characterisation of populations and will be proposed as an international panel for this species.

Use of comparative genomics to develop EST-SSRs for red drum (Sciaenops ocellatus)

Microsatellites physically linked to expressed sequence tags (EST-SSRs) are an important resource for linkage mapping and comparative genomics, and data mining in publicly available EST databases is a common strategy for EST-SSR discovery. At present, many species lack species-specific EST sequence data needed for the efficient characterization of EST-SSRs. This paper describes the discovery and development of EST-SSRs for red drum (Sciaenops ocellatus), an estuarine-dependent sciaenid species of economic importance in the USA and elsewhere, using a phylogenetically informed, comparative genomics approach to primer design. The approach entailed comparing existing genomic resources from species closely allied phylogenetically to red drum, with resources from more distantly related outgroup species. By taking into account the degree to which flanking regions are conserved across taxa, the efficiency of PCR primer design was increased greatly. The amplification success rate for primers designed for red drum was 100 % when using EST libraries from confamilial species and 92 % when using an EST library from a species in the same suborder. The primers developed also amplified EST-SSRs in a wide range of perciform fishes, suggesting potential use in comparative genomics. This study demonstrates that EST-SSRs can be efficiently developed for an organism when limited species-specific data are available by exploiting genomic resources from well-studied species, even those at extended phylogenetic distances.

Fast skeletal muscle transcriptome of the gilthead sea bream (Sparus aurata) determined by next generation sequencing

Background

The gilthead sea bream (Sparus aurata L.) occurs around the Mediterranean and along Eastern Atlantic coasts from Great Britain to Senegal. It is tolerant of a wide range of temperatures and salinities and is often found in brackish coastal lagoons and estuarine areas, particularly early in its life cycle. Gilthead sea bream are extensively cultivated in the Mediterranean with an annual production of 125,000 metric tonnes. Here we present a de novo assembly of the fast skeletal muscle transcriptome of gilthead sea bream using 454 reads and identify gene paralogues, splice variants and microsatellite repeats. An annotated transcriptome of the skeletal muscle will facilitate understanding of the genetic and molecular basis of traits linked to production in this economically important species.

Results

Around 2.7 million reads of mRNA sequence data were generated from the fast myotomal of adult fish (~2 kg) and juvenile fish (~0.09 kg) that had been either fed to satiation, fasted for 3-5d or transferred to low (11°C) or high (33°C) temperatures for 3-5d. Newbler v2.5 assembly resulted in 43,461 isotigs >100 bp. The number of sequences annotated by searching protein and gene ontology databases was 10,465. The average coverage of the annotated isotigs was x40 containing 5655 unique gene IDs and 785 full-length cDNAs coding for proteins containing 58–1536 amino acids. The v2.5 assembly was found to be of good quality based on validation using 200 full-length cDNAs from GenBank. Annotated isotigs from the reference transcriptome were attributable to 344 KEGG pathway maps. We identified 26 gene paralogues (20 of them teleost-specific) and 43 splice variants, of which 12 had functional domains missing that were likely to affect their biological function. Many key transcription factors, signaling molecules and structural proteins necessary for myogenesis and muscle growth have been identified. Physiological status affected the number of reads that mapped to isotigs, reflecting changes in gene expression between treatments.

Conclusions

We have produced a comprehensive fast skeletal muscle transcriptome for the gilthead sea bream, which will provide a resource for SNP discovery in genes with a large effect on production traits of commercial interest and for expression studies of growth and adaptation.

Exploring neutral and adaptive processes in expanding populations of gilthead sea bream, Sparus aurata L., in the North-East Atlantic

Recent studies in empirical population genetics have highlighted the importance of taking into account both neutral and adaptive genetic variation in characterizing microevolutionary dynamics. Here, we explore the genetic population structure and the footprints of selection in four populations of the warm-temperate coastal fish, the gilthead sea bream (Sparus aurata), whose recent northward expansion has been linked to climate change. Samples were collected at four Atlantic locations, including Spain, Portugal, France and the South of Ireland, and genetically assayed using a suite of species-specific markers, including 15 putatively neutral microsatellites and 23 expressed sequence tag-linked markers, as well as a portion of the mitochondrial DNA (mtDNA) control region. Two of the putatively neutral markers, Bld-10 and Ad-10, bore signatures of strong directional selection, particularly in the newly established Irish population, although the potential 'surfing effect' of rare alleles at the edge of the expansion front was also considered. Analyses after the removal of these loci suggest low but significant population structure likely affected by some degree of gene flow counteracting random genetic drift. No signal of historic divergence was detected at mtDNA. BLAST searches conducted with all 38 markers used failed to identify specific genomic regions associated to adaptive functions. However, the availability of genomic resources for this commercially valuable species is rapidly increasing, bringing us closer to the understanding of the interplay between selective and neutral evolutionary forces, shaping population divergence of an expanding species in a heterogeneous milieu.

Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L

Gilthead sea bream (Sparus aurata L.) is an important marine fish in Mediterranean aquaculture. Sex determination by age and/or body weight is a critical life-history trait, the genetic basis for which is largely unknown in this sequential hermaphrodite species. Herein, we performed a partial genome scan to map quantitative trait loci (QTL) affecting body weight and sex using 74 informative microsatellite markers from 10 paternal half-sib families to construct nine linkage groups (LG). In total, four growth-related QTL (two chromosome-wide and two genome-wide) and six QTL related to sex determination (three pairs in three different LGs) were detected (two chromosome-wide and one genome-wide). The proportion of phenotypic variation explained by the body-weight QTL ranged from 9.3% to 17.2%, showing their potential for use in marker-assisted selection. The results obtained offer solid ground to investigate the structure and function of the genomic regions involved in the mechanisms of sex reversal.

In silico mining and characterization of 12 EST-SSRs for the invasive slipper limpet Crepidula fornicata

In silico mining of an expressed sequence tags (ESTs) library was found to be efficient at isolating simple sequence repeats (SSRs) loci in the non-indigenous marine mollusc Crepidula fornicata. Twelve SSR loci were developed for routine genotyping. Cross-species amplification to 8 other Crepidula species showed that the 12 loci are highly specific for C. fornicata. Mendelian inheritance was shown for 11 of them (1 being monomorphic in the analyzed offspring array). The genetic diversity for 88 adults was found to be variable across the 12 loci (2-40 alleles, expected heterozygosity between 0.023 and 0.898) with a high overall exclusion probability of 0.99. The degree of genetic polymorphism found here is similar to that shown for 7 anonymous SSRs previously developed and here used on the same samples. This set of 12 specific loci is relevant to perform reliable population and relatedness analyses in Crepidula fornicata.

Transcriptome analysis in Concholepas concholepas (Gastropoda, Muricidae): mining and characterization of new genomic and molecular markers

The marine gastropod Concholepas concholepas, locally known as the "loco", is the main target species of the benthonic Chilean fisheries. Genetic and genomic tools are necessary to study the genome of this species in order to understand the molecular basis of its development, growth, and other key traits to improve the management strategies and to identify local adaptation to prevent loss of biodiversity. Here, we use pyrosequencing technologies to generate the first transcriptomic database from adult specimens of the loco. After trimming, a total of 140,756 Expressed Sequence Tag sequences were achieved. Clustering and assembly analysis identified 19,219 contigs and 105,435 singleton sequences. BlastN analysis showed a significant identity with Expressed Sequence Tags of different gastropod species available in public databases. Similarly, BlastX results showed that only 895 out of the total 124,654 had significant hits and may represent novel genes for marine gastropods. From this database, simple sequence repeat motifs were also identified and a total of 38 primer pairs were designed and tested to assess their potential as informative markers and to investigate their cross-species amplification in different related gastropod species. This dataset represents the first publicly available 454 data for a marine gastropod endemic to the southeastern Pacific coast, providing a valuable transcriptomic resource for future efforts of gene discovery and development of functional markers in other marine gastropods.