Genetic diversity and population structure of Nibea albiflora in the China Sea revealed by mitochondrial COI sequences

The Complete Mitochondrial Genome of the Luciocyprinus langsoni (Cypriniformes: Cyprinidae): Characterization, Phylogeny, and Genetic Diversity Analysis

Background: Luciocyprinus langsoni is a species belonging to the Cyprinidae family. The objective of this study is to gain a comprehensive understanding of its evolutionary history and genetic characteristics. Methods: The complete mitochondrial genome of L. langsoni was determined using overlapping PCR. A phylogenetic analysis was conducted based on 13 protein-coding genes from 48 species. A population genetic diversity analysis using the COI gene and a selection analysis of 13 protein-coding genes were also performed. Results: The mitogenome is 16,586 base pairs long and consists of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and two control regions. It has a high adenine-thymine (A + T) content. The phylogenetic analysis confirms the placement of L. langsoni within the subfamily Cyprininae. The population genetic diversity analysis reveals low variability in the Hechi Longjiang population. The selection analysis shows that all 13 protein-coding genes have evolved under purifying selection with Ka/Ks ratios below 1. Conclusions: These results enhance our understanding of L. langsoni's evolutionary history and lay a genetic foundation for future studies in population genetics and phylogenetics.

Identification, functional characterization and immune response profiles of interleukin-10 in Nibea albiflora

Interleukin-10 (IL-10) is an immunosuppressive cytokine, which plays a vital role in regulating inflammation for inhibiting the generation and function of pro-inflammatory cytokines in vivo or in vitro. In the present study, the full length cDNA of IL-10 was characterized from Nibea albiflora (named as NaIL-10) of 1238 base pairs (bp), containing a 5'-UTR (untranslated region) of 350 bp, a 3'-UTR of 333 bp and an open reading frame (ORF) of 555 bp (Fig. 1A) to encode 184 amino acid residues with a signal peptide at the N-terminus. The sequence analysis showed that NaIL-10 possessed the typical IL-10 family symbolic motif and conversed cysteine residues, similar to its teleost orthologues. Real-time PCR indicated that NaIL-10 had wide distribution in different healthy tissues, with a relatively high expression in immune-related tissues (head kidney, spleen, kidney, liver and gill). Significantly, up-regulations of NaIL-10 after infection against Vibrio parahaemolyticus, Vibrio alginolyticus and Poly I:C were also observed. Subcellular localization manifested that NaIL-10 mainly distributed in the cytoplasm unevenly and aggregately, and there was also a small amount on the cell membrane, indicating that NaIL-10 was secreted to the extracellular space as the known IL-10 homologous molecules. It could co-locate with IL-10 Rα on the membrane of HEK293T cells for their potential interaction, and GST pull-down and Co-IP studies certified the specific and direct interaction between NaIL-10 and NaIL-10 Rα, confirming that an IL-10 ligand-receptor system existed in N.albiflora. The expression of pro-inflammatory cytokines, including TNF-α, IL-6, IL-1β, were dramatically inhibited in LPS-stimulated RAW264.7 macrophages pre-incubated with recombinant NaIL-10 protein, demonstrating its anti-inflammatory roles. Taken together, the results demonstrated the existence of IL-10 ligand-receptor system in N.albiflora for the first time, and indicated the suppressive function of NaIL-10 on pro-inflammatory cytokine expression in inflammatory response, which would be conducive to better comprehending the role of IL-10 in the immunomodulatory mechanisms of teleost.

Genetic Diversity and Connectivity of Ocypode ceratophthalmus in the East and South China Seas and Its Implications for Conservation

The East and South China Seas are rich in marine resources, but they are also under great pressure from climate change and human activities. Maintaining diversity and connectivity between communities is thought to be effective in mitigating these pressures. To assess the diversity and connectivity among the populations of Ocypode ceratophthalmus in the East and South China Seas, 15 populations from or near 15 marine protected areas in the two seas were studied using COI and D-Loop as genetic markers. The results showed that O. ceratophthalmus populations had high diversity, and the results of a hierarchical analysis of molecular variance and fixation index found that there were no significant genetic structures among these populations. High historical gene flow and high migration rates were further observed among populations by Migrate-n. Furthermore, the COI sequences further showed the asymmetric migration rate with a higher migration rate from south to north than from north to south. This information could provide recommendations for the management of marine protected areas in the East and South China Seas.

Sequence and functional features of a novel scavenger receptor homolog, SCARA5 from Yellow drum (Nibea albiflora)

As an important member in SR-As, member 5 (SCARA5) can swallow apoptotic cells and foreign bodies, and participate multiple signaling pathways to inhibit tumor occurrence, development growth and metastasis. To explore its immune function, SCARA5 was identified from the yellow drum (Nibea albiflora) according to its transcriptome data, and its full-length cDNA was 6968 bp (named as NaSCARA5, GenBank accession no: MW070211) encoding 497 amino acids with a calculated molecular weight of 55.12 kDa, which had the typical motifs of SR family, such as transmembrane helix region, coil region, Pfam collagens region and SR region. BLASTp and the phylogenetic relationship analysis illustrated that the sequences shared high similarity with known SCARA5 of teleosts. Quantitative real time RT-PCR analysis showed that NaSCARA5 was expressed in intestine, stomach, liver, kidney, gill, heart and spleen, with the highest in the spleen (24.42-fold compared with that in heart). After being infected with Polyinosinic:polycytidylic acid (PolyI:C), Vibrio alginolyticus and Vibrio parahaemolyticus, NaSCARA5 mRNA were up-regulated with time dependent mode in spleen, which suggested that NaSCARA5 might play an important role in the immune process of fish. The extracellular domain of NaSCARA5 was successfully expressed in BL21 (DE3), and yielded the target protein of the expected size with many active sites for their conferring protein-protein interaction functions. After being purified by Ni-NAT Superflow resin and renatured, it was found to bind all the tested bacteria (V.parahaemolyticus,V.alginolyticus and Vibrio harveyi). The eukaryotic expression vector of the NaSCARA5-EGFP fusion protein was constructed and transferred into epithelioma papulosum cyprini (EPC) cells, and it was mainly expressed on the cell membrane indicating that NaSCARA5 was a typical transmembrane protein. The aforementioned results indicated that NaSCARA5 played a significant role in the defense against pathogenic bacteria infection as PRRs, which may provide some further understandings of the regulatory mechanisms in the fish innate immune system for SR family.

Molecular phylogenetic and morphometric analysis of population structure and demography of endangered threadfin fish Eleutheronema from Indo-Pacific waters

The threadfin Eleutheronema are the important fishery resources in Indo-Pacific regions and classified as the endangered species with considerable conservation values. Their genetic diversity and population structure remain essentially unknown but are critical for the proper management and sustainable harvests of such important fisheries. Here, the mitochondrial DNA sequences of CO1 and 16s rRNA were determined from 75 individuals of Eleutheronema tetradactylum and 89 individuals of Eleutheronema rhadinum collected from different locations of South China Sea and Thailand coastal waters. Genetic diversity analysis revealed that both E. tetradactylum (Haplotype diversity, H = 0.105–0.211; Nucleotide diversity, π = 0.00017–0.00043) and E. rhadinum (H = 0.074–0.663, π = 0.00013–0.01088) had low diversity. Population structure analysis demonstrated the shallow genetic differentiation among the South China Sea populations. The limited communication between China and Thailand populations caused the high genetic differentiation in all groups due to the low dispersal ability. Reconstruction of CO1 phylogenetic relationships and demographic studies across Indo-West-Pacific regions provided strong evidence for a shared common origin or ancestor of E. tetradactylum and E. rhadinum. Eleutheronema rhadinum were further subdivided into two distinct genetic lineages, with Clade A dominantly distributing in Thailand and Malaysia and Clade B distributing in China coastal waters. Phenotypic divergence, characterized mainly by the depth of caudal peduncle and length of caudal peduncle, was also observed for all populations, which was possibly associated with specific local adaptations to environmental changes. Our study suggested a strong need for the development of proper fishery management strategies and conservation actions for the imperiled Eleutheronema species.

Chromosome-scale assembly and high-density genetic map of the yellow drum, Nibea albiflora

The yellow drum (Nibea albiflora) is an economically important sciaenid fish in East Asian countries. In this study, we sequenced and assembled a near-complete gynogenetic yellow drum genome. We generated 45.63 Gb of Illumina short-reads and 80.27 Gb of PacBio long-reads and assembled them into a 628.01-Mb genome with a contig N50 of 4.42 Mb. Twenty-four chromosomes with a scaffold N50 of 26.73 Mb were obtained using the Hi-C analysis. We predicted a set of 27,069 protein-coding genes, of which 1,581 and 2,583 were expanded and contracted gene families, respectively. The most expanded genes were categorised into the protein binding, zinc-ion binding and ATP binding functional pathways. We built a high-density genetic linkage map that spanned 4,300.2 cM with 24 linkage groups and a resolution of 0.69 cM. The high-quality reference genome and annotated profiles that we produced will not only increase our understanding of the genetic architecture of economic traits in the yellow drum, but also help us explore the evolution and unique biological characteristics of sciaenid fishes.

Genetic Diversity and Population Structure of Gynaephora qinghaiensis in Yushu Prefecture, Qinghai Province Based on the Mitochondrial COI Gene

Gynaephora qinghaiensis (Lepidoptera: Lymantriidae: Gynaephora), a serious economic pest in alpine meadows, is mainly distributed in Yushu prefecture, Qinghai province, China. In this study, we aimed to investigate the genetic diversity and population structure of G. qinghaiensis through analyzing the sequence of 194 mitochondrial cytochrome oxidase subunit (COI) genes (658 bp in length) identified from 10 geographic populations located in three different countries, including Zhiduo, Zaduo, and Chengduo, of Yushu prefecture. Eleven haplotypes were identified from all populations of G. qinghaiensis with high levels of haplotype diversity (0.78500) and low levels of nucleotide diversity (0.00511). High levels of genetic differentiation and low levels of gene flow were also detected among the populations of G. qinghaiensis. Analysis of molecular variance (AMOVA) showed that 90.13% of the variation was attributed to distribution among groups (Chengduo, Zhiduo, and Zaduo), and 5.22% and 4.65% were, respectively, attributed to distribution among populations, within group, and within populations. The result of mantel test showed a highly significant positive correlation (P < 0.01) between F_ST and geographical distance. A maximum likelihood tree showed that most haplotypes were grouped into three clusters corresponding to the three counties, suggesting a significant phylogeographic structure in the populations of G. qinghaiensis. The haplotype networks revealed that H2 may be the most primitive haplotype and the most adaptable in nature. Populations 7# and 8# had haplotype H2 and higher haplotype diversity; therefore, we speculated that the G. qinghaiensis in both populations were more adaptable to the environment and had greater outbreak potential and, therefore, should be focused on in terms of prevention and control. Our findings provide valuable information for further study of the population structure and phylogeny of G. qinghaiensis and provide a theoretical basis for the control of G. qinghaiensis.

Physiological responses to cold and starvation stresses in the liver of yellow drum (Nibea albiflora) revealed by LC-MS metabolomics

As global climate changes, mass mortality in farmed fish associated with the severely cold weather has aroused growing concerns. Yellow drum (Nibea albiflora) is an important maricultured fish in China, whereby its aquaculture suffered from overwinter mortality associated with cold and cold-induced-fasting stresses. Here, by using LC-MS metabolomics combined with transcriptomics, we investigated the physiological responses of yellow drum liver to cold and starvation stresses. The experiment involved four groups: 16 °C fed group (CG1), 16 °C unfed group (CG2), 8 °C fed group (EG1), and 8 °C unfed group (EG2). Under cold stress, a total of 308 and 257 differential metabolites were identified in EG1 vs. CG1 and EG2 vs. CG2, respectively, showing 5 overlapping significant pathways: glutathione metabolism, biosynthesis of unsaturated fatty acids, galactose metabolism, arginine and proline metabolism, and ABC transporters. Intersection analysis identified that glutamate, oxidized glutathione (GSSG), and eicosenoic acid were the common metabolites induced by cold stress. Under starvation stress, a total of 300 and 215 differential metabolites were identified in CG2 vs. CG1 and EG2 vs. EG1, respectively, showing 2 overlapping significant pathways: glutathione metabolism and galactose metabolism. Intersection analysis revealed that glutamate and GSSG were the common metabolites caused by fasting. Under cold and starvation combined stresses, 286 differential metabolites were identified in EG2 vs. CG1, showing 7 influenced pathways: glycerophospholipid metabolism, biosynthesis of unsaturated fatty acids, glutathione metabolism, sphingolipid metabolism, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, autophagy, and purine metabolism. Interestingly, the glutamate and GSSG were induced by both single and combined stresses of cold and starvation treatments. These findings suggest that glutathione metabolism and its related metabolites (glutamate and GSSG) could be potential biomarkers of cold and starvation stresses in yellow drum. Overall, the results of this study provided insights into the physiological regulation in response to cold and starvation stresses in this fish.

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.

Low mtDNA variation and shallow population structure of the Chinese pomfret Pampus chinensis along the China coast

In the present study, population genetic structure and genetic diversity of the Chinese pomfret Pampus chinensis, along the China coast were investigated and compared with that from Indonesia using mitochondrial DNA cytochrome b gene sequences. A total of 28 variable sites (including 18 singleton sites and 10 parsimony information sites) were observed and 23 haplotypes were defined in 330 individuals from 11 localities. The haplotype diversity (H_D ) of the populations ranged from 0·540 to 0·828, the nucleotide diversity (π) ranged from 0·081 to 0·295%. Pairwise F_ST statistics showed that significant genetic divergence occurred among populations from different geographical regions. The high dispersal capabilities, geographic segregation and ocean currents may be responsible for the present population genetic structure in this species. In addition, a population expansion event during the late Pleistocene period was inferred. The time of population expansion was estimated to occur about 117 000-169 000 years ago.

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.

Comparative cytogenetics in three Sciaenid species (Teleostei, Perciformes): evidence of interspecific chromosomal diversification

Background

Species belonging to the Sciaenidae family present a karyotype composed by 48 acrocentric chromosomes and are thus considered a striking example of chromosomal conservation. In this family, three species are extensively studied including Larimichthys crocea, Larimichthys polyactis and Nibea albiflora due to their importance in fishery and aquaculture in East Asia. Despite abundant data of population genetics available for some of them, cytogenetic information on these species is still scarce and obtained by conventional cytogenetic protocols. Therefore, a more detailed cytogenomic investigation was performed in these species to analyze their karyotype differentiation using conventional staining techniques and fluorescence in situ hybridization to map several repetitive DNAs.

Results

The three species showed a slight karyotype differentiation with 4sm + 2st + 42a in L. polyactis, 20st + 28a in L. crocea and 48a in N. albiflora. Additionally, the mapping of repetitive sequences further revealed a number of interspecific differences among them. Particularly, 18S and 5S rDNA sites showed syntenic arrangements in N. albiflora and non-syntenic arrangements in both Larimichthys species. The microsatellites (CA)₁₅ and (GA)₁₅ showed conspicuous terminal clusters in some chromosomes of all species. On the other hand, (CGG)₁₀ repeats, Rex6 elements and U2 snRNA displayed a scattered distribution on the chromosomes.

Conclusions

Although the three Sciaenid species examined displayed a general pattern of karyotypic conservatism, we explored chromosomal diversification among them. The diversificated karyotypic macrostructure is followed by intergeneric evolutionary diversification of the repetitive sequences. The data indicate some degree of intergeneric evolutionary diversification at chromosomal level, and suggest the evolutionary dynamics among Sciaenid species, higher than previously thought. The present cytogenetic data provide new insight into the chromosomal diversification in Sciaenidae, and contribute to inferring the chromosomal rearrangements and trends of karyotype evolution in this fish group.

Population genetic structure and genetic diversity of Chinese pomfret at the coast of the East China Sea and the South China Sea

The Chinese pomfret Pampus chinensis is one of the most economic and ecological important marine fish species in China. In the present study, the population genetic structure and genetic diversity of P. chinensis were evaluated from a total sample size of 180 individuals representing six populations from the East China Sea and the South China Sea using mitochondrial cytochrome c oxidase subunit I (COI) gene. A total of 24 variable sites (including 3 singleton sites and 21 parsimony information sites) were observed, and 18 haplotypes were defined. The haplotype diversity (Hd) of the populations ranged from 0.559 to 0.775, and the nucleotide diversity (π) ranged from 0.330 to 1.090%. Analysis of molecular variance (AMOVA) reveals that the main variation (66.02%) was among individuals within populations. The average pairwise differences and ϕ_ST values indicated significant genetic differentiation between Dongxing population and the other populations. The results of the present study are helpful for the sustainable management and utilization of this species.

Molecular exploration of hidden diversity in the Indo-West Pacific sciaenid clade

The family Sciaenidae, known as croakers or drums, is one of the largest perciform fish families. A recent multi-gene based study investigating the phylogeny and biogeography of global sciaenids revealed that the origin and early diversification of this family occurred in tropical America during the Late Oligocene-Early Miocene before undergoing range expansions to other seas including the Indo-West Pacific, where high species richness is observed. Despite this clarification of the overall evolutionary history of the family, knowledge of the taxonomy and phylogeny of sciaenid genera endemic to the Indo-West Pacific is still limited due to lack of a thorough survey of all taxa. In this study, we used DNA-based approaches to investigate the evolutionary relationships, to explore the species diversity, and to elucidate the taxonomic status of sciaenid species/genera within the Indo-West Pacific clade. Three datasets were herein built for the above objectives: the combined dataset (248 samples from 45 currently recognized species) from one nuclear gene (RAG1) and one mitochondrial gene (COI); the dataset with only RAG1 gene sequences (245 samples from 44 currently recognized species); and the dataset with only COI gene sequences (308 samples from 51 currently recognized species). The latter was primarily used for our biodiversity exploration with two different species delimitation methods (Automatic Barcode Gap Discovery, ABGD and Generalized Mixed Yule Coalescent, GMYC). The results were further evaluated with help of four supplementary criteria for species delimitation (genetic similarity, monophyly inferred from individual gene and combined data trees, geographic distribution, and morphology). Our final results confirmed the validity of 32 currently recognized species and identified several potential new species waiting for formal descriptions. We also reexamined the taxonomic status of the genera, Larimichthys, Nibea, Protonibea and Megalonibea, and suggested a revision of Nibea and proposed a new genus Pseudolarimichthys.