De Novo assembly and annotation of the freshwater crayfish Astacus astacus transcriptome

A novel small non-coding RNA 562 mediates the virulence of Pseudomonas plecoglossicida by regulating the expression of fliP, a key component of flagella T3SS

Pseudomonas plecoglossicida is a vital pathogen that poses a substantial risk to aquaculture. Small RNAs (sRNAs) are non-coding regulatory molecules capable of sensing environmental changes and modulating virulence-associated signaling pathways, such as the assembly of flagella. However, the relevant researches on P. plecoglossicida are an urgent need. Here, we report a novel sRNA, sRNA562, which has potential to regulate the post-transcriptional of fliP, a key component of the lateral flagellar type III secretion system. In this study, the effects of sRNA562 on the virulence of P. plecoglossicida and its role in regulating the pathogenic process were investigated through the use of a constructed sRNA562 deletion strain. The deletion of sRNA562 resulted in an up-regulation of fliP in P. plecoglossicida, and leading to increased swarming motility and enhanced the ability of biofilm formation, adhesion and chemotaxis. Subsequent artificial infection experiment demonstrated that the deletion of sRNA562 increased the virulence of P. plecoglossicida towards hybrid grouper, as evidenced by a reduction in survival rate, elevation of tissue bacterial load, and the exacerbation of histopathological damage. Further studies have found that the deletion of sRNA562 lead to an up-regulation of fliP expression during hybrid grouper infection, thereby enhancing bacterial swarming ability and ultimately heightening pathogenicity, leading to a dysregulated host response to infection, tissue damage and eventually death. Our work revealed a sRNA that exerts negative regulation on the expression of lateral flagella in P. plecoglossicida, thereby impacting its virulence. These findings provide a new perspective on the virulence regulation mechanism of P. plecoglossicida, contributing to a more comprehensive understanding in the field of pathogenicity research.

Protocols for studying the crayfish plague pathogen, Aphanomyces astaci, and its host-pathogen interactions

The crayfish plague caused by the pathogen Aphanomyces astaci has decimated the European and Asian populations of freshwater crayfish and represents an important threat to the other highly susceptible crayfish species in the world, such as the Australian, Madagascar, and South American species. The development and application of molecular methods addressed to the identification of A. astaci has increased exponentially during the last decades in contrast to a slow trend of the pathogen biology and host interaction. There is still a need for a better comprehension of the A. astaci-crayfish interactions, specifically the resistance and tolerance immune mechanism. These types of studies required a robust basic knowledge on the developmental biology of the pathogen in order to reproduce life stages and to perform infection experiments. A great piece of work in this area was carried out during the 1960 s to 80 s in University of Uppsala. Thus, the purpose of this work was to update previous protocols as well as to generate new guidelines to reproduce key developmental biology stages of A. astaci, to eventually identify crayfish populations with higher resistance and tolerance to this pathogen. This work also refers to other methodologies and guidelines for the diagnosis of crayfish plague, the pathogen isolation, and the in vitro production of zoospores.

Transcriptomic analysis reveals potential candidate pathways and genes involved in toxin biosynthesis in true toads

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across ten species of North American true toads (Bufonidae). Upon assembly and annotation of the ten transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g. steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the ten focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

De novo Assembly, Annotation, and Analysis of Transcriptome Data of the Ladakh Ground Skink Provide Genetic Information on High-Altitude Adaptation

The Himalayan Arc is recognized as a global biodiversity hotspot. Among its numerous cryptic and undiscovered organisms, this composite high-mountain ecosystem harbors many taxa with adaptations to life in high elevations. However, evolutionary patterns and genomic features have been relatively rarely studied in Himalayan vertebrates. Here, we provide the first well-annotated transcriptome of a Greater Himalayan reptile species, the Ladakh Ground skink Asymblepharus ladacensis (Squamata: Scincidae). Based on tissues from the brain, an embryonic disc, and pooled organ material, using pair-end Illumina NextSeq 500 RNAseq, we assembled ~77,000 transcripts, which were annotated using seven functional databases. We tested ~1600 genes, known to be under positive selection in anurans and reptiles adapted to high elevations, and potentially detected positive selection for 114 of these genes in Asymblepharus. Even though the strength of these results is limited due to the single-animal approach, our transcriptome resource may be valuable data for further studies on squamate reptile evolution in the Himalayas as a hotspot of biodiversity.

Comprehensive phylogeny of Konosirus punctatus (Clupeiformes: Clupeidae) based on transcriptomic data

Konosirus punctatus is an important species for the structure of marine ecosystems. Meanwhile, it is a native species in the northwest Pacific Ocean and supports important commercial fishery. In the present study, we generated the whole transcriptome of K. punctatus from combined tissues (muscle, liver, gill, heart, kidney, swim bladder and sexual gland) using Illumina RNA-seq technology and a total of 46087110 clean reads were obtained, corresponding to 6531521430 nucleotides. Meanwhile, 10000 clean reads were randomly selected and compared with NT database to examine the possible contamination. Results showed that 6754 clean reads were distributed among some species closely related with K. punctatus, indicating no-pollution. De novo assembly was performed and all clean reads were assembled to produce 71610 longest unigenes with an N50 of 906 bp. Among all the unigenes, 43974 unigenes were annotated in at least one database and 3172 unigenes were annotated in all databases. All unigenes were further analyzed to predict the gene structure and we have obtained a total of 54864 coding sequences (CDSs) and 17326 simple sequence repeats (SSRs). Saturation analyses were applied to evaluate the accuracy of gene expression and we hypothesized that the detection of gene expression might be effective. Finally, single-copy orthologous genes were applied to construct the phylogenetic relationship of K. punctatus. Results showed that K. punctatus diverged from the common ancestor with Alosa alosa, Alosa pseudoharengus and Sardina pilchardus at approx. 61.16–92.52 MYA. The present study will provide a foundational molecular information for the biological research of K. punctatus.

A Comprehensive Review on Crustaceans' Immune System With a Focus on Freshwater Crayfish in Relation to Crayfish Plague Disease

Freshwater crayfish immunity has received great attention due to the need for urgent conservation. This concern has increased the understanding of the cellular and humoral defense systems, although the regulatory mechanisms involved in these processes need updating. There are, however, aspects of the immune response that require clarification and integration. The particular issues addressed in this review include an overall description of the oomycete Aphanomyces astaci, the causative agent of the pandemic plague disease, which affects freshwater crayfish, and an overview of crustaceans' immunity with a focus on freshwater crayfish. It includes a classification system of hemocyte sub-types, the molecular factors involved in hematopoiesis and the differential role of the hemocyte subpopulations in cell-mediated responses, including hemocyte infiltration, inflammation, encapsulation and the link with the extracellular trap cell death pathway (ETosis). In addition, other topics discussed include the identity and functions of hyaline cells, the generation of neoplasia, and the emerging topic of the role of sessile hemocytes in peripheral immunity. Finally, attention is paid to the molecular execution of the immune response, from recognition by the pattern recognition receptors (PRRs), the role of the signaling network in propagating and maintaining the immune signals, to the effector elements such as the putative function of the Down syndrome adhesion molecules (Dscam) in innate immune memory.

Integrated Full-Length Transcriptome and RNA-Seq to Identify Immune System Genes from the Skin of Sperm Whale (Physeter macrocephalus)

Cetaceans are a group of secondary aquatic mammals whose ancestors returned to the ocean from land, and during evolution, their immune systems adapted to the aquatic environment. Their skin, as the primary barrier to environmental pathogens, supposedly evolved to adapt to a new living environment. However, the immune system in the skin of cetaceans and the associated molecular mechanisms are still largely unknown. To better understand the immune system, we extracted RNA from the sperm whale's (Physeter macrocephalus) skin and performed PacBio full-length sequencing and RNA-seq sequencing. We obtained a total of 96,350 full-length transcripts with an average length of 1705 bp and detected 5150 genes that were associated with 21 immune-related pathways by gene annotation enrichment analysis. Moreover, we found 89 encoding genes corresponding to 33 proteins were annotated in the NOD-like receptor (NLR)-signaling pathway, including NOD1, NOD2, RIP2, and NF-kB genes, which were discussed in detail and predicted to play essential roles in the immune system of the sperm whale. Furthermore, NOD1 was highly conservative during evolution by the sequence comparison and phylogenetic tree. These results provide new information about the immune system in the skin of cetaceans, as well as the evolution of immune-related genes.

Iso-seq analysis and functional annotation of the Santa Fe cave crayfish (Procambarus erythrops) transcriptome

The genus Procambarus represents a diverse genus of freshwater crayfish that includes epigean species, stygobitic species, and at least one parthenogenic species. Despite its evolutionary, ecological, and economic importance, most genomic and transcriptomic resources for this genus are limited to a couple of model species. We sequenced the transcriptome of a non-model species, P. erythrops, a geographically restricted stygobitic species from Florida. RNA isolated from gill, muscle and eye tissue was pooled to create a de novo transcriptome assembly using Single Molecule Real-Time sequencing (PacBio), resulting in 19,442 full-length isoforms. The assembly has been deposited in the NCBI (BioProject PRJNA657230). These data will make an important contribution to the comparative study of transcriptome evolution in crayfish and crustaceans.

Exposure to benzo[a]pyrene triggers distinct patterns of microRNA transcriptional profiles in aquatic firefly Aquatica wuhana (Coleoptera: Lampyridae)

Larval aquatic fireflies in fresh water are adversely affected by water pollutants such as benzo(a)pyrene (BaP). However, their response to BaP stress at the microRNA (miRNA)-regulatory level remains unknown. Here, transcriptomes containing 31,872 genes and six miRNA transcriptional profiles were obtained for Aquatica wuhana larvae, and comparative analysis was performed between larvae exposed to BaP (0.01 mg/L) and unexposed controls. Fifteen of 114 miRNAs identified via bioinformatics were detected as differentially expressed (DEMs) upon BaP exposure. Analysis results of predicted target genes of DEM suggests that BaP exposure primarily triggered transcriptional changes of miRNA associated with five major regulatory categories: 1) osmotic balance, 2) energy metabolic efficiency, 3) development, 4) xenobiotic metabolism (oxidative stress), and 5) innate immune response. Based on six innate immune- and xenobiotic metabolism-related pathways enriched by the predicted DEM targets, 11 key BaP-responsive DEMs were further screened to investigate dynamic changes of expression in response to BaP stress at five time points, and also to validate the miRNA sequencing data using quantitative real-time PCR. This study provides valuable information for the protection of firefly resources and supplements the understanding of miRNA regulatory mechanisms in response to water deterioration.

Miocene Diversification and High-Altitude Adaptation of Parnassius Butterflies (Lepidoptera: Papilionidae) in Qinghai-Tibet Plateau Revealed by Large-Scale Transcriptomic Data

Simple Summary

Parnassius butterflies have contributed to fundamental studies in biogeography, insect–plant interactions, and other fields of conservation biology and ecology. However, the early evolutionary pattern and molecular adaptation mechanism of this alpine butterfly group to high altitudes in Qinghai–Tibet Plateau are poorly understood up to now. In this study, we report for the first time, a relatively large-scale transcriptomic dataset of eight Parnassius species and their two closely related papilionid species, a dated phylogeny based on hundreds of gene sequences, and potential genetic mechanisms underlying the high-altitude adaptation by investigating changes in evolutionary rates and positively selected genes. Overall, our findings indicate that the transcriptome data sets reported here can provide some new insights into the spatiotemporally evolutionary pattern and high altitude adaptation of Parnassius butterflies from the extrinsic and intrinsic view, and will support further expressional and functional studies that will help interested researchers to address evolution, biodiversity and conservation questions concerning Parnassius and other butterfly species.

Abstract

The early evolutionary pattern and molecular adaptation mechanism of alpine Parnassius butterflies to high altitudes in Qinghai–Tibet Plateau are poorly understood up to now, due to difficulties in sampling, limited sequence data, and time calibration issues. Here, we present large-scale transcriptomic datasets of eight representative Parnassius species to reveal the phylogenetic timescale and potential genetic basis for high-altitude adaptation with multiple analytic strategies using 476 orthologous genes. Our phylogenetic results strongly supported that the subgenus Parnassius formed a well-resolved basal clade, and the subgenera Tadumia and Kailasius were closely related in the phylogenetic trees. In addition, molecular dating analyses showed that the Parnassius began to diverge at about 13.0 to 14.3 million years ago (middle Miocene), correlated with their hostplant’s spatiotemporal distributions, as well as geological and palaeoenvironmental changes of the Qinghai–Tibet Plateau. Moreover, the accelerated evolutionary rate, candidate positively selected genes and their potentially functional changes were detected, probably contributed to the high-altitude adaptation of Parnassius species. Overall, our study provided some new insights into the spatiotemporally evolutionary pattern and high altitude adaptation of Parnassius butterflies from the extrinsic and intrinsic view, which will help to address evolution, biodiversity, and conservation questions concerning Parnassius and other butterfly species.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

Genetic compensation rather than genetic assimilation drives the evolution of plasticity in response to mild warming across latitudes in a damselfly

Global warming is causing plastic and evolutionary changes in the phenotypes of ectotherms. Yet, we have limited knowledge on how the interplay between plasticity and evolution shapes thermal responses and underlying gene expression patterns. We assessed thermal reaction norm patterns across the transcriptome and identified associated molecular pathways in northern and southern populations of the damselfly Ischnura elegans. Larvae were reared in a common garden experiment at the mean summer water temperatures experienced at the northern (20°C) and southern (24°C) latitudes. This allowed a space-for-time substitution where the current gene expression levels at 24°C in southern larvae are a proxy for the expected responses of northern larvae under gradual thermal evolution to the predicted 4°C warming. Most differentially expressed genes showed fixed differences across temperatures between latitudes, suggesting that thermal genetic adaptation will mainly evolve through changes in constitutive gene expression. Northern populations also frequently showed plastic responses in gene expression to mild warming, while southern populations were much less responsive to temperature. Thermal responsive genes in northern populations showed to a large extent a pattern of genetic compensation, namely gene expression that was induced at 24°C in northern populations remained at a lower constant level in southern populations, and were associated with metabolic and translation pathways. There was instead little evidence for genetic assimilation of an initial plastic response to mild warming. Our data therefore suggest that genetic compensation rather than genetic assimilation may drive the evolution of plasticity in response to mild warming in this damselfly species.

Identification and functional characterization of the astacidin family of proline-rich host defence peptides (PcAst) from the red swamp crayfish (Procambarus clarkii, Girard 1852)

This study reports the identification of four novel proline-rich antimicrobial peptides (PR-AMP) from the transcriptome of the red swamp crayfish Procambarus clarkii. The newly identified putative peptides (PcAst-1b, -1c, -2 and -3), which are related with the previously identified hemocyte-specific PR-AMP astacidin-1, are encoded by the multi-genic astacidin gene family. The screening of available and proprietary transcriptomes allowed to define the taxonomical range of distribution of this gene family to Astacoidea and Parastacoidea. The antimicrobial properties of three synthetic PcAst peptides (PcAst-1a, -1b/c and -2), were characterized against reference bacteria or multidrug resistant clinical isolates, and their cytotoxicity was evaluated towards human transformed cell lines. The antimicrobial activity ranged from potent and broad-spectrum, in low-salt medium, to poor, whereas it was generally low in full nutrient broth. No significant toxic effects were observed on cultured human cells. RNA-seq data from 12 different tissues indicated a strong specificity for haemocytes under naïve physiological condition, with moderate expression (5-fold lower) in gills. Quantitative real time PCR revealed a rapid (within 2 h) and significant up-regulation of PcAst-1a (Astacidin 1) and PcAst-2 expression in response to LPS injection. Due to the variation in antimicrobial potency and inducibility, the roles of the other astacidins (PcAst-1b, -1c and -3) need to be further investigated to determine their significance to the immune responses of the red swamp crayfish.

Disruption of the lactate dehydrogenase and acetate kinase genes in Klebsiella pneumoniae HD79 to enhance 2,3-butanediol production, and related transcriptomics analysis

OBJECTIVES

2,3-Butanediol (2,3-BD) is widely used in several chemical syntheses as well as the manufacture of plastics, solvents, and antifreeze formulations, and can be manufactured by microbial glucose fermentation. Conventional (2,3-BD) fermentation typically has low productivity, yield, and purity, and is expensive for commercial applications. We aimed to delete the lactate dehydrogenase and acetate kinase (ldhA and ack) genes in Klebsiella pneumoniae HD79 by using λRed homologous recombination technology, to eliminate by-products and thereby improve (2,3-BD) production. We also analyzed the resulting gene changes by using transcriptomics.

RESULTS

The yield of (2,3-BD) from the mutant Klebsiella strain was 46.21 g/L, the conversion rate was 0.47 g/g, and the productivity was 0.64 g/L·h, which represented increases of 54.9%, 20.5%, and 106.5% respectively, compared to (WT) strains. Lactate and acetate decreased by 48.2% and 62.8%, respectively. Transcriptomics analysis showed that 4628 genes were differentially expressed (404 significantly up-regulated and 162 significantly down-regulated). Moreover, the (2,3-BD) operon genes were differentially expressed.

CONCLUSION

Our data showed that the biosynthesis of (2,3-BD) was regulated by inducers (lactate and acetate), a regulator (BudR), and carbon flux. Elimination of acidic by-products by ldhA and ack knockdown significantly improved (2,3-BD) production. Our results provide a deeper understanding of the mechanisms underlying (2,3-BD) production, and form a molecular basis for the improvement this process by genetic modification in the future.

New insight into the molecular basis of Fe (III) stress responses of Procambarus clarkii by transcriptome analysis

Iron in excess can have toxic effects on living organisms. In China, the freshwater crayfish Procambarus clarkii is a source of aquatic food with high-quality protein and has significant commercial value. P. clarkii shows oxidative stress on exposure to heavy metals, and antioxidant enzymes, such as ubiquitination enzymes and proteasomes, play important roles in oxidative stress. To understand the antioxidant defense system of P. clarkii, we analyzed the hepatopancreas transcriptomes of P. clarkii after stimulation with FeCl₃. In total, 5199 differentially expressed genes (DEGs) were identified (2747 upregulated and 2452 downregulated). GO analysis revealed that these DEGs belonged to 16 cellular component, 16 molecular function, and 19 biological process subcategories. A total of 1069 DEGs were classified into 25 categories by using COG. Some antioxidant defense pathways, such as "Ubiquitin mediated proteolysis" and "Glutathione metabolism," were identified using KEGG. In addition, quantitative real time-PCR (qRT-PCR) substantiated the up-regulation of a random selection of DEGs including antioxidant and immune defense genes. We obtained information for P. clarkii transcriptome databases and new insights into the responses of P. clarkii hepatopancreas to heavy metals.

Arctic charr brain transcriptome strongly affected by summer seasonal growth but only subtly by feed deprivation

Background

The Arctic charr (Salvelinus alpinus) has a highly seasonal feeding cycle that comprises long periods of voluntary fasting and a short but intense feeding period during summer. Therefore, the charr represents an interesting species for studying appetite-regulating mechanisms in fish.

Results

In this study, we compared the brain transcriptomes of fed and feed deprived charr over a 4 weeks trial during their summer feeding season. Despite prominent differences in body condition between fed and feed deprived charr at the end of the trial, feed deprivation affected the brain transcriptome only slightly. In contrast, the transcriptome differed markedly over time in both fed and feed deprived charr, indicating strong shifts in basic cell metabolic processes possibly due to season, growth, temperature, or combinations thereof. The GO enrichment analysis revealed that many biological processes appeared to change in the same direction in both fed and feed deprived fish. In the feed deprived charr processes linked to oxygen transport and apoptosis were down- and up-regulated, respectively. Known genes encoding for appetite regulators did not respond to feed deprivation. Gene expression of Deiodinase 2 (DIO2), an enzyme implicated in the regulation of seasonal processes in mammals, was lower in response to season and feed deprivation. We further found a higher expression of VGF (non-acronymic) in the feed deprived than in the fed fish. This gene encodes for a neuropeptide associated with the control of energy metabolism in mammals, and has not been studied in relation to regulation of appetite and energy homeostasis in fish.

Conclusions

In the Arctic charr, external and endogenous seasonal factors for example the increase in temperature and their circannual growth cycle, respectively, evoke much stronger responses in the brain than 4 weeks feed deprivation. The absence of a central hunger response in feed deprived charr give support for a strong resilience to the lack of food in this high Arctic species. DIO2 and VGF may play a role in the regulation of energy homeostasis and need to be further studied in seasonal fish.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5874-z) contains supplementary material, which is available to authorized users.

Comparative transcriptomic analysis provides insights into the response to the benzo(a)pyrene stress in aquatic firefly (Luciola leii)

Many studies have reported that behavior and bioluminescence of fireflies could be affected by changes in environment conditions. However, little is known about how the deterioration of the aquatic environment affects aquatic fireflies, particularly with respect to molecular responses following exposure to water pollutants, such as benzo(a)pyrene (BaP), which is a key indicator in environmental risk assessment because of the hazards it poses. Here, whole transcriptome sequencing and gene expression analysis were performed on freshwater fireflies (Luciola leii) exposed to BaP (concentration of 0.01 mg/L). Four transcriptomic libraries were constructed for the control and treatment groups, including two biological replicates. From the mixed pools (each pool contains 60 individuals from three time points), a total of 54,282 unigenes were assembled. Furthermore, 329,337 of Single-nucleotide Polymorphisms (SNPs) and 1324 of Simple Sequence Repeats (SSRs) were predicted using bioinformatics, which is useful for the future development of molecular markers. Subsequently, 2414 differently expressed genes (DEGs) were identified in response to BaP stress in comparison to the control, including 1350 up-regulated and 1064 down-regulated DEGs. Functional enrichment showed that these DEGs are primarily related to innate immunity; xenobiotic biodegradation and response, biomacromolecule metabolism, biosynthesis, and absorption. Eight key BaP-responsive DEGs were screened to survey the dynamic changes of expression in response to BaP stress at different time points, and to validate the RNA sequencing data using quantitative real-time PCR. The results indicate that the expression of genes encoding UGT, CYP3A, CYP9, CYP6AS5 and ADHP were induced, while those encoding UGT2B10L, PTGDS, and ALDH were reduced, to participate in response to the BaP exposure and potentially help counteract the adverse effects of BaP. This investigation provides insight into the toxicological response of fireflies to the occurrence of water deterioration.

De novo characterization of placental transcriptome in the Eurasian beaver (Castor fiber L.)

Our pioneering data provide the first comprehensive view of placental transcriptome of the beaver during single and multiple gestation. RNA-Seq and a de novo approach allowed global pattern identification of C. fiber placental transcriptome. Non-redundant beaver transcriptome comprised 211,802,336 nt of placental transcripts, grouped into 128,459 contigs and clustered into 83,951 unigenes. An Ensembl database search revealed 14,487, 14,994, 15,004, 15,267 and 15,892 non-redundant homologs for Ictidomys tridecemlineatus, Rattus norvegicus, Mus musculus, Homo sapiens and Castor canadensis, respectively. Due to expression levels, the identified transcripts were divided into two sets: non-redundant and highly expressed (FPKM > 2 in at least three examined samples), analysed simultaneously. Among 17,009 highly expressed transcripts, 12,147 had BLASTx hits. GO annotations (175,882) were found for 4301 transcripts that were assigned to biological process (16,386), cellular component (9149) and molecular function (8338) categories; 666 unigenes were also classified into 122 KEGG pathways. Comprehensive analyses were performed for 411 and 3078 highly expressed transcripts annotated with a list of processes linked to ‘placenta’ (31 GO terms) or ‘embryo’ (324 GO terms), respectively. Among transcripts with entire CDS annotation, 281 (placenta) and 34 (embryo) alternative splicing events were identified. A total of 8499 putative SNVs (~ 6.2 SNV/transcript and 1.7 SNV/1 kb) were predicted with 0.1 minimum frequency and maximum variant quality (p value 10e⁻⁹). Our results provide a broad-based characterization of the global expression pattern of the beaver placental transcriptome. Enhancement of transcriptomic resources for C. fiber should improve understanding of crucial pathways relevant to proper placenta development and successful reproduction.

The transcriptome of the marine calanoid copepod Temora longicornis under heat stress and recovery

Understanding the impacts of global change in zooplankton communities is crucial, as alterations in the zooplankton communities can affect entire marine ecosystems. Despite the economic and ecological importance of the calanoid copepod Temora longicornis in the Belgian part of the North Sea, molecular data is still very limited for this species. Using HiSeq Illumina sequencing, we sequenced the whole transcriptome of T. longicornis, after being exposed to realistic temperatures of 14 and 17 °C. After both an acute (1 day) and a more sustained (5 days) thermal exposure to 17 °C, we investigated gene expression differences with animals exposed to 14 °C, which may be critical for the thermal acclimation and resilience of this copepod species. We also studied the possibility of a short term stress recovery of a heat shock. A total of 179,569 transcripts were yielded, of which 44,985 putative ORF transcripts were identified. These transcripts were subsequently annotated into roughly 22,000 genes based on known sequences using Gene Ontology (GO) and KEGG databases. Temora only showed a mild response to both the temperature and the duration of the exposure. We found that the expression of 27 transcripts varied significantly with an increase in temperature of 3 °C, of which eight transcripts were differentially expressed after acute exposure only. Gene set enrichment analysis revealed that, overall, T. longicornis was more impacted by a sustained thermal exposure, rather than an immediate (acute) exposure, with two times as many enriched GO terms in the sustained treatment. We also identified several general stress responses independent of exposure time, such as modified protein synthesis, energy mobilisation, cuticle and chaperone proteins. Finally, we highlighted candidate genes of a possible recovery from heat exposure, identifying similar terms as those enriched in the heat treatments, i.e. related to for example energy metabolism, cuticle genes and extracellular matrix. The data presented in this study provides the first transcriptome available for T. longicornis which can be used for future genomic studies.

De novo annotation of the transcriptome of the Northern Wheatear (Oenanthe oenanthe)

We have sequenced a partial transcriptome of the Northern Wheatear (Oenanthe oenanthe), a species with one of the longest migrations on Earth. The transcriptome was constructed de novo using RNA-Seq sequence data from the pooled mRNA of six different tissues: brain, muscle, intestine, liver, adipose tissue and skin. The samples came from nine captive-bred wheatears collected at three different stages of the endogenous autumn migratory period: (1) lean birds prior the onset of migration, (2) during the fattening stage and (3) individuals at their migratory body mass plateau, when they have almost doubled their lean body mass. The sample structure used to build up the transcriptome of the Northern Wheatears concerning tissue composition and time guarantees the future survey of the regulatory genes involved in the development of the migratory phenotype. Through the pre-migratory period, birds accomplish outstanding physical and behavioural changes that involve all organ systems. Nevertheless, the molecular mechanisms through which birds synchronize and control hyperphagia, fattening, restlessness increase, immunity boosting and tuning the muscles for such endurance flight are still largely unknown. The use of RNA-Seq has emerged as a powerful tool to analyse complex traits on a broad scale, and we believe it can help to characterize the migratory phenotype of wheatears at an unprecedented level. The primary challenge to conduct quantitative transcriptomic studies in non-model species is the availability of a reference transcriptome, which we have constructed and described in this paper. The cDNA was sequenced by pyrosequencing using the Genome Sequencer Roche GS FLX System; with single paired-end reads of about 400 bp. We estimate the total number of genes at 15,640, of which  67% could be annotated using Turkey and Zebra Finch genomes, or protein sequence information from SwissProt and NCBI databases. With our study, we have made a first step towards understanding the migratory phenotype regarding gene expression of a species that has become a model to study birds long-distance migrations.

How crickets become freeze tolerant: The transcriptomic underpinnings of acclimation in Gryllus veletis

Some ectotherms can survive internal ice formation. In temperate regions, freeze tolerance is often induced by decreasing temperature and/or photoperiod during autumn. However, we have limited understanding of how seasonal changes in physiology contribute to freeze tolerance, and how these changes are regulated. During a six week autumn-like acclimation, late-instar juveniles of the spring field cricket Gryllus veletis (Orthoptera: Gryllidae) become freeze tolerant, which is correlated with accumulation of low molecular weight cryoprotectants, elevation of the temperature at which freezing begins, and metabolic rate suppression. We used RNA-Seq to assemble a de novo transcriptome of this emerging laboratory model for freeze tolerance research. We then focused on gene expression during acclimation in fat body tissue due to its role in cryoprotectant production and regulation of energetics. Acclimated G. veletis differentially expressed >3000 transcripts in fat body. This differential expression may contribute to metabolic suppression in acclimated G. veletis, but we did not detect changes in expression that would support cryoprotectant accumulation or enhanced control of ice formation, suggesting that these latter processes are regulated post-transcriptionally. Acclimated G. veletis differentially regulated transcripts that likely coordinate additional freeze tolerance mechanisms, including upregulation of enzymes that may promote membrane and cytoskeletal remodelling, cryoprotectant transporters, cytoprotective proteins, and antioxidants. Thus, while accumulation of cryoprotectants and controlling ice formation are commonly associated with insect freeze tolerance, our results support the hypothesis that many other systems contribute to surviving internal ice formation. Together, this information suggests new avenues for understanding the mechanisms underlying insect freeze tolerance.

Characterization of ladybird Henosepilachna vigintioctopunctata transcriptomes across various life stages

Henosepilachna vigintioctopunctata is a vegetable pest that has spread worldwide. It belongs to the Coccinellidae family, whose members exhibit remarkable diversity, both in terms of their diets and the colored spots that appear on the elytra in the adult stage. Transcriptomic data from H. vigintioctopunctata at different life stages would be useful for further investigating the genetic basis of this dietary diversity and the formation of the colored spots in ladybird beetles, as well as revealing the population dynamics of H. vigintioctopunctata, which could be useful in pest control. Here, we generated a comprehensive RNA-seq data set (a total of ~24 Gb of clean data) for H. vigintioctopunctata by sequencing samples collected at different life stages. We characterized the transcriptomes of each of the four life stages (egg, larva, pupa, adult) and generated a high-coverage pool by combining all the RNA-seq reads. Furthermore, we identified a catalog of simple sequence repeat (SSR) markers. This represents the first study to collect transcriptome data from all life stages of a ladybird beetle.

Crayfish immunity - Recent findings

Freshwater crayfish is an important commodity as well as a successful model for studies on crustacean immunity. Due to the ease with which they are kept and the available methods for hemocyte separation and culture they have proven to be very useful. Here, recent progress regarding pattern recognition, immune effector production and antiviral mechanisms are discussed. Several cases of functional resemblance between vertebrate complement and the crayfish immune reactions are highlighted.

De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species

Background

Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date.

Results

We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species.

Conclusion

New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.

Comparative transcriptomic analysis of Tibetan Gynaephora to explore the genetic basis of insect adaptation to divergent altitude environments

Adaptation of insects to different altitudes remain largely unknown, especially those endemic to the Tibetan Plateau (TP). Here, we generated the transcriptomes of Gynaephora menyuanensis and G. alpherakii, inhabiting different high altitudes on the TP, and used these and the previously available transcriptomic and genomic sequences from low-altitude insects to explore potential genetic basis for divergent high-altitude adaptation in Gynaephora. An analysis of 5,869 orthologous genes among Gynaephora and other three low-altitude insects uncovered that fast-evolving genes and positively selected genes (PSGs) in the two Gynaephora species were enriched in energy metabolism and hypoxia response categories (e.g. mitochondrion, oxidation-reduction process, and response to oxidative stress). Particularly, mTOR signaling pathway involving hypoxia was enriched by PSGs, indicating this well-known pathway in mammal hypoxia adaptation may be an important signaling system in Gynaephora. Furthermore, some PSGs were associated with response to hypoxia (e.g. cytochrome proteins), cold (e.g. dehydrogenase) and DNA repair (e.g. DNA repair proteins). Interestingly, several insect-specific genes that were associated with exoskeleton and cuticle development (e.g. chitinase and ecdysteroids) had experienced positive selection, suggesting the specific adaptive mechanisms in insects. This study is favourable for understanding the adaptive evolution of Gynaephora and even TP insects to divergent altitudes.

Transcriptome analysis of hepatopancraes of Procambarus clarkii challenged with polyriboinosinic polyribocytidylic acid (poly I:C)

Crustacean hepatopancreas regulates metabolic processes, biogenesis and innate immune processes, and the knowledge on its immune genes are crucial to understand antimicrobial mechanisms. In this study, we reported the transcriptomic profile of Procambarus clarkii hepatopancreas after poly I:C administration using high-throughput sequencing. Following de novo assembly 56,716 unigene sequences with an average length of 810 bp was obtained. The unigene sequences were annotated to three ontologies including cellular components, biological processes and molecular functions, further 56,716 unigene sequences were mapped to 25 COG categories. A total of 2497 differentially expressed genes (DEGs) were identified following the comparative analysis between poly I:C treated and control group, and then KEGG enrichment analysis were performed to detect immune related pathways. Quantitative real time polymerase chain reaction showed that the selected DEGs significantly up-regulated following poly I:C administration in comparison to control group. The transcriptomic sequence information will improve the knowledge of this economically important crustacean, and will shed light on its antiviral immune mechanisms.

De novo assembly and annotation of the whole transcriptome of Oratosquilla oratoria

As the representative species of Stomatopoda, Oratosquilla oratoria can tolerate the complex benthic environment. In this study, our goal was to develop the transcriptomic resource of O. oratoria that would support adaptation mechanism studies. We generated the whole transcriptome of O. oratoria from combined tissues (eyestalk, muscle, sexual and viscus) using Hiseq technology. A total of 51,305,284 high-quality clean reads were assembled to produce 59,054 non-redundant transcripts with a mean length of 987nt using Trinity and CORSET software. Among the predictable unigenes, a total of 32,451 unigenes were annotated based on protein databases. Finally, we predicted the coding sequences of 31,822 unigenes and obtained 19,057 SSRs in the present study. The present study will provide an important resource and foundational understanding for future genomic research of O. oratoria.

Sex and tissue specific gene expression patterns identified following de novo transcriptomic analysis of the Norway lobster, Nephrops norvegicus

BACKGROUND

The Norway lobster, Nephrops norvegicus, is economically important in European fisheries and is a key organism in local marine ecosystems. Despite multi-faceted scientific interest in this species, our current knowledge of genetic resources in this species remains very limited. Here, we generated a reference de novo transcriptome for N. norvegicus from multiple tissues in both sexes. Bioinformatic analyses were conducted to detect transcripts that were expressed exclusively in either males or females. Patterns were validated via RT-PCR.

RESULTS

Sixteen N. norvegicus libraries were sequenced from immature and mature ovary, testis and vas deferens (including the masculinizing androgenic gland). In addition, eyestalk, brain, thoracic ganglia and hepatopancreas tissues were screened in males and both immature and mature females. RNA-Sequencing resulted in >600 million reads. De novo assembly that combined the current dataset with two previously published libraries from eyestalk tissue, yielded a reference transcriptome of 333,225 transcripts with an average size of 708 base pairs (bp), with an N50 of 1272 bp. Sex-specific transcripts were detected primarily in gonads followed by hepatopancreas, brain, thoracic ganglia, and eyestalk, respectively. Candidate transcripts that were expressed exclusively either in males or females were highlighted and the 10 most abundant ones were validated via RT-PCR. Among the most highly expressed genes were Serine threonine protein kinase in testis and Vitellogenin in female hepatopancreas. These results align closely with gene annotation results. Moreover, a differential expression heatmap showed that the majority of differentially expressed transcripts were identified in gonad and eyestalk tissues. Results indicate that sex-specific gene expression patterns in Norway lobster are controlled by differences in gene regulation pattern between males and females in somatic tissues.

CONCLUSIONS

The current study presents the first multi-tissue reference transcriptome for the Norway lobster that can be applied to future biological, wild restocking and fisheries studies. Sex-specific markers were mainly expressed in males implying that males may experience stronger selection than females. It is apparent that differential expression is due to sex-specific gene regulatory pathways that are present in somatic tissues and not from effects of genes located on heterogametic sex chromosomes. The N. norvegicus data provide a foundation for future gene-based reproductive studies.

Transcriptome profiling of ontogeny in the acridid grasshopper Chorthippus biguttulus

Acridid grasshoppers (Orthoptera:Acrididae) are widely used model organisms for developmental, evolutionary, and neurobiological research. Although there has been recent influx of orthopteran transcriptomic resources, many use pooled ontogenetic stages obscuring information about changes in gene expression during development. Here we developed a de novo transcriptome spanning 7 stages in the life cycle of the acridid grasshopper Chorthippus biguttulus. Samples from different stages encompassing embryonic development through adults were used for transcriptomic profiling, revealing patterns of differential gene expression that highlight processes in the different life stages. These patterns were validated with semi-quantitative RT-PCR. Embryonic development showed a strongly differentiated expression pattern compared to all of the other stages and genes upregulated in this stage were involved in signaling, cellular differentiation, and organ development. Our study is one of the first to examine gene expression during post-embryonic development in a hemimetabolous insect and we found that only the fourth and fifth instars had clusters of genes upregulated during these stages. These genes are involved in various processes ranging from synthesis of biogenic amines to chitin binding. These observations indicate that post-embryonic ontogeny is not a continuous process and that some instars are differentiated. Finally, genes upregulated in the imago were generally involved in aging and immunity. Our study highlights the importance of looking at ontogeny as a whole and indicates promising directions for future research in orthopteran development.

Effects of cold-acclimation on gene expression in Fall field cricket (Gryllus pennsylvanicus) ionoregulatory tissues

Background

Cold tolerance is a key determinant of temperate insect distribution and performance. Chill-susceptible insects lose ion and water homeostasis during cold exposure, but prior cold acclimation improves both cold tolerance and defense of homeostasis. The mechanisms underlying these processes are mostly unknown; cold acclimation is thought to enhance ion transport in the cold and/or prevent leak of water and ions. To identify candidate mechanisms of cold tolerance plasticity we generated transcriptomes of ionoregulatory tissues (hindgut and Malpighian tubules) from Gryllus pennsylvanicus crickets and compared gene expression in warm- and cold-acclimated individuals.

Results

We assembled a G. pennsylvanicus transcriptome de novo from 286 million 50-bp reads, yielding 70,037 contigs (~44% of which had putative BLAST identities). We compared the transcriptomes of warm- and cold-acclimated hindguts and Malpighian tubules. Cold acclimation led to a ≥ 2-fold change in the expression of 1493 hindgut genes (733 downregulated, 760 upregulated) and 2008 Malpighian tubule genes (1009 downregulated, 999 upregulated). Cold-acclimated crickets had altered expression of genes putatively associated with ion and water balance, including: a downregulation of V-ATPase and carbonic anhydrase in the Malpighian tubules and an upregulation of Na⁺-K⁺ ATPase in the hindgut. We also observed acclimation-related shifts in the expression of cytoskeletal genes in the hindgut, including actin and actin-anchoring/stabilizing proteins, tubulin, α-actinin, and genes involved in adherens junctions organization. In both tissues, cold acclimation led to differential expression of genes encoding cytochrome P450s, glutathione-S-transferases, apoptosis factors, DNA repair, and heat shock proteins.

Conclusions

This is the first G. pennsylvanicus transcriptome, and our tissue-specific approach yielded new candidate mechanisms of cold tolerance plasticity. Cold acclimation may reduce loss of hemolymph volume in the cold by 1) decreasing primary urine production via reduced expression of carbonic anhydrase and V-ATPase in the Malpighian tubules and 2) by increasing Na⁺ (and therefore water) reabsorption across the hindgut via increase in Na⁺-K⁺ ATPase expression. Cold acclimation may reduce chilling injury by remodeling and stabilizing the hindgut epithelial cytoskeleton and cell-to-cell junctions, and by increasing the expression of genes involved in DNA repair, detoxification, and protein chaperones.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3711-9) contains supplementary material, which is available to authorized users.

Loggerhead sea turtle embryos (Caretta caretta) regulate expression of stress response and developmental genes when exposed to a biologically realistic heat stress

Oviparous reptile embryos are expected to breach their critical thermal maxima if temperatures reach those predicted under current climate change models due to the lack of the maternal buffering processes and parental care. Heat-shock proteins (HSPs) are integral in the molecular response to thermal stress, and their expression is heritable, but the roles of other candidate families such as the heat-shock factors (HSFs) have not been determined in reptiles. Here, we subject embryonic sea turtles (Caretta caretta) to a biologically realistic thermal stress and employ de novo transcriptomic profiling of brain tissue to investigate the underlying molecular response. From a reference transcriptome of 302 293 transcripts, 179 were identified as differentially expressed between treatments. As anticipated, genes enriched in the heat-shock treatment were primarily associated with the Hsp families, or were genes whose products play similar protein editing and chaperone functions (e.g. bag3, MYOC and serpinh1). Unexpectedly, genes encoding the HSFs were not significantly upregulated under thermal stress, indicating their presence in unstressed cells in an inactive state. Genes that were downregulated under thermal stress were less well functionally defined but were associated with stress response, development and cellular organization, suggesting that developmental processes may be compromised at realistically high temperatures. These results confirm that genes from the Hsp families play vital roles in the thermal tolerance of developing reptile embryos and, in addition with a number of other genes, should be targets for evaluating the capacity of oviparous reptiles to respond adaptively to the effects of climate change.

De Novo Transcriptome Characterization of a Sterilizing Trematode Parasite (Microphallus sp.) from Two Species of New Zealand Snails

Snail-borne trematodes represent a large, diverse, and evolutionarily, ecologically, and medically important group of parasites, often imposing strong selection on their hosts and causing host morbidity and mortality. Even so, there are very few genomic and transcriptomic resources available for this important animal group. We help to fill this gap by providing transcriptome resources from trematode metacercariae infecting two congeneric snail species, Potamopyrgus antipodarum and P. estuarinus. This genus of New Zealand snails has gained prominence in large part through the development of P. antipodarum and its sterilizing trematode parasite Microphallus livelyi into a textbook model for host–parasite coevolutionary interactions in nature. By contrast, the interactions between Microphallus trematodes and P. estuarinus, an estuary-inhabiting species closely related to the freshwater P. antipodarum, are relatively unstudied. Here, we provide the first annotated transcriptome assemblies from Microphallus isolated from P. antipodarum and P. estuarinus. We also use these transcriptomes to produce genomic resources that will be broadly useful to those interested in host–parasite coevolution, local adaption, and molecular evolution and phylogenetics of this and other snail–trematode systems. Analyses of the two Microphallus transcriptomes revealed that the two trematode types are more genetically differentiated from one another than are the M. livelyi infecting different populations of P. antipodarum, suggesting that the Microphallus infecting P. estuarinus represent a distinct lineage. We also provide a promising set of candidate genes likely involved in parasitic infection and response to salinity stress.

De novo assembly, annotation, and characterization of the whole brain transcriptome of male and female Syrian hamsters

Hamsters are an ideal animal model for a variety of biomedical research areas such as cancer, virology, circadian rhythms, and behavioural neuroscience. The use of hamsters has declined, however, most likely due to the dearth of genetic tools available for these animals. Our laboratory uses hamsters to study acute social stress, and we are beginning to investigate the genetic mechanisms subserving defeat-induced behavioural change. We have been limited, however, by the lack of genetic resources available for hamsters. In this study, we sequenced the brain transcriptome of male and female Syrian hamsters to generate the necessary resources to continue our research. We completed a de novo assembly and after assembly optimization, there were 113,329 transcripts representing 14,530 unique genes. This study is the first to characterize transcript expression in both female and male hamster brains and offers invaluable information to promote understanding of a host of important biomedical research questions for which hamsters are an excellent model.

Deep sequencing of transcriptomes from the nervous systems of two decapod crustaceans to characterize genes important for neural circuit function and modulation

BACKGROUND

Crustaceans have been studied extensively as model systems for nervous system function from single neuron properties to behavior. However, lack of molecular sequence information and tools have slowed the adoption of these physiological systems as molecular model systems. In this study, we sequenced and performed de novo assembly for the nervous system transcriptomes of two decapod crustaceans: the Jonah crab (Cancer borealis) and the American lobster (Homarus americanus).

RESULTS

Forty-two thousand, seven hundred sixty-six and sixty thousand, two hundred seventy-three contigs were assembled from C. borealis and H. americanus respectively, representing 9,489 and 11,061 unique coding sequences. From these transcripts, genes associated with neural function were identified and manually curated to produce a characterization of multiple gene families important for nervous system function. This included genes for 34 distinct ion channel types, 17 biogenic amine and 5 GABA receptors, 28 major transmitter receptor subtypes including glutamate and acetylcholine receptors, and 6 gap junction proteins - the Innexins.

CONCLUSION

With this resource, crustacean model systems are better poised for incorporation of modern genomic and molecular biology technologies to further enhance the interrogation of fundamentals of nervous system function.

A transcriptional blueprint for a spiral-cleaving embryo

Background

The spiral cleavage mode of early development is utilized in over one-third of all animal phyla and generates embryonic cells of different size, position, and fate through a conserved set of stereotypic and invariant asymmetric cell divisions. Despite the widespread use of spiral cleavage, regulatory and molecular features for any spiral-cleaving embryo are largely uncharted. To address this gap we use RNA-sequencing on the spiralian model Platynereis dumerilii to capture and quantify the first complete genome-wide transcriptional landscape of early spiral cleavage.

Results

RNA-sequencing datasets from seven stages in early Platynereis development, from the zygote to the protrochophore, are described here including the de novo assembly and annotation of ~17,200 Platynereis genes. Depth and quality of the RNA-sequencing datasets allow the identification of the temporal onset and level of transcription for each annotated gene, even if the expression is restricted to a single cell. Over 4000 transcripts are maternally contributed and cleared by the end of the early spiral cleavage phase. Small early waves of zygotic expression are followed by major waves of thousands of genes, demarcating the maternal to zygotic transition shortly after the completion of spiral cleavages in this annelid species.

Conclusions

Our comprehensive stage-specific transcriptional analysis of early embryonic stages in Platynereis elucidates the regulatory genome during early spiral embryogenesis and defines the maternal to zygotic transition in Platynereis embryos. This transcriptome assembly provides the first systems-level view of the transcriptional and regulatory landscape for a spiral-cleaving embryo.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2860-6) contains supplementary material, which is available to authorized users.