Search

Evolutionary developmental biology (evo-devo) has been among the most fascinating interdisciplinary fields for decades, which aims to elucidate the origin and evolution of diverse developmental processes. The rapid accumulation of omics data provides unprecedented opportunities to answer many interesting but unresolved evo-devo questions. However, the access and utilization of these resources are hindered by challenges particularly in non-model animals. Here, we establish a comparative multi-omics database for animal evo-devo (EDomics, http://edomics.qnlm.ac) containing comprehensive genomes, bulk transcriptomes, and single-cell data across 40 representative species, many of which are generally used as model organisms for animal evo-devo study. EDomics provides a systematic view of genomic/transcriptomic information from various aspects, including genome assembly statistics, gene features and families, transcription factors, transposable elements, and gene expressional profiles/networks. It also exhibits spatiotemporal gene expression profiles at a single-cell level, such as cell atlas, cell markers, and spatial-map information. Moreover, EDomics provides highly valuable, customized datasets/resources for evo-devo research, including gene family expansion/contraction, inferred core gene repertoires, macrosynteny analysis for karyotype evolution, and cell type evolution analysis. EDomics presents a comprehensive and comparative multi-omics platform for animal evo-devo community to decipher the whole history of developmental evolution across the tree of life.

The pharynx is an endoderm innovation in deuterostome ancestors, the vertebrate descendent structure of which is a pharyngeal developmental organizer involving multi-germ layer and organ derivatives. However, the evolutionary origination of complicated pharynx organs in vertebrates is still largely unknown. Endostyle, a transitional pharyngeal organ exclusively in basal chordates provides an opportunity to reveal the origin of pharyngeal organs. Here, utilizing cutting-edged Stereo-seq and single-cell RNA-seq, we constructed the first spatially-resolved single-cell atlas in the endostyle of urochordate ascidianStyela clava, where the spatial location of Stereo-seq and high capture efficiency of single-cell RNA-seq complement each other and identified 23 highly differentiated cell types. We identified a previously overlooked hemolymphoid region (HLR), which harbors immune and blood cell clusters with enriched stemness capacities, illuminating a mixed rudiment and stem-cell niches for the blood and lymphoid system. More excitingly, we discovered a mechanical-sensitive hair cell candidate in zone 3 homologous to vertebrate acoustico-lateralis system, which was supported by the expression ofin situhybridization-verified inner ear-specific markers, includingPTPRQ,USH2A,WHRN, andADGRV1, ultracellular structure evidence and cross-species comparison. These results thoroughly renewed the comprehension of the basal-chordate pharynx and provides expressional evidence for multiplexed pharyngeal organ evolution.

It has long been known that abundant symbiotic bacteria exist in the tunic and gut of marine ascidians, and that these play crucial roles in host development, physiological metabolism, and environmental adaptation. However, the identity, roles and functions of these symbiotic bacteria are known for only a few strains. In this study, we isolated and cultivated 263 strains of microorganisms from the intestine of the marine ascidian Halocynthia roretzi through a combination of aerobic and anaerobic culture approaches. Most cultivated species, both aerobic and anaerobic, from ascidian stool samples belonged to the genus Bacillus based on 16S rDNA sequencing identification and phylogenetic assays. The distribution of cultured bacteria varied with seasonal changes in environmental conditions. To explore the functions of cultured bacteria, we screened out a strain of Serratia sp. whose extracts showed high antibacterial activity against aquatic pathogens. These findings revealed the potential roles of gut microorganisms in ascidian defense and environmental adaptation, thus providing insights into the interaction and co-evolution between gut bacteria and their hosts.

Zygotic genome activation (ZGA), a universal process in early embryogenesis that occurs during the maternal-to-zygotic transition, involves reprogramming in the zygotic nucleus that initiates global transcription. In recent decades, knowledge of this process has been acquired from research on various model organisms; however, a consensus explanation of the mechanism underlying the process, especially in relation to housekeeping gene reactivation, is lacking. Here, we used hybrids derived from two ascidian species (Ciona robusta and C. savignyi), which diverged >120 Mya with significant divergence among most orthologous genes, to symmetrically document the unique dynamics of ZGA in urochordates. We found two co-ordinated waves of ZGA, representing early developmental and housekeeping gene reactivation, during the 8-cell to 110-cell stage. Comparative analysis revealed the regulatory connection between maternal and zygotic genes as well as allelic-specific expression in a species-rather than parental-related manner, which was attributed to the divergence of cis-regulatory elements. Single-cell RNA sequencing revealed that spatial differential reactivation of paternal housekeeping genes was significantly correlated with the mechanical property of each cell type. These findings potentially provide a new system for understanding the evolution and adaptation of strategies regulating ZGA in basal chordates.

Background Gut microbiota plays important roles in host animal metabolism, homeostasis and environmental adaptation. However, the interplay between the gut microbiome and urochordate ascidian, the most closet relative of vertebrate, remains less explored. In this study, we characterized the gut microbial communities of urochordate ascidian (Halocynthia roretzi) across the changes of season and starvation stress using a comprehensive set of omic approaches including 16S rRNA gene amplicon sequencing, shotgun metagenomics, metabolomic profiling, and transcriptome sequencing. Results The 16S rRNA gene amplicon profiling revealed that ascidians harbor indigenous gut microbiota distinctly different to the marine microbial community and significant variations in composition and abundance of gut bacteria, with predominant bacterial orders representing each season. Depressed alpha-diversities of gut microbiota were observed across starvation stress when compared to the communities in aquafarm condition. Synechococcales involving photosynthesis and its related biosynthesis was reduced in abundance while the enrichments of Xanthomonadales and Legionellales may facilitate bile acid biosynthesis during starvation. Metabolomics analysis found that long chain fatty acids, linolenic acid, cyanoamino acid, and pigments derived from gut bacteria were upregulated, suggesting a beneficial contribution of the gut microbiome to the ascidian under starvation stress. Conclusions Our findings revealed seasonal variation of ascidian gut microbiota. Defense and energy-associated metabolites derived from gut microbiome may provide an adaptive interplay between gut microbiome and ascidian host that maintains a beneficial metabolic system across season and starvation stress. The diversity-generating metabolisms from both microbiota and host might lead to the co-evolution and environmental adaptation. Graphical Abstract

The transient larva-bearing biphasic life cycle is the hallmark of many metazoan phyla, but how metazoan larvae originated remains a major enigma in animal evolution. There are two hypotheses for larval origin. The ‘larva-first’ hypothesis suggests that the first metazoans were similar to extant larvae, with later evolution of the adult-added biphasic life cycle; the ‘adult-first’ hypothesis suggests that the first metazoans were adult forms, with the biphasic life cycle arising later via larval intercalation. Here, we investigate the evolutionary origin of primary larvae by conducting ontogenetic transcriptome profiling for Mollusca—the largest marine phylum characterized by a trochophore larval stage and highly variable adult forms. We reveal that trochophore larvae exhibit rapid transcriptome evolution with extraordinary incorporation of novel genes (potentially contributing to adult shell evolution), and that cell signalling/communication genes (for example, caveolin and innexin) are probably crucial for larval evolution. Transcriptome age analysis of eight metazoan species reveals the wide presence of young larval transcriptomes in both trochozoans and other major metazoan lineages, therefore arguing against the prevailing larva-first hypothesis. Our findings support an adult-first evolutionary scenario with a single metazoan larval intercalation, and suggest that the first appearance of proto-larva probably occurred after the divergence of direct-developing Ctenophora from a metazoan ancestor. Evolutionary analysis of transcriptomes across Metazoa supports a scenario of ancestral adult stages and a single intercalation event at the origin of larvae.

Cell morphology and tissue integrity are essential for embryogenesis. Caveolins are membrane proteins that induce the formation of surface pits called caveolae that serve as membrane reservoirs for cell and tissue protection during development. In vertebrates, caveolin 1 (Cav1) and caveolin 3 (Cav3) are required for caveola formation. However, the formation of caveola and the function of caveolins in invertebrates are largely unknown. In this study, three caveolins, Cav‐a, Cav‐b, and CavY, are identified in the genome of the invertebrate chordate Ciona spp. Based on phylogenetic analysis, Cav‐a is found to be closely related to the vertebrate Cav1 and Cav3. In situ hybridization shows that Cav‐a is expressed in Ciona embryonic notochord and muscle. Cell‐free experiments, model cell culture systems, and in vivo experiments demonstrate that Ciona Cav‐a has the ability to induce membrane curvature at the plasma membrane. Knockdown of Cav‐a in Ciona embryos causes loss of invaginations in the plasma membrane and results in the failure of notochord elongation and lumenogenesis. Expression of a dominant‐negative Cav‐a point mutation causes cells to change shape and become displaced from the muscle and notochord to disrupt tissue integrity. Furthermore, we demonstrate that Cav‐a vesicles show polarized trafficking and localize at the luminal membrane during notochord lumenogenesis. Taken together, these results show that the invertebrate chordate caveolin from Ciona plays crucial roles in tissue integrity and morphology by inducing membrane curvature and intracellular vesicle trafficking during embryogenesis.

The solute carrier 6 (SLC6) gene family, functioning as neurotransmitter transporters, plays the crucial roles in neurotransmission, cellular and organismal homeostasis. In this study, we found an expansion of SLC6 family gene in the genome of chordate invertebrate Ciona savignyi. A total of 40 candidate genes including 29 complete and 11 putative genes were identified as SLC6 family gene homologs. Phylogenetic analysis revealed that most of these Ciona SLC6 genes were highly conserved with the vertebrate ones, although gene duplication and gene losses did exist. Four genes were selected from SLC6 subfamilies to be further investigated for their functional characteristics on cell growth and migration through overexpression approach in cultured cell lines. The results showed both SLC6A7 and SLC6A17 from amino acid transporters AA1 and AA2 sub-families, respectively, significantly suppressed the cell proliferation and migration. While SLC6A1 and SLC6A4, which were classified into GABA and monoamine transporters, respectively, did not affect the cell proliferation and migration in HEK293T, HeLa, and MCF7 cells. The whole set of C. savignyi SLC6 genes identified in this study provides an important genomic resource for future biochemical, physiological, and phylogenetic studies on SLC6 gene family. Our experimental data demonstrated that Ciona amino acid transporters, such as SLC6A7 and SLC6A17, were essential for cell physiology and behaviors, indicating their crucially potential roles in the control of cell proliferation and migration during ascidian embryogenesis.

Interferon regulatory factors (IRFs) are key transcription factors that function in the immune system via the interferon (IFN) pathway. In the current study, we identified and characterized three IRFs (CsIRFL1, CsIRFL2, and CsIRFL3) from ascidian Ciona savignyi. Phylogenetic analysis showed that CsIRFL1 was clustered with two IRFs from Ciona robusta and shrimp IRF apart from the vertebrate IRFs, whereas CsIRFL2 and CsIRFL3 were grouped with an unnamed protein from Oikopleura dioica into a sub-branch highly identifying with the vertebrate IRF4, IRF8, and IRF9. Gene expression analysis revealed that CsIRFL1 and CsIRFL2 expressed in all the examined adult tissues (stomach, intestines, eggs, hemocytes, gonad, heart, and pharynx) and predominantly in hemocytes. However, the expression of CsIRFL3 was undetectable in the tested adult tissues. Furthermore, in situ hybridization showed that CsIRFL1 and CsIRFL2 mainly expressed in immunocytes within hemolymph, including phagocytes, macrophage-like cells, morula cells, and amoebocytes, suggesting CsIRFL1 and CsIRFL2 were involved in ascidian immune responses. We then performed LPS and poly(I:C) challenge assay and found that CsIRFL1 highly expressed in the cultured hemocytes following LPS infection for 24 h. After viral analogue poly(I:C) stimulation, the expression of CsIRFL2 was dramatically upregulated from 12 to 24 h. Meanwhile, two critical components of the IFN signaling pathways, STAT and TBK1, showed the increased expression as well after poly(I:C) induction, indicating that CsIRFL2 and IFN pathways genes were activated under the infection of viral analogue. Thus, our findings suggested that CsIRFL2 was a potential transcriptional regulatory factor that participated in regulating the ascidian anti-virus immune response.

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.

Penaeid shrimps are among the most economically important crustaceans, which provide an important global food source. These species exhibit complex body plans and novelties, such as segments and appendages, which render them interesting organisms for developmental biology study of crustaceans. However, limited genomic resources have been put forward for the researches of them. Here, we report the genome sequencing and draft assembly of two economically important penaeid shrimp species, Marsupenaeus japonicus and Penaeus monodon. A total of 132.86 Gb and 132.83 Gb sequencing data was obtained in the two shrimp species. The genome assembly, a total length of 1.94 Gb and 2.04 Gb in M. japonicus and P. monodon, respectively, covers more than 97% of coding regions. We further identified 626 Mb (34.96%) and 833 Mb (46.68%) repeats, 16,716 and 18,100 genes in these two genomes, respectively. We also identified Hox genes that are important to their body plans. These data will provide valuable resources for the study of selective breeding and some plastic biological characters of penaeid shrimps, including molting, lobstering, brooding eggs and sensitization in humans.

The extracellular matrix (ECM) not only provides essential physical scaffolding for cellular constituents but also initiates crucial biochemical and biomechanical cues that are required for tissue morphogenesis. In this study, we utilized wheat germ agglutinin (WGA) staining to characterize the ECM architecture in ascidian embryos and larvae. The results showed three distinct populations of ECM presenting in Ciona embryogenesis: the outer layer localized at the surface of embryo, an inner layer of notochord sheath and the apical ECM secreted by the notochord. To further elucidate the precise structure of Ciona embryonic ECM, we employed scanning and transmission electron microscopy, and found that the outer membrane was relatively thick with short fibres, whereas the ECM layer in notochord sheath was not as thick as the outer membrane but more regular arranged; the lumen between notochord cells was hydrostatic and sticky. Then, we used the RNA sequencing data from the embryos and larvae of Ciona savignyi to identify ECM genes and acquire their expression patterns. We identified 115 unigenes as 67 ECM genes, and 77 unigenes showed dynamic expression changes between different stages. Our results reveal the architecture, molecular composition and dynamic expression profile of ECM in ascidian embryogenesis, and may increase understanding of the function of the ECM in chordate development., Summary: This study reveals the architecture, molecular composition and dynamic expression profile of the extracellular matrix in ascidian embryos and larvae, providing clues for its function in chordate development.

Tunicates occupy the evolutionary position at the boundary of invertebrates and vertebrates. It exhibits adaptation to broad environmental conditions and is distributed globally. Despite hundreds of years of embryogenesis studies, the genetic basis of the invasive habits of ascidians remains largely unknown. The leathery sea squirt, Styela clava, is an important invasive species. We used the chromosomal‐level genome and transcriptome of S. clava to explore its genomic‐ and molecular‐network‐based mechanisms of adaptation to environments. Compared with Ciona intestinalis type A (C. robusta), the size of the S. clava genome was expanded by 2‐fold, although the gene number was comparable. An increase in transposon number and variation in dominant types were identified as potential expansion mechanisms. In the S. clava genome, the number of genes encoding the heat‐shock protein 70 family and members of the complement system was expanded significantly, and cold‐shock protein genes were transferred horizontally into the S. clava genome from bacteria. The expanded gene families potentially play roles in the adaptation of S. clava to its environments. The loss of key genes in the galactan synthesis pathway might explain the distinct tunic structure and hardness compared with the ascidian Ciona species. We demonstrated further that the integrated thyroid hormone pathway participated in the regulation of larval metamorphosis that provides S. clava with two opportunities for adapting to their environment. Thus, our report of the chromosomal‐level leathery sea squirt genome provides a comprehensive genomic basis for the understanding of environmental adaptation in tunicates.

In penaeid shrimp, mesoderm forms from two sources: naupliar mesoderm founder cells, which invaginate during gastrulation, and posterior mesodermal stem cells called mesoteloblasts, which undergo characteristic teloblastic divisions. The primordial mesoteloblast descends from the ventral mesendoblast, which arrests in cell division at the 32-cell stage and ingresses with its sister dorsal mesendoblast prior to naupliar mesoderm invagination. The naupliar mesoderm forms the muscles of the naupliar appendages (first and second antennae and mandibles), while the mesoteloblasts form the mesoderm, including the muscles, of subsequently formed posterior segments. To better understand the mechanism of mesoderm and muscle formation in penaeid shrimp, twist, snail, and mef2 cDNAs were identified from transcriptomes of Penaeus vannamei, P. japonicus, P. chinensis, and P. monodon. A single Twist ortholog was found, with strong inferred amino acid conservation across all three species. Multiple Snail protein variants were detected, which clustered in a phylogenetic tree with other decapod crustacean Snail sequences. Two closely-related mef2 variants were found in P. vannamei. The developmental mRNA expression of these genes was studied by qPCR in P. vannamei embryos, larvae, and postlarvae. Expression of Pv-twist and Pv-snail began during the limb bud stage and continued through larval stages to the postlarva. Surprisingly, Pv-mef2 expression was found in all stages from the zygote to the postlarva, with the highest expression in the limb bud and protozoeal stages. The results add comparative data on the development of anterior and posterior mesoderm in malacostracan crustaceans, and should stimulate further studies on mesoderm and muscle development in penaeid shrimp.

Background Metamorphosis takes place within the life cycle of most marine invertebrates. The marine ascidian is a classical model to study complex cellular processes and underlying molecular mechanisms involved in its larval metamorphosis. The detailed molecular signaling pathways remain elusive, though extracellular signal-regulated kinases (ERKs) and c-Jun N-terminal kinase (JNK) have been revealed to regulate cell migration, differentiation, and apoptosis in ascidian larval organ regression and juvenile organ development. MicroRNAs (miRNAs) are small non-coding RNAs that modulate gene expression at the post-transcriptional level. Large numbers of miRNAs have been demonstrated to be involved in many developmental and metamorphic processes. However, the identification of miRNAs in ascidian larval metamorphosis has not yet been investigated. Results Totally, 106 known and 59 novel miRNAs were screened out through RNA-sequencing of three small RNA libraries from 18 to 21-h post-fertilization (hpf) tailbud embryos as well as from 42 hpf larvae (after tail regression) in Ciona savignyi. Expression profiling of miRNAs was confirmed by quantitative real-time PCR, showing that the expression levels of csa-miR-4040, csa-miR-4086, csa-miR-4055, csa-miR-4060, csa-miR-216a, csa-miR-216b, csa-miR-217, csa-miR-183, and csa-miR-92c were significantly higher in 42 hpf larvae, whereas those of csa-miR-4018a, csa-miR-4018b, and csa-miR-4000f were higher in 18 and 21 hpf embryos; then, their expression in 42 hpf larvae became significantly low. For these 12 miRNAs, whose expression levels significantly changed, we predicted their target genes through the combination of miRanda and TargetScan. This prediction analysis revealed 332 miRNA-target gene pairs that were associated with the ERK, JNK, and transforming growth factor beta signaling pathways, suggesting that the identified miRNAs are involved in the regulation of C. savignyi larval metamorphosis via controlling the expression of their target genes. Furthermore, we validated the expression of five selected miRNAs by northern blotting. Among the selected miRNAs, the expression patterns of csa-miR-4018a, csa-miR-4018b, and csa-miR-4000f were further examined by whole-mount in situ hybridization. The results showed that all three miRNAs were specifically expressed in a cell population resembling mesenchymal cells at the head and trunk part in swimming larvae but not in metamorphic larvae. Utilizing the luciferase assay, we also confirmed that miR-4000f targeted Mapk1, suggesting that the csa-miR-4018a/csa-miR-4018b/csa-miR-4000f cluster regulates larval metamorphosis through the Mapk1-mediated signaling pathway. Conclusions Totally, 165 miRNAs, including 59 novel ones, were identified from the embryos and larvae of C. savignyi. Twelve of them showed significant changes in expression before and during metamorphosis. In situ hybridization and northern blotting results revealed that three miRNAs are potentially involved in the signaling regulatory network for the migration and differentiation of mesenchymal cells in larval metamorphosis. Furthermore, the luciferase reporter assay revealed that Mapk1 is a target of csa-miR-4000f. Our results not only present a list and profile of miRNAs involved in Ciona metamorphosis but also provide informative cues to further understand their function in ascidian larval metamorphosis. Electronic supplementary material The online version of this article (10.1186/s12864-018-4566-4) contains supplementary material, which is available to authorized users.

Figure S1. Original images of northern blotting presented in Fig. 5. (PDF 1785Â kb)

Table S2. Sequences of identified novel miRNAs in C. savignyi. (DOCX 103Â kb)

Table S1. Sequences of identified known miRNAs in C. savignyi. (DOCX 108Â kb)

Table S4. Sequence of primers for mutant constructs used in overlap extension PCR. (DOCX 60Â kb)

Table S3. Sequences of probes used in northern blotting and in situ hybridization. (DOCX 47Â kb)

Figure S2. Validation of the interaction between csa-miR-4018a or csa-miR-4018b and Mapkk3â 3â ˛-UTR. (PDF 841Â kb)

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.

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.

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.

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.

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.

Peroxiredoxins (Prx) are a family of antioxidant proteins and perform important functions in intracellular signal transduction. Here, we report a Prx gene from Chinese shrimp Fenneropenaeus chinensis. The full-length cDNA of FcPrx gene contained an open reading frame of 735 bp encoding a polypeptide of 275 amino acids. The molecular mass of the deduced amino acid of FcPrx is 27445.43 Da with an estimated pI of 5.71. Sequence comparison showed that the FcPrx shares high identities with Prx IVs and it was named FcPrx4. A real-Time PCR (qRT-PCR) assay was developed to assess the mRNA expression of FcPrx4 in different tissues and temporal expression in hemocytes and hepatopancreas of F. chinensis challenged by white spot syndrome virus (WSSV). Transcripts of FcPrx4 can be detected in all tissues examined. The expression of FcPrx4 showed significant up-regulation in shrimp hemocytes and hepatopancreas after artificial infection with WSSV. A fusion protein containing FcPrx4 was produced in vitro and was confirmed by Matrix-assisted laser desorption ionization mass spectrometry (MALDI-MS) assay. And activity analysis indicated that the recombinant FcPrx4 proteins can reduce H2O2 in the presence of dithiothreitol.

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.

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.