The analyses of the complete mitochondrial genomes of three crabs revealed novel gene rearrangements and phylogenetic relationships of Brachyura.

BACKGROUND: Brachyura crab is the largest branch of Decapoda crustacean. Phylogenetic relationships within Brachyura remain controversial to be investigated. The mitochondrial genome (mitogenome) is an important molecular marker for studying the phylogenetic relationships of Brachyura. METHODS AND RESULTS: To understand the phylogeny of Brachyura, the three complete mitogenomes from Charybdis annulata, Leptodius exaratus, and Spider crab were sequenced and annotated. Their full length was 15,747, 15,716, and 16,608 bp long, respectively. The first two crabs both contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. However, Spider crab contained 13 PCGs, two rRNA genes, 25 tRNA genes and a control region. The mitogenomes of each of the three crabs exhibited high AT content (67.8%, 69.1%, and 70.8%), with negative AT skews (-0.014, - 0.028, and - 0.017) and GC skews (-0.269, - 0.286, and - 0.341). The gene order of C. annulata was identical to the ancestor of Brachyura. Compared with the ancestor of Brachyura, L. exaratus exhibited the gene rearrangements of Val (V)-rrnS-control region, and Spider crab had the four copies of Lys (K). Phylogenetic analyses indicated that C. annulata belonged to Portunidae family, Portunoidea superfamilies, L. exaratus belonged to Xanthidae family, Xanthoidea superfamilies, and Spider crab belonged to Mithracidae family, Majoidea superfamilies. Phylogenetic analyses showed that the two species (Somanniathelphusa boyangensis and Huananpotamon lichuanense) belonging to the Potamoidea were sister groups to the Thoracotremata, thus supporting the conclusion that Heterotremata is polyphyletic. CONCLUSION: The results of this study enriched the crab mitogenome database and enabled us to better understand the phylogenetic relationships of Brachyura.

Phylogenetic relationships of Grapsoidea and insights into the higher phylogeny of Brachyuran.

Decapoda is one of the most diverse crustacean orders, and has become an important research subject. However, the phylogenetic relationships among the main lineages of Decapoda remain uncertain, especially in the order Brachyura. Herein, we sequenced the whole mitochondrial genome of V. litterata and constructed a phylogenetic tree to understand its phylogenetic relationships with other species. The results showed that the mitochondrial genome of V. litterata was generally similar to mitogenomes of Metazoa reported in the literature, with a size of 16,247 bp, 37 genes, and a control region. Both AT-skew and GC-skew were negative, indicating more abundant Cs and Ts than Gs and As. The gene arrangement of V. litterata is identical to those of Eriocheir hepuensis, Cyclograpsus granulosus, Hemigrapsus sanguineus, Helicana wuana, and Helice tientsinensis but differs from the pancrustacean ground pattern and typical arrangement of Brachyuran crabs. Phylogenetic reconstruction showed that V. litterata belongs to the Varunidae.

Transcriptome analysis of immune-related genes in Sesarmops sinensis hepatopancreas in reaction to peptidoglycan challenge.

Sesarmops sinensis is a dominant omnivorous crab species, which plays an important ecological function in salt marsh ecosystems. To better understand its immune system and immune related genes under pathogen infection, the transcriptome was analyzed by comparing the data of S. sinensis hepatopancreas stimulated by PBS and PGN. A set of assembly and annotation identified 39,039 unigenes with an average length of 1105 bp, obtaining 1300 differentially expressed genes (DEGs) in all, which included 466 remarkably up-regulated unigenes and 834 remarkably down-regulated unigenes. In addition, based on mensurable real time-polymerase chain reaction and high-throughput sequencing, several immune responsive genes were found to be markedly up-regulated under PGN stimulation. In conclusion, in addition to enriching the existing transcriptome data of S. sinensis, this study also clarified the immune response of S. sinensis to PGN stimulation, which will help us to further understand the crustacean's immune system.

A complement factor I (CFI) gene mediates innate immune responses in yellow catfish Pelteobagrus fulvidraco.

This study isolated CFI gene from Pelteobagrus fulvidraco and named it PfCFI. The cDNA of PfCFI is 2374 bp long, including a 52 bp 5' untranslated sequence, a 222 bp 3' untranslated sequence, and an open reading frame (ORF) of 2100 bp encoding polypeptide consisting of 699 amino acids. Phylogenetic analysis revealed that the PfCFI was closely related to CFI of Ictalurus punctatus. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis indicate that there is the PfCFI gene which expressed in all the rest of tested tissues in varied levels, and mainly distributed in liver and least in heart. The reseachers induce the expressions level of PfCFI gene in liver, spleen, head kidney and blood at different points in time after challenged with lipopolysaccharide (LPS), and polyriboinosinic polyribocytidylic acid (poly I:C), respectively. Together these results suggested that CFI gene plays an important role in resistance to pathogens in yellow catfish immunity.

Chitinase involved in immune regulation by mediated the toll pathway of crustacea Procambarus clarkii.

Chitinase can degrade chitin and play an essential role in animal immunity and plant defense. The immune functions of Chitinase in Procambarus clarkii (P. clarkii) remain to elucidate. Here, we identified PcChitinase 2 gene sequence from P. clarkii and studied its spatial and temporal expression profiles. The PcChitinase 2 transcribed unequally in different tissues; however, its expression was highest in those of stomach, gut, and hepatopancreas. The challenge with lipolysaccharide or peptidoglycan significantly up-regulated the expression of PcChitinase 2 in hepatopancreas. The knockdown of the PcChitinase 2 gene by double-stranded RNA suppressed most of the Toll-pathway-related immune genes (phospholipase, lectin, sptazle Cactus, serine proteikinase, anti-lipopolysaccharide factor, and Toll) production were significantly increased. Our results suggest PcChitinase 2 may be involved in the innate immune responses of P. clarkii by modulating the toll pathway.

Transcriptome analysis of differentially expressed genes in the red swamp crayfish Procambarus clarkii challenged with Aeromonas hydrophila.

As an important economic species in China, aquaculture of the crayfish Procambarus clarkii has suffered huge losses due to infection by pathogenic bacteria, mainly by Aeromonas hydrophila, which leads to high mortality and huge economic loss. To better understand the immune response of crayfish against bacterial infection, we compared and analyzed transcriptome data of hepatopancreatic tissue from P. clarkii that were either challenged with A. hydrophila or treated with PBS. After assembly and annotation of the data, 32,041 unigenes with an average length of 1512 base pairs were identified. Compared to control group, Differential gene expression (DEG) analysis revealed 608 DEGs were obtained, of which 274 unigenes were upregulated and 334 were downregulated in the A. hydrophila group. Furthermore, the expression levels of eight selected immune-related DEGs were validated by qRT-PCR, substantiating the reliability of RNA-seq results. This study not only provides effective data support for immune defense strategies of P. clarkii in response to bacterial infections, but also provides new information about the P. clarkii immune system and defense mechanisms, and a valuable basis for further studies to elucidate the molecular immune mechanisms of this species.

Bmcas-1 plays an important role in response against BmNPV infection in vitro.

Apoptosis, as one kind of innate immune system, is involved in host response against pathogens innovation. Caspases play a vital role in the execution stage of host cell apoptosis. It has been reported that Bmcaspase-1 (Bmcas-1) has a close relationship with Bombyx mori nucleopolyhedrovirus (BmNPV) infection for its differentially expressed patterns after viral infection. However, its underlying response mechanism is still unclear. The significant differential expression of Bmcas-1 in different tissues of differentially resistant strains revealed its vital role in BmNPV infection. To further validate its role in BmNPV infection, budded virus (BV)-eGFP was analyzed after knockdown and overexpression of Bmcas-1 by small interfering RNA and the pIZT-mCherry vector, respectively. The reproduction of BV-eGFP obviously increased at 72 h after knockdown of Bmcas-1, and decreased after overexpression in BmN cells. Moreover, the conserved functional domain of Cas-1 among different species and the closed evolutionary relationship of Cas-1 in Lepidoptera hinted that Bmcas-1 might be associated with apoptosis, and this was also validated by the apoptosis inducer, Silvestrol, and the inhibitor, Z-DEVD-FMK. Therefore, Bmcas-1 plays an essential antiviral role by activating apoptosis, and this result lays a fundament for clarifying the molecular mechanism of silkworm in response against BmNPV infection and breeding of resistant strains.

Transcriptome analysis reveals antioxidant defense mechanisms in the red swamp crayfish Procambarus clarkia after exposure to chromium.

Chromium (Cr) as a chromate anion has a strong redox capacity that seriously threatens the ecological environment and human health. Cr can contaminate water and impart toxicity to aquatic species. Procambarus clarkii is an important food source that once represented a large proportion of the aquaculture industry due to its rapid reproduction and high economic value. However, there have been reports on the death of P. clarkii due to heavy metal pollution. The underlying mechanism regarding heavy metal toxicity was studied in this paper. The transcriptome data of hemocytes extracted from P. clarkii injected with Cr were analyzed by high-throughput sequencing and compared to the control group. In total, 48,128,748 clean reads were obtained in the treatment group and 56,480,556 clean reads were obtained in the control group. The reads were assembled using Trinity and the identified unigenes were then annotated. Then, 421 differentially-expressed genes (DEGs) were found, 170 of which were upregulated and 251 downregulated. Many of these genes were found to be related to glutathione metabolism and transportation. The glutathione metabolic pathway of P. clarkii was thus activated by Cr exposure to detoxify and maintain body function. Validation of DEGs with quantitative real-time PCR confirms the changes in gene expression. Thus, this study provides data supporting a glutathione-focused response of P. clarkii to exposure to heavy metals.

Characterization of the complete mitochondrial genome of Helice latimera and its phylogenetic implications in Brachyura.

Mitochondrial genomes (mitogenomes) help advance our learning of molecular evolution and phylogenetic relationships. The mitogenome of H. latimera is 16,246 bp in length, which typically contains 37 animal mitogenome genes consisting of 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes, as well as a control region. The AT content of H. latimera is 69.1%. The A + T skew of the mitogenome of H. latimera was slightly negative (-0.017). The size of Thirteen PCGs is from 162 bp to 1731 bp. Twenty-two tRNA genes ranged from 62 to 73 bp and were highly A + T biased. All tRNA genes owed a typical cloverleaf structure, not including the trnS1 gene lacking a dihydroxyuridine arm. One PCG, two rRNAs, and 12 of the tRNAs were rearranged compared to the pancrustacean gene order. Phylogenetic analysis revealed the locationt of H. latimera among the Varunidae family.

The red swamp crayfish, Procambarus clarkii cathepsin C, participates in the innate immune response to the viral and bacterial pathogens.

The cathepsin C, a lysosomal cysteine protease, involves the modulation of immune and inflammatory responses in living organisms. However, the knowledge on cathepsin C in red swamp crayfish (Procambarus clarkii), a freshwater crustacean with economic values, remained unclear. In the present study, we provide identification and molecular characterization of cathepsin C from P. clarkii. (Hereafter Pc-cathepsin C). The Pc-cathepsin C cDNA contained a 1356 bp open reading frame that encoded a protein of 451 amino acid residues. The deduced amino acid sequence comprised of cathepsin C exclusion domain and pept_C1 domain, and also catalytic residues (Cys248, His395 and Asn417). Analysis of the transcriptional patterns of the Pc-cathepsin C gene revealed that it was broadly distributed in various tissues of P. clarkii, and it was more abundant in the hepatopancreas and gut. Following a challenge with viral and bacterial pathogen-associated molecular patterns, the expression of Pc-cathepsin C was strongly enhanced at different time points. The knockdown of Pc-cathepsin C, altered the expression of immune-responsive genes, suggesting its immunoregulatory role in P. clarkii. This study has identified and provided the immunoregulatory function of Pc-cathepsin C, which will contribute to further investigation of the molecular mechanism of cathepsin C in crustaceans.

Proteomic analysis of differentially expressed proteins in the lipopolysaccharide-stimulated hepatopancreas of the freshwater crayfish Procambarus clarkii.

Procambarus clarkii is one of the most important aquatic invertebrates in China and has high commercial value. However, aquaculture has suffered great economic loss due to outbreaks of infectious diseases in P. clarkii. To identify red swamp crayfish related proteins involved in the response to bacterial infection, we analysed immune-related proteins following lipopolysaccharide (LPS) stimulation by quantitative proteomics. The proteome of the hepatopancreas of P. clarkii challenged with LPS and phosphate-buffered saline was analysed to evaluate the immune response. Based on liquid chromatography coupled with tandem mass spectrometry, 16 upregulated and 29 downregulated proteins were identified. A Gene Ontology analysis demonstrated 5 biological process, 11 cellular component, and 6 molecular function subcategories. The Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that the identified proteins were mainly involved in metabolism, phagosome, and ribosome. Real-time quantitative reverse transcription-PCR revealed that eight immune-related genes were upregulated after LPS stimulation compared to the control. Taken together, the data enhance our understanding of the immune response of crayfish to LPS.

Mitochondrial genome of the yellow catfish Pelteobagrus fulvidraco and insights into Bagridae phylogenetics.

The mitochondrial genome (mitogenome) can provide important information for understanding phylogenetic analysis and molecular evolution. Herein, we amplified the complete mitogenome sequence of Pelteobagrus fulvidraco. The mitogenome was 16,526 bp in length and included 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and a non-coding control region (D-loop). Both the organization and location of genes in the mitogenome were consistent with those from Siluriformes fishes previously published in GenBank. The phylogenetic relationships based on Bayesian inference (BI) and Maximum likelihood (ML) methods showed that P. fulvidraco has close relationships with Pelteobagrus eupogon and Tachysurus intermedius, suggesting that P. fulvidraco belongs to Tachysurus. This study provides evidence that Tachysurus, Pseudobagrus and Leiocassis do not form monophyly, but that these three genera form a monophyletic group. Our results provide reference for further phylogenetic research of the Bagridae species.

A non-mammalian Toll-like receptor 26 (TLR26)gene mediates innate immune responses in yellow catfish Pelteobagrus fulvidraco.

In this study, we identified a fish-specific Toll-like receptor (TLR) in Pelteobagrus fulvidraco, an economically important freshwater fish in China. This TLR, PfTLR26, was shown to be encoded by a 3084 bp open reading frame (ORF), producing a polypeptide 1027 amino acids in length. The PfTLR26 protein contains a signal peptide, eight leucine-rich repeat (LRR) domains, two LRR_TYP domains in the extracellular region, and a Toll/interleukin (IL)-1 receptor (TIR) domain in the cytoplasmic region, consistent with the characteristic TLR domain architecture. This predicted 117.1 kDa protein was highly homologous to those of other fish, with phylogenetic analysis revealing the closest relation to TLR26 of Ictalurus punctatus. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis showed that the PfTLR26 gene was expressed in all tissues tested, with the highest expression levels seen in the head kidney and blood, and the lowest seen in muscle. PfTLR26 exhibited significant upregulation in liver, spleen, head kidney, and blood at different time points following challenge with the common TLR agonists lipopolysaccharide (LPS) and polyriboinosinic polyribocytidylic acid (Poly I:C). Taken together, these results suggest that PfTLR26 may be an important component of the P. fulvidraco innate immune system, participating in the transduction of TLR signaling under pathogen stimulation.

Characterization and expression analysis of immune-related genes in the red swamp crayfish, Procambarus clarkii in response to lipopolysaccharide challenge.

To learn more about red swamp crayfish related genes in response to bacterial infections, we investigated immune-related genes induced by lipopolysaccharide (LPS) in the hepatopancreas using high-throughput sequencing method. In present the study, a total of 55,107 unigenes were identified, with an average length of 678 bp. A total of 2215 differentially expressed genes (DEGs) were found, including 669 up-regulated genes and 1546 down-regulated genes. The result of Gene ontology (GO) analysis revealed that 3017 DEGs were enriched in 19 biological process subcategories, 17 cellular component subcategories and 15 molecular function subcategories. The top 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that "ribosome" was the most abundant group, which had 34 DEGs. KEGG enrichment analysis identified several immune response pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results exhibited that several immune responsive genes were greatly up-regulated following LPS stimulation as observed in the results of high-throughput sequencing. Overall, this study provides new insight into the immune defense mechanisms of P. clarkii against LPS infection.

Transcriptome analysis of differential expressed genes in hepatopancreas of Procambarus clarkii challenged with peptidoglycan.

Procambarus clarkii is one of the most economically important species in Chinese aquaculture, and is widely cultured. Infection of P. clarkii populations with bacterial pathogens causes high mortality and great economic loss, therefore disease control is of significant economic importance. P. clarkii is a model system for studying immune responses in invertebrates, and its immune system consists solely of the innate response. In the present study, we examined gene expression related to immune function in P. clarkii in response to pathogen challenge. The transcriptome of hepatopancreas tissue from P. clarkii challenged with peptidoclycan (PGN) was analyzed and compared to control specimens. After assembly and annotation, 48,661 unigenes were identified with an average length of 671.54 bp. A total of 2533 differentially expressed genes (DEGs) were obtained, including 765 significantly up-regulated unigenes and 1757 significantly down-regulated unigenes. Gene ontology (GO) analysis demonstrated 19 biological process subcategories, 16 cellular component subcategories, and 17 molecular function subcategories that were enriched among these DEGs. Enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database revealed enrichment among immune responses pathways. Taken together, this study not only enriches the existing P. clarkii transcriptome database, but also elucidates immune responses of crayfish that are activated in response to PGN challenge.

Transcriptome-wide identification of differentially expressed genes in Procambarus clarkii in response to chromium challenge.

Because of the high protein content and rich meat quality of crayfish Procambarus clarkii, it has become widely popular in China in recent years and has a high economic value. When P. clarkii is stimulated by heavy metals, it reacts to oxidation. P. clarkii has evolved antioxidant defense systems, including antioxidant enzymes such as catalase (CAT). The hexavalent form of Cr (VI) is a pathogenic factor that is of particular concern in aqueous systems because of its great toxicity to living organisms. In this study, we characterized the transcriptome of P. clarkii using a RNA sequencing method and performed a comparison between K2Cr2O7-treated samples and controls. In total, 34,237 unigenes were annotated. We identified 5098 significantly differentially expressed genes (DEGs), including 2536 and 2562 were significantly up-regulated and down-regulated, respectively. In addition, quantitative real time-PCR (qRT-PCR) confirmed the up-regulation of a random selection of DEGs. Our results contribute to a more comprehensive understanding of the antioxidant defense system used by P. clarkii in response to heavy metal stress.

New insight into the molecular basis of Fe (III) stress responses of Procambarus clarkii by transcriptome analysis.

Iron in excess can have toxic effects on living organisms. In China, the freshwater crayfish Procambarus clarkii is a source of aquatic food with high-quality protein and has significant commercial value. P. clarkii shows oxidative stress on exposure to heavy metals, and antioxidant enzymes, such as ubiquitination enzymes and proteasomes, play important roles in oxidative stress. To understand the antioxidant defense system of P. clarkii, we analyzed the hepatopancreas transcriptomes of P. clarkii after stimulation with FeCl3. In total, 5199 differentially expressed genes (DEGs) were identified (2747 upregulated and 2452 downregulated). GO analysis revealed that these DEGs belonged to 16 cellular component, 16 molecular function, and 19 biological process subcategories. A total of 1069 DEGs were classified into 25 categories by using COG. Some antioxidant defense pathways, such as "Ubiquitin mediated proteolysis" and "Glutathione metabolism," were identified using KEGG. In addition, quantitative real time-PCR (qRT-PCR) substantiated the up-regulation of a random selection of DEGs including antioxidant and immune defense genes. We obtained information for P. clarkii transcriptome databases and new insights into the responses of P. clarkii hepatopancreas to heavy metals.

Transcriptomic analysis of immune-related genes in the lipopolysaccharide-stimulated hepatopancreas of the mudflat crab Helice tientsinensis.

The mudflat crab Helice tientsinensis is one of the most economically important aquaculture species in China. Nevertheless, it is susceptible to various diseases caused by viruses, bacteria and rickettsia-like organisms. A better understanding of the immune system and genes related to the responses to bacterial and viral infection is required. Herein, the hepatopancreas transcriptome of H. tientsinensis was analyzed by comparing control and lipopolysaccharide (LPS)-stimulated RNA-Seq data, yielding 91,885,038 bp and 13.78 Gb of clean reads. Following assembly and annotation, 93,207 unigenes with an average length of 883 bp were identified, of which 31,674 and 13,700 were annotated in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Following LPS, 4845 differentially expressed genes (DEGs) were identified, of which 2491 and 2354 were up- and down-regulated, respectively. To further investigate immune-related DEGs, KEGG enrichment analysis identified immune response pathways, most notably the peroxisome and Toll-like receptor signaling pathways. Quantitative real time-PCR (qRT-PCR) confirmed the up-regulation of a random selection of DEGs. This systematic transcriptomic analysis of the innate immune pathway in H. tientsinensis expands our understanding of the immune system in crabs.

De novo transcriptome assembly and analysis of differential gene expression following lipopolysaccharide challenge in Pelteobagrus fulvidraco.

The yellow catfish, Pelteobagrus fulvidraco, has been recognized as an important freshwater aquaculture species in Eastern and Southeast Asia. To gain a better understanding of the immune response in P. fulvidraco, we analyzed its transcriptome following stimulation with lipopolysaccharide (LPS). Phosphate buffer saline (PBS) was used as control. Following assembly and annotation, 72,152 unigenes with an average length of 1090 bp were identified. A total of 370 differentially expressed genes (DEGs) in the P. fulvidraco were observed at 12 h post LPS treatment, including 197 up-regulated genes and 173 down-regulated genes. Clusters of Orthologous Groups of proteins (KOG/COG) annotation demonstrated that a total of 18,819 unigenes classified into 26 categories. Gene ontology (GO) analysis revealed 20 biological process subcategories, 7 cellular component subcategories and 20 molecular function subcategories. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified immune responses pathways. Quantitative reverse transcription polymerase chain reaction measured the expression of 18 genes involved in the immune response. CXCL2-like chemokine (CXCL2), goose-type lysozyme (LYZ G), and cathepsin K (CTSK) were significantly up-regulated. This study enriches the P. fulvidraco transcriptome database and provides insight into the immune response of P. fulvidraco against infection.

The complete mitochondrial genome of Clostera anachoreta (Lepidoptera: Notodontidae) and phylogenetic implications for Noctuoidea species.

In this study, the complete mitochondrial genome (mitogenome) of Clostera anachoreta (Lepidoptera: Notodontidae) was sequenced. It comprises 15,456 base pairs (bp), including 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and one non-coding control region (CR), as found in other lepidopterans. Gene order is identical to that of typical lepidopterans. There are 15 intergenic spacers ranging from 2 to 49bp, and 9 overlapping regions ranging from 1 to 8bp, occurring throughout the genome. The CR is 347bp long. All PCGs are initiated by ATN codons. We found a typical gene rearrangement in C. anachoreta (tRNAMet-tRNAIle-tRNAGln), which is different from ancestral insects (tRNAIle-tRNAGln-tRNAMet). The gene rearrangement can be explained by a duplication/random loss model. Phylogenetic analyses indicate that C. anachoreta belongs to Notodontidae, and that the monophyly of Lepidopteran families is well supported.