Phytohormonal networks promote differentiation of fiber initials on pre-anthesis cotton ovules grown in vitro and in planta.

The number of cotton (Gossypium sp.) ovule epidermal cells differentiating into fiber initials is an important factor affecting cotton yield and fiber quality. Despite extensive efforts in determining the molecular mechanisms regulating fiber initial differentiation, only a few genes responsible for fiber initial differentiation have been discovered. To identify putative genes directly involved in the fiber initiation process, we used a cotton ovule culture technique that controls the timing of fiber initial differentiation by exogenous phytohormone application in combination with comparative expression analyses between wild type and three fiberless mutants. The addition of exogenous auxin and gibberellins to pre-anthesis wild type ovules that did not have visible fiber initials increased the expression of genes affecting auxin, ethylene, ABA and jasmonic acid signaling pathways within 1 h after treatment. Most transcripts expressed differentially by the phytohormone treatment in vitro were also differentially expressed in the ovules of wild type and fiberless mutants that were grown in planta. In addition to MYB25-like, a gene that was previously shown to be associated with the differentiation of fiber initials, several other differentially expressed genes, including auxin/indole-3-acetic acid (AUX/IAA) involved in auxin signaling, ACC oxidase involved in ethylene biosynthesis, and abscisic acid (ABA) 8'-hydroxylase an enzyme that controls the rate of ABA catabolism, were co-regulated in the pre-anthesis ovules of both wild type and fiberless mutants. These results support the hypothesis that phytohormonal signaling networks regulate the temporal expression of genes responsible for differentiation of cotton fiber initials in vitro and in planta.

Construction of a plant-transformation-competent BIBAC library and genome sequence analysis of polyploid Upland cotton (Gossypium hirsutum L.).

BACKGROUND: Cotton, one of the world's leading crops, is important to the world's textile and energy industries, and is a model species for studies of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. Here, we report the construction of a plant-transformation-competent binary bacterial artificial chromosome (BIBAC) library and comparative genome sequence analysis of polyploid Upland cotton (Gossypium hirsutum L.) with one of its diploid putative progenitor species, G. raimondii Ulbr. RESULTS: We constructed the cotton BIBAC library in a vector competent for high-molecular-weight DNA transformation in different plant species through either Agrobacterium or particle bombardment. The library contains 76,800 clones with an average insert size of 135 kb, providing an approximate 99% probability of obtaining at least one positive clone from the library using a single-copy probe. The quality and utility of the library were verified by identifying BIBACs containing genes important for fiber development, fiber cellulose biosynthesis, seed fatty acid metabolism, cotton-nematode interaction, and bacterial blight resistance. In order to gain an insight into the Upland cotton genome and its relationship with G. raimondii, we sequenced nearly 10,000 BIBAC ends (BESs) randomly selected from the library, generating approximately one BES for every 250 kb along the Upland cotton genome. The retroelement Gypsy/DIRS1 family predominates in the Upland cotton genome, accounting for over 77% of all transposable elements. From the BESs, we identified 1,269 simple sequence repeats (SSRs), of which 1,006 were new, thus providing additional markers for cotton genome research. Surprisingly, comparative sequence analysis showed that Upland cotton is much more diverged from G. raimondii at the genomic sequence level than expected. There seems to be no significant difference between the relationships of the Upland cotton D- and A-subgenomes with the G. raimondii genome, even though G. raimondii contains a D genome (D5). CONCLUSIONS: The library represents the first BIBAC library in cotton and related species, thus providing tools useful for integrative physical mapping, large-scale genome sequencing and large-scale functional analysis of the Upland cotton genome. Comparative sequence analysis provides insights into the Upland cotton genome, and a possible mechanism underlying the divergence and evolution of polyploid Upland cotton from its diploid putative progenitor species, G. raimondii.

Cloning and characterization of homeologous cellulose synthase catalytic subunit 2 genes from allotetraploid cotton (Gossypium hirsutum L.).

Cellulose synthase catalytic subunits (CesAs) are the catalytic sites within a multisubunit complex for cellulose biosynthesis in plants. CesAs have been extensively studied in diploid plants, but are not well characterized in polyploid plants. Gossypium hirsutum is an allotetraploid cotton species producing over 90% of the world's cotton fibers. Although G. hirsutum CesAs (GhCesAs) are responsible for cellulose production in cotton fiber, very limited numbers of GhCesA genes have been identified. Here, we report isolating and characterizing a pair of homeologous CesA2 genes and their full-length cDNAs from allotetraploid cotton. The GhCesA2-A(T) gene from the A-subgenome and GhCesA2-D(T) gene from the D-subgenome were screened from a G. hirsutum BAC library. These genes shared 92% sequence similarity throughout the entire sequence. The coding sequences were nearly identical, and the deduced amino acid sequences from GhCesA2-A(T) (1,039 amino acids) and GhCesA2-D(T) (1,040 amino acids) were identical except four amino acids, whereas the noncoding sequences showed divergence. Sequence analyses showed that all exons of GhCesA2-A(T) contained consensus splice donor dinucleotides, but one exon in GhCesA2-D(T) contained nonconsensus splice donor dinucleotides. Although the nonconsensus splice donor dinucleotides were previously suggested to be involved in alternative splice or pseudogenization, our results showed that a majority of GhCesA2-A(T) and GhCesA2-D(T) transcripts consisted of functional and full-length transcripts with little evidence for alternative mRNA isoforms in developing cotton fibers. Expression analyses showed that GhCesA2-A(T) and GhCesA2-D(T) shared common temporal and spatial expression patterns, and they were highly and preferentially expressed during the cellulose biosynthesis stage in developing cotton fibers. The observations of higher expression levels of both GhCesA2-A(T) and GhCesA2-D(T) in developing fibers of one near-isogenic line (NIL) with higher fiber bundle strength over the other NIL with lower fiber bundle strength suggested that the differential expression of genes associated with secondary cell wall cellulose biosynthesis in developing fiber might affect cotton fiber properties.

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.

Near-isogenic cotton germplasm lines that differ in fiber-bundle strength have temporal differences in fiber gene expression patterns as revealed by comparative high-throughput profiling.

Gene expression profiles of developing cotton (Gossypium hirsutum L.) fibers from two near-isogenic lines (NILs) that differ in fiber-bundle strength, short-fiber content, and in fewer than two genetic loci were compared using an oligonucleotide microarray. Fiber gene expression was compared at five time points spanning fiber elongation and secondary cell wall (SCW) biosynthesis. Fiber samples were collected from field plots in a randomized, complete block design, with three spatially distinct biological replications for each NIL at each time point. Microarray hybridizations were performed in a loop experimental design that allowed comparisons of fiber gene expression profiles as a function of time between the two NILs. Overall, developmental expression patterns revealed by the microarray experiment agreed with previously reported cotton fiber gene expression patterns for specific genes. Additionally, genes expressed coordinately with the onset of SCW biosynthesis in cotton fiber correlated with gene expression patterns of other SCW-producing plant tissues. Functional classification and enrichment analysis of differentially expressed genes between the two NILs revealed that genes associated with SCW biosynthesis were significantly up-regulated in fibers of the high-fiber quality line at the transition stage of cotton fiber development. For independent corroboration of the microarray results, 15 genes were selected for quantitative reverse transcription PCR analysis of fiber gene expression. These analyses, conducted over multiple field years, confirmed the temporal difference in fiber gene expression between the two NILs. We hypothesize that the loci conferring temporal differences in fiber gene expression between the NILs are important regulatory sequences that offer the potential for more targeted manipulation of cotton fiber quality.

Genome-wide analysis reveals rapid and dynamic changes in miRNA and siRNA sequence and expression during ovule and fiber development in allotetraploid cotton (Gossypium hirsutum L.).

BACKGROUND: Cotton fiber development undergoes rapid and dynamic changes in a single cell type, from fiber initiation, elongation, primary and secondary wall biosynthesis, to fiber maturation. Previous studies showed that cotton genes encoding putative MYB transcription factors and phytohormone responsive factors were induced during early stages of ovule and fiber development. Many of these factors are targets of microRNAs (miRNAs) that mediate target gene regulation by mRNA degradation or translational repression. RESULTS: Here we sequenced and analyzed over 4 million small RNAs derived from fiber and non-fiber tissues in cotton. The 24-nucleotide small interfering RNAs (siRNAs) were more abundant and highly enriched in ovules and fiber-bearing ovules relative to leaves. A total of 31 miRNA families, including 27 conserved, 4 novel miRNA families and a candidate-novel miRNA, were identified in at least one of the cotton tissues examined. Among 32 miRNA precursors representing 19 unique miRNA families identified, 7 were previously reported, and 25 new miRNA precursors were found in this study. Sequencing, miRNA microarray, and small RNA blot analyses showed a trend of repression of miRNAs, including novel miRNAs, during ovule and fiber development, which correlated with upregulation of several target genes tested. Moreover, 223 targets of cotton miRNAs were predicted from the expressed sequence tags derived from cotton tissues, including ovules and fibers. The cotton miRNAs examined triggered cleavage in the predicted sites of the putative cotton targets in ovules and fibers. CONCLUSIONS: Enrichment of siRNAs in ovules and fibers suggests active small RNA metabolism and chromatin modifications during fiber development, whereas general repression of miRNAs in fibers correlates with upregulation of a dozen validated miRNA targets encoding transcription and phytohormone response factors, including the genes found to be highly expressed in cotton fibers. Rapid and dynamic changes in siRNAs and miRNAs may contribute to ovule and fiber development in allotetraploid cotton.

Phenotypic and molecular evaluation of cotton hairy roots as a model system for studying nematode resistance.

Agrobacterium rhizogenes-induced cotton (Gossypium hirsutum L.) hairy roots were evaluated as a model system for studying molecular cotton-nematode interactions. Hairy root cultures were developed from the root-knot nematode (RKN) (Meloidogyne incognita [Kofoid and White] Chitwood, race 3)-resistant breeding line M315 and from the reniform nematode (RN) (Rotylenchulus reniformis Linford & Oliveira)-resistant accession GB713 (G. barbadense L.) and compared to a nematode-susceptible culture derived from the obsolete cultivar DPL90. M315, GB713, and DPL90 hairy roots differed significantly in their appearance and growth potential; however, these differences were not correlated with transcript levels of the A. rhizogenes T-DNA genes rolB and aux2 which help regulate hairy root initiation and proliferation. DPL90 hairy roots were found to support both RKN and RN reproduction in tissue culture, whereas M315 and GB713 hairy roots were resistant to RKN and RN, respectively. M315 hairy roots showed constitutive up-regulation of the defense gene MIC3 (Meloidogyne Induced Cotton3) compared to M315 whole-plant roots and DPL90 hairy roots. Our data show the potential use of cotton hairy roots in maintaining monoxenic RKN and RN cultures and suggest hairy roots may be useful in evaluating the effect of manipulated host gene expression on nematode resistance in cotton.

Salt tolerance of Arabidopsis thaliana requires maturation of N-glycosylated proteins in the Golgi apparatus.

Protein N-glycosylation in the endoplasmic reticulum (ER) and in the Golgi apparatus is an essential process in eukaryotic cells. Although the N-glycosylation pathway in the ER has been shown to regulate protein quality control, salt tolerance, and cellulose biosynthesis in plants, no biological roles have been linked functionally to N-glycan modifications that occur in the Golgi apparatus. Herein, we provide evidence that mutants defective in N-glycan maturation, such as complex glycan 1 (cgl1), are more salt-sensitive than wild type. Salt stress caused growth inhibition, aberrant root-tip morphology, and callose accumulation in cgl1, which were also observed in an ER oligosaccharyltransferase mutant, staurosporin and temperature sensitive 3a (stt3a). Unlike stt3a, cgl1 did not cause constitutive activation of the unfolded protein response. Instead, aberrant modification of the plasma membrane glycoprotein KORRIGAN 1/RADIALLY SWOLLEN 2 (KOR1/RSW2) that is necessary for cellulose biosynthesis occurred in cgl1 and stt3a. Genetic analyses identified specific interactions among rsw2, stt3a, and cgl1 mutations, indicating that the function of KOR1/RSW2 protein depends on complex N-glycans. Furthermore, cellulose deficient rsw1-1 and rsw2-1 plants were also salt-sensitive. These results establish that plant protein N-glycosylation functions beyond protein folding in the ER and is necessary for sufficient cell-wall formation under salt stress.

Cu/Zn superoxide dismutases in developing cotton fibers: evidence for an extracellular form.

Hydrogen peroxide and other reactive oxygen species are important signaling molecules in diverse physiological processes. Previously, we discovered superoxide dismutase (SOD) activity in extracellular protein preparations from fiber-bearing cotton (Gossypium hirsutum L.) seeds. We show here, based on immunoreactivity, that the enzyme is a Cu/Zn-SOD (CSD). Immunogold localization shows that CSD localizes to secondary cell walls of developing cotton fibers. Five cotton CSD cDNAs were cloned from cotton fiber and classified into three subfamilies (Group 1: GhCSD1; Group 2: GhCSD2a and GhCSD2b; Group 3: GhCSD3 and GhCSD3s). Members of Group 1 and 2 are expressed throughout fiber development, but predominant during the elongation stage. Group 3 CSDs are also expressed throughout fiber development, but transiently increase in abundance at the transition period between cell elongation and secondary cell wall synthesis. Each of the three GhCSDs also has distinct patterns of expression in tissues other than fiber. Overexpression of cotton CSDs fused to green fluorescent protein in transgenic Arabidopsis demonstrated that GhCSD1 localizes to the cytosol, GhCSD2a localizes to plastids, and GhCSD3 is translocated to the cell wall. Subcellular fractionation of proteins from transgenic Arabidopsis seedlings confirmed that only c-myc epitope-tagged GhCSD3 co-purifies with cell wall proteins. Extracellular CSDs have been suggested to be involved in lignin formation in secondary cell walls of other plants. Since cotton fibers are not lignified, we suggest that extracellular CSDs may be involved in other plant cell wall growth and development processes.

Involvement of extracellular Cu/Zn superoxide dismutase in cotton fiber primary and secondary cell wall biosynthesis.

Extracellular Cu/Zn superoxide dismutases (CSDs) that catalyze the conversion of superoxide to hydrogen peroxide have been suggested to be involved in lignification of secondary walls in spinach, pine and aspen. In cotton fibers, hydrogen peroxide was proposed to be involved in the induction of secondary cell wall biosynthesis. Recently, we identified extracellular CSDs from developing cotton fibers using both immunological and epitope tagging techniques. Since cotton fibers are not lignified, we suggested that extracellular CSDs may be involved in plant cell wall growth and development processes other than lignification. In this addendum, we have further characterized the extracellular CSD in cotton fiber. Immunoblots, enzyme activity assays, and transcript levels show that an extracellular CSD is present in elongating primary walls as well as thickening secondary walls of cotton fibers. Our working model proposes that extracellular hydrogen peroxide levels, regulated by redox status-related enzymes including extracellular CSDs and peroxidases, may affect the processes of wall loosening and wall tightening.

Accumulation of genome-specific transcripts, transcription factors and phytohormonal regulators during early stages of fiber cell development in allotetraploid cotton.

Gene expression during the early stages of fiber cell development and in allopolyploid crops is poorly understood. Here we report computational and expression analyses of 32 789 high-quality ESTs derived from Gossypium hirsutum L. Texas Marker-1 (TM-1) immature ovules (GH_TMO). The ESTs were assembled into 8540 unique sequences including 4036 tentative consensus sequences (TCs) and 4504 singletons, representing approximately 15% of the unique sequences in the cotton EST collection. Compared with approximately 178 000 existing ESTs derived from elongating fibers and non-fiber tissues, GH_TMO ESTs showed a significant increase in the percentage of genes encoding putative transcription factors such as MYB and WRKY and genes encoding predicted proteins involved in auxin, brassinosteroid (BR), gibberellic acid (GA), abscisic acid (ABA) and ethylene signaling pathways. Cotton homologs related to MIXTA, MYB5, GL2 and eight genes in the auxin, BR, GA and ethylene pathways were induced during fiber cell initiation but repressed in the naked seed mutant (N1N1) that is impaired in fiber formation. The data agree with the known roles of MYB and WRKY transcription factors in Arabidopsis leaf trichome development and the well-documented phytohormonal effects on fiber cell development in immature cotton ovules cultured in vitro. Moreover, the phytohormonal pathway-related genes were induced prior to the activation of MYB-like genes, suggesting an important role of phytohormones in cell fate determination. Significantly, AA sub-genome ESTs of all functional classifications including cell-cycle control and transcription factor activity were selectively enriched in G. hirsutum L., an allotetraploid derived from polyploidization between AA and DD genome species, a result consistent with the production of long lint fibers in AA genome species. These results suggest general roles for genome-specific, phytohormonal and transcriptional gene regulation during the early stages of fiber cell development in cotton allopolyploids.

Characterization of GhRac1 GTPase expressed in developing cotton (Gossypium hirsutum L.) fibers.

Cytoskeleton assembly plays an important role in determining cotton fiber cell length and morphology and is developmentally regulated. As in other plant cells, it is not clear how cytoskeletal assembly in fibers is regulated. Recently, several Rac/Rop GTPases in Arabidopsis were shown to regulate isotropic and polar cell growth of root hairs and pollen tubes by controlling assembly of the cytoskeleton. GhRac1, isolated from cottonseeds, is a member of the Rac/Rop GTPase family and is abundantly expressed in rapidly growing cotton tissues. GhRac1 shows the greatest sequence similarity to the group IV subfamily of Arabidopsis Rac/Rop genes. Overexpression of GhRac1 in E. coli led to the production of a functional GTPase as shown by in vitro enzyme activity assay. In contrast to other Rac/Rop GTPases found in cotton fiber, GhRac1 is highly expressed during the elongation stage of fiber development with expression decreasing dramatically when the rate of fiber elongation declines. The association of highest GhRac1 expression during stages of maximal cotton fiber elongation suggests that GhRac1 GTPase may be a potential regulator of fiber elongation by controlling cytoskeletal assembly.

Cotton-fiber germin-like protein. II: Immunolocalization, purification, and functional analysis.

Cotton (Gossypium hirsutum L.) contains a germin-like protein (GLP), GhGLP1, that shows tissue-specific accumulation in fiber. The fiber GLP is an oligomeric, glycosylated protein with a subunit size of approximately 25.5 kDa. Accumulation of GhGLP1 occurs during the period of fiber elongation [4-14 days post-anthesis (DPA)]. During early phases of fiber development (2-4 DPA), GhGLP1 localizes to cytoplasmic vesicles as shown by confocal immunofluorescent microscopy. In slightly older fibers (7-10 DPA), GhGLP1 localizes to the apoplast. In other plants, germins and GLPs have been reported to have enzymatic activities including oxalate oxidase (OxO), superoxide dismutase, and ADP-glucose pyrophosphatase. Cotton fiber extracts did not contain OxO activity, nor did intact fibers stain for OxO activity. A four-step purification protocol involving ammonium sulfate precipitation of a 1.0 M NaCl extract, ion-exchange chromatography on DEAE-Trisacryl M, lectin-affinity chromatography, and gel filtration chromatography resulted in electrophoretically pure GhGLP1. While 1.0 M NaCl extracts from 10-14 DPA fiber contained superoxide dismutase and phosphodiesterase activities, GhGLP1 could be separated from both enzyme activities by the purification protocol. Although a GLP accumulates in the cotton fiber apoplast during cell elongation, the function of this protein in fiber growth and development remains unknown.

Cotton fiber germin-like protein. I. Molecular cloning and gene expression.

The presence of cotton ( Gossypium hirsutum L.) fiber transcripts coding for a germin-like protein (GLP) was revealed by differential display analysis in which early stages of cotton fiber development between a wild type line, Texas Marker-1 (TM1) and a near isogenic mutant, Naked Seed (N1) were compared. Transcripts of the cotton GLP ( GhGLP1) accumulated specifically in TM1, but did not accumulate in the mutant although the GhGLP1 gene was present in both lines. The deduced protein sequence of GhGLP1 is similar to Prunus persica auxin-binding proteins, a barley ADP-glucose pyrophosphatase/phosphodiesterase and two different classes of hydrogen peroxide-producing enzymes: wheat germin oxalate oxidase and moss extracellular Mn-superoxide dismutase. Cotton GLPs constitute a multigene family like those of Arabidopsis, rice, soybean, and barley. GhGLP1 transcripts accumulated to their highest levels during the period of fiber expansion, followed by a sharp decline when the rate of cell expansion decreased. While germins and GLPs appear to be involved in defense mechanisms in some plants, both biotic and abiotic stress down-regulated the expression of GhGLP1. Numerous functions have been proposed for dicot GLPs. However, to date, there is little direct evidence for how these proteins function in vivo. The association of maximal GhGLP1 expression with stages of maximal cotton fiber elongation suggests that some GLPs may be important for cell wall expansion.