Control of plant trichome development by a cotton fiber MYB gene

Control of trichome formation in Arabidopsis by poplar single-repeat R3 MYB transcription factors

In Arabidopsis, trichome formation is regulated by the interplay of R3 MYBs and several others transcription factors including the WD40-repeat protein TRANSPARENT TESTA GLABRA1 (TTG1), the R2R3 MYB transcription factor GLABRA1 (GL1), the bHLH transcription factor GLABRA3 (GL3) or ENHANCER OF GLABRA3 (EGL3), and the homeodomain protein GLABRA2 (GL2). R3 MYBs including TRICHOMELESS1 (TCL1), TCL2, TRYPTICHON (TRY), CAPRICE (CPC), ENHANCER OF TRY AND CPC1 (ETC1), ETC2 and ETC3 negatively regulate trichome formation by competing with GL1 for binding GL3 or EGL3, thus blocking the formation of TTG1-GL3/EGL3-GL1, an activator complex required for the activation of the trichome positive regulator gene GL2. However, it is largely unknown if R3 MYBs in other plant species especially woody plants have similar functions. By BLASTing the Populus trichocarpa protein database using the entire amino acid sequence of TCL1, an Arabidopsis R3 MYB transcription factor, we identified a total of eight R3 MYB transcription factor genes in poplar, namely P. trichocarpa TRICHOMELESS1 through 8 (PtrTCL1-PtrTCL8). The amino acid signature required for interacting with bHLH transcription factors and the amino acids required for cell-to-cell movement of R3 MYBs are not fully conserved in all PtrTCLs. When tested in Arabidopsis protoplasts, however, all PtrTCLs interacted with GL3. Expressing each of the eight PtrTCL genes in Arabidopsis resulted in either glabrous phenotypes or plants with reduced trichome numbers, and expression levels of GL2 in all transgenic plants tested were greatly reduced. Expression of PtrTCL1 under the control of TCL1 native promoter almost completely complemented the mutant phenotype of tcl. In contrast, expression of PtrTCL1 under the control of TRY native promoter in the try mutant, or under the control of CPC native promoter in the cpc mutant resulted in glabrous phenotypes, suggesting that PtrTCL1 functions similarly to TCL1, but not TRY and CPC.

Functional analysis of the seed coat-specific gene GbMYB2 from cotton

MYB transcription factors are essential for cotton fiber development. We isolated the R2R3-MYB gene GbMYB2 from Gossypium barbadense. RNA in situ hybridization analysis showed that GbMYB2 is mainly expressed in the outer integuments of cotton ovules and in elongating fibers. GbMYB2 expression increased throughout fiber initiation and during the elongation stage. The expression level of GbMYB2 was higher in the Gossypium hirsutum cultivar Xu142 than in the Xu142 (fl) mutant. Overexpression of GbMYB2 in Arabidopsis caused thicker leaf trichomes and longer roots to develop due to the activation of trichome development-related genes such as GL2. These results indicate that GbMYB2 is an R2R3-MYB gene that is involved in fiber development.

Induction of cotton ovule culture fibre branching by co-expression of cotton BTL, cotton SIM, and Arabidopsis STI genes

The highly elongated single-celled cotton fibre consists of lint and fuzz, similar to the Arabidopsis trichome. Endoreduplication is an important determinant in Arabidopsis trichome initiation and morphogenesis. Fibre development is also controlled by functional homologues of Arabidopsis trichome patterning genes, although fibre cells do not have a branched shape like trichomes. The identification and characterization of the homologues of 10 key Arabidopsis trichome branching genes in Gossypium arboreum are reported here. Nuclear ploidy of fibres was determined, and gene function in cotton callus and fibre cells was investigated. The results revealed that the nuclear DNA content was constant in fuzz, whereas a limited and reversible change occurred in lint after initiation. Gossypeum arboreum branchless trichomes (GaBLT) was not transcribed in fibres. The homologue of STICHEL (STI), which is essential for trichome branching, was a pseudogene in Gossypium. Targeted expression of GaBLT, Arabidopsis STI, and the cytokinesis-repressing GaSIAMESE in G. hirsutum fibre cells cultured in vitro resulted in branching. The findings suggest that the distinctive developmental mechanism of cotton fibres does not depend on endoreduplication. This important component may be a relic function that can be activated in fibre cells.

Medicago glucosyltransferase UGT72L1: potential roles in proanthocyanidin biosynthesis

In the first reaction specific for proanthocyanidin (PA) biosynthesis in Arabidopsis thaliana and Medicago truncatula, anthocyanidin reductase (ANR) converts cyanidin to (-)-epicatechin. The glucosyltransferase UGT72L1 catalyzes formation of epicatechin 3'-O-glucoside (E3'OG), the preferred substrate for MATE transporters implicated in PA biosynthesis in both species. The mechanism of PA polymerization is still unclear, but may involve the laccase-like polyphenol oxidase TRANSPARENT TESTA 10 (TT10). We have employed a combination of cell biological, biochemical and genetic approaches to evaluate this PA pathway model. The promoter regions of UGT72L1 and MtANR share common cis-acting elements and direct overlapping, but partially distinct, expression patterns. UGT72L1 and MtANR are localized in the cytosol, whereas TT10 is localized to the vacuole. Over-expression of UGT72L1 in M. truncatula hairy roots results in increased accumulation of PA-like compounds, and loss of function of UGT72L1 partially reduces epicatechin, E3'OG and extractable PA levels in M. truncatula seeds. Expression of UGT72L1 in A. thaliana leads to a massive increase in E3'OG in immature seed, but reduced levels of extractable PAs. However, when UGT72L1 was expressed in the Arabidopsis tt10 mutant, extractable PA levels increased and seed coat browning was delayed. Our results suggest that glycosylation of epicatechin is important for both PA precursor transport and assembly, but that additional redundant pathways may exist.

The cotton transcription factor TCP14 functions in auxin-mediated epidermal cell differentiation and elongation

Plant-specific TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors play crucial roles in development, but their functional mechanisms remain largely unknown. Here, we characterized the cellular functions of the class I TCP transcription factor GhTCP14 from upland cotton (Gossypium hirsutum). GhTCP14 is expressed predominantly in fiber cells, especially at the initiation and elongation stages of development, and its expression increased in response to exogenous auxin. Induced heterologous overexpression of GhTCP14 in Arabidopsis (Arabidopsis thaliana) enhanced initiation and elongation of trichomes and root hairs. In addition, root gravitropism was severely affected, similar to mutant of the auxin efflux carrier PIN-FORMED2 (PIN2) gene. Examination of auxin distribution in GhTCP14-expressing Arabidopsis by observation of auxin-responsive reporters revealed substantial alterations in auxin distribution in sepal trichomes and root cortical regions. Consistent with these changes, expression of the auxin uptake carrier AUXIN1 (AUX1) was up-regulated and PIN2 expression was down-regulated in the GhTCP14-expressing plants. The association of GhTCP14 with auxin responses was also evidenced by the enhanced expression of auxin response gene IAA3, a gene in the AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) family. Electrophoretic mobility shift assays showed that GhTCP14 bound the promoters of PIN2, IAA3, and AUX1, and transactivation assays indicated that GhTCP14 had transcription activation activity. Taken together, these results demonstrate that GhTCP14 is a dual-function transcription factor able to positively or negatively regulate expression of auxin response and transporter genes, thus potentially acting as a crucial regulator in auxin-mediated differentiation and elongation of cotton fiber cells.

Identification and characterization of cotton genes involved in fuzz-fiber development

Cotton fuzz fibers, like Arabidopsis trichomes, are elongated unicells. It is postulated that a transcriptional complex of GLABRA1 (GL1), GLABRA3 (GL3), and TRANSPARENT TESTAGLABRA1 (TTG1) might be in existence in Arabidopsis as evidenced by their physical interaction in yeast, and the complex regulates expression of GLABRA2 (GL2) controlling trichome cell differentiation; it is also assumed that TRIPTYCHON (TRY) and CAPRICE (CPC) counteract the complex formation in neighboring cells. Here, the homologs GaMYB23 (a homolog of GL1), GaDEL65 (a homolog of GL3), GaTTG1, GaCPC and GaTRY were identified in Gossypium arboreum. We show that GaMYB23 can bind to and activate the promoters of GaCPC, GaGL2 and GaTRY, and that GaMYB23, GaTRY and GaTTG1 could interact with GaDEL65 in yeast and in planta. In situ analysis showed that GaMYB23, GaGL2, GaDEL65, and GaTRY were predominantly expressed in fuzz fiber, but GaTRY proteins were primarily found in undeveloped epidermal cells. A G. arboreum fuzzless mutant with consistently high level GaMYB23 transcript has lost the detectable GaMYB23-promoter of GaGL2 complex, corresponding to sharply reduced transcription of GaGL2. Our results support that cotton homologs to the genetic molecules regulating Arabidopsis trichome differentiation interacted in the epidermis of ovules and the redundant GaMYB23 serves as a negative regulator in fuzz-fiber patterning.

Using genome-referenced expressed sequence tag assembly to analyze the origin and expression patterns of Gossypium hirsutum transcripts

We assembled a total of 297,239 Gossypium hirsutum (Gh, a tetraploid cotton, AADD) expressed sequence tag (EST) sequences that were available in the National Center for Biotechnology Information database, with reference to the recently published G. raimondii (Gr, a diploid cotton, DD) genome, and obtained 49,125 UniGenes. The average lengths of the UniGenes were increased from 804 and 791 bp in two previous EST assemblies to 1,019 bp in the current analysis. The number of putative cotton UniGenes with lengths of 3 kb or more increased from 25 or 34 to 1,223. As a result, thousands of originally independent G. hirsutum ESTs were aligned to produce large contigs encoding transcripts with very long open reading frames, indicating that the G. raimondii genome sequence provided remarkable advantages to assemble the tetraploid cotton transcriptome. Significant different distribution patterns within several GO terms, including transcription factor activity, were observed between D- and A-derived assemblies. Transcriptome analysis showed that, in a tetraploid cotton cell, 29,547 UniGenes were possibly derived from the D subgenome while another 19,578 may come from the A subgenome. Finally, some of the in silico data were confirmed by reverse transcription polymerase chain reaction experiments to show the changes in transcript levels for several gene families known to play key role in cotton fiber development. We believe that our work provides a useful platform for functional and evolutionary genomic studies in cotton.

Trichomes as models for studying plant cell differentiation

Trichomes, originating from epidermal cells, are present on nearly all terrestrial plants. They exist in diverse forms, are readily accessible, and serve as an excellent model system for analyzing the molecular mechanisms in plant cell differentiation, including cell fate choices, cell cycle control, and cell morphogenesis. In Arabidopsis, two regulatory models have been identified that function in parallel in trichome formation; the activator-inhibitor model and the activator-depletion model. Cotton fiber, a similar unicellular structure, is controlled by some functional homologues of Arabidopsis trichome-patterning genes. Multicellular trichomes, as in tobacco and tomato, may form through a distinct pathway from unicellular trichomes. Recent research has shown that cell cycle control participates in trichome formation. In this review, we summarize the molecular mechanisms involved in the formation of unicellular and multicellular trichomes, and discuss the integration of the cell cycle in its initiation and morphogenesis.

A cotton BURP domain protein interacts with α-expansin and their co-expression promotes plant growth and fruit production

Plant growth requires cell wall extension. The cotton AtRD22-Like 1 gene GhRDL1, predominately expressed in elongating fiber cells, encodes a BURP domain-containing protein. Here, we show that GhRDL1 is localized in cell wall and interacts with GhEXPA1, an α-expansin functioning in wall loosening. Transgenic cotton overexpressing GhRDL1 showed an increase in fiber length and seed mass, and an enlargement of endopleura cells of ovules. Expression of either GhRDL1 or GhEXPA1 alone in Arabidopsis led to a substantial increase in seed size; interestingly, their co-expression resulted in the increased number of siliques, the nearly doubled seed mass, and the enhanced biomass production. Cotton plants overexpressing GhRDL1 and GhEXPA1 proteins produced strikingly more fruits (bolls), leading to up to 40% higher fiber yield per plant without adverse effects on fiber quality and vegetative growth. We demonstrate that engineering cell wall protein partners has a great potential in promoting plant growth and crop yield.

LuFLA1PRO and LuBGAL1PRO promote gene expression in the phloem fibres of flax (Linum usitatissimum)

Cell type-specific promoters were identified that drive gene expression in an industrially important product. To identify flax (Linum usitatissimum) gene promoters, we analyzed the genomic regions upstream of a fasciclin-like arabinogalactan protein (LuFLA1) and a beta-galactosidase (LuBGAL1). Both of these genes encode transcripts that have been found to be highly enriched in tissues bearing phloem fibres. Using a beta-glucuronidase (GUS) reporter construct, we found that a 908-bp genomic sequence upstream of LuFLA1 (LuFLA1PRO) directed GUS expression with high specificity to phloem fibres undergoing secondary cell wall development. The DNA sequence upstream of LuBGAL1 (LuBGAL1PRO) likewise produced GUS staining in phloem fibres with developing secondary walls, as well as in tissues of developing flowers and seed bolls. These data provide further evidence of a specific role for LuFLA1 in phloem fibre development, and demonstrate the utility of LuFLA1PRO and LuBGAL1PRO as tools for biotechnology and further investigations of phloem fibre development.

Construction of genetic linkage map with chromosomal assigment and quantitative trait loci associated with some important agronomic traits in cotton

Cotton is the world's leading natural fiber and second most important oilseed crop and has been a focus of genetic, systematic and breeding research. The genetic and physiological bases of some important agronomic traits in cotton were investigated by QTL mapping through constructing of genetic map with chromosomal assignment. A segregating F2 population derived from an interspecific cross (G. barbadense x G. hirsutum) between two genotypes, cvs. "Giza 83" and "Deltapine" was used in this study. Different molecular markers including SSR, EST, EST-SSR, AFLP and RAPD were employed to identify markers that reveal differences between the parents. In total 42 new markers were merged with 140 previously mapped markers to produce a new map with 182 loci covering a total length of 2370.5 cM. Among these new markers, some of them were used to assign chromosomes to the produced 26 linkage groups. The LG2, LG3, LG11 and LG26 were assigned to chromosomes 1, 6, 5 and 20 respectively. Single point analysis was used to identify genomic regions controlling traits for plant height, number of nodes at flowering time, bolling date, days to flowering and number of bolls. In total 40 significant QTL were identified for the five traits on 11 linkage groups (1, 2, 3, 4, 5, 10, 11, 12, 18, 19 and 23). This work represents an improvement of the previously constructed genetic map in addition to chromosomal assignment and detection of new significant QTL for the five traits in Egyptian cotton. The Significant QTLs detected in this study can be employed in marker assisted selection for molecular breeding programs aiming at developing cotton cultivars with improved agronomic traits.

Agrobacterium rhizogenes-induced cotton hairy root culture as an alternative tool for cotton functional genomics

Although well-accepted as the ultimate method for cotton functional genomics, Agrobacterium tumefaciens-mediated cotton transformation is not widely used for functional analyses of cotton genes and their promoters since regeneration of cotton in tissue culture is lengthy and labor intensive. In certain cases, A. rhizogenes-induced hairy root culture has been a suitable molecular tool for functional analyses of genes and promoters for plants that are difficult to regenerate by A. tumefaciens-mediated transformation. Similarly, A. rhizogenes-induced hairy root cultures are an alternative tool for cotton functional genomics. In this chapter, the advantages and disadvantages of using A. rhizogenes-induced cotton hairy root culture over A. tumefaciens-mediated cotton transformation are discussed. The procedures for transformation, generation, selection, and molecular analyses of transgenic cotton hairy roots are introduced by describing the functional analysis of a cotton promoter in cotton hairy roots generated by A. rhizogenes-mediated transformation.

Plant glandular trichomes as targets for breeding or engineering of resistance to herbivores

Glandular trichomes are specialized hairs found on the surface of about 30% of all vascular plants and are responsible for a significant portion of a plant's secondary chemistry. Glandular trichomes are an important source of essential oils, i.e., natural fragrances or products that can be used by the pharmaceutical industry, although many of these substances have evolved to provide the plant with protection against herbivores and pathogens. The storage compartment of glandular trichomes usually is located on the tip of the hair and is part of the glandular cell, or cells, which are metabolically active. Trichomes and their exudates can be harvested relatively easily, and this has permitted a detailed study of their metabolites, as well as the genes and proteins responsible for them. This knowledge now assists classical breeding programs, as well as targeted genetic engineering, aimed to optimize trichome density and physiology to facilitate customization of essential oil production or to tune biocide activity to enhance crop protection. We will provide an overview of the metabolic diversity found within plant glandular trichomes, with the emphasis on those of the Solanaceae, and of the tools available to manipulate their activities for enhancing the plant's resistance to pests.

A WUSCHEL-related homeobox 3B gene, depilous (dep), confers glabrousness of rice leaves and glumes

BACKGROUND

Glabrousness is an important agricultural trait for the practical breeding of rice. In this study, depilous (dep), the gene responsible for glabrous leaves and glumes of rice was identified by map-based cloning.

RESULTS

The dep gene encodes a WUSCHEL-related homeobox 3B that was fine-mapped to a 22-kb region on the short arm of chromosome 5 using progenies derived from crosses between Koshihikari (pubescent) and GLSL15, an Oryza glaberrima chromosome segment substitution line (glabrous). Complementation tests confirmed the conditioning of the glabrous phenotype by the dep gene. Phylogenetic analysis showed that dep groups with the WOX3 family of plant-specific homeobox transcription factors that are involved in regulating lateral organ development. Localization of dep in the nucleus indicates the function of the gene as a transcription factor. Spatial expression of the gene was observed in the base of young shoots, the leaf sheath, midrib, young roots and nodal structures.

CONCLUSION

The identification and cloning of dep will not only provide basis for future research on the elucidation of the molecular mechanisms underlying trichome formation in rice but will also aid in breeding programs for the development of glabrous varieties.

Comparative transcriptomic analysis of roots of contrasting Gossypium herbaceum genotypes revealing adaptation to drought

Background

Root length and its architecture govern the adaptability of plants to various stress conditions, including drought stress. Genetic variations in root growth, length, and architecture are genotypes dependent. In this study, we compared the drought-induced transcriptome of four genotypes of Gossypium herbaceum that differed in their drought tolerance adaptability. Three different methodologies, namely, microarray, pyrosequencing, and qRT–PCR, were used for transcriptome analysis and validation.

Results

The variations in root length and growth were found among four genotypes of G.herbaceum when exposed to mannitol-induced osmotic stress. Under osmotic stress, the drought tolerant genotypes Vagad and GujCot-21 showed a longer root length than did by drought sensitive RAHS-14 and RAHS-IPS-187. Further, the gene expression patterns in the root tissue of all genotypes were analyzed. We obtained a total of 794 differentially expressed genes by microarray and 104928 high-quality reads representing 53195 unigenes from the root transcriptome. The Vagad and GujCot-21 respond to water stress by inducing various genes and pathways such as response to stresses, response to water deprivation, and flavonoid pathways. Some key regulatory genes involved in abiotic stress such as AP2 EREBP, MYB, WRKY, ERF, ERD9, and LEA were highly expressed in Vagad and GujCot-21. The genes RHD3, NAP1, LBD, and transcription factor WRKY75, known for root development under various stress conditions, were expressed specifically in Vagad and GujCot-21. The genes related to peroxidases, transporters, cell wall-modifying enzymes, and compatible solutes (amino acids, amino sugars, betaine, sugars, or sugar alcohols) were also highly expressed in Vagad and Gujcot-21.

Conclusion

Our analysis highlights changes in the expression pattern of genes and depicts a small but highly specific set of drought responsive genes induced in response to drought stress. Some of these genes were very likely to be involved in drought stress signaling and adaptation, such as transmembrane nitrate transporter, alcohol dehydrogenase, pyruvate decarboxylase, sucrose synthase, and LEA. These results might serve as the basis for an in-depth genomics study of Gossypium herbaceum, including a comparative transcriptome analysis and the selection of genes for root traits and drought tolerance.

Deep sequencing reveals differences in the transcriptional landscapes of fibers from two cultivated species of cotton

Cotton (Gossypium) fiber is the most prevalent natural product used in the textile industry. The two major cultivated species, G. hirsutum (Gh) and G. barbadense (Gb), are allotetraploids with contrasting fiber quality properties. To better understand the molecular basis for their fiber differences, EST pyrosequencing was used to document the fiber transcriptomes at two key development stages, 10 days post anthesis (dpa), representing the peak of fiber elongation, and 22 dpa, representing the transition to secondary cell wall synthesis. The 617,000 high quality reads (89% of the total 692,000 reads) from 4 libraries were assembled into 46,072 unigenes, comprising 38,297 contigs and 7,775 singletons. Functional annotation of the unigenes together with comparative digital gene expression (DGE) revealed a diverse set of functions and processes that were partly linked to specific fiber stages. Globally, 2,770 contigs (7%) showed differential expression (>2-fold) between 10 and 22 dpa (irrespective of genotype), with 70% more highly expressed at 10 dpa, while 2,248 (6%) were differentially expressed between the genotypes (irrespective of stage). The most significant genes with differential DGE at 10 dpa included expansins and lipid transfer proteins (higher in Gb), while at 22 dpa tubulins, cellulose, and sucrose synthases showed higher expression in Gb. DGE was compared with expression data of 10 dpa-old fibers from Affymetrix microarrays. Among 543 contigs showing differential expression on both platforms, 74% were consistent in being either over-expressed in Gh (242 genes) or in Gb (161 genes). Furthermore, the unigene set served to identify 339 new SSRs and close to 21,000 inter-genotypic SNPs. Subsets of 88 SSRs and 48 SNPs were validated through mapping and added 65 new loci to a RIL genetic map. The new set of fiber ESTs and the gene-based markers complement existing available resources useful in basic and applied research for crop improvement in cotton.

Composition and physiological function of the wax layers coating Arabidopsis leaves: β-amyrin negatively affects the intracuticular water barrier

Plants prevent dehydration by coating their aerial, primary organs with waxes. Wax compositions frequently differ between species, organs, and developmental stages, probably to balance limiting nonstomatal water loss with various other ecophysiological roles of surface waxes. To establish structure-function relationships, we quantified the composition and transpiration barrier properties of the gl1 mutant leaf waxes of Arabidopsis (Arabidopsis thaliana) to the necessary spatial resolution. The waxes coating the upper and lower leaf surfaces had distinct compositions. Moreover, within the adaxial wax, the epicuticular layer contained more wax and a higher relative quantity of alkanes, whereas the intracuticular wax had a higher percentage of alcohols. The wax formed a barrier against nonstomatal water loss, where the outer layer contributed twice as much resistance as the inner layer. Based on this detailed description of Arabidopsis leaf waxes, structure-function relationships can now be established by manipulating one cuticle component and assessing the effect on cuticle functions. Next, we ectopically expressed the triterpenoid synthase gene AtLUP4 (for lupeol synthase4 or β-amyrin synthase) to compare water loss with and without added cuticular triterpenoids in Arabidopsis leaf waxes. β-Amyrin accumulated solely in the intracuticular wax, constituting up to 4% of this wax layer, without other concomitant changes of wax composition. This triterpenoid accumulation caused a significant reduction in the water barrier effectiveness of the intracuticular wax.

GbTCP, a cotton TCP transcription factor, confers fibre elongation and root hair development by a complex regulating system

As the most important natural raw material for textile industry, cotton fibres are an excellent model for studying single-cell development. Although expression profiling and functional genomics have provided some data, the mechanism of fibre development is still not well known. A class I TCP transcription factor (designated GbTCP), encoding 344 amino acids, was isolated from the normalized cDNA library of sea-island cotton fibre (from -2 to 25 days post anthesis). GbTCP was preferentially expressed in the elongating cotton fibre from 5 to 15 days post anthesis. Some expression was also observed in stems, apical buds, and petals. RNAi silencing of GbTCP produced shorter fibre, a reduced lint percentage, and a lower fibre quality than the wild-type plants. Overexpression of GbTCP enhanced root hair initiation and elongation in Arabidopsis and regulated branching. Solexa sequencing and Affymetrix GeneChip analysis indicated that GbTCP positively regulates the level of jasmonic acid (JA) and, as a result, activates downstream genes (reactive oxygen species, calcium signalling, ethylene biosynthesis and response, and several NAC and WRKY transcription factors) necessary for elongation of fibres and root hairs. JA content analysis in cotton also confirmed that GbTCP has a profound effect on JA biosynthesis. In vitro ovule culture showed that an appropriate concentration of JA promoted fibre elongation. The results suggest that GbTCP is an important transcription factor for fibre and root hair development by regulating JA biosynthesis and response and other pathways, including reactive oxygen species, calcium channel and ethylene signalling.

AtMYB2 transcription factor can interact with the CMO promoter and regulate its downstream gene expression

The pC5 promoter, a region of the choline monooxygenase (CMO) promoter, contains an AtMYB2 transcription factor recognition sequence, TAACCA, and we have examined the interaction between AtMYB2 and the pC5 promoter. The AtMYB2 gene was cloned from Arabidopsis, expressed in Escherichia coli and transferred into pC5-GUS transgenic tobacco plants. Using an electrophoretic mobility shift assay, the AtMYB2 fusion protein binds to the TAACCA sequence in the pC5 promoter. As GUS activity was higher in pC5-GUS/AtMYB2 transgenic tobacco than in pC5-GUS plants, the AtMYB2 protein can bind to the CMO promoter in vitro and activate the transcription of the GUS reporter gene in vivo. The AtMYB2 transcription factor can therefore interact with the CMO promoter directly and regulate its downstream gene expression.

The draft genome of a diploid cotton Gossypium raimondii

We have sequenced and assembled a draft genome of G. raimondii, whose progenitor is the putative contributor of the D subgenome to the economically important fiber-producing cotton species Gossypium hirsutum and Gossypium barbadense. Over 73% of the assembled sequences were anchored on 13 G. raimondii chromosomes. The genome contains 40,976 protein-coding genes, with 92.2% of these further confirmed by transcriptome data. Evidence of the hexaploidization event shared by the eudicots as well as of a cotton-specific whole-genome duplication approximately 13-20 million years ago was observed. We identified 2,355 syntenic blocks in the G. raimondii genome, and we found that approximately 40% of the paralogous genes were present in more than 1 block, which suggests that this genome has undergone substantial chromosome rearrangement during its evolution. Cotton, and probably Theobroma cacao, are the only sequenced plant species that possess an authentic CDN1 gene family for gossypol biosynthesis, as revealed by phylogenetic analysis.

Epidermal cell differentiation in cotton mediated by the homeodomain leucine zipper gene, GhHD-1

Gossypium hirsutum L. (cotton) fibres are specialized trichomes a few centimetres in length that grow from the seed coat. Few genes directly involved in the differentiation of these epidermal cells have been identified. These include GhMYB25-like and GhMYB25, two related MYB transcription factors that regulate fibre cell initiation and expansion. We have also identified a putative homeodomain leucine zipper (HD-ZIP) transcription factor, GhHD-1, expressed in trichomes and early fibres that might play a role in cotton fibre initiation. Here, we characterize GhHD-1 homoeologues from tetraploid G. hirsutum and show, using reporter constructs and quantitative real-time PCR (qRT-PCR), that they are expressed predominantly in epidermal tissues during early fibre development, and in other tissues bearing epidermal trichomes. Silencing of GhHD-1 reduced trichome formation and delayed the timing of fibre initiation. Constitutive overexpression of GhHD-1 increased the number of fibres initiating on the seed, but did not affect leaf trichomes. Expression of GhHD-1 in cotton silenced for different fibre MYBs suggest that in ovules it acts downstream of GhMYB25-like, but is unaffected in GhMYB25- or GhMYB109-silenced plants. Microarray analysis of silencing and overexpression lines of GhHD-1 indicated that it potentially regulates the levels of ethylene and reactive oxidation species (ROS) through a WRKY transcription factor and calcium-signalling pathway genes to activate downstream genes necessary for cell expansion and elongation.

Overexpression of cotton (Gossypium hirsutum) dirigent1 gene enhances lignification that blocks the spread of Verticillium dahliae

Dirigent super-family abounds throughout the plant kingdom, especially vascular plants. To elucidate the function of cotton (Gossypium hirsutum) DIR genes in lignification, two cDNAs (designated GhDIR1 and GhDIR2) encoding putative dirigent proteins were isolated from cotton cDNA libraries. Real-time quantitative reverse transcription-polymerase chain reaction analysis revealed that GhDIR1 transcript was preferentially accumulated in cotton hypocotyls, whereas GhDIR2 was predominantly expressed in cotton fibers. Overexpression of GhDIR1 gene resulted in an increase in lignin content in transgenic cotton plants, compared with that of wild type. Histochemical assay revealed that the transgenic plants displayed more widespread lignification than that of wild type in epidermis and vascular bundle. Furthermore, the transgenic cotton plants displayed more tolerance to the infection of Verticillium dahliae. Our data suggest that GhDIR1 may be involved in cotton lignification which can block the spread of fungal pathogen V. dahliae.

Hyperspectral imaging techniques for rapid identification of Arabidopsis mutants with altered leaf pigment status

The spectral reflectance signature of living organisms provides information that closely reflects their physiological status. Because of its high potential for the estimation of geomorphic biological parameters, particularly of gross photosynthesis of plants, two-dimensional spectroscopy, via the use of hyperspectral instruments, has been widely used in remote sensing applications. In genetics research, in contrast, the reflectance phenotype has rarely been the subject of quantitative analysis; its potential for illuminating the pathway leading from the gene to phenotype remains largely unexplored. In this study, we employed hyperspectral imaging techniques to identify Arabidopsis mutants with altered leaf pigment status. The techniques are comprised of two modes; the first is referred to as the 'targeted mode' and the second as the 'non-targeted mode'. The 'targeted' mode is aimed at visualizing individual concentrations and compositional parameters of leaf pigments based on reflectance indices (RIs) developed for Chls a and b, carotenoids and anthocyanins. The 'non-targeted' mode highlights differences in reflectance spectra of leaf samples relative to reference spectra from the wild-type leaves. Through the latter approach, three mutant lines with weak irregular reflectance phenotypes, that are hardly identifiable by simple observation, were isolated. Analysis of these and other mutants revealed that the RI-based targeted pigment estimation was robust at least against changes in trichome density, but was confounded by genetic defects in chloroplast photorelocation movement. Notwithstanding such a limitation, the techniques presented here provide rapid and high-sensitive means to identify genetic mechanisms that coordinate leaf pigment status with developmental stages and/or environmental stress conditions.

Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana

The ultimate understanding of how biological diversity arises, is maintained, and lost depends on identifying the genes responsible. Although a good deal has been discovered about gene function over the past few decades, far less is understood about gene effects, that is, how natural variation in a gene contributes to natural variation in phenotypes. Trichome density in Arabidopsis thaliana is an ideal trait for studies of natural molecular and phenotypic variation, as trichome initiation is genetically well-characterized and trichome density is highly variable in and among natural populations. Here, we show that variation at GLABRA1 (GL1), an R2R3 MYB transcription factor gene, which has a role in trichome initiation, has qualitative and likely quantitative effects on trichome density in natural accessions of A. thaliana. Specifically, we characterize four independent loss-of-function alleles for GL1, each of which yields a glabrous phenotype. Further, we find that a pattern of common polymorphisms confined to the GL1 locus is associated with quantitative variation for trichome density. While mutations resulting in a glabrous phenotype are primarily coding changes, the pattern resulting in quantitative variation spans both coding and regulatory regions. These data show that GL1 is an important source of trichome density variation within A. thaliana and, along with recent reports, suggest that the TTG1 epidermal cell fate pathway generally may be the key molecular genetic source of natural trichome density variation and an important model for the study of molecular evolution.

GbPDF1 is involved in cotton fiber initiation via the core cis-element HDZIP2ATATHB2

Cotton (Gossypium spp.) fiber cells are seed trichomes derived from the epidermal layer of the cotton seed coat. The molecular components responsible for regulating fiber cell differentiation have not been fully elucidated. A cotton PROTODERMAL FACTOR1 gene (GbPDF1) was found to be expressed preferentially during fiber initiation and early elongation, with highest accumulation in fiber cells 5 d post anthesis. PDF1 silencing caused retardation of fiber initiation and produced shorter fibers and lower lint percentage compared with the wild type, indicating that the gene is required for cotton fiber development. Further analysis showed that a higher accumulation of hydrogen peroxide occurred in the RNA interference transgenic cotton lines. Meanwhile, the expression of several genes related to ethylene and pectin synthesis or sugar transport during cotton fiber growth was found to be significantly reduced in the PDF1-suppressed cotton. Three proteins interacting with GbPDF1 in yeast and in planta might involve cellular signaling or metabolism. GbPDF1 promoter::GUS constructs in transgenic cotton were predominantly expressed in the epidermis of ovules and developing fibers. Progressive deletions of the GbPDF1 promoter showed that a 236-bp promoter fragment was sufficient for basal GbPDF1 transcription in cotton. Mutation of putative regulatory sequences showed that HDZIP2ATATHB2, an element within the fragment, was essential for PGbPDF1-1 expression. The binding activity between this cis-element and nuclear extracts from fiber-bearing cotton ovules at 5 d post anthesis was specific. We conclude that GbPDF1 plays a critical role together with interaction partners in hydrogen peroxide homeostasis and steady biosynthesis of ethylene and pectin during fiber development via the core cis-element HDZIP2ATATHB2.

Improving peppermint essential oil yield and composition by metabolic engineering

Peppermint (Mentha × piperita L.) was transformed with various gene constructs to evaluate the utility of metabolic engineering for improving essential oil yield and composition. Oil yield increases were achieved by overexpressing genes involved in the supply of precursors through the 2C-methyl-D-erythritol 4-phosphate (MEP) pathway. Two-gene combinations to enhance both oil yield and composition in a single transgenic line were assessed as well. The most promising results were obtained by transforming plants expressing an antisense version of (+)-menthofuran synthase, which is critical for adjusting the levels of specific undesirable oil constituents, with a construct for the overexpression of the MEP pathway gene 1-deoxy-D-xylulose 5-phosphate reductoisomerase (up to 61% oil yield increase over wild-type controls with low levels of the undesirable side-product (+)-menthofuran and its intermediate (+)-pulegone). Elite transgenic lines were advanced to multiyear field trials, which demonstrated consistent oil yield increases of up to 78% over wild-type controls and desirable effects on oil composition under commercial growth conditions. The transgenic expression of a gene encoding (+)-limonene synthase was used to accumulate elevated levels of (+)-limonene, which allows oil derived from transgenic plants to be recognized during the processing of commercial formulations containing peppermint oil. Our study illustrates the utility of metabolic engineering for the sustainable agricultural production of high quality essential oils at a competitive cost.

A cotton gene encoding novel MADS-box protein is preferentially expressed in fibers and functions in cell elongation

Cotton fibers, as natural fibers, are widely used in the textile industry in the world. In order to find genes involved in fiber development, a cDNA (designated as GhMADS11) encoding a novel MADS protein with 151 amino acid residues was isolated from cotton fiber cDNA library. The deduced protein shares high similarity with Arabidopsis AP1 and AGL8 in MADS domain. However, the GhMADS11 protein (being absent of the partial K-domain and normal C-terminus) is shorter than AP1 and AGL8 by the reason of gene frameshift mutation during evolution. The experimental results revealed that GhMADS11 was not a transcriptional activator, and it did not form homodimer. GhMADS11 transcripts were specifically accumulated in elongating fibers, but no or very low signals of its expression were detected in other tissues of cotton. Overexpression of GhMADS11 in fission yeast promotes atypical cell elongation by 1.4-2.0-fold. Furthermore, morphological analysis indicated that the transformed cells expressing GhMADS11m, a MIKC-type derivative of GhMADS11 by the site-directed mutation, displayed the same phenotype as that of the transformed cells with GhMADS11. The concurrence of these data sets suggested that GhMADS11 protein may function in fiber cell elongation, and its MADS domain and partial K-domain are sufficient for this function.

Fine-mapping of the woolly gene controlling multicellular trichome formation and embryonic development in tomato

Trichomes are small hairs that originate from the epidermal cells of nearly all land plants, and they exist in unicellular and multicellular forms. The regulatory pathway of unicellular trichomes in Arabidopsis is well characterized. However, little is known about the multicellular trichome formation in tomato (Solanum lycopersicum). The woolly (Wo) gene controls multicellular trichome initiation and leads to embryonic lethality when homozygous in tomato. To clone and characterize Wo, the gene was fine-mapped to a DNA fragment of ~200 kb using the map-based cloning strategy. A series of sequence-based molecular markers, including simple sequence repeat, sequence characterized amplified region, and cleaved amplified polymorphic sequence were utilized in this study. Analysis of the sequence indicated that this region carries 19 putative open reading frames. These results will provide not only the important information for the isolation and characterization of Wo but also the starting point for studying the regulatory pathway responsible for trichome formation and embryonic lethality in tomato.

Activation of Arabidopsis seed hair development by cotton fiber-related genes

Each cotton fiber is a single-celled seed trichome or hair, and over 20,000 fibers may develop semi-synchronously on each seed. The molecular basis for seed hair development is unknown but is likely to share many similarities with leaf trichome development in Arabidopsis. Leaf trichome initiation in Arabidopsis thaliana is activated by GLABROUS1 (GL1) that is negatively regulated by TRIPTYCHON (TRY). Using laser capture microdissection and microarray analysis, we found that many putative MYB transcription factor and structural protein genes were differentially expressed in fiber and non-fiber tissues. Gossypium hirsutum MYB2 (GhMYB2), a putative GL1 homolog, and its downstream gene, GhRDL1, were highly expressed during fiber cell initiation. GhRDL1, a fiber-related gene with unknown function, was predominately localized around cell walls in stems, sepals, seed coats, and pollen grains. GFP:GhRDL1 and GhMYB2:YFP were co-localized in the nuclei of ectopic trichomes in siliques. Overexpressing GhRDL1 or GhMYB2 in A. thaliana Columbia-0 (Col-0) activated fiber-like hair production in 4–6% of seeds and had on obvious effects on trichome development in leaves or siliques. Co-overexpressing GhRDL1 and GhMYB2 in A. thaliana Col-0 plants increased hair formation in ∼8% of seeds. Overexpressing both GhRDL1 and GhMYB2 in A. thaliana Col-0 try mutant plants produced seed hair in ∼10% of seeds as well as dense trichomes inside and outside siliques, suggesting synergistic effects of GhRDL1 and GhMYB2 with try on development of trichomes inside and outside of siliques and seed hair in A. thaliana. These data suggest that a different combination of factors is required for the full development of trichomes (hairs) in leaves, siliques, and seeds. A. thaliana can be developed as a model a system for discovering additional genes that control seed hair development in general and cotton fiber in particular.

A regulatory gene induces trichome formation and embryo lethality in tomato

Trichomes are universal biological structures originating from the aerial epidermis, which serve as an excellent model to study plant differentiation at the cell level. Although the pathway regulating trichome formation in the Rosids has been well characterized, only very recently a few genes were identified for trichome initiation in the Asterids. In this study, we cloned Woolly (Wo), essential for trichome formation in tomato. Transgenic experiments revealed that the woolly phenotype is caused by the mutation in Wo which encodes a homeodomain protein containing a bZIP motif and a START domain. We identified three alleles of Wo and found that each allele contains a missense mutation, which respectively results in an amino acid substitution at the C terminus. Microarray and expression analysis showed that the expression of a B-type cyclin gene, SlCycB2, is possibly regulated by Wo, which also participates in trichome formation. Suppression of Wo or SlCycB2 expression by RNAi decreased the number of type I trichomes, and direct protein-protein interaction was detected between them, implying that both proteins may work together in the regulation of this type of trichome formation. Cytological observation and Wo transcript analysis in the developing seeds showed that embryo development was also correlated with Wo.

Differential regulation of a MYB transcription factor is correlated with transgenerational epigenetic inheritance of trichome density in Mimulus guttatus

• Epigenetic inheritance, transgenerational transmission of traits not proximally determined by DNA sequence, has been linked to transmission of chromatin modifications and gene regulation, which are known to be sensitive to environmental factors. Mimulus guttatus increases trichome (plant hair) density in response to simulated herbivore damage. Increased density is expressed in progeny even if progeny do not experience damage. To better understand epigenetic inheritance of trichome production, we tested the hypothesis that candidate gene expression states are inherited in response to parental damage. • Using M. guttatus recombinant inbred lines, offspring of leaf-damaged and control plants were raised without damage. Relative expression of candidate trichome development genes was measured in offspring. Line and parental damage effects on trichome density were measured. Associations between gene expression, trichome density, and response to parental damage were determined. • We identified M. guttatus MYB MIXTA-like 8 as a possible negative regulator of trichome development. We found that parental leaf damage induces down-regulation of MYB MIXTA-like 8 in progeny, which is associated with epigenetically inherited increased trichome density. • Our results link epigenetic transmission of an ecologically important trait with differential gene expression states - providing insight into a mechanism underlying environmentally induced 'soft inheritance'.

Evolutionary and comparative analysis of MYB and bHLH plant transcription factors

The expansion of gene families encoding regulatory proteins is typically associated with the increase in complexity characteristic of multi-cellular organisms. The MYB and basic helix-loop-helix (bHLH) families provide excellent examples of how gene duplication and divergence within particular groups of transcription factors are associated with, if not driven by, the morphological and metabolic diversity that characterize the higher plants. These gene families expanded dramatically in higher plants; for example, there are approximately 339 and 162 MYB and bHLH genes, respectively, in Arabidopsis, and approximately 230 and 111, respectively, in rice. In contrast, the Chlamydomonas genome has only 38 MYB genes and eight bHLH genes. In this review, we compare the MYB and bHLH gene families from structural, evolutionary and functional perspectives. The knowledge acquired on the role of many of these factors in Arabidopsis provides an excellent reference to explore sequence-function relationships in crops and other plants. The physical interaction and regulatory synergy between particular sub-classes of MYB and bHLH factors is perhaps one of the best examples of combinatorial plant gene regulation. However, members of the MYB and bHLH families also interact with a number of other regulatory proteins, forming complexes that either activate or repress the expression of sets of target genes that are increasingly being identified through a diversity of high-throughput genomic approaches. The next few years are likely to witness an increasing understanding of the extent to which conserved transcription factors participate at similar positions in gene regulatory networks across plant species.

GhMYB25-like: a key factor in early cotton fibre development

MYB transcription factors have been implicated in regulation of the development of ovule epidermal cells into the elongated seed fibres of cotton. An R2R3 MYB, GhMYB25-like, identified from its reduced expression in a fibreless mutant of cotton (Xu142 fl), is here shown to play a key role in the very early stages of fibre cell differentiation. A GhMYB25-like promoter-GUS construct was expressed predominantly in the epidermal layers of cotton ovules before anthesis (-3days post-anthesis, dpa), increasing in expression in 0-dpa ovules, primarily in those epidermal cells expanding into fibres, and then in elongating fibres at +3dpa, declining thereafter. This was consistent with GhMYB25-like transcript abundance during fibre development. RNA interference suppression of GhMYB25-like resulted in cotton plants with fibreless seeds, but normal trichomes elsewhere, phenocopying the Xu142 fl mutant. Like Xu142 fl these plants had reduced expression of the fibre-expressed MYBs, GhMYB25 and GhMYB109, indicating that GhMYB25-like is upstream from those MYBs. This hierarchy was supported by the absence of any change in transcript level of GhMYB25-like in GhMYB25- and GhMYB109-silenced transgenic lines. Transgenic cotton with an additional copy of the native gene had elevated expression of GhMYB25-like in ovules, but no obvious increase in fibre initials, suggesting that there may be other factors that interact with GhMYB25-like to differentiate epidermal cells into fibre cells.

Functional analysis of Gossypium hirsutum cellulose synthase catalytic subunit 4 promoter in transgenic Arabidopsis and cotton tissues

Gossypium hirsutum cellulose synthase catalytic subunit 4 (GhCesA4) plays an important role in cellulose biosynthesis during cotton fiber development. The transcript levels of GhCesA4 are significantly up-regulated as secondary cell wall cellulose is produced in developing cotton fibers. To understand the molecular mechanisms involved in transcriptional regulation of GhCesA4, β-glucuronidase (GUS) activity regulated by a GhCesA4 promoter (-2574/+56) or progressively deleted promoters were determined in both cotton tissues and transgenic Arabidopsis. The spatial regulation of GhCesA4 expression was similar between cotton tissues and transgenic Arabidopsis. GUS activity regulated by the GhCesA4 promoter (-2574/+56) was found in trichomes and root vascular tissues in both cotton and transgenic Arabidopsis. The -2574/-1824 region was responsible for up-regulation of GhCesA4 expression in trichomes and root vascular tissues in transgenic Arabidopsis. The -1824/-1355 region negatively regulated GhCesA4 expression in most Arabidopsis vascular tissues. For vascular expression in stems and leaves, the -898/-693 region was required. The -693/-320 region of the GhCesA4 promoter was necessary for basal expression of GhCesA4 in cotton roots as well as Arabidopsis roots. Exogenous phytohormonal treatments on transgenic Arabidopsis revealed that phytohormones may be involved in the differential regulation of GhCesA4 during cotton fiber development.

Cotton BCP genes encoding putative blue copper-binding proteins are functionally expressed in fiber development and involved in response to high-salinity and heavy metal stresses

Copper is vitally required for plants at low concentrations but extremely toxic for plants at elevated concentrations. Plants have evolved a series of mechanisms to prevent the consequences of the excess or deficit of copper. These mechanisms require copper-interacting proteins involved in copper trafficking. Blue copper-binding proteins (BCPs) are a class of copper proteins containing one blue copper-binding domain binding a single type I copper. To investigate the role of BCPs in plant development and in response to stresses, we isolated nine cDNAs encoding the putative blue copper-binding proteins (GhBCPs) from cotton (Gossypium hirsutum). Meanwhile, four corresponding genes (including GhBCP1-GhBCP4), which contain a single intron inserted in their conserved position, were isolated from cotton genome. Quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) analysis indicated that the nine GhBCP genes are differentially expressed in cotton tissues. Among them, GhBCP1 and GhBCP4 were predominantly expressed in fibers, while the transcripts of GhBCP2 and GhBCP3 were accumulated at relatively high levels in fibers. These four genes were strongly expressed in early fiber elongation, but dramatically declined with further fiber development. In addition, these GhBCP genes were upregulated in fibers by Cu(2+) , Zn(2+) , high-salinity and drought stresses, but downregulated in fibers by Al(3+) treatment. Overexpression of GhBCP1 and GhBCP4 in yeast (Schizosaccharomyces pombe) significantly increased the cell growth rate under Cu(2+) , Zn(2+) and high-salinity stresses. These results suggested that these GhBCPs may participate in the regulation of fiber development and in response to high-salinity and heavy metal stresses in cotton.

Polyploidization altered gene functions in cotton (Gossypium spp.)

Cotton (Gossypium spp.) is an important crop plant that is widely grown to produce both natural textile fibers and cottonseed oil. Cotton fibers, the economically more important product of the cotton plant, are seed trichomes derived from individual cells of the epidermal layer of the seed coat. It has been known for a long time that large numbers of genes determine the development of cotton fiber, and more recently it has been determined that these genes are distributed across At and Dt subgenomes of tetraploid AD cottons. In the present study, the organization and evolution of the fiber development genes were investigated through the construction of an integrated genetic and physical map of fiber development genes whose functions have been verified and confirmed. A total of 535 cotton fiber development genes, including 103 fiber transcription factors, 259 fiber development genes, and 173 SSR-contained fiber ESTs, were analyzed at the subgenome level. A total of 499 fiber related contigs were selected and assembled. Together these contigs covered about 151 Mb in physical length, or about 6.7% of the tetraploid cotton genome. Among the 499 contigs, 397 were anchored onto individual chromosomes. Results from our studies on the distribution patterns of the fiber development genes and transcription factors between the At and Dt subgenomes showed that more transcription factors were from Dt subgenome than At, whereas more fiber development genes were from At subgenome than Dt. Combining our mapping results with previous reports that more fiber QTLs were mapped in Dt subgenome than At subgenome, the results suggested a new functional hypothesis for tetraploid cotton. After the merging of the two diploid Gossypium genomes, the At subgenome has provided most of the genes for fiber development, because it continues to function similar to its fiber producing diploid A genome ancestor. On the other hand, the Dt subgenome, with its non-fiber producing D genome ancestor, provides more transcription factors that regulate the expression of the fiber genes in the At subgenome. This hypothesis would explain previously published mapping results. At the same time, this integrated map of fiber development genes would provide a framework to clone individual full-length fiber genes, to elucidate the physiological mechanisms of the fiber differentiation, elongation, and maturation, and to systematically study the functional network of these genes that interact during the process of fiber development in the tetraploid cottons.

Molecular cloning and expression analysis of a novel SANT/MYB gene from Gossypium barbadense

MYB family transcription factors are implicated in multiple developmental processes. Herein, a new full-length cDNA encoding a SANT/MYB transcription factor (designated as GbRL2) was cloned and characterized from cotton (Gossypium barbadense L.) for the first time. The full-length cDNA of GbRL2 was 573 bp with a 240 bp open reading frame (ORF) encoding a deduced protein of 80 amino acid polypeptide with a calculated molecular mass of 8.96 kDa and an isoelectric point of 8.96. Sequence alignment revealed that GbRL2 had high homology with other single SANT/MYB domain containing genes, including the RADIALIS genes in Antirrhinum majus and Bournea leiophylla. Semi-quantitative reverse transcript polymerase chain reaction (RT-PCR) revealed that at seedling stage, GbRL2 was strongly expressed in leaves but merely in stems. In opening flowers, the expression of GbRL2 was moderate in the petals but could not be detected in stamens. In ovules, the expression of GbRL2 could not be detected at -3 days post-anthesis (DPA) but increased during early elongation stage (0 DPA, +3 DPA, +5 DPA and +8 DPA). The transcripts of GbRL2 could also be detected at +8 DPA elongating fibers. We also examined the expression of RL2 gene in Gossypium hirstum cultivar Xu-142 and its fuzzless-lintless-seed mutant fl plants. The GhRL2 gene was ectopically expressed at -3 DPA in the fl mutant while the expression of GhRL2 in WT could not be detected. The expression of GhRL2 decreased early (+5 DPA) while that of WT was still strong. Our results suggest that GbRL2 may participate in development of various organs and may be a target for genetic improvement of cotton fiber.

Evidence that high activity of vacuolar invertase is required for cotton fiber and Arabidopsis root elongation through osmotic dependent and independent pathways, respectively

Vacuolar invertase (VIN) has long been considered as a major player in cell expansion. However, direct evidence for this view is lacking due, in part, to the complexity of multicellular plant tissues. Here, we used cotton (Gossypium spp.) fibers, fast-growing single-celled seed trichomes, to address this issue. VIN activity in elongating fibers was approximately 4-6-fold higher than that in leaves, stems, and roots. It was undetectable in fiberless cotton seed epidermis but became evident in initiating fibers and remained high during their fast elongation and dropped when elongation slowed. Furthermore, a genotype with faster fiber elongation had significantly higher fiber VIN activity and hexose levels than a slow-elongating genotype. By contrast, cell wall or cytoplasmic invertase activities did not show correlation with fiber elongation. To unravel the molecular basis of VIN-mediated fiber elongation, we cloned GhVIN1, which displayed VIN sequence features and localized to the vacuole. Once introduced to Arabidopsis (Arabidopsis thaliana), GhVIN1 complemented the short-root phenotype of a VIN T-DNA mutant and enhanced the elongation of root cells in the wild type. This demonstrates that GhVIN1 functions as VIN in vivo. In cotton fiber, GhVIN1 expression level matched closely with VIN activity and fiber elongation rate. Indeed, transformation of cotton fiber with GhVIN1 RNA interference or overexpression constructs reduced or enhanced fiber elongation, respectively. Together, these analyses provide evidence on the role of VIN in cotton fiber elongation mediated by GhVIN1. Based on the relative contributions of sugars to sap osmolality in cotton fiber and Arabidopsis root, we conclude that VIN regulates their elongation in an osmotic dependent and independent manner, respectively.

Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant

Cotton fiber as a single-celled trichome is a biological model system for studying cell differentiation and elongation. However, the complexity of its gene expression and regulatory mechanism allows only marginal progress. Here, we report the high-throughput tag-sequencing (Tag-seq) analysis using Solexa Genome Analyzer platform on transcriptome of -2 to 1 (fiber initiation, stage I) and 2-8 (fiber elongation, stage II) days post anthesis (DPA) cotton (Gossypium hirsutum) ovules (wild type: WT; Xuzhou 142 and its mutant: fuzzless/lintless or flM, in the same background). To this end, we sequenced 3.5-3.8 million tags representing 0.7-1.0 million unique transcripts for each library (WT1, WT2, M1, and M2). After removal of low quality tags, we obtained a total of 2,973,104, 3,139,306, 2,943,654, and 3,392,103 clean sequences that corresponded to 357,852, 280,787, 372,952, and 382,503 distinct tags for WT1, WT2, M1, and M2, respectively. All clean tags were aligned to the publicly available cotton transcript database (TIGR, http://www.tigr.org). About 15% of the distinct tags were uniquely mapped to the reference genes, and 31.4% of existing genes were matched by tags. The tag mapping to the database sequences generated 23,854, 24,442, 23,497, and 19,957 annotated genes for WT1, WT2, M1, and M2 libraries, respectively. Analyses of differentially expressed genes revealed the substantial changes in gene type and abundance between the wild type and mutant libraries. Among the 20 most differentially expressed genes in WT1/M1 and WT2/M2 libraries were cellulose synthase, phosphatase, and dehydrogenase, all of which are involved in the fiber cell development. Overall, the deep-sequencing analyses demonstrate the high degree of transcriptional complexity in early developing fibers and represent a major improvement over the microarrays for analyzing transcriptional changes on a large scale.

An L1 box binding protein, GbML1, interacts with GbMYB25 to control cotton fibre development

Transcription factors play key roles in plant development through their interaction with cis-elements and/or other transcription factors. A HD-Zip IV family transcription factor, Gossypium barbadense Meristem Layer 1 (GbML1) has been identified and characterized here. GbML1 specifically bound to the L1 box and the promoters of GbML1 and GbRDL1. GbML1 physically interacted with a key regulator of cotton fibre development, GbMYB25. Truncated and point mutation assays indicated the START-SAD domain was required for the binding to the C terminal domain (CTD) of GbMYB25. GbML1 overexpression in Arabidopsis increased the number of trichomes on stems and leaves and increased the accumulation of anthocyanin in leaves. Taken together, the L1 box binding protein, GbML1 was identified as the first partner for GbMYB25 and the role of START domain was discovered to be a protein binding domain in plants. Our findings will help the improvement of cotton fibre production and the understanding of the key role of HD-Zip family and MYB family in plants.

Temporal control of trichome distribution by microRNA156-targeted SPL genes in Arabidopsis thaliana

The production and distribution of plant trichomes is temporally and spatially regulated. After entering into the flowering stage, Arabidopsis thaliana plants have progressively reduced numbers of trichomes on the inflorescence stem, and the floral organs are nearly glabrous. We show here that SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) genes, which define an endogenous flowering pathway and are targeted by microRNA 156 (miR156), temporally control the trichome distribution during flowering. Plants overexpressing miR156 developed ectopic trichomes on the stem and floral organs. By contrast, plants with elevated levels of SPLs produced fewer trichomes. During plant development, the increase in SPL transcript levels is coordinated with the gradual loss of trichome cells on the stem. The MYB transcription factor genes TRICHOMELESS1 (TCL1) and TRIPTYCHON (TRY) are negative regulators of trichome development. We show that SPL9 directly activates TCL1 and TRY expression through binding to their promoters and that this activation is independent of GLABROUS1 (GL1). The phytohormones cytokinin and gibberellin were reported to induce trichome formation on the stem and inflorescence via the C2H2 transcription factors GIS, GIS2, and ZFP8, which promote GL1 expression. We show that the GIS-dependent pathway does not affect the regulation of TCL1 and TRY by miR156-targeted SPLs, represented by SPL9. These results demonstrate that the miR156-regulated SPLs establish a direct link between developmental programming and trichome distribution.

Characterization and fine mapping of the glabrous leaf and hull mutants (gl1) in rice (Oryza sativa L.)

The glabrous leaf and hull (gl1) mutants were isolated from M(2) generation of indica cultivar 93-11. These mutants produced smooth leaves and hairless glumes under normal growth conditions. By analyzing through scanning electron microscope, it was revealed that the leaf trichomes, including macro and micro hairs, were deficient in these mutants. Genetic analysis indicated that the mutation was controlled by a single recessive gene. Using nine SSR markers and one InDel marker, the gl1 gene was mapped between RM1200 and RM2010 at the short arm of chromosome 5, which was consistent with the mapping of gl1 in previous studies. To facilitate the map-based cloning of the gl1 gene, 12 new InDel markers were developed. A high-resolution genetic and physical map was constructed by using 1,396 mutant individuals of F(2) mapping population. Finally, the gl1 was fine mapped in 54-kb region containing 10 annotated genes. Cloning and sequencing of the target region from four gl1 mutants (gl1-1, gl1-2, gl1-3 and gl1-4) and four glabrous rice varieties (Jackson, Jefferson, Katy and Lemont) all showed that the same single point mutation (A-->T) occurred in the 5'-untranslated region (UTR) of the locus Os05g0118900 (corresponding to the 3'-UTR of STAR2). RT-PCR analysis of the locus Os05g0118900 revealed that its mRNA expression level was normal in gl1 mutant. RNA secondary structure prediction showed that the single point mutation resulted in a striking RNA conformational change. These results suggest that the single point mutation is most likely responsible for the glabrous leaf and hull phenotypes in rice.

Disruption of the 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) gene results in albino, dwarf and defects in trichome initiation and stomata closure in Arabidopsis

1-Deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) is an important enzyme involved in the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway which provides the basic five-carbon units for isoprenoid biosynthesis. To investigate the role of the MEP pathway in plant development and metabolism, we carried out detailed analyses on a dxr mutant (GK_215C01) and two DXR transgenic co-suppression lines, OX-DXR-L2 and OX-DXR-L7. We found that the dxr mutant was albino and dwarf. It never bolted, had significantly reduced number of trichomes and most of the stomata could not close normally in the leaves. The two co-suppression lines produced more yellow inflorescences and albino sepals with no trichomes. The transcription levels of genes involved in trichome initiation were found to be strongly affected, including GLABRA1, TRANSPARENT TESTA GLABROUS 1, TRIPTYCHON and SPINDLY, expression of which is regulated by gibberellic acids (GAs). Exogenous application of GA(3) could partially rescue the dwarf phenotype and the trichome initiation of dxr, whereas exogenous application of abscisic acid (ABA) could rescue the stomata closure defect, suggesting that lower levels of both GA and ABA contribute to the phenotype in the dxr mutants. We further found that genes involved in the biosynthetic pathways of GA and ABA were coordinately regulated. These results indicate that disruption of the plastidial MEP pathway leads to biosynthetic deficiency of photosynthetic pigments, GAs and ABA, and thus the developmental abnormalities, and that the flux from the cytoplasmic mevalonate pathway is not sufficient to rescue the deficiency caused by the blockage of the plastidial MEP pathway. These results reveal a critical role for the MEP biosynthetic pathway in controlling the biosynthesis of isoprenoids.

Identification and evaluation of new reference genes in Gossypium hirsutum for accurate normalization of real-time quantitative RT-PCR data

BACKGROUND

Normalizing through reference genes, or housekeeping genes, can make more accurate and reliable results from reverse transcription real-time quantitative polymerase chain reaction (qPCR). Recent studies have shown that no single housekeeping gene is universal for all experiments. Thus, suitable reference genes should be the first step of any qPCR analysis. Only a few studies on the identification of housekeeping gene have been carried on plants. Therefore qPCR studies on important crops such as cotton has been hampered by the lack of suitable reference genes.

RESULTS

By the use of two distinct algorithms, implemented by geNorm and NormFinder, we have assessed the gene expression of nine candidate reference genes in cotton: GhACT4, GhEF1alpha5, GhFBX6, GhPP2A1, GhMZA, GhPTB, GhGAPC2, GhbetaTUB3 and GhUBQ14. The candidate reference genes were evaluated in 23 experimental samples consisting of six distinct plant organs, eight stages of flower development, four stages of fruit development and in flower verticils. The expression of GhPP2A1 and GhUBQ14 genes were the most stable across all samples and also when distinct plants organs are examined. GhACT4 and GhUBQ14 present more stable expression during flower development, GhACT4 and GhFBX6 in the floral verticils and GhMZA and GhPTB during fruit development. Our analysis provided the most suitable combination of reference genes for each experimental set tested as internal control for reliable qPCR data normalization. In addition, to illustrate the use of cotton reference genes we checked the expression of two cotton MADS-box genes in distinct plant and floral organs and also during flower development.

CONCLUSION

We have tested the expression stabilities of nine candidate genes in a set of 23 tissue samples from cotton plants divided into five different experimental sets. As a result of this evaluation, we recommend the use of GhUBQ14 and GhPP2A1 housekeeping genes as superior references for normalization of gene expression measures in different cotton plant organs; GhACT4 and GhUBQ14 for flower development, GhACT4 and GhFBX6 for the floral organs and GhMZA and GhPTB for fruit development. We also provide the primer sequences whose performance in qPCR experiments is demonstrated. These genes will enable more accurate and reliable normalization of qPCR results for gene expression studies in this important crop, the major source of natural fiber and also an important source of edible oil. The use of bona fide reference genes allowed a detailed and accurate characterization of the temporal and spatial expression pattern of two MADS-box genes in cotton.

The essential role of GhPEL gene, encoding a pectate lyase, in cell wall loosening by depolymerization of the de-esterified pectin during fiber elongation in cotton

Cotton fiber elongation, largely achieved by cell wall loosening, is an important stage during cotton fiber development. In this present research, a fiber preferential cDNA encoding a pectate lyase (PEL) which could exclusively degrade the de-esterified pectin was isolated from a cotton (Gossypium hirsutum) fiber cDNA library. Subsequently, the corresponding PEL genes were isolated from four different cotton species and characterized. In vitro enzyme assays indicated that GhPEL really exhibited cleavage-activity against de-esterified pectin. The temporal-spatial expression analyses revealed that the GhPEL gene was preferentially expressed in fibers at 10 days-post anthesis (DPA). Antisense GhPEL transgenic cotton plants were generated by Agrobacterium-mediated transformation. Six homozygous lines, each with one or two copies of the transgene inserted as determined by southern blot analysis of the NPTII gene, were selected for further functional analysis. The GhPEL expression during fiber elongation in these transgenic lines was significantly suppressed in various degrees. Furthermore, the reduction of GhPEL enzymatic activity by decreasing GhPEL transcripts severely affected the degradation of de-esterified pectin in primary cell walls of transgenic cotton fibers, which consequently blocked cell wall loosening in early fiber development. Ultimately, the fiber elongation of all these transgenic lines was repressed. These results suggested that GhPEL may play an important role in the process of normal fiber elongation in cotton.

Functional analysis of GbAGL1, a D-lineage gene from cotton (Gossypium barbadense)

Cotton fibres originate from the outer ovule integument and D-lineage genes are essential for ovule development and their roles can be described by the 'ABCDE' model of flower development. To investigate the role of D-lineage genes during ovule and fibre development, GbAGL1 (GenBank accession number: FJ198049) was isolated from G. barbadense by using the SMART RACE strategy. Sequence and phylogenetic analyses revealed that GbAGL1 was a member of the D-lineage gene family. Southern blot analysis showed that GbAGL1 belonged to a low-copy gene family. Semi-quantitative RT-PCR and RNA in situ hybridization analyses revealed that the GbAGL1 gene in G. barbadense was highly expressed in whole floral bud primordia and the floral organs including ovules and fibres, but the signals were barely observed in vegetative tissues. GbAGL1 expression increased gradually with the ovule developmental stages. Over-expression of GbAGL1 in Arabidopsis caused obvious homeotic alternations in the floral organs, such as early flowering, and an extruded stigma, which were the typical phenotypes of the D-lineage gene family. In addition, a complementation test revealed that GbAGL1 could rescue the phenotypes of the stk mutant. Our study indicated that GbAGL1 was a D-lineage gene that was involved in ovule development and might play key roles in fibres development.

Apyrase (nucleoside triphosphate-diphosphohydrolase) and extracellular nucleotides regulate cotton fiber elongation in cultured ovules

Ectoapyrase enzymes remove the terminal phosphate from extracellular nucleoside tri- and diphosphates. In Arabidopsis (Arabidopsis thaliana), two ectoapyrases, AtAPY1 and AtAPY2, have been implicated as key modulators of growth. In fibers of cotton (Gossypium hirsutum), transcript levels for GhAPY1 and GhAPY2, two closely related ectoapyrases that have high sequence similarity to AtAPY1 and AtAPY2, are up-regulated when fibers enter their rapid growth phase. In an ovule culture system, fibers release ATP as they grow, and when their ectoapyrase activity is blocked by the addition of polyclonal anti-apyrase antibodies or by two different small molecule inhibitors, the medium ATP level rises and fiber growth is suppressed. High concentrations of the poorly hydrolyzable nucleotides ATPgammaS and ADPbetaS applied to the medium inhibit fiber growth, and low concentrations of them stimulate growth, but treatment with adenosine 5'-O-thiomonophosphate causes no change in the growth rate. Both the inhibition and stimulation of growth by applied nucleotides can be blocked by an antagonist that blocks purinoceptors in animal cells, and by adenosine. Treatment of cotton ovule cultures with ATPgammaS induces increased levels of ethylene, and two ethylene antagonists, aminovinylglycine and silver nitrate, block both the growth stimulatory and growth inhibitory effects of applied nucleotides. In addition, the ethylene precursor, 1-aminocyclopropane-1-carboxylic acid, lowers the concentration of nucleotide needed to promote fiber growth. These data indicate that ectoapyrases and extracellular nucleotides play a significant role in regulating cotton fiber growth and that ethylene is a likely downstream component of the signaling pathway.