Sequencing Medicago truncatula expressed sequenced tags using 454 Life Sciences technology

More than Three Decades of Bm86: What We Know and Where to Go

Tick and tick-borne disease control have been a serious research focus for many decades. In a global climate of increasing acaricide resistance, host immunity against tick infestation has become a much-needed complementary strategy to common chemical control. From the earliest acquired resistance studies in small animal models to proof of concept in large production animals, it was the isolation, characterization, and final recombinant protein production of the midgut antigen Bm86 from the Australian cattle tick strain of Rhipicephalus (Boophilus) microplus (later reinstated as R. (B.) australis) that established tick subunit vaccines as a viable alternative in tick and tick-borne disease control. In the past 37 years, this antigen has spawned numerous tick subunit vaccines (either Bm86-based or novel), and though we are still describing its molecular structure and function, this antigen remains the gold standard for all tick vaccines. In this paper, advances in tick vaccine development over the past three decades are discussed alongside the development of biotechnology, where existing gaps and future directives in the field are highlighted.

The Prospects of gene introgression from crop wild relatives into cultivated lentil for climate change mitigation

Crop wild relatives (CWRs), landraces and exotic germplasm are important sources of genetic variability, alien alleles, and useful crop traits that can help mitigate a plethora of abiotic and biotic stresses and crop yield reduction arising due to global climatic changes. In the pulse crop genus Lens, the cultivated varieties have a narrow genetic base due to recurrent selections, genetic bottleneck and linkage drag. The collection and characterization of wild Lens germplasm resources have offered new avenues for the genetic improvement and development of stress-tolerant, climate-resilient lentil varieties with sustainable yield gains to meet future food and nutritional requirements. Most of the lentil breeding traits such as high-yield, adaptation to abiotic stresses and resistance to diseases are quantitative and require the identification of quantitative trait loci (QTLs) for marker assisted selection and breeding. Advances in genetic diversity studies, genome mapping and advanced high-throughput sequencing technologies have helped identify many stress-responsive adaptive genes, quantitative trait loci (QTLs) and other useful crop traits in the CWRs. The recent integration of genomics technologies with plant breeding has resulted in the generation of dense genomic linkage maps, massive global genotyping, large transcriptomic datasets, single nucleotide polymorphisms (SNPs), expressed sequence tags (ESTs) that have advanced lentil genomic research substantially and allowed for the identification of QTLs for marker-assisted selection (MAS) and breeding. Assembly of lentil and its wild species genomes (~4Gbp) opens up newer possibilities for understanding genomic architecture and evolution of this important legume crop. This review highlights the recent strides in the characterization of wild genetic resources for useful alleles, development of high-density genetic maps, high-resolution QTL mapping, genome-wide studies, MAS, genomic selections, new databases and genome assemblies in traditionally bred genus Lens for future crop improvement amidst the impending global climate change.

Developing New Tools to Fight Human Pathogens: A Journey through the Advances in RNA Technologies

A long scientific journey has led to prominent technological advances in the RNA field, and several new types of molecules have been discovered, from non-coding RNAs (ncRNAs) to riboswitches, small interfering RNAs (siRNAs) and CRISPR systems. Such findings, together with the recognition of the advantages of RNA in terms of its functional performance, have attracted the attention of synthetic biologists to create potent RNA-based tools for biotechnological and medical applications. In this review, we have gathered the knowledge on the connection between RNA metabolism and pathogenesis in Gram-positive and Gram-negative bacteria. We further discuss how RNA techniques have contributed to the building of this knowledge and the development of new tools in synthetic biology for the diagnosis and treatment of diseases caused by pathogenic microorganisms. Infectious diseases are still a world-leading cause of death and morbidity, and RNA-based therapeutics have arisen as an alternative way to achieve success. There are still obstacles to overcome in its application, but much progress has been made in a fast and effective manner, paving the way for the solid establishment of RNA-based therapies in the future.

Selection and validation of reference genes for measuring gene expression in Piper species at different life stages using RT-qPCR analysis

The secondary metabolism of Piper species is known to produce a myriad of natural products from various biosynthetic pathways which, represent a rich source of previously uncharacterized chemical compounds. The determination of gene expression profiles in multiple tissue/organ samples could provide valuable clues towards understanding the potential biological functions of chemical changes in these plants. Studies on gene expression by RT-qPCR require particularly careful selection of suitable reference genes as a control for normalization. Here, we provide a study for the identification of reliable reference genes in P. arboreum, P. gaudichaudianum, P. malacophyllum, and P. tuberculatum, at two different life stages: 2-month-old seedlings and adult plants. To do this, annotated sequences were recovered from transcriptome datasets of the above listed Piper spp. These sequences were subjected to expression analysis using RT-qPCR, followed by analysis using the geNorm and NormFinder algorithms. A set of five genes were identified showing stable expression: ACT7 (Actin-7), Cyclophilin (Peptidyl-prolyl cis-trans isomerase), EF1α (Elongation factor 1-alpha), RNABP (RNA-binding protein), and UBCE (Ubiquitin conjugating enzyme). The universality of these genes was then validated using two target genes, ADC (arginine decarboxylase) and SAMDC (S-adenosylmethionine decarboxylase), which are involved in the biosynthesis of polyamines. We showed that normalization genes varied according to Piper spp., and we provide a list of recommended pairs of the best combination for each species. This study provides the first set of suitable candidate genes for gene expression studies in the four Piper spp. assayed, and the findings will facilitate subsequent transcriptomic and functional gene research.

Protocol for RNA-seq Expression Analysis in Yeast

Genome-wide sequencing of RNA (RNA-seq) has become an inexpensive tool to gain key insights into cellular and disease mechanisms. Sample preparation and sequencing are streamlined and allow the acquisition of hundreds of gene expression profiles in a few days; however, in particular, data processing, curation, and analysis involve numerous steps that can be overwhelming to non-experts. Here, the sample preparation, sequencing, and data processing workflow for RNA-seq expression analysis in yeast is described. While this protocol covers only a small portion of the RNA-seq landscape, the principal workflow common to such experiments is described, allowing the reader to adapt the protocol where necessary. Graphic abstract: Basic workflow of RNA-seq expression analysis.

Non-Coding RNAs in Cancer Diagnosis and Therapy: Focus on Lung Cancer

Over the last several decades, clinical evaluation and treatment of lung cancers have largely improved with the classification of genetic drivers of the disease, such as EGFR, ALK, and ROS1. There are numerous regulatory factors that exert cellular control over key oncogenic pathways involved in lung cancers. In particular, non-coding RNAs (ncRNAs) have a diversity of regulatory roles in lung cancers such that they have been shown to be involved in inducing proliferation, suppressing apoptotic pathways, increasing metastatic potential of cancer cells, and acquiring drug resistance. The dysregulation of various ncRNAs in human cancers has prompted preclinical studies examining the therapeutic potential of restoring and/or inhibiting these ncRNAs. Furthermore, ncRNAs demonstrate tissue-specific expression in addition to high stability within biological fluids. This makes them excellent candidates as cancer biomarkers. This review aims to discuss the relevance of ncRNAs in cancer pathology, diagnosis, and therapy, with a focus on lung cancer.

Transcriptome analysis of Plantago major as a phytoremediator to identify some genes related to cypermethrin detoxification

Cypermethrin (CYP) is a toxic manmade chemical compound belonging to pyrethroid insecticides contaminating the environment. Plantago major (PM) has numerous excellent advantages like high biomass yield and great stress tolerance, which make it able to increase the efficacy of phytoremediation. So far, no study has directly or indirectly made a transcriptome analysis (RNA-seq) of PM under CYP stress. The aim of this study is to identify the genes in PM related to CYP detoxification (10 μg mL^-1) and compared with control. In this study, BGISEQ-500 high-throughput sequencing technology independently developed by BGI was used to sequence the transcriptome of P. major. Six libraries were constructed including (CK_1, CK_2, and CK_3) and (CYP_1, CYP_2, and CYP_3) were sequenced for transcripts involved in CYP detoxification. Our data showed that de novo assembly generated 138,806 unigenes with an average length of 1129 bp. Analyzing the annotation results of the KEGG database between the samples revealed 37,177 differentially expressed genes (DEGs), 18,062 down- and 19,115 upregulated under CYP treatment compared with control. A set of 107 genes of cytochrome P450 (Cyt P450), 43 genes of glutathione S-transferases (GST), 25 genes of glycosyltransferases (GTs), 113 genes from ABC transporters, 21 genes from multidrug and toxin efflux (MATE), 11 genes from oligopeptide transporter (OPT), and 3 genes of metallothioneins (MT) were upregulated notably. By using quantitative real-time PCR (qRT-PCR), the results of gene expression for 12 randomly selected DEGs were confirmed, showing the different patterns of response to CYP in PM tissues. Furthermore, the enzyme activity of Cyt P450 and GST in PM under CYP stress was significantly increased in roots and leaves than in control. This study introduces a clue to understand the metabolic pathways of plants used in phytoremediation by identifying the highly expressed genes related to phytoremediation which would be utilized to enhance pesticide detoxification and reduce pollution problem.

Optimized design of antisense oligomers for targeted rRNA depletion

RNA sequencing (RNA-seq) is extensively used to quantify gene expression transcriptome-wide. Although often paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA), many applications require alternate approaches to counteract the high proportion of ribosomal RNA (rRNA) in total RNA. Recently, digestion using RNaseH and antisense DNA oligomers tiling target rRNAs has emerged as an alternative to commercial rRNA depletion kits. Here, we present a streamlined, more economical RNaseH-mediated rRNA depletion with substantially lower up-front costs, using shorter antisense oligos only sparsely tiled along the target RNA in a 5-min digestion reaction. We introduce a novel Web tool, Oligo-ASST, that simplifies oligo design to target regions with optimal thermodynamic properties, and additionally can generate compact, common oligo pools that simultaneously target divergent RNAs, e.g. across different species. We demonstrate the efficacy of these strategies by generating rRNA-depletion oligos for Xenopus laevis and for zebrafish, which expresses two distinct versions of rRNAs during embryogenesis. The resulting RNA-seq libraries reduce rRNA to <5% of aligned reads, on par with poly(A) selection, and also reveal expression of many non-adenylated RNA species. Oligo-ASST is freely available at https://mtleelab.pitt.edu/oligo to design antisense oligos for any taxon or to target any abundant RNA for depletion.

Biotechnological Perspectives of Omics and Genetic Engineering Methods in Alfalfa

For several decades, researchers are working to develop improved major crops with better adaptability and tolerance to environmental stresses. Forage legumes have been widely spread in the world due to their great ecological and economic values. Abiotic and biotic stresses are main factors limiting legume production, however, alfalfa (Medicago sativa L.) shows relatively high level of tolerance to drought and salt stress. Efforts focused on alfalfa improvements have led to the release of cultivars with new traits of agronomic importance such as high yield, better stress tolerance or forage quality. Alfalfa has very high nutritional value due to its efficient symbiotic association with nitrogen-fixing bacteria, while deep root system can help to prevent soil water loss in dry lands. The use of modern biotechnology tools is challenging in alfalfa since full genome, unlike to its close relative barrel medic (Medicago truncatula Gaertn.), was not released yet. Identification, isolation, and improvement of genes involved in abiotic or biotic stress response significantly contributed to the progress of our understanding how crop plants cope with these environmental challenges. In this review, we provide an overview of the progress that has been made in high-throughput sequencing, characterization of genes for abiotic or biotic stress tolerance, gene editing, as well as proteomic and metabolomics techniques bearing biotechnological potential for alfalfa improvement.

Elucidation of Galactomannan Biosynthesis Pathway Genes through Transcriptome Sequencing of Seeds Collected at Different Developmental Stages of Commercially Important Indian Varieties of Cluster Bean (Cyamopsis tetragonoloba L.)

Cyamopsis tetragonoloba (L) endosperm predominantly contains guar gum a polysaccharide, which has tremendous industrial applications in food, textile, paper, oil drilling and water treatment. In order to understand the genes controlling galactomannan biosynthesis, mRNA was isolated from seeds collected at different developmental stages; young pods, mature pods and young leaf from two guar varieties, HG365 and HG870 and subjected to Illumina sequencing. De novo assembly of fourteen individual read files from two varieties of guar representing seven developmental stages gave a total of 1,13,607 contigs with an N50 of 1,244 bases. Annotation of assemblies with GO mapping revealed three levels of distribution, namely, Biological Processes, Molecular Functions and Cellular Components. GO studies identified major genes involved in galactomannan biosynthesis: Cellulose synthase D1 (CS D1) and GAUT-like gene families. Among the polysaccharide biosynthetic process (GO:0000271) genes the transcript abundance for CS was found to be predominantly more in leaf samples, whereas, the transcript abundance for GAUT-like steadily increased from 65% to 90% and above from stage1 to stage5 indicating accumulation of galactomannan in developing seeds; and validated by qRT-PCR analysis. Galactomannan quantification by HPLC showed HG365 (12.98–20.66%) and HG870 (7.035–41.2%) gradually increasing from stage1 to stage 5 (10–50 DAA) and highest accumulation occurred in mature and dry seeds with 3.8 to 7.1 fold increase, respectively. This is the first report of transcriptome sequencing and complete profiling of guar seeds at different developmental stages, young pods, mature pods and young leaf material from two commercially important Indian varieties and elucidation of galactomannan biosynthesis pathway. It is envisaged that the data presented herein will be very useful for improvement of guar through biotechnological interventions in future.

Legume genetic resources and transcriptome dynamics under abiotic stress conditions

Grain legumes are an important source of nutrition and income for billions of consumers and farmers around the world. However, the low productivity of new legume varieties, due to the limited genetic diversity available for legume breeding programmes and poor policymaker support, combined with an increasingly unpredictable global climate is resulting in a large gap between current yields and the increasing demand for legumes as food. Hence, there is a need for novel approaches to develop new high-yielding legume cultivars that are able to cope with a range of environmental stressors. Next-generation technologies are providing the tools that could enable the more rapid and cost-effective genomic and transcriptomic studies for most major crops, allowing the identification of key functional and regulatory genes involved in abiotic stress resistance. In this review, we provide an overview of the recent achievements regarding abiotic stress resistance in a wide range of legume crops and highlight the transcriptomic and miRNA approaches that have been used. In addition, we critically evaluate the availability and importance of legume genetic resources with desirable abiotic stress resistance traits.

Computational functional genomics-based approaches in analgesic drug discovery and repurposing

Persistent pain is a major healthcare problem affecting a fifth of adults worldwide with still limited treatment options. The search for new analgesics increasingly includes the novel research area of functional genomics, which combines data derived from various processes related to DNA sequence, gene expression or protein function and uses advanced methods of data mining and knowledge discovery with the goal of understanding the relationship between the genome and the phenotype. Its use in drug discovery and repurposing for analgesic indications has so far been performed using knowledge discovery in gene function and drug target-related databases; next-generation sequencing; and functional proteomics-based approaches. Here, we discuss recent efforts in functional genomics-based approaches to analgesic drug discovery and repurposing and highlight the potential of computational functional genomics in this field including a demonstration of the workflow using a novel R library ‘dbtORA’.

De Novo Transcriptome Sequencing and the Hypothetical Cold Response Mode of Saussurea involucrata in Extreme Cold Environments

Saussurea involucrata grows in high mountain areas covered by snow throughout the year. The temperature of this habitat can change drastically in one day. To gain a better understanding of the cold response signaling pathways and molecular metabolic reactions involved in cold stress tolerance, genome-wide transcriptional analyses were performed using RNA-Seq technologies. A total of 199,758 transcripts were assembled, producing 138,540 unigenes with 46.8 Gb clean data. Overall, 184,416 (92.32%) transcripts were successfully annotated. The 365 transcription factors identified (292 unigenes) belonged to 49 transcription factor families associated with cold stress responses. A total of 343 transcripts on the signal transduction (132 upregulated and 212 downregulated in at least any one of the conditions) were strongly affected by cold temperature, such as the CBL-interacting serine/threonine-protein kinase (CIPKs), receptor-like protein kinases, and protein kinases. The circadian rhythm pathway was activated by cold adaptation, which was necessary to endure the severe temperature changes within a day. There were 346 differentially expressed genes (DEGs) related to transport, of which 138 were upregulated and 22 were downregulated in at least any one of the conditions. Under cold stress conditions, transcriptional regulation, molecular transport, and signal transduction were involved in the adaptation to low temperature in S. involucrata. These findings contribute to our understanding of the adaptation of plants to harsh environments and the survival traits of S. involucrata. In addition, the present study provides insight into the molecular mechanisms of chilling and freezing tolerance.

De novo transcriptome sequencing and comparative analysis to discover genes related to floral development in Cymbidium faberi Rolfe

Cymbidium faberi is a traditional orchid flower in China that is highly appreciated for its fragrant aroma from its zygomorphic flowers. One bottleneck of the commercial production of C. faberi is the long vegetative growth phase of the orchid and the difficulty of the regulation of its flowering time. Moreover, its flower size, shape and color are often targeting traits for orchid breeders. Understanding the molecular mechanisms of floral development in C. faberi will ultimately benefit the genetic improvement of this orchid plant. The goal of this study is to identify potential genes and regulatory networks related to the floral development in C. faberi by using transcriptome sequencing, de novo assembly and computational analyses. The vegetative and flower buds of C. faberi were sampled for such comparisons. The RNA-seq yielded about 189,300 contigs that were assembled into 172,959 unigenes. Furthermore, a total of 13,484 differentially expressed unigenes (DEGs) were identified between the vegetative and flower buds. There were 7683 down-regulated and 5801 up-regulated DEGs in the flower buds compared to those in the vegetative buds, among which 3430 and 6556 DEGs were specifically enriched in the flower or vegetative buds, respectively. A total of 173 DEGs orthologous to known genes associated with the floral organ development, floral symmetry and flowering time were identified, including 12 TCP transcription factors, 34 MADS-box genes and 28 flowering time related genes. Furthermore, expression levels of ten genes potentially involved in floral development and flowering time were verified by quantitative real-time PCR. The identified DEGs will facilitate the functional genetic studies for further understanding the flower development of C. faberi.

Dissecting the Root Nodule Transcriptome of Chickpea (Cicer arietinum L.)

A hallmark trait of chickpea (Cicer arietinum L.), like other legumes, is the capability to convert atmospheric nitrogen (N2) into ammonia (NH3) in symbiotic association with Mesorhizobium ciceri. However, the complexity of molecular networks associated with the dynamics of nodule development in chickpea need to be analyzed in depth. Hence, in order to gain insights into the chickpea nodule development, the transcriptomes of nodules at early, middle and late stages of development were sequenced using the Roche 454 platform. This generated 490.84 Mb sequence data comprising 1,360,251 reads which were assembled into 83,405 unigenes. Transcripts were annotated using Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways analysis. Differential expression analysis revealed that a total of 3760 transcripts were differentially expressed in at least one of three stages, whereas 935, 117 and 2707 transcripts were found to be differentially expressed in the early, middle and late stages of nodule development respectively. MapMan analysis revealed enrichment of metabolic pathways such as transport, protein synthesis, signaling and carbohydrate metabolism during root nodulation. Transcription factors were predicted and analyzed for their differential expression during nodule development. Putative nodule specific transcripts were identified and enriched for GO categories using BiNGO which revealed many categories to be enriched during nodule development, including transcription regulators and transporters. Further, the assembled transcriptome was also used to mine for genic SSR markers. In conclusion, this study will help in enriching the transcriptomic resources implicated in understanding of root nodulation events in chickpea.

Emerging Genomic Tools for Legume Breeding: Current Status and Future Prospects

Legumes play a vital role in ensuring global nutritional food security and improving soil quality through nitrogen fixation. Accelerated higher genetic gains is required to meet the demand of ever increasing global population. In recent years, speedy developments have been witnessed in legume genomics due to advancements in next-generation sequencing (NGS) and high-throughput genotyping technologies. Reference genome sequences for many legume crops have been reported in the last 5 years. The availability of the draft genome sequences and re-sequencing of elite genotypes for several important legume crops have made it possible to identify structural variations at large scale. Availability of large-scale genomic resources and low-cost and high-throughput genotyping technologies are enhancing the efficiency and resolution of genetic mapping and marker-trait association studies. Most importantly, deployment of molecular breeding approaches has resulted in development of improved lines in some legume crops such as chickpea and groundnut. In order to support genomics-driven crop improvement at a fast pace, the deployment of breeder-friendly genomics and decision support tools seems appear to be critical in breeding programs in developing countries. This review provides an overview of emerging genomics and informatics tools/approaches that will be the key driving force for accelerating genomics-assisted breeding and ultimately ensuring nutritional and food security in developing countries.

Discovering New Biology through Sequencing of RNA

Sequencing of RNA (RNA-Seq) was invented approximately 1 decade ago and has since revolutionized biological research. This update provides a brief historic perspective on the development of RNA-Seq and then focuses on the application of RNA-Seq in qualitative and quantitative analyses of transcriptomes. Particular emphasis is given to aspects of data analysis. Since the wet-lab and data analysis aspects of RNA-Seq are still rapidly evolving and novel applications are continuously reported, a printed review will be rapidly outdated and can only serve to provide some examples and general guidelines for planning and conducting RNA-Seq studies. Hence, selected references to frequently update online resources are given.

Role of bioinformatics in establishing microRNAs as modulators of abiotic stress responses: the new revolution

microRNAs (miRs) are a class of 21-24 nucleotide long non-coding RNAs responsible for regulating the expression of associated genes mainly by cleavage or translational inhibition of the target transcripts. With this characteristic of silencing, miRs act as an important component in regulation of plant responses in various stress conditions. In recent years, with drastic change in environmental and soil conditions different type of stresses have emerged as a major challenge for plants growth and productivity. The identification and profiling of miRs has itself been a challenge for research workers given their small size and large number of many probable sequences in the genome. Application of computational approaches has expedited the process of identification of miRs and their expression profiling in different conditions. The development of High-Throughput Sequencing (HTS) techniques has facilitated to gain access to the global profiles of the miRs for understanding their mode of action in plants. Introduction of various bioinformatics databases and tools have revolutionized the study of miRs and other small RNAs. This review focuses the role of bioinformatics approaches in the identification and study of the regulatory roles of plant miRs in the adaptive response to stresses.

De novo transcriptome analysis of Medicago falcata reveals novel insights about the mechanisms underlying abiotic stress-responsive pathway

Background

The entire world is facing a deteriorating environment. Understanding the mechanisms underlying plant responses to external abiotic stresses is important for breeding stress-tolerant crops and herbages. Phytohormones play critical regulatory roles in plants in the response to external and internal cues to regulate growth and development. Medicago falcata is one of the stress-tolerant candidate leguminous species and is able to fix atmospheric nitrogen. This ability allows leguminous plants to grow in nitrogen deficient soils.

Methods

We performed Illumina sequencing of cDNA prepared from abiotic stress treated M. falcata. Sequencedreads were assembled to provide a transcriptome resource. Transcripts were annotated using BLASTsearches against the NCBI non-redundant database and gene ontology definitions were assigned. Acomparison among the three abiotic stress treated samples was carried out. The expression of transcriptswas confirmed with qRT-PCR.

Results

We present an abiotic stress-responsive M. falcata transcriptome using next-generation sequencing data from samples grown under standard, dehydration, high salinity, and cold conditions. We combined reads from all samples and de novo assembled 98,515 transcripts to build the M. falcata gene index. A comprehensive analysis of the transcriptome revealed abiotic stress-responsive mechanisms underlying the metabolism and core signalling components of major phytohormones. We identified nod factor signalling pathways during early symbiotic nodulation that are modified by abiotic stresses. Additionally, a global comparison of homology between the M. falcata and M. truncatula transcriptomes, along with five other leguminous species, revealed a high level of global sequence conservation within the family.

Conclusions

M. falcata is shown to be a model candidate for studying abiotic stress-responsive mechanisms in legumes. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to abiotic stresses. Our data provides many gene candidates that might be used for herbage and crop breeding. Additionally, FalcataBase (http://bioinformatics.cau.edu.cn/falcata/) was built for storing these data.

Electronic supplementary material

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

Comparative transcriptome analysis of developmental stages of the Limonium bicolor leaf generates insights into salt gland differentiation

With the expansion of saline land worldwide, it is essential to establish a model halophyte to study the salt-tolerance mechanism. The salt glands in the epidermis of Limonium bicolor (a recretohalophyte) play a pivotal role in salt tolerance by secreting excess salts from tissues. Despite the importance of salt secretion, nothing is known about the molecular mechanisms of salt gland development. In this study, we applied RNA sequencing to profile early leaf development using five distinct developmental stages, which were quantified by successive collections of the first true leaves of L. bicolor with precise spatial and temporal resolution. Specific gene expression patterns were identified for each developmental stage. In particular, we found that genes controlling salt gland differentiation in L. bicolor may evolve in a trichome formation, which was also confirmed by mutants with increased salt gland densities. Genes involved in the special ultrastructure of salt glands were also elucidated. Twenty-six genes were proposed to participate in salt gland differentiation. Our dataset sheds light on the molecular processes underpinning salt gland development and thus represents a first step towards the bioengineering of active salt-secretion capacity in crops.

Novel and Stress Relevant EST Derived SSR Markers Developed and Validated in Peanut

With the aim to increase the number of functional markers in resource poor crop like cultivated peanut (Arachis hypogaea), large numbers of available expressed sequence tags (ESTs) in the public databases, were employed for the development of novel EST derived simple sequence repeat (SSR) markers. From 16424 unigenes, 2784 (16.95%) SSRs containing unigenes having 3373 SSR motifs were identified. Of these, 2027 (72.81%) sequences were annotated and 4124 gene ontology terms were assigned. Among different SSR motif-classes, tri-nucleotide repeats (33.86%) were the most abundant followed by di-nucleotide repeats (27.51%) while AG/CT (20.7%) and AAG/CTT (13.25%) were the most abundant repeat-motifs. A total of 2456 EST-SSR novel primer pairs were designed, of which 366 unigenes having relevance to various stresses and other functions, were PCR validated using a set of 11 diverse peanut genotypes. Of these, 340 (92.62%) primer pairs yielded clear and scorable PCR products and 39 (10.66%) primer pairs exhibited polymorphisms. Overall, the number of alleles per marker ranged from 1-12 with an average of 3.77 and the PIC ranged from 0.028 to 0.375 with an average of 0.325. The identified EST-SSRs not only enriched the existing molecular markers kitty, but would also facilitate the targeted research in marker-trait association for various stresses, inter-specific studies and genetic diversity analysis in peanut.

De novo sequencing analysis of the Rosa roxburghii fruit transcriptome reveals putative ascorbate biosynthetic genes and EST-SSR markers

Rosa roxburghii Tratt. is a well-known ornamental rose species native to China. In addition, the fruits of this species are valued for their nutritional and medicinal characteristics, especially their high ascorbic acid (AsA) levels. Nevertheless, AsA biosynthesis in R. roxburghii fruit has not been explored in detail because of a lack of genomic resources for this species. High-throughput transcriptomic sequencing generating large volumes of transcript sequence data can aid in gene discovery and molecular marker development. In this study, we generated more than 53 million clean reads using Illumina paired-end sequencing technology. De novo assembly yielded 106,590 unigenes, with an average length of 343 bp. On the basis of sequence similarity to known proteins, 9301 and 2393 unigenes were classified into Gene Ontology and Clusters of Orthologous Group categories, respectively. There were 7480 unigenes assigned to 124 pathways in the Kyoto Encyclopedia of Gene and Genome pathway database. BLASTx searches identified 498 unique putative transcripts encoding various transcription factors, some known to regulate fruit development. qRT-PCR validated the expressions of most of the genes encoding the main enzymes involved in ascorbate biosynthesis. In addition, 9131 potential simple sequence repeat (SSR) loci were identified among the unigenes. One hundred and two primer pairs were synthesized and 71 pairs produced an amplification product during initial screening. Among the amplified products, 30 were polymorphic in the 16 R. roxburghii germplasms tested. Our study was the first to produce a large volume of transcriptome data from R. roxburghii. The resulting sequence collection is a valuable resource for gene discovery and marker-assisted selective breeding in this rose species.

De novo transcriptome assembly and the identification of gene-associated single-nucleotide polymorphism markers in Asian and American ginseng roots

We performed de novo transcriptome sequencing for Panax ginseng and Panax quinquefolius accessions using the 454 GS FLX Titanium System and discovered annotation-based genome-wide single-nucleotide polymorphism (SNPs) using next-generation ginseng transcriptome data without reference genome sequence. The comprehensive transcriptome characterization with the mature roots of four ginseng accessions generated 297,170 reads for 'Cheonryang' cultivar, 305,673 reads for 'Yunpoong' cultivar, 311,861 reads for the G03080 breeding line, and 308,313 reads for P. quinquefolius. In transcriptome assembly, the lengths of the sample read were 156.42 Mb for 'Cheonryang', 161.95 Mb for 'Yunpoong', 165.07 Mb for G03080 breeding line, and 166.48 Mb for P. quinquefolius. A total of 97 primer pairs were designed with the homozygous SNP presented in all four accessions. SNP genotyping using high-resolution melting (HRM) analysis was performed to validate the putative SNP markers of 97 primer pairs. Out of the 73 primer pairs, 73 primer pairs amplified the target sequence and 34 primer pairs showed polymorphic melting curves in samples from 11 P. ginseng cultivars and one P. quinquefolius accession. Among the 34 polymorphic HRM-SNP primers, four primers were useful to distinguish ginseng cultivars. In the present study, we demonstrated that de novo transcriptome assembly and mapping analyses are useful in providing four HRM-SNP primer pairs that reliably show a high degree of polymorphism among ginseng cultivars.

Bovine serum albumin in saliva mediates grazing response in Leymus chinensis revealed by RNA sequencing

Background

Sheepgrass (Leymus chinensis) is an important perennial forage grass across the Eurasian Steppe and is adaptable to various environmental conditions, but little is known about its molecular mechanism responding to grazing and BSA deposition. Because it has a large genome, RNA sequencing is expensive and impractical except for the next-generation sequencing (NGS) technology.

Results

In this study, NGS technology was employed to characterize de novo the transcriptome of sheepgrass after defoliation and grazing treatments and to identify differentially expressed genes (DEGs) responding to grazing and BSA deposition. We assembled more than 47 M high-quality reads into 120,426 contigs from seven sequenced libraries. Based on the assembled transcriptome, we detected 2,002 DEGs responding to BSA deposition during grazing. Enrichment analysis of Gene ontology (GO), EuKaryotic Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways revealed that the effects of grazing and BSA deposition involved more apoptosis and cell oxidative changes compared to defoliation. Analysis of DNA fragments, cell oxidative factors and the lengths of leaf scars after grazing provided physiological and morphological evidence that BSA deposition during grazing alters the oxidative and apoptotic status of cells.

Conclusions

This research greatly enriches sheepgrass transcriptome resources and grazing-stress-related genes, helping us to better understand the molecular mechanism of grazing in sheepgrass. The grazing-stress-related genes and pathways will be a valuable resource for further gene-phenotype studies.

Electronic supplementary material

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

De novo assembly of the common bean transcriptome using short reads for the discovery of drought-responsive genes

The common bean (Phaseolus vulgaris L.) is one of the most important food legumes, far ahead of other legumes. The average grain yield of the common bean worldwide is much lower than its potential yields, primarily due to drought in the field. However, the gene network that mediates plant responses to drought stress remains largely unknown in this species. The major goals of our study are to identify a large scale of genes involved in drought stress using RNA-seq. First, we assembled 270 million high-quality trimmed reads into a non-redundant set of 62,828 unigenes, representing approximately 49 Mb of unique transcriptome sequences. Of these unigenes, 26,501 (42.2%) common bean unigenes had significant similarity with unigenes/predicted proteins from other legumes or sequenced plants. All unigenes were functionally annotated within the GO, COG and KEGG pathways. The strategy for de novo assembly of transcriptome data generated here will be useful in other legume plant transcriptome studies. Second, we identified 10,482 SSRs and 4,099 SNPs in transcripts. The large number of genetic markers provides a resource for gene discovery and development of functional molecular markers. Finally, we found differential expression genes (DEGs) between terminal drought and optimal irrigation treatments and between the two different genotypes Long 22-0579 (drought tolerant) and Naihua (drought sensitive). DEGs were confirmed by quantitative real-time PCR assays, which indicated that these genes are functionally associated with the drought-stress response. These resources will be helpful for basic and applied research for genome analysis and crop drought resistance improvement in the common bean.

Transcriptome analysis of the oil-rich tea plant, Camellia oleifera, reveals candidate genes related to lipid metabolism

Background

Rapidly driven by the need for developing sustainable sources of nutritionally important fatty acids and the rising concerns about environmental impacts after using fossil oil, oil-plants have received increasing awareness nowadays. As an important oil-rich plant in China, Camellia oleifera has played a vital role in providing nutritional applications, biofuel productions and chemical feedstocks. However, the lack of C. oleifera genome sequences and little genetic information have largely hampered the urgent needs for efficient utilization of the abundant germplasms towards modern breeding efforts of this woody oil-plant.

Results

Here, using the 454 GS-FLX sequencing platform, we generated approximately 600,000 RNA-Seq reads from four tissues of C. oleifera. These reads were trimmed and assembled into 104,842 non-redundant putative transcripts with a total length of ∼38.9 Mb, representing more than 218-fold of all the C. oleifera sequences currently deposited in the GenBank (as of March 2014). Based on the BLAST similarity searches, nearly 42.6% transcripts could be annotated with known genes, conserved domains, or Gene Ontology (GO) terms. Comparisons with the cultivated tea tree, C. sinensis, identified 3,022 pairs of orthologs, of which 211 exhibited the evidence under positive selection. Pathway analysis detected the majority of genes potentially related to lipid metabolism. Evolutionary analysis of omega-6 fatty acid desaturase (FAD2) genes among 20 oil-plants unexpectedly suggests that a parallel evolution may occur between C. oleifera and Olea oleifera. Additionally, more than 2,300 simple sequence repeats (SSRs) and 20,200 single-nucleotide polymorphisms (SNPs) were detected in the C. oleifera transcriptome.

Conclusions

The generated transcriptome represents a considerable increase in the number of sequences deposited in the public databases, providing an unprecedented opportunity to discover all related-genes associated with lipid metabolic pathway in C. oleifera. It will greatly enhance the generation of new varieties of C. oleifera with increased yields and high quality.

Contrasting soil microbial community functional structures in two major landscapes of the Tibetan alpine meadow

The grassland and shrubland are two major landscapes of the Tibetan alpine meadow, a region very sensitive to the impact of global warming and anthropogenic perturbation. Herein, we report a study showing that a majority of differences in soil microbial community functional structures, measured by a functional gene array named GeoChip 4.0, in two adjacent shrubland and grassland areas, were explainable by environmental properties, suggesting that the harsh environments in the alpine grassland rendered niche adaptation important. Furthermore, genes involved in labile carbon degradation were more abundant in the shrubland than those of the grassland but genes involved in recalcitrant carbon degradation were less abundant, which was conducive to long-term carbon storage and sequestration in the shrubland despite low soil organic carbon content. In addition, genes of anerobic nitrogen cycling processes such as denitrification and dissimilatory nitrogen reduction were more abundant, shifting soil nitrogen cycling toward ammonium biosynthesis and consequently leading to higher soil ammonium contents. We also noted higher abundances of stress genes responsive to nitrogen limitation and oxygen limitation, which might be attributed to low total nitrogen and higher water contents in the shrubland. Together, these results provide mechanistic knowledge about microbial linkages to soil carbon and nitrogen storage and potential consequences of vegetation shifts in the Tibetan alpine meadow.

Transcriptomic analysis of candidate osmoregulatory genes in the eastern oyster Crassostrea virginica

BACKGROUND

The eastern oyster, Crassostrea virginica, is a euryhaline species that can thrive across a wide range of salinities (5-35). As with all estuarine species, individual oysters must be able to regulate their osmotic balance in response to constant temporal variation in salinity. At the population level, recurrent viability selection may be an additional mechanism shaping adaptive osmoregulatory phenotypes at the margins of oyster salinity tolerance. To identify candidate genes for osmoregulation, we sequenced, assembled, and annotated the transcriptome of wild juvenile eastern oysters from 'high' and 'low' salinity regimes. Annotations and candidates were mostly based on the Pacific oyster (Crassostrea gigas) genome sequence so osmoregulatory relevance in C. virginica was explored by testing functional enrichment of genes showing spatially discrete patterns of expression and by quantifying coding sequence divergence.

RESULTS

The assembly of sequence reads and permissive clustering of potentially oversplit alleles resulted in 98,729 reftigs (contigs and singletons). Of these, 50,736 were annotated with 9,307 belonging to a set of candidate osmoregulatory genes identified from the C. gigas genome. A total of 218,777 SNPs (0.0185 SNPs/bp) were identified in annotated reftigs of C. virginica. Amino acid divergence between translations of C. virginica annotated reftigs and C. gigas coding sequence averaged 23.2 % with an average dN/dS ratio of 0.074, suggesting purifying selection on protein sequences. The high and low salinity source oysters each expressed a subset of genes unique to that group, and the functions for these annotated genes were consistent with known molecular mechanisms for osmotic regulation in molluscs.

CONCLUSIONS

Most of the osmoregulatory gene candidates experimentally identified in C. gigas are present in this C. virginica transcriptome. In general these congeners show coding sequence divergence too high to make the C. gigas genome a useful reference for C. virginica bioinformatics. However, strong purifying selection is characteristic of the osmoregulatory candidates so functional annotations are likely to correspond. An initial examination of C. virginica presence/absence expression patterns across the salinity gradient in a single estuary suggests that many of these candidates have expression patterns that co-vary with salinity, consistent with osmoregulatory function in C. virginica.

Transcriptomic analysis of Asiatic lily in the process of vernalization via RNA-seq

The new transcriptomes provided comprehensive sequence profiling data of transcriptomic variation during vernalization in Lily Asiatic Hybrids 'Tiny ghost'. A number of 52,277,184 sequencing raw reads totaling 5.11 Gbp of the chilling treatment (4 °C) sample and 39,466,176 sequencing raw reads totaling 3.85 Gbp of room temperature control (25 °C) sample were assembled de novo into 68,718 unigenes with a mean length of 984 bp, and a total of 33,208 (45.6%) unigenes were annotated by using public protein databases with a cut-off E value about 10(-5). There are 6,153 unigenes of which were assigned to specific metabolic pathways by the Kyoto encyclopedia of genes and genomes. Gene Ontology analysis of the annotated unigenes revealed that the majority of sequenced genes were associated with signal transduction mechanisms, posttranslational modification, protein turnover and chaperones. In addition, the genes expression levels were compared just after vernalization completion between the cold treatment and room temperature control. There are 68,116 unigenes were differentially expressed, and hierarchical clustering analysis arranged 7,301 significantly differentially-expressed unigenes into 56 groups. Six genes related to the vernalization were selected to confirm their expression levels by using quantitative real-time polymerase chain reaction. Furthermore, typical vernalization unigenes VRN1 and VRN2 were identified, and also some vernalization-associated unigenes, such as CBF, SOC, TaAGL, AP2, LEA, LIM et al. were also annotated in the present study. As for VRN1 and VRN2, their expressions were consistent with some previous related studies. Also, this was the first time the vernalization genes VRN1 and VRN2 were founded in lily. According to the results of the present studies, we predicted that they would play an important role during vernalization in Lily Asiatic Hybrids; these data provided the foundation for the future studies of metabolism during vernalization of Asiatic lily.

Global transcriptome analysis of developing chickpea (Cicer arietinum L.) seeds

Understanding developmental processes, especially in non-model crop plants, is extremely important in order to unravel unique mechanisms regulating development. Chickpea (C. arietinum L.) seeds are especially valued for their high carbohydrate and protein content. Therefore, in order to elucidate the mechanisms underlying seed development in chickpea, deep sequencing of transcriptomes from four developmental stages was undertaken. In this study, next generation sequencing platform was utilized to sequence the transcriptome of four distinct stages of seed development in chickpea. About 1.3 million reads were generated which were assembled into 51,099 unigenes by merging the de novo and reference assemblies. Functional annotation of the unigenes was carried out using the Uniprot, COG and KEGG databases. RPKM based digital expression analysis revealed specific gene activities at different stages of development which was validated using Real time PCR analysis. More than 90% of the unigenes were found to be expressed in at least one of the four seed tissues. DEGseq was used to determine differentially expressing genes which revealed that only 6.75% of the unigenes were differentially expressed at various stages. Homology based comparison revealed 17.5% of the unigenes to be putatively seed specific. Transcription factors were predicted based on HMM profiles built using TF sequences from five legume plants and analyzed for their differential expression during progression of seed development. Expression analysis of genes involved in biosynthesis of important secondary metabolites suggested that chickpea seeds can serve as a good source of antioxidants. Since transcriptomes are a valuable source of molecular markers like simple sequence repeats (SSRs), about 12,000 SSRs were mined in chickpea seed transcriptome and few of them were validated. In conclusion, this study will serve as a valuable resource for improved chickpea breeding.

Next generation sequencing analysis of human platelet PolyA+ mRNAs and rRNA-depleted total RNA

Background

Platelets are small anucleate cells circulating in the blood vessels where they play a key role in hemostasis and thrombosis. Here, we compared platelet RNA-Seq results obtained from polyA+ mRNA and rRNA-depleted total RNA.

Materials and Methods

We used purified, CD45 depleted, human blood platelets collected by apheresis from three male and one female healthy blood donors. The Illumina HiSeq 2000 platform was employed to sequence cDNA converted either from oligo(dT) isolated polyA+ RNA or from rRNA-depleted total RNA. The reads were aligned to the GRCh37 reference assembly with the TopHat/Cufflinks alignment package using Ensembl annotations. A de novo assembly of the platelet transcriptome using the Trinity software package and RSEM was also performed. The bioinformatic tools HTSeq and DESeq from Bioconductor were employed for further statistical analyses of read counts.

Results

Consistent with previous findings our data suggests that mitochondrially expressed genes comprise a substantial fraction of the platelet transcriptome. We also identified high transcript levels for protein coding genes related to the cytoskeleton function, chemokine signaling, cell adhesion, aggregation, as well as receptor interaction between cells. Certain transcripts were particularly abundant in platelets compared with other cell and tissue types represented by RNA-Seq data from the Illumina Human Body Map 2.0 project. Irrespective of the different library preparation and sequencing protocols, there was good agreement between samples from the 4 individuals. Eighteen differentially expressed genes were identified in the two sexes at 10% false discovery rate using DESeq.

Conclusion

The present data suggests that platelets may have a unique transcriptome profile characterized by a relative over-expression of mitochondrially encoded genes and also of genomic transcripts related to the cytoskeleton function, chemokine signaling and surface components compared with other cell and tissue types. The in vivo functional significance of the non-mitochondrial transcripts remains to be shown.

Functional genomics to study stress responses in crop legumes: progress and prospects

Legumes are important food crops worldwide, contributing to more than 33% of human dietary protein. The production of crop legumes is frequently impacted by abiotic and biotic stresses. It is therefore important to identify genes conferring resistance to biotic stresses and tolerance to abiotic stresses that can be used to both understand molecular mechanisms of plant response to the environment and to accelerate crop improvement. Recent advances in genomics offer a range of approaches such as the sequencing of genomes and transcriptomes, gene expression microarray as well as RNA-seq based gene expression profiling, and map-based cloning for the identification and isolation of biotic and abiotic stress-responsive genes in several crop legumes. These candidate stress associated genes should provide insights into the molecular mechanisms of stress tolerance and ultimately help to develop legume varieties with improved stress tolerance and productivity under adverse conditions. This review provides an overview on recent advances in the functional genomics of crop legumes that includes the discovery as well as validation of candidate genes.

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 sequencing of sunflower genome for SNP discovery using RAD (Restriction site Associated DNA) approach

Background

Application of Single Nucleotide Polymorphism (SNP) marker technology as a tool in sunflower breeding programs offers enormous potential to improve sunflower genetics, and facilitate faster release of sunflower hybrids to the market place. Through a National Sunflower Association (NSA) funded initiative, we report on the process of SNP discovery through reductive genome sequencing and local assembly of six diverse sunflower inbred lines that represent oil as well as confection types.

Results

A combination of Restriction site Associated DNA Sequencing (RAD-Seq) protocols and Illumina paired-end sequencing chemistry generated high quality 89.4 M paired end reads from the six lines which represent 5.3 GB of the sequencing data. Raw reads from the sunflower line, RHA 464 were assembled de novo to serve as a framework reference genome. About 15.2 Mb of sunflower genome distributed over 42,267 contigs were obtained upon assembly of RHA 464 sequencing data, the contig lengths ranged from 200 to 950 bp with an N₅₀ length of 393 bp. SNP calling was performed by aligning sequencing data from the six sunflower lines to the assembled reference RHA 464. On average, 1 SNP was located every 143 bp of the sunflower genome sequence. Based on several filtering criteria, a final set of 16,467 putative sequence variants with characteristics favorable for Illumina Infinium Genotyping Technology (IGT) were mined from the sequence data generated across six diverse sunflower lines.

Conclusion

Here we report the molecular and computational methodology involved in SNP development for a complex genome like sunflower lacking reference assembly, offering an attractive tool for molecular breeding purposes in sunflower.

Tracing the transcriptomic changes in synthetic Trigenomic allohexaploids of Brassica using an RNA-Seq approach

Polyploidization has played an important role in plant evolution and speciation, and newly formed allopolyploids have experienced rapid transcriptomic changes. Here, we compared the transcriptomic differences between a synthetic Brassica allohexaploid and its parents using a high-throughput RNA-Seq method. A total of 35,644,409 sequence reads were generated, and 32,642 genes were aligned from the data. Totals of 29,260, 29,060, and 29,697 genes were identified in Brassicarapa, Brassicacarinata, and Brassica allohexaploid, respectively. We compared 7,397 differentially expressed genes (DEGs) between Brassica hexaploid and its parents, as well as 2,545 nonadditive genes of Brassica hexaploid. We hypothesized that the higher ploidy level as well as secondary polyploidy might have influenced these changes. The majority of the 3,184 DEGs between Brassica hexaploid and its paternal parent, B. rapa, were involved in the biosynthesis of secondary metabolites, plant–pathogen interactions, photosynthesis, and circadian rhythm. Among the 2,233 DEGs between Brassica hexaploid and its maternal parent, B. carinata, several played roles in plant–pathogen interactions, plant hormone signal transduction, ribosomes, limonene and pinene degradation, photosynthesis, and biosynthesis of secondary metabolites. There were more significant differences in gene expression between the allohexaploid and its paternal parent than between it and its maternal parent, possibly partly because of cytoplasmic and maternal effects. Specific functional categories were enriched among the 2,545 nonadditive genes of Brassica hexaploid compared with the additive genes; the categories included response to stimulus, immune system process, cellular process, metabolic process, rhythmic process, and pigmentation. Many transcription factor genes, methyltransferases, and methylation genes showed differential expression between Brassica hexaploid and its parents. Our results demonstrate that the Brassica allohexaploid can generate extensive transcriptomic diversity compared with its parents. These changes may contribute to the normal growth and reproduction of allohexaploids.

Transcriptome analysis in sheepgrass (Leymus chinensis): a dominant perennial grass of the Eurasian Steppe

BACKGROUND

Sheepgrass [Leymus chinensis (Trin.) Tzvel.] is an important perennial forage grass across the Eurasian Steppe and is known for its adaptability to various environmental conditions. However, insufficient data resources in public databases for sheepgrass limited our understanding of the mechanism of environmental adaptations, gene discovery and molecular marker development.

RESULTS

The transcriptome of sheepgrass was sequenced using Roche 454 pyrosequencing technology. We assembled 952,328 high-quality reads into 87,214 unigenes, including 32,416 contigs and 54,798 singletons. There were 15,450 contigs over 500 bp in length. BLAST searches of our database against Swiss-Prot and NCBI non-redundant protein sequences (nr) databases resulted in the annotation of 54,584 (62.6%) of the unigenes. Gene Ontology (GO) analysis assigned 89,129 GO term annotations for 17,463 unigenes. We identified 11,675 core Poaceae-specific and 12,811 putative sheepgrass-specific unigenes by BLAST searches against all plant genome and transcriptome databases. A total of 2,979 specific freezing-responsive unigenes were found from this RNAseq dataset. We identified 3,818 EST-SSRs in 3,597 unigenes, and some SSRs contained unigenes that were also candidates for freezing-response genes. Characterizations of nucleotide repeats and dominant motifs of SSRs in sheepgrass were also performed. Similarity and phylogenetic analysis indicated that sheepgrass is closely related to barley and wheat.

CONCLUSIONS

This research has greatly enriched sheepgrass transcriptome resources. The identified stress-related genes will help us to decipher the genetic basis of the environmental and ecological adaptations of this species and will be used to improve wheat and barley crops through hybridization or genetic transformation. The EST-SSRs reported here will be a valuable resource for future gene-phenotype studies and for the molecular breeding of sheepgrass and other Poaceae species.

Sequencing breakthroughs for genomic ecology and evolutionary biology

Techniques involving whole-genome sequencing and whole-population sequencing (metagenomics) are beginning to revolutionize the study of ecology and evolution. This revolution is furthest advanced in the Bacteria and Archaea, and more sequence data are required for genomic ecology to be fully applied to the majority of eukaryotes. Recently developed next-generation sequencing technologies provide practical, massively parallel sequencing at lower cost and without the requirement for large, automated facilities, making genome and transcriptome sequencing and resequencing possible for more projects and more species. These sequencing methods include the 454 implementation of pyrosequencing, Solexa/Illumina reversible terminator technologies, polony sequencing and AB SOLiD. All of these methods use nanotechnology to generate hundreds of thousands of small sequence reads at one time. These technologies have the potential to bring the genomics revolution to whole populations, and to organisms such as endangered species or species of ecological and evolutionary interest. A future is now foreseeable where ecologists may resequence entire genomes from wild populations and perform population genetic studies at a genome, rather than gene, level. The new technologies for high throughput sequencing, their limitations and their applicability to evolutionary and environmental studies, are discussed in this review.

De novo sequence assembly and characterization of Lycoris aurea transcriptome using GS FLX titanium platform of 454 pyrosequencing

Background

Lycoris aurea, also called Golden Magic Lily, is an ornamentally and medicinally important species of the Amaryllidaceae family. To date, the sequencing of its whole genome is unavailable as a non-model organism. Transcriptomic information is also scarce for this species. In this study, we performed de novo transcriptome sequencing to produce the first comprehensive expressed sequence tag (EST) dataset for L. aurea using high-throughput sequencing technology.

Methodology and Principal Findings

Total RNA was isolated from leaves with sodium nitroprusside (SNP), salicylic acid (SA), or methyl jasmonate (MeJA) treatment, stems, and flowers at the bud, blooming, and wilting stages. Equal quantities of RNA from each tissue and stage were pooled to construct a cDNA library. Using 454 pyrosequencing technology, a total of 937,990 high quality reads (308.63 Mb) with an average read length of 329 bp were generated. Clustering and assembly of these reads produced a non-redundant set of 141,111 unique sequences, comprising 24,604 contigs and 116,507 singletons. All of the unique sequences were involved in the biological process, cellular component and molecular function categories by GO analysis. Potential genes and their functions were predicted by KEGG pathway mapping and COG analysis. Based on our sequence analysis and published literatures, many putative genes involved in Amaryllidaceae alkaloids synthesis, including PAL, TYDC OMT, NMT, P450, and other potentially important candidate genes, were identified for the first time in this Lycoris. Furthermore, 6,386 SSRs and 18,107 high-confidence SNPs were identified in this EST dataset.

Conclusions

The transcriptome provides an invaluable new data for a functional genomics resource and future biological research in L. aurea. The molecular markers identified in this study will provide a material basis for future genetic linkage and quantitative trait loci analyses, and will provide useful information for functional genomic research in future.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa

BACKGROUND

Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology.

RESULTS

Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology.

CONCLUSIONS

Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC

BACKGROUND

The MHC, which is regarded as the most polymorphic region in the genomes of jawed vertebrates, plays a central role in the immune system by encoding various proteins involved in the immune response. The chicken MHC-B genomic region has a highly streamlined gene content compared to mammalian MHCs. Its core region includes genes encoding Class I and Class IIB molecules but is only ~92Kb in length. Sequences of other galliform MHCs show varying degrees of similarity as that of chicken. The black grouse (Tetrao tetrix) is a wild galliform bird species which is an important model in conservation genetics and ecology. We sequenced the black grouse core MHC-B region and combined this with available data from related species (chicken, turkey, gold pheasant and quail) to perform a comparative genomics study of the galliform MHC. This kind of analysis has previously been severely hampered by the lack of genomic information on avian MHC regions, and the galliformes is still the only bird lineage where such a comparison is possible.

RESULTS

In this study, we present the complete genomic sequence of the MHC-B locus of black grouse, which is 88,390 bp long and contains 19 genes. It shows the same simplicity as, and almost perfect synteny with, the corresponding genomic region of chicken. We also use 454-transcriptome sequencing to verify expression in 17 of the black grouse MHC-B genes. Multiple sequence inversions of the TAPBP gene and TAP1-TAP2 gene block identify the recombination breakpoints near the BF and BLB genes. Some of the genes in the galliform MHC-B region also seem to have been affected by selective forces, as inferred from deviating phylogenetic signals and elevated rates of non-synonymous nucleotide substitutions.

CONCLUSIONS

We conclude that there is large synteny between the MHC-B region of the black grouse and that of other galliform birds, but that some duplications and rearrangements have occurred within this lineage. The MHC-B sequence reported here will provide a valuable resource for future studies on the evolution of the avian MHC genes and on links between immunogenetics and ecology of black grouse.

De novo sequencing and comparative analysis of the blueberry transcriptome to discover putative genes related to antioxidants

Blueberry (Vaccinium spp.) is an important small fruit crop rich in antioxidants. However, tissue-specific transcriptome and genomic data in public databases are not sufficient for an understanding of the molecular mechanisms associated with antioxidants, especially the biosynthesis of anthocyanins. Here, we obtained more than 64 million sequencing reads from blueberry skin and pulp using Illumina sequencing technology. De novo assemblies yielded 34,464 unigenes, among them 1236 transcripts and 862 putative transcription factors involved in the main antioxidant biosynthesis pathway were identified. Comparative transcript profiling allowed the identification of 92 differentially expressed genes with potential relevance in regulating the fruit metabolism and anthocyanin content during ripening. A series of qRT-PCR confirmed the high expression level of the anthocyanin pathway genes in the skin of the blue fruit from the in silico study. This sequence collection provides a significant resource for the blueberry research and breeding work.

Yellow lupin (Lupinus luteus L.) transcriptome sequencing: molecular marker development and comparative studies

Background

Yellow lupin (Lupinus luteus L.) is a minor legume crop characterized by its high seed protein content. Although grown in several temperate countries, its orphan condition has limited the generation of genomic tools to aid breeding efforts to improve yield and nutritional quality. In this study, we report the construction of 454-expresed sequence tag (EST) libraries, carried out comparative studies between L. luteus and model legume species, developed a comprehensive set of EST-simple sequence repeat (SSR) markers, and validated their utility on diversity studies and transferability to related species.

Results

Two runs of 454 pyrosequencing yielded 205 Mb and 530 Mb of sequence data for L1 (young leaves, buds and flowers) and L2 (immature seeds) EST- libraries. A combined assembly (L1L2) yielded 71,655 contigs with an average contig length of 632 nucleotides. L1L2 contigs were clustered into 55,309 isotigs. 38,200 isotigs translated into proteins and 8,741 of them were full length. Around 57% of L. luteus sequences had significant similarity with at least one sequence of Medicago, Lotus, Arabidopsis, or Glycine, and 40.17% showed positive matches with all of these species. L. luteus isotigs were also screened for the presence of SSR sequences. A total of 2,572 isotigs contained at least one EST-SSR, with a frequency of one SSR per 17.75 kbp. Empirical evaluation of the EST-SSR candidate markers resulted in 222 polymorphic EST-SSRs. Two hundred and fifty four (65.7%) and 113 (30%) SSR primer pairs were able to amplify fragments from L. hispanicus and L. mutabilis DNA, respectively. Fifty polymorphic EST-SSRs were used to genotype a sample of 64 L. luteus accessions. Neighbor-joining distance analysis detected the existence of several clusters among L. luteus accessions, strongly suggesting the existence of population subdivisions. However, no clear clustering patterns followed the accession’s origin.

Conclusion

L. luteus deep transcriptome sequencing will facilitate the further development of genomic tools and lupin germplasm. Massive sequencing of cDNA libraries will continue to produce raw materials for gene discovery, identification of polymorphisms (SNPs, EST-SSRs, INDELs, etc.) for marker development, anchoring sequences for genome comparisons and putative gene candidates for QTL detection.

Quantitative RNA-Seq analysis in non-model species: assessing transcriptome assemblies as a scaffold and the utility of evolutionary divergent genomic reference species

Background

How well does RNA-Seq data perform for quantitative whole gene expression analysis in the absence of a genome? This is one unanswered question facing the rapidly growing number of researchers studying non-model species. Using Homo sapiens data and resources, we compared the direct mapping of sequencing reads to predicted genes from the genome with mapping to de novo transcriptomes assembled from RNA-Seq data. Gene coverage and expression analysis was further investigated in the non-model context by using increasingly divergent genomic reference species to group assembled contigs by unique genes.

Results

Eight transcriptome sets, composed of varying amounts of Illumina and 454 data, were assembled and assessed. Hybrid 454/Illumina assemblies had the highest transcriptome and individual gene coverage. Quantitative whole gene expression levels were highly similar between using a de novo hybrid assembly and the predicted genes as a scaffold, although mapping to the de novo transcriptome assembly provided data on fewer genes. Using non-target species as reference scaffolds does result in some loss of sequence and expression data, and bias and error increase with evolutionary distance. However, within a 100 million year window these effect sizes are relatively small.

Conclusions

Predicted gene sets from sequenced genomes of related species can provide a powerful method for grouping RNA-Seq reads and annotating contigs. Gene expression results can be produced that are similar to results obtained using gene models derived from a high quality genome, though biased towards conserved genes. Our results demonstrate the power and limitations of conducting RNA-Seq in non-model species.

Genomic comparison between a virulent type A1 strain of Francisella tularensis and its attenuated O-antigen mutant

We report the complete genome sequences of TI0902, a highly virulent type A1 strain, and TIGB03, a related, attenuated chemical mutant strain. Compared to the wild type, the mutant strain had 45 point mutations and a 75.9-kb duplicated region that had not been previously observed in Francisella species.

Molecular characterization, immune responses and DNA protection activity of rock bream (Oplegnathus fasciatus), peroxiredoxin 6 (Prx6)

In this study, we describe the molecular characterization, immune responses of rock bream, Oplegnathus fasciatus peroxiredoxin 6 cDNA (RbPrx6) and DNA protection activity of its recombinant protein. The full-length cDNA sequence of RbPrx6 was identified after pyrosequencing of rock bream cDNA library. RbPrx6 consists of 663 bp open reading frame (ORF) that codes for a putative protein of 221 amino acids with predicted molecular mass of 27 kDa. It showed characteristic peroxiredoxin super-family domain similar to vertebrate Prx counterparts. In the pair-wise comparison, RbPrx6 showed the highest amino acid identity (92.8%) to Scophthalmus maximus Prx6. Real-time RT-PCR analysis revealed that constitutive expression of RbPrx6 transcripts in eleven tissues selected from un-challenged fish showing the highest level in liver. Synthetic polyinosinic:polycytidylic acid (poly I:C) and iridovirus containing supernatant, up-regulated the RbPrx6 mRNA in liver. Purified recombinant RbPrx6 protein was able to protect supercoiled plasmid DNA from damages that is induced by metal-catalyzed generation of reactive oxygen species. Our results suggest that RbPrx6 may play an important role in regulating oxidative stress by scavenging of ROS, involving immune reactions and minimizing the DNA damage in rock bream.

The Physalis peruviana leaf transcriptome: assembly, annotation and gene model prediction

Background

Physalis peruviana commonly known as Cape gooseberry is a member of the Solanaceae family that has an increasing popularity due to its nutritional and medicinal values. A broad range of genomic tools is available for other Solanaceae, including tomato and potato. However, limited genomic resources are currently available for Cape gooseberry.

Results

We report the generation of a total of 652,614 P. peruviana Expressed Sequence Tags (ESTs), using 454 GS FLX Titanium technology. ESTs, with an average length of 371 bp, were obtained from a normalized leaf cDNA library prepared using a Colombian commercial variety. De novo assembling was performed to generate a collection of 24,014 isotigs and 110,921 singletons, with an average length of 1,638 bp and 354 bp, respectively. Functional annotation was performed using NCBI’s BLAST tools and Blast2GO, which identified putative functions for 21,191 assembled sequences, including gene families involved in all the major biological processes and molecular functions as well as defense response and amino acid metabolism pathways. Gene model predictions in P. peruviana were obtained by using the genomes of Solanum lycopersicum (tomato) and Solanum tuberosum (potato). We predict 9,436 P. peruviana sequences with multiple-exon models and conserved intron positions with respect to the potato and tomato genomes. Additionally, to study species diversity we developed 5,971 SSR markers from assembled ESTs.

Conclusions

We present the first comprehensive analysis of the Physalis peruviana leaf transcriptome, which will provide valuable resources for development of genetic tools in the species. Assembled transcripts with gene models could serve as potential candidates for marker discovery with a variety of applications including: functional diversity, conservation and improvement to increase productivity and fruit quality. P. peruviana was estimated to be phylogenetically branched out before the divergence of five other Solanaceae family members, S. lycopersicum, S. tuberosum, Capsicum spp, S. melongena and Petunia spp.