Genome-wide QTL and eQTL mapping reveal genes associated with growth rate trait of the Pacific white shrimp (Litopenaeus vannamei).

BACKGROUND: Growth rate is a crucial economic trait for farmed animals, but the genetic regulation of this trait is largely unknown in non-model organisms such as shrimp. RESULTS: In this study, we performed genome-wide phenotypic quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) mapping analyses to identify genes affecting the growth rate of Pacific white shrimp (Litopenaeus vannamei), which is the most commercially-farmed crustacean worldwide. We used RNA-sequencing of 268 individuals in a mapping population, and subsequently validated our findings through gene silencing and shrimp growth experiments. We constructed a high-density genetic linkage map comprising 5533 markers spanning 44 linkage groups, with a total distance of 6205.75 cM and an average marker interval of 1.12 cM. Our analyses identified 11 QTLs significantly correlated with growth rate, and 117,525 eQTLs. By integrating QTL and eQTL data, we identified a gene (metalloreductase STEAP4) highly associated with shrimp growth rate. RNA interference (RNAi) analysis and growth experiments confirmed that STEAP4 was significantly correlated with growth rate in L. vannamei. CONCLUSIONS: Our results indicate that the comprehensive analysis of QTL and eQTL can effectively identify genes involved in complex animal traits. This is important for marker-assisted selection (MAS) of animals. Our work contributes to the development of shrimp breeding and available genetic resources.

Genome-wide analysis of ATP-binding cassette (ABC) transporter in Penaeus vannamei and identification of two ABC genes involved in immune defense against Vibrio parahaemolyticus by affecting NF-κB signaling pathway.

The ATP-binding cassette (ABC) transporters have crucial roles in various biological processes such as growth, development and immune defense in eukaryotes. However, the roles of ABC transporters in the immune system of crustaceans remain elusive. In this study, 38 ABC genes were systematically identified and characterized in Penaeus vannamei. Bioinformation analysis revealed that PvABC genes were categorized into ABC A-H eight subfamilies with 17 full-transporters, 11 half transporters and 10 soluble proteins, and multiple immunity-related cis-elements were found in gene promoter regions. Expression analysis showed that most PvABC genes were widely and highly expressed in immune-related tissues and responded to the stimulation of Vibrio parahaemolyticus. To investigate whether PvABC genes mediated innate immunity, PvABCC5, PvABCF1 and PvABCB4 were selected for dsRNA interference experiment. Knockdown of PvABCF1 and PvABCC5 not PvABCB4 increased the cumulative mortality of P. vannamei and bacterial loads in hepatopancreas after infection with V. parahaemolyticus. Further analysis showed that the PvABCF1 and PvABCC5 knockdown decreased expression levels of NF-κB pathway genes and antimicrobial peptides (AMPs). Collectively, these findings indicated that PvABCF1 and PvABCC5 might restrict V. parahaemolyticus challenge by positively regulating NF-κB pathway and then promoting the expression of AMPs, which would contribute to overall understand the function of ABC genes in innate immunity of invertebrates.

Cloning, Identification, and Functional Analysis of the Foxl2 Gene in Procambarus clarkii.

Procambarus clarkii is the most widely distributed freshwater shrimp in China, with important economic value and great potential for development. The forkheadboxL2 (Foxl2) gene has been found to be involved in the reproductive development of many crustaceans. To understand the role of the Foxl2 gene in the gonad development of P. clarkii, we designed CDS-specific primers for the P. clarkii Foxl2 (PcFoxl2) gene and cloned its CDS sequence using RT-PCR. The nucleotide and protein sequence information was then analyzed through bioinformatics analysis. The expression and subcellular localization of PcFoxl2 in various tissues were detected using qRT-PCR and in situ hybridization. The effects of PcFoxl2 knockdown on gonad development were investigated using RNA interference. The results showed that the CDS length of the PcFoxl2 gene was 1614 bp and encoded 537 amino acids. Protein sequence comparison and phylogenetic analysis showed that PcFoxl2 was the closest relative to Crayfish. qRT-PCR analysis indicated that the expression level of PcFoxl2 in the testis was significantly higher (>40 fold) than that in the ovary (p < 0.01). The in situ hybridization results showed that PcFoxl2 was expressed in both the cytoplasm and the nucleus of egg cells, and that the expression was strongest in egg cells at the early stage of yolk synthesis, while weak in the secondary oocytes. The positive signal was strongest in the spermatocyte nucleolus, while only a trace signal was observed in the cytoplasm. After interfering with the PcFoxl2 gene using dsRNA, the expression of PcFoxl2 in the RNA interference group was significantly lower than that in the control group, and this interference effect lasted for one week. Moreover, the gonad index of the experimental group was significantly lower than that of the control group (p < 0.05) after 10 days of P. clarkii cultivation following PcFoxl2 knockdown. The expression levels of the nanos and S3a genes, which are related to gonad development, decreased significantly after PcFoxl2 gene interference. The results suggest that the Foxl2 gene is involved in the growth and development of gonads, particularly in the development of testis, and is related to the early development of oocytes. This study provides a theoretical basis for the artificial breeding of P. clarkii.

Integrated Analysis of DNA Methylome and Transcriptome Reveals Epigenetic Regulation of Cold Tolerance in Litopenaeus vannamei.

DNA methylation is an important epigenetic modification that has been shown to be associated with responses to non-biological stressors. However, there is currently no research on DNA methylation in response to environmental signals in shrimp. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles and differentially expressed genes between two strains of Litopenaeus vannamei with significantly different cold tolerance through whole genome bisulfite sequencing (WGBS) and transcriptome sequencing. Between Lv-C and Lv-T (constant temperature of 28 °C and low temperatures of 18 °C and 10 °C) under cytosine-guanine (CG) environments, 39,100 differentially methylated regions (DMRs) were identified, corresponding to 9302 DMR-related genes (DMRGs). The DMRs were mainly located in the gene body (exons and introns). Gene Ontology (GO) analysis showed that these DMRGs were significantly enriched in cell parts, catalytic activity, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed significant enrichment of these DMRGs in pathways such as proteasome (ko03050), oxidative phosphorylation (ko00190), mTOR signaling pathway (ko04150), fatty acid metabolism (ko01212), and fatty acid degradation (ko00071). The comprehensive results suggested that L. vannamei mainly regulates gene expression in response to low temperatures through hypermethylation or demethylation of some genes involved in thermogenesis, glycolysis, the autophagy pathway, the peroxisome, and drug metabolism pathways. These results provide important clues for studying DNA methylation patterns and identifying cold tolerance genes in shrimp.

The complete mitochondrial genome of Scatophagus argus (Perciformes, Scatophagidae) from Beibu Bay and phylogenetic analysis.

The spotted scat (Scatophagus argus, Linnaeus, 1766) is a subtropical fish that is widely distributed in the coastal waters of Indo-Pacific. Here, we report the complete mitochondrial genome of S. argus. The mitogenome is 16,783 base pairs (56.0% A + T content) in length and consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a 1007 bp D-loop region. Phylogenetic analysis showed that the relationship between S. argus and Selenotoca multifasciata was close.

The complete mitochondrial genome of the mudskipper, Boleophthalmus pectinirostris (Gobiiformes, Oxudercidae) from Beibu Bay.

The mudskipper, Boleophthalmus pectinirostris (B. pectinirostris), is an amphibious fish that lives in the intertidal mudflats. It is a cultured economic fish with nutritional and pharmacological value. Here, we report the complete mitochondrial genome sequence of B. pectinirostris, which is 17,111 base pairs (55.3% A + T content) in length and consists of 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a 1453 bp D-loop region. The complete mitochondrial genome of B. pectinirostris will provide useful genetic information for future phylogenetic and taxonomic classification of B. pectinirostris.

The complete mitochondrial genome of yellowfin seabream, Acanthopagrus latus (Percoiformes, Sparidae) from Beibu Bay.

The yellowfin seabream, Acanthopagrus latus Houttuyn 1782, is a commercially and ecologically important species and a good model for studies of sexual differentiation. In this study, the complete mitochondrial genome of A. latus has been determined, which is 16,635 base pairs (54.3% A + T content) in length and consists of 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and a 948 bp D-loop region. The phylogenetic analyses showed that A. latus has a close relationship with Acanthopagrus schlegelii Bleeker 1854.

Single-Cell Ribonucleic Acid Sequencing Clarifies Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus Vannamei).

To characterize the cold tolerance mechanism of the Pacific white shrimp (Litopenaeus vannamei), we performed single-cell RNA sequencing (scRNA-seq) of ∼5185 hepatopancreas cells from cold-tolerant (Lv-T) and common (Lv-C) L. vannamei at preferred and low temperatures (28°C and 10°C, respectively). The cells fell into 10 clusters and 4 cell types: embryonic, resorptive, blister-like, and fibrillar. We identified differentially expressed genes between Lv-T and Lv-C, which were mainly associated with the terms "immune system," "cytoskeleton," "antioxidant system," "digestive enzyme," and "detoxification," as well as the pathways "metabolic pathways of oxidative phosphorylation," "metabolism of xenobiotics by cytochrome P450," "chemical carcinogenesis," "drug metabolism-cytochrome P450," and "fatty acid metabolism." Reconstruction of fibrillar cell trajectories showed that, under low temperature stress, hepatopancreas cells had two distinct fates, cell fate 1 and cell fate 2. Cell fate 1 was mainly involved in signal transduction and sensory organ development. Cell fate 2 was mainly involved in metabolic processes. This study preliminarily clarifies the molecular mechanisms underlying cold tolerance in L. vannamei, which will be useful for the breeding of shrimp with greater cold tolerance.

Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei).

BACKGROUND: Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available. RESULTS: In this study, we constructed a high-density genetic map of the Pacific white shrimp (Litopenaeus vannamei) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41-8.46% of the phenotypic variance in L. vannamei survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct L. vannamei families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance. CONCLUSIONS: In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.

Construction of a High-Density Genetic Map and Identification of Quantitative Trait Loci for Nitrite Tolerance in the Pacific White Shrimp (Litopenaeus vannamei).

Nitrite is a major environmental toxin in aquaculture systems that disrupts multiple physiological functions in aquatic animals. Although nitrite tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available. In this study, we constructed a high-density genetic map of Litopenaeus vannamei with 17,242 single nucleotide polymorphism markers spanning 6,828.06 centimorgans (cM), with an average distance of 0.4 cM between adjacent markers on 44 linkage groups (LGs). Using this genetic map, we identified two markers associated with nitrite tolerance. We then sequenced the transcriptomes of the most nitrite-tolerant and nitrite-sensitive individuals from each of four genetically distinct L. vannamei families (LV-I-4). We found 2,002, 1,983, 1,954, and 1,867 differentially expressed genes in families LV-1, LV-2, LV-3, and LV-4, respectively. By integrating QTL and transcriptomics analyses, we identified a candidate gene associated with nitrite tolerance. This gene was annotated as solute carrier family 26 member 6 (SLC26A6). RNA interference (RNAi) analysis demonstrated that SLC26A6 was critical for nitrite tolerance in L. vannamei. The present study increases our understanding of the molecular mechanisms underlying nitrite tolerance in shrimp and provides a basis for molecular-marker-assisted shrimp breeding.

Penaeus stylirostris densovirus proteins CP and NS1 interact with peritrophin of Litopenaeus vannamei.

The Penaeus stylirostris densovirus (PstDNV) is a major virus of shrimps that severely harms the shrimp farming industry. Peritrophin is a peritrophic membrane protein with chitin binding activity. To examine the roles of peritrophin in viral infection, we used yeast two-hybrid to analyze the interaction between the Pacific white shrimp (Litopenaeus vannamei) peritrophin and PstDNV proteins (CP, NS1 and NS2). The yeast two-hybrid results showed that NS1 and peritrophin had an interaction, CP and peritrophin had an interaction as well, and NS2 had no interaction with peritrophin. We validated the interactions with GST pull-down assays. We then conducted RNA interference and qRT-PCR. The results showed that when pre-injection of dsRNA-peritrophin, the quantity of PstDNV in the shrimps injected with viruses was significantly lower than in the control group (P < 0.01), indicating the viral infection was decreased when the peritrophin gene expression was inhibited. The results indicated that peritrophin of L. vannamei participated in the PstDNV infection.

The Penaeus stylirostris densovirus capsid protein interacts with the Litopenaeus vannamei BCCIP protein.

Viral capsid proteins play an important role in the viral infection process. To identify the cellular proteins in shrimp that interact with the Penaeus stylirostris densovirus capsid protein (PstDNV-CP), we constructed a yeast two-hybrid (Y2H) cDNA library of the muscle tissue of Litopenaeus vannamei, and hybridized the bait vector pGBKT7-CP with this library. Cloning and sequencing showed that the shrimp protein interacting with PstDNV-CP was a homolog of BRCA2 and CDKN1A(p21)-interacting protein (BCCIP). We named this protein L. vannamei BCCIP (LvBCCIP). Further analysis showed that LvBCCIP interacted with L. vannamei calmodulin (LvCaM). We validated the interactions between PstDNV-CP and LvBCCIP, and between LvBCCIP and LvCaM, with GST pulldown assays. The gene expression of LvBCCIP increased significantly after PstDNV challenge. In addition, the PstDNV titer of PstDNV-challenged shrimp was significantly reduced after LvBCCIP expression was inhibited using double-stranded RNA (dsRNA) interference. These results indicated that LvBCCIP is critical to PstDNV pathogenesis in L. vannamei. Interestingly, the growth rate of L. vannamei was significantly reduced when LvBCCIP gene expression was silenced, indicating that LvBCCIP may also be associated with growth regulation in L. vannamei. Thus, the interaction between PstDNV-CP and LvBCCIP might explain why PstDNV infection leads to runt-deformity syndrome in shrimp.

mRNA and microRNA transcriptomics analyses in intermuscular bones of two carp species, rice flower carp (Cyprinus carpio var. Quanzhounensis) and Jian carp (Cyprinus carpio var. Jian).

The rice flower carp (Cyprinus carpio var. Quanzhounensis) is a bony fish (superclass Osteichthyes), with very soft bones that is a significant feature unlike the other carp species. In this study, we analyzed the mRNA and microRNA (miRNA) transcriptomes in the intermuscular bones of rice flower carp and Jian carp (Cyprinus carpio var. Jian) (a typical member of common carp in China), using Illumina RNA sequencing. We identified 55,340 genes (including 47,541 known genes and 8231 predicted new genes) and 662 miRNAs (including 595 known miRNAs and 67 novel miRNAs) in the two species. By comparing the transcriptomes of the two species, we identified 1523 differentially expressed genes (DEGs) (including 576 up - and 947 downregulated DEGs) and 352 differentially expressed miRNAs (DEMs) (including 85 up- and 267 downregulated DEMs). According to the Gene Ontology (GO) annotation, 7 DEGs and 12 DEMs were found to be involved in the regulation of bone mineralization. The results of this study improve our understanding of the mRNA and miRNA profiles of carp bones, particularly those of the rice flower carp.

The Penaeus stylirostris densovirus capsid interacts with Litopenaeus vannamei troponin I.

The Penaeus stylirostris densovirus (PstDNV) (also known as infectious hypodermal and hematopoietic necrosis virus, IHHNV), a very small DNA virus, is a major shrimp pathogen. The PstDNV genome encodes only two nonstructural proteins and one capsid protein. This virus is thus an ideal, simple model for the investigation of virus-host interactions. To explore the role of the PstDNV capsid in viral infections, a yeast two-hybrid (Y2H) cDNA library was constructed based on Pacific white shrimp, Litopenaeus vannamei mRNA. The Y2H library was then screened, using the PstDNV capsid protein as bait. We identified a host protein that interacted strongly with the PstDNV capsid as L. vannamei troponin I (LvTnI). An in vitro co-immunoprecipitation experiment further supported this interaction. In addition, an in vivo neutralization experiment showed that the vaccination with anti-LvTnI significantly reduced PstDNV copies in PstDNV-challenged shrimp, indicating that the interaction between the PstDNV capsid and cellular LvTnI is essential for PstDNV infection. This result has important implications for our understanding of the mechanisms by which PstDNV infects shrimp.

Single-molecule long-read sequencing facilitates shrimp transcriptome research.

Although shrimp are of great economic importance, few full-length shrimp transcriptomes are available. Here, we used Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology to generate transcripts from the Pacific white shrimp (Litopenaeus vannamei). We obtained 322,600 full-length non-chimeric reads, from which we generated 51,367 high-quality unique full-length transcripts. We corrected errors in the SMRT sequences by comparison with Illumina-produced short reads. We successfully annotated 81.72% of all unique SMRT transcripts against the NCBI non-redundant database, 58.63% against Swiss-Prot, 45.38% against Gene Ontology, 32.57% against Clusters of Orthologous Groups of proteins (COG), and 47.83% against Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Across all transcripts, we identified 3,958 long non-coding RNAs (lncRNAs) and 80,650 simple sequence repeats (SSRs). Our study provides a rich set of full-length cDNA sequences for L. vannamei, which will greatly facilitate shrimp transcriptome research.

Proteomic Responses Under Cold Stress Reveal Unique Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus vannamei).

The Pacific white shrimp (Litopenaeus vannamei), one of the most widely cultured shrimp species in the world, often suffers from cold stress. To understand the molecular mechanism of cold tolerance in Pacific white shrimp, we conducted a proteomic analysis on two contrasting shrimp cultivars, namely, cold-tolerant Guihai2 (GH2) and cold-sensitive Guihai1 (GH1), under normal temperature (28°C), under cold stress (16°C), and during recovery to 28°C. In total, 3,349 proteins were identified, among which 2,736 proteins were quantified. Based on gene ontology annotations, differentially expressed proteins largely belonged to biological processes, cellular components, and molecular functions. KEGG pathway annotations indicated that the main changes were observed in the lysosome, ribosomes, and oxidative phosphorylation. Subcellular localization analysis showed a significant increase in proteins present in cytosol, extracellular regions, and mitochondria. Combining enrichment-based clustering analysis and qRT-PCR analysis, we found that glutathione S-transferase, zinc proteinase, m7GpppX diphosphatase, AP2 transcription complex, and zinc-finger transcription factors played a major role in the cold stress response in Pacific white shrimp. Moreover, structure proteins, including different types of lectin and DAPPUDRAFT, were indispensable for cold stress tolerance of the Pacific white shrimp. Results indicate the molecular mechanisms of the Pacific white shrimp in response to cold stress and provide new insight into breeding new cultivars with increased cold tolerance.

Identification of microRNAs involved in cold adaptation of Litopenaeus vannamei by high-throughput sequencing.

The Litopenaeus vannamei (L. vannamei) is one of the most widely cultured shrimp species in the world, with low temperature being one of the most serious threats to its growth and survival. To examine the potential regulatory mechanism of cold adaptation, we conducted a microRNAs (miRNAs) analysis on the hepatopancreas of L. vannamei under normal temperature 28 °C (M28), cold acclimation 16 °C for 6 days (M16), and recovered under normal temperature (MR). In total 14,754,823, 14,945,246 and 15,880,093 raw reads representing 10,690,259, 8,587,144, and 11,512,941 unique sequences of 18-32 nt length were obtained from the M28, M16 and MR libraries, respectively. After comparing the miRNA sequences with the miRBase database, 68 known mature miRNAs and 47 novel miRNAs were identified. Expression analysis showed that 34 miRNAs were significantly differential expressed in response to cold adaptation. Compared to the M28 library, 21 miRNAs were upregulated and 13 miRNAs were downregulated significantly in the M16 library. After recovery to normal temperature, there are 16 miRNAs upregulated and 15 miRNAs downregulated significantly compared to M28 library. Then, five significantly differential expressed miRNAs under cold acclimation including three known miRNAs (mja-miR-6491, mja-miR-6494, and Bta-miR-2478) and two newly-identified miRNAs (novel_68 and novel_5) were selected for validation by RT-qPCR in the hepatopancreas and muscle tissues of cold treated shrimps. The expression trend of most the miRNAs from RT-qPCR were consistent with the next-generation sequencing data. Further, the Gene Ontology (GO) annotation showed that the metabolic process GO term was significantly enriched with target genes of the differentially expressed miRNAs. Additionally, KEGG pathway analysis suggested that the fatty acid degradation and glycerolipid metabolism pathways etc. are significantly enriched with the target genes. These findings may contribute to a better understanding of the molecular mechanisms governing the responses to low temperature in L. vannamei.

Molecular cloning, expression, promoter analysis and functional characterization of a new Crustin from Litopenaeus vannamei.

Antimicrobial peptides (AMPs) are the most important players in the innate immune system, providing a principal first-line of defense against the invading pathogens. Crustin, a type of whey acidic protein (WAP) domain-containing and cationic cysteine-rich AMP, can function in a protease inhibition or an effector molecule manner. In the present study, a new Crustin was cloned and identified from Pacific white shrimp Litopenaeus vannamei and designated as LvCrustinA. The full-length cDNA of LvCrustinA was 687 bp, with a 519 bp open reading frame (ORF) that encoded a peptide of 172 amino acids. Domain analysis indicated that LvCrustinA contained a Glycine-rich region in the N-terminal and a single WAP domain within eight cysteines in the C-terminal. The 5' upstream regulatory sequence of 1249 bp (promoter) was obtained using a genome walking method, and it contained several conserved transcription factors binding motifs including NF-κB, AP-1 and STAT (Signal transducers and activators of transcription). Dual-reporter assay showed that NF-κB transcription factors LvDorsal and LvRelish, and AP-1 transcription factor Lvc-Jun could up-regulate the promoter activity of LvCrustinA, suggesting that NF-κB and JNK-c-Jun pathways could be involved in regulating the expression of LvCrustinA. Moreover, LvCrustinA was abundantly expressed in immune related tissues such as gill, hemocyte and epithelium, and its expression was up-regulated in response to Vibrio parahaemolyticus and White spot syndrome virus (WSSV) challenges in gill tissue, suggesting that LvCrustinA could be involved in the host defense against bacterial and viral infection. Additionally, RNAi mediated knockdown of LvCrustinA resulted in shrimps with the higher cumulative mortality during V. parahaemolyticus and WSSV infection. Taken together, these results provided some insight into the expression and transcriptional regulatory role of LvCrustinA, and its defensive role against pathogenic infection.

Identification of cold responsive genes in Pacific white shrimp (Litopenaeus vannamei) by suppression subtractive hybridization.

The Pacific white shrimp (Litopenaeus vannamei) is one of the most widely cultured shrimp species in the world. Despite L. vannamei having tropical origins, it is being reared subtropically, with low temperature stress being one of the most severe threats to its growth, survival and distribution. To unravel the molecular basis of cold tolerance in L. vannamei, the suppression subtractive hybridization (SSH) platform was employed to identify cold responsive genes in the hepatopancreas of L. vannamei. Both forward and reverse cDNA libraries were constructed, followed by dot blot hybridization, cloning, sequence analysis and quantitative real-time PCR. These approaches identified 92 cold induced and 48 cold inhibited ESTs to give a total of 37 cold induced and 17 cold inhibited contigs. Some of the identified genes related to stress response or cell defense, such as tetraspanins (TSPANs), DEAD-box helicase, heat shock proteins (HSPs) and metallothionein (MT), which were more abundant in the forward SSH library than in the reverse SSH library. The most abundant Est was a tetraspanin-8 (TSPAN8) homolog dubbed LvTSPAN8. A multiple sequence alignment and transmembrane domain prediction was also performed for LvTSPAN8. LvTSPAN8 expression was also examined in the gills, muscle, heart and hepatopancreas following cold exposure and showed the highest expression levels in the hepatopancreas. Overall, this study was able to identify several known genes and novel genes via SSH that appear to be associated with cold stress and will help to provide further insights into the molecular mechanisms regulating cold tolerance in L. vannamei.

Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei).

BACKGROUND: The Pacific white shrimp (Litopenaeus vannamei) is the world's most prevalent cultured crustacean species. However, the supply of high-quality broodstocks is limited and baseline information related to its reproductive activity and molecular issues related to gonad development are scarce. In this study, we performed transcriptome sequencing on the gonads of adult male and female L. vannamei to identify sex-related genes. RESULTS: A total of 25.16 gigabases (Gb) of sequences were generated from four L. vannamei gonadal tissue libraries. After quality control, 24.11 Gb of clean reads were selected from the gonadal libraries. De-novo assembly of all the clean reads generated a total of 65,218 unigenes with a mean size of 1021 bp and a N50 of 2000 bp. A search of all-unigene against Nr, SwissProt, KEGG, COG and NT databases resulted in 26,482, 23,062, 20,659, 11,935 and 14,626 annotations, respectively, providing a total of 30,304 annotated unigenes. Among annotated unigenes, 12,320 unigenes were assigned to gene ontology categories and 20,659 unigenes were mapped to 258 KEGG pathways. By comparing the ovary and testis libraries, 19,279 testicular up-regulated and 3,529 ovarian up-regulated unigenes were identified. Enrichment analysis of differentially expressed unigenes resulted in 1060 significantly enriched GO terms and 34 significantly enriched KEGG pathways. Nine ovary-specific, 6 testis-specific, 45 testicular up-regulated and 39 ovarian up-regulated unigenes were then confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, using all-unigenes as a reference, a total of 13,233 simple sequence repeats (SSRs) were identified in 10,411 unigene sequences. CONCLUSIONS: The present study depicts the first large-scale RNA sequencing of shrimp gonads. We have identified many important sex-related functional genes, GO terms and pathways, all of which will facilitate future research into the reproductive biology of shrimp. We expect that the SSRs detected in this study can then be used as genetic markers for germplasm evaluation of breeding and imported populations.