Comparative analysis of sequences expressed during the liquid-cultured mycelia and fruit body stages of Pleurotus ostreatus

To characterize genes involved in fruit body development, two complementary DNA (cDNA) libraries were constructed from RNA isolated from liquid-cultured mycelia and fruit bodies of Pleurotus ostreatus. Using single-pass sequencing of cDNA clones, 952 and 1069 expressed sequence tags (ESTs) were generated from liquid-cultured mycelia and fruit body cDNA library, respectively. A BLASTX search revealed that 390 of the liquid-cultured mycelia ESTs (41%) and 531 of the fruit body ESTs (50%) showed significant similarity to protein sequences described in the nonredudant database (E values < or =1 x 10(-5)). When liquid-cultured mycelia and fruit body ESTs were compared by the SeqMan II program, among the total of 2021 ESTs, 1256 ESTs were unigenes, and 66 unigenes (5.3%) were commonly expressed during both stages. The functional catalogs of the ESTs were made by comparison with functionally identified Saccharomyces cerevisiae genes. Liquid-cultured mycelium ESTs were compared with fruit body ESTs and changes of the expressed genes during fruit body development were analyzed.     

Expressed sequence tag analysis in tef (Eragrostis tef (Zucc) Trotter).

Tef (Eragrostis tef (Zucc.) Trotter) is the most important cereal crop in Ethiopia; however, there is very little DNA sequence information available for this species. Expressed sequence tags (ESTs) were generated from 4 cDNA libraries: seedling leaf, seedling root, and inflorescence of E. tef and seedling leaf of Eragrostis pilosa, a wild relative of E. tef. Clustering of 3603 sequences produced 530 clusters and 1890 singletons, resulting in 2420 tef unigenes. Approximately 3/4 of tef unigenes matched protein or nucleotide sequences in public databases. Annotation of unigenes associated 68% of the putative tef genes with gene ontology categories. Identification of the translated unigenes for conserved protein domains revealed 389 protein family domains (Pfam), the most frequent of which was protein kinase. A total of 170 ESTs containing simple sequence repeats (EST-SSRs) were identified and 80 EST-SSR markers were developed. In addition, 19 single-nucleotide polymorphism (SNP) and (or) insertion-deletion (indel) and 34 intron fragment length polymorphism (IFLP) markers were developed. The EST database and molecular markers generated in this study will be valuable resources for further tef genetic research.     

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.     

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.     

An ordered EST catalogue and gene expression profiles of cassava (<Emphasis Type="Italic">Manihot esculenta</Emphasis>) at key growth stages

A cDNA library was constructed from the root tissues of cassava variety Huanan 124 at the root bulking stage. A total of 9,600 cDNA clones from the library were sequenced with single-pass from the 5′-terminus to establish a catalogue of expressed sequence tags (ESTs). Assembly of the resulting EST sequences resulted in 2,878 putative unigenes. Blastn analysis showed that 62.6% of the unigenes matched with known cassava ESTs and the rest had no ‘hits’ against the cassava database in the integrative PlantGDB database. Blastx analysis showed that 1,715 (59.59%) of the unigenes matched with one or more GenBank protein entries and 1,163 (40.41%) had no ‘hits’. A cDNA microarray with 2,878 unigenes was developed and used to analyze gene expression profiling of Huanan 124 at key growth stages including seedling, formation of root system, root bulking, and starch maturity. Array data analysis revealed that (1) the higher ratio of up-regulated ribosome-related genes was accompanied by a high ratio of up-regulated ubiquitin, proteasome-related and protease genes in cassava roots; (2) starch formation and degradation simultaneously occur at the early stages of root development but starch degradation is declined partially due to decrease in UDP-glucose dehydrogenase activity with root maturity; (3) starch may also be synthesized in situ in roots; (4) starch synthesis, translocation, and accumulation are also associated probably with signaling pathways that parallel Wnt, LAM, TCS and ErbB signaling pathways in animals; (5) constitutive expression of stress-responsive genes may be due to the adaptation of cassava to harsh environments during long-term evolution.     

Analysis of expressed sequence tags from Gibberella zeae (anamorph Fusarium graminearum)

Gibberella zeae is a broad host range pathogen that infects many crop plants, including wheat and barley, and causes head blight and rot diseases throughout the world. To better understand fungal development and pathogenicity, we have generated 7996 ESTs from three cDNA libraries. Two libraries were generated from carbon-(C-) and nitrogen- (N-) starved mycelia and one library was generated from cultures of maturing perithecia (P). In other fungal pathogens, starvation conditions have been shown to act as cues to induce infection-related gene expression. To assign putative function to cDNAs, sequences were initially assembled using StackPack. The estimated total number of genes identified from the three EST databases was 2110: 1088 contigs and 1022 singleton sequences. These 2110 sequences were compared to a yeast protein sequence reference set and to the GenBank nonredundant database using BLASTX. Based on presumptive gene function identified by this process, we found that the two starved cultures had similar, but not identical, patterns of gene expression, whereas the developmental cultures were distinct in their pattern of expression. Of the three libraries, the perithecium library had the greatest percentage (46%) of ESTS falling into the "unclassified" category. Homologues of some known fungal virulence or pathogenicity factors were found primarily in the N- and C-libraries. Comparisons also were made with ESTs from the related fungi, Neurospora crassa and Magnaporthe grisea and the genomic sequence of N. crassa.     

Analysis of expressed sequence tags from grapevine flower and fruit and development of simple sequence repeat markers

A total of 6,230 EST sequences were produced from 7,561 clones in a cDNA library generated from grapevine (Vitis vinifera cv. ‘Summer Black’) flower and fruit tissues in this study. After cluster and assembly analysis of the datasets, 3,582 unigenes (GenBank accession numbers GW836604–GW840185) were established, among which 381 were new grapevine EST sequences. Out of the 381 new ESTs, 289 could be mapped on the 19 grapevine chromosomes. 913 unique ESTs with known or putative functions were assigned to 11 putative cellular roles. 540 potentially workable grapevine EST-SSRs were developed from 3,582 unigenes and about 42.6% of these unigenes were identified as true-to-type SSR loci and could amplify polymorphic bands from 22 individual plants of V. vinifera L, indicating that grapevine EST datasets are a valuable source for the development of functional simple sequence repeat (SSR) markers.     

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.