Transcriptomic analysis of differentially expressed genes during anther development in genetic male sterile and wild type cotton by digital gene-expression profiling
© Wei et al.; licensee BioMed Central Ltd. 2013
Received: 2 July 2012
Accepted: 1 February 2013
Published: 12 February 2013
Cotton (Gossypium hirsutum) anther development involves a diverse range of gene interactions between sporophytic and gametophytic tissues. However, only a small number of genes are known to be specifically involved in this developmental process and the molecular mechanism of the genetic male sterility (GMS) is still poorly understand. To fully explore the global gene expression during cotton anther development and identify genes related to male sterility, a digital gene expression (DGE) analysis was adopted.
Six DGE libraries were constructed from the cotton anthers of the wild type (WT) and GMS mutant (in the WT background) in three stages of anther development, resulting in 21,503 to 37,352 genes detected in WT and GMS mutant anthers. Compared with the fertile isogenic WT, 9,595 (30% of the expressed genes), 10,407 (25%), and 3,139 (10%) genes were differentially expressed at the meiosis, tetrad, and uninucleate microspore stages of GMS mutant anthers, respectively. Using both DGE experiments and real-time quantitative RT-PCR, the expression of many key genes required for anther development were suppressed in the meiosis stage and the uninucleate microspore stage in anthers of the mutant, but these genes were activated in the tetrad stage of anthers in the mutant. These genes were associated predominantly with hormone synthesis, sucrose and starch metabolism, the pentose phosphate pathway, glycolysis, flavonoid metabolism, and histone protein synthesis. In addition, several genes that participate in DNA methylation, cell wall loosening, programmed cell death, and reactive oxygen species generation/scavenging were activated during the three anther developmental stages in the mutant.
Compared to the same anther developmental stage of the WT, many key genes involved in various aspects of anther development show a reverse gene expression pattern in the GMS mutant, which indicates that diverse gene regulation pathways are involved in the GMS mutant anther development. These findings provide the first insights into the mechanism that leads to genetic male sterility in cotton and contributes to a better understanding of the regulatory network involved in anther development in cotton.
Male sterility is a simple and efficient pollination control system that is widely exploited in hybrid breeding. In cotton breeding, cytoplasmic male sterility (CMS) and genetic male sterility (GMS) have been used to produce hybrid seeds. Both types of lines have a maternally (former) or nuclear (later) inherited trait and each line is unable to produce or release functional pollen. Such lines are suitable as maternal plants for the utilization of hybrid vigor.
Although the CMS system can economically generate a completely male-sterile population, hybrid seed production using this system involves development of three lines, it usually takes years to develop the maintainer and restorer lines because most CMS systems have stringent restoring-maintaining relationships, and any CMS system potentially causes negative cytoplasmic effects and shows unstable sterility . The advantage of GMS system is that it involves only two lines and the genes responsible for male sterility are relatively easy to transfer to any desired genetic background. For example, in genetic male sterility controlled by a recessive gene(s) (RGMS), most breeding lines can serve as restorers, thus it is easy to combine elite lines to produce hybrids that show high heterosis. The molecular mechanism of GMS is currently a research hotspot in plant science.
GMS-sterile lines are mainly controlled by alleles of nuclear genes designated as ‘ms’ that affect male reproduction. These alleles are usually recessive, but some are dominant (Ms) and both types are typically expressed in specific sporophytic tissues at different stages. In G. hirsutum, the GMS mutant in ‘Dong A’ is controlled by one pair of recessive genes , and it has the same genetic background with its wild type (WT). Therefore, they are ideal genetic materials for studying cotton anther development and male sterility.
In plants, anther is a bilaterally symmetrical structure with four lobes, each lobe containing meiotic cells at the center surrounded by four somatic cell layers, which are the epidermis, the endothecium, the middle layer, and the tapetum from the surface to the interior . The meiotic cells produce microspore in anther lodes. Among the four somatic cell layers, the tapetum cells are of considerable physiological significance because all nutritional materials entering the sporogenous cells passes through or originates from the tapetum . Tapetal aberrations are frequently observed in male sterile mutants, with premature or delayed degradation of the tapetum resulting in male sterility .
The anther formation and maturation is a critical phase in the plant life cycle, which commences at the end of meiosis with the formation of a tetrad and ends at the dehiscence of anthers when the mature pollen grains are released. This process involves a diverse range of gene interactions between sporophytic and gametophytic tissues [6, 7]. Using forward and reverse genetic approaches, a growing number of genes have been identified to have vital roles in anther development . For instance, SPOROCYTELESS (SPL)/NOZZLE (NZZ) gene required for cell division and differentiation is essential for the earliest anther development in Arabidopsis[9, 10]. ABORTED MICROSPORES (AMS) is involved in tapetal development and microspore development in rice and Arabidopsis. TAPETAL DETERMINANT1 (TPD1), DYSFUNCTIONAL TAPETUM (DYT1), and UNDEVELOPED TAPETUM1 are involved in tapetal and microsporocyte determination [12–14]. MALE STERILITY1 (MS1) is required for tapetal development and pollen wall formation in Arabidopsis. MALE STERILITY2 (MS2) is involved in pollen wall biosynthesis in rice . Cytochrome P450 family member CYP704B1 and CYP704B2 are essential for the formation of both cuticle and extine during plant male reproductive and spore development in Arabidopsis and rice [17, 18]. LAP5 and LAP6 encode anther-specific proteins with similarity to chalcone synthase essential for pollen extine development in Arabidopsis . However, the regulatory mechanism underlying anther development as well as the genetic network that governs male sterility is yet not to be fully understood.
In recent decades, most studies on male sterility focused on differential display methods, such as amplified fragment length polymorphisms [20, 21]. These low-throughput methods are inadequate to fully detect gene expression in different samples. Although transcriptome profiling studies based on microarray data can detect detailed gene expression in different samples, this approach also has limitations because genes are represented by unspecific probe sets and cannot be reliably detected at low expression levels. High-throughput tag-sequencing for digital gene expression (DGE) analysis is a powerful, recently developed tool that allows the concomitant sequencing of millions of signatures to the genome, and identification of specific genes and the abundance of gene expression in a sample tissue . This method provides a more qualitative and quantitative description of gene expression than previous microarray-based assays . In a direct comparison with high-throughput mRNA sequencing (RNA-seq), both methods provide similar assessments of relative transcript abundance, but DGE better detects expression differences for poorly expressed genes and does not exhibit transcript length bias .
In the present study, DGE based on the Solexa Genome Analyzer platform was applied to analyze the complex regulatory network underlying anther development and to investigate the differences in gene expression between the ‘Dong A’ male sterile mutant and its wild type (WT). After high-throughput tag-sequencing, an integrated bioinformatic analysis was performed to identify the expression patterns of genes and critical pathways in the male-sterile mutant and WT at three stages of male gametophyte development. The aim of the study was to gain increased insight into the molecular mechanisms of cotton male sterility. The results yielded sets of up-regulated and down-regulated genes associated with male sterility, and candidate genes associated with cotton male sterility are discussed. Comparison of the gene expression patterns at three developmental stages of the WT and GMS mutant anther provides an improved understanding of the molecular mechanisms of cotton anther development and GMS.
Phenotypic analysis of impaired anthers in cotton male-sterile mutant
From meiosis to the uninucleate microspore stage, the anther epidermal cell walls of the WT were strongly stained and gradually thickened, and decreased in electron density; the anther cuticle also became thicker (Figure 2M, N and O). Surprisingly, although the epidermis of the GMS mutant anther developed a cell wall, the wall was weakly stained, which indicated that decreased amounts of lipophilic materials were deposited in the outer epidermal cell wall. Moreover, the thickness of the cell wall and cuticle were greatly reduced in the GMS mutant (Figure 2P, Q and R). Thus, abnormal synthesis of the pollen wall and defective formation of the anther epidermal cell wall and cuticle in the GMS mutant anthers were indicated.
Analysis of digital gene expression (DGE) libraries
Previous studies have demonstrated that the peak of male sterility mainly occurs in the uninucleate microspore stage of anthers in ‘Dong A’ GMS mutant . Therefore, not only the uninucleate microspore stage anthers, but also anthers before the uninucleate microspore stage should be collected to analyze the differentially expressed genes between the WT and the GMS mutant anther. Based on this fact, anthers were harvested at the meiosis stage (WT: F-1; mutant: S-1), tetrad stage (WT: F-2; mutant: S-2), and uninucleate microspore stage (WT: F-3; mutant: S-3) to construct six DGE libraries. The Solexa Genome Analyzer platform was used to perform high-throughput tag-sequencing analysis of cotton anther libraries. The transcriptome during cotton anther development was characterized and differences in the gene regulatory pathways were analyzed at the three stages of WT and GMS mutant anther development.
Categorization and abundance of tags
Distinct tag numbers
All tag mapping to gene
Total % of clean tags
Distinct tag numbers
Distinct tag % of clean tags
Unambiguous Tag mapping to gene
Total % of clean tags
Distinct tag numbers
Distinct Tag % of clean tags
All tag-mapped genes
% of ref genes
Unambiguous tag-mapped genes
% of ref genes
Total % of clean tags
Distinct tag numbers
Distinct tag % of clean tags
Heterogeneity and redundancy are two significant characteristics of mRNA expression. A small number of mRNAs are present at a very high abundance, whereas expression of the majority of mRNAs remains at a very low level. The distribution of clean tag expression can be used to evaluate the normality of the whole DGE data. In the present study, the distribution of the distinct clean tag copy numbers showed extremely similar tendencies. Among the distinct clean tags in the six libraries, only 3.2–5.2% possessed more than 100 copies. The majority of distinct clean tags (28.64–60.88%) had 2–5 copies (Additional file 1: Figure S1), which indicated that the whole DGE data among the six libraries was normally distributed.
Analysis of tag mapping
A reference gene database that included 116, 520 sequences of the G. hirsutum unigenes was preprocessed for tag mapping (http://occams.dfci.harvard.edu/pub/bio/tgi/data/Gossypium). Among the reference sequences, the genes with a CATG site accounted for 87.52% of the total. To obtain the reference tags, all of the CATG +17 tags in the gene were taken as gene reference tags. Finally, 284,576 total reference tag sequences with 236,167 unambiguous reference tags were obtained in the six libraries. On this basis, 44,968 unambiguous tag-mapped genes were detected in the six DGE libraries.
To estimate whether the sequencing depth was sufficient for the transcriptome coverage, the sequencing saturation in the six libraries was analyzed. The genes that were mapped by all clean tags and unambiguous clean tags increased with the total number of tags. However, when the sequencing counts reached 2 million tags or higher, the number of detected genes was saturated (Additional file 2). Given that Solexa sequencing can distinguish transcripts that originate from both DNA strands, using the strand-specific nature of the sequencing tags obtained, we found that 26.25–33.96% distinct clean tags showed a perfect match to sense strand-specific transcripts, and 7.43–12.98% distinct clean tags showed a perfect match to antisense strand-specific transcripts in the six libraries (Additional file 3). These results indicated that antisense genes also play important roles in the transcriptional regulation of cotton anther development.
Transcriptome diversity quantified by DGE profiles in the WT and GMS mutant anthers
The fertile anthers exhibit transcript complexity during three anther stages: of a total of 42,488 genes that were expressed over three stages of WT anthers, 22,131 were constitutively expressed, 13,269 were stage-specific and 7,088 were expressed at two stages (Figure 3b; see Additional file 4 for the gene list for each category). These dynamic patterns of gene expression reinforce the conclusion that male reproductive development is a highly complex process in plants. Strikingly, during the three stages of WT anthers, the transcriptome in the tetrad stage were significantly more than other two stages. However, the transcriptome in the tetrad stage is lower than the same two stages in maize anther development . This reversely increased transcript diversity in the tetrad stage of cotton anthers indicates that cotton anther development may have some unique gene regulation processes and pathways.
During the transition from the meiosis to tetrad stage, there was a major transcriptome change: 10,022 stage-specific transcripts expressed in the tetrad stage anther (fufous checked bar in Figure 3d, see Additional file 4) and 3,764 transcripts expressed in the meiosis stage anther were not detectable at the tetrad stage (green checked bar in Figure 3d, see Additional file 4). Notably, the stage-specific transcripts in the tetrad stage of WT anthers were also significantly more than other two stages (Figure 3b), indicated that the anther in the tetrad stage might express a distinctive suite of genes for core cellular functions, which may be vital for tetrad generation or microspore separation during the cotton anther development. During the subsequent transition from the tetrad to uninucleate microspore stage, 1,372 transcripts were specific (fufous checked bar in Figure 3d), while 2,525 transcripts (blue bar in Figure 3d, see Additional file 4) that were present at both the meiosis and tetrad microspore stages were no longer detectable.
Although mutant anthers were defective, the transcriptome of mutant anther was also highly dynamic. Of a total of 34,932 genes that are expressed over the three stages of mutant anthers, 16,145 transcripts were shared across all the three mutant anther stages, 9,737 were stage-specific transcripts and 9,050 genes were expressed at two stages (Figure 3c, see Additional file 5). During the meiosis to tetrad stage transition, 4,376 new transcripts appeared (fufous checked bar in Figure 3d) and 3,794 transcripts that were expressed in the meiosis stage anther were not detectable at the tetrad stage (green checked bar in Figure 3d, see Additional file 5). At the subsequent tetrad to uninucleate microspore transition, 2,305 transcripts appeared (fufous checked bar in Figure 3d), and 1,564 transcripts (blue bar in Figure 3d, see Additional file 5) present at both the meiosis and tetrad microspore stages were no longer detectable.
In comparison with the three stages of WT anthers, transcripts have differentially degrees of reduction in three corresponding stages of mutant anthers (Figure 3a). Especially, the number of stage-specific transcripts decreased more than twofold in the tetrad stage of mutant anthers (10,022 and 4,627 stage-specific transcripts detected in the tetrad stage of WT and mutant anthers, Figure 3b and c). Considering the observations of tapetal cells underwent acute and abnormal degradation in the tetrad stage of mutant anthers, so many stage-specific transcripts missing in the tetrad stage of mutant anthers were likely to disrupt the temporally coordinated growth and differentiation of anther cells, resulting in the observed defects in tapetal cell division.
Genes associated with male sterility in GMS mutant anthers
To characterize the functional consequences of gene expression changes associated with male sterility, a pathway analysis of the differentially expressed genes based on the KEGG database was performed. The results indicated that several key branch-point genes that participate in histone modification and DNA methylation, hormone signaling, sucrose and starch metabolism, the pentose phosphate pathway, glycolysis, flavonoid metabolism, programmed cell death (PCD), cell wall loosening, and reactive oxygen species (ROS) generation/scavenging, were differentially expressed in GMS mutant anthers (Additional file 7).
Selected genes with altered expression in anthers of the GMS mutant
TPM (transcript per million clean tag)
Histone modification and DNA methylation
Gibberellic acid receptor
Gibberellin 3-hydroxylase 1
Ethylene-responsive transcription factor
ACC oxidase 4
Sucrose and starch metabolism
Sucrose transporter 1
alpha-1,4-glucan phosphorylase L isozyme
Granule bound starch synthase
Hexose transporter 6
Starch branching enzyme II-1
Pentose phosphate pathway
Fructose bisphosphate aldolase
Leucoanthocyanidin reductase 1
Program cell death (PCD)
Cysteine proteinase inhibitor
Cell wall development
Actin depolymerizing factor 1
Actin depolymerizing factor 2
Reactive oxygen species (ROS) generation/scavenging
Class III peroxidase
Cytosolic ascorbate peroxidase 1
Catalase isozyme 1
Several genes implicated in ethylene synthesis and GA signaling, the expression of which changed significantly in the GMS mutant anther, were identified. These differentially expressed genes included those that encode an ethylene receptor, ethylene-responsive transcription factor, ACC oxidase 4, and GA receptor. These genes were also distinctly up-regulated at the tetrad stage in GMS mutant anthers, but exhibited lower expression at the uninucleate microspore stage in GMS mutant anthers (Table 2).
Changes in expression of carbon synthesis genes were apparent in the GMS mutant anthers. Starch and sucrose production declined as transcript levels of the starch-branching enzyme and sucrose synthetase decreased. Most genes that participated in starch and sucrose metabolism showed much lower expression at the meiosis and uninucleate microspore stages, but were up-regulated at the tetrad stage in GMS mutant anthers. Notably, we detected much higher expression of sugar transferase at two anther developmental stages in the GMS mutant, especially at the meiosis stage, and expression of sugar transferase was twenty-fold higher than that at the same stages of WT anther development. These findings indicated that decreased synthesis and transporting of sucrose and starch occurred in GMS mutant anthers.
Similar to most of the genes involved in starch and sucrose metabolism, many genes associated with glycolysis and the pentose phosphate pathway were up-regulated at the tetrad stage, but down-regulated at the meiosis and uninucleate microspore stages, of GMS mutant anthers. These genes included those that encoded hexokinase, phosphofructokinase, phosphoglycerate kinase, pyruvate kinase, and alcohol dehydrogenase, as well as genes that encoded fructose-bisphosphate aldolase and 6-phosphogluconate dehydrogenase.
Genes involved in flavonoid metabolism also have the same expression patterns in the GMS mutant anthers, such as chalcone synthase and flavonoid 3′, 5′-hydroxylase, anthocyanidin reductase, and leucoanthocyanidin reductase were down-regulated at the meiosis and uninucleate microspore stages, but were up-regulated at the tetrad stage, in GMS mutant anthers (Table 2). However, genes involved in cell wall loosening, such as a beta-galactosidase, katanin-like protein, actin depolymerizing factor 1, and actin depolymerizing factor 2, were up-regulated at the three developmental stages in the GMS mutant anthers (Table 2). These observations indicated that the actin cytoskeleton balance may be disturbed in the GMS mutant anther, which would directly affect pollen cell wall development.
In active oxygen metabolism and PCD pathways, expression of several ROS generation/scavenging-related genes and PCD-related genes were changed in the GMS mutant anthers. Up-regulated genes included those that encode cysteine proteinase and cytochrome c, and down-regulated genes included those that encode a cysteine proteinase inhibitor, peroxidase precursor, class III peroxidise, Cu/Zn superoxide dismutase (SOD), and catalase isozyme, at the uninucleate microspore stage of GMS mutant anthers (Table 2).
Tag-mapped genes were confirmed by qRT-PCR
In flowering plants, the production of functional pollen grains is a prerequisite for the propagation, which relies on cooperative functional interactions between gametophytic and sporophytic tissues within the anther. As a non-photosynthetic male reproductive organ, anther must obtain photosynthetic assimilates predominantly from source organs to support pollen development and maturation . The sink strength of anthers is highest at early stages of flower development, with the anthers requiring large amounts of sugar to support their development . At advanced stages, pollen maturation requires the accumulation of starch, which is converted from sugar and functions as an energy reserve for germination ; thus, the presence of sufficient levels of sucrose is of vital importance for growth of the male reproductive cells in plants. Disturbances in sugar metabolism and unloading in the anther can significantly impair pollen development and cause male sterility .
To ensure reproductive success, flowering plants established two barriers for protecting the pollen grains from various environmental stresses, such as bacterial and fungal attacks. One barrier is the anther wall, which covered by a cuticle can protect the anther from various stresses. Another barrier is the cell wall of the pollen grain, which can be divided into three layers: an outer extine, an inner intine, and a lipid- and protein-rich pollen coat in the crevices of the extine. The extine is a multi-layered structure, primarily consisting of sporopollenin . Sporopollenin is a complex polymer composed of fatty acids and phenolic compounds . The sporopollenin polymer is responsible for the extine’s unparalleled physical strength, chemical inertness, and elasticity to protect pollen against various stresses. After anther dehiscence, the extine forms the outer wall of the pollen grain. These two major barriers are vital for the development of viable pollen grains and male sterility .
More importantly, anther development is apparently influenced by ethylene and gibberellic acid (GA). Ethylene is a key gaseous phytohormone that regulates many aspects of plant growth and development [39, 40]. The mature pollen is characterized by a high content of ethylene. During the anther development, fertile male gametophyte development is accompanied by two peaks of ethylene production by anther tissues . The first peak occurs during microspore development simultaneously with degeneration of both tapetal tissues and the middle layers of the anther wall. The second peak coincides with maturation and dispersal of the pollen grains. GA also plays an essential role in the regulation of reproductive development, especially of anther development. In A. thaliana, GAs can accelerate flowering, even under short-day conditions. In maize (Zea mays) and castor bean (Ricinus communis), GAs can promote the development of female flowers [42, 43]. Conversely, GA deficiency or insensitivity in several plant species causes abnormal development of anthers, which often leads to male sterility .
To explore the key underlying molecular switches resulting in the male sterility of the cotton GMS mutant, six DGE libraries were set up during the anther development of Dong A WT and GMS mutant in the current work. To the best of our knowledge, the present study is the first to attempt to perform a deep sequencing of the Dong A WT and GMS mutant in cotton anther development, which may facilitate identification of systemic gene expressions and regulatory mechanisms for male sterility.
After analyzing the differences in gene expressions of the WT and the GMS mutant during the anther development, several opposite gene expression were discovered to exist in the anther development of GMS mutant. These genes involved in histone modification and DNA methylation, carbon and energy metabolism, hormone signaling, pollen cell wall development, and ROS metabolism. For example, three histone constitution-related genes and two methyltransferase genes were differentially expressed in the GMS mutant anthers (Table 2). Both histone modification and DNA methylation play essential roles in genome management, and control gene expression or silence . As a main enzyme for DNA methylation, DNA methyltransferase is not only associated with DNA methylation, but also link to many important biological activities, including cell proliferation, senescence . In the present study, three histone constitution-related genes encoding histone H3 and histone H2B were up-regulated at the tetrad stage, but down-regulated at the meiosis and uninucleate microspore stages in the GMS mutant anther. The much lower expression of these histone constitution-related genes at the meiosis and uninucleate microspore stage of GMS mutant anther may affect the normal chromosome structure and influence the functional gene expression in the GMS mutant anther. Besides, two methyltransferase genes were gradually up-regulated in the three stages of GMS mutant anthers, but they were gradually decreased in the three corresponding stages of WT anthers (Table 2), indicating that the expression of these methyltransferase genes were activated in the anther development of GMS mutant, but they were suppressed in the anther development of WT. The increased expression of methyltransferas genes in the GMS mutant anther development may change the level of DNA methylation and repress the functional gene expressions in the anther development of GMS mutant, which may explain the results of the down-regulated genes outnumbered up-regulated genes at the three stages of GMS mutant anther development (Figure 4).
The present research also found that expression of many key branch-point genes associated with the sucrose and starch metabolism, the pentose phosphate pathway, and the glycolysis pathways were affected in GMS mutant anthers. All of these genes with potential roles in carbon and energy metabolism were down-regulated at the uninucleate microspore stage in GMS mutant anthers, which indicated that energy metabolism was disturbed or suppressed during the later anther developmental stage of GMS mutant. In facts, pollen requires the accumulation of sucrose and starch as an energy reserve for maturation in later anther developmental stages. These opposite gene expression patterns in the GMS mutant anthers may directly reduce sucrose transportation and accelerate sucrose depletion in anthers, leading to nutrient substance deficiency in the GMS mutant and related to male sterility.
More interestingly, five genes involved in ethylene metabolism were differential expressed in WT and GMS mutant anthers. Previous studies on male sterility noted changes in ethylene and gibberellic acid (GA) content in GMS lines [47, 48]. These genes were up-regulated by almost two-fold at the tetrad stage in the GMS mutant anther, but down-regulated at the meiosis and uninucleate microspore stages in the WT anther. Considering the TEM observations at the tetrad stage in GMS mutant anthers, surprisingly we found that the higher expression of genes involved in ethylene metabolism coincided with the peak in tapetal tissue degeneration in the GMS mutant anther, which indicated that the higher expression of these ethylene metabolism-related genes may directly lead to the premature degeneration of the tapetal layer in GMS mutant anthers. During the uninucleate microspore stage in WT and GMS mutant anthers (the stage of pollen prematuration), these genes showed lower expression in the GMS mutant anther but higher expression in the WT anther. Given that the pollen grains in the GMS mutant are not dispersed, this opposing pattern of ethylene metabolism expression in the GMS mutant anther suggests that ethylene also may act as an important signal mediating the response to tapetal degeneration and pollen maturation.
The present research also found three GA metabolism-related genes (GA receptor, gibberellin 3-hydroxylase 1, and DELLA protein) that were up-regulated at the tetrad stage, but down-regulated at the meiosis and uninucleate microspore stages in the GMS mutant anther (Table 2). Given the essential roles of the GA receptor, and involvement of the DELLA protein in GA signaling for pollen extine production and male sterility in A. thaliana and rice [49, 50], the much lower expression of these GA metabolism-related genes at the uninucleate microspore stage of GMS mutant anther development could affect pollen extine formation.
Interestingly, the expression pattern of these GA metabolism-related genes was similar with the expression patterns of the ethylene metabolism-related genes at the same stages of WT and GMS mutant anther development. Ethylene treatment can markedly affect GA content; however, GA does not significantly affect ethylene-related gene expression in relation to regulation of hypocotyl and root length, which suggests that ethylene acts upstream via GA to regulate hypocotyl and root development . On the basis of the similar expression pattern of the ethylene and GA metabolism-related genes in the WT and GMS mutant anthers, we reasoned that ethylene also may act upstream via GA to regulate cotton anther development, thereby leading to the higher expression of GA metabolism-related genes at the tetrad stage, but lower expression at the meiosis and uninucleate microspore stages in GMS mutant anthers. The opposing expression pattern of these key branch-points genes involved in ethylene and GA metabolism during GMS mutant anther development may act as an important hormone signal-mediating response to male sterility.
More importantly, expression of key branch-point genes involved in flavonoid metabolism, which include chalcone synthase, flavonoid 3′, 5′-hydroxylase, anthocyanidin reductase, and leucoanthocyanidin reductase increased at the tetrad stage but were down-regulated at the uninucleate microspore stage in GMS mutant anthers, which indicated that flavonoid metabolism was initially activated, then suppressed at advanced stages of GMS mutant anther development (Figure 6). In particular, at the uninucleate microspore stage in GMS anthers, the expression of chalcone synthase and flavonoid 3′, 5′-hydroxylase was down-regulated ten-fold and two-fold, respectively. Flavonoids are important for pollen germination and fertility in several plant species, but not for A. thaliana. Chalcone synthase and flavonoid 3′, 5′-hydroxylase are two key branch-point genes for flavonoid biosynthesis, which is vital important for pollen extine formation . Mutation in the key branch-point genes involved in flavonoid metabolism in the former plant species results in male sterility and the absence of flavonoids in the mature stamens. Exogenous application of flavonols can completely restore their fertility . The much lower expression of these key branch-point genes may lead to the excessive reduction of flavonoid synthesis and accumulation of phenolic compounds at the uninucleate microspore stage in GMS mutant anthers, which are considered to be vital for extine formation and related to male sterility in the GMS mutant.
The present research also found that several genes involved in pollen cell wall development were affected in the GMS mutant anther development. For example, at the uninucleate microspore stage, lower expression of the alpha-tubulin and tubulin alpha-2 chain genes in GMS mutant anthers was detected. Reduction of the tubulin level could affect tube elongation and result in aberrant cytoskeletal structures in microspores . In thermo-sensitive genetic male sterile line of wheat, anthers show a shrunken plasma membrane and aberrant actin cytoskeletal structures, which indicate that the cytoskeletal organization is altered during anther development . In contrast, actin depolymerizing factor genes were highly expressed at the three anther developmental stages in the GMS mutant (Figure 6). In G. hirsutum, over-expression of the actin depolymerizing factor (GhADF7) can alter the balance of actin depolymerization and polymerization, thus leading to defective cytokinesis and partial male sterility . Surprisingly, in the present study beta-D-glucosidase (an important cell wall-degrading enzyme) exhibited higher expression at the meiosis and uninucleate microspore stages in the GMS mutant anthers. The abnormal expression of alpha-tubulin, tubulin alpha-2 chain, actin depolymerizing factor, and beta-D-glucosidase genes in the GMS mutant anther, may directly influence the synthesis of F-actin cytoskeleton and further alter the cell wall developmental pattern in the GMS mutant anther development.
Notably, compared with the WT anther, much lower expression of the peroxidase precursor, class III peroxidise, Cu/Zn SOD, and catalase isozyme1 was detected at the uninucleate microspore stage in the GMS mutant anther, which indicated that not only various ROS-scavenging enzymes were suppressed, but also their generation may be blocked in the GMS mutant anther. In plants ROS have prominent roles in the induction, signaling, and execution of PCD . In the cotton CMS line, at the peak of anther abortion, excessive accumulation of ROS and a significant down-regulation of ROS-scavenging enzyme expression coincide with male cell death in male sterility . In present study, the decreased ROS-scavenging enzyme activity in the GMS mutant anther may lead to a transient oxidative burst and significantly increased ROS accumulation at the uninucleate microspore stage in the GMS mutant anther. Furthermore, at the uninucleate microspore stage (the peak stage of anther abortion in the GMS mutant anther), we detected two PCD-related genes (cysteine proteinase and cytochrome c) that were up-regulated in the GMS mutant anther. It is considered that ROS can trigger and promote the expression of cysteine proteinase and cytochrome c, which lead to PCD in plants . These results are consistent with recent findings that male cell death is coincides with excessive formation of ROS at the peak of anther abortion in a cotton CMS line. Surprisingly, a cysteine proteinase inhibitor is down-regulated at the three stages of GMS mutant anther development. On the basis of these results, it is possible that the expression balance of the cysteine proteinase and cysteine proteinase inhibitor might be disturbed in later GMS mutant anther development, thus resulting in the imbalance in oxidative metabolism and male cell death in the GMS mutant.
The present results indicate that expression of genes participated in many diverse molecular functions were altered during anther development in the GMS mutant. Several key branch-point genes involved in histone modification and DNA methylation, hormone signaling, carbon and energy metabolism, pollen wall development, and ROS generation or scavenging were differentially expressed in the GMS mutant anther. To some degree, compared to the WT anther at the same developmental stage, these differentially expressed genes exhibit opposite expression patterns in the GMS mutant anther, which indicates that the hormone signals and energy metabolism are disturbed or blocked during anther development in the GMS mutant. These changes may be the major factors that related to male sterility. These findings provide systemic insights into the mechanism of male sterility and contribute to an improved understanding of the molecular mechanism of male sterility in cotton. Some key branch-point genes involved in cotton anther development are good candidate genes for functional analysis of anther development in the future.
Plants of upland cotton (G. hirsutum) cv. ‘Dong A’ (WT) and the GMS mutant in the ‘Dong A’ background were grown in an experimental field at the China Agricultural Academy of Science Cotton Research Institute under standard field conditions during the spring and summer of 2010. Previous study revealed that when the longitudinal length of buds reach 5.0 mm, 6.5 mm, and 9 mm, respectively, the pollen mother cell of the GMS mutant enter the meiosis, tetrad and uninucleate stages . According to these sampling criterions and combined with microscopic examination of pollen mother cells, developing anthers at these three stages of development were harvested during early morning on the basis of floral bud length. The excised anthers were frozen in liquid nitrogen and stored at −70°C prior to examination. In addition, WT and GMS mutant anthers at the meiosis, tetrad, and uninucleate stages were harvested for total RNA isolation.
Total sugar content measurement
Anthers were harvested and frozen at −70°C. The samples were ground into a powder with a pestle and mortar. Twenty milliliters of water were added to glass tubes containing 1 g of ground anther tissue, incubated at 100°C for 10 min, then centrifuged at 2500 g for 5 min. A total of 2 mL of a solution containing glucose, fructose, or galactose was prepared. To this solution, 2 mL of a glucose solution was added to achieve final concentrations of 0%, 2%, 4%, 6%, 8%, and 10% glucose for optimization. An anthrone colorimetric method was adopted to determine the total sugar content in the WT and GMS mutant anthers .
Sequencing and library construction
Total RNA was extracted from anthers using the pBiozol Total RNA Extraction Reagent (BioFlux) in accordance with the manufacturer’s instructions. RNA was precipitated with ethanol, dissolved in diethypyrocarbonate-treated water (DEPC) and stored at −70°C. All RNA samples were examined for protein contamination (as indicated by the A260/A280 ratio) and reagent contamination (indicated by the A260/A230 ratio) with a Nanodrop ND 1000 spectrophotometer (NanoDrop, Wilmington, DE).
Total RNA purity and degradation were checked with a 1% agarose gel before proceeding. The samples for transcriptomic analysis were prepared using the Illumina kit following the manufacturer’s recommendations. The extracted total RNAs were resolved on a denatured 15% polyacrylamide gel. Briefly, mRNA was purified from 6 μg total RNA using oligo (dT) magnetic beads. Following purification, the mRNA was fragmented into small pieces using divalent cations under an elevated temperature and the cleaved RNA fragments were used for first-strand cDNA synthesis using reverse transcriptase and random primers.
Quantitative RT-PCR analysis
Quantitative RT-PCR analysis was used to verify the DGE results. The RNA samples used for the qRT-PCR assays were identical to those used for the DGE experiments. Gene-specific primers were designed on the basis of the reference unigene sequences with Primer Premier 5.0 (see Additional file 8). The qRT-PCR assay was performed in accordance with the manufacturer’s specifications. The reactions were incubated in a 96-well plate at 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 60 s. The cotton 18S RNA gene (forward primer: 5′-ATCAGCTCGCGTTGACTACGT-3′; reverse primer: 5′- ACACTTCACCGGACCATTCAAT-3′) was used to normalize the amount of gene-specific RT-PCR products, and the relative expression levels of genes were calculated with the 2−ΔΔCt method.
The anthers wish to thank the National Basic Research Program of China (grant no. 2010CB126006) and the 863 Project of China (grant no. 2011AA10A102) for the financial support provided to this project.
- Schnable PS, Wise RP: The molecular basis of cytoplasmic male sterility and fertility restoration. Trends Plant Sci. 1998, 3: 175-180. 10.1016/S1360-1385(98)01235-7.Google Scholar
- Sheng TZ: Thesises on male sterile of cotton. 1989, Chengdu: Sichuan Science and Technology PressGoogle Scholar
- Goldberg RB, Beals TP, Sanders PM: Anther development: basic principles and practical applications. Plant Cell. 1993, 5: 1217-1229.PubMed CentralPubMedGoogle Scholar
- Pacini E, Franchi GG, Hesse M: The tapetum: its form, function and possible phylogeny in Embryophyta. Plant Syst Evol. 1985, 149: 155-185. 10.1007/BF00983304.Google Scholar
- Kaul ML: Male sterility in higher plants: Monograph on Theoretical and Applied Genetics. 1988, Berlin:Google Scholar
- Scott RJ, Spielman M, Dickinson HG: Stamen structure and function. Plant Cell. 2004, 16: S46-S60. 10.1105/tpc.017012.PubMed CentralPubMedGoogle Scholar
- Ma H: Molecular genetic analysis of microsporogenesis and microgametogenesis in flowering plants. Annu Rev Plant Biol. 2005, 56: 393-434. 10.1146/annurev.arplant.55.031903.141717.PubMedGoogle Scholar
- Bhatt AM, Canales C, Dickinson HG: Plant meiosis: the means to 1N. Trends Plant Sci. 2001, 6: 114-121. 10.1016/S1360-1385(00)01861-6.PubMedGoogle Scholar
- Schiefthaler U, Balasubramanian S, Sieber P, Chevalier D, Wisman E, Schneitz K: Molecular analysis of NOZZLE, a gene involved in pattern formation and early sporogenesis during sex organ development in Arabidopsis thaliana. Proc Natl Acad Sci USA. 1999, 96: 1664-1669.Google Scholar
- Yang WC, Ye D, Xu J, Sundaresan V: The SPOROCYTELESS gene of Arabidopsis is required for initiation of sporogenesis and encodes a novel nuclear protein. Genes Dev. 1999, 13: 2108-2117. 10.1101/gad.13.16.2108.PubMed CentralPubMedGoogle Scholar
- Sorensen AM, Krober S, Unte US, Huijser P, Dekker K, Saedler H: The Arabidopsis ABORTED MICROSPORES (AMS) gene encodes a MYC class transcription factor. Plant J. 2003, 33: 413-423. 10.1046/j.1365-313X.2003.01644.x.PubMedGoogle Scholar
- Yang SL, Xie L, Mao HZ, Puah CS, Yang WC, Jiang L, Sundaresan V, Ye D: TAPETUM DETERMINANT1 is MS1 Is Required for Tapetal Development 3547 required for cell specialization in the Arabidopsis anther. Plant Cell. 2003, 15: 2792-2804. 10.1105/tpc.016618.PubMed CentralPubMedGoogle Scholar
- Zhang W, Sun YL, Timofejeva L, Chen C, Grossniklaus U, Ma H: Regulation of Arabidopsis tapetum development and function by DYSFUNCTIONAL TAPETUM (DYT1) encoding a putative bHLH transcription factor. Developmen t. 2006, 133: 3085-3095.Google Scholar
- Jung KH, Han MJ, Lee YS, Kim YW, Hwang I, Kim MJ, Kim YK, Nahm BH, An G: Rice Undeveloped Tapetum1 is a major regulator of early tapetum development. Plant Cell. 2005, 17: 2705-2722. 10.1105/tpc.105.034090.PubMed CentralPubMedGoogle Scholar
- Yang CY, Vizcay-Barrena G, Conner K, Wilson ZA: MALE STERILITY1 is required for tapetal development and pollen wall biosynthesis. Plant Cell. 2007, 19: 3530-3548. 10.1105/tpc.107.054981.PubMed CentralPubMedGoogle Scholar
- Aarts MG, Hodge R, Kalantidis K, Florack D, Wilson ZA, Mulligan BJ, Stiekema WJ, Scott R, Pereira A: The Arabidopsis MALE STERILITY 2 protein shares similarity with reductases in elongation/condensation complexes. Plant J. 1997, 12: 615-623. 10.1046/j.1365-313X.1997.d01-8.x.PubMedGoogle Scholar
- Anna AD, Jay S, Marc M, Franck P, Michiyo M, Robert S, Birger LM, Daphne P: CYP704B1 is a long-chain fatty acid v-Hydroxylase essential for sporopollenin synthesis in pollen of arabidopsis. Plant Phy. 2009, 151: 574-589. 10.1104/pp.109.144469.Google Scholar
- Li H, Pinot F, Sauveplane V, Werck-Reichhart D, Diehl P, Schreiber L, Franke R, Zhang P, Chen L, Gao YW, Liang WQ, Zhang DB: Cytochrome P450 family member CYP704B2 catalyzes the v-hydroxylation of fatty acids and is required for anther Cutin biosynthesis and pollen exine formation in rice. Plant Cell. 2010, 22: 173-190. 10.1105/tpc.109.070326.PubMed CentralPubMedGoogle Scholar
- Dobritsa AA, Lei ZT, Nishikawa S, Urbanczyk-Wochniak E, Huhman DV, Preuss D, Sumner LW: LAP5 and LAP6 encode anther-specific proteins with similarity to chalcone synthase essential for pollen exine Development in Arabidopsis. Plant Phy. 2010, 153: 937-955. 10.1104/pp.110.157446.Google Scholar
- Zhang JF, Turley RBMD, Stewart J: Comparative analysis of gene expression between CMS-D8 restored plants and normal non-restoring fertile plants in cotton by differential display. Plant Cell Rep. 2008, 27: 553-561. 10.1007/s00299-007-0492-7.PubMedGoogle Scholar
- He JP, Ke LP, Hong DF, Xie YZ, Wang GC, Liu PW, Yang GS: Fine mapping of a recessive genic male sterility gene (Bnms3) in rapeseed (Brassica napus) with AFLP and Arabidopsis-derived PCR markers. Theor Appl Genet. 2008, 117: 11-18. 10.1007/s00122-008-0747-x.PubMedGoogle Scholar
- Bentley DR: Whole-genome re-sequencing. Curr Opin Genet Dev. 2006, 16: 545-552. 10.1016/j.gde.2006.10.009.PubMedGoogle Scholar
- Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, Menezes RX, Boer JM, Ommen GJ, Dunnen JT: Deep sequencing- based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36: 141-145. 10.1093/nar/gkn705.Google Scholar
- Lewis Z, Li HJ, Schmidt-Küntzel A, Warren WC, Gregory SB: Digital gene expression for non-model organisms. Genome Res. 2011, 21: 1905-1915. 10.1101/gr.122135.111.Google Scholar
- Hang GW, Zhang DM, Huang YX, Guo YG, Shi SC, Li ZG, Gao DK, Feng GB: The utilization of male sterile recessive genes in hybrid seed production of cotton (G. hirsutum ). Scientia Agricultura Sinica. 1981, 1: 5–1-Google Scholar
- Honys D, Twell D: Transcriptome analysis of haploid male gametophyte development in Arabidopsis. Genome Bio. 2004, 5: R85-10.1186/gb-2004-5-11-r85.Google Scholar
- Wang DX, Oses-Prieto JA, Li KH, Fernandes JF, Burlingame AL, Walbot V: The male sterile 8 mutation of maize disrupts the temporal progression of the transcriptome and results in the mis-regulation of metabolic functions. Plant J. 2010, 63: 939-951. 10.1111/j.1365-313X.2010.04294.x.PubMed CentralPubMedGoogle Scholar
- Ma J, Skibbe DS, Fernandes J, Walbot V: Male reproductive development: gene expression profiling of maize anther and pollen ontogeny. Genome Bio. 2008, 9: R181-10.1186/gb-2008-9-12-r181.Google Scholar
- Xie FL, Sun GL, Stiller JW, Zhang BH: Genome-wide functional analysis of the cotton transcriptome by creating an integrated EST database. PLoS One. 2011, 6: 11-Google Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 10: 986-995.Google Scholar
- Coll NS, Epple P, Dang JL: Programmed cell death in the plant immune system. Cell Death Differ. 2011, 11: 1-10.Google Scholar
- Goetz M, Godt DE, Guivarch A, Kahmann U, Chriqui D, Roitsch T: Induction of male sterility in plants by metabolic engineering of the carbohydrate supply. Proc Natl Acad Sci USA. 2001, 98: 6522-6527. 10.1073/pnas.091097998.PubMed CentralPubMedGoogle Scholar
- Oliver SN, Dennis ES, Dolferus R: ABA regulates apoplastic sugar transport and is a potential signal for cold-induced pollen sterility in rice. Plant Cell Physiol. 2007, 48: 1319-1330. 10.1093/pcp/pcm100.PubMedGoogle Scholar
- Datta R, Chamusco KC, Chourey PS: Starch biosynthesis during pollen maturation is associated with altered patterns of gene expression in maize. Plant Physiol. 2002, 130: 1645-1656. 10.1104/pp.006908.PubMed CentralPubMedGoogle Scholar
- Mamun EA, Alfred S, Cantrill LC, Overall RL, Sutton BG: Effects of chilling on male gametophyte development in rice. Cell Biol. 2006, 30: 583-591.Google Scholar
- Ahlers F, Lambert J, Wiermann R: Acetylation and silylation of piperidine solubilized sporopollenin from pollen of Typha angustifolia L. Z Naturforsch. 2003, 58: 807-811.Google Scholar
- Wiermann R, Ahlers F, Schmitz-Thom I: Sporopollenin: 209–227. 2001, Weinheim, Germany: Wiley-VCH VerlagGoogle Scholar
- Piffanelli P, Ross JH, Murphy DJ: Biogenesis and function of the lipidic structures of pollen grains. Sex Plant Reprod. 1998, 11: 65-80. 10.1007/s004970050122.Google Scholar
- Johnson PR, Ecker JR: The ethylene gas signal transduction pathway: a molecular perspective. Annu Rev Genet. 1998, 32: 227-254. 10.1146/annurev.genet.32.1.227.PubMedGoogle Scholar
- Bleecker AB, Kende H: Ethylene: a gaseous signal molecule in plants. Annu Rev Cell Dev Biol. 2000, 16: 1-18. 10.1146/annurev.cellbio.16.1.1.PubMedGoogle Scholar
- Kovaleva LV, Alla D, Alexander V, Viktor R: Ethylene is involved in the control of male gametophyte development and germination in Petunia. J Plant Growth Regul. 2011, 30: 64-73. 10.1007/s00344-010-9168-6.Google Scholar
- Bagnall DJ: Control of flowering in Arabidopsis thaliana by light, vernalisation and gibberellins. Aust J Plant Physiol. 1992, 19: 401-409. 10.1071/PP9920401.Google Scholar
- Pharis RP, King RW: Gibberellins and reproductive development in seed plants. Annu Rev Plant Physiol. 1985, 36: 517-568. 10.1146/annurev.pp.36.060185.002505.Google Scholar
- Thornsberry JM, Goodman MM, Doebley J, Kresovich S, Nielsen D, Buckler ES: Dwarf8 polymorphisms associate with variation in flowering time. Nat Genet. 2001, 28: 286-289. 10.1038/90135.PubMedGoogle Scholar
- Cheng X, Blumenthal RM: Mammalian DNA methyltransferases: a structural perspective. Structure. 2008, 16: 341-350. 10.1016/j.str.2008.01.004.PubMed CentralPubMedGoogle Scholar
- Berger SL: The complex language of chromatin regulation during transcription. Nature. 2007, 447: 407-412. 10.1038/nature05915.PubMedGoogle Scholar
- Ishimaru K, Takada K, Watanabe S, Kamada H, Ezura H: Stable male sterility induced by the expression of mutated melon ethylene receptor genes in Nicotiana tabacum. Plant Sci. 2006, 3: 355-359.Google Scholar
- He YC, Wei Q, Ge J, Jiang AM, Gan L, Song ZJ, Cai DT: Genome duplication effects on pollen development and the interrelated physiological substances in tetraploid rice with polyploid meiosis stability. Planta. 2010, 232: 1219-1228. 10.1007/s00425-010-1249-z.PubMedGoogle Scholar
- Feng S: Coordinated regulation of Arabidopsis thaliana development by light and gibberellins. Nature. 2008, 451: 475-479. 10.1038/nature06448.PubMed CentralPubMedGoogle Scholar
- Ueguchi-Tanaka M, Hirano K, Hasegawa Y, Kitano H, Matsuoka M: Release of the repressive activity of rice DELLA protein SLR1by Gibberellin does not require SLR1 degradation in the gid2 mutant. Plant Cell. 2004, 16: 2001-2019. 10.1105/tpc.104.022830.Google Scholar
- Liu XY, Yang XD, Zhao X, Yang JH, Zhang MF: Reduced expression of CTR1 gene modulated by mitochondria causes enhanced ethylene response in cytoplasmic male-sterile Brassica juncea. PhysiologiaPlantarum. 2012, 145: 332-340.Google Scholar
- Meer IM, Stam ME, Tunen AJ, Mol JN, Stuitje AR: Antisense inhibition of flavonoid biosynthesis in petunia anthers results in male sterility. Plant Cell. 1992, 4: 253-262.PubMed CentralPubMedGoogle Scholar
- Stefan M, Anja P, Ulrich M: Multifunctional flavonoid dioxygenases: flavonol and anthocyanin biosynthesis in Arabidopsis thaliana L. Phytochemistry. 2010, 71: 1040-1049. 10.1016/j.phytochem.2010.04.016.Google Scholar
- Mo Y, Nagel C, Taylor LP: Biochemical complementation of chalcone synthase mutants defines a role for flavonoids in functional pollen. Proc Natl Acad Sci USA. 1992, 89: 7213-7217. 10.1073/pnas.89.15.7213.PubMed CentralPubMedGoogle Scholar
- Christopher J, Staiger NS, Poulter JL, Henty VE, Franklin T: Regulation of actin dynamics by actin-binding proteins in pollen. J Exp Bot. 2010, 16: 1-18.Google Scholar
- Xu CG, Liu ZT, Zhang LP, Zhao CP, Yuan SH, Zhang FT: Organization of actin cytoskeleton during meiosis I in a wheat thermo-sensitive genic male sterile line. Protoplasma. 2012, 10: 1-8.Google Scholar
- Li XB, Xu D, Wang XL, Huang GQ, Luo J, Li DD, Zhang ZT, Xu WL: Three cotton genes preferentially expressed in flower tissues encode actin-depolymerizing factors which are involved in F-actin dynamics in cells. J Exp Bot. 2010, 61: 41-53. 10.1093/jxb/erp280.PubMed CentralPubMedGoogle Scholar
- Frank VB, James FD: Reactive oxygen species in plant cell death. Plant Physiol. 2006, 141: 384-390. 10.1104/pp.106.078295.Google Scholar
- Jiang PD, Zhang XQ, Zhu YG, Zhu W, Xie HY, Wang XD: Metabolism of reactive oxygen species in cotton cytoplasmic male sterility and its restoration. Plant Cell Rep. 2007, 26: 1627-1634. 10.1007/s00299-007-0351-6.PubMedGoogle Scholar
- Lorrain S, Vailleau F, Balague C, Roby D: Lesion mimic mutants: keys for deciphering cell death and defense pathways in plants. Trends Plant Sci. 2003, 8: 263-271. 10.1016/S1360-1385(03)00108-0.PubMedGoogle Scholar
- Hou L, Xiao YH, Li XB, Wang WF, Luo XY, Pei Y: The cDNA-AFLP differential display in developing anthers between cotton male sterile and fertile line of “Dong A”. Acta Genetica Sinica. 2002, 29: 359-363.PubMedGoogle Scholar
- Fu ZD, Zhang ZL, Qu WQ: Metabolism: Experiments of plant physiology. 2004, Higher Education Harbor: Higher Education PressGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.