Population genetic diversity of mud crab (Scylla paramamosain) in Hainan Island of China based on mitochondrial DNA

Genetic Variations and Expansion of the Blue Swimmer Crab (Portunus pelagicus) in Southeast Asia

Blue swimmer crabs (Portunus pelagicus) primarily inhabit the coastal regions of Southeast Asia. This study aimed to examine their population genetic structure, phylogeography, and historical demography by applying partial sequences to the control region in the mitochondrial DNA. The analysis focused on 6 populations comprising 166 individual samples collected from Taiwan (Kezailiao Kaohsiung (KK) and Penghu (PH)), mainland China (Xiamen (XM) and Hong Kong (HK)), Vietnam (Hanoi (VN)), and Singapore (SGP). The estimated nucleotide diversity (π) for all of the samples was 0.062, with values ranging between 0.018 (PH) and 0.045 (HK). Our study collected 116 haplotypes and classified them into lineages A and B. The XM, HK, VN, and SGP populations made up lineage A, and the PH and KK populations comprised lineage B. For each lineage, the results highlighted indications of demographic expansion. Lineage A can be traced back to the second-to-last interglacial period, while lineage B possibly originated in the last glacial period. There were significant differences in the FST values among the six populations, except for the pairs HK–XM and PH–KK. In conclusion, the phylogeography and genetic structure of P. pelagicus in Southeast Asia were potentially affected by the Pleistocene glacial cycles and human introduction. These findings further highlight the possible dispersal routes of P. pelagicus throughout Southeast Asia.

Transcriptome and metabolome integration analysis of mud crab Scylla paramamosain challenged to Vibrio parahaemolyticus infection

Vibrio parahaemolyticus (V. parahaemolyticus) is a common pathogen for marine crustacean, which causes severe illnesses in aquatic animals. Therefore, it is meaningful to explore the mechanism during V. parahaemolyticus infection. In this study, to investigate the immune responses of mud crab Scylla paramamosain (S. paramamosain) to V. parahaemolyticus, we established the metabolic and transcriptional profiles of mud crab hemocytes challenged with V. parahaemolyticus. The results indicated that V. parahaemolyticus infection could induce a series of metabolism alterations at both metabolome and transcriptome levels, including biosynthesis of amino acids and Aminoacyl-tRNA, Purine and pyrimidine metabolism, TCA cycle and glutamine metabolism. In this context, through the integration of metabolomics and transcriptomics, our study provided a more comprehensive understanding of the biological process in mud crab against pathogen infection.

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.

Genetic diversity and population structure of native mitten crab (Eriocheir sensu stricto) by microsatellite markers and mitochondrial COI gene sequence

Mitten crab (Eriocheir sensu stricto) is an indigenous and economically important species in East Asia. The genetic diversity and population structure of four wild native Eriocheir s. s. populations, from Yangtze (YZ), Oujiang (OJ), Minjiang (MJ), and Nanliujiang (NLJ) River basin, were studied. Genetic diversity of the four populations was evaluated by eight microsatellite and mitochondrial cytochrome oxidase subunits I (COI) markers. The polymorphism information content (PIC) values ranged from 0.7827 to 0.8580 and Shannon Wiener index (I) values varied from 2.0722 to 2.4088 respectively, according to microsatellite markers analysis. The haplotype diversity (h) values ranged from 0.52101 to 0.87097, while the values of nucleotide diversity (π) varied from 0.00139 to 0.02796 based on mitochondrial COI gene sequence. Four wild Eriocheir s. s. populations all showed high genetic diversity. As microsatellite analysis results, pair-wise F_ST values of YZ, OJ, MJ and NLJ populations were 0.1347, 0.1103 and 0.0943 respectively, which indicated that they were in medium genetic differentiation (0.05 < F_ST < 0.15). This was consistent with mitochondrial COI gene's conclusion. The bottleneck analysis was a hint that Eriocheir s. s. experienced bottleneck effect under SMM models of Sign Test and Wilcoxon Sign-rank Test. Based on genetic structure, YZ, OJ and MJ populations were in the same cluster separated from NLJ population. In summary, by microsatellite and mitochondrial COI markers analysis, all four Eriocheir s. s. populations had high genetic diversity, and genetic differentiation was consistent with geographical distance. Populations OJ and MJ were closer to population YZ genetically, and they had large genetic distance with NLJ population. This study provides a particularly important basis for mitten crab resources conservation and utilization in south China.

Cloning, tissue distribution and nutritional regulation of a fatty acyl Elovl4-like elongase in mud crab, Scylla paramamosain (Estampador, 1949)

In this report, the full-length cDNA of fatty acyl Elovl4-like elongase was cloned from the hepatopancreas of Scylla paramamosain by rapid-amplification of cDNA ends (RACE). To the best of our knowledge, this is the first report of Elovl4-like elongase in crustaceans. The full-length cDNA of Elovl4-like was 1119bp, which included a 5'-terminal untranslated region (UTR) of 58bp, a 3'-terminal UTR of 44bp and an open reading frame (ORF) of 1017bp encoding a polypeptide of 338 amino acids. Tissue distribution analysis revealed that Elovl4-like transcripts are widely distributed in various organs, with high mRNA levels in the hepatopancreas and cranial ganglia. Further, Elovl4-like transcriptional levels in hepatopancreas were up-regulated in proportion to the replacement of dietary fish oil (FO) with soybean oil (SO). The result showed that Elovl4-like transcripts increased about 0.83 and 1.12-fold respectively when SO constituted 80% and 100% of total oil (P<0.05). These results may contribute to better understanding of the long-chain polyunsaturated fatty acids (LC-PUFA) biosynthetic pathway and regulation mechanism in this species.

Genetic diversity and genetic structure of farmed and wild Chinese mitten crab (Eriocheir sinensis) populations from three major basins by mitochondrial DNA COI and Cyt b gene sequences

The Chinese mitten crab, Eriocheir sinensis, is one of the important native crab species in East Asian region, which has been widely cultured throughout China, particularly in river basins of Yangtze, Huanghe and Liaohe. This study was designed to evaluate the genetic diversity and genetic structure of cultured and wild E. sinensis populations from the three river basins based on mitochondrial DNA (mtDNA) cytochrome oxidase subunit I (COI) and cytochrome b (Cyt b). The results showed that there were 62 variable sites and 30 parsimony informative sites in the 647 bp of sequenced mtDNA COI from 335 samples. Similarly, a 637 bp segment of Cyt b provided 59 variable sites and 26 parsimony informative sites. AMOVA showed that the levels of genetic differentiation were low among six populations. Although the haplotype diversity and nucleotide diversity of Huanghe wild population had slightly higher than the other populations, there were no significant differences. There was no significant differentiation between the genetic and geographic distance of the six populations, and haplotype network diagram indicated that there may exist genetic hybrids of E. sinensis from different river basins. The results of clustering and neutrality tests revealed that the distance of geographical locations were not completely related to their genetic distance values for the six populations. In conclusion, these results have great significance for the evaluation and exploitation of germplasm resources of E. sinensis.

Complete sequence and characterization of the Ectropis oblique mitochondrial genome and its phylogenetic implications

In the present study, we sequenced the entire mitochondrial genome (mitogenome) of Ectropis oblique using PCR amplification and sequencing methods The entire mitogenome is 15,356bp long, including 13 protein coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and A+T rich element. The base composition and gene arrangement are identical to those of other lepidopterans. Of these 37 genes, twenty three are resided on heavy strand, while fourteen located on light strand. The newly sequenced mitogenome displayed a biased Adenine/thymine 84.94% usage versus guanine/cytosine. The AT skewness is 0.030 and the GC skewness is -0.180. Twelve out of 13 PCGs initiated with canonical start codon (ATN), while cox1 started with CGA. The A+T rich element of E. oblique is 363bp long and contains of many features common to lepidopteran insects, including the 'ATAGA' motif, a 22bp poly (T) stretch and a microsatellite-like (AT)₈ element upstream of trnM. Further, this mitogenome comprises 15 intergenic spacer and overlapping regions. Phylogenetic analyses showed that E. oblique belongs to Geometridae, and that the monophyly of Lepidoptera superfamilies is well supported.

Genetic structure and historical demography of the blue swimming crab (Portunus pelagicus) from southeastern sea of China based on mitochondrial COI gene

In this study, the population genetic structure and historical demography of the blue swimming crab, Portunus pelagicus, from southeastern sea of China were investigated using cytochrome c oxidase subunit I (COI) gene of mitochondrion. A total of 889 bp segment of COI gene was sequenced, which showed a high haplotype diversity (0.6833-0.8142) and low nucleotide diversity (0.0021-0.0034). Among 30 haplotypes defined in this study, one (H1) was the most dominant (47.7%) and shared by each locality, while the majority (23) were rare and only existed in one individual. The AMOVA analysis revealed a limited population genetic structure, which suggested a high level of gene flow along the distribution areas of China. This conclusion was supported by the pairwise F_ST comparison values. The topology of the neighbour-joining tree constructed using 30 haplotypes from four localities presented two distinct clades (clade A and clade B). Meanwhile, three sequences of P. pelagicus downloaded from NCBI database showed a high-level divergence with the individuals collected in our study, which might form a new cryptical species. The individuals of clade B were cryptically embedded in the whole population, with a low frequency (7.7-24.2%), while clade A accounted for 75.8-92.3%. Neutrality tests and mismatch analyses suggested a late Pleistocene population expansion for both clade A (47,000-66,000 years ago) and clade B (74,000-100,000 years ago). This study should provide insight into phylogeny, population genetic structure, conservation genetics, and sustainable management of P. pelagicus.

Genetic diversity and population structure of Portunus sanguinolentus (Herbst, 1783) revealed by mtDNA COI sequences

In this paper, population genetic diversity and genetic structure of Portunus sanguinolentus distributed along China coast was investigated based on mitochondrial DNA cytochrome c oxidase subunit I (COI) gene. Sixty-eight haplotypes and 71 polymorphic sites were detected in 110 individuals sampled from four localities. The average haplotype diversity (h) and nucleotide diversity (π) were 0.9576 and 0.0051, respectively. Pairwise F_ST statistics among different localities were negative or low and insignificant, suggesting a widespread gene flow. The AMOVA analysis also indicated a broad genetic exchange among four localities (F_ST= -0.0039, p = 0.559). The high dispersal capabilities of larvae of P. sanguinolentus coupled with the limited physical barriers in marine environment may be responsible for the high gene flow among P. sanguinolentus populations along southeastern coast of China. Moreover, neutrality test and mismatch analysis implied that P. sanguinolentus population might have undergone a rapid population expansion event. The time of population expansion was estimated to occur about 229,000-321,000 years ago.

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.

Mitochondrial DNA diversity of mud crab Scylla olivacea (Portunidae) in Peninsular Malaysia: a preliminary assessment

A primary factor in population management and wildlife conservation is the delineation of population units derived from descriptions of population genetic structure. Yet, predicting factors that influence the patterns of gene flow in a population particularly at landscape scales remains a major challenge in evolutionary biology. Here we report a population genetic study of the mud crab Scylla olivacea examined based on a 542 bp segment of the mitochondrial DNA cytochrome c oxidase I gene among 91 individuals from six localities in the west and east coast of Peninsular Malaysia. In total 55 unique haplotypes were distinguished with 45 private haplotypes and a single common haplotype shared among all populations studied. The other ten haplotypes were shared among various populations. The sharing of this haplotype reflects the connection of the mangrove areas between east and west coast of Peninsular Malaysia. High haplotype diversity (h = 0.968 ± 0.021; mean ± SD) and low nucleotide diversity (π = 0.120 ± 0.015; mean ± SD) were displayed, which may be indicative of genetic bottleneck events. No significant phylogenetic lineages were recognized using neighbour-joining and maximum parsimony methods. Hierarchical AMOVA analysis indicated that 99.33 % of the genetic variation was contained within populations and 0.67 % occurred among populations, suggesting no geographical patterning among populations studied, supported by F st test. Mismatch distribution analysis showed that the observed distribution of the pairwise mutation differences among haplotypes was multimodal, which is not concordant with a sudden range expansion scenario. However, neutrality tests showed non-significant negative values suggesting that the populations studied may have experienced past population growth, but the expansion may have been restricted to separate local areas that resulted in the non-significant negative Fu's Fs and Tajima's D value. Overall, this present preliminary study was able to be a reference on the phylogenetic relationships and assessment of genetic structure of Scylla sp. in Malaysia.

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.

Coupling genetic and species distribution models to examine the response of the Hainan Partridge (Arborophila ardens) to late quaternary climate

Understanding the historical dynamics of animal species is critical for accurate prediction of their response to climate changes. During the late Quaternary period, Southeast Asia had a larger land area than today due to lower sea levels, and its terrestrial landscape was covered by extensive forests and savanna. To date, however, the distribution fluctuation of vegetation and its impacts on genetic structure and demographic history of local animals during the Last Glacial Maximum (LGM) are still disputed. In addition, the responses of animal species on Hainan Island, located in northern Southeast Asia, to climate changes during the LGM are poorly understood. Here, we combined phylogeographic analysis, paleoclimatic evidence, and species distribution models to examine the response of the flightless Hainan Partridge (Arborophila ardens) to climate change. We concluded that A. ardens survived through LGM climate changes, and its current distribution on Hainan Island was its in situ refuge. Range model results indicated that A. ardens once covered a much larger area than its current distribution. Demographic history described a relatively stable pattern during and following the LGM. In addition, weak population genetic structure suggests a role in promoting gene flow between populations with climate-induced elevation shifts. Human activities must be considered in conservation planning due to their impact on fragmented habitats. These first combined data for Hainan Partridge demonstrate the value of paired genetic and SDMs study. More related works that might deepen our understanding of the responses of the species in Southeast Asia to late Quaternary Climate are needed.

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

The mud crab (Scylla paramamosain) is a carnivorous portunid crab, mainly distributed along the southeastern coast of China. Mitochondrial DNA analysis in a previous study indicated a high level of genetic diversity and a low level of genetic differentiation. In this study, population genetic diversity and differentiation of S. paramamosain were investigated using nine microsatellite markers. In total, 397 wild specimens from 11 locations on the southeastern coast of China were sampled and genotyped. A high level of genetic diversity was observed, with the number of alleles, and the observed and expected heterozygosity per location in the range 7.8-9.6, 0.62-0.77 and 0.66-0.76, respectively. AMOVA analysis indicated a low level of genetic differentiation among the 11 locations, despite the fact that a statistically significant fixation index (F(ST)) value was found (F(ST)=0.0183, P<0.05). Out of 55 pairwise location comparisons, 39 showed significant F(ST) values (P<0.05), but all of them were lower than 0.05, except for one between Sanmen and Shantou locations. No significant deficiency of heterozygotes (inbreeding coefficient F(IS)=0.0007, P>0.05) was detected for all locations except Sanmen and Zhanjiang. Cluster analysis using UPGMA showed that all locations fell into one group except Sanmen. Significant association was found between genetic differentiation in terms of F(ST)/(1-F(ST)) and the natural logarithm of geographical distance (r(2)=0.1139, P=0.02), indicating that the genetic variation pattern closely resembled an isolation by distance model. This study supports the proposal of high genetic diversity and low genetic differentiation in S. paramamosain along the southeastern coast of China.