High genetic diversity and low differentiation in mud crab (Scylla paramamosain) along the southeastern coast of China revealed by microsatellite markers

Population genetic diversity of mud crab (Scylla paramamosain) from southeast coastal regions of China based on mitochondrial COI gene sequence

Mud crab (Scylla paramamosain) is a native and economically important species in East Asia. In order to provide a comprehensive genetic background for the resource protection and management in this species, the genetic diversity and population structure of S. paramamosain were analyzed based on mitochondrial cytochrome oxidase subunit I (COI) gene sequence. Totally, 599 individuals were sampled from 20 populations, including 18 wild and two cultured populations from five provinces along southeastern coast of China. After the sequencing of a 638 bp fragment of COI gene, 84 variable sites were found and no insertion or deletion was detected. The H2 (haplotype 2) was the dominant hapolotye shared by 301 individuals (50.3% of all individuals) and existed in all localities. In addition, a high percent of unique haplotypes (53 of 93 haplotyoes) was found. The average nucleotide diversity (π) of all populations was 0.00194, ranging from 0.00010 (NHF1) to 0.00305 (SHCM). The haplotype diversity (h) ranged from 0.067 (NHF1) to 0.876 (CM) with an average 0.738. All of the populations showed high h (>0.5) except NHF1. Genetic distance ranged from 0.00063 to 0.00337 between populations and from 0.00010 to 0.00374 within populations. The molecular variance analysis (AMOVA) showed that the total genetic variation mainly occurred within populations (99.68%) and only 0.32% was contributed by among populations variation. No significant genetic differentiation was observed among all wild populations except that between BA and SHCM (Fst = 0.0.05707, P < 0.05), indicating a low level of genetic differentiation among localities. It is worth to note that two progeny population (HNF1 and NHF1) showed significantly different genetic variation, which suggested that a large quantity of parents could help to improve the genetic diversity of progeny populations.

Genome survey and development of polymorphic microsatellite loci for Sillago sihama based on Illumina sequencing technology

In this study, we first conducted a genome survey assay for Sillago sihama by Illumina sequencing platform, and then developed 15 polymorphic microsatellite loci in a wild population. A total of 129.46 Gb raw data were obtained, of which 115.07 Gb were clean data, with a sequencing depth of 179.3-folds. This genome was estimated to be 522.6 Mb in size, with the heterozygosity, repeat content and GC content being 0.63%, 21% and 44%. A total of 630,028 microsatellites were identified from the genome, of which, dinucleotide repeat was the most abundant (56.80%), followed by mononucleotide repeat (30.23%). Furthermore, 60 pairs of primers were designed and synthesized based on microsatellite sequences, of which 15 were polymorphic in a wild population. A total of 91 alleles were found, with an average of 6.07 per locus. Number of alleles, observed and expected heterozygosity per locus ranged from two to 13, from 0.250 to 0.862, and from 0.396 to 0.901, respectively. Twelve loci were highly informative (PIC > 0.5), and the others were medium informative (0.25 < PIC < 0.5). Seven loci deviated from Hardy-Weinberg equilibrium after Bonferroni correction (P < 0.0033). No significant linkage disequilibrium was detected between loci pairs. This study provided a large number of genomic resources and 15 polymorphic microsatellite loci that should be helpful for the further genetic studies in S. sihama.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

High Connectivity among Blue Crab (Callinectes sapidus) Populations in the Western South Atlantic

Population connectivity in the blue crab Callinectes sapidus was evaluated along 740 km of the Western South Atlantic coast. Blue crabs are the most exploited portunid in Brazil. Despite their economic importance, few studies report their ecology or population structure. Here we sampled four estuarine areas in southern Brazil during winter 2013 and summer 2014 in order to evaluate diversity, gene flow and structure of these populations. Nine microsatellite markers were evaluated for 213 adult crabs, with identification of seven polymorphic loci and 183 alleles. Pairwise FST values indicated low population structure ranging from -0.00023 to 0.01755. A Mantel test revealed that the geographic distance does not influence genetic (r = -0.48), and structure/migration rates confirmed this, showing that even the populations located at the opposite extremities of our covered region presented low FST and exchanged migrants. These findings show that there is a significant amount of gene flow between blue crab populations in South Brazil, likely influenced by local current dynamics that allow the transport of a high number of larvae between estuaries. Considering the elevated gene flow, the populations can be considered a single genetic stock. However, further information on population size and dynamics, as well as fishery demands and impacts at different regions, are necessary for harvest management purposes.

Transcriptome analysis of the mud crab (Scylla paramamosain) by 454 deep sequencing: assembly, annotation, and marker discovery

In this study, we reported the characterization of the first transcriptome of the mud crab (Scylla paramamosain). Pooled cDNAs of four tissue types from twelve wild individuals were sequenced using the Roche 454 FLX platform. Analysis performed included de novo assembly of transcriptome sequences, functional annotation, and molecular marker discovery. A total of 1,314,101 high quality reads with an average length of 411 bp were generated by 454 sequencing on a mixed cDNA library. De novo assembly of these 1,314,101 reads produced 76,778 contigs (consisting of 818,154 reads) with 5.4-fold average sequencing coverage. The remaining 495,947 reads were singletons. A total of 78,268 unigenes were identified based on sequence similarity with known proteins (E≤0.00001) in UniProt and non-redundant protein databases. Meanwhile, 44,433 sequences were identified (E≤0.00001) using a BLASTN search against the NCBI nucleotide database. Gene Ontology (GO) analysis indicated that biosynthetic process, cell part, and ion binding were the most abundant terms in biological process, cellular component, and molecular function categories, respectively. Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analysis revealed that 4,878 unigenes distributed in 281 different pathways. In addition, 19,011 microsatellites and 37,063 potential single nucleotide polymorphisms were detected from the transcriptome of S. paramamosain. Finally, thirty polymorphic microsatellite markers were developed and used to assess genetic diversity of a wild population of S. paramamosain. So far, existing sequence resources for S. paramamosain are extremely limited. The present study provides a characterization of transcriptome from multiple tissues and individuals, as well as an assessment of genetic diversity of a wild population. These sequence resources will facilitate the investigation of population genetic diversity, the development of genetic maps, and the conduct of molecular marker-assisted breeding in S. paramamosain and related crab species.

Characterization of 40 single nucleotide polymorphism (SNP) via Tm-shift assay in the mud crab (Scylla paramamosain)

In this study, single nucleotide polymorphism (SNP) were identified, confirmed and genotyped in the mud crab (Scylla paramamosain) using Tm-shift assay. High quality sequences (13, 311 bp long) were obtained by re-sequencing that contained 91 SNPs, with a density of one SNP every 146 bp. Of all 91 SNPs, 40 were successfully genotyped and characterized using 30 wild specimens by Tm-shift assay. The minor allele frequency per locus ranged from 0.017 to 0.500. The observed and expected heterozygosity, and polymorphism information content (PIC) ranged from 0.000 to 0.600, from 0.033 to 0.509, and from 0.033 to 0.375, respectively, with an average of 0.142, 0.239 and 0.198 per locus. Seventeen SNPs were significantly deviated from Hardy-Weinberg equilibrium. No significant linkage disequilibrium between pairs of loci was detected after sequential Bonferroni correction (P > 0.00125). Seventeen SNPs were related with known function genes. This study provided new molecular markers for investigation of population genetic diversity, construction of genetic linkage maps and molecular marker-assisted selection in this important crustacean species.

Identification of transcriptome-derived microsatellite markers and their association with the growth performance of the mud crab (Scylla paramamosain)

Microsatellite markers from a transcriptome sequence library were initially isolated, and their genetic variation was characterized in a wild population of the mud crab (Scylla paramamosain). We then tested the association between these microsatellite markers and the growth performance of S. paramamosain. A total of 129 polymorphic microsatellite markers were identified, with an observed heterozygosity ranging from 0.19 to 1.00 per locus, an expected heterozygosity ranging from 0.23 to 0.96 per locus, and a polymorphism information content (PIC) ranging from 0.21 to 0.95 per locus. Of these microsatellite markers, 30 showed polymorphism in 96 full-sib individuals of a first generation family. Statistical analysis indicated that three microsatellite markers were significantly associated with 12 growth traits of S. paramamosain. Of these three markers, locus Scpa36 was significantly associated with eight growth traits, namely, carapace length, abdomen width (AW), body height (BH), fixed finger length of the claw, fixed finger width of the claw, fixed finger height of the claw, meropodite length of pereiopod 2, and meropodite length of pereiopod 3 (MLP3) (P<0.05). Locus Scpa75 was significantly associated with five growth traits, namely, internal carapace width, AW, carapace width at spine 8, distance between lateral spine 2 (DLS2), and MLP3 (P<0.05). Locus Spm30 was significantly associated with BH, DLS2, and body weight (P<0.05). Further analysis suggested a set of genotypes (BC at Scpa36, BC and BD at Scpa75, and AC at Spm30) that have great potential in the selection of S. paramamosain for growth traits. These findings will facilitate the development of population conservation genetics and molecular marker-assisted selective breeding of S. paramamosain and other closely related species.

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA(His) differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.