Rapid transcriptome and proteome profiling of a non-model marine invertebrate, Bugula neritina

Differential Gene Expression Between Polymorphic Zooids of the Marine Bryozoan Bugulina stolonifera

Bryozoans are a diverse phylum of marine and freshwater colonial invertebrates containing approximately 6,300 described living species. Bryozoans grow by budding new physiologically connected colony members (zooids) from a founding individual that forms from a metamorphosed larva. In some species these zooids come in different shapes and sizes and are specialized to serve different tasks within the colony. A complex interaction of genotype, environment, and developmental pathway shapes zooid fate, however, the specific mechanisms underlying the establishment of this division of labor remain unknown. Here, the first characterization of differential gene expression between polymorphic zooids of a bryozoan colony is presented. The development of different zooid types of lab-cultured Bugulina stolonifera colonies including feeding autozooids, avicularia (derived non-feeding zooids that are homologous to feeding autozooids but shaped like a bird's beak), and rhizoids (a branching network of non-feeding anchoring zooids) was explored using RNA sequencing, de novo transcriptome assembly, and differential gene expression analyses. High throughput sequencing of cDNA libraries yielded an average of 14.9 ± 1.3 (SE) million high-quality paired-end reads per sample. Data for the first de novo transcriptome assemblies of B. stolonifera and the first characterization of genes involved in the formation and maintenance of zooid types within a bryozoan colony are presented. In a comparison between autozooid and avicularium tissues, 1,097 significant differentially expressed genes were uncovered. This work provides a much-needed foundation for understanding the mechanisms involved in the development of polymorphic zooids and the establishment of division of labor in bryozoans.

Proteome of larval metamorphosis induced by epinephrine in the Fujian oyster Crassostrea angulata

Background

The Fujian oyster Crassostrea angulata is an economically important species that has typical settlement and metamorphosis stages. The development of the oyster involves complex morphological and physiological changes, the molecular mechanisms of which are as yet unclear.

Results

In this study, changes in proteins were investigated during larval settlement and metamorphosis of Crassostrea angulata using epinephrine induction. Protein abundance and identity were characterized using label-free quantitative proteomics, tandem mass spectrometry (MS/ MS), and Mascot methods. The results showed that more than 50% (764 out of 1471) of the quantified proteins were characterized as differentially expressed. Notably, more than two-thirds of the differentially expressed proteins were down-regulated in epinephrine-induced larvae. The results showed that “metabolic process” was closely related to the development of settlement and metamorphosis; 5 × 10^− 4 M epinephrine induced direct metamorphosis of larvae and was non-toxic. Calmodulin and MAPK pathways were involved in the regulation of settlement of the oyster. Expression levels of immune-related proteins increased during metamorphosis. Hepatic lectin-like proteins, cadherins, calmodulin, calreticulin, and cytoskeletal proteins were involved in metamorphosis. The nervous system may be remodeled in larval metamorphosis induced by epinephrine. Expression levels of proteins that were enriched in the epinephrine signaling pathway may reflect the developmental stage of the larvae, that may reflect whether or not larvae were directly involved in metamorphosis when the larvae were treated with epinephrine.

Conclusion

The study provides insight into proteins that function in energy metabolism, immune responses, settlement and metamorphosis, and shell formation in C. angulata. The results contribute valuable information for further research on larval settlement and metamorphosis.

Graphical abstract

First transcriptome analysis of bryozoan Fredericella sultana, the primary host of myxozoan parasite Tetracapsuloides bryosalmonae

Bryozoans are aquatic invertebrate moss animals that are found worldwide. Fredericella sultana is a freshwater bryozoan and is the most common primary host of myxozoan parasite, Tetracapsuloides bryosalmonae. However, limited genomic resources are available for this bryozoan, which hampers investigations into the molecular mechanisms of host-parasite interactions. To better understand these interactions, there is a need to build a transcriptome dataset of F. sultana, for functional genomics analysis by large-scale RNA sequencing. Total RNA was extracted from zooids of F. sultana cultivated under controlled laboratory conditions. cDNA libraries were prepared and were analyzed by the Illumina paired-ends sequencing. The sequencing data were used for de novo transcriptome assembly and functional annotation. Approximately 118 million clean reads were obtained, and assembled into 85,544 contigs with an average length of 852 bp, an N50 of 1,085 bp, and an average GC content 51.4%. A total of 23,978 (28%) contigs were annotated using BLASTX analysis. Of these transcripts, 4,400 contigs had highest similarity to brachiopod species Lingula anatina. Based on Gene ontology (GO) annotation, the most highly scored categories of biological process were categorized into cellular process (27%), metabolic process (24%), and biological regulation (8%) in the transcriptome of F. sultana. This study gives first insights into the transcriptome of F. sultana and provides comprehensive genetic resources for the species. We believe that the transcriptome of F. sultana will serve as a useful genomic dataset to accelerate research of functional genomics and will help facilitate whole genome sequencing and annotation. Candidate genes potentially involved in growth, proteolysis, and stress/immunity-response were identified, and are worthy of further investigation.

Neurocognitive Disorders in Heart Failure: Novel Pathophysiological Mechanisms Underpinning Memory Loss and Learning Impairment

Heart failure (HF) is a major public health issue affecting more than 26 million people worldwide. HF is the most common cardiovascular disease in elder population; and it is associated with neurocognitive function decline, which represent underlying brain pathology diminishing learning and memory faculties. Both HF and neurocognitive impairment are associated with recurrent hospitalization episodes and increased mortality rate in older people, but particularly when they occur simultaneously. Overall, the published studies seem to confirm that HF patients display functional impairments relating to attention, memory, concentration, learning, and executive functioning compared with age-matched controls. However, little is known about the molecular mechanisms underpinning neurocognitive decline in HF. The present review round step recent evidence related to the possible molecular mechanism involved in the establishment of neurocognitive disorders during HF. We will make a special focus on cerebral ischemia, neuroinflammation and oxidative stress, Wnt signaling, and mitochondrial DNA alterations as possible mechanisms associated with cognitive decline in HF. Also, we provide an integrative mechanism linking pathophysiological hallmarks of altered cardiorespiratory control and the development of cognitive dysfunction in HF patients. Graphical Abstract Main molecular mechanisms involved in the establishment of cognitive impairment during heart failure. Heart failure is characterized by chronic activation of brain areas responsible for increasing cardiac sympathetic load. In addition, HF patients also show neurocognitive impairment, suggesting that the overall mechanisms that underpin cardiac sympathoexcitation may be related to the development of cognitive disorders in HF. In low cardiac output, HF cerebral infarction due to cardiac mural emboli and cerebral ischemia due to chronic or intermittent cerebral hypoperfusion has been described as a major mechanism related to the development of CI. In addition, while acute norepinephrine (NE) release may be relevant to induce neural plasticity in the hippocampus, chronic or tonic release of NE may exert the opposite effects due to desensitization of the adrenergic signaling pathway due to receptor internalization. Enhanced chemoreflex drive is a major source of sympathoexcitation in HF, and this phenomenon elevates brain ROS levels and induces neuroinflammation through breathing instability. Importantly, both oxidative stress and neuroinflammation can induce mitochondrial dysfunction and vice versa. Then, this ROS inflammatory pathway may propagate within the brain and potentially contribute to the development of cognitive impairment in HF through the activation/inhibition of key molecular pathways involved in neurocognitive decline such as the Wnt signaling pathway.

Exploring the regulatory role of nitric oxide (NO) and the NO-p38MAPK/cGMP pathway in larval settlement of the bryozoan Bugula neritina

The bryozoan Bugula neritina is a cosmopolitan marine fouling species that causes major fouling problems in sub-tropical waters. Settlement of B. neritina larvae can be triggered without an obvious external cue. Here, the negative regulatory role of nitric oxide (NO) during larval settlement of B. neritina was demonstrated to be mediated by cyclic guanosine monophosphate (cGMP). Although the regulatory role of the NO-p38 MAPK signaling axis in larval settlement was not evident, inhibition of nitric oxide synthase (NOS) led to the deactivation of p38 MAPK. Exclusive localization of NO and NO signaling components in sensory-related organs of the larvae is consistent with its signal transduction function in metamorphosis. Overall, this study provides new insights into the regulatory roles of the NO-p38MAPK/cGMP pathway in B. neritina settlement.

Metatranscriptomics profile of the gill microbial community during Bathymodiolus azoricus aquarium acclimatization at atmospheric pressure

Background

The deep-sea mussels Bathymodiolus azoricus (Bivalvia: Mytilidae) are the dominant macrofauna subsisting at the hydrothermal vents site Menez Gwen in the Mid-Atlantic Ridge (MAR). Their adaptive success in such challenging environments is largely due to their gill symbiotic association with chemosynthetic bacteria. We examined the response of vent mussels as they adapt to sea-level environmental conditions, through an assessment of the relative abundance of host-symbiont related RNA transcripts to better understand how the gill microbiome may drive host-symbiont interactions in vent mussels during hypothetical venting inactivity.

Results

The metatranscriptome of B. azoricus was sequenced from gill tissues sampled at different time-points during a five-week acclimatization experiment, using Next-Generation-Sequencing. After Illumina sequencing, a total of 181,985,262 paired-end reads of 150 bp were generated with an average of 16,544,115 read per sample. Metatranscriptome analysis confirmed that experimental acclimatization in aquaria accounted for global gill transcript variation. Additionally, the analysis of 16S and 18S rRNA sequences data allowed for a comprehensive characterization of host-symbiont interactions, which included the gradual loss of gill endosymbionts and signaling pathways, associated with stress responses and energy metabolism, under experimental acclimatization. Dominant active transcripts were assigned to the following KEGG categories: “Ribosome”, “Oxidative phosphorylation” and “Chaperones and folding catalysts” suggesting specific metabolic responses to physiological adaptations in aquarium environment.

Conclusions

Gill metagenomics analyses highlighted microbial diversity shifts and a clear pattern of varying mRNA transcript abundancies and expression during acclimatization to aquarium conditions which indicate change in bacterial community activity. This approach holds potential for the discovery of new host-symbiont associations, evidencing new functional transcripts and a clearer picture of methane metabolism during loss of endosymbionts. Towards the end of acclimatization, we observed trends in three major functional subsystems, as evidenced by an increment of transcripts related to genetic information processes; the decrease of chaperone and folding catalysts and oxidative phosphorylation transcripts; but no change in transcripts of gluconeogenesis and co-factors-vitamins.

A novel assessment of the traction forces upon settlement of two typical marine fouling invertebrates using PDMS micropost arrays

Marine biofouling poses a severe threat to maritime and aquaculture industries. To prevent the attachment of marine biofouling organisms on man-made structures, countless cost and effort was spent annually. In particular, most attention has been paid on the development of efficient and environmentally friendly fouling-resistant coatings, as well as larval settlement mechanism of several major biofouling invertebrates. In this study, polydimethylsiloxane (PDMS) micropost arrays were utilized as the settlement substrata and opposite tractions were identified during early settlement of the barnacle Amphibalanus amphitrite and the bryozoan Bugula neritina The settling A. amphitrite pushed the periphery microposts with an average traction force of 376.2 nN, while settling B. neritina pulled the periphery microposts with an average traction force of 205.9 nN. These micropost displacements are consistent with the body expansion of A. amphitrite during early post-settlement metamorphosis stage and elevation of wall epithelium of B. neritina during early pre-ancestrula stage, respectively. As such, the usage of micropost array may supplement the traditional histological approach to indicate the early settlement stages or even the initiation of larval settlement of marine fouling organisms, and could finally aid in the development of automatic monitoring platform for the real-time analysis on this complex biological process.

HSP90 regulates larval settlement of the bryozoan Bugula neritina through NO pathway

The larvae of many sessile marine invertebrates go through a settlement process, during which the planktonic larvae attach to a substrate and metamorphose into sessile juveniles. Larval attachment and metamorphosis (herein defined as “settlement”) are complex processes mediated by many signaling pathways. Nitric oxide (NO) signaling is one of the pathways that inhibits larval settlement in marine invertebrates across different phyla. NO is synthesized by NO synthase (NOS), which is a client of molecular chaperon heat shock protein 90 (HSP90). In the present study, we provide evidence that NO, a gaseous messenger, regulates larval settlement of B. neritina. By using pharmacological bioassays and western blotting, we demonstrated that NO inhibits larval settlement of B. neritina and that NO signals occur mainly in the sensory organ of swimming larvae. The settlement rate of B. neritina larvae decreased after heat shock treatment. Inhibition of HSP90 induced larval settlement, and attenuated the inhibition of NO donors during larval settlement. In addition, the expression level of both HSP90 and NOS declined upon settlement. These results demonstrate that HSP90 regulates the larval settlement of B. neritina by interacting with the NO pathway.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

Transcriptomic alterations in Sitophilus zeamais in response to allyl isothiocyanate fumigation

To study the fumigation mechanisms of Allyl isothiocyanate (AITC) a promising biorational alternative to present fumigants (phosphine and methyl bromide), and provide theoretical basis for its further development in the control of stored grain pests, this research presents a transcriptome analysis of Sitophilus zeamais fumigated with AITC at the concentration of LC₅₀ (5.69μg/mL) and control over 8h. 21,869,022 and 23,873,110 clean reads in insects fumigated with AITC and control were gained, respectively. The results of RNA-seq were confirmed by qRT-PCR determination of the expression levels of NADH dehydrogenase subunit 6 and Vacuolar ATP synthase subunit B in the insects fumigated with AITC at different concentrations. After enrichment analysis of differentially expressed genes, 117 over-expressed and 271 down-regulated transcripts were gained. Following GO enrichment, these transcripts were classified into 38 GO subgroups (at level 2), and the majority enriched GO terms were "Binding" "Cell process" and "metabolic". KEGG enrichment analysis showed that the majority enriched pathway were "Folding, sorting and degradation", "Transport and catabolism", "Energy metabolism", and "Carbohydrate metabolism". Connected with previous researches on mechanisms of isothiocyanates, cytoskeleton collapse and mitochondria dysfunction are proposed to be significant lethal mechanisms of AITC.

Transcriptome Analysis in Domesticated Species: Challenges and Strategies

Domesticated species occupy a special place in the human world due to their economic and cultural value. In the era of genomic research, domesticated species provide unique advantages for investigation of diseases and complex phenotypes. RNA sequencing, or RNA-seq, has recently emerged as a new approach for studying transcriptional activity of the whole genome, changing the focus from individual genes to gene networks. RNA-seq analysis in domesticated species may complement genome-wide association studies of complex traits with economic importance or direct relevance to biomedical research. However, RNA-seq studies are more challenging in domesticated species than in model organisms. These challenges are at least in part associated with the lack of quality genome assemblies for some domesticated species and the absence of genome assemblies for others. In this review, we discuss strategies for analyzing RNA-seq data, focusing particularly on questions and examples relevant to domesticated species.

De Novo Assembly and Characterization of the Invasive Northern Pacific Seastar Transcriptome

Invasive species are a major threat to global biodiversity but can also serve as valuable model systems to examine important evolutionary processes. While the ecological aspects of invasions have been well documented, the genetic basis of adaptive change during the invasion process has been hampered by a lack of genomic resources for the majority of invasive species. Here we report the first larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. Approximately 117.5 million 100 base-pair (bp) paired-end reads were sequenced from a single RNA-Seq library from a pooled set of full-sibling A. amurensis bipinnaria larvae. We evaluated the efficacy of a pre-assembly error correction pipeline on subsequent de novo assembly. Error correction resulted in small but important improvements to the final assembly in terms of mapping statistics and core eukaryotic genes representation. The error-corrected de novo assembly resulted in 115,654 contigs after redundancy clustering. 41,667 assembled contigs were homologous to sequences from NCBI’s non-redundant protein and UniProt databases. We assigned Gene Ontology, KEGG Orthology, Pfam protein domain terms and predicted protein-coding sequences to > 36,000 contigs. The final transcriptome dataset generated here provides functional information for 18,319 unique proteins, comprising at least 11,355 expressed genes. Furthermore, we identified 9,739 orthologs to P. miniata proteins, evaluated our annotation pipeline and generated a list of 150 candidate genes for responses to several environmental stressors that may be important for adaptation of A. amurensis in the invasive range. Our study has produced a large set of A. amurensis RNA contigs with functional annotations that can serve as a resource for future comparisons to other echinoderm transcriptomes and gene expression studies. Our data can be used to study the genetic basis of adaptive change and other important evolutionary processes during a successful invasion.

Using proteomic analysis to investigate uniconazole-induced phytohormone variation and starch accumulation in duckweed (Landoltia punctata)

BACKGROUND

Duckweed (Landoltia punctata) has the potential to remediate wastewater and accumulate enormous amounts of starch for bioethanol production. Using systematical screening, we determined that the highest biomass and starch percentage of duckweed was obtained after uniconazole application. Uniconazole contributes to starch accumulation of duckweed, but the molecular mechanism is still unclear.

RESULTS

To elucidate the mechanisms of high starch accumulation, in the study, the responses of L. punctata to uniconazole were investigated using a quantitative proteomic approach combined with physiological and biochemical analysis. A total of 3327 proteins were identified. Among these identified proteins, a large number of enzymes involved in endogenous hormone synthetic and starch metabolic pathways were affected. Notably, most of the enzymes involved in abscisic acid (ABA) biosynthesis showed up-regulated expression, which was consistent with the content variation. The increased endogenous ABA may up-regulate expression of ADP-glucose pyrophosphorylase to promote starch biosynthesis. Importantly, the expression levels of several key enzymes in the starch biosynthetic pathway were up-regulated, which supported the enzymatic assay results and may explain why there is increased starch accumulation.

CONCLUSIONS

These generated data linked uniconazole with changes in expression of enzymes involved in hormone biosynthesis and starch metabolic pathways and elucidated the effect of hormones on starch accumulation. Thus, this study not only provided insights into the molecular mechanisms of uniconazole-induced hormone variation and starch accumulation but also highlighted the potential for duckweed to be feedstock for biofuel as well as for sewage treatment.

RNA-Seq analysis and gene discovery of Andrias davidianus using Illumina short read sequencing

The Chinese giant salamander, Andrias davidianus, is an important species in the course of evolution; however, there is insufficient genomic data in public databases for understanding its immunologic mechanisms. High-throughput transcriptome sequencing is necessary to generate an enormous number of transcript sequences from A. davidianus for gene discovery. In this study, we generated more than 40 million reads from samples of spleen and skin tissue using the Illumina paired-end sequencing technology. De novo assembly yielded 87,297 transcripts with a mean length of 734 base pairs (bp). Based on the sequence similarities, searching with known proteins, 38,916 genes were identified. Gene enrichment analysis determined that 981 transcripts were assigned to the immune system. Tissue-specific expression analysis indicated that 443 of transcripts were specifically expressed in the spleen and skin. Among these transcripts, 147 transcripts were found to be involved in immune responses and inflammatory reactions, such as fucolectin, β-defensins and lymphotoxin beta. Eight tissue-specific genes were selected for validation using real time reverse transcription quantitative PCR (qRT-PCR). The results showed that these genes were significantly more expressed in spleen and skin than in other tissues, suggesting that these genes have vital roles in the immune response. This work provides a comprehensive genomic sequence resource for A. davidianus and lays the foundation for future research on the immunologic and disease resistance mechanisms of A. davidianus and other amphibians.

Natural antifouling compounds: Effectiveness in preventing invertebrate settlement and adhesion

Biofouling represents a major economic issue regarding maritime industries and also raise important environmental concern. International legislation is restricting the use of biocidal-based antifouling (AF) coatings, and increasing efforts have been applied in the search for environmentally friendly AF agents. A wide diversity of natural AF compounds has been described for their ability to inhibit the settlement of macrofouling species. However poor information on the specific AF targets was available before the application of different molecular approaches both on invertebrate settlement strategies and bioadhesive characterization and also on the mechanistic effects of natural AF compounds. This review focuses on the relevant information about the main invertebrate macrofouler species settlement and bioadhesive mechanisms, which might help in the understanding of the reported effects, attributed to effective and non-toxic natural AF compounds towards this macrofouling species. It also aims to contribute to the elucidation of promising biotechnological strategies in the development of natural effective environmentally friendly AF paints.

Estrogen alters the profile of the transcriptome in river snail Bellamya aeruginosa

We evaluated the transcriptome dynamics of the freshwater river snail Bellamya aeruginosa exposed to 17β-estradiol (E2) using the Roche/454 GS-FLX platform. In total, 41,869 unigenes, with an average length of 586 bp, representing 36,181 contigs and 5,688 singlets were obtained. Among them, 18.08, 36.85, and 25.47 % matched sequences in the GenBank non-redundant nucleic acid database, non-redundant protein database, and Swiss protein database, respectively. Annotation of the unigenes with gene ontology, and then mapping them to biological pathways, revealed large groups of genes related to growth, development, reproduction, signal transduction, and defense mechanisms. Significant differences were found in gene expression in both liver and testicular tissues between control and E2-exposed organisms. These changes in gene expression will help in understanding the molecular mechanisms of the response to physiological stress in the river snail exposed to estrogen, and will facilitate research into biological processes and underlying physiological adaptations to xenoestrogen exposure in gastropods.

De novo transcriptome characterization of Lilium 'Sorbonne' and key enzymes related to the flavonoid biosynthesis

Lily is an important cut-flower and bulb crop in the commercial market. Here, transcriptome profiling of Lilium 'Sorbonne' was conducted through de novo sequencing based on Illumina platform. This research aims at revealing basic information and data that can be used for applied purposes especially the molecular regulatory information on flower color formation in lily. In total, 36,920,680 short reads which corresponded to 3.32 GB of total nucleotides, were produced through transcriptome sequencing. These reads were assembled into 39,636 Unigenes, of which 30,986 were annotated in Nr, Nt, Swiss-Prot, KEGG, COG, GO databases. Based on the three public protein databases, a total of 32,601 coding sequences were obtained. Meanwhile, 19,242 Unigenes were assigned to 128 KEGG pathways. Those with the greatest representation by unique sequences were for ''metabolic pathways'' (5,406 counts, 28.09 %). Our transcriptome revealed 156 Unigenes that encode key enzymes in the flavonoid biosynthesis pathway including CHS, CHI, F3H, FLS, DFR, etc. MISA software identified 2,762 simple sequence repeats, from which 1,975 primers pairs were designed. Over 2,762 motifs were identified, of which the most frequent was AG/CT (659, 23.86 %), followed by A/T (615, 22.27 %) and CCG/CGG (416, 15.06 %). Based on the results, we believe that the color formation of the Lilium 'Sorbonne' flower was mainly controlled by the flavonoid biosynthesis pathway. Additionally, this research provides initial genetic resources that will be valuable to the lily community for other molecular biology research, and the SSRs will facilitate marker-assisted selection in lily breeding.

Transcriptome analysis elucidates key developmental components of bryozoan lophophore development

The most recent phylogenomic study suggested that Bryozoa (Ectoprocta), Brachiopoda, and Phoronida are monophyletic, implying that the lophophore of bryozoans, phoronids and brachiopods is a synapomorphy. Understanding the molecular mechanisms of the lophophore development of the Lophophorata clade can therefore provide us a new insight into the formation of the diverse morphological traits in metazoans. In the present study, we profiled the transcriptome of the Bryozoan (Ectoproct) Bugula neritina during the swimming larval stage (SW) and the early (4 h) and late (24 h) metamorphic stages using the Illumina HiSeq2000 platform. Various genes that function in development, the immune response and neurogenesis showed differential expression levels during metamorphosis. In situ hybridization of 23 genes that participate in the Wnt, BMP, Notch, and Hedgehog signaling pathways revealed their regulatory roles in the development of the lophophore and the ancestrula digestive tract. Our findings support the hypothesis that developmental precursors of the lophophore and the ancestrula digestive tract are pre-patterned by the differential expression of key developmental genes according to their fate. This study provides a foundation to better understand the developmental divergence and/or convergence among developmental precursors of the lophophore of bryozoans, branchiopods and phoronids.

Rapid development of microsatellite markers for Callosobruchus chinensis using Illumina paired-end sequencing

Background

The adzuki bean weevil, Callosobruchus chinensis L., is one of the most destructive pests of stored legume seeds such as mungbean, cowpea, and adzuki bean, which usually cause considerable loss in the quantity and quality of stored seeds during transportation and storage. However, a lack of genetic information of this pest results in a series of genetic questions remain largely unknown, including population genetic structure, kinship, biotype abundance, and so on. Co-dominant microsatellite markers offer a great resolving power to determine these events. Here, we report rapid microsatellite isolation from C. chinensis via high-throughput sequencing.

Principal Findings

In this study, 94,560,852 quality-filtered and trimmed reads were obtained for the assembly of genome using Illumina paired-end sequencing technology. In total, the genome with total length of 497,124,785 bp, comprising 403,113 high quality contigs was generated with de novo assembly. More than 6800 SSR loci were detected and a suit of 6303 primer pair sequences were designed and 500 of them were randomly selected for validation. Of these, 196 pair of primers, i.e. 39.2%, produced reproducible amplicons that were polymorphic among 8 C. chinensis genotypes collected from different geographical regions. Twenty out of 196 polymorphic SSR markers were used to analyze the genetic diversity of 18 C. chinensis populations. The results showed the twenty SSR loci were highly polymorphic among these populations.

Conclusions

This study presents a first report of genome sequencing and de novo assembly for C. chinensis and demonstrates the feasibility of generating a large scale of sequence information and SSR loci isolation by Illumina paired-end sequencing. Our results provide a valuable resource for C. chinensis research. These novel markers are valuable for future genetic mapping, trait association, genetic structure and kinship among C. chinensis.

Quantitative proteomics study of larval settlement in the Barnacle Balanus amphitrite

Barnacles are major sessile components of the intertidal areas worldwide, and also one of the most dominant fouling organisms in fouling communities. Larval settlement has a crucial ecological effect not only on the distribution of the barnacle population but also intertidal community structures. However, the molecular mechanisms involved in the transition process from the larval to the juvenile stage remain largely unclear. In this study, we carried out comparative proteomic profiles of stage II nauplii, stage VI nauplii, cyprids, and juveniles of the barnacle Balanus amphitrite using label-free quantitative proteomics, followed by the measurement of the gene expression levels of candidate proteins. More than 700 proteins were identified at each stage; 80 were significantly up-regulated in cyprids and 95 in juveniles vs other stages. Specifically, proteins involved in energy and metabolism, the nervous system and signal transduction were significantly up-regulated in cyprids, whereas proteins involved in cytoskeletal remodeling, transcription and translation, cell proliferation and differentiation, and biomineralization were up-regulated in juveniles, consistent with changes associated with larval metamorphosis and tissue remodeling in juveniles. These findings provided molecular evidence for the morphological, physiological and biological changes that occur during the transition process from the larval to the juvenile stages in B. amphitrite.

Transcriptome analysis reveals unique C4-like photosynthesis and oil body formation in an arachidonic acid-rich microalga Myrmecia incisa Reisigl H4301

BACKGROUND

Arachidonic acid (ArA) is important for human health because it is one of the major components of mammalian brain membrane phospholipids. The interest in ArA inspired the search for a new sustainable source, and the green microalga Myrmecia incisa Reisigl H4301 has been found a potential ArA-producer due to a high content of intracellular ArA. To gain more molecular information about metabolism pathways, including the biosynthesis of ArA in the non-model microalga, a transcriptomic analysis was performed.

RESULTS

The 454 pyrosequencing generated 371,740 high-quality reads, which were assembled into 51,908 unique sequences consisting of 22,749 contigs and 29,159 singletons. A total of 11,873 unique sequences were annotated through BLAST analysis, and 3,733 were assigned to Gene Ontology (GO) categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis uncovered a C4-like photosynthesis pathway in M. incisa. The biosynthesis pathways of lipid particularly those of ArA and triacylglycerol (TAG) were analyzed in detail, and TAG was proposed to be accumulated in oil bodies in the cytosol with the help of caleosin or oil globule-associated proteins. In addition, the carotenoid biosynthesis pathways are discussed.

CONCLUSION

This transcriptomic analysis of M. incisa enabled a global understanding of mechanisms involved in photosynthesis, de novo biosynthesis of ArA, metabolism of carotenoids, and accumulation of TAG in M. incisa. These findings provided a molecular basis for the research and possibly economic exploitation of this ArA-rich microalga.

Transcriptome and proteome data reveal candidate genes for pollinator attraction in sexually deceptive orchids

BACKGROUND

Sexually deceptive orchids of the genus Ophrys mimic the mating signals of their pollinator females to attract males as pollinators. This mode of pollination is highly specific and leads to strong reproductive isolation between species. This study aims to identify candidate genes responsible for pollinator attraction and reproductive isolation between three closely related species, O. exaltata, O. sphegodes and O. garganica. Floral traits such as odour, colour and morphology are necessary for successful pollinator attraction. In particular, different odour hydrocarbon profiles have been linked to differences in specific pollinator attraction among these species. Therefore, the identification of genes involved in these traits is important for understanding the molecular basis of pollinator attraction by sexually deceptive orchids.

RESULTS

We have created floral reference transcriptomes and proteomes for these three Ophrys species using a combination of next-generation sequencing (454 and Solexa), Sanger sequencing, and shotgun proteomics (tandem mass spectrometry). In total, 121 917 unique transcripts and 3531 proteins were identified. This represents the first orchid proteome and transcriptome from the orchid subfamily Orchidoideae. Proteome data revealed proteins corresponding to 2644 transcripts and 887 proteins not observed in the transcriptome. Candidate genes for hydrocarbon and anthocyanin biosynthesis were represented by 156 and 61 unique transcripts in 20 and 7 genes classes, respectively. Moreover, transcription factors putatively involved in the regulation of flower odour, colour and morphology were annotated, including Myb, MADS and TCP factors.

CONCLUSION

Our comprehensive data set generated by combining transcriptome and proteome technologies allowed identification of candidate genes for pollinator attraction and reproductive isolation among sexually deceptive orchids. This includes genes for hydrocarbon and anthocyanin biosynthesis and regulation, and the development of floral morphology. These data will serve as an invaluable resource for research in orchid floral biology, enabling studies into the molecular mechanisms of pollinator attraction and speciation.

De novo assembly, gene annotation and marker development using Illumina paired-end transcriptome sequences in celery (Apium graveolens L.)

Background

Celery is an increasing popular vegetable species, but limited transcriptome and genomic data hinder the research to it. In addition, a lack of celery molecular markers limits the process of molecular genetic breeding. High-throughput transcriptome sequencing is an efficient method to generate a large transcriptome sequence dataset for gene discovery, molecular marker development and marker-assisted selection breeding.

Principal Findings

Celery transcriptomes from four tissues were sequenced using Illumina paired-end sequencing technology. De novo assembling was performed to generate a collection of 42,280 unigenes (average length of 502.6 bp) that represent the first transcriptome of the species. 78.43% and 48.93% of the unigenes had significant similarity with proteins in the National Center for Biotechnology Information (NCBI) non-redundant protein database (Nr) and Swiss-Prot database respectively, and 10,473 (24.77%) unigenes were assigned to Clusters of Orthologous Groups (COG). 21,126 (49.97%) unigenes harboring Interpro domains were annotated, in which 15,409 (36.45%) were assigned to Gene Ontology(GO) categories. Additionally, 7,478 unigenes were mapped onto 228 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Large numbers of simple sequence repeats (SSRs) were indentified, and then the rate of successful amplication and polymorphism were investigated among 31 celery accessions.

Conclusions

This study demonstrates the feasibility of generating a large scale of sequence information by Illumina paired-end sequencing and efficient assembling. Our results provide a valuable resource for celery research. The developed molecular markers are the foundation of further genetic linkage analysis and gene localization, and they will be essential to accelerate the process of breeding.

Mini-review: Molecular mechanisms of antifouling compounds

Various antifouling (AF) coatings have been developed to protect submerged surfaces by deterring the settlement of the colonizing stages of fouling organisms. A review of the literature shows that effective AF compounds with specific targets are ones often considered non-toxic. Such compounds act variously on ion channels, quorum sensing systems, neurotransmitters, production/release of adhesive, and specific enzymes that regulate energy production or primary metabolism. In contrast, AF compounds with general targets may or may not act through toxic mechanisms. These compounds affect a variety of biological activities including algal photosynthesis, energy production, stress responses, genotoxic damage, immunosuppressed protein expression, oxidation, neurotransmission, surface chemistry, the formation of biofilms, and adhesive production/release. Among all the targets, adhesive production/release is the most common, possibly due to a more extensive research effort in this area. Overall, the specific molecular targets and the molecular mechanisms of most AF compounds have not been identified. Thus, the information available is insufficient to draw firm conclusions about the types of molecular targets to be used as sensitive biomarkers for future design and screening of compounds with AF potential. In this review, the relevant advantages and disadvantages of the molecular tools available for studying the molecular targets of AF compounds are highlighted briefly and the molecular mechanisms of the AF compounds, which are largely a source of speculation in the literature, are discussed.

Marine proteomics: a critical assessment of an emerging technology

The application of proteomics to marine sciences has increased in recent years because the proteome represents the interface between genotypic and phenotypic variability and, thus, corresponds to the broadest possible biomarker for eco-physiological responses and adaptations. Likewise, proteomics can provide important functional information regarding biosynthetic pathways, as well as insights into mechanism of action, of novel marine natural products. The goal of this review is to (1) explore the application of proteomics methodologies to marine systems, (2) assess the technical approaches that have been used, and (3) evaluate the pros and cons of this proteomic research, with the intent of providing a critical analysis of its future roles in marine sciences. To date, proteomics techniques have been utilized to investigate marine microbe, plant, invertebrate, and vertebrate physiology, developmental biology, seafood safety, susceptibility to disease, and responses to environmental change. However, marine proteomics studies often suffer from poor experimental design, sample processing/optimization difficulties, and data analysis/interpretation issues. Moreover, a major limitation is the lack of available annotated genomes and proteomes for most marine organisms, including several "model species". Even with these challenges in mind, there is no doubt that marine proteomics is a rapidly expanding and powerful integrative molecular research tool from which our knowledge of the marine environment, and the natural products from this resource, will be significantly expanded.

Gel-based and gel-free identification of proteins and phosphopeptides during egg-to-larva transition in polychaete Neanthes arenaceodentata

The polychaete Neanthes arenaceodentata- is cosmopolitan in distribution-, has been used as a laboratory test animal. Life history of this species has several unique features; the female dies after spawning and the male incubates the fertilized eggs through the 21-segmented stage. The larvae leave the tube and commence feeding. Changes in protein abundance and phosphorylation were examined during early development of N. arenaceodentata. A gel-based approach and gel-free enrichment of phosphopeptides coupled with mass spectrometry were used to identify proteins and phosphopeptides in fertilized ova and larval stages. Patterns of proteins and phosphoproteins changed from fertilized ova to larval stages. Twelve proteins occurred in phosphorylated form and nine as stage specific proteins. Cytoskeletal proteins have exhibited differential phosphorylation from ova to larval stages; whereas, other proteins exhibited stage-specific phosphorylation patterns. Ten phosphopeptides were identified that showed phosphorylation sites on serine or threonine residues. Sixty percent of the identified proteins were related to structural reorganization and others with protein synthesis, stress response and attachment. The abundance and distribution of two cytoskeleton proteins were examined further by 2-DE Western blot analysis. This is the first report on changes in protein expression and phosphorylation sites at Thr/Ser in early development of N. arenaceodentata. The 2-DE proteome maps and identified phosphoproteins contributes toward understanding the state of fertilized ova and early larval stages and serves as a basis for further studies on proteomics changes under different developmental conditions in this and other polychaete species.

Butenolide inhibits marine fouling by altering the primary metabolism of three target organisms

Butenolide is a very promising antifouling compound that inhibits ship hull fouling by a variety of marine organisms, but its antifouling mechanism was previously unknown. Here we report the first study of butenolide's molecular targets in three representative fouling organisms. In the barnacle Balanus (=Amphibalanus) amphitrite, butenolide bound to acetyl-CoA acetyltransferase 1 (ACAT1), which is involved in ketone body metabolism. Both the substrate and the product of ACAT1 increased larval settlement under butenolide treatment, suggesting its functional involvement. In the bryozoan Bugula neritina, butenolide bound to very long chain acyl-CoA dehydrogenase (ACADVL), actin, and glutathione S-transferases (GSTs). ACADVL is the first enzyme in the very long chain fatty acid β-oxidation pathway. The inhibition of this primary pathway for energy production in larvae by butenolide was supported by the finding that alternative energy sources (acetoacetate and pyruvate) increased larval attachment under butenolide treatment. In marine bacterium Vibrio sp. UST020129-010, butenolide bound to succinyl-CoA synthetase β subunit (SCSβ) and inhibited bacterial growth. ACAT1, ACADVL, and SCSβ are all involved in primary metabolism for energy production. These findings suggest that butenolide inhibits fouling by influencing the primary metabolism of target organisms.

Involvement of Wnt signaling pathways in the metamorphosis of the bryozoan Bugula neritina

In this study, we analyzed the metamorphosis of the marine bryozoan Bugula neritina. We observed the morphogenesis of the ancestrula. We defined three distinct pre-ancestrula stages based on the anatomy of the developing polypide and the overall morphology of pre-ancestrula. We then used an annotation based enrichment analysis tool to analyze the B. neritina transcriptome and identified over-representation of genes related to Wnt signaling pathways, suggesting its involvement in metamorphosis. Finally, we studied the temporal-spatial gene expression studies of several Wnt pathway genes. We found that one of the Wnt ligand, BnWnt10, was expressed spatially opposite to the Wnt antagonist BnsFRP within the blastemas, which is the presumptive polypide. Down-stream components of the canonical Wnt signaling pathway were exclusively expressed in the blastemas. Bnβcatenin and BnFz5/8 were exclusively expressed in the blastemas throughout the metamorphosis. Based on the genes expression patterns, we propose that BnWnt10 and BnsFRP may relate to the patterning of the polypide, in which the two genes served as positional signals and contributed to the polarization of the blastemas. Another Wnt ligand, BnWnt6, was expressed in the apical part of the pre-ancestrula epidermis. Overall, our findings suggest that the Wnt signaling pathway may be important to the pattern formation of polypide and the development of epidermis.

High through-put sequencing of the Parhyale hawaiensis mRNAs and microRNAs to aid comparative developmental studies

Understanding the genetic and evolutionary basis of animal morphological diversity will require comparative developmental studies that use new model organisms. This necessitates development of tools for the study of genetics and also the generation of sequence information of the organism to be studied. The development of next generation sequencing technology has enabled quick and cost effective generation of sequence information. Parhyale hawaiensis has emerged as a model organism of choice due to the development of advanced molecular tools, thus P. hawaiensis genetic information will help drive functional studies in this organism.

Here we present a transcriptome and miRNA collection generated using next generation sequencing platforms. We generated approximately 1.7 million reads from a P. hawaiensis cDNA library constructed from embryos up to the germ band stage. These reads were assembled into a dataset comprising 163,501 transcripts.

Using the combined annotation of Annot8r and pfam2go, Gene Ontology classifications was assigned to 20,597 transcripts. Annot8r was used to provide KEGG orthology to our transcript dataset. A total of 25,292 KEGG pathway assignments were defined and further confirmed with reciprocal blast against the NCBI nr protein database. This has identified many P. hawaiensis gene orthologs of key conserved signalling pathways involved in development. We also generated small RNA sequences from P. hawaiensis, identifying 55 conserved miRNAs. Sequenced small RNAs that were not annotated by stringent comparison to mirBase were used to search the Daphnia pulex for possible novel miRNAs. Using a conservative approach, we have identified 51 possible miRNA candidates conserved in the Daphnia pulex genome, which could be potential crustacean/arthropod specific miRNAs. Our study presents gene and miRNA discovery in a new model organism that does not have a sequenced genome. The data provided by our work will be valuable for the P. hawaiensis community as well as the wider evolutionary developmental biology community.

Transcriptome sequencing and characterization for the sea cucumber Apostichopus japonicus (Selenka, 1867)

Background

Sea cucumbers are a special group of marine invertebrates. They occupy a taxonomic position that is believed to be important for understanding the origin and evolution of deuterostomes. Some of them such as Apostichopus japonicus represent commercially important aquaculture species in Asian countries. Many efforts have been devoted to increasing the number of expressed sequence tags (ESTs) for A. japonicus, but a comprehensive characterization of its transcriptome remains lacking. Here, we performed the large-scale transcriptome profiling and characterization by pyrosequencing diverse cDNA libraries from A. japonicus.

Results

In total, 1,061,078 reads were obtained by 454 sequencing of eight cDNA libraries representing different developmental stages and adult tissues in A. japonicus. These reads were assembled into 29,666 isotigs, which were further clustered into 21,071 isogroups. Nearly 40% of the isogroups showed significant matches to known proteins based on sequence similarity. Gene ontology (GO) and KEGG pathway analyses recovered diverse biological functions and processes. Candidate genes that were potentially involved in aestivation were identified. Transcriptome comparison with the sea urchin Strongylocentrotus purpuratus revealed similar patterns of GO term representation. In addition, 4,882 putative orthologous genes were identified, of which 202 were not present in the non-echinoderm organisms. More than 700 simple sequence repeats (SSRs) and 54,000 single nucleotide polymorphisms (SNPs) were detected in the A. japonicus transcriptome.

Conclusion

Pyrosequencing was proven to be efficient in rapidly identifying a large set of genes for the sea cucumber A. japonicus. Through the large-scale transcriptome sequencing as well as public EST data integration, we performed a comprehensive characterization of the A. japonicus transcriptome and identified candidate aestivation-related genes. A large number of potential genetic markers were also identified from the A. japonicus transcriptome. This transcriptome resource would lay an important foundation for future genetic or genomic studies on this species.

Proteomics in evolutionary ecology: linking the genotype with the phenotype

The study of the proteome (proteomics), which includes the dynamics of protein expression, regulation, interactions and its function, has played a less prominent role in evolutionary and ecological investigations in comparison with the study of the genome and transcriptome. There are, however, a number of arguments suggesting that this situation should change. First, the proteome is closer to the phenotype than the genome or the transcriptome, and as such may be more directly responsive to natural selection, and thus closely linked to adaptation. Second, there is evidence of a low correlation between protein and transcript expression levels across genes in many different organisms. Finally, there have been some recent important technological improvements in proteomics methods that make them feasible, practical and useful to address a wide range of evolutionary questions even in nonmodel organisms. The different proteomic methods, their limitations and problems when interpreting empirical data are described and discussed. In addition, the proteomic literature pertaining to evolutionary ecology is reviewed with examples, and potential applications of proteomics in a variety of evolutionary contexts are outlined. New proteomic research trends such as the study of posttranslational modifications and protein-protein interactions, as well as the combined use of the different -omics approaches, are discussed in relation to the development of a more functional and integrated perspective, needed for achieving a more comprehensive knowledge of evolutionary change.

Toward an understanding of the molecular mechanisms of barnacle larval settlement: a comparative transcriptomic approach

BACKGROUND

The barnacle Balanus amphitrite is a globally distributed biofouler and a model species in intertidal ecology and larval settlement studies. However, a lack of genomic information has hindered the comprehensive elucidation of the molecular mechanisms coordinating its larval settlement. The pyrosequencing-based transcriptomic approach is thought to be useful to identify key molecular changes during larval settlement.

METHODOLOGY AND PRINCIPAL FINDINGS

Using 454 pyrosequencing, we collected totally 630,845 reads including 215,308 from the larval stages and 415,537 from the adults; 23,451 contigs were generated while 77,785 remained as singletons. We annotated 31,720 of the 92,322 predicted open reading frames, which matched hits in the NCBI NR database, and identified 7,954 putative genes that were differentially expressed between the larval and adult stages. Of these, several genes were further characterized with quantitative real-time PCR and in situ hybridization, revealing some key findings: 1) vitellogenin was uniquely expressed in late nauplius stage, suggesting it may be an energy source for the subsequent non-feeding cyprid stage; 2) the locations of mannose receptors suggested they may be involved in the sensory system of cyprids; 3) 20 kDa-cement protein homologues were expressed in the cyprid cement gland and probably function during attachment; and 4) receptor tyrosine kinases were expressed higher in cyprid stage and may be involved in signal perception during larval settlement.

CONCLUSIONS

Our results provide not only the basis of several new hypotheses about gene functions during larval settlement, but also the availability of this large transcriptome dataset in B. amphitrite for further exploration of larval settlement and developmental pathways in this important marine species.

Phosphoproteome analysis during larval development and metamorphosis in the spionid polychaete Pseudopolydora vexillosa

Background

The metamorphosis of the spionid polychaete Pseudopolydora vexillosa includes spontaneous settlement onto soft-bottom habitats and morphogenesis that can be completed in a very short time. A previous study on the total changes to the proteome during the various developmental stages of P. vexillosa suggested that little or no de novo protein synthesis occurs during metamorphosis. In this study, we used multicolor fluorescence detection of proteins in 2-D gels for differential analysis of proteins and phosphoproteins to reveal the dynamics of post-translational modification proteins in this species. A combination of affinity chromatography, 2D-PAGE, and mass spectrometry was used to identify the phosphoproteins in pre-competent larvae, competent larvae, and newly metamorphosed juveniles.

Results

We reproducibly detected 210, 492, and 172 phosphoproteins in pre-competent larvae, competent larvae, and newly metamorphosed juveniles, respectively. The highest percentage of phosphorylation was observed during the competent larval stage. About 64 stage-specific phosphoprotein spots were detected in the competent stage, and 32 phosphoproteins were found to be significantly differentially expressed in the three stages. We identified 38 phosphoproteins, 10 of which were differentially expressed during metamorphosis. These phosphoproteins belonged to six categories of biological processes: (1) development, (2) cell differentiation and integrity, (3) transcription and translation, (4) metabolism, (5) protein-protein interaction and proteolysis, and (6) receptors and enzymes.

Conclusion

This is the first study to report changes in phosphoprotein expression patterns during the metamorphosis of the marine polychaete P. vexillosa. The higher degree of phosphorylation during the process of attaining competence to settle and metamorphose may be due to fast morphological transitions regulated by various mechanisms. Our data are consistent with previous studies showing a high percentage of phosphorylation during competency in the barnacle Balanus amphitrite and the bryozoan Bugula neritina. The identified phosphoproteins may play an important role during metamorphosis, and further studies on the location and functions of important proteins during metamorphosis are warranted.

Transcriptome sequencing and annotation of the microalgae Dunaliella tertiolecta: pathway description and gene discovery for production of next-generation biofuels

Background

Biodiesel or ethanol derived from lipids or starch produced by microalgae may overcome many of the sustainability challenges previously ascribed to petroleum-based fuels and first generation plant-based biofuels. The paucity of microalgae genome sequences, however, limits gene-based biofuel feedstock optimization studies. Here we describe the sequencing and de novo transcriptome assembly for the non-model microalgae species, Dunaliella tertiolecta, and identify pathways and genes of importance related to biofuel production.

Results

Next generation DNA pyrosequencing technology applied to D. tertiolecta transcripts produced 1,363,336 high quality reads with an average length of 400 bases. Following quality and size trimming, ~ 45% of the high quality reads were assembled into 33,307 isotigs with a 31-fold coverage and 376,482 singletons. Assembled sequences and singletons were subjected to BLAST similarity searches and annotated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology (KO) identifiers. These analyses identified the majority of lipid and starch biosynthesis and catabolism pathways in D. tertiolecta.

Conclusions

The construction of metabolic pathways involved in the biosynthesis and catabolism of fatty acids, triacylglycrols, and starch in D. tertiolecta as well as the assembled transcriptome provide a foundation for the molecular genetics and functional genomics required to direct metabolic engineering efforts that seek to enhance the quantity and character of microalgae-based biofuel feedstock.

Comparing de novo assemblers for 454 transcriptome data

BACKGROUND

Roche 454 pyrosequencing has become a method of choice for generating transcriptome data from non-model organisms. Once the tens to hundreds of thousands of short (250-450 base) reads have been produced, it is important to correctly assemble these to estimate the sequence of all the transcripts. Most transcriptome assembly projects use only one program for assembling 454 pyrosequencing reads, but there is no evidence that the programs used to date are optimal. We have carried out a systematic comparison of five assemblers (CAP3, MIRA, Newbler, SeqMan and CLC) to establish best practices for transcriptome assemblies, using a new dataset from the parasitic nematode Litomosoides sigmodontis.

RESULTS

Although no single assembler performed best on all our criteria, Newbler 2.5 gave longer contigs, better alignments to some reference sequences, and was fast and easy to use. SeqMan assemblies performed best on the criterion of recapitulating known transcripts, and had more novel sequence than the other assemblers, but generated an excess of small, redundant contigs. The remaining assemblers all performed almost as well, with the exception of Newbler 2.3 (the version currently used by most assembly projects), which generated assemblies that had significantly lower total length. As different assemblers use different underlying algorithms to generate contigs, we also explored merging of assemblies and found that the merged datasets not only aligned better to reference sequences than individual assemblies, but were also more consistent in the number and size of contigs.

CONCLUSIONS

Transcriptome assemblies are smaller than genome assemblies and thus should be more computationally tractable, but are often harder because individual contigs can have highly variable read coverage. Comparing single assemblers, Newbler 2.5 performed best on our trial data set, but other assemblers were closely comparable. Combining differently optimal assemblies from different programs however gave a more credible final product, and this strategy is recommended.

Dependency on de novo protein synthesis and proteomic changes during metamorphosis of the marine bryozoan Bugula neritina

Background

Metamorphosis in the bryozoan Bugula neritina (Linne) includes an initial phase of rapid morphological rearrangement followed by a gradual phase of morphogenesis. We hypothesized that the first phase may be independent of de novo synthesis of proteins and, instead, involves post-translational modifications of existing proteins, providing a simple mechanism to quickly initiate metamorphosis. To test our hypothesis, we challenged B. neritina larvae with transcription and translation inhibitors. Furthermore, we employed 2D gel electrophoresis to characterize changes in the phosphoproteome and proteome during early metamorphosis. Differentially expressed proteins were identified by liquid chromatography tandem mass spectrometry and their gene expression patterns were profiled using semi-quantitative real time PCR.

Results

When larvae were incubated with transcription and translation inhibitors, metamorphosis initiated through the first phase but did not complete. We found a significant down-regulation of 60 protein spots and the percentage of phosphoprotein spots decreased from 15% in the larval stage to12% during early metamorphosis. Two proteins--the mitochondrial processing peptidase beta subunit and severin--were abundantly expressed and phosphorylated in the larval stage, but down-regulated during metamorphosis. MPPbeta and severin were also down-regulated on the gene expression level.

Conclusions

The initial morphogenetic changes that led to attachment of B. neritina did not depend on de novo protein synthesis, but the subsequent gradual morphogenesis did. This is the first time that the mitochondrial processing peptidase beta subunit or severin have been shown to be down-regulated on both gene and protein expression levels during the metamorphosis of B. neritina. Future studies employing immunohistochemistry to reveal the expression locality of these two proteins during metamorphosis should provide further evidence of the involvement of these two proteins in the morphogenetic rearrangement of B. neritina.