Phylogenetically distinct cellulose synthase genes support secondary wall thickening in arabidopsis shoot trichomes and cotton fiber

Genome-wide analysis of the Sus gene family in cotton

Sucrose synthase (Sus) is a key enzyme in plant sucrose metabolism. In cotton, Sus (EC 2.4.1.13) is the main enzyme that degrades sucrose imported into cotton fibers from the phloem of the seed coat. This study demonstrated that the genomes of Gossypium arboreum L., G. raimondii Ulbr., and G. hirsutum L., contained 8, 8, and 15 Sus genes, respectively. Their structural organizations, phylogenetic relationships, and expression profiles were characterized. Comparisons of genomic and coding sequences identified multiple introns, the number and positions of which were highly conserved between diploid and allotetraploid cotton species. Most of the phylogenetic clades contained sequences from all three species, suggesting that the Sus genes of tetraploid G. hirsutum derived from those of its diploid ancestors. One Sus group (Sus I) underwent expansion during cotton evolution. Expression analyses indicated that most Sus genes were differentially expressed in various tissues and had development-dependent expression profiles in cotton fiber cells. Members of the same orthologous group had very similar expression patterns in all three species. These results provide new insights into the evolution of the cotton Sus gene family, and insight into its members' physiological functions during fiber growth and development.

Tertiary model of a plant cellulose synthase

A 3D atomistic model of a plant cellulose synthase (CESA) has remained elusive despite over forty years of experimental effort. Here, we report a computationally predicted 3D structure of 506 amino acids of cotton CESA within the cytosolic region. Comparison of the predicted plant CESA structure with the solved structure of a bacterial cellulose-synthesizing protein validates the overall fold of the modeled glycosyltransferase (GT) domain. The coaligned plant and bacterial GT domains share a six-stranded β-sheet, five α-helices, and conserved motifs similar to those required for catalysis in other GT-2 glycosyltransferases. Extending beyond the cross-kingdom similarities related to cellulose polymerization, the predicted structure of cotton CESA reveals that plant-specific modules (plant-conserved region and class-specific region) fold into distinct subdomains on the periphery of the catalytic region. Computational results support the importance of the plant-conserved region and/or class-specific region in CESA oligomerization to form the multimeric cellulose-synthesis complexes that are characteristic of plants. Relatively high sequence conservation between plant CESAs allowed mapping of known mutations and two previously undescribed mutations that perturb cellulose synthesis in Arabidopsis thaliana to their analogous positions in the modeled structure. Most of these mutation sites are near the predicted catalytic region, and the confluence of other mutation sites supports the existence of previously undefined functional nodes within the catalytic core of CESA. Overall, the predicted tertiary structure provides a platform for the biochemical engineering of plant CESAs.

Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species

Background

Cotton is the leading fiber crop worldwide. Gossypium barbadense is an important species of cotton because of its extra-long staple fibers with superior luster and silkiness. However, a systematic analysis and utilization of cDNA sequences from G. barbadense fiber development remains understudied.

Results

A total of 21,079 high quality sequences were generated from two non-normalized cDNA libraries prepared by using a mixture of G. barbadense Hai7124 fibers and ovules. After assembly processing, a set of 8,653 unigenes were obtained. Of those, 7,786 were matched to known proteins and 7,316 were assigned to functional categories. The molecular functions of these unigenes were mostly related to binding and catalytic activity, and carbohydrate, amino acid, and energy metabolisms were major contributors among the subsets of metabolism. Sequences comparison between G. barbadense and G. hirsutum revealed that 8,245 unigenes from G. barbadense were detected the similarity with those released publicly in G. hirsutum, however, the remaining 408 sequences had no hits against G. hirsutum unigenes database. Furthermore, 13,275 putative ESTs InDels loci involved in the orthologous and/or homoeologous differences between/within G. barbadense and G. hirsutum were discovered by in silico analyses, and 2,160 InDel markers were developed by ESTs with more than five insertions or deletions. By gel electrophoresis combined with sequencing verification, 71.11% candidate InDel loci were reconfirmed orthologous and/or homoeologous loci polymorphisms using G. hirsutum acc TM-1 and G. barbadense cv Hai7124. Blastx result showed among 2,160 InDel loci, 81 with significant function similarity with known genes associated with secondary wall synthesis process, indicating the important roles in fiber quality in tetraploid cultivated cotton species.

Conclusion

Sequence comparisons and InDel markers development will lay the groundwork for promoting the identification of genes related to superior agronomic traits, genetic differentiation and comparative genomic studies between G. hirsutum and G. barbadense.

Functional analysis of complexes with mixed primary and secondary cellulose synthases

In higher plants, cellulose is synthesized by cellulose synthase complexes, which contain multiple isoforms of cellulose synthases (CESAs). Among the total 10 CESA genes in Arabidopsis, recessive mutations at three of them cause the collapse of mature xylem cells in inflorescence stems of Arabidopsis (irx1(cesa8), irx3(cesa7) and irx5(cesa4)). These CESA genes are considered secondary cell wall CESAs. The others (the function CESA10 is still unknown) are thought to be specialized for cellulose synthesis in the primary cell wall. A split-ubiquitin membrane yeast two-hybrid system was used to assess interactions among four primary CESAs (CESA1, CESA2, CESA3, CESA6) and three secondary CESAs (CESA4, CESA7, CESA8). Our results showed that primary CESAs could physically interact with secondary CESAs in a limited fashion. Analysis of transgenic lines showed that CESA1 could partially rescue irx1(cesa8) null mutants, resulting in complementation of the plant growth defect, collapsed xylem and cellulose content deficiency. These results suggest that mixed primary and secondary CESA complexes are functional using experimental set-ups.

Cotton fiber cell walls of Gossypium hirsutum and Gossypium barbadense have differences related to loosely-bound xyloglucan

Cotton fiber is an important natural textile fiber due to its exceptional length and thickness. These properties arise largely through primary and secondary cell wall synthesis. The cotton fiber of commerce is a cellulosic secondary wall surrounded by a thin cuticulated primary wall, but there were only sparse details available about the polysaccharides in the fiber cell wall of any cotton species. In addition, Gossypium hirsutum (Gh) fiber was known to have an adhesive cotton fiber middle lamella (CFML) that joins adjacent fibers into tissue-like bundles, but it was unknown whether a CFML existed in other commercially important cotton fibers. We compared the cell wall chemistry over the time course of fiber development in Gh and Gossypium barbadense (Gb), the two most important commercial cotton species, when plants were grown in parallel in a highly controlled greenhouse. Under these growing conditions, the rate of early fiber elongation and the time of onset of secondary wall deposition were similar in fibers of the two species, but as expected the Gb fiber had a prolonged elongation period and developed higher quality compared to Gh fiber. The Gb fibers had a CFML, but it was not directly required for fiber elongation because Gb fiber continued to elongate rapidly after CFML hydrolysis. For both species, fiber at seven ages was extracted with four increasingly strong solvents, followed by analysis of cell wall matrix polysaccharide epitopes using antibody-based Glycome Profiling. Together with immunohistochemistry of fiber cross-sections, the data show that the CFML of Gb fiber contained lower levels of xyloglucan compared to Gh fiber. Xyloglucan endo-hydrolase activity was also higher in Gb fiber. In general, the data provide a rich picture of the similarities and differences in the cell wall structure of the two most important commercial cotton species.

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.

Complexes with mixed primary and secondary cellulose synthases are functional in Arabidopsis plants

In higher plants, cellulose is synthesized by so-called rosette protein complexes with cellulose synthases (CESAs) as catalytic subunits of the complex. The CESAs are divided into two distinct families, three of which are thought to be specialized for the primary cell wall and three for the secondary cell wall. In this article, the potential of primary and secondary CESAs forming a functional rosette complex has been investigated. The membrane-based yeast two-hybrid and biomolecular fluorescence systems were used to assess the interactions between three primary (CESA1, CESA3, CESA6), and three secondary (CESA4, CESA7, CESA8) Arabidopsis (Arabidopsis thaliana) CESAs. The results showed that all primary CESAs can physically interact both in vitro and in planta with all secondary CESAs. Although CESAs are broadly capable of interacting in pairwise combinations, they are not all able to form functional complexes in planta. Analysis of transgenic lines showed that CESA7 can partially rescue defects in the primary cell wall biosynthesis in a weak cesa3 mutant. Green fluorescent protein-CESA protein fusions revealed that when CESA3 was replaced by CESA7 in the primary rosette, the velocity of the mixed complexes was slightly faster than the native primary complexes. CESA1 in turn can partly rescue defects in secondary cell wall biosynthesis in a cesa8ko mutant, resulting in an increase of cellulose content relative to cesa8ko. These results demonstrate that sufficient parallels exist between the primary and secondary complexes for cross-functionality and open the possibility that mixed complexes of primary and secondary CESAs may occur at particular times.

Genome-wide transcriptomic analysis of cotton under drought stress reveal significant down-regulation of genes and pathways involved in fibre elongation and up-regulation of defense responsive genes

Cotton is an important source of natural fibre used in the textile industry and the productivity of the crop is adversely affected by drought stress. High throughput transcriptomic analyses were used to identify genes involved in fibre development. However, not much information is available on cotton genome response in developing fibres under drought stress. In the present study a genome wide transcriptome analysis was carried out to identify differentially expressed genes at various stages of fibre growth under drought stress. Our study identified a number of genes differentially expressed during fibre elongation as compared to other stages. High level up-regulation of genes encoding for enzymes involved in pectin modification and cytoskeleton proteins was observed at fibre initiation stage. While a large number of genes encoding transcription factors (AP2-EREBP, WRKY, NAC and C2H2), osmoprotectants, ion transporters and heat shock proteins and pathways involved in hormone (ABA, ethylene and JA) biosynthesis and signal transduction were up-regulated and genes involved in phenylpropanoid and flavonoid biosynthesis, pentose and glucuronate interconversions and starch and sucrose metabolism pathways were down-regulated during fibre elongation. This study showed that drought has relatively less impact on fibre initiation but has profound effect on fibre elongation by down-regulating important genes involved in cell wall loosening and expansion process. The comprehensive transcriptome analysis under drought stress has provided valuable information on differentially expressed genes and pathways during fibre development that will be useful in developing drought tolerant cotton cultivars without compromising fibre quality.

Cellulose synthase complexes: composition and regulation

Live cell imaging has greatly advanced our knowledge on the molecular mechanism by which cellulose is deposited. Both the actin and microtubule cytoskeleton are involved in assuring the proper distribution, organization, and dynamics of cellulose synthase complexes (CSCs). This review is an update on the most recent progress on the characterization of the composition, regulation, and trafficking of CSCs. With the newly identified cellulose synthase interactive protein 1 (CSI1) on hand, we begin to unveil the mystery of an intimate relationship between cellulose microfibrils and microtubules.

Cellulose synthesis in two secondary cell wall processes in a single cell type

Plant cells have a rigid cell wall that constrains internal turgor pressure yet extends in a regulated and organized manner to allow the cell to acquire shape. The primary load-bearing macromolecule of a plant cell wall is cellulose, which forms crystalline microfibrils that are organized with respect to a cell's function and shape requirements. A primary cell wall is deposited during expansion whereas secondary cell wall is synthesized post expansion during differentiation. A complex form of asymmetrical cellular differentiation occurs in Arabidopsis seed coat epidermal cells, where we have recently shown that two secondary cell wall processes occur that utilize different cellulose synthase (CESA) proteins. One process is to produce pectinaceous mucilage that expands upon hydration and the other is a radial wall thickening that reinforced the epidermal cell structure. Our data illustrate polarized specialization of CESA5 in facilitating mucilage attachment to the parent seed and CESA2, CESA5 and CESA9 in radial cell wall thickening and formation of the columella. Herein, we present a model for the complexity of cellulose biosynthesis in this highly differentiated cell type with further evidence supporting each cellulosic secondary cell wall process.

Perturbation of wood cellulose synthesis causes pleiotropic effects in transgenic aspen

Genetic manipulation of cellulose biosynthesis in trees may provide novel insights into the growth and development of trees. To explore this possibility, the overexpression of an aspen secondary wall-associated cellulose synthase (PtdCesA8) gene was attempted in transgenic aspen (Populus tremuloides L.) and unexpectedly resulted in silencing of the transgene as well as its endogenous counterparts. The main axis of the transgenic aspen plants quickly stopped growing, and weak branches adopted a weeping growth habit. Furthermore, transgenic plants initially developed smaller leaves and a less extensive root system. Secondary xylem (wood) of transgenic aspen plants contained as little as 10% cellulose normalized to dry weight compared to 41% cellulose typically found in normal aspen wood. This massive reduction in cellulose was accompanied by proportional increases in lignin (35%) and non-cellulosic polysaccharides (55%) compared to the 22% lignin and 36% non-cellulosic polysaccharides in control plants. The transgenic stems produced typical collapsed or 'irregular' xylem vessels that had altered secondary wall morphology and contained greatly reduced amounts of crystalline cellulose. These results demonstrate the fundamental role of secondary wall cellulose within the secondary xylem in maintaining the strength and structural integrity required to establish the vertical growth habit in trees.

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.

De novo assembled expressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq

Background

De novo assembly of transcript sequences produced by short-read DNA sequencing technologies offers a rapid approach to obtain expressed gene catalogs for non-model organisms. A draft genome sequence will be produced in 2010 for a Eucalyptus tree species (E. grandis) representing the most important hardwood fibre crop in the world. Genome annotation of this valuable woody plant and genetic dissection of its superior growth and productivity will be greatly facilitated by the availability of a comprehensive collection of expressed gene sequences from multiple tissues and organs.

Results

We present an extensive expressed gene catalog for a commercially grown E. grandis × E. urophylla hybrid clone constructed using only Illumina mRNA-Seq technology and de novo assembly. A total of 18,894 transcript-derived contigs, a large proportion of which represent full-length protein coding genes were assembled and annotated. Analysis of assembly quality, length and diversity show that this dataset represent the most comprehensive expressed gene catalog for any Eucalyptus tree. mRNA-Seq analysis furthermore allowed digital expression profiling of all of the assembled transcripts across diverse xylogenic and non-xylogenic tissues, which is invaluable for ascribing putative gene functions.

Conclusions

De novo assembly of Illumina mRNA-Seq reads is an efficient approach for transcriptome sequencing and profiling in Eucalyptus and other non-model organisms. The transcriptome resource (Eucspresso, http://eucspresso.bi.up.ac.za/) generated by this study will be of value for genomic analysis of woody biomass production in Eucalyptus and for comparative genomic analysis of growth and development in woody and herbaceous plants.

Regulated trafficking of cellulose synthases

New findings reveal that proteins involved in cellulose biosynthesis undergo regulated trafficking between intracellular compartments and the plasma membrane. The coordinated secretion and internalization of these proteins involve both the actin and cortical microtubule cytoskeletons. This regulated trafficking allows the dynamic remodeling of cellulose synthase complex (CSC) secretion during cell expansion and differentiation. Several new actors of the cellulose synthesis machinery have been recently identified.

Gene expression in developing fibres of Upland cotton (Gossypium hirsutum L.) was massively altered by domestication

Background

Understanding the evolutionary genetics of modern crop phenotypes has a dual relevance to evolutionary biology and crop improvement. Modern upland cotton (Gossypium hirsutum L.) was developed following thousands of years of artificial selection from a wild form, G. hirsutum var. yucatanense, which bears a shorter, sparser, layer of single-celled, ovular trichomes ('fibre'). In order to gain an insight into the nature of the developmental genetic transformations that accompanied domestication and crop improvement, we studied the transcriptomes of cotton fibres from wild and domesticated accessions over a developmental time course.

Results

Fibre cells were harvested between 2 and 25 days post-anthesis and encompassed the primary and secondary wall synthesis stages. Using amplified messenger RNA and a custom microarray platform designed to interrogate expression for 40,430 genes, we determined global patterns of expression during fibre development. The fibre transcriptome of domesticated cotton is far more dynamic than that of wild cotton, with over twice as many genes being differentially expressed during development (12,626 versus 5273). Remarkably, a total of 9465 genes were diagnosed as differentially expressed between wild and domesticated fibres when summed across five key developmental time points. Human selection during the initial domestication and subsequent crop improvement has resulted in a biased upregulation of components of the transcriptional network that are important for agronomically advanced fibre, especially in the early stages of development. About 15% of the differentially expressed genes in wild versus domesticated cotton fibre have no homology to the genes in databases.

Conclusions

We show that artificial selection during crop domestication can radically alter the transcriptional developmental network of even a single-celled structure, affecting nearly a quarter of the genes in the genome. Gene expression during fibre development within accessions and expression alteration arising from evolutionary change appears to be 'modular' - complex genic networks have been simultaneously and similarly transformed, in a coordinated fashion, as a consequence of human-mediated selection. These results highlight the complex alteration of the global gene expression machinery that resulted from human selection for a longer, stronger and finer fibre, as well as other aspects of fibre physiology that were not consciously selected. We illustrate how the data can be mined for genes that were unwittingly targeted by aboriginal and/or modern domesticators during crop improvement and/or which potentially control the improved qualities of domesticated cotton fibre.

See Commentary: http://www.biomedcentral.com/1741-7007/8/137