Jianbo Yuan

Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific White Shrimp Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is the dominant crustacean species in global seafood mariculture. Understanding the genome and genetic architecture is useful for deciphering complex traits and accelerating the breeding program in shrimp. In this study, a genome survey was conducted and a high-density linkage map was constructed using a next-generation sequencing approach. The genome survey was used to identify preliminary genome characteristics and to generate a rough reference for linkage map construction. De novo SNP discovery resulted in 25,140 polymorphic markers. A total of 6,359 high-quality markers were selected for linkage map construction based on marker coverage among individuals and read depths. For the linkage map, a total of 6,146 markers spanning 4,271.43 cM were mapped to 44 sex-averaged linkage groups, with an average marker distance of 0.7 cM. An integration analysis linked 5,885 genome scaffolds and 1,504 BAC clones to the linkage map. Based on the high-density linkage map, several QTLs for body weight and body length were detected. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of L. vannamei and other penaeid shrimp species.

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.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Virus-derived small RNAs in the penaeid shrimp Fenneropenaeus chinensis during acute infection of the DNA virus WSSV

Small interfering RNAs (siRNAs) and microRNAs (miRNAs) are two classes of small RNAs (sRNAs) that are critical for virus-host interplay via the RNA interference (RNAi) pathway. One virus-derived siRNA and numerous miRNAs has been reported for the double-stranded DNA virus white spot syndrome virus (WSSV), however, the expression profiles of these different types of sRNAs have not been assessed. Here, by sequencing the sRNAs and mRNAs of WSSV-infected Chinese shrimp (Fenneropenaeus chinensis), we found that the viral transcripts were universally targeted by WSSV-derived siRNAs, supporting a pivotal role for RNAi in the anti-viral immunity of shrimp. The genesis of WSSV-derived siRNAs was associated with long RNA structures. Moreover, by separating miRNAs from siRNAs, 12 WSSV miRNAs were identified. Investigation of conserved viral miRNA targets in different host species indicated the involvement of viral miRNAs in host immune responses. Collectively, our data provide new insights into the role of the RNAi pathway in the interplay between DNA viruses and crustaceans.

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.

Wnt gene family members and their expression profiling in Litopenaeus vannamei

ABSTRACT The Wnt gene family encodes secreted glycoproteins involved in a wide variety of biological processes, including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt pathway exists in all metazoan animals, however, the relevant research is rare in crustaceans. Here we described 12 Wnt genes representing 12 Wnt gene subfamilies in the Pacific white shrimp, Litopenaeus vannamei. Based on homolog annotations and phylogenetic analyses, we named these 12 Wnt genes as LvWnt1, LvWnt2, LvWnt4‐11, LvWnt16, and LvWntA. All the corresponding LvWnt proteins shared a conserved Wnt1 domain and 22 conserved cysteine residues. LvWnt1 and LvWnt6 were adjacent in a scaffold in the shrimp genome. Furthermore, we performed expression analyses of LvWnt genes at different developmental stages, during the molting process, in different tissues and after different pathogenic infection. We showed that each LvWnt gene had a unique expression pattern at different developmental stages but only a few of them expressed in adult shrimp. All the investigated LvWnt genes were initially expressed at the gastrula or limb bud embryo stages. Among them, LvWnt8 was specifically high expressed only in early embryos. LvWntA and LvWnt5 displayed high and similar expression profiles during the molting process, and LvWnt6 and LvWnt16 were specifically expressed in the thoracic ganglion, ventral nerve, intestines and gill tissues, respectively. We also found the expression of LvWntA, LvWnt5, LvWnt6, LvWnt9, and LvWnt16 were varied in the different tissues after infected with Staphylococcus aureus, Vibrio parahaemolyticus and white spot syndrome virus (WSSV), which indicated that they might participate in immune response in L. vannamei. This study provided an insight into the repertoire of the Wnt gene structure and expression in shrimps, and furthermore, might promote the understanding of development, growth and immune response of shrimps and crustaceans. HighlightsWe first identified and analyzed 12 Wnt gens in Litopenaeus vannamei.LvWnt genes play important role in development, molting and immunity of L. vannamei.LvWnt5 was up‐regulated when treated with S. aureus, V. parahaemolyticus and WSSV, and LvWnt9 was down‐regulated.LvWnt genes may regulate immune response when treated with different pathogens.

Analysis of genomic characters reveals that four distinct gene clusters are correlated with different functions in Burkholderia cenocepacia AU 1054

Possessing three circular chromosomes is a distinct genomic characteristic of Burkholderia cenocepacia AU 1054, a clinically important pathogen in cystic fibrosis. In this study, base composition, codon usage and functional role category were analyzed in the B. cenocepacia AU 1054 genome. Although no bias in the base and codon usage was detected between any two chromosomes, function differences did exist in the genes of each chromosome. Similar base composition and differential functional role categories indicated that genes on these three chromosomes were relatively stable and that a proper division of labor was established. Based on variations in the base or codon usage, four small gene clusters were observed in all of the genes. Multivariate analysis revealed that protein hydrophobicity played a predominant role in shaping base usage bias, while horizontal gene transfer and the gene expression level were the two most important factors that affected the codon usage bias. Interestingly, we also found that these gene clusters were correlated with different biological functions: (i) 45 pyrimidine-leading-codon preferred genes were predominantly involved in regulatory function; (ii) most drug resistance-related genes involved in 826 genes that coding for hydrophobic proteins; (iii) most of the 111 horizontal transfer genes were responsible for genomic plasticity; and (iv) 73 highly expressed genes (predicted by their codon adaptation index values) showed environmental adaptation to cystic fibrosis. Our results showed that genes with base or codon usage bias were affected by mutational pressure and natural selection, and their functions could contribute to drug assistance and transmissible activity in B. cenocepacia.

Actin genes and their expression in pacific white shrimp, Litopenaeus vannamei

Actin is a multi-functional gene family that can be divided into muscle-type actins and non-muscle-type actins. In this study, 37 unigenes encoding actins were identified from RNA-Seq data of Pacific white shrimp, Litopenaeus vannamei. According to phylogenetic analysis, four and three cDNAs belong to cytoplasmic- and heart-type actins and were named LvActinCT and LvActinHT, respectively. 10 cDNAs belong to the slow-type skeletal muscle actins, and 18 belong to the fast-type skeletal muscle actins; they were designated LvActinSSK and LvActinFSK, respectively. Some muscle actin genes formed gene clusters in the genome. Multiple alternative transcription starts sites (ATSSs) were found for LvActinCT1. Based on the early developmental expression profile, almost all LvActins were highly expressed between the early limb bud and post-larval stages. Using LvActinSSK5 as probes, slow-type muscle was localized in pleopod muscle and superficial ventral muscle. We also found three actin genes that were down-regulated in the hemocytes of white spot syndrome virus (WSSV)- and Vibrio parahaemolyticus-infected L. vannamei. This study provides valuable information on the actin gene structure of shrimp, furthers our understanding of the shrimp muscle system and helps us develop strategies for disease control and sustainable shrimp farming.