A comprehensive assessment of the transcriptome of cork oak (Quercus suber) through EST sequencing

Background

Cork oak (Quercus suber) is one of the rare trees with the ability to produce cork, a material widely used to make wine bottle stoppers, flooring and insulation materials, among many other uses. The molecular mechanisms of cork formation are still poorly understood, in great part due to the difficulty in studying a species with a long life-cycle and for which there is scarce molecular/genomic information. Cork oak forests are of great ecological importance and represent a major economic and social resource in Southern Europe and Northern Africa. However, global warming is threatening the cork oak forests by imposing thermal, hydric and many types of novel biotic stresses. Despite the economic and social value of the Q. suber species, few genomic resources have been developed, useful for biotechnological applications and improved forest management.

Results

We generated in excess of 7 million sequence reads, by pyrosequencing 21 normalized cDNA libraries derived from multiple Q. suber tissues and organs, developmental stages and physiological conditions. We deployed a stringent sequence processing and assembly pipeline that resulted in the identification of ~159,000 unigenes. These were annotated according to their similarity to known plant genes, to known Interpro domains, GO classes and E.C. numbers. The phylogenetic extent of this ESTs set was investigated, and we found that cork oak revealed a significant new gene space that is not covered by other model species or EST sequencing projects. The raw data, as well as the full annotated assembly, are now available to the community in a dedicated web portal at http://www.corkoakdb.org.

Conclusions

This genomic resource represents the first trancriptome study in a cork producing species. It can be explored to develop new tools and approaches to understand stress responses and developmental processes in forest trees, as well as the molecular cascades underlying cork differentiation and disease response.     

Identification and comparative analysis of the Eriocheir sinensis microRNA transcriptome response to Spiroplasma eriocheiris infection using a deep sequencing approach

The Chinese mitten crab Eriocheir sinensis is one of the most important freshwater aquaculture crustacean species in China. MicroRNAs (miRNAs) are small non-coding RNAs that are important effectors in the intricate host-pathogen interaction network. To increase the repertoire of miRNAs characterized in crustaceans and to examine the relationship between host miRNA expression and pathogen infection, we used the Illumina/Solexa deep sequencing technology to sequence two small RNA libraries prepared from haemocytes of E. sinensis under normal conditions and during infection with Spiroplasma eriocheiris. The high-throughput sequencing resulted in approximately 30,975,151 and 30,826,277 raw reads corresponding to 12,077,088 and 16,271,545 high-quality mappable reads for the normal and infected haemocyte samples, respectively. Bioinformatic analyses identified 735 unique miRNAs, including 36 that are conserved in crustaceans, 134 that are novel to crabs but are present in other arthropods (PN-type), and 565 that are completely new (PC-type). Two hundred twenty-eight unique miRNAs displayed significant differential expression between the normal and infected haemocyte samples (p < 0.0001). Of these, 133 (58%) were significantly up-regulated and 95 (42%) were significantly down-regulated upon challenge with S. eriocheiris. Real-time quantitative PCR (RT-qPCR) experiments were preformed for 10 miRNAs of the two samples, and agreement was found between the sequencing and RT-qPCR data. To our knowledge, this is the first report of comprehensive identification of E. sinensis miRNAs and of expression analysis of E. sinensis miRNAs after exposure to S. eriocheiris. Many miRNAs were differentially regulated when exposed to the pathogen, and these findings support the hypothesis that certain miRNAs might be essential in host-pathogen interactions. Our results suggest that elucidation of the molecular mechanisms responsible for miRNA regulation of the host’s innate immune system should help with the development of new control strategies to prevent or treat S. eriocheiris infections in crustaceans.     

Osprey: a comprehensive tool employing novel methods for the design of oligonucleotides for DNA sequencing and microarrays

We have developed a software package called Osprey for the calculation of optimal oligonucleotides for DNA sequencing and the creation of microarrays based on either PCR-products or directly spotted oligomers. It incorporates a novel use of position-specific scoring matrices, for the sensitive and specific identification of secondary binding sites anywhere in the target sequence. Using accelerated hardware is faster and more efficient than the traditional pairwise alignments used in most oligo-design software. Osprey consists of a module for target site selection based on user input, novel utilities for dealing with problematic sequences such as repeats, and a common code base for the identification of optimal oligonucleotides from the target list. Overall, these improvements provide a program that, without major increases in run time, reflects current DNA thermodynamics models, improves specificity and reduces the user's data preprocessing and parameterization requirements. Using a TimeLogic™ hardware accelerator, we report up to 50-fold reduction in search time versus a linear search strategy. Target sites may be derived from computer analysis of DNA sequence assemblies in the case of sequencing efforts, or genome or EST analysis in the case of microarray development in both prokaryotes and eukaryotes.