Characterization of the Zoarces viviparus liver transcriptome using massively parallel pyrosequencing

Physiological analysis and transcriptome sequencing of a delayed-green leaf mutant 'Duojiao' of ornamental crabapple (Malus sp.)

Malus spectabilis 'Duojiao' is a spontaneous delayed-green leaf color mutant of M. spectabilis 'Riversii' and has chloroplasts with irregularly arranged vesicles and indistinct stromal lamellae. The yellow leaves of mutant have less chlorophyll (Chl), carotenoids, and flavonoids. Measurement of photosynthetic gas exchange indicated that the mutant has lower photosynthetic activity than 'Riversii' plants. Transcriptome sequencing with the Illumina platform was used to characterize differences in gene expression between the leaves of plants with yellow and green colors and elucidate the molecular mechanisms responsible for variation in leaf color in ornamental crabapple. In the comparison group of mutant yellow leaves and the maternal green leaves, 1848 differentially significant expressed genes (DEGs) were annotated by transcriptomic analysis. Many DEGs and transcription factors were identified related to chloroplast development, Chl synthesis and degradation, photosynthesis, carotenoid biosynthesis, flavonoid biosynthesis and other pathways related to plant leaf color formation. Among these, the Chl biosynthesis-related coproporphyrinogen gene, oxidative decarboxylase gene, and Chl a oxygenase gene were down-regulated, indicating that Chl biosynthesis was blocked. GLK1, which regulates chloroplast development, was down-regulated in yellow leaves. Parallel experiments showed that the content of the Chl synthesis precursors, protoporphyrinogen IX, chlorophyllide a, and chlorophyllide b and the activity of chlorophyllogen III oxidase and chlorophyllide a oxygenase in the yellow leaves of 'Duojiao' were lower than those in the green leaves of 'Riversii'. Thus, leaf color formation is greatly affected by Chl metabolism and chloroplast development. The reliability of the RNA-sequencing data was confirmed by quantitative real-time PCR analysis with 12 selected DEGs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01248-7.

Antifreeze protein dispersion in eelpouts and related fishes reveals migration and climate alteration within the last 20 Ma

Antifreeze proteins inhibit ice growth and are crucial for the survival of supercooled fish living in icy seawater. Of the four antifreeze protein types found in fishes, the globular type III from eelpouts is the one restricted to a single infraorder (Zoarcales), which is the only clade know to have antifreeze protein-producing species at both poles. Our analysis of over 60 unique antifreeze protein gene sequences from several Zoarcales species indicates this gene family arose around 18 Ma ago, in the Northern Hemisphere, supporting recent data suggesting that the Arctic Seas were ice-laden earlier than originally thought. The Antarctic was subject to widespread glaciation over 30 Ma and the Notothenioid fishes that produce an unrelated antifreeze glycoprotein extensively exploited the adjoining seas. We show that species from one Zoarcales family only encroached on this niche in the last few Ma, entering an environment already dominated by ice-resistant fishes, long after the onset of glaciation. As eelpouts are one of the dominant benthic fish groups of the deep ocean, they likely migrated from the north to Antarctica via the cold depths, losing all but the fully active isoform gene along the way. In contrast, northern species have retained both the fully active (QAE) and partially active (SP) isoforms for at least 15 Ma, which suggests that the combination of isoforms is functionally advantageous.

De novo assembly and characterization of the liver transcriptome of Mugil incilis (lisa) using next generation sequencing

Mugil incilis (lisa) is an important commercial fish species in many countries, living along the coasts of the western Atlantic Ocean. It has been used as a model organism for environmental monitoring and ecotoxicological investigations. Nevertheless, available genomic and transcriptomic information for this organism is extremely deficient. The aim of this study was to characterize M. incilis hepatic transcriptome using Illumina paired-end sequencing. A total of 32,082,124 RNA-Seq read pairs were generated utilizing the HiSeq platform and subsequently cleaned and assembled into 93,912 contigs (N50 = 2,019 bp). The analysis of species distribution revealed that M. incilis contigs had the highest number of hits to Stegastes partitus (13.4%). Using a sequence similarity search against the public databases GO and KEGG, a total of 7,301 and 16,967 contigs were annotated, respectively. KEGG database showed genes related to environmental information, metabolism and organismal system pathways were highly annotated. Complete or partial coding DNA sequences for several candidate genes associated with stress responses/detoxification of xenobiotics, as well as housekeeping genes, were employed to design primers that were successfully tested and validated by RT-qPCR. This study presents the first transcriptome resources for Mugil incilis and provides basic information for the development of genomic tools, such as the identification of RNA markers, useful to analyze environmental impacts on this fish Caribbean species.

Biomarker responses in eelpouts from four coastal areas in Sweden, Denmark and Germany

To increase our understanding of possible chemical impacts on coastal fish populations in the Baltic Sea, Kattegat and Skagerrak, the viviparous eelpout (Zoarces viviparus) was used as sentinel species in two major sampling campaigns (spring and autumn) in 16 different coastal sites. Condition factor (CF), liver somatic index (LSI), gonad somatic index (GSI) were measured and the activity of the hepatic enzymes ethoxyresorufin-O-deethylase (EROD), glutathione reductase GR), glutathione S-transferase (GST), catalase (CAT) and muscular activity of acetylcholinesterase (AChE) were assessed. PAH metabolites in bile were also analyzed. The most notable finding in the data set was the low EROD activity in eelpouts collected at the relatively polluted region in Germany compared to the other regions, which could be due to an inhibition of the CYP1A-system or to adaptation to chronic exposure of pollutants in this area. Additionally, low AChE activity was noted in the German region in the autumn campaign and low AChE activity detected in the Danish region in the spring campaign. These differences suggest possible season-specific differences in the use and release of AChE-inhibiting chemicals in the Danish and German regions. Clustering of biomarkers on site level indicated a relationship between CF and GSI and suggested that sites with a high CF contained eelpout that put a larger effort into their larvae development. Clustering of the oxidative stress markers GR, GST and CAT on the individual level reflected a possible coordinated regulation of these enzymes. Overall, the results support the importance of taking into account general regional differences and seasonal variation in biomarker activity when monitoring and assessing the effects of pollution. Despite the expected seasonal variation for most of the measured endpoint, several markers (GSI, EROD and CF) vary similarly between all selected sites in both spring and autumn. This suggests that the differences between sites for these endpoints are independent of season.

The transcriptome of metamorphosing flatfish

BACKGROUND

Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. Metamorphosis in vertebrates is driven by thyroid hormones (THs), but how they orchestrate the cellular, morphological and functional modifications associated with maturation to juvenile/adult states in flatfish is an enigma. Since THs act via thyroid receptors that are ligand activated transcription factors, we hypothesized that the maturation of tissues during metamorphosis should be preceded by significant modifications in the transcriptome. Targeting the unique metamorphosis of flatfish and taking advantage of the large size of Atlantic halibut (Hippoglossus hippoglossus) larvae, we determined the molecular basis of TH action using RNA sequencing.

RESULTS

De novo assembly of sequences for larval head, skin and gastrointestinal tract (GI-tract) yielded 90,676, 65,530 and 38,426 contigs, respectively. More than 57 % of the assembled sequences were successfully annotated using a multi-step Blast approach. A unique set of biological processes and candidate genes were identified specifically associated with changes in morphology and function of the head, skin and GI-tract. Transcriptome dynamics during metamorphosis were mapped with SOLiD sequencing of whole larvae and revealed greater than 8,000 differentially expressed (DE) genes significantly (p < 0.05) up- or down-regulated in comparison with the juvenile stage. Candidate transcripts quantified by SOLiD and qPCR analysis were significantly (r = 0.843; p < 0.05) correlated. The majority (98 %) of DE genes during metamorphosis were not TH-responsive. TH-responsive transcripts clustered into 6 groups based on their expression pattern during metamorphosis and the majority of the 145 DE TH-responsive genes were down-regulated.

CONCLUSIONS

A transcriptome resource has been generated for metamorphosing Atlantic halibut and over 8,000 DE transcripts per stage were identified. Unique sets of biological processes and candidate genes were associated with changes in the head, skin and GI-tract during metamorphosis. A small proportion of DE transcripts were TH-responsive, suggesting that they trigger gene networks, signalling cascades and transcription factors, leading to the overt changes in tissue occurring during metamorphosis.

Sequencing and characterization of a multi-organ Arctic charr transcriptome: A toolbox for investigating polymorphism and seasonal life in a high Arctic fish

The Arctic charr (Salvelinus alpinus L.) inhabits fresh water ecosystems of the high North. The species has developed a strong phenotypic plasticity and variability in life history characteristics which has made this species an attractive model for investigations on phenotype plasticity, morph formation and ecological speciation. Further, the extreme seasonal variations in environmental conditions (e.g. food availability) in the high North induce seasonal changes in phenotype, which require precise timing mechanisms and physiological preparations. Individual gating of life-history strategies (e.g. formation of resident and sea-migrating morphs) and transitions (e.g. maturation) depends on conditional traits (size/energy status) at specific assessment time windows, and complex neuroendocrine regulation, which so far is poorly understood. In the absence of a reference genome, and in order to facilitate the investigation of the complex biological mechanisms of this unique fish model, the present study reveals a reference transcriptome for the Arctic charr. Using Roche 454 GS FLX+, we targeted various organs being either at the crossroads of many key pathways (neuroendocrine, metabolic, behavioral), of different ontological origins or displaying complementary physiological functions. The assemblage yielded 34,690 contigs greater than 1000bp with an average length (1690bp) and annotation rate (52%) within the range, or even higher, than what has been previously obtained with other teleost de novo transcriptomes. We dramatically improve the publically available transcript data on this species that may indeed be useful for various disciplines, from basic research to applied aspects related to conservation issues and aquaculture.

A gene to organism approach--assessing the impact of environmental pollution in eelpout (Zoarces viviparus) females and larvae

A broad biomarker approach was applied to study the effects of marine pollution along the Swedish west coast using the teleost eelpout (Zoarces viviparus) as the sentinel species. Measurements were performed on different biological levels, from the molecular to the organismal, including measurements of messenger RNA (mRNA), proteins, cellular and tissue changes, and reproductive success. Results revealed that eelpout captured in Stenungsund had significantly higher hepatic ethoxyresorufin O-deethylase activity, high levels of both cytochrome P4501A and diablo homolog mRNA, and high prevalence of dead larvae and nuclear damage in erythrocytes. Eelpout collected in Göteborg harbor displayed extensive macrovesicular steatosis, whereby the majority of hepatocytes were affected throughout the liver, which could indicate an effect on lipid metabolism. Results also indicate that eelpouts collected at polluted sites might have an affected immune system, with lower mRNA expression of genes involved in the innate immune system and a higher number of lymphocytes. Biomarker assessment also was performed on livers dissected from unborn eelpout larvae collected from the ovary of the females. No significant differences were noted, which might indicate that the larvae to some extent are protected from effects of environmental pollutants. In conclusion, usage of the selected set of biological markers, covering responses from gene to organism, has demonstrated site-specific biomarker patterns that provided a broad and comprehensive picture of the impact of environmental stressors.

Gene set by de novo assembly of Perilla species and expression profiling between P. frutescens (L.) var. frutescens and var. crispa

Perilla frutescens (L.) Britt. is a self-pollinating annual species and is widely cultivated in China, Korea and Japan as an economic crop and a source of medicine and spices. In this study, we sequenced one cultivar variety (PF98095) of P. frutescens (L.) var. frutescens Britt., which was assembled as reference and other three varieties (PF11109, weedy of var. frutescens, PF06336 and PF06353, cultivars of varieties crispa) in order to carry out comparative expression profiling within cultivar and weedy in varieties frutescens and between varieties frutescens and varieties crispa of cultivar type in P. frutescens. Assembly of PF98095, annotation mapping, DEG (differentially expressed gene) profiling, and comparative analysis were performed. We found that more than 65% of the reads were mapped to the reference of P. frutescens gene set. Moreover, we detected 22,962 DEGs in the weedy variety compared to the cultivar, and also, 22,138 and 23,845 DEGs were identified in two cultivars according to the reference, respectively. The DEGs and functional classification were developed to analyze the differences between weedy and cultivar and between varieties frutescens and varieties crispa of Perilla. Furthermore, candidate genes for the different color and seed size of Perilla were identified that could be further investigated in future study. The herein results may play a significant role, and contribute in functional transcriptome studies of Perilla.

Transcript profiling reveals auxin and cytokinin signaling pathways and transcription regulation during in vitro organogenesis of Ramie (Boehmeria nivea L. Gaud)

In vitro organogenesis, one of the most common pathways leading to in vitro plant regeneration, is widely used in biotechnology and the fundamental study of plant biology. Although previous studies have constructed a complex regulatory network model for Arabidopsis in vitro organogenesis, no related study has been reported in ramie. To generate more complete observations of transcriptome content and dynamics during ramie in vitro organogenesis, we constructed a reference transcriptome library and ten digital gene expression (DGE) libraries for illumina sequencing. Approximately 111.34 million clean reads were obtained for transcriptome and the DGE libraries generated between 13.5 and 18.8 million clean reads. De novo assembly produced 43,222 unigenes and a total of 5,760 differentially expressed genes (DEGs) were filtered. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database, 26 auxin related and 11 cytokinin related DEGs were selected for qRT-PCR validation of two ramie cultivars, which had high (Huazhu No. 5) or extremely low (Dazhuhuangbaima) shoot regeneration abilities. The results revealed differing regulation patterns of auxin and cytokinin in different genotypes. Here we report the first genome-wide gene expression profiling of in vitro organogenesis in ramie and provide an overview of transcription and phytohormone regulation during the process. Furthermore, the auxin and cytokinin related genes have distinct expression patterns in two ramie cultivars with high or extremely low shoot regeneration ability, which has given us a better understanding of the in vitro organogenesis mechanism. This result will provide a foundation for future phytohormone research and lead to improvements of the ramie regeneration system.

Generation and characterization of the sea bass Dicentrarchus labrax brain and liver transcriptomes

The sea bass Dicentrarchus labrax is the center of interest of an increasing number of basic or applied research investigations, even though few genomic or transcriptomic data is available. Current public data only represent a very partial view of its transcriptome. To fill this need, we characterized brain and liver transcriptomes in a generalist manner that would benefit the entire scientific community. We also tackled some bioinformatics questions, related to the effect of RNA fragment size on the assembly quality. Using Illumina RNA-seq, we sequenced organ pools from both wild and farmed Atlantic and Mediterranean fishes. We built two distinct cDNA libraries per organ that only differed by the length of the selected mRNA fragments. Efficiency of assemblies performed on either or both fragments size differed depending on the organ, but remained very close reflecting the quality of the technical replication. We generated more than 19,538Mbp of data. Over 193million reads were assembled into 35,073 contigs (average length=2374bp; N50=3257). 59% contigs were annotated with SwissProt, which corresponded to 12,517 unique genes. We compared the Gene Ontology (GO) contig distribution between the sea bass and the tilapia. We also looked for brain and liver GO specific signatures as well as KEGG pathway coverage. 23,050 putative micro-satellites and 134,890 putative SNPs were identified. Our sampling strategy and assembly pipeline provided a reliable and broad reference transcriptome for the sea bass. It constitutes an indisputable quantitative and qualitative improvement of the public data, as it provides 5 times more base pairs with fewer and longer contigs. Both organs present unique signatures consistent with their specific physiological functions. The discrepancy in fragment size effect on assembly quality between organs lies in their difference in complexity and thus does not allow prescribing any general strategy. This information on two key organs will facilitate further functional approaches.

Transcriptome profiling of fruit development and maturation in Chinese white pear (Pyrus bretschneideri Rehd)

BACKGROUND

Pear (Pyrus spp) is an important fruit species worldwide; however, its genetics and genomic information is limited. Combining the Solexa/Illumina RNA-seq high-throughput sequencing approach (RNA-seq) with Digital Gene Expression (DGE) analysis would be a powerful tool for transcriptomic study. This paper reports the transcriptome profiling analysis of Chinese white pear (P. bretschneideri) using RNA-seq and DGE to better understand the molecular mechanisms in fruit development and maturation of Chinese white pear.

RESULTS

De novo transcriptome assembly and gene expression analysis of Chinese white pear were performed in an unprecedented depth (5.47 gigabase pairs) using high-throughput Illumina RNA-seq combined with a tag-based Digital Gene Expression (DGE) system. Approximately, 60.77 million reads were sequenced, trimmed, and assembled into 90,227 unigenes. These unigenes comprised 17,619 contigs and 72,608 singletons with an average length of 508 bp and had an N50 of 635 bp. Sequence similarity analyses against six public databases (Uniprot, NR, and COGs at NCBI, Pfam, InterPro, and KEGG) found that 61,636 unigenes can be annotated with gene descriptions, conserved protein domains, or gene ontology terms. By BLASTing all 61,636 unigenes in KEGG, a total of 31,215 unigenes were annotated into 121 known metabolic or signaling pathways in which a few primary, intermediate, and secondary metabolic pathways are directly related to pear fruit quality. DGE libraries were constructed for each of the five fruit developmental stages. Variations in gene expression among all developmental stages of pear fruit were significantly different in a large amount of unigenes.

CONCLUSION

Extensive transcriptome and DGE profiling data at five fruit developmental stages of Chinese white pear have been obtained from a deep sequencing, which provides comprehensive gene expression information at the transcriptional level. This could facilitate understanding of the molecular mechanisms in fruit development and maturation. Such a database can also be used as a public information platform for research on molecular biology and functional genomics in pear and other related species.

Model of gene expression in extreme cold - reference transcriptome for the high-Antarctic cryopelagic notothenioid fish Pagothenia borchgrevinki

BACKGROUND

Among the cold-adapted Antarctic notothenioid fishes, the high-latitude bald notothen Pagothenia borchgrevinki is particularly notable as the sole cryopelagic species, exploiting the coldest and iciest waters of the Southern Ocean. Because P. borchgrevinki is a frequent model for investigating notothenioid cold-adaptation and specialization, it is imperative that "omic" tools be developed for this species. In the absence of a sequenced genome, a well annotated reference transcriptome of the bald notothen will serve as a model of gene expression in the coldest and harshest of all polar marine environments, useful for future comparative studies of cold adaptation and thermal responses in polar teleosts and ectotherms.

RESULTS

We sequenced and annotated a reference transcriptome for P. borchgrevinki, with added attention to capturing the transcriptional responses to acute and chronic heat exposures. We sequenced by Roche 454 a normalized cDNA library constructed from pooled mRNA encompassing multiple tissues taken from environmental, warm acclimating, and acute heat stressed specimens. The resulting reads were assembled into 42,620 contigs, 17,951 of which could be annotated. We utilized this annotated portion of the reference transcriptome to map short Illumina reads sequenced from the gill and liver of environmental specimens, and also compared the gene expression profiles of these two tissue transcriptomes with those from the temperate model fish Danio rerio. From this, we identified a conserved group of 58 GO terms, in which terms related to transcription and its regulation, ubiquitin-protein ligase activity, protein ubiquitination, and protein binding among others are more prevalent in the bald notothen, suggesting the pertinent genes play essential roles in cold temperature functioning.

CONCLUSION

We sequenced multiple tissue transcriptomes from native and heat-exposed experimental specimens of the high Antarctic, cryopelagic notothenioid P. borchgrevinki to construct a reference transcriptome. In a proof of concept, we utilized the annotated reference transcriptome to profile the gene expression patterns of gill and liver, and identified a suite of over and under-represented GO terms when compared to the tropical water zebrafish suggesting these functions may be important for surviving in freezing waters. The transcriptome resource from this study will aid future investigations of cold adaptation and thermal response of polar ectothermic species.

Illumina-based de novo transcriptome sequencing and analysis of Amanita exitialis basidiocarps

Amanita exitialis is a lethal mushroom that was first discovered in Guangdong Province, China. The high content of amanitin in its basidiocarps makes it lethal to humans. To comprehensively characterize the A. exitialis transcriptome and analyze the Amanita toxins as well as their related gene family, transcriptome sequencing of A. exitialis was performed using Illumina HiSeq 2000 technology. A total of 25,563,688 clean reads were collected and assembled into 62,137 cDNA contigs with an average length of 481 bp and N50 length of 788 bp. A total of 27,826 proteins and 39,661 unigenes were identified among the assembled contigs. All of the unigenes were classified into 166 functional categories for understanding the gene functions and regulation pathways. The genes contributing to toxic peptide biosynthesis were analyzed. From this set, eleven gene sequences encoding the toxins or related cyclic peptides were discovered in the transcriptome. Three of these sequences matched the peptide toxins α-amanitin, β-amanitin, and phallacidin, while others matched amanexitide and seven matched unknown peptides. All of the genes encoding peptide toxins were confirmed by polymerase chain reaction (PCR) in A. exitialis, and the phylogenetic relationships among these proprotein sequences were discussed. The gene polymorphism and degeneracy of the toxin encoding sequences were found and analyzed. This study provides the first primary transcriptome of A. exitialis, which provided comprehensive gene expression information on the lethal amanitas at the transcriptional level, and could lay a strong foundation for functional genomics studies in those fungi.

Hepatic transcriptome profiling indicates differential mRNA expression of apoptosis and immune related genes in eelpout (Zoarces viviparus) caught at Göteborg harbor, Sweden

The physiology and reproductive performance of eelpout (Zoarces viviparus) have been monitored along the Swedish coast for more than three decades. In this study, transcriptomic profiling was applied for the first time as an exploratory tool to search for new potential candidate biomarkers and to investigate possible stress responses in fish collected from a chronically polluted area. An oligonucleotide microarray with more than 15,000 sequences was used to assess differentially expressed hepatic mRNA levels in female eelpout collected from the contaminated area at Göteborg harbor compared to fish from a national reference site, Fjällbacka. Genes involved in apoptosis and DNA damage (e.g., SMAC/diablo homolog and DDIT4/DNA-damage-inducible protein transcript 4) had higher mRNA expression levels in eelpout from the harbor compared to the reference site, whereas mRNA expression of genes involved in the innate immune system (e.g., complement components and hepcidin) and protein transport/folding (e.g., signal recognition particle and protein disulfide-isomerase) were expressed at lower levels. Gene Ontology enrichment analysis revealed that genes involved biological processes associated with protein folding, immune responses and complement activation were differentially expressed in the harbor eelpout compared to the reference site. The differential mRNA expression of selected genes involved in apoptosis/DNA damage and in the innate immune system was verified by quantitative PCR, using the same fish in addition to eelpout captured four years later. Thus, our approach has identified new potential biomarkers of pollutant exposure and has generated hypotheses on disturbed physiological processes in eelpout. Despite a higher mRNA expression of genes related to apoptosis (e.g., diablo homolog) in eelpout captured in the harbor there were no significant differences in the number of TUNEL-positive apoptotic cells between sites. The mRNA level of genes involved in apoptosis/DNA damage and the status of the innate immune system in fish species captured in polluted environments should be studied in more detail to lay the groundwork for future biomonitoring studies.

Characterization of the pearl oyster (Pinctada martensii) mantle transcriptome unravels biomineralization genes

Pearl oyster, Pinctada martensii, is a marine bivalve species widely distributed in tropic and subtropic marine coasts. Mantle is the special tissue of P. martensii that secretes biomineralization proteins inducing shell deposition as well as iridescent nacre both in the inner shell and artificial nucleus. The pearl oyster is very efficient for artificial pearl production and is therefore an ideal organism for studies into the processes of biomineralization. However, deficiency of transcriptome information limits the insight into biomineralization mechanisms and pearl formation. In this study, we sequenced and characterized the P. martensii mantle transcriptome using 454 pyrosequencing. A total of 25,723 unique transcripts were assembled from 220,824 quality reads, followed by annotation and Gene Ontology classification analysis. A total of 146 unique transcript segments homologous to 49 reference biomineralization genes were identified, including calcineurin-binding protein, amorphous calcium carbonate binding protein 1, calmodulin, calponin-like protein, carbonic anhydrase 1, glycine-rich shell matrix protein, lysine-rich matrix protein, mantle gene or protein, nacrein, pearlin, PIF, regucalcin, and shematrin. The sequence data enabled the identification of 10,285 potential single nucleotide polymorphism loci and 7,836 putative indels, providing a resource for molecular biomarker, population genetics, and functional genomic studies. A large number of candidate genes for biomineralization were identified, considerably enriching resources for the study of shell formation. These sequence data will notably advance biomineralization and transcriptome study in pearl oyster and other Pinctada species.

De novo sequencing and analysis of the Ulva linza transcriptome to discover putative mechanisms associated with its successful colonization of coastal ecosystems

BACKGROUND

The green algal genus Ulva Linnaeus (Ulvaceae, Ulvales, Chlorophyta) is well known for its wide distribution in marine, freshwater, and brackish environments throughout the world. The Ulva species are also highly tolerant of variations in salinity, temperature, and irradiance and are the main cause of green tides, which can have deleterious ecological effects. However, limited genomic information is currently available in this non-model and ecologically important species. Ulva linza is a species that inhabits bedrock in the mid to low intertidal zone, and it is a major contributor to biofouling. Here, we presented the global characterization of the U. linza transcriptome using the Roche GS FLX Titanium platform, with the aim of uncovering the genomic mechanisms underlying rapid and successful colonization of the coastal ecosystems.

RESULTS

De novo assembly of 382,884 reads generated 13,426 contigs with an average length of 1,000 bases. Contiguous sequences were further assembled into 10,784 isotigs with an average length of 1,515 bases. A total of 304,101 reads were nominally identified by BLAST; 4,368 isotigs were functionally annotated with 13,550 GO terms, and 2,404 isotigs having enzyme commission (EC) numbers were assigned to 262 KEGG pathways. When compared with four other full sequenced green algae, 3,457 unique isotigs were found in U. linza and 18 conserved in land plants. In addition, a specific photoprotective mechanism based on both LhcSR and PsbS proteins and a C4-like carbon-concentrating mechanism were found, which may help U. linza survive stress conditions. At least 19 transporters for essential inorganic nutrients (i.e., nitrogen, phosphorus, and sulphur) were responsible for its ability to take up inorganic nutrients, and at least 25 eukaryotic cytochrome P450s, which is a higher number than that found in other algae, may be related to their strong allelopathy. Multi-origination of the stress related proteins, such as glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase and heat-shock proteins, may also contribute to colonization of U. linza under stress conditions.

CONCLUSIONS

The transcriptome of U. linza uncovers some potential genomic mechanisms that might explain its ability to rapidly and successfully colonize coastal ecosystems, including the land-specific genes; special photoprotective mechanism based on both LhcSR and PsbS; development of C4-like carbon-concentrating mechanisms; muti-origin transporters for essential inorganic nutrients; multiple and complex P450s; and glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase, and heat-shock proteins that are related to stress resistance.

Evolutionary force in confamiliar marine vertebrates of different temperature realms: adaptive trends in zoarcid fish transcriptomes

BACKGROUND

Studies of temperature-induced adaptation on the basis of genomic sequence data were mainly done in extremophiles. Although the general hypothesis of an increased molecular flexibility in the cold is widely accepted, the results of thermal adaptation are still difficult to detect at proteomic down to the genomic sequence level. Approaches towards a more detailed picture emerge with the advent of new sequencing technologies. Only small changes in primary protein structure have been shown to modify kinetic and thermal properties of enzymes, but likewise for interspecies comparisons a high genetic identity is still essential to specify common principles. The present study uses comprehensive transcriptomic sequence information to uncover general patterns of thermal adaptation on the RNA as well as protein primary structure.

RESULTS

By comparing orthologous sequences of two closely related zoarcid fish inhabiting different latitudinal zones (Antarctica: Pachycara brachycephalum, temperate zone: Zoarces viviparus) we were able to detect significant differences in the codon usage. In the cold-adapted species a lower GC content in the wobble position prevailed for preserved amino acids. We were able to estimate 40-60% coverage of the functions represented within the two compared zoarcid cDNA-libraries on the basis of a reference genome of the phylogenetically closely related fish Gasterosteus aculeatus. A distinct pattern of amino acid substitutions could be identified for the non-synonymous codon exchanges, with a remarkable surplus of serine and reduction of glutamic acid and asparagine for the Antarctic species.

CONCLUSION

Based on the differences between orthologous sequences from confamiliar species, distinguished mainly by the temperature regimes of their habitats, we hypothesize that temperature leaves a signature on the composition of biological macromolecules (RNA, proteins) with implications for the transcription and translation level. As the observed pattern of amino acid substitutions only partly support the flexibility hypothesis further evolutionary forces may be effective at the global transcriptome level.

Carbonyl reductase mRNA abundance and enzymatic activity as potential biomarkers of oxidative stress in marine fish

Carbonyl reductase (CBR) is an enzyme involved in protection from oxidative stress. In rainbow trout (Oncorhynchus mykiss), the hepatic mRNA abundance of the two isoforms (A and B) is increased after exposure to treated sewage effluents, as well as after exposure with β-naphthoflavone (β-NF) and the pro-oxidant paraquat. In this study, we show that the same chemicals similarly increase the single known hepatic CBR mRNA level and CBR catalytic activity in the coastal living eelpout (Zoarces viviparus). Hepatic CBR mRNA abundance and catalytic activity were also compared between eelpout collected at contaminated and reference sites on the Swedish west coast, but no differences were observed. In conclusion, CBR is a potential biomarker candidate for monitoring the exposure and effects of AhR agonists and/or pro-oxidants in the marine environment, but more research is needed to investigate temporal regulation as well as dose dependency for different chemicals. The mRNA and enzymatic assays presented in this study provide two additional tools for researchers interested in expanding their biomarker battery.

Does ketoprofen or diclofenac pose the lowest risk to fish?

Ketoprofen and diclofenac are non-steroidal anti-inflammatory drugs (NSAIDs) often used for similar indications, and both are frequently found in surface waters. Diclofenac affects organ histology and gene expression in fish at around 1 μg/L. Here, we exposed rainbow trout to ketoprofen (1, 10 and 100 μg/L) to investigate if this alternative causes less risk for pharmacological responses in fish. The bioconcentration factor from water to fish blood plasma was <0.05 (4 for diclofenac based on previous studies). Ketoprofen only reached up to 0.6 ‰ of the human therapeutic plasma concentration, thus the probability of target-related effects was estimated to be fairly low. Accordingly, a comprehensive analysis of hepatic gene expression revealed no consistent responses. In some contrast, trout exposed to undiluted, treated sewage effluents bioconcentrated ketoprofen and other NSAIDs much more efficiently, according to a meta-analysis of recent studies. Neither of the setups is however an ideal representation of the field situation. If a controlled exposure system with a single chemical in pure water is a reasonable representation of the environment, then the use of ketoprofen is likely to pose a lower risk for wild fish than diclofenac, but if bioconcentration factors from effluent-exposed fish are applied, the risks may be more similar.

Transcriptomic characterization of the larval stage in gilthead seabream (Sparus aurata) by 454 pyrosequencing

Gilthead seabream (Sparus aurata) is a teleost belonging to the family Sparidae with a high economical relevance in the Mediterranean countries. Although genomic tools have been developed in this species in order to investigate its physiology at the molecular level and consequently its culture, genomic information on post-embryonic development is still scarce. In this study, we have investigated the transcriptome of a marine teleost during the larval stage (from hatching to 60 days after hatching) by the use of 454 pyrosequencing technology. We obtained a total of 68,289 assembled contigs, representing putative transcripts, belonging to 54,606 different clusters. Comparison against all S. aurata expressed sequenced tags (ESTs) from the NCBI database revealed that up to 34,722 contigs, belonging to about 61% of gene clusters, are sequences previously not described. Contigs were annotated through an iterative Blast pipeline by comparison against databases such as NCBI RefSeq from Danio rerio, SwissProt or NCBI teleost ESTs. Our results indicate that we have enriched the number of annotated sequences for this species by more than 50% compared with previously existing databases for the gilthead seabream. Gene Ontology analysis of these novel sequences revealed that there is a statistically significant number of transcripts with key roles in larval development, differentiation, morphology, and growth. Finally, all information has been made available online through user-friendly interfaces such as GBrowse and a Blast server with a graphical frontend.

Global hepatic gene expression in rainbow trout exposed to sewage effluents: a comparison of different sewage treatment technologies

Effluents from sewage treatment plants contain a mixture of micropollutants with the potential of harming aquatic organisms. Thus, addition of advanced treatment techniques to complement existing conventional methods has been proposed. Some of the advanced techniques could, however, potentially produce additional compounds affecting exposed organisms by unknown modes of action. In the present study the aim was to improve our understanding of how exposure to different sewage effluents affects fish. This was achieved by explorative microarray and quantitative PCR analyses of hepatic gene expression, as well as relative organ sizes of rainbow trout exposed to different sewage effluents (conventionally treated, granular activated carbon, ozonation (5 or 15 mg/L), 5 mg/L ozone plus a moving bed biofilm reactor, or UV-light treatment in combination with hydrogen peroxide). Exposure to the conventionally treated effluent caused a significant increase in liver and heart somatic indexes, an effect removed by all other treatments. Genes connected to xenobiotic metabolism, including cytochrome p450 1A, were differentially expressed in the fish exposed to the conventionally treated effluents, though only effluent treatment with granular activated carbon or ozone at 15 mg/L completely removed this response. The mRNA expression of heat shock protein 70 kDa was induced in all three groups exposed to ozone-treated effluents, suggesting some form of added stress in these fish. The induction of estrogen-responsive genes in the fish exposed to the conventionally treated effluent was effectively reduced by all investigated advanced treatment technologies, although the moving bed biofilm reactor was least efficient. Taken together, granular activated carbon showed the highest potential of reducing responses in fish induced by exposure to sewage effluents.

Applications for next-generation sequencing in fish ecotoxicogenomics

The new technologies for next-generation sequencing (NGS) and global gene expression analyses that are widely used in molecular medicine are increasingly applied to the field of fish biology. This has facilitated new directions to address research areas that could not be previously considered due to the lack of molecular information for ecologically relevant species. Over the past decade, the cost of NGS has decreased significantly, making it possible to use non-model fish species to investigate emerging environmental issues. NGS technologies have permitted researchers to obtain large amounts of raw data in short periods of time. There have also been significant improvements in bioinformatics to assemble the sequences and annotate the genes, thus facilitating the management of these large datasets.The combination of DNA sequencing and bioinformatics has improved our abilities to design custom microarrays and study the genome and transcriptome of a wide variety of organisms. Despite the promising results obtained using these techniques in fish studies, NGS technologies are currently underused in ecotoxicogenomics and few studies have employed these methods. These issues should be addressed in order to exploit the full potential of NGS in ecotoxicological studies and expand our understanding of the biology of non-model organisms.

De novo assembly and characterization of the root transcriptome of Aegilops variabilis during an interaction with the cereal cyst nematode

BACKGROUND

Aegilops variabilis No.1 is highly resistant to cereal cyst nematode (CCN). However, a lack of genomic information has restricted studies on CCN resistance genes in Ae. variabilis and has limited genetic applications in wheat breeding.

RESULTS

Using RNA-Seq technology, we generated a root transcriptome at a sequencing depth of 4.69 gigabases of Ae. variabilis No. 1 from a pooled RNA sample. The sample contained equal amounts of RNA extracted from CCN-infected and untreated control plants at three time-points. Using the Trinity method, nearly 52,081,238 high-quality trimmed reads were assembled into a non-redundant set of 118,064 unigenes with an average length of 500 bp and an N50 of 599 bp. The total assembly was 59.09 Mb of unique transcriptome sequences with average read-depth coverage of 33.25×. In BLAST searches of our database against public databases, 66.46% (78,467) of the unigenes were annotated with gene descriptions, conserved protein domains, or gene ontology terms. Functional categorization further revealed 7,408 individual unigenes and three pathways related to plant stress resistance.

CONCLUSIONS

We conducted high-resolution transcriptome profiling related to root development and the response to CCN infection in Ae. variabilis No.1. This research facilitates further studies on gene discovery and on the molecular mechanisms related to CCN resistance.

Transcriptome sequencing of black grouse (Tetrao tetrix) for immune gene discovery and microsatellite development

The black grouse (Tetrao tetrix) is a galliform bird species that is important for both ecological studies and conservation genetics. Here, we report the sequencing of the spleen transcriptome of black grouse using 454 GS FLX Titanium sequencing. We performed a large-scale gene discovery analysis with a focus on genes that might be related to fitness in this species and also identified a large set of microsatellites. In total, we obtained 182 179 quality-filtered sequencing reads that we assembled into 9035 contigs. Using these contigs and 15 794 length-filtered (greater than 200 bp) singletons, we identified 7762 transcripts that appear to be homologues of chicken genes. A specific BLAST search with an emphasis on immune genes found 308 homologous chicken genes that have immune function, including ten major histocompatibility complex-related genes located on chicken chromosome 16. We also identified 1300 expressed sequence tag microsatellites and were able to design suitable flanking primers for 526 of these. A preliminary test of the polymorphism of the microsatellites found 10 polymorphic microsatellites of the 102 tested. Genomic resources generated in this study should greatly benefit future ecological, evolutionary and conservation genetic studies on this species.

De novo characterization of the Timema cristinae transcriptome facilitates marker discovery and inference of genetic divergence

Adaptation to different ecological environments can promote speciation. Although numerous examples of such 'ecological speciation' now exist, the genomic basis of the process, and the role of gene flow in it, remains less understood. This is, at least in part, because systems that are well characterized in terms of their ecology often lack genomic resources. In this study, we characterize the transcriptome of Timema cristinae stick insects, a system that has been researched intensively in terms of ecological speciation, but for which genomic resources have not been previously developed. Specifically, we obtained >1 million 454 sequencing reads that assembled into 84,937 contigs representing approximately 18,282 unique genes and tens of thousands of potential molecular markers. Second, as an illustration of their utility, we used these genomic resources to assess multilocus genetic divergence within both an ecotype pair and a species pair of Timema stick insects. The results suggest variable levels of genetic divergence and gene flow among taxon pairs and genes and illustrate a first step towards future genomic work in Timema.

Comparison between Normalised and Unnormalised 454-Sequencing Libraries for Small-Scale RNA-Seq Studies

Next-generation sequencing of transcriptomes (RNA-Seq) is being used increasingly in studies of nonmodel organisms. Here, we evaluate the effectiveness of normalising cDNA libraries prior to sequencing in a small-scale study of the zebra finch. We find that assemblies produced from normalised libraries had a larger number of contigs but used fewer reads compared to unnormalised libraries. Considerably more genes were also detected using the contigs produced from normalised cDNA, and microsatellite discovery was up to 73% more efficient in these. There was a positive correlation between the detected expression level of genes in normalised and unnormalised cDNA, and there was no difference in the number of genes identified as being differentially expressed between blood and spleen for the normalised and unnormalised libraries. We conclude that normalised cDNA libraries are preferable for many applications of RNA-Seq and that these can also be used in quantitative gene expression studies.

The polychaete, Perinereis nuntia ESTs and its use to uncover potential biomarker genes for molecular ecotoxicological studies

The polychaete, Perinereis nuntia, has been used as an indicator species to assess the environmental condition of benthic communities in coastal marine environments. Recently, high-throughput sequencing technology has been proven to be a useful method for analyzing expressed sequence tags (ESTs) in non-model species. Thus, we have obtained extensive cDNA information by the pyrosequencing method, to utilize the polychaete species as a test organism for sediment quality monitoring studies. From the total RNA of P. nuntia, cDNA was reversely synthesized and randomly sequenced using a GS-FLX sequencer. In the assembly stage 1, 40,379 transcripts (13,666 contigs and 26,713 singletons) were acquired and showed 47% hitting rate compared with the GenBank non-redundant (NR) amino acid sequence database using BLASTX. After the stage-2 assembly, 21,657 transcripts were identified and showed 28% hitting rate. Finally, we obtained 6 064 unigenes that corresponded to the GenBank NR amino acid sequence database using BLASTX. Of the transcripts obtained in this species, we found a number of stress- and cell defense-related genes (e.g. heat shock protein family, antioxidant-related genes, cytochrome P450 genes) that are potentially useful for sediment monitoring at the molecular level, indicating that the pyrosequencing method is an effective approach to uncover gene families of potential biomarker genes simultaneously, and thus make transcriptomic studies possible. To confirm the usefulness of those potential biomarker genes, we analyzed the comparative profiling of P. nuntia mRNA transcripts between the samples collected from the polychaete aquaculture farm and the southern coast field of South Korea. In this paper, we summarize the expressed cDNA information of P. nuntia and discussed its potential use in environmental genomics and ecotoxicological studies for uncovering the potential molecular mechanisms of environmental stresses and chemical toxicity to the indicator species, P. nuntia in marine sediments.

Genomic toolboxes for conservation biologists

Conservation genetics is expanding its research horizon with a genomic approach, by incorporating the modern techniques of next-generation sequencing (NGS). Application of NGS overcomes many limitations of conservation genetics. First, NGS allows for genome-wide screening of markers, which may lead to a more representative estimation of genetic variation within and between populations. Second, NGS allows for distinction between neutral and non-neutral markers. By screening populations on thousands of single nucleotide polymorphism markers, signals of selection can be found for some markers. Variation in these markers will give insight into functional rather than neutral genetic variation. Third, NGS facilitates the study of gene expression. Conservation genomics will increase our insight in how the environment and genes interact to affect phenotype and fitness. In addition, the NGS approach opens a way to study processes such as inbreeding depression and local adaptation mechanistically. Conservation genetics programs are directed to a fundamental understanding of the processes involved in conservation genetics and should preferably be started in species for which large databases on ecology, demography and genetics are available. Here, we describe and illustrate the connection between the application of NGS technologies and the research questions in conservation. The perspectives of conservation genomics programs are also discussed.

Eelpout (Zoarces viviparus) in marine environmental monitoring

The implementation of the EU Marine Strategy Framework Directive necessitates the development of common criteria and methodological standards for marine environmental monitoring and assessment across Europe. Eelpout (Zoarces viviparus) is proposed as a key indicator organism in the Baltic and North Sea regions. This benthic fish species is widely used in ecotoxicological studies and as a bioindicator of local pollution due to its stationary behavior. Eelpout is included in the environmental monitoring program of several Baltic States, covering both chemical and biological effects measurements, and samples have been archived in environmental specimen banks for >15 years. A method for evaluating the frequency of larval aberrations has been suggested as a standardized assessment tool. The large scientific knowledge-base and considerable experience of long-term chemical and biological effects monitoring and specimen banking, make eelpout a suitable species for the assessment of Good Environmental Status in the Baltic and North Seas.

Physiology and mRNA expression in rainbow trout (Oncorhynchus mykiss) after long-term exposure to the new antifoulant medetomidine

Medetomidine is under evaluation for use as an antifouling agent, and its effects on non-target aquatic organisms are therefore of interest. In this study, rainbow trout was exposed to low (0.5 and 5.0nM) concentrations of medetomidine for up to 54 days. Recently we have reported on effects on paleness and melanophore aggregation of medetomidine in these fish. Here, specific growth rates were investigated together with a broad set of physiological parameters including plasma levels of growth hormone (GH), insulin-like growth factor-I (IGF-I) and leptin, glucose and haemoglobin (Hb), hematocrit (Ht), condition factor, liver and heart somatic indexes (LSI, HSI). Hepatic enzyme activities of CYP1A (EROD activity), glutathione S-transferases (GST) and glutathione reductase (GR) were also measured. Additionally, hepatic mRNA expression was analysed through microarray and quantitative PCR in fish sampled after 31 days of exposure. Medetomidine at both concentrations significantly lowered blood glucose levels and the higher concentration significantly reduced the LSI. The mRNA expression analysis revealed few differentially expressed genes in the liver and the false discovery rate was high. Taken together, the results suggest that medetomidine at investigated concentrations could interfere with carbohydrate metabolism of exposed fish but without any clear consequences for growth.

Diclofenac in fish: blood plasma levels similar to human therapeutic levels affect global hepatic gene expression

Diclofenac is a nonsteroidal anti-inflammatory drug frequently found in the aquatic environment. Previous studies have reported histological changes in the liver, kidney, and gills of fish at concentrations similar to those measured in treated sewage effluents (approximately 1 µg/L). Analyses or predictions of blood plasma levels in fish allow a direct comparison with human therapeutic plasma levels and may therefore be used to indicate a risk for pharmacological effects in fish. To relate internal exposure to a pharmacological interaction, we investigated global hepatic gene expression together with bioconcentration in blood plasma and liver of rainbow trout (Oncorhynchus mykiss) exposed to waterborne diclofenac. At the highest exposure concentration (81.5 µg/L), the fish plasma concentration reached approximately 88% of the human therapeutic levels (C(max) ) after two weeks. Using an oligonucleotide microarray followed by quantitative PCR, we found extensive effects on hepatic gene expression at this concentration, and some genes were found to be regulated down to the lowest exposure concentration tested (1.6 µg/L), corresponding to a plasma concentration approximately 1.5% of the human C(max) . Thus, at concentrations detected in European surface waters, diclofenac can affect the expression of multiple genes in exposed fish. Functional analysis of differentially expressed genes revealed effects on biological processes such as inflammation and the immune response, in agreement with the mode of action of diclofenac in mammals. In contrast to some previously reported results, the bioconcentration factor was found to be stable (4.02 ± 0.75 for blood plasma and 2.54 ± 0.36 for liver) regardless of the water concentration.

Gene discovery and tissue-specific transcriptome analysis in chickpea with massively parallel pyrosequencing and web resource development

Chickpea (Cicer arietinum) is an important food legume crop but lags in the availability of genomic resources. In this study, we have generated about 2 million high-quality sequences of average length of 372 bp using pyrosequencing technology. The optimization of de novo assembly clearly indicated that hybrid assembly of long-read and short-read primary assemblies gave better results. The hybrid assembly generated a set of 34,760 transcripts with an average length of 1,020 bp representing about 4.8% (35.5 Mb) of the total chickpea genome. We identified more than 4,000 simple sequence repeats, which can be developed as functional molecular markers in chickpea. Putative function and Gene Ontology terms were assigned to at least 73.2% and 71.0% of chickpea transcripts, respectively. We have also identified several chickpea transcripts that showed tissue-specific expression and validated the results using real-time polymerase chain reaction analysis. Based on sequence comparison with other species within the plant kingdom, we identified two sets of lineage-specific genes, including those conserved in the Fabaceae family (legume specific) and those lacking significant similarity with any non chickpea species (chickpea specific). Finally, we have developed a Web resource, Chickpea Transcriptome Database, which provides public access to the data and results reported in this study. The strategy for optimization of de novo assembly presented here may further facilitate the transcriptome sequencing and characterization in other organisms. Most importantly, the data and results reported in this study will help to accelerate research in various areas of genomics and implementing breeding programs in chickpea.

Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.

Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance

Background

Until recently, read lengths on the Solexa/Illumina system were too short to reliably assemble transcriptomes without a reference sequence, especially for non-model organisms. However, with read lengths up to 100 nucleotides available in the current version, an assembly without reference genome should be possible. For this study we created an EST data set for the common pond snail Radix balthica by Illumina sequencing of a normalized transcriptome. Performance of three different short read assemblers was compared with respect to: the number of contigs, their length, depth of coverage, their quality in various BLAST searches and the alignment to mitochondrial genes.

Results

A single sequencing run of a normalized RNA pool resulted in 16,923,850 paired end reads with median read length of 61 bases. The assemblies generated by VELVET, OASES, and SeqMan NGEN differed in the total number of contigs, contig length, the number and quality of gene hits obtained by BLAST searches against various databases, and contig performance in the mt genome comparison. While VELVET produced the highest overall number of contigs, a large fraction of these were of small size (< 200bp), and gave redundant hits in BLAST searches and the mt genome alignment. The best overall contig performance resulted from the NGEN assembly. It produced the second largest number of contigs, which on average were comparable to the OASES contigs but gave the highest number of gene hits in two out of four BLAST searches against different reference databases. A subsequent meta-assembly of the four contig sets resulted in larger contigs, less redundancy and a higher number of BLAST hits.

Conclusion

Our results document the first de novo transcriptome assembly of a non-model species using Illumina sequencing data. We show that de novo transcriptome assembly using this approach yields results useful for downstream applications, in particular if a meta-assembly of contig sets is used to increase contig quality. These results highlight the ongoing need for improvements in assembly methodology.

Transcriptomic analysis of the housefly (Musca domestica) larva using massively parallel pyrosequencing

To explore the transcriptome of Musca domestica larvae and to identify unique sequences, we used massively parallel pyrosequencing on the Roche 454-FLX platform to generate a substantial EST dataset of this fly. As a result, we obtained a total of 249,555 ESTs with an average read length of 373 bp. These reads were assembled into 13,206 contigs and 20,556 singletons. Using BlastX searches of the Swissprot and Nr databases, we were able to identify 4,814 contigs and 8,166 singletons as unique sequences. Subsequently, the annotated sequences were subjected to GO analysis and the search results showed a majority of the query sequences were assignable to certain gene ontology terms. In addition, functional classification and pathway assignment were performed by KEGG and 2,164 unique sequences were mapped into 184 KEGG pathways in total. As the first attempt on large-scale RNA sequencing of M. domestica, this general picture of the transcriptome can establish a fundamental resource for further research on functional genomics.

Transcriptome sequencing of lentil based on second-generation technology permits large-scale unigene assembly and SSR marker discovery

Background

Lentil (Lens culinaris Medik.) is a cool-season grain legume which provides a rich source of protein for human consumption. In terms of genomic resources, lentil is relatively underdeveloped, in comparison to other Fabaceae species, with limited available data. There is hence a significant need to enhance such resources in order to identify novel genes and alleles for molecular breeding to increase crop productivity and quality.

Results

Tissue-specific cDNA samples from six distinct lentil genotypes were sequenced using Roche 454 GS-FLX Titanium technology, generating c. 1.38 × 10⁶ expressed sequence tags (ESTs). De novo assembly generated a total of 15,354 contigs and 68,715 singletons. The complete unigene set was sequence-analysed against genome drafts of the model legume species Medicago truncatula and Arabidopsis thaliana to identify 12,639, and 7,476 unique matches, respectively. When compared to the genome of Glycine max, a total of 20,419 unique hits were observed corresponding to c. 31% of the known gene space. A total of 25,592 lentil unigenes were subsequently annoated from GenBank. Simple sequence repeat (SSR)-containing ESTs were identified from consensus sequences and a total of 2,393 primer pairs were designed. A subset of 192 EST-SSR markers was screened for validation across a panel 12 cultivated lentil genotypes and one wild relative species. A total of 166 primer pairs obtained successful amplification, of which 47.5% detected genetic polymorphism.

Conclusions

A substantial collection of ESTs has been developed from sequence analysis of lentil genotypes using second-generation technology, permitting unigene definition across a broad range of functional categories. As well as providing resources for functional genomics studies, the unigene set has permitted significant enhancement of the number of publicly-available molecular genetic markers as tools for improvement of this species.

Sequencing and characterization of the guppy (Poecilia reticulata) transcriptome

Background

Next-generation sequencing is providing researchers with a relatively fast and affordable option for developing genomic resources for organisms that are not among the traditional genetic models. Here we present a de novo assembly of the guppy (Poecilia reticulata) transcriptome using 454 sequence reads, and we evaluate potential uses of this transcriptome, including detection of sex-specific transcripts and deployment as a reference for gene expression analysis in guppies and a related species. Guppies have been model organisms in ecology, evolutionary biology, and animal behaviour for over 100 years. An annotated transcriptome and other genomic tools will facilitate understanding the genetic and molecular bases of adaptation and variation in a vertebrate species with a uniquely well known natural history.

Results

We generated approximately 336 Mbp of mRNA sequence data from male brain, male body, female brain, and female body. The resulting 1,162,670 reads assembled into 54,921 contigs, creating a reference transcriptome for the guppy with an average read depth of 28×. We annotated nearly 40% of this reference transcriptome by searching protein and gene ontology databases. Using this annotated transcriptome database, we identified candidate genes of interest to the guppy research community, putative single nucleotide polymorphisms (SNPs), and male-specific expressed genes. We also showed that our reference transcriptome can be used for RNA-sequencing-based analysis of differential gene expression. We identified transcripts that, in juveniles, are regulated differently in the presence and absence of an important predator, Rivulus hartii, including two genes implicated in stress response. For each sample in the RNA-seq study, >50% of high-quality reads mapped to unique sequences in the reference database with high confidence. In addition, we evaluated the use of the guppy reference transcriptome for gene expression analyses in a congeneric species, the sailfin molly (Poecilia latipinna). Over 40% of reads from the sailfin molly sample aligned to the guppy transcriptome.

Conclusions

We show that next-generation sequencing provided a reliable and broad reference transcriptome. This resource allowed us to identify candidate gene variants, SNPs in coding regions, and sex-specific gene expression, and permitted quantitative analysis of differential gene expression.

Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds

BACKGROUND

Tea is one of the most popular non-alcoholic beverages worldwide. However, the tea plant, Camellia sinensis, is difficult to culture in vitro, to transform, and has a large genome, rendering little genomic information available. Recent advances in large-scale RNA sequencing (RNA-seq) provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes.

RESULTS

Using high-throughput Illumina RNA-seq, the transcriptome from poly (A)+ RNA of C. sinensis was analyzed at an unprecedented depth (2.59 gigabase pairs). Approximate 34.5 million reads were obtained, trimmed, and assembled into 127,094 unigenes, with an average length of 355 bp and an N50 of 506 bp, which consisted of 788 contig clusters and 126,306 singletons. This number of unigenes was 10-fold higher than existing C. sinensis sequences deposited in GenBank (as of August 2010). Sequence similarity analyses against six public databases (Uniprot, NR and COGs at NCBI, Pfam, InterPro and KEGG) found 55,088 unigenes that could be annotated with gene descriptions, conserved protein domains, or gene ontology terms. Some of the unigenes were assigned to putative metabolic pathways. Targeted searches using these annotations identified the majority of genes associated with several primary metabolic pathways and natural product pathways that are important to tea quality, such as flavonoid, theanine and caffeine biosynthesis pathways. Novel candidate genes of these secondary pathways were discovered. Comparisons with four previously prepared cDNA libraries revealed that this transcriptome dataset has both a high degree of consistency with previous EST data and an approximate 20 times increase in coverage. Thirteen unigenes related to theanine and flavonoid synthesis were validated. Their expression patterns in different organs of the tea plant were analyzed by RT-PCR and quantitative real time PCR (qRT-PCR).

CONCLUSIONS

An extensive transcriptome dataset has been obtained from the deep sequencing of tea plant. The coverage of the transcriptome is comprehensive enough to discover all known genes of several major metabolic pathways. This transcriptome dataset can serve as an important public information platform for gene expression, genomics, and functional genomic studies in C. sinensis.

Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds

Background

Tea is one of the most popular non-alcoholic beverages worldwide. However, the tea plant, Camellia sinensis, is difficult to culture in vitro, to transform, and has a large genome, rendering little genomic information available. Recent advances in large-scale RNA sequencing (RNA-seq) provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes.

Results

Using high-throughput Illumina RNA-seq, the transcriptome from poly (A)⁺ RNA of C. sinensis was analyzed at an unprecedented depth (2.59 gigabase pairs). Approximate 34.5 million reads were obtained, trimmed, and assembled into 127,094 unigenes, with an average length of 355 bp and an N50 of 506 bp, which consisted of 788 contig clusters and 126,306 singletons. This number of unigenes was 10-fold higher than existing C. sinensis sequences deposited in GenBank (as of August 2010). Sequence similarity analyses against six public databases (Uniprot, NR and COGs at NCBI, Pfam, InterPro and KEGG) found 55,088 unigenes that could be annotated with gene descriptions, conserved protein domains, or gene ontology terms. Some of the unigenes were assigned to putative metabolic pathways. Targeted searches using these annotations identified the majority of genes associated with several primary metabolic pathways and natural product pathways that are important to tea quality, such as flavonoid, theanine and caffeine biosynthesis pathways. Novel candidate genes of these secondary pathways were discovered. Comparisons with four previously prepared cDNA libraries revealed that this transcriptome dataset has both a high degree of consistency with previous EST data and an approximate 20 times increase in coverage. Thirteen unigenes related to theanine and flavonoid synthesis were validated. Their expression patterns in different organs of the tea plant were analyzed by RT-PCR and quantitative real time PCR (qRT-PCR).

Conclusions

An extensive transcriptome dataset has been obtained from the deep sequencing of tea plant. The coverage of the transcriptome is comprehensive enough to discover all known genes of several major metabolic pathways. This transcriptome dataset can serve as an important public information platform for gene expression, genomics, and functional genomic studies in C. sinensis.

De novo transcriptome characterization and development of genomic tools for Scabiosa columbaria L. using next-generation sequencing techniques

Next-generation sequencing (NGS) technologies are increasingly applied in many organisms, including nonmodel organisms that are important for ecological and conservation purposes. Illumina and 454 sequencing are among the most used NGS technologies and have been shown to produce optimal results at reasonable costs when used together. Here, we describe the combined application of these two NGS technologies to characterize the transcriptome of a plant species of ecological and conservation relevance for which no genomic resource is available, Scabiosa columbaria. We obtained 528,557 reads from a 454 GS-FLX run and a total of 28,993,627 reads from two lanes of an Illumina GAII single run. After read trimming, the de novo assembly of both types of reads produced 109,630 contigs. Both the contigs and the >75 bp remaining singletons were blasted against the Uniprot/Swissprot database, resulting in 29,676 and 10,515 significant hits, respectively. Based on sequence similarity with known gene products, these sequences represent at least 12,516 unique genes, most of which are well covered by contig sequences. In addition, we identified 4320 microsatellite loci, of which 856 had flanking sequences suitable for PCR primer design. We also identified 75,054 putative SNPs. This annotated sequence collection and the relative molecular markers represent a main genomic resource for S. columbaria which should contribute to future research in conservation and population biology studies. Our results demonstrate the utility of NGS technologies as starting point for the development of genomic tools in nonmodel but ecologically important species.

De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum

Background

Because of their phylogenetic position and unique characteristics of their biology and life cycle, ferns represent an important lineage for studying the evolution of land plants. Large and complex genomes in ferns combined with the absence of economically important species have been a barrier to the development of genomic resources. However, high throughput sequencing technologies are now being widely applied to non-model species. We leveraged the Roche 454 GS-FLX Titanium pyrosequencing platform in sequencing the gametophyte transcriptome of bracken fern (Pteridium aquilinum) to develop genomic resources for evolutionary studies.

Results

681,722 quality and adapter trimmed reads totaling 254 Mbp were assembled de novo into 56,256 unique sequences (i.e. unigenes) with a mean length of 547.2 bp and a total assembly size of 30.8 Mbp with an average read-depth coverage of 7.0×. We estimate that 87% of the complete transcriptome has been sequenced and that all transcripts have been tagged. 61.8% of the unigenes had blastx hits in the NCBI nr protein database, representing 22,596 unique best hits. The longest open reading frame in 52.2% of the unigenes had positive domain matches in InterProScan searches. We assigned 46.2% of the unigenes with a GO functional annotation and 16.0% with an enzyme code annotation. Enzyme codes were used to retrieve and color KEGG pathway maps. A comparative genomics approach revealed a substantial proportion of genes expressed in bracken gametophytes to be shared across the genomes of Arabidopsis, Selaginella and Physcomitrella, and identified a substantial number of potentially novel fern genes. By comparing the list of Arabidopsis genes identified by blast with a list of gametophyte-specific Arabidopsis genes taken from the literature, we identified a set of potentially conserved gametophyte specific genes. We screened unigenes for repetitive sequences to identify 548 potentially-amplifiable simple sequence repeat loci and 689 expressed transposable elements.

Conclusions

This study is the first comprehensive transcriptome analysis for a fern and represents an important scientific resource for comparative evolutionary and functional genomics studies in land plants. We demonstrate the utility of high-throughput sequencing of a normalized cDNA library for de novo transcriptome characterization and gene discovery in a non-model plant.

The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus

BACKGROUND

Most evolutionary developmental biology ("evo-devo") studies of emerging model organisms focus on small numbers of candidate genes cloned individually using degenerate PCR. However, newly available sequencing technologies such as 454 pyrosequencing have recently begun to allow for massive gene discovery in animals without sequenced genomes. Within insects, although large volumes of sequence data are available for holometabolous insects, developmental studies of basally branching hemimetabolous insects typically suffer from low rates of gene discovery.

RESULTS

We used 454 pyrosequencing to sequence over 500 million bases of cDNA from the ovaries and embryos of the milkweed bug Oncopeltus fasciatus, which lacks a sequenced genome. This indirectly developing insect occupies an important phylogenetic position, branching basal to Diptera (including fruit flies) and Hymenoptera (including honeybees), and is an experimentally tractable model for short-germ development. 2,087,410 reads from both normalized and non-normalized cDNA assembled into 21,097 sequences (isotigs) and 112,531 singletons. The assembled sequences fell into 16,617 unique gene models, and included predictions of splicing isoforms, which we examined experimentally. Discovery of new genes plateaued after assembly of ~1.5 million reads, suggesting that we have sequenced nearly all transcripts present in the cDNA sampled. Many transcripts have been assembled at close to full length, and there is a net gain of sequence data for over half of the pre-existing O. fasciatus accessions for developmental genes in GenBank. We identified 10,775 unique genes, including members of all major conserved metazoan signaling pathways and genes involved in several major categories of early developmental processes. We also specifically address the effects of cDNA normalization on gene discovery in de novo transcriptome analyses.

CONCLUSIONS

Our sequencing, assembly and annotation framework provide a simple and effective way to achieve high-throughput gene discovery for organisms lacking a sequenced genome. These data will have applications to the study of the evolution of arthropod genes and genetic pathways, and to the wider evolution, development and genomics communities working with emerging model organisms.[The sequence data from this study have been submitted to GenBank under study accession number SRP002610 (http://www.ncbi.nlm.nih.gov/sra?term=SRP002610). Custom scripts generated are available at http://www.extavourlab.com/protocols/index.html. Seven Additional files are available.].

Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus

The accumulation of genetic incompatibilities between isolated populations is thought to lead to the evolution of intrinsic postzygotic isolation. The molecular basis for these mechanisms, however, remains poorly understood. The intertidal copepod Tigriopus californicus provides unique opportunities for addressing mechanistic questions regarding the early stages of speciation; hybrids between highly divergent populations are fertile and viable, but exhibit reduced fitness at the F(2) or later generations. Given the current scarcity of genomic information in taxa at incipient stages of reproductive isolation, we utilize high-throughout 454 pyrosequencing to characterize a substantial fraction of protein-coding regions (the transcriptome) of T. californicus. Our sequencing effort was divided equally between two divergent populations in order to estimate levels of divergence and to reveal patterns of selection across the transcriptome. Assembly of sequences generated over 40,000 putatively unique transcripts (unigenes) for each population, 19,622 of which were orthologous between populations. BLAST searches of public databases determined protein identity and functional features for 15,402 and 12,670 unigenes, respectively. Based on rates of nonsynonymous and synonymous substitutions in 5897 interpopulation orthologs (those >150 bp and with at least 2X coverage), we identified 229 potential targets of positive selection. Many of these genes are predicted to be involved in several metabolic processes, and to function in hydrolase, peptidase and binding activities. The library of T. californicus coding regions, annotated with their predicted functions and level of divergence, will serve as an invaluable resource for elucidating molecular mechanisms underlying the early stages of speciation.

De novo assembled expressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq

Background

De novo assembly of transcript sequences produced by short-read DNA sequencing technologies offers a rapid approach to obtain expressed gene catalogs for non-model organisms. A draft genome sequence will be produced in 2010 for a Eucalyptus tree species (E. grandis) representing the most important hardwood fibre crop in the world. Genome annotation of this valuable woody plant and genetic dissection of its superior growth and productivity will be greatly facilitated by the availability of a comprehensive collection of expressed gene sequences from multiple tissues and organs.

Results

We present an extensive expressed gene catalog for a commercially grown E. grandis × E. urophylla hybrid clone constructed using only Illumina mRNA-Seq technology and de novo assembly. A total of 18,894 transcript-derived contigs, a large proportion of which represent full-length protein coding genes were assembled and annotated. Analysis of assembly quality, length and diversity show that this dataset represent the most comprehensive expressed gene catalog for any Eucalyptus tree. mRNA-Seq analysis furthermore allowed digital expression profiling of all of the assembled transcripts across diverse xylogenic and non-xylogenic tissues, which is invaluable for ascribing putative gene functions.

Conclusions

De novo assembly of Illumina mRNA-Seq reads is an efficient approach for transcriptome sequencing and profiling in Eucalyptus and other non-model organisms. The transcriptome resource (Eucspresso, http://eucspresso.bi.up.ac.za/) generated by this study will be of value for genomic analysis of woody biomass production in Eucalyptus and for comparative genomic analysis of growth and development in woody and herbaceous plants.

The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing

Gene expression divergence is one of the mechanisms thought to be involved in the emergence of incipient species. Next-generation sequencing has become an extremely valuable tool for the study of this process by allowing whole transcriptome sequencing, or RNA-Seq. We have conducted a 454 GS-FLX pyrosequencing experiment to refine our understanding of adaptive divergence between dwarf and normal lake whitefish species (Coregonus clupeaformis spp.). The objectives were to: (i) investigate transcriptomic divergence as measured by liver RNA-Seq; (ii) test the correlation between divergence in expression and sequence polymorphism; and (iii) investigate the extent of allelic imbalance. We also compared the results of RNA-seq with those of a previous microarray study performed on the same fish. Following de novo assembly, results showed that normal whitefish overexpressed more contigs associated with protein synthesis while dwarf fish overexpressed more contigs related to energy metabolism, immunity and DNA replication and repair. Moreover, 63 SNPs showed significant allelic imbalance, and this phenomenon prevailed in the recently diverged dwarf whitefish. Results also showed an absence of correlation between gene expression divergence as measured by RNA-Seq and either polymorphism rate or sequence divergence between normal and dwarf whitefish. This study reiterates an important role for gene expression divergence, and provides evidence for allele-specific expression divergence as well as evolutionary decoupling of regulatory and coding sequences in the adaptive divergence of normal and dwarf whitefish. It also demonstrates how next-generation sequencing can lead to a more comprehensive understanding of transcriptomic divergence in a young species pair.

Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered European eel

Background

Once highly abundant, the European eel (Anguilla anguilla L.; Anguillidae; Teleostei) is considered to be critically endangered and on the verge of extinction, as the stock has declined by 90-99% since the 1980s. Yet, the species is poorly characterized at molecular level with little sequence information available in public databases.

Results

The first European eel transcriptome was obtained by 454 FLX Titanium sequencing of a normalized cDNA library, produced from a pool of 18 glass eels (juveniles) from the French Atlantic coast and two sites in the Mediterranean coast. Over 310,000 reads were assembled in a total of 19,631 transcribed contigs, with an average length of 531 nucleotides. Overall 36% of the contigs were annotated to known protein/nucleotide sequences and 35 putative miRNA identified.

Conclusions

This study represents the first transcriptome analysis for a critically endangered species. EeelBase, a dedicated database of annotated transcriptome sequences of the European eel is freely available at http://compgen.bio.unipd.it/eeelbase. Considering the multiple factors potentially involved in the decline of the European eel, including anthropogenic factors such as pollution and human-introduced diseases, our results will provide a rich source of data to discover and identify new genes, characterize gene expression, as well as for identification of genetic markers scattered across the genome to be used in various applications.

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.