Abiotic stress reduces plant growth and crop productivity. However, the mechanism underlying posttranscriptional regulations of stress response remains elusive. Herein, we report the posttranscriptional mechanism of nucleocytoplasmic RNA transport of stress-responsive transcripts mediated by EgRBP42, a heterogeneous nuclear ribonucleoprotein-like RNA-binding protein from oil palm, which could be necessary for rapid protein translation to confer abiotic stress tolerance in plants. Transgenic Arabidopsis overexpressing EgRBP42 showed early flowering through alteration of gene expression of flowering regulators and exhibited tolerance towards heat, cold, drought, flood, and salinity stresses with enhanced poststress recovery response by increasing the expression of its target stress-responsive genes. EgRBP42 harbours nucleocytoplasmic shuttling activity mediated by the nuclear localization signal and the M9-like domain of EgRBP42 and interacts directly with regulators in the nucleus, membrane, and the cytoplasm. EgRBP42 regulates the nucleocytoplasmic RNA transport of target stress-responsive transcripts through direct binding to their AG-rich motifs. Additionally, EgRBP42 transcript and protein induction by environmental stimuli are regulated at the transcriptional and posttranscriptional levels. Taken together, the posttranscriptional regulation of RNA transport mediated by EgRBP42 may change the stress-responsive protein profiles under abiotic stress conditions leading to a better adaptation of plants to environmental changes.
The mycoheterotrophic lifestyle has enabled some plant lineages to obtain carbon from their mycorrhizal symbionts. The mycoheterotrophic genus Epirixanthes (Polygalaceae) consists of six species from tropical Asia. Although it is probably closely related to the chlorophyllous genus Salomonia and linked to arbuscular mycorrhizal fungi, lack of DNA sequence data has thus far prevented these hypotheses from being tested. Therefore, the evolutionary history of Epirixanthes remains largely unknown.
Bioinformatics analyses of caleosin/peroxygenases (CLO/PXG) demonstrated that these genes are present in the vast majority of Viridiplantae taxa for which sequence data are available. Functionally active CLO/PXG proteins with roles in abiotic stress tolerance and lipid droplet storage are present in some Trebouxiophycean and Chlorophycean green algae but are absent from the small number of sequenced Prasinophyceaen genomes. CLO/PXG-like genes are expressed during dehydration stress in Charophyte algae, a sister clade of the land plants (Embryophyta). CLO/PXG-like sequences are also present in all of the >300 sequenced Embryophyte genomes, where some species contain as many as 10-12 genes that have arisen via selective gene duplication. Angiosperm genomes harbour at least one copy each of two distinct CLO/PX isoforms, termed H (high) and L (low), where H-forms contain an additional C-terminal motif of about 30-50 residues that is absent from L-forms. In contrast, species in other Viridiplantae taxa, including green algae, non-vascular plants, ferns and gymnosperms, contain only one (or occasionally both) of these isoforms per genome. Transcriptome and biochemical data show that CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression. CLO/PXG proteins can associate with cytosolic lipid droplets and/or bilayer membranes. Many of the analysed isoforms also have peroxygenase activity and are involved in oxylipin metabolism. The distribution of CLO/PXG-like genes is consistent with an origin >1 billion years ago in at least two of the earliest diverging groups of the Viridiplantae, namely the Chlorophyta and the Streptophyta, after the Viridiplantae had already diverged from other Archaeplastidal groups such as the Rhodophyta and Glaucophyta. While algal CLO/PXGs have roles in lipid packaging and stress responses, the Embryophyte proteins have a much wider spectrum of roles and may have been instrumental in the colonisation of terrestrial habitats and the subsequent diversification as the major land flora.
The possibility of controlling Fusarium wilt--caused by Fusarium oxysporum sp. cubensec (race 4)--was investigated by genetic engineering of banana plants for constitutive expression of rice thaumatin-like protein (tlp) gene. Transgene was introduced to cauliflower-like bodies' cluster, induced from meristemic parts of male inflorescences, using particle bombardment with plasmid carrying a rice tlp gene driving by the CaMV 35S promoter. Hygromycin B was used as the selection reagent. The presence and integration of rice tlp gene in genomic DNA confirmed by PCR and Southern blot analyses. RT-PCR revealed the expression of transgene in leaf and root tissues in transformants. Bioassay of transgenic banana plants challenged with Fusarium wilt pathogen showed that expression of TLP enhanced resistance to F. oxysporum sp. cubensec (race 4) compared to control plants.
SAG21/LEA5 is an unusual late embryogenesis abundant protein in Arabidopsis thaliana, that is primarily mitochondrially located and may be important in regulating translation in both chloroplasts and mitochondria. SAG21 expression is regulated by a plethora of abiotic and biotic stresses and plant growth regulators indicating a complex regulatory network. To identify key transcription factors regulating SAG21 expression, yeast-1-hybrid screens were used to identify transcription factors that bind the 1685 bp upstream of the SAG21 translational start site. Thirty-three transcription factors from nine different families bound to the SAG21 promoter, including members of the ERF, WRKY and NAC families. Key binding sites for both NAC and WRKY transcription factors were tested through site directed mutagenesis indicating the presence of cryptic binding sites for both these transcription factor families. Co-expression in protoplasts confirmed the activation of SAG21 by WRKY63/ABO3, and SAG21 upregulation elicited by oligogalacturonide elicitors was partially dependent on WRKY63, indicating its role in SAG21 pathogen responses. SAG21 upregulation by ethylene was abolished in the erf1 mutant, while wound-induced SAG21 expression was abolished in anac71 mutants, indicating SAG21 expression can be regulated by several distinct transcription factors depending on the stress condition.
Basal stem rot (BSR) is a major disease of oil palm caused by a pathogenic fungus, Ganoderma boninense. However, the interaction between the host plant and its pathogen is not well characterized. To better understand the response of oil palm to G. boninense, transcript profiles of eleven putative defence-related genes from oil palm were measured by quantitative reverse-transcription (qRT)-PCR in the roots of oil palms treated with G. boninense from 3 to 12 weeks post infection (wpi). These transcripts encode putative Bowman-Birk serine protease inhibitors (EgBBI1 and 2), defensin (EgDFS), dehydrin (EgDHN), early methionine-labeled polypeptides (EgEMLP1 and 2), glycine-rich RNA binding protein (EgGRRBP), isoflavone reductase (EgIFR), metallothionein-like protein (EgMT), pathogenesis-related-1 protein (EgPRP), and type 2 ribosome-inactivating protein (EgT2RIP). The transcript abundance of EgBBI2 increased in G. boninense-treated roots at 3 and 6wpi compared to those of controls; while the transcript abundance of EgBBI1, EgDFS, EgEMLP1, EgMT, and EgT2RIP increased in G. boninense-treated roots at 6 or 12wpi. Meanwhile, the gene expression of EgDHN was up-regulated at all three time points in G. boninense-treated roots. The expression profiles of the eleven transcripts were also studied in leaf samples upon inoculation of G. boninense and Trichoderma harzianum to identify potential biomarkers for early detection of BSR. Two candidate genes (EgEMLP1 and EgMT) that have different profiles in G. boninense-treated leaves compared to those infected by T. harzianum may have the potential to be developed as biomarkers for early detection of G. boninense infection.
Polygonum minus is an aromatic plant, which contains high abundance of terpenoids, especially the sesquiterpenes C15H24. Sesquiterpenes were believed to contribute to the many useful biological properties in plants. This study aimed to functionally characterize a full length sesquiterpene synthase gene from P. minus. P. minus sesquiterpene synthase (PmSTS) has a complete open reading frame (ORF) of 1689 base pairs encoding a 562 amino acid protein. Similar to other sesquiterpene synthases, PmSTS has two large domains: the N-terminal domain and the C-terminal metal-binding domain. It also consists of three conserved motifs: the DDXXD, NSE/DTE, and RXR. A three-dimensional protein model for PmSTS built clearly distinguished the two main domains, where conserved motifs were highlighted. We also constructed a phylogenetic tree, which showed that PmSTS belongs to the angiosperm sesquiterpene synthase subfamily Tps-a. To examine the function of PmSTS, we expressed this gene in Arabidopsis thaliana. Two transgenic lines, designated as OE3 and OE7, were further characterized, both molecularly and functionally. The transgenic plants demonstrated smaller basal rosette leaves, shorter and fewer flowering stems, and fewer seeds compared to wild type plants. Gas chromatography-mass spectrometry analysis of the transgenic plants showed that PmSTS was responsible for the production of β -sesquiphellandrene.
Vanda Mimi Palmer (VMP), an orchid hybrid of Vanda tesselata and Vanda Tan Chay Yan is a highly scented tropical orchid which blooms all year round. Previous studies revealed that VMP produces a variety of isoprenoid volatiles during daylight. Isoprenoids are well known to contribute significantly to the scent of most fragrant plants. They are a large group of secondary metabolites which may possess valuable characteristics such as flavor, fragrance and toxicity and are produced via two pathways, the mevalonate (MVA) pathway or/and the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway. In this study, a sesquiterpene synthase gene denoted VMPSTS, previously isolated from a floral cDNA library of VMP was cloned and expressed in Lactococcus lactis to characterize the functionality of the protein. L. lactis, a food grade bacterium which utilizes the mevalonate pathway for isoprenoid production was found to be a suitable host for the characterization of plant terpene synthases. Through recombinant expression of VMPSTS, it was revealed that VMPSTS produced multiple sesquiterpenes and germacrene D dominates its profile.
Transgenic technologies have been applied to a wide range of biological research. However, information on the potential epigenetic effects of transgenic technology is still lacking. Here, we show that the transgenic process can simultaneously induce both genetic and epigenetic changes in rice. We analyzed genetic, epigenetic, and phenotypic changes in plants subjected to tissue culture regeneration, using transgenic lines expressing the same coding sequence from two different promoters in transgenic lines of two rice cultivars: Wuyunjing7 (WYJ7) and Nipponbare (NP). We determined the expression of OsNAR2.1 in two overexpression lines generated from the two cultivars, and in the RNA interference (RNAi) OsNAR2.1 line in NP. DNA methylation analyses were performed on wild-type cultivars (WYJ7 and NP), regenerated lines (CK, T0 plants), segregation-derived wild-type from pOsNAR2.1-OsNAR2.1 (SDWT), pOsNAR2.1-OsNAR2.1, pUbi-OsNAR2.1, and RNAi lines. Interestingly, we observed global methylation decreased in the T0 regenerated line of WYJ7 (CK-WJY7) and pOsNAR2.1-OsNAR2.1 lines but increased in pUbi-OsNAR2.1 and RNAi lines of NP. Furthermore, the methylation pattern in SDWT returned to the WYJ7 level after four generations. Phenotypic changes were detected in all the generated lines except for SDWT. Global methylation was found to decrease by 13% in pOsNAR2.1-OsNAR2.1 with an increase in plant height of 4.69% compared with WYJ7, and increased by 18% in pUbi-OsNAR2.1 with an increase of 17.36% in plant height compared with NP. This suggests an absence of a necessary link between global methylation and the phenotype of transgenic plants with OsNAR2.1 gene over-expression. However, epigenetic changes can influence phenotype during tissue culture, as seen in the massive methylation in CK-WYJ7, T0 regenerated lines, resulting in decreased plant height compared with the wild-type, in the absence of a transformed gene. We conclude that in the transgenic lines the phenotype is mainly determined by the nature and function of the transgene after four generations of transformation, while the global epigenetic modification is dependent on the genetic background. Our research suggests an innovative insight in explaining the reason behind the occurrence of transgenic plants with random and undesirable phenotypes.
Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca(2+) transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca(2+) transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.
Oil palm is one of the most productive oil-producing crops and can store up to 90% oil in its fruit mesocarp. Oil palm fruit is a sessile drupe consisting of a fleshy mesocarp from which palm oil is extracted. Biochemical changes in the mesocarp cell walls, polyamines, and hormones at different ripening stages of oil palm fruits were studied, and the relationship between the structural and the biochemical metabolism of oil palm fruits during ripening is discussed. Time-course analysis of the changes in expression of polyamines, hormones, and cell-wall-related genes and metabolites provided insights into the complex processes and interactions involved in fruit development. Overall, a strong reduction in auxin-responsive gene expression was observed from 18 to 22 weeks after pollination. High polyamine concentrations coincided with fruit enlargement during lipid accumulation and latter stages of maturation. The trend of abscisic acid (ABA) concentration was concordant with GA₄ but opposite to the GA₃ profile such that as ABA levels increase the resulting elevated ABA/GA₃ ratio clearly coincides with maturation. Polygalacturonase, expansin, and actin gene expressions were also observed to increase during fruit maturation. The identification of the master regulators of these coordinated processes may allow screening for oil palm variants with altered ripening profiles.
Chlorella is a popular microalga with robust physiological and biochemical characteristics, which can be cultured under various conditions. The exploration of the small RNA content of Chlorella could improve strategies for the enhancement of metabolite production from this microalga. In this study, stress was introduced to the Chlorella sorokiniana culture to produce high-value metabolites such as carotenoids and phenolic content. The small RNA transcriptome of C. sorokiniana was sequenced, focusing on microRNA (miRNA) content. From the analysis, 98 miRNAs were identified in cultures subjected to normal and stress conditions. The functional analysis result showed that the miRNA targets found were most often involved in the biosynthesis of secondary metabolites, followed by protein metabolism, cell cycle, and porphyrin and chlorophyll metabolism. Furthermore, the biosynthesis of secondary metabolites such as carotenoids, terpenoids, and lipids was found mostly in stress conditions. These results may help to improve our understanding of regulatory mechanisms of miRNA in the biological and metabolic process of Chlorella species. It is important and timely to determine the true potential of this microalga species and to support the potential for genetic engineering of microalgae as they receive increasing focus for their development as an alternative source of biofuel, food, and health supplements.
Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.
Bambara groundnut (Vigna subterranea (L.) Verdc.) is an African legume and is a promising underutilized crop with good seed nutritional values. Low temperature stress in a number of African countries at night, such as Botswana, can effect the growth and development of bambara groundnut, leading to losses in potential crop yield. Therefore, in this study we developed a computational pipeline to identify and analyze the genes and gene modules associated with low temperature stress responses in bambara groundnut using the cross-species microarray technique (as bambara groundnut has no microarray chip) coupled with network-based analysis. Analyses of the bambara groundnut transcriptome using cross-species gene expression data resulted in the identification of 375 and 659 differentially expressed genes (p<0.01) under the sub-optimal (23°C) and very sub-optimal (18°C) temperatures, respectively, of which 110 genes are commonly shared between the two stress conditions. The construction of a Highest Reciprocal Rank-based gene co-expression network, followed by its partition using a Heuristic Cluster Chiseling Algorithm resulted in 6 and 7 gene modules in sub-optimal and very sub-optimal temperature stresses being identified, respectively. Modules of sub-optimal temperature stress are principally enriched with carbohydrate and lipid metabolic processes, while most of the modules of very sub-optimal temperature stress are significantly enriched with responses to stimuli and various metabolic processes. Several transcription factors (from MYB, NAC, WRKY, WHIRLY & GATA classes) that may regulate the downstream genes involved in response to stimulus in order for the plant to withstand very sub-optimal temperature stress were highlighted. The identified gene modules could be useful in breeding for low-temperature stress tolerant bambara groundnut varieties.
Basal stem rot (BSR) of oil palm is a disastrous disease caused by a white-rot fungus Ganoderma boninense Pat. Non-ribosomal peptides (NRPs) synthesized by non-ribosomal peptide synthetases (NRPSs) are a group of secondary metabolites that act as fungal virulent factors during pathogenesis in the host. In this study, we aimed to isolate NRPS gene of G. boninense strain UPMGB001 and investigate the role of this gene during G. boninense-oil palm interaction. The isolated NRPS DNA fragment of 8322 bp was used to predict the putative peptide sequence of different domains and showed similarity with G. sinense (85%) at conserved motifs of three main NRPS domains. Phylogenetic analysis of NRPS peptide sequences demonstrated that NRPS of G. boninense belongs to the type VI siderophore family. The roots of 6-month-old oil palm seedlings were artificially inoculated for studying NRPS gene expression and disease severity in the greenhouse. The correlation between high disease severity (50%) and high expression (67-fold) of G. boninense NRPS gene at 4 months after inoculation and above indicated that this gene played a significant role in the advancement of BSR disease. Overall, these findings increase our knowledge on the gene structure of NRPS in G. boninense and its involvement in BSR pathogenesis as an effector gene.
Hevea brasiliensis is the most widely cultivated species for commercial production of natural rubber (cis-polyisoprene). In this study, 10,040 expressed sequence tags (ESTs) were generated from the latex of the rubber tree, which represents the cytoplasmic content of a single cell type, in order to analyse the latex transcription profile with emphasis on rubber biosynthesis-related genes. A total of 3,441 unique transcripts (UTs) were obtained after quality editing and assembly of EST sequences. Functional classification of UTs according to the Gene Ontology convention showed that 73.8% were related to genes of unknown function. Among highly expressed ESTs, a significant proportion encoded proteins related to rubber biosynthesis and stress or defence responses. Sequences encoding rubber particle membrane proteins (RPMPs) belonging to three protein families accounted for 12% of the ESTs. Characterization of these ESTs revealed nine RPMP variants (7.9-27 kDa) including the 14 kDa REF (rubber elongation factor) and 22 kDa SRPP (small rubber particle protein). The expression of multiple RPMP isoforms in latex was shown using antibodies against REF and SRPP. Both EST and quantitative reverse transcription-PCR (QRT-PCR) analyses demonstrated REF and SRPP to be the most abundant transcripts in latex. Besides rubber biosynthesis, comparative sequence analysis showed that the RPMPs are highly similar to sequences in the plant kingdom having stress-related functions. Implications of the RPMP function in cis-polyisoprene biosynthesis in the context of transcript abundance and differential gene expression are discussed.
We have characterized an oil palm (Elaeis guineensis Jacq.) constitutive promoter that is derived from a translationally control tumor protein (TCTP) gene. The TCTP promoter was fused transcriptionally with the gusA reporter gene and transferred to monocot and dicot systems in order to study its regulatory role in a transient expression study. It was found that the 5' region of TCTP was capable of driving the gusA expression in all the oil palm tissues tested, including immature embryo, embryogenic callus, embryoid, young leaflet from mature palm, green leaf, mesocarp and stem. It could also be used in dicot systems as it was also capable of driving gusA expression in tobacco leaves. The results indicate that the TCTP promoter could be used for the production of recombinant proteins that require constitutive expression in the plant system.
Oil palm suspension cultures were initiated by transferring the gel-like friable embryogenic tissue onto liquid medium supplemented with auxins. In this study, transcripts that were differentially expressed in oil palm suspension cells cultured at different auxin concentrations were examined using suppression subtractive hybridization. Total RNA was first isolated from oil palm suspension cells proliferated in liquid medium with different hormone concentrations for 6 months. Four different hormone combinations: T1 (0.1 mg/l 2,4-D and 1.0 mg/l NAA), T2 (0.4 mg/l 2,4-D and 1.0 mg/l NAA), T3 (1.0 mg/l NAA), and T4 (0.4 mg/l 2,4-D) were used for the treatments. The first and second subtractions were performed using samples T1 and T2 in forward and reverse order. The other two subtractions were forward and reverse subtractions of T3 and T4, respectively. Reverse northern analyses showed that 14.13% of these clones were preferentially expressed in T1, 13.70% in T2, 14.75% in T3, and 15.70% in T4. Among the 294 cDNA clones that were sequenced, 61 contigs (assembled from 165 sequences) and 129 singletons were obtained. Among the 61 contigs, 10 contigs consist of sequences from treatment T1, 8 contigs were from treatment T2, 10 contigs were contains sequences of treatment T3 and 13 contigs contains sequences of treatment T4. Northern analyses of five transcripts that were shown to be differentially expressed in the oil palm suspension cells by reverse northern analysis revealed that transcripts 16A1 (a putative lignostilbene-alpha,beta-dioxygenase, EgLSD) and 16H12 (a putative ethylene responsive 6, EgER6) were differentially expressed in oil palm suspension cells treated with different levels of auxin.
Abscisic acid (ABA) is an important phytohormone involved in the abiotic stress resistance in plants. The ABA-responsive element (ABRE) binding factors play significant roles in the plant development and response to abiotic stresses, but none so far have been isolated and characterized from the oil palm. Two ABA-responsive cDNA clones, named EABF and EABF1, were isolated from the oil palm fruits using yeast one-hybrid system. The EABF had a conserved AP2/EREBP DNA-binding domain (DNA-BD) and a potential nuclear localization sequence (NLS). No previously known DNA-BD was identified from the EABF1 sequence. The EABF and EABF1 proteins were classified as DREB/CBF and bZIP family members based on the multiple sequence alignment and phylogenetic analysis. Both proteins showed ABRE-binding and transcriptional activation properties in yeast. Furthermore, both proteins were able to trans-activate the down-stream expression of the LacZ reporter gene in yeast. An electrophoretic mobility shift assay revealed that in addition to the ABRE sequence, both proteins could bind to the DRE sequence as well. Transcriptional analysis revealed that the expression of EABF was induced in response to the ABA in the oil palm fruits and leaves, but not in roots, while the EABF1 was constitutively induced in all tissues. The expressions of both genes were strongly induced in fruits in response to the ABA, ethylene, methyl jasmonate, drought, cold and high-salinity treatments, indicating that the EABF and EABF1 might act as connectors among different stress signal transduction pathways. Our results indicate that the EABF and EABF1 are novel stress-responsive transcription factors, which are involved in the abiotic stress response and ABA signaling in the oil palm and could be used for production of stress-tolerant transgenic crops.
Silicon (Si) is the second most abundant element in soil after oxygen. It is not an essential element for plant growth and formation but plays an important role in increasing plant tolerance towards different kinds of abiotic and biotic stresses. The molecular mechanism of Si absorption and accumulation may differ between plants, such as monocotyledons and dicotyledons. Silicon absorption and accumulation in mangrove plants are affected indirectly by some proteins rich in serine and proline amino acids. The expression level of the genes responsible for Si absorption varies in different parts of plants. In this study, Si is mainly observed in the epidermal roots' cell walls of mangrove plants compared to other parts. The present work was carried out to discover further information on Si stress responsive genes in Rhizophora apiculata, using the suppression subtractive hybridization technique. To construct the cDNA library, two-month-old seedlings were exposed to 0.5, 1, and 1.5 mM SiO2 for 15 hrs and for 1 to 6 days resulting in a total of 360 high quality ESTs gained. Further examination by RT-PCR and real-time qRT-PCR showed the expression of a candidate gene of serine-rich protein.