Complete mitochondrial genome of Parasa sinica: New insights into the phylogeny of Limacodidae.

In this study, the whole mitochondrial genome (mitogenome) of Parasa sinica was sequenced. It contains 15,306 base pairs (bp), 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and one non-coding regulatory area (CR), all of which are shared by other lepidopterans. It follows the same gene order as ordinary lepidopterans. Further, out of these 37 genes, 23 are present on the heavy strand whereas the remaining 14 are located on the light strand. The A + T composition of the mitogenome is relatively high. Although P. sinica has a negative AT-skew and GC-skew, the GC-skew value is significantly lower than the AT-skew value. All PCGs, with the exception of CO1, carry the same start codon (ATN). All tRNAs exhibit the usual cloverleaf secondary structure. We identified the conserved motif "ATAGA + poly-T″ found in other lepidopteran insects at the beginning of the CR. We collected the concatenated PCGs sequences in the mitochondrial genome of 15 species of Zygaenoidea, with the sequences of Geometridae as outgroups, including P. sinica, and constructed phylogenetic trees using Bayesian inference (BI) and maximum likelihood (ML) methods. The monolineage of each superfamily is usually well supported. Based on phylogenetic analysis, P. sinica is a member of family Limacodidae, strongly supporting the monophyly of the Zygaenoidea groups.

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.

Comparative transcriptome analysis of eyes reveals the adaptive mechanism of mantis shrimp (oratosquilla oratoria) induced by a dark environment.

The light-dark cycle significantly impacts the growth and development of animals. Mantis shrimps (Oratosquilla oratoria) receive light through their complex photoreceptors. To reveal the adaptive expression mechanism of the mantis shrimp induced in a dark environment, we performed comparative transcriptome analysis with O. oratoria cultured in a light environment (Oo-L) as the control group and O. oratoria cultured in a dark environment (Oo-D) as the experimental group. In the screening of differentially expressed genes (DEGs) between the Oo-L and Oo-D groups, a total of 88 DEGs with |log2FC| > 1 and FDR < 0.05 were identified, of which 78 were upregulated and 10 were downregulated. Then, FBP1 and Pepck were downregulated in the gluconeogenesis pathway, and MKNK2 was upregulated in the MAPK classical pathway, which promoted cell proliferation and differentiation, indicating that the activity of mantis shrimp was slowed and the metabolic rate decreases in the dark environment. As a result, the energy was saved for its growth and development. At the same time, we performed gene set enrichment analysis (GSEA) on all DEGs. In the KEGG pathway analysis, each metabolic pathway in the dark environment showed a slowing trend. GO was enriched in biological processes such as eye development, sensory perception and sensory organ development. The study showed that mantis shrimp slowed down metabolism in the dark, while the role of sensory organs prominent. It provides important information for further understanding the energy metabolism and has great significance to study the physiology of mantis shrimp in dark environment.

Whole genome evaluation analysis and preliminary Assembly of Oratosquilla oratoria (Stomatopoda: Squillidae).

BACKGROUND: As the dominant species of Stomatopoda, Oratosquilla oratoria has not been fully cultivated artificially, and the fishery production mainly depends on marine fishing. Due to the lack of stomatopod genome, the development of molecular breeding of mantis shrimps still lags behind. METHODS AND RESULTS: A survey analysis was performed to obtain the genome size, GC content and heterozygosity ratio in order to provide a fundation for subsequent whole-genome sequencing. The results showed that the estimated genome size of the O. oratoria was about 2.56 G, and the heterozygosity ratio was 1.81%, indicating that it is a complex genome. Then the sequencing data was preliminarily assembled with k-mer = 51 by SOAPdenovo software to obtain a genome size of 3.01G and GC content of 40.37%. According to ReapeatMasker and RepeatModerler analysis, the percentage of repeats in O. oratoria was 45.23% in the total genome, similar to 44% in Survey analysis. The MISA tool was used to analyze the simple sequence repeat (SSR) characteristics of genome sequences including Oratosquilla oratoria, Macrobrachium nipponense, Fenneropenaeus chinensis, Eriocheir japonica sinensis, Scylla paramamosain and Paralithodes platypus. All crustacean genomes showed similar SSRs characteristics, with the highest proportion of di-nucleotide repeat sequences. And AC/GT and AGG/CCT repeats were the main types of di-nucleotide and tri-nucleotide repeats in O. oratoria. CONCLUSION: This study provided a reference for the genome assembly and annotation of the O. oratoria, and also provided a theoretical basis for the development of molecular markers of O. oratoria.

Time-course transcriptome analysis reveals regulation of Arabidopsis seed dormancy by the transcription factors WOX11/12.

The induction of seed dormancy and its release involve a finely regulated genetic program controlled by various environmental and developmental cues that are critical for plant survival and population expansion. Light plays a key role in seed dormancy and germination, but the molecular mechanisms underlying the control of dormancy are unclear. In the present study, high-resolution temporal RNA-seq in Arabidopsis identified WOX11 as encoding a hub transcription factor during the seed dormancy induction and release stages. This gene might have evolved from gymnosperms and expanded in angiosperms with highly conserved expression patterns in seeds. WOX11 and its homolog WOX12 were highly expressed from 2 d after pollination, and mRNA abundance was greatly increased during the seed dormancy induction and release stages. Further, we found that WOX11 plays a role in the regulation of seed dormancy downstream of phytochrome B (PHYB)-mediated red-light signaling during the induction stage, indicating that WOX11/12 are newly identified components of red-light signaling transduction. Taken together, our results suggest that WOX11/12-mediated PHYB signaling regulates seed dormancy in Arabidopsis, and provide insights into the developmental regulation and evolutionary adaptation of plants to changes in the light environment.

Differentially expressed genes in head kidney of Pelteobagrus fulvidraco following Vibrio cholerae challenge.

The yellow catfish (Pelteobagrus fulvidraco) is a freshwater fish with high economic value in eastern China. Nevertheless, pathogens causing bacterial diseases in P. fulvidraco have brought about huge economic loss and high mortality in artificial aquaculture. For disease control, it is critical to further understand the immune system of yellow catfish and immune-related genes with which they respond to pathogenic infections. In this study, high-throughput sequencing methods were used to analyze the transcriptomic spectrum of the head kidney from P. fulvidraco challenged by Vibrio cholera. A total of 45,544 unique transcript fragments (unigenes) were acquired after assembly and annotation, with an average length of 1,373 bp. Additionally, 674 differentially expressed genes (DEGs) were identified after stimulation with V. cholerae, 353 and 321 genes were identified as remarkably up- or downregulated, respectively. To further study the immune-related DEGs, we performed KEGG enrichment and GO enrichment. The results showed gene regulation of response to stimulus, immune response, immune system progress, response to external stimuli and cellular response to stimuli. Analysis of KEGG enrichment is important to identify chief immune related pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results indicated 10 immune response genes that were found to be upregulated compared to a control group after 6 h of V. cholerae challenging. In summary, the results of our study are helpful to determine the defense mechanisms and immune system responses of yellow catfish in reaction to bacterial challenges.

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.

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.

Differentially expressed genes involved in immune pathways from yellowhead catfish (Tachysurus fulvidraco) after poly (I:C) challenge.

Yellowhead catfish (Tachysurus fulvidraco) is an important aquaculture fish species in China with a high market value. Infectious diseases pose serious threats in farmed fish species, and although vaccines can prevent certain infections, they rely on potent adjuvants. In this study, we analyzed the transcriptomic profiles of spleens from poly (I:C)-treated T. fulvidraco. We obtained 46,362,922 reads corresponding to 490,926 transcripts and 318,059 genes. Gene annotation using different databases and subsequent differential gene expression analyses led to the identification of 5587 differentially expressed genes (DEGs), of which 2473 were up-regulated and 3114 were down-regulated in poly (I:C)-treated fish. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs revealed the significant dysregulation of immune- and cancer-related genes in the spleens of poly (I:C)-treated fish. Notably, several components of JAK-STAT, MAPK, and p53 signaling pathways were significantly dysregulated in response to poly (I:C) treatment. Quantitative real-time PCR (qRT-PCR) analysis of 11 randomly selected immune response genes confirmed the reliability of our findings. In conclusion, our findings provide novel insight into the immune responses of T. fulvidraco and suggest that poly (I:C) may represent a promising adjuvant of fish vaccines.

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.

Chromosome-level genome assembly of Paralithodes platypus provides insights into evolution and adaptation of king crabs.

The blue king crab, Paralithodes platypus, which belongs to the family Lithodidae, is a commercially and ecologically important species. However, a high-quality reference genome for the king crab has not yet been reported. Here, we assembled the first chromosome-level blue king crab genome, which contains 104 chromosomes and an N50 length of 51.15 Mb. Furthermore, we determined that the large genome size can be attributed to the insertion of long interspersed nuclear elements and long tandem repeats. Genome assembly assessment showed that 96.54% of the assembled transcripts could be aligned to the assembled genome. Phylogenetic analysis showed the blue king crab to have a close relationship with the Eubrachyura crabs, from which it diverged 272.5 million years ago. Population history analyses indicated that the effective population of the blue king crab declined sharply and then gradually increased from the Cretaceous and Neogene periods, respectively. Furthermore, gene families related to developmental pathways, steroid and thyroid hormone synthesis, and inflammatory regulation were expanded in the genome, suggesting that these genes contributed substantially to the environmental adaptation and unique body plan evolution of the blue king crab. The high-quality reference genome reported here provides a solid molecular basis for further study of the blue king crab's development and environmental adaptation.

Comparative Mitochondrial Genome Analyses of Sesarmid and Other Brachyuran Crabs Reveal Gene Rearrangements and Phylogeny.

Mitochondrial genomes (mitogenomes) are important for understanding molecular evolution and phylogenetic relationships. The complete mitogenome of Perisesarma bidens was determined, which is 15,641 bp in length. The A + T content of P. bidens mitogenome was 74.81%. The AT skew was slightly negative (-0.021). The 22 tRNAs ranged from 65 to 73 bp and were highly A + T biased. All tRNA genes had typical cloverleaf structures, except for the trnS1 gene, which lacked a dihydrouridine (DHU) arm. The gene order within the P. bidens mitogenome was identical to the pancrustacean ground pattern, except for the translocation of the trnH. Additionally, the gene order of trnI-trnQ-trnM in pancrustacean ground pattern became trnQ-trnI-trnM in P. bidens. Phylogenetic analyses supported the inclusion of P. bidens in Sesarmidae and the promotion of Sesarminae to Sesarmidae. The results will help us to better understand the status and evolutionary history of Grapsoidea crabs.

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.

A novel modulation of physiological regulation in cultured Chinese mitten crab (Eriocheir japonica sinensis) in response to consistent salinity changes.

The life history of the Chinese mitten crab (Eriocheir japonica sinensis) includes two migrations: a feeding migration and a reproductive migration. Ambient salinity is one of the most critical factors during migration. In this study, the salinity adaptation mechanism of Chinese mitten crabs was simulated using continuous salinity changes. The expression of six key genes [Na+/K+-ATPase α subunit (NAK-α), V-type H+-ATPase subunit A (VHA-A), Zinc transporter (ZnT), Cl- channel protein 2 (CLCN2), ubiquitin/ribosomal S27 fusionprotein (S27), and glutathione S-transferase (GST)] and the activities of three enzymes [Na+/K+-ATPase (NAK), V-type H+-ATPase (VHA), and glutathione S-transferase (GST)] were evaluated in ten groups exposed to a range of salinity changes during mariculture based on the transcriptome data obtained from short term salinity-induced crabs (ES) compared to control group in freshwater crabs (EF). The results revealed that different genes exhibited different roles in physiological regulation. In total, 3,599 unigenes were significantly and differentially expressed in a comparison between the EF and ES treatments. A novel modulation of gene expression and the corresponding enzyme activity of NAK and VHA exhibited similar patterns. As genes related to osmoregulation, NAK and VHA showed similar patterns of both gene expression and enzyme activity in mariculture. During the gradual change in salinity from 0‰ to 25‰ and back to 0‰, the gene expression and enzyme activities of NAK and VHA initially increased (0‰ â†’ 10‰), weakened (10‰ â†’ 20‰) and then increased again (20‰ â†’ 25‰ â†’ 0‰). S27 could serve as a reference gene in the expression analysis of Chinese mitten crabs under salinity stress. ZnT and CLCN2 were involved in osmoregulation as functional proteins. Our findings provide insights into the regulation mechanisms employed during the migration of the Chinese mitten crab.

Transcriptome Analysis Reveals the Tolerance Mechanism of Mantis Shrimp (Oratosquilla oratoria) under a Lipopolysaccharide Challenge.

Lipopolysaccharide (LPS), a major cell wall component of Gram-negative bacteria, is considered to lead to some disease development in commercial crustaceans. However, mantis shrimps Oratosquilla oratoria (Crustacea: Stomatopoda) have a strong vitality and ability to resist disease. To study the tolerance mechanism of mantis shrimp, transcriptome analyses were conducted in hepatopancreas of O. oratoria under LPS challenge investigation. Totally, 84 547 044 clean reads were obtained from transcriptomes (43 159 230 in OP (control), 41 387 814 in OL (treatment), respectively). Unigenes, the longest transcript of each gene, with a total length of 68 318 880 bp and the total number of 100 978 were obtained. 8369 (8.28%) of unigenes were successfully annotated in all databases and 54 888 (54.35%) were annotated in at least one database. Finally, 1012 differentially expressed genes (DEGs) including 439 and 573 showed significantly upregulated and downregulated were determined between OL and OP, respectively. Moreover, those DEGs only expressed in OL or OP accounted for 8.99%. The functional classification based on GO and KEGG indicated that the common enrichment categories for the DEGs are "amino sugar metabolic" and "cellular homeostasis" and that the progress of nutrient metabolic and homeostasis in cells is important in facing variable environmental conditions. Protein-protein interaction analysis elucidated proteins, ß-actin (ACTB_G1), T-complex protein subunits (TCPs), heat shock proteins (HSPs), hydroxysteroid dehydrogenase-like protein 2 (HSDL2), kinesin family member 5 (KIF5), methylglutaconyl-CoA hydratase (AUH), and myosin heavy chain (MYH) may play key roles in response to an LPS challenge. This study laid a foundation to further investigate the possible adaptation way that O. oratoria survives in a bacterial challenge.

Chromosome-level genome assembly reveals the unique genome evolution of the swimming crab (Portunus trituberculatus).

BACKGROUND: The swimming crab, Portunus trituberculatus, is an important commercial species in China and is widely distributed in the coastal waters of Asia-Pacific countries. Despite increasing interest in swimming crab research, a high-quality chromosome-level genome is still lacking. FINDINGS: Here, we assembled the first chromosome-level reference genome of P. trituberculatus by combining the short reads, Nanopore long reads, and Hi-C data. The genome assembly size was 1.00 Gb with a contig N50 length of 4.12 Mb. In addition, BUSCO assessment indicated that 94.7% of core eukaryotic genes were present in the genome assembly. Approximately 54.52% of the genome was identified as repetitive sequences, with a total of 16,796 annotated protein-coding genes. In addition, we anchored contigs into chromosomes and identified 50 chromosomes with an N50 length of 21.80 Mb by Hi-C technology. CONCLUSIONS: We anticipate that this chromosome-level assembly of the P. trituberculatus genome will not only promote study of basic development and evolution but also provide important resources for swimming crab reproduction.

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.

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.