A cDNA microarray analysis of the response to heat stress in hepatopancreas tissue of the porcelain crab Petrolisthes cinctipes

Intertidal zone organisms experience thermal stress during periods of low tide, and much work has shown that induction of heat shock proteins and ubiquitination occurs in response to this stress. However, less is known of other cellular pathways that are regulated following thermal stress in these organisms. Here, we used a functional genomics approach to identify genes that were up- and downregulated following heat stress in the intertidal porcelain crab, Petrolisthes cinctipes using custom cDNA microarrays made from 13,824 cloned P. cinctipes ESTs representing 6717 unique consensus sequences. Statistically significant differences in gene expression between heat stressed and control groups were determined with R/maanova. Genes upregulated following heat stress were involved with protein folding, protein degradation, protein synthesis and gluconeogenesis, suggesting that heat stress accelerated protein turnover. Genes downregulated following heat stress were involved with detoxification, oxygen transport, oxidative phosphorylation, and lipid metabolism, suggesting that the animals were avoiding the generation of reactive oxygen species. ESTs matching hypothetical proteins and ESTs that had no GenBank match were also found to have been both upregulated and downregulated following heat stress, suggesting that novel genes may be involved in the heat stress response.     

Comparative analysis of expressed sequence tags (ESTs) from Triticum monococcum shoot apical meristem at vegetative and reproductive stages

Triticum monococcum has recently drawn the attention of biologists to discover and utilize novel genes and alleles. To explore the molecular features of the genetic network governing floral transition in shoot apical meristem (SAM) of spring growth habit T. monococcum, two expressed sequence tag (EST) libraries containing 3,031 ESTs from vegetative SAM (VS) and 2,647 ESTs from early reproductive SAM (RS) were analyzed. Assembly of ESTs resulted in 2,303 unigenes for VS library (368 contigs and 1,935 singletons) and 1,890 unigenes (337 contigs and 1,553 singletons) for RS library. The 67.05 % of VS unigenes and 66.30 % of RS unigenes showed significant similarity with genes of known, putative and or unknown function, whereas the remaining 32.95 % of the VS unigenes and 33.7 % of RS unigenes displayed no significant match with the public protein database. The 1,064 and 866 unigenes of VS and RS libraries were assigned to functional categories using Pageman ontology tool. Further analysis revealed that the switch from VS to RS caused significant changes in the abundance of unigenes assigned to some functional categories. A total of 37 genes were identified which were significantly differentially expressed between vegetative and reproductive stages of T. monococcum SAM. Investigation of the differentially expressed genes revealed the importance of the genes involved in energy metabolism, ubiquitin/26S proteasome system, polyamines biosynthesis and signaling of reactive oxygen species in SAM differentiation towards floral transition in T. monococcum.     

EST-based in silico identification and in vitro test of antimicrobial peptides in Brassica napus

Background

Brassica napus is the third leading source of vegetable oil in the world after soybean and oil palm. The accumulation of gene sequences, especially expressed sequence tags (ESTs) from plant cDNA libraries, has provided a rich resource for genes discovery including potential antimicrobial peptides (AMPs). In this study, we used ESTs including those generated from B. napus cDNA libraries of seeds, pathogen-challenged leaves and deposited in the public databases, as a model, to perform in silico identification and consequently in vitro confirmation of putative AMP activities through a highly efficient system of recombinant AMP prokaryotic expression.

Results

In total, 35,788 were generated from cDNA libraries of pathogen-challenged leaves and 187,272 ESTs from seeds of B. napus, and the 644,998 ESTs of B. napus were downloaded from the EST database of PlantGDB. They formed 201,200 unigenes. First, all the known AMPs from the AMP databank (APD2 database) were individually queried against all the unigenes using the BLASTX program. A total of 972 unigenes that matched the 27 known AMP sequences in APD2 database were extracted and annotated using Blast2GO program. Among these unigenes, 237 unigenes from B. napus pathogen-challenged leaves had the highest ratio (1.15 %) in this unigene dataset, which is 13 times that of the unigene datasets of B. napus seeds (0.09 %) and 2.3 times that of the public EST dataset. About 87 % of each EST library was lipid-transfer protein (LTP) (32 % of total unigenes), defensin, histone, endochitinase, and gibberellin-regulated proteins. The most abundant unigenes in the leaf library were endochitinase and defensin, and LTP and histone in the pub EST library. After masking of the repeat sequence, 606 peptides that were orthologous matched to different AMP families were found. The phylogeny and conserved structural motifs of seven AMPs families were also analysed. To investigate the antimicrobial activities of the predicted peptides, 31 potential AMP genes belonging to different AMP families were selected to test their antimicrobial activities after bioinformatics identification. The AMP genes were all optimized according to Escherichia coli codon usage and synthetized through one-step polymerase chain reaction method. The results showed that 28 recombinant AMPs displayed expected antimicrobial activities against E. coli and Micrococcus luteus and Sclerotinia sclerotiorum strains.

Conclusion

The study not only significantly expanded the number of known/predicted peptides, but also contributed to long-term plant genetic improvement for increased resistance to diverse pathogens of B.napus. These results proved that the high-throughput method developed that combined an in silico procedure with a recombinant AMP prokaryotic expression system is considerably efficient for identification of new AMPs from genome or EST sequence databases.

Electronic supplementary material

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

Use of a large-scale Triticeae expressed sequence tag resource to reveal gene expression profiles in hexaploid wheat (Triticum aestivum L.).

The US Wheat Genome Project, funded by the National Science Foundation, developed the first large public Triticeae expressed sequence tag (EST) resource. Altogether, 116,272 ESTs were produced, comprising 100,674 5' ESTs and 15 598 3' ESTs. These ESTs were derived from 42 cDNA libraries, which were created from hexaploid bread wheat (Triticum aestivum L.) and its close relatives, including diploid wheat (T. monococcum L. and Aegilops speltoides L.), tetraploid wheat (T. turgidum L.), and rye (Secale cereale L.), using tissues collected from various stages of plant growth and development and under diverse regimes of abiotic and biotic stress treatments. ESTs were assembled into 18,876 contigs and 23,034 singletons, or 41,910 wheat unigenes. Over 90% of the contigs contained fewer than 10 EST members, implying that the ESTs represented a diverse selection of genes and that genes expressed at low and moderate to high levels were well sampled. Statistical methods were used to study the correlation of gene expression patterns, based on the ESTs clustered in the 1536 contigs that contained at least 10 5' EST members and thus representing the most abundant genes expressed in wheat. Analysis further identified genes in wheat that were significantly upregulated (p < 0.05) in tissues under various abiotic stresses when compared with control tissues. Though the function annotation cannot be assigned for many of these genes, it is likely that they play a role associated with the stress response. This study predicted the possible functionality for 4% of total wheat unigenes, which leaves the remaining 96% with their functional roles and expression patterns largely unknown. Nonetheless, the EST data generated in this project provide a diverse and rich source for gene discovery in wheat.     

Development and annotation of perennial Triticeae ESTs and SSR markers

Triticeae contains hundreds of species of both annual and perennial types. Although substantial genomic tools are available for annual Triticeae cereals such as wheat and barley, the perennial Triticeae lack sufficient genomic resources for genetic mapping or diversity research. To increase the amount of sequence information available in the perennial Triticeae, three expressed sequence tag (EST) libraries were developed and annotated for Pseudoroegneria spicata, a mixture of both Elymus wawawaiensis and E. lanceolatus, and a Leymus cinereus x L. triticoides interspecific hybrid. The ESTs were combined into unigene sets of 8 780 unigenes for P. spicata, 11 281 unigenes for Leymus, and 7 212 unigenes for Elymus. Unigenes were annotated based on putative orthology to genes from rice, wheat, barley, other Poaceae, Arabidopsis, and the non-redundant database of the NCBI. Simple sequence repeat (SSR) markers were developed, tested for amplification and polymorphism, and aligned to the rice genome. Leymus EST markers homologous to rice chromosome 2 genes were syntenous on Leymus homeologous groups 6a and 6b (previously 1b), demonstrating promise for in silico comparative mapping. All ESTs and SSR markers are available on an EST information management and annotation database (http://titan.biotec.uiuc.edu/triticeae/).     

Suppression Subtractive Hybridization Analysis of Genes Regulated by Application of Exogenous Abscisic Acid in Pepper Plant (Capsicum annuum L.) Leaves under Chilling Stress

Low temperature is one of the major factors limiting pepper (Capsicum annuum L.) production during winter and early spring in non-tropical regions. Application of exogenous abscisic acid (ABA) effectively alleviates the symptoms of chilling injury, such as wilting and formation of necrotic lesions on pepper leaves; however, the underlying molecular mechanism is not understood. The aim of this study was to identify genes that are differentially up- or downregulated in ABA-pretreated hot pepper seedlings incubated at 6°C for 48 h, using a suppression subtractive hybridization (SSH) method. A total of 235 high-quality ESTs were isolated, clustered and assembled into a collection of 73 unigenes including 18 contigs and 55 singletons. A total of 37 unigenes (50.68%) showed similarities to genes with known functions in the non-redundant database; the other 36 unigenes (49.32%) showed low similarities or unknown functions. Gene ontology analysis revealed that the 37 unigenes could be classified into nine functional categories. The expression profiles of 18 selected genes were analyzed using quantitative RT-PCR; the expression levels of 10 of these genes were at least two-fold higher in the ABA-pretreated seedlings under chilling stress than water-pretreated (control) plants under chilling stress. In contrast, the other eight genes were downregulated in ABA-pretreated seedlings under chilling stress, with expression levels that were one-third or less of the levels observed in control seedlings under chilling stress. These results suggest that ABA can positively and negatively regulate genes in pepper plants under chilling stress.     

Identification and analysis of expressed sequence tags present in xylem tissues of kelampayan (Neolamarckia cadamba (Roxb.) Bosser)

The large-scale genomic resource for kelampayan was generated from a developing xylem cDNA library. A total of 6,622 high quality expressed sequence tags (ESTs) were generated through high-throughput 5’ EST sequencing of cDNA clones. The ESTs were analyzed and assembled to generate 4,728 xylogenesis unigenes distributed in 2,100 contigs and 2,628 singletons. About 59.3 % of the ESTs were assigned with putative identifications whereas 40.7 % of the sequences showed no significant similarity to any sequences in GenBank. Interestingly, most genes involved in lignin biosynthesis and several other cell wall biosynthesis genes were identified in the kelampayan EST database. The identified genes in this study will be candidates for functional genomics and association genetic studies in kelampayan aiming at the production of high value forests.     

Transcriptome generation and analysis from spleen of Indian catfish, Clarias batrachus (Linnaeus, 1758) through normalized cDNA library

Catfishes are commercially important fish for both the fisheries and aquaculture industry. Clarias batrachus, an Indian catfish species is economically important owing to its high demand. A normalized cDNA library was constructed from spleen of the Indian catfish to identify genes associated with immune function. One thousand nine hundred thirty seven ESTs were submitted to the GenBank with an average read length of approximately 700 bp. Clustering analysis of ESTs yielded 1,698 unique sequences, including 184 contigs and 1,514 singletons. Significant homology to known genes was found by homology searches against data in GenBank in 576 (34 %) ESTs, including similarity to functionally annotated unigenes for 158 ESTs. Additionally, 433 ESTs revealed similarity to unigenes and ESTs in the dbEST but the remaining 658 EST sequences (39 %) did not match any sequence in GenBank. Of a total of 1,698 ESTs generated, 65 ESTs were found to be associated with immune functions. Gene Ontology and KEGG pathway analyses of C. batrachus ESTs collectively revealed a preponderance of immune relevant pathways apart from the presence of pathways involved in protein processing, localization, folding and protein degradation. This study constitutes first EST analysis of lymphoid organ in aquaculturally important Indian catfish species and could pave the way for further research of immune-related genes and functional genomics in this catfish.     

Expressed sequence tags from phytophthora sojae reveal genes specific to development and infection

Six unique expressed sequence tag (EST) libraries were generated from four developmental stages of Phytophthora sojae P6497. RNA was extracted from mycelia, swimming zoospores, germinating cysts, and soybean (Glycine max (L.) Merr.) cv. Harosoy tissues heavily infected with P. sojae. Three libraries were created from mycelia growing on defined medium, complex medium, and nutrient-limited medium. The 26,943 high-quality sequences obtained clustered into 7,863 unigenes composed of 2,845 contigs and 5,018 singletons. The total number of P. sojae unigenes matching sequences in the genome assembly was 7,412 (94%). Of these unigenes, 7,088 (90%) matched gene models predicted from the P. sojae sequence assembly, but only 2,047 (26%) matched P. ramorum gene models. Analysis of EST frequency from different growth conditions and morphological stages revealed genes that were specific to or highly represented in particular growth conditions and life stages. Additionally, our results indicate that, during infection, the pathogen derives most of its carbon and energy via glycolysis of sugars in the plant. Sequences identified with putative roles in pathogenesis included avirulence homologs possessing the RxLR motif, elicitins, and hydrolytic enzymes. This large collection of P. sojae ESTs will serve as a valuable public genomic resource.