Dataset Information

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx.

ABSTRACT: BACKGROUND: New methods are needed for genomic-scale analysis of emerging model organisms that exemplify important biological questions but lack fully sequenced genomes. For example, there is an urgent need to understand the potential for corals to adapt to climate change, but few molecular resources are available for studying these processes in reef-building corals. To facilitate genomics studies in corals and other non-model systems, we describe methods for transcriptome sequencing using 454, as well as strategies for assembling a useful catalog of genes from the output. We have applied these methods to sequence the transcriptome of planulae larvae from the coral Acropora millepora. RESULTS: More than 600,000 reads produced in a single 454 sequencing run were assembled into approximately 40,000 contigs with five-fold average sequencing coverage. Based on sequence similarity with known proteins, these analyses identified approximately 11,000 different genes expressed in a range of conditions including thermal stress and settlement induction. Assembled sequences were annotated with gene names, conserved domains, and Gene Ontology terms. Targeted searches using these annotations identified the majority of genes associated with essential metabolic pathways and conserved signaling pathways, as well as novel candidate genes for stress-related processes. Comparisons with the genome of the anemone Nematostella vectensis revealed approximately 8,500 pairs of orthologs and approximately 100 candidate coral-specific genes. More than 30,000 SNPs were detected in the coral sequences, and a subset of these validated by re-sequencing. CONCLUSION: The methods described here for deep sequencing of the transcriptome should be widely applicable to generate catalogs of genes and genetic markers in emerging model organisms. Our data provide the most comprehensive sequence resource currently available for reef-building corals, and include an extensive collection of potential genetic markers for association and population connectivity studies. The characterization of the larval transcriptome for this widely-studied coral will enable research into the biological processes underlying stress responses in corals and evolutionary adaptation to global climate change.

SUBMITTER: Meyer E

PROVIDER: S-EPMC2689275 | biostudies-literature | 2009

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx.

Meyer Eli E Aglyamova Galina V GV Wang Shi S Buchanan-Carter Jade J Abrego David D Colbourne John K JK Willis Bette L BL Matz Mikhail V MV

BMC genomics 20090512

<h4>Background</h4>New methods are needed for genomic-scale analysis of emerging model organisms that exemplify important biological questions but lack fully sequenced genomes. For example, there is an urgent need to understand the potential for corals to adapt to climate change, but few molecular resources are available for studying these processes in reef-building corals. To facilitate genomics studies in corals and other non-model systems, we describe methods for transcriptome sequencing usin ...[more]

PMID: 19435504

Similar Datasets

Project description:BACKGROUND: The Himalayan or Indian Mayapple (Podophyllum hexandrum Royle) produces podophyllotoxin, which is used in the production of semisynthetic anticancer drugs. High throughput transcriptome sequences or genomic sequence data from the Indian Mayapple are essential for further understanding of the podophyllotoxin biosynthetic pathway. RESULTS: 454 pyrosequencing of a P. hexandrum cell culture normalized cDNA library generated 2,667,207 raw reads and 1,503,232 high quality reads, with an average read length of 138 bp. The denovo assembly was performed by Newbler using default and optimized parameters. The optimized parameter generated 40, 380 assembled sequences, comprising 12,940 contigs and 27,440 singlets which resulted in better assembly as compared to default parameters. BLASTX analysis resulted in the annotation of 40,380 contigs/singlet using a cut-off value of ? 1E-03. High similarity to Medicago truncatula using optimized parameters and to Populus trichocarpa using default parameters was noted. The Kyoto encyclopedia of genes and genomes (KEGG) analysis using KEGG Automatic Annotation Server (KAAS) combined with domain analysis of the assembled transcripts revealed putative members of secondary metabolism pathways that may be involved in podophyllotoxin biosynthesis. A proposed schematic pathway for phenylpropanoids and podophyllotoxin biosynthesis was generated. Expression profiling was carried out based on fragments per kilobase of exon per million fragments (FPKM). 1036 simple sequence repeats were predicted in the P. hexandrum sequences. Sixty-nine transcripts were mapped to 99 mature and precursor microRNAs from the plant microRNA database. Around 961 transcripts containing transcription factor domains were noted. High performance liquid chromatography analysis showed the peak accumulation of podophyllotoxin in 12-day cell suspension cultures. A comparative qRT-PCR analysis of phenylpropanoid pathway genes identified in the present data was performed to analyze their expression patterns in 12-day cell culture, callus and rhizome. CONCLUSIONS: The present data will help the identification of the potential genes and transcription factors involved in podophyllotoxin biosynthesis in P. hexandrum. The assembled transcripts could serve as potential candidates for marker discovery and conservation, which should form the foundations for future endeavors.

Project description: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.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.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.

Dataset Information

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx.

Publications

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets