Expressed sequence tags: an overview

Label-Free Sample Preparation for Single-Cell Proteomics

In recent years, single-cell proteomics (SCP) has become a valuable addition to other single-cell omics technologies for studying cellular heterogeneity. The amount of protein in a single cell is very limited, and in contrast to sequencing techniques, there are currently no means for protein amplification. Therefore, most single-cell proteomics approaches aim to maximize sample preparation efficiency while minimizing peptide loss. By reducing processing volumes to sub-microliters and avoiding manual transfer steps that could lead to peptide loss, peptide recovery, and the robustness of SCP workflows have been significantly improved. In this chapter, we describe a protocol for label-free SCP sample preparation using the cellenONE® platform and the proteoCHIP LF 48 substrate prior to analysis with high-performance liquid chromatography-mass spectrometry.

Housekeeping gene stability in Pesudomonas aeruginosa PAO1 under the pressure of commonly used antibiotics in molecular microbiology assays

Pseudomonas aeruginosa is an opportunistic human pathogen notorious for its remarkable capacity of multi-drug resistance, and has become one of the most important model bacteria in clinical bacteriology research. Quantitative real-time PCR is a reliable method widely used in gene expression analysis, for which the selection of a set of appropriate housekeeping genes is a key prerequisite for the accuracy of the results. However, it is easy to overlook that the expression level of housekeeping gene may vary in different conditions, especially in the condition of molecular microbiology assays, where tested strains are generally cultured under the pre-set antibiotic selection pressures, and how this affects the stability of commonly used housekeeping genes remains unclear. In this study, the expression stability of ten classic housekeeping genes (algD, gyrA, anr, nadB, recA, fabD, proC, ampC, rpoS, and rpsL) under the pressure of eight laboratory commonly used antibiotics (kanamycin, gentamycin, tetracycline, chloramphenicol, hygromycin B, apramycin, tellurite, and zeocin) were tested. Results showed that the stability of housekeeping gene expression was indeed affected by the types of antibiotics added, and of course the best reference gene set varied for different antibiotics. This study provides a comprehensive summary of the effects of laboratory antibiotics on the stability of housekeeping genes in P. aeruginosa, highlighting the necessity to select housekeeping genes according to the type of antibiotics used in the initial stage of experiment.

Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger

More over half of the world's population depends on rice as a major food crop. Rice (Oryza sativa L.) is vulnerable to abiotic challenges including drought, cold, and salinity since it grown in semi-aquatic, tropical, or subtropical settings. Abiotic stress resistance has bred into rice plants since the earliest rice cultivation techniques. Prior to the discovery of the genome, abiotic stress-related genes were identified using forward genetic methods, and abiotic stress-tolerant lines have developed using traditional breeding methods. Dynamic transcriptome expression represents the degree of gene expression in a specific cell, tissue, or organ of an individual organism at a specific point in its growth and development. Transcriptomics can reveal the expression at the entire genome level during stressful conditions from the entire transcriptional level, which can be helpful in understanding the intricate regulatory network relating to the stress tolerance and adaptability of plants. Rice (Oryza sativa L.) gene families found comparatively using the reference genome sequences of other plant species, allowing for genome-wide identification. Transcriptomics via gene expression profiling which have recently dominated by RNA-seq complements genomic techniques. The identification of numerous important qtl,s genes, promoter elements, transcription factors and miRNAs involved in rice response to abiotic stress was made possible by all of these genomic and transcriptomic techniques. The use of several genomes and transcriptome methodologies to comprehend rice (Oryza sativa, L.) ability to withstand abiotic stress have been discussed in this review.

Investigation of genetic relationships within three Miscanthus species using SNP markers identified with SLAF-seq

BACKGROUND

Miscanthus, which is a leading dedicated-energy grass in Europe and in parts of Asia, is expected to play a key role in the development of the future bioeconomy. However, due to its complex genetic background, it is difficult to investigate phylogenetic relationships in this genus. Here, we investigated 50 Miscanthus germplasms: 1 female parent (M. lutarioriparius), 30 candidate male parents (M. lutarioriparius, M. sinensis, and M. sacchariflorus), and 19 offspring. We used high-throughput Specific-Locus Amplified Fragment sequencing (SLAF-seq) to identify informative single nucleotide polymorphisms (SNPs) in all germplasms.

RESULTS

We identified 257,889 SLAF tags, of which 87,162 were polymorphic. Each tag was 264-364 bp long. The obtained 724,773 population SNPs were used to investigate genetic relationships within three species of Miscanthus. We constructed a phylogenetic tree of the 50 germplasms using the obtained SNPs and grouped them into two clades: one clade comprised of M. sinensis alone and the other one included the offspring, M. lutarioriparius, and M. sacchariflorus. Genetic cluster analysis had revealed that M. lutarioriparius germplasm C3 was the most likely male parent of the offspring.

CONCLUSIONS

As a high-throughput sequencing method, SLAF-seq can be used to identify informative SNPs in Miscanthus germplasms and to rapidly characterize genetic relationships within this genus. Our results will support the development of breeding programs with the focus on utilizing Miscanthus cultivars with elite biomass- or fiber-production potential for the developing bioeconomy.

Unraveling the versatility of CCD4: Metabolic engineering, transcriptomic and computational approaches

Carotenoids are economically valuable isoprenoids synthesized by plants and microorganisms, which play a paramount role in their overall growth and development. Carotenoid cleavage dioxygenases are a vast group of enzymes that specifically cleave thecarotenoids to produce apocarotenoids. Recently, CCDs are a subject of talk because of their contributions to different aspects of plant growth and due to their significance in the production of economically valuable apocarotenoids. Among them, CCD4 stands unique because of its versatility in performing metabolic roles. This review focuses on the multiple functionalities of CCD4 like pigmentation, volatile apocarotenoid production, stress responses, etc. Interestingly, through our literature survey we arrived at a conclusion that CCD4 could perform functions of other carotenoid cleaving enzymes.The metabolic engineering, transcriptomic, and computational approaches adopted to reveal the contributions of CCD4 were also considered here for the study.Phylogenetic analysis was performed to delve into the evolutionary relationships of CCD4 in different plant groups. A tree of 81CCD genes from 64 plant species was constructed, signifying the presence of well-conserved families. Gene structures were illustrated and the difference in the number and position of exons could be considered as a factor behind functional versatility and substrate tolerance of CCD4 in different plants.

Guidelines for Setting Up a mRNA Sequencing Experiment and Best Practices for Bioinformatic Data Analysis

RNA-sequencing, commonly referred to as RNA-seq, is the most recently developed method for the analysis of transcriptomes. It uses high-throughput next-generation sequencing technologies and has revolutionized our understanding of the complexity and dynamics of whole transcriptomes.In this chapter, we recall the key developments in transcriptome analysis and dissect the different steps of the general workflow that can be run by users to design and perform a mRNA-seq experiment as well as to process mRNA-seq data obtained by the Illumina technology. The chapter proposes guidelines for completing a mRNA-seq study properly and makes available recommendations for best practices based on recent literature and on the latest developments in technology and algorithms. We also remark the large number of choices available (especially for bioinformatic data analysis) in front of which the scientist may be in trouble.In the last part of the chapter we discuss the new frontiers of single-cell RNA-seq and isoform sequencing by long read technology.

TagSeq for gene expression in non-model plants: A pilot study at the Santa Rita Experimental Range NEON core site

PREMISE

TagSeq is a cost-effective approach for gene expression studies requiring a large number of samples. To date, TagSeq studies in plants have been limited to those with a high-quality reference genome. We tested the suitability of reference transcriptomes for TagSeq in non-model plants, as part of a study of natural gene expression variation at the Santa Rita Experimental Range National Ecological Observatory Network (NEON) core site.

METHODS

Tissue for TagSeq was sampled from multiple individuals of four species (Bouteloua aristidoides and Eragrostis lehmanniana [Poaceae], Tidestromia lanuginosa [Amaranthaceae], and Parkinsonia florida [Fabaceae]) at two locations on three dates (56 samples total). One sample per species was used to create a reference transcriptome via standard RNA-seq. TagSeq performance was assessed by recovery of reference loci, specificity of tag alignments, and variation among samples.

RESULTS

A high fraction of tags aligned to each reference and mapped uniquely. Expression patterns were quantifiable for tens of thousands of loci, which revealed consistent spatial differentiation in expression for all species.

DISCUSSION

TagSeq using de novo reference transcriptomes was an effective approach to quantifying gene expression in this study. Tags were highly locus specific and generated biologically informative profiles for four non-model plant species.

Cancer-testis antigens in canine histiocytic sarcoma and other malignancies

Cancer-testis antigens (CTAs) are a category of self proteins aberrantly expressed in diverse malignancies, mostly solid tumours, due to epigenetic de-repression. Normally expressed only in fetal or gametogenic tissues, CTAs are tantalizing immunotherapy targets, since autoimmunity risks appear minimal. Few prevalent CTAs have been identified in human hematologic cancers, and just two in their veterinary counterparts. We sought to discover new CTAs in canine hematologic cancers such as histiocytic sarcoma (HS) and lymphoma to foster immunotherapy development. To accomplish this, the ligandome binding the dog leukocyte antigen (DLA)-88*508:01 class I allele overexpressed in an HS line was searched by mass spectrometry to identify possible CTA-derived peptides, which could serve as CD8+ T-cell epitopes. Twenty-two peptides mapped to 5 human CTAs and 12 additional proteins with CTA characteristics. Expression of five promising candidates was then evaluated in tumour and normal tissue by quantitative and end-point RT-PCR. The ortholog of an established CTA, IGF2BP3, had unexpectedly high expression in peripheral blood mononuclear cells (PBMCs). Four other testis-enhanced proteins were also assessed. AKR1E2, SPECC1 and TPX2 were expressed variably in HS and T-cell lymphoma biopsies, but also at high levels in critical tissues, including kidney, brain and marrow, diminishing their utility. A more tissue-restricted candidate, NT5C1B, was detected in T-cell lymphomas, but also at low levels in some normal dog tissues. These results illustrate the feasibility of discovering canine CTAs by a reverse approach, proceeding from identification of MHC class I-presented peptides to a comparative RNA expression survey of tumours and normal tissues.

Identification of Expressed miRNAs in Human Rheumatoid Arthritis Using Computational Approach - Discovery of a New miR-7167 from Human

BACKGROUND

Rheumatoid Arthritis (RA) is a chronic inflammatory and autoimmune disease leading to bones and joints destruction. It is one of the major causes of lifetime disability and mortality among humans in the developing and developed countries. It was evident that epigenetic dysregulation is related to the pathogenesis of RA. MicroRNAs (miRNAs) are small non-coding RNAs that are epigenetic regulators for diverse biological processes and also provided novel molecular insights in the formation of arthritis.

OBJECTIVE

The influences of miRNAs in the alteration of gene regulation during the pathogenesis of arthritis were exposed in recent years.

METHOD

The computational approach to identify miRNA through EST-based homology is more powerful, economical and time-efficient. In this study, we applied EST-based homology search to identify miRNAs responsible for the development of arthritis in human beings.

RESULTS

Our study on 36519 ESTs in human RA condition revealed the expression of four miRNAs, HSA-miR-198, HSA-miR-4647, has-miR-7167-5p and has-miR-7167-3p. The present study is the first report about has-miR-7167 that was homologous to Macaca mulatta.

CONCLUSION

The predicted targets of these identified miRNAs revealed many biological functions in the pathogenesis of RA. Further elaborated studies on these miRNAs will help to understand their function in the development of RA and the use of miRNAs as therapeutic targets in the future.

Ocimum basilicum miRNOME revisited: A cross kingdom approach

O. basilicum is medicinally important herb having inevitable role in human health. However, the mechanism of action is largely unknown. Present study aims to understand the mechanism of regulation of key human target genes that could plausibly modulated by O. basilicum miRNAs in cross kingdom manner using computational and system biology approach. O. basilicum miRNA sequences were retrieved and their corresponding human target genes were identified using psRNA target and interaction analysis of hub nodes. Six O. basilicum derived miRNAs were found to modulate 26 human target genes which were associated `with PI3K-AKTand MAPK signaling pathways with PTPN11, EIF2S2, NOS1, IRS1 and USO1 as top 5 Hub nodes. O. basilicum miRNAs not only regulate key human target genes having a significance in various diseases but also paves the path for future studies that might explore potential of miRNA mediated cross-kingdom regulation, prevention and treatment of various human diseases including cancer.

Advances in Fasciola hepatica research using 'omics' technologies

The liver fluke Fasciola hepatica is an economically important pathogen of livestock worldwide, as well as being an important neglected zoonosis. Parasite control is reliant on the use of drugs, particularly triclabendazole, which is effective against multiple parasite stages. However, the spread of parasites resistant to triclabendazole has intensified the pursuit for novel control strategies. Emerging 'omics' technologies are helping advance our understanding of liver fluke biology, specifically the molecules that act at the host-parasite interface and are central to infection, virulence and long-term survival within the definitive host. This review discusses the technological sequencing advances that have facilitated the unbiased analysis of liver fluke biology, resulting in an extensive range of 'omics' datasets. In addition, we highlight the 'omics' studies of host responses to F. hepatica infection that, when combined with the parasite datasets, provide the opportunity for integrated analyses of host-parasite interactions. These extensive datasets will form the foundation for future in-depth analysis of F. hepatica biology and development, and the search for new drug or vaccine interventions.

Differential gene expression analysis in ageing muscle and drug discovery perspectives

Identifying therapeutic target genes represents the key step in functional genomics-based therapies. Within this context, the disease heterogeneity, the exogenous factors and the complexity of genomic structure and function represent important challenges. The functional genomics aims to overcome such obstacles via identifying the gene functions and therefore highlight disease-causing genes as therapeutic targets. Genomic technologies promise to reshape the research on ageing muscle, exercise response and drug discovery. Herein, we describe the functional genomics strategies, mainly differential gene expression methods microarray, serial analysis of gene expression (SAGE), massively parallel signature sequence (MPSS), RNA sequencing (RNA seq), representational difference analysis (RDA), and suppression subtractive hybridization (SSH). Furthermore, we review these illustrative approaches that have been used to discover new therapeutic targets for some complex diseases along with the application of these tools to study the modulation of the skeletal muscle transcriptome.

Difficulty in obtaining the complete mRNA coding sequence at 5' region (5' end mRNA artifact): Causes, consequences in biology and medicine and possible solutions for obtaining the actual amino acid sequence of proteins (Review)

The known difficulty in obtaining the actual full length, complete sequence of a messenger RNA (mRNA) may lead to the erroneous determination of its coding sequence at the 5' region (5' end mRNA artifact), and consequently to the wrong assignment of the translation start codon, leading to the inaccurate prediction of the encoded polypeptide at its amino terminus. Among the known human genes whose study was affected by this artifact, we can include disco interacting protein 2 homolog A (DIP2A; KIAA0184), Down syndrome critical region 1 (DSCR1), SON DNA binding protein (SON), trefoil factor 3 (TFF3) and URB1 ribosome biogenesis 1 homolog (URB1; KIAA0539) on chromosome 21, as well as receptor for activated C kinase 1 (RACK1, also known as GNB2L1), glutaminyl‑tRNA synthetase (QARS) and tyrosyl-DNA phosphodiesterase 2 (TDP2) along with another 474 loci, including interleukin 16 (IL16). In this review, we discuss the causes of this issue, its quantitative incidence in biomedical research, the consequences in biology and medicine, and the possible solutions for obtaining the actual amino acid sequence of proteins in the post-genomics era.

Characterization of the gene encoding component C3 of the complement system from the spider Loxosceles laeta venom glands: Phylogenetic implications

A transcriptome analysis of the venom glands of the spider Loxosceles laeta, performed by our group, in a previous study (Fernandes-Pedrosa et al., 2008), revealed a transcript with a sequence similar to the human complement component C3. Here we present the analysis of this transcript. cDNA fragments encoding the C3 homologue (Lox-C3) were amplified from total RNA isolated from the venom glands of L. laeta by RACE-PCR. Lox-C3 is a 5178 bps cDNA sequence encoding a 190kDa protein, with a domain configuration similar to human C3. Multiple alignments of C3-like proteins revealed two processing sites, suggesting that Lox-C3 is composed of three chains. Furthermore, the amino acids consensus sequences for the thioester was found, in addition to putative sequences responsible for FB binding. The phylogenetic analysis showed that Lox-C3 belongs to the same group as two C3 isoforms from the spider Hasarius adansoni (Family Salcitidae), showing 53% homology with these. This is the first characterization of a Loxosceles cDNA sequence encoding a human C3 homologue, and this finding, together with our previous finding of the expression of a FB-like molecule, suggests that this spider species also has a complement system. This work will help to improve our understanding of the innate immune system in these spiders and the ancestral structure of C3.

Understanding Haemonchus contortus Better Through Genomics and Transcriptomics

Parasitic roundworms (nematodes) cause substantial mortality and morbidity in animals globally. The barber's pole worm, Haemonchus contortus, is one of the most economically significant parasitic nematodes of small ruminants worldwide. Although this and related nematodes can be controlled relatively well using anthelmintics, resistance against most drugs in common use has become a major problem. Until recently, almost nothing was known about the molecular biology of H. contortus on a global scale. This chapter gives a brief background on H. contortus and haemonchosis, immune responses, vaccine research, chemotherapeutics and current problems associated with drug resistance. It also describes progress in transcriptomics before the availability of H. contortus genomes and the challenges associated with such work. It then reviews major progress on the two draft genomes and developmental transcriptomes of H. contortus, and summarizes their implications for the molecular biology of this worm in both the free-living and the parasitic stages of its life cycle. The chapter concludes by considering how genomics and transcriptomics can accelerate research on Haemonchus and related parasites, and can enable the development of new interventions against haemonchosis.

Transcriptome sequencing of a thalloid bryophyte; Dumortiera hirsuta (Sw) Nees: assembly, annotation, and marker discovery

Bryophytes are the first land plants but are scarcely studied at the molecular level. Here, we report transcriptome sequencing and functional annotation of Dumortiera hirsuta, as a representative bryophyte. Approximately 0.5 million reads with ~195 Mb data were generated by sequencing of mRNA using 454 pyrosequencer. De novo assembly of reads yielded 85,240 unigenes (12,439 contigs and 72,801 singletons). BlastX search at NCBI-NR database showed similarity of 33,662 unigenes with 10-(10) e-value. A total of 23,685 unigenes were annotated at TAIR10 protein database. The annotated unigenes were further classified using the Gene Ontology. Analysis at Kyoto Encyclopedia of Genes and Genomes pathway database identified 95 pathways with significant scores, among which metabolic and biosynthesis of secondary metabolite were the major ones. Phenylpropanoid pathway was elucidated and selected genes were characterized by real time qPCR. We identified 447 transcription factors belonging to 41 families and 1594 eSSRs in 1479 unigenes. D. hirsuta unigenes showed homology across the taxa from algae to angiosperm indicating their role as the connecting link between aquatic and terrestrial plants. This could be a valuable genomic resource for molecular and evolutionary studies. Further, it sheds light for the isolation and characterization of new genes with unique functions.

De novo assembly and characterization of the transcriptome of the pancreatic fluke Eurytrema pancreaticum (trematoda: Dicrocoeliidae) using Illumina paired-end sequencing

Eurytrema pancreaticum is one of the most common trematodes living in the pancreatic and bile ducts of ruminants and also occasionally infects humans, causing eurytremiasis. In spite of its economic and medical importance, very little is known about the genomic resources of this parasite. Herein, we performed de novo sequencing, assembly and characterization of the transcriptome of adult E. pancreaticum. Approximately 36.4 million high-quality clean reads were obtained, and the length of the transcript contigs ranged from 66 to 19,968 nt with mean length of 479 nt and N50 length of 1094 nt, and then 23,573 unigenes were assembled. Of these unigenes, 15,353 (65.1%) were annotated by blast searches against the NCBI non-redundant protein database. Among these, 15,267 (64.8%), 2732 (11.6%) and 10,354 (43.9%) of the unigenes had significant similarity with proteins in the NR, NT and Swiss-Prot databases, respectively. 5510 (23.4%) and 4567 (19.4%) unigenes were assigned to GO and COG, respectively. 8886 (37.7%) unigenes were identified and mapped onto 254 pathways in the KEGG Pathway database. Furthermore, we found that 105 (1.18%) unigenes were related to pancreatic secretion and 61 (0.7%) to pancreatic cancer. The present study represents the first transcriptome of any members of the family Dicrocoeliidae, which has little genomic information available in the public databases. The novel transcriptome of E. pancreaticum should provide a useful resource for designing new strategies against pancreatic flukes and other trematodes of human and animal health significance.

De Novo Transcriptome Analysis of Allium cepa L. (Onion) Bulb to Identify Allergens and Epitopes

Allium cepa (onion) is a diploid plant with one of the largest nuclear genomes among all diploids. Onion is an example of an under-researched crop which has a complex heterozygous genome. There are no allergenic proteins and genomic data available for onions. This study was conducted to establish a transcriptome catalogue of onion bulb that will enable us to study onion related genes involved in medicinal use and allergies. Transcriptome dataset generated from onion bulb using the Illumina HiSeq 2000 technology showed a total of 99,074,309 high quality raw reads (~20 Gb). Based on sequence homology onion genes were categorized into 49 different functional groups. Most of the genes however, were classified under 'unknown' in all three gene ontology categories. Of the categorized genes, 61.2% showed metabolic functions followed by cellular components such as binding, cellular processes; catalytic activity and cell part. With BLASTx top hit analysis, a total of 2,511 homologous allergenic sequences were found, which had 37–100% similarity with 46 different types of allergens existing in the database. From the 46 contigs or allergens, 521 B-cell linear epitopes were identified using BepiPred linear epitope prediction tool. This is the first comprehensive insight into the transcriptome of onion bulb tissue using the NGS technology, which can be used to map IgE epitopes and prediction of structures and functions of various proteins.

Insect phylogenomics

Phylogenomics, the integration of phylogenetics with genome data, has emerged as a powerful approach to study the evolution and systematics of species. Recently, several studies employing phylogenomic tools have provided better insights into insect evolution. Next-generation sequencing methods are now increasingly used by entomologists to generate genomic and transcript sequences of various insect species and strains. These data provide opportunities for comparative genomics and large-scale multigene phylogenies of diverse lineages of insects. Phy-logenomic investigations help us to better understand systematic and evolutionary relationships of insect species that play important roles as herbivores, predators, detritivores, pollinators and disease vectors. It is important that we critically assess the prospects and limitations of phylogenomic methods. In this review, I describe the current status, outline the major challenges and remark on potential future applications of phylogenomic tools in studying insect systematics and evolution.

Transcriptomic analysis of Siberian ginseng (Eleutherococcus senticosus) to discover genes involved in saponin biosynthesis

Background

Eleutherococcus senticosus, Siberian ginseng, is a highly valued woody medicinal plant belonging to the family Araliaceae. E. senticosus produces a rich variety of saponins such as oleanane-type, noroleanane-type, 29-hydroxyoleanan-type, and lupane-type saponins. Genomic or transcriptomic approaches have not been used to investigate the saponin biosynthetic pathway in this plant.

Result

In this study, de novo sequencing was performed to select candidate genes involved in the saponin biosynthetic pathway. A half-plate 454 pyrosequencing run produced 627,923 high-quality reads with an average sequence length of 422 bases. De novo assembly generated 72,811 unique sequences, including 15,217 contigs and 57,594 singletons. Approximately 48,300 (66.3%) unique sequences were annotated using BLAST similarity searches. All of the mevalonate pathway genes for saponin biosynthesis starting from acetyl-CoA were isolated. Moreover, 206 reads of cytochrome P450 (CYP) and 145 reads of uridine diphosphate glycosyltransferase (UGT) sequences were isolated. Based on methyl jasmonate (MeJA) treatment and real-time PCR (qPCR) analysis, 3 CYPs and 3 UGTs were finally selected as candidate genes involved in the saponin biosynthetic pathway.

Conclusions

The identified sequences associated with saponin biosynthesis will facilitate the study of the functional genomics of saponin biosynthesis and genetic engineering of E. senticosus.

Electronic supplementary material

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

Whole transcriptome analysis using next-generation sequencing of sterile-cultured Eisenia andrei for immune system research

Recently, earthworms have become a useful model for research into the immune system, and it is expected that results obtained using this model will shed light on the sophisticated vertebrate immune system and the evolution of the immune response, and additionally help identify new biomolecules with therapeutic applications. However, for earthworms to be used as a genetic model of the invertebrate immune system, basic molecular and genetic resources, such as an expressed sequence tag (EST) database, must be developed for this organism. Next-generation sequencing technologies have generated EST libraries by RNA-seq in many model species. In this study, we used Illumina RNA-sequence technology to perform a comprehensive transcriptome analysis using an RNA sample pooled from sterile-cultured Eisenia andrei. All clean reads were assembled de novo into 41,423 unigenes using the Trinity program. Using this transcriptome data, we performed BLAST analysis against the GenBank non-redundant (NR) database and obtained a total of 12,285 significant BLAST hits. Furthermore, gene ontology (GO) analysis assigned 78 unigenes to 24 immune class GO terms. In addition, we detected a unigene with high similarity to beta-1,3-glucuronyltransferase 1 (GlcAT-P), which mediates a glucuronyl transfer reaction during the biosynthesis of the carbohydrate epitope HNK-1 (human natural killer-1, also known as CD57), a marker of NK cells. The identified transcripts will be used to facilitate future research into the immune system using E. andrei.

Impact of next-generation technologies on exploring socioeconomically important parasites and developing new interventions

High-throughput molecular and computer technologies have become instrumental for systems biological explorations of pathogens, including parasites. For instance, investigations of the transcriptomes of different developmental stages of parasitic nematodes give insights into gene expression, regulation and function in a parasite, which is a significant step to understanding their biology, as well as interactions with their host(s) and disease. This chapter (1) gives a background on some key parasitic nematodes of socioeconomic importance, (2) describes sequencing and bioinformatic technologies for large-scale studies of the transcriptomes and genomes of these parasites, (3) provides some recent examples of applications and (4) emphasizes the prospects of fundamental biological explorations of parasites using these technologies for the development of new interventions to combat parasitic diseases.

Genome-wide transcriptome analysis of genes involved in flavonoid biosynthesis between red and white strains of Magnolia sprengeri pamp

BACKGROUND

Magnolia sprengeri Pamp is one of the most highly valuable medicinal and ornamental plants of the Magnolia Family. The natural color of M. sprengeri is variable. The complete genome sequence of M. sprengeri is not available; therefore we sequenced the transcriptome of white and red petals of M. sprengeri using Illumina technology. We focused on the identity of structural and regulatory genes encoding the enzymes involved in the determination of flower color.

RESULTS

We sequenced and annotated a reference transcriptome for M. sprengeri, and aimed to capture the transcriptional determinanats of flower color. We sequenced a normalized cDNA library of white and red petals using Illumina technology. The resulting reads were assembled into 77,048 unique sequences, of which 28,243 could be annotated by Gene Ontology (GO) analysis, while 48,805 transcripts lacked GO annotation. The main enzymes involved in the flavonoid biosynthesis, such as phenylalanine ammonia-Lyase, cinnamat-4-Hydroxylase, dihydroflavonol-4-reductase, flavanone 3-hydroxylase, flavonoid-3'-hydroxylase, flavonol synthase, chalcone synthase and anthocyanidin synthase, were identified in the transcriptome. A total of 270 transcription factors were sorted into three families, including MYB, bHLH and WD40 types. Among these transcription factors, eight showed 4-fold or greater changes in transcript abundance in red petals compared with white petals. High-performance liquid chromatography analysis of anthocyanin compositions showed that the main anthocyanin in the petals of M. sprengeri is cyanidin-3-O-glucoside chloride and its content in red petals was 26-fold higher than that in white petals.

CONCLUSION

This study presents the first next-generation sequencing effort and transcriptome analysis of a non-model plant from the Family Magnoliaceae. Genes encoding key enzymes were identified and the metabolic pathways involved in biosynthesis and catabolism of M. sprengeri flavonoids were reconstructed. Identification of these genes and pathways adds to the current knowledge of the molecular biology and biochemistry of their production in plant. Such insights into the mechanisms supporting metabolic processes could be used to genetically to enhance flower color among members of the Magnoliaceae.

De novo sequencing and comparative analysis of holy and sweet basil transcriptomes

Background

Ocimum L. of family Lamiaceae is a well known genus for its ethnobotanical, medicinal and aromatic properties, which are attributed to innumerable phenylpropanoid and terpenoid compounds produced by the plant. To enrich genomic resources for understanding various pathways, de novo transcriptome sequencing of two important species, O. sanctum and O. basilicum, was carried out by Illumina paired-end sequencing.

Results

The sequence assembly resulted in 69117 and 130043 transcripts with an average length of 1646 ± 1210.1 bp and 1363 ± 1139.3 bp for O. sanctum and O. basilicum, respectively. Out of the total transcripts, 59648 (86.30%) and 105470 (81.10%) from O. sanctum and O. basilicum, and respectively were annotated by uniprot blastx against Arabidopsis, rice and lamiaceae. KEGG analysis identified 501 and 952 transcripts from O. sanctum and O. basilicum, respectively, related to secondary metabolism with higher percentage of transcripts for biosynthesis of terpenoids in O. sanctum and phenylpropanoids in O. basilicum. Higher digital gene expression in O. basilicum was validated through qPCR and correlated to higher essential oil content and chromosome number (O. sanctum, 2n = 16; and O. basilicum, 2n = 48). Several CYP450 (26) and TF (40) families were identified having probable roles in primary and secondary metabolism. Also SSR and SNP markers were identified in the transcriptomes of both species with many SSRs linked to phenylpropanoid and terpenoid pathway genes.

Conclusion

This is the first report of a comparative transcriptome analysis of Ocimum species and can be utilized to characterize genes related to secondary metabolism, their regulation, and breeding special chemotypes with unique essential oil composition in Ocimum.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-588) contains supplementary material, which is available to authorized users.

The first set of expressed sequence tags (EST) from the medicinal mushroom Agaricus subrufescens delivers resource for gene discovery and marker development

Agaricus subrufescens is one of the most important culinary-medicinal cultivable mushrooms with potentially high-added-value products and extended agronomical valorization. The development of A. subrufescens-related technologies is hampered by, among others, the lack of suitable molecular tools. Thus, this mushroom is considered as a genomic orphan species with a very limited number of available molecular markers or sequences. To fill this gap, this study reports the generation and analysis of the first set of expressed sequence tags (EST) for A. subrufescens. cDNA fragments obtained from young sporophores (SP) and vegetative mycelium in liquid culture (CL) were sequenced using 454 pyrosequencing technology. After assembly process, 4,989 and 5,125 sequences were obtained in SP and CL libraries, respectively. About 87% of the EST had significant similarity with Agaricus bisporus-predicted proteins, and 79% correspond to known proteins. Functional categorization according to Gene Ontology could be assigned to 49% of the sequences. Some gene families potentially involved in bioactive compound biosynthesis could be identified. A total of 232 simple sequence repeats (SSRs) were identified, and a set of 40 EST-SSR polymorphic markers were successfully developed. This EST dataset provides a new resource for gene discovery and molecular marker development. It constitutes a solid basis for further genetic and genomic studies in A. subrufescens.

The origins and progress of genomics research on Tef (Eragrostis tef)

Tef, Eragrostis tef (Zucc.) Trotter, is the most important cereal in Ethiopia. Tef is cultivated by more than five million small-scale farmers annually and constitutes the staple food for more than half of the population of 80 million. The crop is preferred by both farmers and consumers due to its beneficial traits associated with its agronomy and utilization. The genetic and phenotypic diversity of tef in Ethiopia is a national treasure of potentially global importance. In order for this diversity to be effectively conserved and utilized, a better understanding at the genomic level is necessary. In the recent years, tef has become the subject of genomic research in Ethiopia and abroad. Genomic-assisted tef improvement holds tremendous potential for improving productivity, thereby benefiting the smallholder farmers who have cultivated and relied on the crop for thousands of years. It is hoped that such research endeavours will provide solutions to some of the age-old problems of tef's husbandry. In this review, we provide a brief description of the genesis and progress of tef genomic research to date, suggest ways to utilize the genomic tools developed so far, discuss the potential of genomics to enable sustainable conservation and use of tef genetic diversity and suggest opportunities for the future research.

Genomic resources in fruit plants: an assessment of current status

The availability of many genomic resources such as genome sequences, functional genomics resources including microarrays and RNA-seq, sufficient numbers of molecular markers, express sequence tags (ESTs) and high-density genetic maps is causing a rapid acceleration of genetics and genomic research of many fruit plants. This is leading to an increase in our knowledge of the genes that are linked to many horticultural and agronomically important traits. Recently, some progress has also been made on the identification and functional analysis of miRNAs in some fruit plants. This is one of the most active research fields in plant sciences. The last decade has witnessed development of genomic resources in many fruit plants such as apple, banana, citrus, grapes, papaya, pears, strawberry etc.; however, many of them are still not being exploited. Furthermore, owing to lack of resources, infrastructure and research facilities in many lesser-developed countries, development of genomic resources in many underutilized or less-studied fruit crops, which grow in these countries, is limited. Thus, research emphasis should be given to those fruit crops for which genomic resources are relatively scarce. The development of genomic databases of these less-studied fruit crops will enable biotechnologists to identify target genes that underlie key horticultural and agronomical traits. This review presents an overview of the current status of the development of genomic resources in fruit plants with the main emphasis being on genome sequencing, EST resources, functional genomics resources including microarray and RNA-seq, identification of quantitative trait loci and construction of genetic maps as well as efforts made on the identification and functional analysis of miRNAs in fruit plants.

Structural and functional characterization of the chaperone Hsp70 from sugarcane. Insights into conformational changes during cycling from cross-linking/mass spectrometry assays

Hsp70 cycles from an ATP-bound state, in which the affinity for unfolded polypeptides is low, to an ADP-bound state, in which the affinity for unfolded polypeptides is high, to assist with cell proteostasis. Such cycling also depends on co-chaperones because these proteins control both the Hsp70 ATPase activity and the delivery of unfolded polypeptide chains. Although it is very important, structural information on the entire protein is still scarce. This work describes the first cloning of a cDNA predicted to code for a cytosolic Saccharum spp. (sugarcane) Hsp70, named SsHsp70 here, the purification of the recombinant protein and the characterization of its structural conformation in solution by chemical cross-linking coupled to mass spectrometry. The in vivo expression of SsHsp70 in sugarcane extracts was confirmed by Western blot. Recombinant SsHsp70 was monomeric, both ADP and ATP binding increased its stability and it was efficient in cooperating with co-chaperones: ATPase activity was stimulated by Hsp40s, and it aided the refolding of an unfolded polypeptide delivered by a member of the small Hsp family. The structural conformation results favor a model in which nucleotide-free SsHsp70 is highly dynamic and may fluctuate among different conformations that may resemble those in which nucleotide is bound.

BIOLOGICAL SIGNIFICANCE

Validation of a sugarcane EST as a true mRNA that encodes a cytosolic Hsp70 (SsHsp70) as confirmed by in vivo expression and characterization of the structure and function of the recombinant protein. SsHsp70 was monomeric, both ADP and ATP binding increased its stability and was efficient in interacting and cooperating with co-chaperones to enhance ATPase activity and refold unfolded proteins. The conformation of nucleotide-free SsHsp70 in solution was much more dynamic than suggested by crystal structures of other Hsp70s. This article is part of a Special Issue entitled: Environmental and structural proteomics.

Understanding non-coding DNA regions in yeast

Non-coding transcripts play an important role in gene expression regulation in all species, including budding and fission yeast. Such regulatory transcripts include intergenic ncRNA (non-coding RNA), 5′ and 3′ UTRs, introns and antisense transcripts. In the present review, we discuss advantages and limitations of recently developed sequencing techniques, such as ESTs, DNA microarrays, RNA-Seq (RNA sequencing), DRS (direct RNA sequencing) and TIF-Seq (transcript isoform sequencing). We provide an overview of methods applied in yeast and how each of them has contributed to our knowledge of gene expression regulation and transcription.

Transcriptome sequencing and microarray design for functional genomics in the extremophile Arabidopsis relative Thellungiella salsuginea (Eutrema salsugineum)

Background

Most molecular studies of plant stress tolerance have been performed with Arabidopsis thaliana, although it is not particularly stress tolerant and may lack protective mechanisms required to survive extreme environmental conditions. Thellungiella salsuginea has attracted interest as an alternative plant model species with high tolerance of various abiotic stresses. While the T. salsuginea genome has recently been sequenced, its annotation is still incomplete and transcriptomic information is scarce. In addition, functional genomics investigations in this species are severely hampered by a lack of affordable tools for genome-wide gene expression studies.

Results

Here, we report the results of Thellungiella de novo transcriptome assembly and annotation based on 454 pyrosequencing and development and validation of a T. salsuginea microarray. ESTs were generated from a non-normalized and a normalized library synthesized from RNA pooled from samples covering different tissues and abiotic stress conditions. Both libraries yielded partially unique sequences, indicating their necessity to obtain comprehensive transcriptome coverage. More than 1 million sequence reads were assembled into 42,810 unigenes, approximately 50% of which could be functionally annotated. These unigenes were compared to all available Thellungiella genome sequence information. In addition, the groups of Late Embryogenesis Abundant (LEA) proteins, Mitogen Activated Protein (MAP) kinases and protein phosphatases were annotated in detail. We also predicted the target genes for 384 putative miRNAs. From the sequence information, we constructed a 44 k Agilent oligonucleotide microarray. Comparison of same-species and cross-species hybridization results showed superior performance of the newly designed array for T. salsuginea samples. The developed microarrays were used to investigate transcriptional responses of T. salsuginea and Arabidopsis during cold acclimation using the MapMan software.

Conclusions

This study provides the first comprehensive transcriptome information for the extremophile Arabidopsis relative T. salsuginea. The data constitute a more than three-fold increase in the number of publicly available unigene sequences and will greatly facilitate genome annotation. In addition, we have designed and validated the first genome-wide microarray for T. salsuginea, which will be commercially available. Together with the publicly available MapMan software this will become an important tool for functional genomics of plant stress tolerance.

Whole transcriptome analysis using next-generation sequencing of model species Setaria viridis to support C4 photosynthesis research

Setaria viridis is an emerging model species for genetic studies of C4 photosynthesis. Many basic molecular resources need to be developed to support for this species. In this paper, we performed a comprehensive transcriptome analysis from multiple developmental stages and tissues of S. viridis using next-generation sequencing technologies. Sequencing of the transcriptome from multiple tissues across three developmental stages (seed germination, vegetative growth, and reproduction) yielded a total of 71 million single end 100 bp long reads. Reference-based assembly using Setaria italica genome as a reference generated 42,754 transcripts. De novo assembly generated 60,751 transcripts. In addition, 9,576 and 7,056 potential simple sequence repeats (SSRs) covering S. viridis genome were identified when using the reference based assembled transcripts and the de novo assembled transcripts, respectively. This identified transcripts and SSR provided by this study can be used for both reverse and forward genetic studies based on S. viridis.

De novo assembly, functional annotation and comparative analysis of Withania somnifera leaf and root transcriptomes to identify putative genes involved in the withanolides biosynthesis

Withania somnifera is one of the most valuable medicinal plants used in Ayurvedic and other indigenous medicine systems due to bioactive molecules known as withanolides. As genomic information regarding this plant is very limited, little information is available about biosynthesis of withanolides. To facilitate the basic understanding about the withanolide biosynthesis pathways, we performed transcriptome sequencing for Withania leaf (101L) and root (101R) which specifically synthesize withaferin A and withanolide A, respectively. Pyrosequencing yielded 8,34,068 and 7,21,755 reads which got assembled into 89,548 and 1,14,814 unique sequences from 101L and 101R, respectively. A total of 47,885 (101L) and 54,123 (101R) could be annotated using TAIR10, NR, tomato and potato databases. Gene Ontology and KEGG analyses provided a detailed view of all the enzymes involved in withanolide backbone synthesis. Our analysis identified members of cytochrome P450, glycosyltransferase and methyltransferase gene families with unique presence or differential expression in leaf and root and might be involved in synthesis of tissue-specific withanolides. We also detected simple sequence repeats (SSRs) in transcriptome data for use in future genetic studies. Comprehensive sequence resource developed for Withania, in this study, will help to elucidate biosynthetic pathway for tissue-specific synthesis of secondary plant products in non-model plant organisms as well as will be helpful in developing strategies for enhanced biosynthesis of withanolides through biotechnological approaches.

Transcriptome analysis using RNA-Seq

Transcriptome analysis by next-generation sequencing (RNA-seq) allows investigation of a transcriptome at unsurpassed resolution. One major benefit is that RNA-seq is independent of a priori knowledge on the sequence under investigation, thereby also allowing analysis of poorly characterized Plasmodium species. Here we provide a detailed protocol for RNA isolation and fragmentation, ribosomal RNA depletion, and cDNA synthesis that enables the preparation of a sequencing library from 1 to 2 μg of total RNA. Although we focus our discussion on the quantitative measurement of gene expression, this protocol is suited for many applications of RNA-seq and allows analysis of most RNA species.

Transcriptomics and comparative analysis of three antarctic notothenioid fishes

For the past 10 to 13 million years, Antarctic notothenioid fish have undergone extraordinary periods of evolution and have adapted to a cold and highly oxygenated Antarctic marine environment. While these species are considered an attractive model with which to study physiology and evolutionary adaptation, they are poorly characterized at the molecular level, and sequence information is lacking. The transcriptomes of the Antarctic fishes Notothenia coriiceps, Chaenocephalus aceratus, and Pleuragramma antarcticum were obtained by 454 FLX Titanium sequencing of a normalized cDNA library. More than 1,900,000 reads were assembled in a total of 71,539 contigs. Overall, 40% of the contigs were annotated based on similarity to known protein or nucleotide sequences, and more than 50% of the predicted transcripts were validated as full-length or putative full-length cDNAs. These three Antarctic fishes shared 663 genes expressed in the brain and 1,557 genes expressed in the liver. In addition, these cold-adapted fish expressed more Ub-conjugated proteins compared to temperate fish; Ub-conjugated proteins are involved in maintaining proteins in their native state in the cold and thermally stable Antarctic environments. Our transcriptome analysis of Antarctic notothenioid fish provides an archive for future studies in molecular mechanisms of fundamental genetic questions, and can be used in evolution studies comparing other fish.

De novo sequencing and analysis of root transcriptome using 454 pyrosequencing to discover putative genes associated with drought tolerance in Ammopiptanthus mongolicus

Background

De novo assembly of transcript sequences produced by next-generation sequencing technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Ammopiptanthus mongolicus, a super-xerophytic broadleaf evergreen wood, is an ecologically important foundation species in desert ecosystems and exhibits substantial drought tolerance in Mid-Asia desert. Root plays an important role in water absorption of plant. There are insufficient transcriptomic and genomic data in public databases for understanding of the molecular mechanism underlying the drought tolerance of A. mongolicus. Thus, high throughput transcriptome sequencing from A. mongolicus root is helpful to generate a large amount of transcript sequences for gene discovery and molecular marker development.

Results

A total of 672,002 sequencing reads were obtained from a 454 GS XLR70 Titanium pyrosequencer with a mean length of 279 bp. These reads were assembled into 29,056 unique sequences including 15,173 contigs and 13,883 singlets. In our assembled sequences, 1,827 potential simple sequence repeats (SSR) molecular markers were discovered. Based on sequence similarity with known plant proteins, the assembled sequences represent approximately 9,771 proteins in PlantGDB. Based on the Gene ontology (GO) analysis, hundreds of drought stress-related genes were found. We further analyzed the gene expression profiles of 27 putative genes involved in drought tolerance using quantitative real-time PCR (qRT-PCR) assay.

Conclusions

Our sequence collection represents a major transcriptomic resource for A. mongolicus, and the large number of genetic markers predicted should contribute to future research in Ammopiptanthus genus. The potential drought stress related transcripts identified in this study provide a good start for further investigation into the drought adaptation in Ammopiptanthus.

The barley EST DNA Replication and Repair Database (bEST-DRRD) as a tool for the identification of the genes involved in DNA replication and repair

BACKGROUND

The high level of conservation of genes that regulate DNA replication and repair indicates that they may serve as a source of information on the origin and evolution of the species and makes them a reliable system for the identification of cross-species homologs. Studies that had been conducted to date shed light on the processes of DNA replication and repair in bacteria, yeast and mammals. However, there is still much to be learned about the process of DNA damage repair in plants.

DESCRIPTION

These studies, which were conducted mainly using bioinformatics tools, enabled the list of genes that participate in various pathways of DNA repair in Arabidopsis thaliana (L.) Heynh to be outlined; however, information regarding these mechanisms in crop plants is still very limited. A similar, functional approach is particularly difficult for a species whose complete genomic sequences are still unavailable. One of the solutions is to apply ESTs (Expressed Sequence Tags) as the basis for gene identification. For the construction of the barley EST DNA Replication and Repair Database (bEST-DRRD), presented here, the Arabidopsis nucleotide and protein sequences involved in DNA replication and repair were used to browse for and retrieve the deposited sequences, derived from four barley (Hordeum vulgare L.) sequence databases, including the "Barley Genome version 0.05" database (encompassing ca. 90% of barley coding sequences) and from two databases covering the complete genomes of two monocot models: Oryza sativa L. and Brachypodium distachyon L. in order to identify homologous genes. Sequences of the categorised Arabidopsis queries are used for browsing the repositories, which are located on the ViroBLAST platform. The bEST-DRRD is currently used in our project during the identification and validation of the barley genes involved in DNA repair.

CONCLUSIONS

The presented database provides information about the Arabidopsis genes involved in DNA replication and repair, their expression patterns and models of protein interactions. It was designed and established to provide an open-access tool for the identification of monocot homologs of known Arabidopsis genes that are responsible for DNA-related processes. The barley genes identified in the project are currently being analysed to validate their function.

Pyrosequencing the Bemisia tabaci transcriptome reveals a highly diverse bacterial community and a robust system for insecticide resistance

Background

Bemisia tabaci (Gennadius) is a phloem-feeding insect poised to become one of the major insect pests in open field and greenhouse production systems throughout the world. The high level of resistance to insecticides is a main factor that hinders continued use of insecticides for suppression of B. tabaci. Despite its prevalence, little is known about B. tabaci at the genome level. To fill this gap, an invasive B. tabaci B biotype was subjected to pyrosequencing-based transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes.

Methodology and Principal Findings

Using Roche 454 pyrosequencing, 857,205 reads containing approximately 340 megabases were obtained from the B. tabaci transcriptome. De novo assembly generated 178,669 unigenes including 30,980 from insects, 17,881 from bacteria, and 129,808 from the nohit. A total of 50,835 (28.45%) unigenes showed similarity to the non-redundant database in GenBank with a cut-off E-value of 10–5. Among them, 40,611 unigenes were assigned to one or more GO terms and 6,917 unigenes were assigned to 288 known pathways. De novo metatranscriptome analysis revealed highly diverse bacterial symbionts in B. tabaci, and demonstrated the host-symbiont cooperation in amino acid production. In-depth transcriptome analysis indentified putative molecular markers, and genes potentially involved in insecticide resistance and nutrient digestion. The utility of this transcriptome was validated by a thiamethoxam resistance study, in which annotated cytochrome P450 genes were significantly overexpressed in the resistant B. tabaci in comparison to its susceptible counterparts.

Conclusions

This transcriptome/metatranscriptome analysis sheds light on the molecular understanding of symbiosis and insecticide resistance in an agriculturally important phloem-feeding insect pest, and lays the foundation for future functional genomics research of the B. tabaci complex. Moreover, current pyrosequencing effort greatly enriched the existing whitefly EST database, and makes RNAseq a viable option for future genomic analysis.

Insect omics research coming of age1This review is part of a virtual symposium on recent advances in understanding a variety of complex regulatory processes in insect physiology and endocrinology, including development, metabolism, cold hardiness, food intake and digestion, and diuresis, through the use of omics technologies in the postgenomic era

As more and more insect genomes are fully sequenced and annotated, omics technologies, including transcriptomic, proteomic, peptidomics, and metobolomic profiling, as well as bioinformatics, can be used to exploit this huge amount of sequence information for the study of different biological aspects of insect model organisms. Omics experiments are an elegant way to deliver candidate genes, the function of which can be further explored by genetic tools for functional inactivation or overexpression of the genes of interest. Such tools include mainly RNA interference and are currently being developed in diverse insect species. In this manuscript, we have reviewed how omics technologies were integrated and applied in insect biology.

EST and transcriptome analysis of cephalochordate amphioxus--past, present and future

The cephalochordates, commonly known as amphioxus or lancelets, are now considered the most basal chordate group, and the studies of these organisms therefore offer important insights into various levels of evolutionary biology. In the past two decades, the investigation of amphioxus developmental biology has provided key knowledge for understanding the basic patterning mechanisms of chordates. Comparative genome studies of vertebrates and amphioxus have uncovered clear evidence supporting the hypothesis of two-round whole-genome duplication thought to have occurred early in vertebrate evolution and have shed light on the evolution of morphological novelties in the complex vertebrate body plan. Complementary to the amphioxus genome-sequencing project, a large collection of expressed sequence tags (ESTs) has been generated for amphioxus in recent years; this valuable collection represents a rich resource for gene discovery, expression profiling and molecular developmental studies in the amphioxus model. Here, we review previous EST analyses and available cDNA resources in amphioxus and discuss their value for use in evolutionary and developmental studies. We also discuss the potential advantages of applying high-throughput, next-generation sequencing (NGS) technologies to the field of amphioxus research.

Analysis of the transcriptome of Panax notoginseng root uncovers putative triterpene saponin-biosynthetic genes and genetic markers

Background

Panax notoginseng (Burk) F.H. Chen is important medicinal plant of the Araliacease family. Triterpene saponins are the bioactive constituents in P. notoginseng. However, available genomic information regarding this plant is limited. Moreover, details of triterpene saponin biosynthesis in the Panax species are largely unknown.

Results

Using the 454 pyrosequencing technology, a one-quarter GS FLX titanium run resulted in 188,185 reads with an average length of 410 bases for P. notoginseng root. These reads were processed and assembled by 454 GS De Novo Assembler software into 30,852 unique sequences. A total of 70.2% of unique sequences were annotated by Basic Local Alignment Search Tool (BLAST) similarity searches against public sequence databases. The Kyoto Encyclopedia of Genes and Genomes (KEGG) assignment discovered 41 unique sequences representing 11 genes involved in triterpene saponin backbone biosynthesis in the 454-EST dataset. In particular, the transcript encoding dammarenediol synthase (DS), which is the first committed enzyme in the biosynthetic pathway of major triterpene saponins, is highly expressed in the root of four-year-old P. notoginseng. It is worth emphasizing that the candidate cytochrome P450 (Pn02132 and Pn00158) and UDP-glycosyltransferase (Pn00082) gene most likely to be involved in hydroxylation or glycosylation of aglycones for triterpene saponin biosynthesis were discovered from 174 cytochrome P450s and 242 glycosyltransferases by phylogenetic analysis, respectively. Putative transcription factors were detected in 906 unique sequences, including Myb, homeobox, WRKY, basic helix-loop-helix (bHLH), and other family proteins. Additionally, a total of 2,772 simple sequence repeat (SSR) were identified from 2,361 unique sequences, of which, di-nucleotide motifs were the most abundant motif.

Conclusion

This study is the first to present a large-scale EST dataset for P. notoginseng root acquired by next-generation sequencing (NGS) technology. The candidate genes involved in triterpene saponin biosynthesis, including the putative CYP450s and UGTs, were obtained in this study. Additionally, the identification of SSRs provided plenty of genetic makers for molecular breeding and genetics applications in this species. These data will provide information on gene discovery, transcriptional regulation and marker-assisted selection for P. notoginseng. The dataset establishes an important foundation for the study with the purpose of ensuring adequate drug resources for this species.

Functional genomics of a living fossil tree, Ginkgo, based on next-generation sequencing technology

Ginkgo biloba is monotypic species native to China and has old, dioecious, medicinally important characteristics. The functional genes related to these characteristics have not been effectively explored due to a limited number of expressed sequence tags (ESTs) from Ginkgo. To discover novel functional genes efficiently and to understand the development of a living fossil tree, Ginkgo, we used massive parallel pyrosequencing on the Roche 454 GS FLX Titanium platform to generate 64 057 ESTs. The ESTs combined with the 21 590 Ginkgo ESTs in genbank were assembled into 22 304 unique putative transcripts, in which 13 922 novel unique putative transcripts were identified by 454 sequencing. After being assigned to putative functions with Gene Ontology terms, a detailed view of the Ginkgo biological systems was displayed, including characterization of unique putative transcripts with homology to known key enzymes and transcription factors involved in ginkgolide/bilobalide and flavonoid biosynthetic pathways, as well as unique putative transcripts related to development, response to disease and defence. The fact that three full-length Ginkgo genes encoding key enzymes were found and cloned, suggests that high-throughput sequencing technology is superior to traditional gene-by-gene approach in discovery of genes. Additionally, a total of 204 simple sequence repeat motifs were detected. Our study not only lays the foundations for transcriptome-led studies in biosynthetic mechanisms, but also contributes significantly to the understanding of functional genomics and development in non-model plants.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

454 EST analysis detects genes putatively involved in ginsenoside biosynthesis in Panax ginseng

Panax ginseng C.A. Meyer is one of the most highly valued medicinal plants in the world. To analyze the transcriptome of P. ginseng and discover the genes involved in ginsenoside biosynthesis, cDNAs derived from the total RNA of 11-year-old, wood-grown P. ginseng roots were analyzed by 454 sequencing. A total of 217,529 high quality reads (expressed sequence tags, ESTs), with an average length of 409 bases, were generated from a one-quarter run to yield 31,741 unique sequences. The majority (20,198; 63.6%) of the unique sequences were annotated using BLAST similarity searches. A total of 16,810 and 16,577 unique sequences were assigned to functional classifications and biochemical pathways based on Gene Ontology analysis and the Kyoto Encyclopedia of Genes and Genomes assignment, respectively. Nine genes involved in the biosynthesis of ginsenoside skeletons and many candidate genes putatively responsible for modification of the skeletons, including 133 cytochrome P450s and 235 glycosyltransferases, were identified. From these candidates, six transcripts encoding UDP-glycosyltransferases that were most likely to be involved in ginsenoside biosynthesis were selected. These results open a new avenue by which to explore and exploit biosynthetic and biochemical properties that may lead to drug improvement. These 454 ESTs will provide the foundation for further functional genomic research into the traditional herb P. ginseng or its closely related species.

The mining of toxin-like polypeptides from EST database by single residue distribution analysis

Background

Novel high throughput sequencing technologies require permanent development of bioinformatics data processing methods. Among them, rapid and reliable identification of encoded proteins plays a pivotal role. To search for particular protein families, the amino acid sequence motifs suitable for selective screening of nucleotide sequence databases may be used. In this work, we suggest a novel method for simplified representation of protein amino acid sequences named Single Residue Distribution Analysis, which is applicable both for homology search and database screening.

Results

Using the procedure developed, a search for amino acid sequence motifs in sea anemone polypeptides was performed, and 14 different motifs with broad and low specificity were discriminated. The adequacy of motifs for mining toxin-like sequences was confirmed by their ability to identify 100% toxin-like anemone polypeptides in the reference polypeptide database. The employment of novel motifs for the search of polypeptide toxins in Anemonia viridis EST dataset allowed us to identify 89 putative toxin precursors. The translated and modified ESTs were scanned using a special algorithm. In addition to direct comparison with the motifs developed, the putative signal peptides were predicted and homology with known structures was examined.

Conclusions

The suggested method may be used to retrieve structures of interest from the EST databases using simple amino acid sequence motifs as templates. The efficiency of the procedure for directed search of polypeptides is higher than that of most currently used methods. Analysis of 39939 ESTs of sea anemone Anemonia viridis resulted in identification of five protein precursors of earlier described toxins, discovery of 43 novel polypeptide toxins, and prediction of 39 putative polypeptide toxin sequences. In addition, two precursors of novel peptides presumably displaying neuronal function were disclosed.

Tangible benefits of the aphid Acyrthosiphon pisum genome sequencing for aphid proteomics: Enhancements in protein identification and data validation for homology-based proteomics

Homology-driven proteomics promises to reveal functional biology in insects with sparse genome sequence information. A proteomics study comparing plant virus transmission competent and refractive genotypes of the aphid Schizaphis graminum isolated numerous candidate proteins involved in virus transmission, but limited genome sequence information hampered their identification. The complete genome of the pea aphid, Acyrthosiphon pisum, released in 2008, enabled us to double the number of protein identifications beyond what was possible using available EST libraries and other insect sequences. This was concomitant with a dramatic increase of the number of MS and MS/MS peptide spectra matching the genome-derived protein sequence. LC-MS/MS proved to be the most robust method of peptide detection. Cross-matching spectral data to multiple EST sequences and error tolerant searching to identify amino acid substitutions enhanced the percent coverage of the Schizaphis graminum proteins. 2-D electrophoresis provided the protein pI and MW which enabled the refinement of the candidate protein selection and provided a measure of protein abundance when coupled to the spectral data. Thus, the homology-based proteomics pipeline for insects should include efforts to maximize the number of peptide matches to the protein to increase certainty in protein identification and relative protein abundance.

TESTLoc: protein subcellular localization prediction from EST data

Background

The eukaryotic cell has an intricate architecture with compartments and substructures dedicated to particular biological processes. Knowing the subcellular location of proteins not only indicates how bio-processes are organized in different cellular compartments, but also contributes to unravelling the function of individual proteins. Computational localization prediction is possible based on sequence information alone, and has been successfully applied to proteins from virtually all subcellular compartments and all domains of life. However, we realized that current prediction tools do not perform well on partial protein sequences such as those inferred from Expressed Sequence Tag (EST) data, limiting the exploitation of the large and taxonomically most comprehensive body of sequence information from eukaryotes.

Results

We developed a new predictor, TESTLoc, suited for subcellular localization prediction of proteins based on their partial sequence conceptually translated from ESTs (EST-peptides). Support Vector Machine (SVM) is used as computational method and EST-peptides are represented by different features such as amino acid composition and physicochemical properties. When TESTLoc was applied to the most challenging test case (plant data), it yielded high accuracy (~85%).

Conclusions

TESTLoc is a localization prediction tool tailored for EST data. It provides a variety of models for the users to choose from, and is available for download at http://megasun.bch.umontreal.ca/~shenyq/TESTLoc/TESTLoc.html