Zhaoxia Cui

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): Unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L. crocea individuals and two L. polyactis individuals. Interestingly, all seven CR sequences from L. crocea were identical while the three sequences from L. polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L. crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L. crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys.

Molecular cloning and characterization of a heat shock protein 70 gene in swimming crab (Portunus trituberculatus).

Hsp70 can stimulate cells of the innate immune system directly by acting as danger-signaling molecules. To understand the immune defense mechanisms of swimming crab Portunus trituberculatus (Decapoda: Brachyura: Portunidae), the cDNA of Hsp70 (designated PtHsp70) was cloned by the combination of homology cloning and rapid amplification of cDNA ends (RACE). The full-length PtHsp70 cDNA was 2195 bp, including an open reading frame (ORF) of 1950 bp encoding a polypeptide of 650 amino acids with estimated molecular mass of 71.1514 kDa and theoretical isoelectric point of 5.38. FASTA and BLAST analysis indicated that PtHsp70 should be an inducible cytosolic member of the Hsp70 family. The coding region of PtHsp70 was uninterrupted and four SNPs with 1133C/T, 1311C/T, 1551C/T and 1809 A/G were detected by direct sequencing of 20 genomic samples. Using fluorescent real-time quantitative PCR, the transcriptional expression of PtHsp70 showed a clear time-dependent response after challenge by Vibrio alginolyticus, the main causative agent of emulsification disease causing large mortality in P. trituberculatus. This is the first report on the expression of Hsp70 induced by pathogen stimulation in Brachyura. Phylogenetic analysis revealed that the inducible Hsp70s were divided into two groups in crab and PtHsp70 was clustered into the Hsp/Hsc group (Clade I) by maximum-likelihood (ML) and Bayesian inference (BI) methods. GAP repeat and GGMP motif of inducible Hsp70 gene in the crab species were only found in Clade I.

Transcriptome Analysis and Discovery of Genes Involved in Immune Pathways from Hepatopancreas of Microbial Challenged Mitten Crab Eriocheir sinensis

Background The Chinese mitten crab Eriocheir sinensis is an important economic crustacean and has been seriously attacked by various diseases, which requires more and more information for immune relevant genes on genome background. Recently, high-throughput RNA sequencing (RNA-seq) technology provides a powerful and efficient method for transcript analysis and immune gene discovery. Methods/Principal Findings A cDNA library from hepatopancreas of E. sinensis challenged by a mixture of three pathogen strains (Gram-positive bacteria Micrococcus luteus, Gram-negative bacteria Vibrio alginolyticus and fungi Pichia pastoris; 108 cfu·mL−1) was constructed and randomly sequenced using Illumina technique. Totally 39.76 million clean reads were assembled to 70,300 unigenes. After ruling out short-length and low-quality sequences, 52,074 non-redundant unigenes were compared to public databases for homology searching and 17,617 of them showed high similarity to sequences in NCBI non-redundant protein (Nr) database. For function classification and pathway assignment, 18,734 (36.00%) unigenes were categorized to three Gene Ontology (GO) categories, 12,243 (23.51%) were classified to 25 Clusters of Orthologous Groups (COG), and 8,983 (17.25%) were assigned to six Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Potentially, 24, 14, 47 and 132 unigenes were characterized to be involved in Toll, IMD, JAK-STAT and MAPK pathways, respectively. Conclusions/Significance This is the first systematical transcriptome analysis of components relating to innate immune pathways in E. sinensis. Functional genes and putative pathways identified here will contribute to better understand immune system and prevent various diseases in crab.

Evolution and phylogeny of the mud shrimps (Crustacea: Decapoda) revealed from complete mitochondrial genomes

BackgroundThe evolutionary history and relationships of the mud shrimps (Crustacea: Decapoda: Gebiidea and Axiidea) are contentious, with previous attempts revealing mixed results. The mud shrimps were once classified in the infraorder Thalassinidea. Recent molecular phylogenetic analyses, however, suggest separation of the group into two individual infraorders, Gebiidea and Axiidea. Mitochondrial (mt) genome sequence and structure can be especially powerful in resolving higher systematic relationships that may offer new insights into the phylogeny of the mud shrimps and the other decapod infraorders, and test the hypothesis of dividing the mud shrimps into two infraorders.ResultsWe present the complete mitochondrial genome sequences of five mud shrimps, Austinogebia edulis, Upogebia major, Thalassina kelanang (Gebiidea), Nihonotrypaea thermophilus and Neaxius glyptocercus (Axiidea). All five genomes encode a standard set of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a putative control region. Except for T. kelanang, mud shrimp mitochondrial genomes exhibited rearrangements and novel patterns compared to the pancrustacean ground pattern. Each of the two Gebiidea species (A. edulis and U. major) and two Axiidea species (N. glyptocercus and N. thermophiles) share unique gene order specific to their infraorders and analyses further suggest these two derived gene orders have evolved independently. Phylogenetic analyses based on the concatenated nucleotide and amino acid sequences of 13 protein-coding genes indicate the possible polyphyly of mud shrimps, supporting the division of the group into two infraorders. However, the infraordinal relationships among the Gebiidea and Axiidea, and other reptants are poorly resolved. The inclusion of mt genome from more taxa, in particular the reptant infraorders Polychelida and Glypheidea is required in further analysis.ConclusionsPhylogenetic analyses on the mt genome sequences and the distinct gene orders provide further evidences for the divergence between the two mud shrimp infraorders, Gebiidea and Axiidea, corroborating previous molecular phylogeny and justifying their infraordinal status. Mitochondrial genome sequences appear to be promising markers for resolving phylogenetic issues concerning decapod crustaceans that warrant further investigations and our present study has also provided further information concerning the mt genome evolution of the Decapoda.

Complete mitochondrial genome of the Asian paddle crab Charybdis japonica (Crustacea: Decapoda: Portunidae): gene rearrangement of the marine brachyurans and phylogenetic considerations of the decapods.

Given the commercial and ecological importance of the Asian paddle crab, Charybdis japonica, there is a clearly need for genetic and molecular research on this species. Here, we present the complete mitochondrial genome sequence of C. japonica, determined by the long-polymerase chain reaction and primer walking sequencing method. The entire genome is 15,738xa0bp in length, encoding a standard set of 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes, plus the putative control region, which is typical for metazoans. The total A+T content of the genome is 69.2%, lower than the other brachyuran crabs except for Callinectes sapidus. The gene order is identical to the published marine brachyurans and differs from the ancestral pancrustacean order by only the position of the tRNAHis gene. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of 13 protein-coding genes strongly support the monophyly of Dendrobranchiata and Pleocyemata, which is consistent with the previous taxonomic classification. However, the systematic status of Charybdis within subfamily Thalamitinae of family Portunidae is not supported. C. japonica, as the first species of Charybdis with complete mitochondrial genome available, will provide important information on both genomics and molecular ecology of the group.

A thioredoxin with antioxidant activity identified from Eriocheir sinensis

Thioredoxin, with a redox-active disulfide/dithiol in the active site, is the major ubiquitous disulfide reductase responsible for maintaining proteins in their reduced state. In the present study, the cDNA encoding thioredoxin-1 (designated EsTrx1) was cloned from Chinese mitten crab Eriocheir sinensis by using rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of EsTrx1 was of 641 bp, containing a 5 untranslated region (UTR) of 17 bp, a 3 UTR of 306 bp with a poly (A) tail, and an open reading frame (ORF) of 318 bp encoding a polypeptide of 105 amino acids. The high similarity of EsTrx1 with Trx1s from other animals indicated that EsTrx1 should be a new member of the Trx1 sub-family. Quantitative real-time PCR analysis revealed the presence of EsTrx1 transcripts in gill, gonad, hepatopancreas, muscle, heart and haemocytes. The expression of EsTrx1 mRNA in haemocytes was up-regulated after Listonella anguillarum challenge, reached the maximum level at 6h post-stimulation, and then dropped back to the original level gradually. In order to elucidate its biological functions, EsTrx1 was recombined and expressed in E. coli BL21 (DE3). The rEsTrx1 was demonstrated to possess the expected redox activity in enzymatic analysis, and to be more potent than GSH in antioxidant capacity. These results together indicated that EsTrx1 could function as an important antioxidant in a physiological context, and perhaps is involved in the responses to bacterial challenge.

High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis

The sex determination system in crabs is believed to be XY-XX from karyotypy, but centromeres could not be identified in some chromosomes and their morphology is not completely clear. Using quantitative trait locus mapping of the gender phenotype, we revealed a ZW-ZZ sex determination system in Eriocheir sinensis and presented a high-density linkage map covering ~98.5% of the genome, with 73 linkage groups corresponding to the haploid chromosome number. All sex-linked markers in the family we used were located on a single linkage group, LG60, and sex linkage was confirmed by genome-wide association studies (GWAS). Forty-six markers detected by GWAS were heterozygous and segregated only in the female parent. The female LG60 was thus the putative W chromosome, with the homologous male LG60 as the Z chromosome. The putative Z and W sex chromosomes were identical in size and carried many homologous loci. Sex ratio (5:1) skewing towards females in induced triploids using unrelated animals also supported a ZW-ZZ system. Transcriptome data were used to search for candidate sex-determining loci, but only one LG60 gene was identified as an ankyrin-2 gene. Double sex- and mab3-related transcription factor 1 (Dmrt1), a Z-linked gene in birds, was located on a putative autosome. With complete genome sequencing and transcriptomic data, more genes on putative sex chromosomes will be characterised, thus leading towards a comprehensive understanding of the sex determination and differentiation mechanisms of E. sinensis, and decapod crustaceans in general.

Three isoforms of anti-lipopolysaccharide factor identified from eyestalk cDNA library of swimming crab Portunus trituberculatus.

Anti-lipopolysaccharide factors (ALFs), as the potent antimicrobial peptides, can bind and neutralize lipopolysaccharide (LPS) and exhibit broad spectrum antimicrobial activities. In this study, three isoforms of the ALF homologues (PtesALF1-3) were identified from eyestalk cDNA library of swimming crab Portunus trituberculatus. The full-length cDNA sequences of PtesALF1, 2 and 3 were 1138, 1052 and 1057 bp encoding 92, 108 and 123 amino acids, respectively. PtesALF1-3 contained two conserved cysteine residues and shared high similarity with other reported ALFs. Predicted tertiary structures of PtesALF2 and 3 containing four β-strands and three α-helix were similar to that described in Limulus polyphemus, while PtesALF1 had only one α-helix in its spatial structure. Sequence analysis revealed PtesALF1-3 were encoded by the same genomic locus and generated by alternative splicing of the pre-mRNA. Totally 89 SNPs including 18 in coding region and 71 in noncoding region were detected by direct sequencing of 30 genomic samples. The mRNA expression of PtesALF1 and PtesALF1-3 transcripts was mainly detected in haemocytes but showed different expression pattern in other tissues including hepatopancreas, gill, eyestalk and muscle. After challenge with Vibrio alginolyticus, the temporal expression level of PtesALF1-3 transcripts in haemocytes showed a clear time-dependent response expression pattern with two peaks within the experimental period of 32 h, while PtesALF1 was up-regulated only once with obvious decrease at 6 h and significant increase at 24 h. These results suggest that the PtesALF isoforms have different tissue specificity and might provide multiple protective functions against invading bacteria in P. trituberculatus.

Multiple isoforms of immune-related genes from hemocytes and eyestalk cDNA libraries of swimming crab Portunus trituberculatus

Expressed sequence tags (ESTs) analysis has been shown to be an efficient approach not only for gene discovery, but also for gene expression profiles performance. Two full-length enriched cDNA libraries were constructed from hemocytes and eyestalk of Portunus trituberculatus, respectively, and randomly sequenced to collect genomic information and identify genes involved in immune defense response. A total of 99 unigenes including 64 unigenes (6.00% of 1066 unigenes) in hemocytes library and 35 unigens (6.86% of 510 unigenes) in eyestalk library are identified to be immune genes. These genes are categorized into six classes, viz. antimicrobial peptides, redox proteins, melanization related proteins, chaperone proteins, clottable proteins and other immune factors. The content and category of immune genes in eyestalk library indicate eyestalk might have unrecognized role in crab immunity. Five immune genes containing multiple protein isoforms are identified and characterized, including anti-lipopolysaccharide factor (PtALF1-7), crustin (PtCrustin1-3), thioredoxin (PtTrx1-2), clip domain serine proteinase (PtcSP1-5) and kazal-type proteinase inhibitor (PtKPI1-4). Sequence alignment and phylogenetic analysis reveal PtALF1-7 contain two conserved cysteine residues and might be encoded by multiple genomic loci. PtCrustin1-3 share the consensus cysteine motif and are considered as Type I crustins. PtTrx1 possesses the critical structural cysteine residue C⁷³ of Trx-1, while PtTrx2 has the N-terminal mitochondrial translocation signal of Trx-2. Sequence analysis shows PtcSP1-5 contain one clip domain and one partial SP catalytic triad domain. PtKPI1-4 present one typical Kazal domain consisting of six conserved cysteine residues. Some protein isoforms are tissue-specific, which might suggest they have different origins and perform diverse functions. Except PtALF1-3 and PtCrustin1, the other isoformes in this study are firstly identified from P. trituberculatus. Especially, PtTrx2 are firstly identified from crustaceans. Our research will provide useful genomic information of P. trituberculatus and be helpful in understanding the molecular mechanisms of crab immunity.

Molecular coordinated regulation of gene expression during ovarian development in the penaeid shrimp

To understand the molecular events of ovarian development in penaeid shrimp, RNA arbitrarily primed polymerase chain reaction (RAP-PCR) was used to identify differentially expressed genes during ovarian maturation in Metapenaeus ensis. From a screening of 700 clones in a cDNA library of the shrimp ovary by the products of RAP-PCR of different maturation stages, 91 fragments with differentially expressed pattern as revealed by dot-blot hybridization were isolated and sequenced. Forty-two of these fragments show significant sequence similarity to known gene products and the differentially expressed pattern of 10 putative genes were further characterized via Northern hybridization. Putative glyceraldehyde–3–phosphate dehydrogenase and arginine kinase are related to provision of energy for active cellular function in oocyte development. Translationally controlled tumor protein, actin, and keratin are related to the organization of cytoskeleton to accomplish growth and development of oocytes. High mobility group protein DSP1, heat shock protein 70, and nucleoside diphosphate kinase may act as repressors before the onset of ovarian maturation. Peptidyl-prolyl cis-trans isomerase and glutathione peroxidase are related to the stabilization of proteins and oocytes. This study provides new insights on the molecular events in the ovarian development in the shrimp.