Jiankai Wei

SNP discovery in the transcriptome of white Pacific shrimp Litopenaeus vannamei by next generation sequencing.

The application of next generation sequencing technology has greatly facilitated high throughput single nucleotide polymorphism (SNP) discovery and genotyping in genetic research. In the present study, SNPs were discovered based on two transcriptomes of Litopenaeus vannamei (L. vannamei) generated from Illumina sequencing platform HiSeq 2000. One transcriptome of L. vannamei was obtained through sequencing on the RNA from larvae at mysis stage and its reference sequence was de novo assembled. The data from another transcriptome were downloaded from NCBI and the reads of the two transcriptomes were mapped separately to the assembled reference by BWA. SNP calling was performed using SAMtools. A total of 58,717 and 36,277 SNPs with high quality were predicted from the two transcriptomes, respectively. SNP calling was also performed using the reads of two transcriptomes together, and a total of 96,040 SNPs with high quality were predicted. Among these 96,040 SNPs, 5,242 and 29,129 were predicted as non-synonymous and synonymous SNPs respectively. Characterization analysis of the predicted SNPs in L. vannamei showed that the estimated SNP frequency was 0.21% (one SNP per 476 bp) and the estimated ratio for transition to transversion was 2.0. Fifty SNPs were randomly selected for validation by Sanger sequencing after PCR amplification and 76% of SNPs were confirmed, which indicated that the SNPs predicted in this study were reliable. These SNPs will be very useful for genetic study in L. vannamei, especially for the high density linkage map construction and genome-wide association studies.

Comparative Transcriptomic Characterization of the Early Development in Pacific White Shrimp Litopenaeus vannamei

Penaeid shrimp has a distinctive metamorphosis stage during early development. Although morphological and biochemical studies about this ontogeny have been developed for decades, researches on gene expression level are still scarce. In this study, we have investigated the transcriptomes of five continuous developmental stages in Pacific white shrimp (Litopenaeus vannamei) with high throughput Illumina sequencing technology. The reads were assembled and clustered into 66,815 unigenes, of which 32,398 have putative homologues in nr database, 14,981 have been classified into diverse functional categories by Gene Ontology (GO) annotation and 26,257 have been associated with 255 pathways by KEGG pathway mapping. Meanwhile, the differentially expressed genes (DEGs) between adjacent developmental stages were identified and gene expression patterns were clustered. By GO term enrichment analysis, KEGG pathway enrichment analysis and functional gene profiling, the physiological changes during shrimp metamorphosis could be better understood, especially histogenesis, diet transition, muscle development and exoskeleton reconstruction. In conclusion, this is the first study that characterized the integrated transcriptomic profiles during early development of penaeid shrimp, and these findings will serve as significant references for shrimp developmental biology and aquaculture research.

Whole Transcriptome Analysis Provides Insights into Molecular Mechanisms for Molting in Litopenaeus vannamei

Molting is one of the most important biological processes in shrimp growth and development. All shrimp undergo cyclic molting periodically to shed and replace their exoskeletons. This process is essential for growth, metamorphosis, and reproduction in shrimp. However, the molecular mechanisms underlying shrimp molting remain poorly understood. In this study, we investigated global expression changes in the transcriptomes of the Pacific white shrimp, Litopenaeus vannamei, the most commonly cultured shrimp species worldwide. The transcriptome of whole L. vannamei was investigated by RNA-sequencing (RNA-seq) throughout the molting cycle, including the inter-molt (C), pre-molt (D0, D1, D2, D3, D4), and post-molt (P1 and P2) stages, and 93,756 unigenes were identified. Among these genes, we identified 5,117 genes differentially expressed (log2ratio ≥1 and FDR ≤0.001) in adjacent molt stages. The results were compared against the National Center for Biotechnology Information (NCBI) non-redundant protein/nucleotide sequence database, Swiss-Prot, PFAM database, the Gene Ontology database, and the Kyoto Encyclopedia of Genes and Genomes database in order to annotate gene descriptions, associate them with gene ontology terms, and assign them to pathways. The expression patterns for genes involved in several molecular events critical for molting, such as hormone regulation, triggering events, implementation phases, skelemin, immune responses were characterized and considered as mechanisms underlying molting in L. vannamei. Comparisons with transcriptomic analyses in other arthropods were also performed. The characterization of major transcriptional changes in genes involved in the molting cycle provides candidates for future investigation of the molecular mechanisms. The data generated in this study will serve as an important transcriptomic resource for the shrimp research community to facilitate gene and genome annotation and to characterize key molecular processes underlying shrimp development.

Horizontally transferred genes in the genome of Pacific white shrimp, Litopenaeus vannamei.

BackgroundIn recent years, as the development of next-generation sequencing technology, a growing number of genes have been reported as being horizontally transferred from prokaryotes to eukaryotes, most of them involving arthropods. As a member of the phylum Arthropoda, the Pacific white shrimp Litopenaeus vannamei has to adapt to the complex water environments with various symbiotic or parasitic microorganisms, which provide a platform for horizontal gene transfer (HGT).ResultsIn this study, we analyzed the genome-wide HGT events in L. vannamei. Through homology search and phylogenetic analysis, followed by experimental PCR confirmation, 14 genes with HGT event were identified: 12 of them were transferred from bacteria and two from fungi. Structure analysis of these genes showed that the introns of the two fungi-originated genes were substituted by shrimp DNA fragment, two genes transferred from bacteria had shrimp specific introns inserted in them. Furthermore, around other three bacteria-originated genes, there were three large DNA segments inserted into the shrimp genome. One segment was a transposon that fully transferred, and the other two segments contained only coding regions of bacteria. Functional prediction of these 14 genes showed that 6 of them might be related to energy metabolism, and 4 others related to defense of the organism.ConclusionsHGT events from bacteria or fungi were happened in the genome of L. vannamei, and these horizontally transferred genes can be transcribed in shrimp. This is the first time to report the existence of horizontally transferred genes in shrimp. Importantly, most of these genes are exposed to a negative selection pressure and appeared to be functional.

RNA-Seq reveals the dynamic and diverse features of digestive enzymes during early development of Pacific white shrimp Litopenaeus vannamei

The Pacific white shrimp (Litopenaeus vannamei), with high commercial value, has a typical metamorphosis pattern by going through embryo, nauplius, zoea, mysis and postlarvae during early development. Its diets change continually in this period, and a high mortality of larvae also occurs in this period. Since there is a close relationship between diets and digestive enzymes, a comprehensive investigation about the types and expression patterns of all digestive enzyme genes during early development of L. vannamei is of considerable significance for shrimp diets and larvae culture. Using RNA-Seq data, the types and expression characteristics of the digestive enzyme genes were analyzed during five different development stages (embryo, nauplius, zoea, mysis and postlarvae) in L. vannamei. Among the obtained 66,815 unigenes, 296 were annotated as 16 different digestive enzymes including five types of carbohydrase, seven types of peptidase and four types of lipase. Such a diverse suite of enzymes illustrated the capacity of L. vannamei to exploit varied diets to fit their nutritional requirements. The analysis of their dynamic expression patterns during development also indicated the importance of transcriptional regulation to adapt to the diet transition. Our study revealed the diverse and dynamic features of digestive enzymes during early development of L. vannamei. These results would provide support to better understand the physiological changes during diet transition.

Genome Sequences of Marine Shrimp Exopalaemon carinicauda Holthuis Provide Insights into Genome Size Evolution of Caridea

Crustacea, particularly Decapoda, contains many economically important species, such as shrimps and crabs. Crustaceans exhibit enormous (nearly 500-fold) variability in genome size. However, limited genome resources are available for investigating these species. Exopalaemon carinicauda Holthuis, an economical caridean shrimp, is a potential ideal experimental animal for research on crustaceans. In this study, we performed low-coverage sequencing and de novo assembly of the E. carinicauda genome. The assembly covers more than 95% of coding regions. E. carinicauda possesses a large complex genome (5.73 Gb), with size twice higher than those of many decapod shrimps. As such, comparative genomic analyses were implied to investigate factors affecting genome size evolution of decapods. However, clues associated with genome duplication were not identified, and few horizontally transferred sequences were detected. Ultimately, the burst of transposable elements, especially retrotransposons, was determined as the major factor influencing genome expansion. A total of 2 Gb repeats were identified, and RTE-BovB, Jockey, Gypsy, and DIRS were the four major retrotransposons that significantly expanded. Both recent (Jockey and Gypsy) and ancestral (DIRS) originated retrotransposons responsible for the genome evolution. The E. carinicauda genome also exhibited potential for the genomic and experimental research of shrimps.

Transcriptome analysis on the exoskeleton formation in early developmetal stages and reconstruction scenario in growth-moulting in Litopenaeus vannamei

Exoskeleton construction is an important issue in shrimp. To better understand the molecular mechanism of exoskeleton formation, development and reconstruction, the transcriptome of the entire developmental process in Litopenaeus vannamei, including nine early developmental stages and eight adult-moulting stages, was sequenced and analysed using Illumina RNA-seq technology. A total of 117,539 unigenes were obtained, with 41.2% unigenes predicting the full-length coding sequence. Gene Ontology, Clusters of Orthologous Group (COG), the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and functional annotation of all unigenes gave a better understanding of the exoskeleton developmental process in L. vannamei. As a result, more than six hundred unigenes related to exoskeleton development were identified both in the early developmental stages and adult-moulting. A cascade of sequential expression events of exoskeleton-related genes were summarized, including exoskeleton formation, regulation, synthesis, degradation, mineral absorption/reabsorption, calcification and hardening. This new insight on major transcriptional events provide a deep understanding for exoskeleton formation and reconstruction in L. vannamei. In conclusion, this is the first study that characterized the integrated transcriptomic profiles cover the entire exoskeleton development from zygote to adult-moulting in a crustacean, and these findings will serve as significant references for exoskeleton developmental biology and aquaculture research.

Identification and expression analysis of long noncoding RNAs in embryogenesis and larval metamorphosis of Ciona savignyi

Long noncoding RNAs (lncRNAs) play important roles in diverse developmental and pathological processes through chromatin reprogramming, cis regulation and posttranscriptional modification. They have been extensively studied in both vertebrates and invertebrates. However, the information of lncRNAs in urochordate is still lacking. In this study, we used the RNA-Seq data from three developmental stages (18, 21 and 42 hours post fertilization, hpf) of embryos and larvae in Ciona savignyi to identify candidate lncRNAs and analyze their expression profiles. A total of 29,944 unigenes were predicted as lncRNAs, five of which had hits with lncRNAs in NONCODE database. The acquired lncRNAs had an average length of 466 nt. The peaks of length, GC content and minimum free energy of the lncRNAs were significantly lower than that of the message RNAs (mRNAs). The average expression levels of lncRNAs were also lower than those of mRNAs. Among the three developmental stages, highly expressed lncRNAs concentrated in 18 hpf embryos. While, for those lncRNAs specifically up-regulated in 21 hpf embryos, their co-expressed mRNAs were enriched in GO terms of membrane, indicating these lncRNAs are involved in the regulation of luminal membrane biogenesis, and extracellular matrix secretion through membrane localized proteins during Ciona notochord tubulogenesis. The lncRNAs in 42 hpf larvae were distinct from those in 18 and 21 hpf embryos. This result is associated with the fact that swimming larvae are transiting into metamorphic juveniles at this stage, indicating lncRNAs are involved in the regulation of larval metamorphosis. Overall, our study identified a large number of lncRNAs in C. savignyi and revealed their expression characteristics and dynamics during Ciona embryogenesis and larval metamorphosis. The results will help to further understand the function of lncRNAs in chordate development and the evolution of lncRNAs.

Penaeid shrimp brachyury: sequence analysis and expression during gastrulation

Gastrulation occurs by a variety of morphogenetic movements, often correlated with diverse expression of the T-box transcription factor Brachyury (Bra). Bra may be expressed in ectoderm, mesoderm, or endoderm, but its role in cell fate specification or regulation of gastrulation movements has not been studied in the development of crustaceans. Penaeid shrimp (Decapoda: Dendrobranchiata: Penaeidae) develop by complete cleavage and gastrulation by invagination to a free-swimming nauplius larva. Penaeid gastrulation diverges from other decapods and from insects, occurring early at a low cell number with the formation of a radial invagination. Toward a better understanding of gastrulation movements in penaeid shrimp, bra was identified from newly available penaeid shrimp genomes and transcriptomes of Litopenaeus vannamei, Marsupenaeus japonicus, and Penaeus monodon. Additional bra homologs were obtained from the outgroups Sicyonia ingentis (Decapoda: Dendrobranchiata: Sicyoniidae) and the caridean shrimp Caridina multidentata (Decapoda: Pleocymata). The genes encoded penaeid shrimp Bra proteins of 551–552 amino acids, containing the highly conserved T-box DNA-binding region. The N-terminal Smad1-binding domain, conserved in most animals, was absent in shrimp Bra. The R1 repressor domain was the best conserved of the C-terminal regulatory domains, which were widely divergent compared to other species. The penaeid shrimp bra gene consisted of six exons, with splice sites conserved with other phyla across the animal kingdom. Real-time qPCR and FPKM analysis showed that shrimp bra mRNA was strongly expressed during gastrulation. These findings begin to address the evolution of gastrulation in shrimp at the molecular level.

Differential gene expression analysis based on expressed sequence tags(EST) from different tissues of Fenneropenaeus chinensis

Bioinformatics analysis was conducted based on the acquired expressed sequence tags(ESTs)data from different tissues of Fenneropenaeus chinensis.A total of 10 446,2 690,1 067 and 1 282 original ESTs were obtained from cephalothorax,blood,eyestalk,and ovary of adult F.chinensis respectively.By clustering and assembling,3 454,1 053,406 and 544 unigenes were generated,and then annotated by searching in NR,GO,KEGG databases.The tissue specific transcripts were identified through sequence homology analysis and classified by GO annotation.The result indicated that specific transcripts in blood group were specifically enriched in the GO term of virion part,locomotion,rhythmic process,cell killing and multi-organism process.Specific transcripts in eyestalk group were specifically enriched in the GO term of molecular transducer activity,response to stimulus,multicellular organismal process,pigmentation and biological regulation.Specific transcripts in ovary group were specifically enriched in the GO term of nutrient reservoir activity,reproduction,anatomical structure formation,transporter activity and establishment of localization.We also analyzed the highly expressed genes in each tissue according to the number of ESTs of each unigene.The result indicated that genes encoding peritrophin,elongation factor 1-alpha,thrombospondin and arginine kinase were widespread and highly expressed in tissues,suggesting that they were involved in a variety of important biological processes in shrimp.In addition,genes encoding penaeidin,cytochrome c oxidase and 14-3-3-like protein were highly expressed in blood group and genes encoding arrestin and rhodopsin were highly expressed in eyestalk group.For understanding the expressions of common genes,we analyzed genes related to peritrophin and peroxiredoxin,revealing that genes of peritrophin in F.chinensis were multiform and highly expressed.This result may provide reference for further functional study of peritrophin.