Tissue expression map of a large number of expressed sequence tags and its application to in silico screening of stress response genes in common wheat

In order to assess global changes in gene expression patterns in stress-induced tissues, we conducted large-scale analysis of expressed sequence tags (ESTs) in common wheat. Twenty-one cDNA libraries derived from stress-induced tissues, such as callus, as well as liquid cultures and abiotic stress conditions (temperature treatment, desiccation, photoperiod, moisture and ABA) were constructed. Several thousand colonies were randomly selected from each of these 21 cDNA libraries and sequenced from both the 5′ and 3′ ends. By computing abundantly expressed ESTs, correlated expression patterns of genes across the tissues were monitored. Furthermore, the relationships between gene expression profiles among the stress-induced tissues were inferred from the gene expression patterns. Multi-dimensional analysis of EST data is analogous to microarray experiments. As an example, genes specifically induced and/or suppressed by cold acclimation and heat-shock treatments were selected in silico. Four hundred and ninety genes showing fivefold induction or 218 genes for suppression in comparison to the control expression level were selected. These selected genes were annotated with the BLAST search. Furthermore, gene ontology was conducted for these genes with the InterPro search. Because genes regulated in response to temperature treatment were successfully selected, this method can be applied to other stress-treated tissues. Then, the method was applied to screen genes in response to abiotic stresses such as drought and ABA treatments. In silico selection of screened genes from virtual display should provide a powerful tool for functional plant genomics.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

EST sequencing and time course microarray hybridizations identify more than 700 <Emphasis Type="Italic">Medicago truncatula</Emphasis> genes with developmental expression regulation in flowers and pods

To evaluate the molecular mechanisms during pod and seed formation in legumes, starting with the development of reproductive organs, we constructed two cDNA libraries from developing flowers (MtFLOW) and pods including seeds (MtPOSE) of the model plant Medicago truncatula Gaertner. A total of 2,516 expressed sequence tags (ESTs) clustered into 1,776 nonredundant sequences (2k-set), which were annotated and assigned to functional classes. While about 30% of the ESTs encoded proteins of yet unknown function, typical annotations pointed to seed storage proteins, LTPs and lipoxygenases. The 2k-set was used to upgrade Mt6k-RIT microarrays (Küster et al. in J Biotechnol 108: 95, 2004) to Mt8k versions representing approximately 6,300 nonredundant M. truncatula genes. These were used to perform time course expression profiling studies based on hybridizations of samples that covered eight different developmental stages from flower buds to almost mature pods versus leaves as a common reference. About 180 up- and 70 downregulated genes were typically found for each stage and in total, 782 genes were either twofold up- or downregulated in at least one of the eight stages investigated. Based on this set, a combination of self-organizing map and hierarchical clustering revealed genes displaying expression regulation during characteristic stages of M. truncatula flower and pod development. Amongst those, several genes encoded proteins related to seed metabolism and development including novel regulators and proteins involved in signaling.Electronic supplementary material Supplementary material is available for this article at     

Large-scale identification of expressed sequence tags involved in rice and rice blast fungus interaction

To better understand the molecular basis of the defense response against the rice blast fungus (Magnaporthe grisea), a large-scale expressed sequence tag (EST) sequencing approach was used to identify genes involved in the early infection stages in rice (Oryza sativa). Six cDNA libraries were constructed using infected leaf tissues harvested from 6 conditions: resistant, partially resistant, and susceptible reactions at both 6 and 24 h after inoculation. Two additional libraries were constructed using uninoculated leaves and leaves from the lesion mimic mutant spl11. A total of 68,920 ESTs were generated from 8 libraries. Clustering and assembly analyses resulted in 13,570 unique sequences from 10,934 contigs and 2,636 singletons. Gene function classification showed that 42% of the ESTs were predicted to have putative gene function. Comparison of the pathogen-challenged libraries with the uninoculated control library revealed an increase in the percentage of genes in the functional categories of defense and signal transduction mechanisms and cell cycle control, cell division, and chromosome partitioning. In addition, hierarchical clustering analysis grouped the eight libraries based on their disease reactions. A total of 7,748 new and unique ESTs were identified from our collection compared with the KOME full-length cDNA collection. Interestingly, we found that rice ESTs are more closely related to sorghum (Sorghum bicolor) ESTs than to barley (Hordeum vulgare), wheat (Triticum aestivum), and maize (Zea mays) ESTs. The large cataloged collection of rice ESTs in this study provides a solid foundation for further characterization of the rice defense response and is a useful public genomic resource for rice functional genomics studies.     

Gene expression profiling of the plant pathogenic basidiomycetous fungus Rhizoctonia solani AG 4 reveals putative virulence factors

Rhizoctonia solani is a ubiquitous basidiomycetous soilborne fungal pathogen causing damping-off of seedlings, aerial blights and postharvest diseases. To gain insight into the molecular mechanisms of pathogenesis a global approach based on analysis of expressed sequence tags (ESTs) was undertaken. To get broad gene-expression coverage, two normalized EST libraries were developed from mycelia grown under high nitrogen-induced virulent and low nitrogen/methylglucose-induced hypovirulent conditions. A pilot-scale assessment of gene diversity was made from the sequence analyses of the two libraries. A total of 2280 cDNA clones was sequenced that corresponded to 220 unique sequence sets or clusters (contigs) and 805 singlets, making up a total of 1025 unique genes identified from the two virulence-differentiated cDNA libraries. From the total sequences, 295 genes (38.7%) exhibited strong similarities with genes in public databases and were categorized into 11 functional groups. Approximately 61.3% of the R. solani ESTs have no apparent homologs in publicly available fungal genome databases and are considered unique genes. We have identified several cDNAs with potential roles in fungal pathogenicity, virulence, signal transduction, vegetative incompatibility and mating, drug resistance, lignin degradation, bioremediation and morphological differentiation. A codon-usage table has been formulated based on 14694 R. solani EST codons. Further analysis of ESTs might provide insights into virulence mechanisms of R. solani AG 4 as well as roles of these genes in development, saprophytic colonization and ecological adaptation of this important fungal plant pathogen.     

CpG island libraries from human Chromosomes 18 and 22: landmarks for novel genes

CpG islands are found at the 5′ end of approximately 60% of human genes and so are important genomic landmarks. They are concentrated in early-replicating, highly acetylated gene-rich regions. With respect to CpG island content, human Chrs 18 and 22 are very different from each other: Chr 18 appears to be CpG island poor, whereas Chr 22 appears to be CpG island rich. We have constructed and validated CpG island libraries from flow-sorted Chrs 18 and 22 and used these to estimate the difference in number of CpG islands found on these two chromosomes. These libraries contain normalized collections of sequences from the 5′ end of genes. Clones from the libraries were sequenced and compared with the sequence databases; one third matched ESTs, thus anchoring these ESTs at the 5′ end of their gene. However, it was striking that many clones either had no match or matched only existing CpG island clones. This suggests that a significant proportion of 5′ gene sequences are absent from databases, presumably either because they are difficult to clone or the gene is poorly expressed and/or has a restricted expression pattern. This point should be taken into consideration if the currently available libraries are those used for the elucidation of complete, as opposed to partial, gene sequences. The Chr 18 and 22 CpG island libraries are a sequence resource for the isolation of such 5′ gene sequences from specific human chromosomes. Received: 15 November 1999 / Accepted: 31 January 2000     

Gene expression profiling in porcine mammary gland during lactation and identification of breed-and developmental-stage-specific genes

A total of 28941 ESTs were sequenced from five 5′-directed non-normalized cDNA libraries, which were assembled into 2212 contigs and 5642 singlets using CAP3. These sequences were annotated and clustered into 6857 unique genes, 2072 of which having no functional annotations were considered as novel genes. These genes were further classified into Gene Ontology categories. By comparing the expression profiles, we identified some breed-and developmental-stage-specific gene groups. These genes may be relative to reproductive performance or play important roles in milk synthesis, secretion and mammary involution. The unknown EST sequences and expression profiles at different developmental stages and breeds are very important resources for further research.     

Genome-wide analysis of defense-responsive genes in bacterial blight resistance of rice mediated by the recessive<Emphasis Type="Italic"> R</Emphasis> gene<Emphasis Type="Italic"> xa13</Emphasis>

Defense responses triggered by dominant and recessive disease resistance ( R) genes are presumed to be regulated by different molecular mechanisms. In order to characterize the genes activated in defense responses against bacterial blight mediated by the recessive R gene xa13, two pathogen-induced subtraction cDNA libraries were constructed using the resistant rice line IRBB13—which carries xa13 —and its susceptible, near-isogenic, parental line IR24. Clustering analysis of expressed sequence tags (ESTs) identified 702 unique expressed sequences as being involved in the defense responses triggered by xa13; 16% of these are new rice ESTs. These sequences define 702 genes, putatively encoding a wide range of products, including defense-responsive genes commonly involved in different host-pathogen interactions, genes that have not previously been reported to be associated with pathogen-induced defense responses, and genes (38%) with no homology to previously described functional genes. In addition, R -like genes putatively encoding nucleotide-binding site/leucine rich repeat (NBS-LRR) and LRR receptor kinase proteins were observed to be induced in the disease resistance activated by xa13. A total of 568 defense-responsive ESTs were mapped to 588 loci on the rice molecular linkage map through bioinformatic analysis. About 48% of the mapped ESTs co-localized with quantitative trait loci (QTLs) for resistance to various rice diseases, including bacterial blight, rice blast, sheath blight and yellow mottle virus. Furthermore, some defense-responsive sequences were conserved at similar locations on different chromosomes. These results reveal the complexity of xa13 -mediated resistance. The information obtained in this study provides a large source of candidate genes for understanding the molecular bases of defense responses activated by recessive R genes and of quantitative disease resistance.Electronic Supplementary Material Supplementary material is available in the online version of this article at The first two authors contributed equally to this workCommunicated by R. Hagemann     

Genes expressed in zoospores of<Emphasis Type="Italic"> Phytophthora nicotianae</Emphasis>

The genus Phytophthora includes many highly destructive plant pathogens. In many Phytophthora species, pathogen dispersal and initiation of plant infection are achieved by motile, biflagellate zoospores that are chemotactically attracted to suitable infection sites. In order to study gene expression in zoospores, we have constructed a cDNA library using mRNA from zoospores of Phytophthora nicotianae. The library was arrayed and screened using probes derived from mycelium or zoospore mRNA. More than 400 clones representing genes preferentially expressed in zoospores were identified and sequenced from the 5 end of the insert. The expressed sequence tags (ESTs) generated were found to represent 240 genes. The ESTs were compared to sequences in GenBank and in the Phytophthora Genome Consortium database, and classified according to putative function based on homology to known proteins. To further characterize the identified genes, a colony array was created on replicate nylon filters and screened with probes derived from four Phytophthora developmental stages including zoospores, germinating cysts, vegetative mycelium and sporulating hyphae, and from inoculated and uninoculated tobacco seedlings. Data from sequence analysis and colony array screening were compiled into a local database, and searched to identify genes that are preferentially expressed in zoospores for future functional analysis.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by C. A. M. J. J. van den Hondel     

Full-length enriched cDNA library construction from tissues related to energy metabolism in pigs

Genome sequencing of the pig is being accelerated because of its importance as an evolutionary and biomedical model animal as well as a major livestock animal. However, information on expressed porcine genes is insufficient to allow annotation and use of the genomic information. A series of expressed sequence tags of 5′ ends of five full-length enriched cDNA libraries (SUSFLECKs) were functionally characterized. SUSFLECKs were constructed from porcine abdominal fat, induced fat cells, loin muscle, liver, and pituitary gland, and were composed of non-normalized and normalized libraries. A total of 55,658 ESTs that were sequenced once from the 5′ ends of clones were produced and assembled into 17,684 unique sequences with 7,736 contigs and 9,948 singletons. In Gene Ontology analysis, two significant biological process leaf nodes were found: gluconeogenesis and translation elongation. In functional domain analysis based on the Pfam database, the beta transducin repeat domain of WD40 protein was the most frequently occurring domain. Twelve genes, including SLC25A6, EEF1G, EEF1A1, COX1, ACTA1, SLA, and ANXA2, were significantly more abundant in fat tissues than in loin muscle, liver, and pituitary gland in the SUSFLECKs. These characteristics of SUSFLECKs determined by EST analysis can provide important insight to discover the functional pathways in gene networks and to expand our understanding of energy metabolism in the pig.     

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.     

Effects of postharvest storage and dormancy status on ABA content, metabolism, and expression of genes involved in ABA biosynthesis and metabolism in potato tuber tissues

At harvest, and for an indeterminate period thereafter, potato tubers will not sprout and are physiologically dormant. Abscisic acid (ABA) has been shown to play a critical role in tuber dormancy control but the mechanisms controlling ABA content during dormancy as well as the sites of ABA synthesis and catabolism are unknown. As a first step in defining the sites of synthesis and cognate processes regulating ABA turnover during storage and dormancy progression, gene sequences encoding the ABA biosynthetic enzymes zeaxanthin epoxidase (ZEP) and 9-cis-epoxycarotenoid dioxygenase (NCED) and three catabolism-related genes were used to quantify changes in their relative mRNA abundances in three specific tuber tissues (meristems, their surrounding periderm and underlying cortex) by qRT-PCR. During storage, StZEP expression was relatively constant in meristems, exhibited a biphasic pattern in periderm with transient increases during early and mid-to-late-storage, and peaked during mid-storage in cortex. Expression of two members of the potato NCED gene family was found to correlate with changes in ABA content in meristems (StNCED2) and cortex (StNCED1). Conversely, expression patterns of three putative ABA-8′-hydroxylase (CYP707A) genes during storage varied in a tissue-specific manner with expression of two of these genes rising in meristems and periderm and declining in cortex during storage. These results suggest that ABA synthesis and metabolism occur in all tuber tissues examined and that tuber ABA content during dormancy is the result of a balance of synthesis and metabolism that increasingly favors catabolism as dormancy ends and may be controlled at the level of StNCED and StCYP707A gene activities Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Characterization of gene expression during potato tuber development in individuals and populations using the luciferase reporter system

Analysis of gene expression and enzyme activity in pooled tuber samples has previously indicated different developmental events occurring in a fixed sequential order during tuber development, starting with the up-regulation of starch synthesis then induction of protein storage followed by cell division and cell enlargement. In this report we analysed in vivo promoter activity of genes related to cell division and storage of reserves during tuber development in individual in vitro tubers, using the non invasive firefly luciferase reporter system. The average activity of the storage related promoters (AGPaseS and Pat21) was up-regulated prior to visible swelling, while the average activity of both cell cycle genes (cycB1;1 and CDC2a) showed an up-regulation after the onset of swelling. However, this novel system allowed expression analysis in individual tubers, which showed a variable up-regulation of both storage genes in relation to the moment of swelling, from 4 days before to 10 days after the onset of swelling. We conclude that during the first stages of tuber development, the moment of storage gene induction is independent from swelling. These results indicate that the developmental program of potato tubers does not consist of a fixed sequential order of events, but consists of independent developmental programs (storage and swelling), together resulting in the formation of a potato tuber. It is concluded that analysis of developmental programs by studying individuals may result in new insights, possibly obscured when using pooled samples.     

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.