Molecular and biochemical evidence for phenylpropanoid synthesis and presence of wall-linked phenolics in cotton fibers

Systematic characterization of cinnamyl alcohol dehydrogenase members revealed classification and function divergence in Haplomitrium mnioides

Cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.195) is considered to be a key enzyme in lignin biosynthesis, which can catalyze cinnamyl aldehyde to produce cinnamyl alcohol. In this study, three putative CADs were characterized from the liverwort Haplomitrium mnioides. The sequence alignment and phylogenetic analysis revealed that HmCADs belonged to a multigene family, with three HmCADs belonging to class II, class III, and class IV, respectively. In vitro enzymatic studies demonstrated that HmCAD2 exhibited high affinity and catalytic activity towards five cinnamyl aldehydes, followed by HmCAD3 with poor catalytic activity, and HmCAD1 catalyzed only the reaction of p-coumaryl aldehyde and coniferyl aldehyde with extremely low catalytic capacity. Protein-substrate binding simulations were performed to investigate the differences in catalytic activity exhibited when proteins catalyzed different substrates. Furthermore, distinct expression patterns of three HmCADs were identified in different plant tissues. Subcellular localization tests confirmed that HmCAD1/2/3 was located in the cytoplasm. The simulated responses of HmCADs to different stresses showed that HmCAD1 played a positive role in coping with each stress, while HmCAD2/3 was weak. These findings demonstrate the diversity of CADs in liverwort, highlight the divergent role of HmCAD1/2/3 in substrate catalysis, and also suggest their possible involvement in stress response, thereby providing new insights into CAD evolution while emphasizing their potential distinctive and collaborative contributions to the normal growth of primitive liverworts.

Regulation of lignin biosynthesis by GhCAD37 affects fiber quality and anther vitality in upland cotton

Cotton stands as a pillar in the textile industry due to its superior natural fibers. Lignin, a complex polymer synthesized from phenylalanine and deposited in mature cotton fibers, is believed to be essential for fiber quality, although the precise effects remain largely unclear. In this study, we characterized two ubiquitously expressed cinnamyl alcohol dehydrogenases (CAD), GhCAD37A and GhCAD37D (GhCAD37A/D), in Gossypium hirsutum. GhCAD37A/D possess CAD enzymatic activities, to catalyze the generation of monolignol products during lignin biosynthesis. Analysis of transgenic cotton knockout and overexpressing plants revealed that GhCAD37A/D are important regulators of fiber quality, positively impacting breaking strength but negatively affecting fiber length and elongation percentage by modulating lignin biosynthesis in fiber cells. Moreover, GhCAD37A/D are shown to modulate anther vitality and affect stem lodging trait in cotton by influencing lignin biosynthesis in the vascular bundles of anther and stem, respectively. Additionally, our study revealed that Ghcad37A/D knockout plants displayed red stem xylem, likely due to the overaccumulation of aldehyde intermediates in the phenylpropanoid metabolism pathway, as indicated by metabolomics analysis. Thus, our work illustrates that GhCAD37A/D are two important enzymes of lignin biosynthesis in different cotton organs, influencing fiber quality, anther vitality, and stem lodging.

GhMYB52 Like: A Key Factor That Enhances Lint Yield by Negatively Regulating the Lignin Biosynthesis Pathway in Fibers of Upland Cotton (Gossypium hirsutum L.)

In the context of sustainable agriculture and biomaterial development, understanding and enhancing plant secondary cell wall formation are crucial for improving crop fiber quality and biomass conversion efficiency. This is especially critical for economically important crops like upland cotton (Gossypium hirsutum L.), for which fiber quality and its processing properties are essential. Through comprehensive genome-wide screening and analysis of expression patterns, we identified a particularly high expression of an R2R3 MYB transcription factor, GhMYB52 Like, in the development of the secondary cell wall in cotton fiber cells. Utilizing gene-editing technology to generate a loss-of-function mutant to clarify the role of GhMYB52 Like, we revealed that GhMYB52 Like does not directly contribute to cellulose synthesis in cotton fibers but instead represses a subset of lignin biosynthesis genes, establishing it as a lignin biosynthesis inhibitor. Concurrently, a substantial decrease in the lint index, a critical measure of cotton yield, was noted in parallel with an elevation in lignin levels. This study not only deepens our understanding of the molecular mechanisms underlying cotton fiber development but also offers new perspectives for the molecular improvement of other economically important crops and the enhancement of biomass energy utilization.

Effect of the HCT Gene on Lignin Synthesis and Fiber Development in Gossypium barbadense

As one of the key enzymes in the metabolic pathway of phenylpropane, shikimate hydroxycinnamoyl transferase (HCT) is mainly involved in the biosynthesis of the plant secondary cell wall, which is closely related to cotton fiber quality. In this study, whole-genome identification and bioinformatics analysis of the HCT gene family were performed in G. barbadense. In the whole genome, we identified 136 GbHCT genes encoding 309-504 amino acids. Phylogenetic analysis divided the genome into 5 subfamilies, which were located on 25 chromosomes. Collinear analysis of polyploidization and tandem duplication events were the main driving forces for the rapid expansion and evolution of this family, and the genes underwent loose purifying selection constraints after duplication. Gene promoters identified a variety of cis-acting elements related to plant hormones and the stress response. Several members of the GbHCT family were highly expressed during the development of cotton fiber, and different members had different expression patterns in cotton fiber. After GbHCT114 gene silencing in cotton, the amount of stem surface trichomes and lignin content decreased, and the cell morphology and arrangement changed. After the GbHCT114 gene was overexpressed in Arabidopsis thaliana (L.) Heynh., the number of stem and leaf surface trichomes and the cross-sectional area of the secondary xylem duct cell wall increased. In addition, utilizing transcriptomic analysis, differentially expressed genes associated with lignin synthesis and fiber development were identified. Taken together, the results obtained in this study confirm that the GbHCT114 gene regulates plant trichome development, which lays a theoretical foundation for future studies on the function of GbHCT114 in cotton.

Comparative proteomic analysis identified proteins and the phenylpropanoid biosynthesis pathway involved in the response to ABA treatment in cotton fiber development

Abscisic acid (ABA) is a plant hormone that plays an important role in cotton fiber development. In this study, the physiological changes and proteomic profiles of cotton (Gossypium hirsutum) ovules were analyzed after 20 days of ABA or ABA inhibitor (ABAI) treatment. The results showed that compared to the control (CK), the fiber length was significantly decreased under ABA treatment and increased under ABAI treatment. Using a tandem mass tags-based quantitative technique, the proteomes of cotton ovules were comprehensively analyzed. A total of 7321 proteins were identified, of which 365 and 69 differentially accumulated proteins (DAPs) were identified in ABA versus CK and ABAI versus CK, respectively. Specifically, 345 and 20 DAPs were up- and down-regulated in the ABA group, and 65 and 4 DAPs were up- and down-regulated in the ABAI group, respectively. The DAPs in the ABA group were mainly enriched in the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis and flavonoid secondary metabolism, whereas the DAPs in the ABAI group were mainly enriched in the indole alkaloid biosynthesis and phenylpropanoid biosynthesis pathways. Moreover, 9 proteins involved in phenylpropanoid biosynthesis were upregulated after ABA treatment, suggesting that this pathway might play important roles in the response to ABA, and 3 auxin-related proteins were upregulated, indicating that auxin might participate in the regulation of fiber development under ABAI treatment.

Phenomics and transcriptomics analyses reveal deposition of suberin and lignin in the short fiber cell walls produced from a wild cotton species and two mutants

Fiber length is one of the major properties determining the quality and commercial value of cotton. To understand the mechanisms regulating fiber length, genetic variations of cotton species and mutants producing short fibers have been compared with cultivated cottons generating long and normal fibers. However, their phenomic variation other than fiber length has not been well characterized. Therefore, we compared physical and chemical properties of the short fibers with the long fibers. Fiber characteristics were compared in two sets: 1) wild diploid Gossypium raimondii Ulbrich (short fibers) with cultivated diploid G. arboreum L and tetraploid G. hirsutum L. (long fibers); 2) G. hirsutum short fiber mutants, Ligon-lintless 1 (Li1) and 2 (Li2) with their near isogenic line (NIL), DP-5690 (long fibers). Chemical analyses showed that the short fibers commonly consisted of greater non-cellulosic components, including lignin and suberin, than the long fibers. Transcriptomic analyses also identified up-regulation of the genes related to suberin and lignin biosynthesis in the short fibers. Our results may provide insight on how high levels of suberin and lignin in cell walls can affect cotton fiber length. The approaches combining phenomic and transcriptomic analyses of multiple sets of cotton fibers sharing a common phenotype would facilitate identifying genes and common pathways that significantly influence cotton fiber properties.

Phylogeny and functional characterization of the cinnamyl alcohol dehydrogenase gene family in Phryma leptostachya

Phryma leptostachya has attracted increasing attention because it is rich in furofuran lignans with a wide range of biological activities. Biosynthesis of furofuran lignans begins with the dimerization of coniferyl alcohol, one of the monolignol. Cinnamyl alcohol dehydrogenase (CAD) catalyzes the final step of monolignol biosynthesis, reducing cinnamyl aldehydes to cinnamyl alcohol. As it is in the terminal position of monolignol biosynthesis, its type and activity can cause significant changes in the total amount and composition of lignans. Herein, combined with bioinformatics analysis and in vitro enzyme assays, we clarified that CAD in P. leptostachya belonged to a multigene family, and identified nearly the entire CAD gene family. Our in-depth characterization about the functions and structures of two major CAD isoforms, PlCAD2 and PlCAD3, showed that PlCAD2 exhibited the highest catalytic activity, and coniferyl aldehyde was its preferred substrate, followed by PlCAD3, and sinapyl aldehyde was its preferred substrate. Considering the accumulation patterns of furofuran lignans and expression patterns of PlCADs, we speculated that PlCAD2 was the predominant CAD isoform responsible for furofuran lignans biosynthesis in P. leptostachya. Moreover, these CADs found here can also provide effective biological parts for lignans and lignins biosynthesis.

Polymyxin B1 and E2 From Paenibacillus polymyxa Y-1 for Controlling Rice Bacterial Disease

To discover novel microbial pesticide for controlling rice bacterial disease, polymyxin B1 and E2 were firstly isolated from the supernatant of fermentation broth of Paenibacillus polymyxa Y-1 by bioactivity tracking separation. It is shown that polymyxin B1 and E2 had remarkable in vitro inhibitory activities to Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas oryzae pv. oryzicola (Xoc) with the EC50 values of 0.19 μg/ml and 0.21 μg/ml against Xoo, and 0.32 μg/ml and 0.41 μg/ml against Xoc, respectively, which were better than those of Zhongshengmycin (0.31 μg/ml and 0.73 μg/ml) and Bismerthiazol (77.48 μg/ml and 85.30 μg/ml). Polymyxins B1 and E2 had good protection and curative activities against rice bacterial leaf blight (BLB) and rice bacterial leaf streak (BLS) in vivo. The protection and curative activities of polymyxins B1 (45.8 and 35.8%, respectively) and E2 (41.2 and 37.0%, respectively) to BLB were superior to those of Zhongshengmycin (34.8 and 29.8%, respectively) and Bismerthiazol (38.0 and 33.5%, respectively). Meanwhile, the protection and curative activities of polymyxins B1 (44.8 and 39.8%, respectively) and E2 (42.9 and 39.9%, respectively) to BLS were also superior to those of Zhongshengmycin (39.7 and 32.0%, respectively) and Bismerthiazol (41.5 and 34.3%, respectively). Polymyxin B1 exerted the anti-pesticide properties via destroying the cell integrity of Xoo, reducing its infectivity and enhancing rice resistance against pathogens through activating the phenylpropanoid biosynthesis pathway of rice. It is indicated that polymyxin B1 and E2 were potential microbial pesticides for controlling rice bacterial disease.

Genome-Wide Identification and Characterization of GhCOMT Gene Family during Fiber Development and Verticillium Wilt Resistance in Cotton

Caffeic acid O-methyltransferases (COMTs) play an essential role in lignin synthesis procession, especially in the plant’s phenylalanine metabolic pathway. The content of COMT genes in cotton and the relationship between their expression patterns have not been studied clearly in cotton. In this study, we have identified 190 COMT genes in cotton, which were classified into three groups (I, II and III), and mapped on the cotton chromosomes. In addition, we found that 135 of the 190 COMT genes result from dispersed duplication (DSD) and whole-genome duplication (WGD), indicating that DSD and WGD were the main forces driving COMT gene expansion. The Ka/Ks analysis showed that GhCOMT43 and GhCOMT41 evolved from GaCOMT27 and GrCOMT14 through positive selection. The results of qRT-PCR showed that GhCOMT13, GhCOMT28, GhCOMT39 and GhCOMT55 were related to lignin content during the cotton fiber development. GhCOMT28, GhCOMT39, GhCOMT55, GhCOMT56 and GhCOMT57 responded to Verticillium Wilt (VW) and maybe related to VW resistance through lignin synthesis. Conclusively, this study found that GhCOMTs were highly expressed in the secondary wall thickening stage and VW. These results provide a clue for studying the functions of GhCOMTs in the development of cotton fiber and VW resistance and could lay a foundation for breeding cotton cultivates with higher quantity and high resistance to VW.

Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes

BACKGROUND

Dirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton.

RESULTS

In this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls.

CONCLUSIONS

These findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

Multi-responses of O-methyltransferase genes to salt stress and fiber development of Gossypium species

BACKGROUND

O-methyltransferases (OMTs) are an important group of enzymes that catalyze the transfer of a methyl group from S-adenosyl-L-methionine to their acceptor substrates. OMTs are divided into several groups according to their structural features. In Gossypium species, they are involved in phenolics and flavonoid pathways. Phenolics defend the cellulose fiber from dreadful external conditions of biotic and abiotic stresses, promoting strength and growth of plant cell wall.

RESULTS

An OMT gene family, containing a total of 192 members, has been identified and characterized in three main Gossypium species, G. hirsutum, G. arboreum and G. raimondii. Cis-regulatory elements analysis suggested important roles of OMT genes in growth, development, and defense against stresses. Transcriptome data of different fiber developmental stages in Chromosome Substitution Segment Lines (CSSLs), Recombination Inbred Lines (RILs) with excellent fiber quality, and standard genetic cotton cultivar TM-1 demonstrate that up-regulation of OMT genes at different fiber developmental stages, and abiotic stress treatments have some significant correlations with fiber quality formation, and with salt stress response. Quantitative RT-PCR results revealed that GhOMT10_Dt and GhOMT70_At genes had a specific expression in response to salt stress while GhOMT49_At, GhOMT49_Dt, and GhOMT48_At in fiber elongation and secondary cell wall stages.

CONCLUSIONS

Our results indicate that O-methyltransferase genes have multi-responses to salt stress and fiber development in Gossypium species and that they may contribute to salt tolerance or fiber quality formation in Gossypium.

Bacillus pumilus LZP02 Promotes Rice Root Growth by Improving Carbohydrate Metabolism and Phenylpropanoid Biosynthesis

Elucidation of the underlying mechanisms of plant growth promotion of rhizobacteria is very important. This study explored the mechanism by which Bacillus pumilus LZP02 promotes growth in rice roots through proteomic, transcriptomic, and metabolomic techniques. The results showed that B. pumilus LZP02 promoted the absorption of phosphorous, calcium, and magnesium ions by colonization of rice roots and enhanced peroxidase, catalase, superoxide dismutase, and Ca²⁺Mg²⁺ adenosine triphosphatase activities and chlorophyll contents in rice. The proteomic results showed that most of the differentially expressed proteins were involved in carbohydrate metabolism and that the biosynthesis of other secondary metabolites was also increased. According to RNA-seq and reverse transcription-quantitative PCR analyses, expression of some genes involved in carbohydrate metabolism and phenylpropanoid biosynthesis was upregulated in rice roots. Regarding metabolomics, phenylpropanoid biosynthesis, starch and sucrose metabolism, the pentose phosphate pathway, and glyoxylate and dicarboxylate metabolism were increased. The results indicated that B. pumilus LZP02 promoted the growth of rice roots by enhancing carbohydrate metabolism and phenylpropanoid biosynthesis.

Proteomic profiling of cotton fiber developmental transition from cell elongation to secondary wall deposition

Cotton fiber developmental transition from elongation to secondary cell wall biosynthesis is a critical growth shifting phase that affects fiber final length, strength, and other properties. Morphological dynamic analysis indicated an asynchronous fiber developmental pattern between two most important commercial cotton species, Gossypium hirsutum (Gh) and G. barbadense (Gb). Using isobaric tags for relative and absolute quantitation techniques, we examined the temporal changes of protein expression at three representative development periods (15-19, 19-23, and 23-27 dpa) in both species. Strikingly, a large proportion of differentially expressed proteins (DEPs) were identified at 19-23 dpa in Gh and at 23-27 dpa in Gb, corresponding to their fiber developmental transition timing. To better understand fiber transitional development, we comparatively analyzed those DEPs in 19-23 dpa of Gh vs. in 23-27 dpa of Gb, and noted that these cotton species indeed share fundamentally similar fiber developmental features under the biological processes. We also showed that there are limited overlaps in both specific upregulated and downregulated proteins between the two species, suggesting species-specific protein regulations in the development process. Proteomic profiling results revealed dynamic changes of several key proteins and biological processes that are potentially correlated with fiber developmental transition. During the transition, upregulated proteins are mainly involved in carbohydrate/energy metabolism, oxidation-reduction, cytoskeleton, protein turnover, Ca2+ signaling, etc., whereas important downregulated proteins are mostly involved in phenylpropanoid and flavonoid secondary metabolism pathways. The gene expressions of several changed proteins in this key stage were also examined by quantitative reverse transcription polymerase chain reaction. Overall, the present study provides accurate pictures of the regulatory networks of functional proteins during the fiber developmental transition, therefore highlighting candidate genes/proteins and related pathways for the cotton fiber improvement.

GhbHLH18 negatively regulates fiber strength and length by enhancing lignin biosynthesis in cotton fibers

Cotton fibers are developed epidermal cells of the seed coat and contain large amounts of cellulose and minor lignin-like components. Lignin in the cell walls of cotton fibers effectively provides mechanical strength and is also presumed to restrict fiber elongation and secondary cell wall synthesis. To analyze the effect of lignin and lignin-like phenolics on fiber quality and the transcriptional regulation of lignin synthesis in cotton fibers, we characterized the function of a bHLH transcription factor, GhbHLH18, during fiber elongation stage. GhbHLH18 knock-down plants have longer and stronger fibers, and accumulate less lignin-like phenolics in mature cotton fibers than control plants. By mining public transcriptomic data for developing fibers, we discovered that GhbHLH18 is coexpressed with most lignin synthesis pathway genes. Furthermore, we showed that GhbHLH18 strongly binds to the E-box in the promoter region of GhPER8 and activates its expression. Transient over expression of GhPER8 protein in tobacco leaves significantly decreased the content of coniferyl alcohol and sinapic alcohol-the substrate respectively for G-lignin and S-lignin biosynthesis. These results suggest that GhbHLH18 is negatively associated with fiber quality by activating peroxidase-mediated lignin metabolism, thus the paper represents an alternative strategy to improve fiber quality.

Genome-Wide Identification of R2R3-MYB Transcription Factors Regulating Secondary Cell Wall Thickening in Cotton Fiber Development

MYB proteins represent one of the largest transcription factor (TF) families in plants, some of which act as key transcriptional regulators of secondary cell wall (SCW) biosynthesis. Cotton (Gossypium hirsutum) fiber is thought to be an ideal single-cell model to study cell elongation and SCW biosynthesis. However, little knowledge regarding the TFs controlling fiber SCW biosynthesis, particularly for R2R3-MYBs is known. By far, no comprehensive genome-wide analysis of the secondary wall-associated R2R3-MYBs has been reported in cultivated tetraploid upland cotton. In this study, we identified 419 R2R3-MYB genes by systematically examining the cotton genome. A combination of phylogenetic, RNA-seq and co-expression analyses indicated that 36 R2R3-MYBs were either preferentially or highly expressed in 20 day post anthesis (dpa) fibers and are putative SCW regulators. Among these MYB genes, 22 MYBs are homologs of known SCW MYB proteins and the other 14 MYBs are novel proteins without prior reported SCW biosynthesis-related functions. Finally, we highlighted on the roles of two MYBs named GhMYB46_D13 and GhMYB46_D9, both of which displayed the highest expression in 20 dpa fibers. Expression of GhMYB46_D13 or GhMYB46_D9 individually in Arabidopsis resulted in ectopic SCW deposition in transgenic plants. Furthermore, both GhMYB46_D13 and GhMYB46_D9 were able to activate the cotton fiber SCW cellulose synthase gene promoters. Thus, we have identified 36 R2R3-MYBs as potential SCW regulators in cotton fibers that represent strong candidates for further functional studies during fiber development and SCW thickening.

GhHUB2, a ubiquitin ligase, is involved in cotton fiber development via the ubiquitin-26S proteasome pathway

Cotton fibers, which are extremely elongated single cells of epidermal seed trichomes and have highly thickened cell walls, constitute the most important natural textile material worldwide. However, the regulation of fiber development is not well understood. Here, we report that GhHUB2, a functional homolog of AtHUB2, controls fiber elongation and secondary cell wall (SCW) deposition. GhHUB2 is ubiquitously expressed, including within fibers. Overexpression of GhHUB2 in cotton increased fiber length and SCW thickness, while RNAi knockdown of GhHUB2 resulted in shortened fibers and thinner cell walls. We found that GhHUB2 interacted with GhKNL1, a transcriptional repressor predominantly expressed in developing fibers, and that GhHUB2 ubiquitinated and degraded GhKNL1 via the ubiquitin-26S proteasome pathway. GhHUB2 negatively regulated GhKNL1 protein levels and lead to the disinhibition of genes such as GhXTH1, Gh1,3-β-G, GhCesA4, GhAGP4, GhCTL1, and GhCOBL4, thus promoting fiber elongation and enhancing SCW biosynthesis. We found that GhREV-08, a transcription factor that participates in SCW deposition and auxin signaling pathway, was a direct target of GhKNL1. In conclusion, our study uncovers a novel function of HUB2 in plants in addition to its monoubiquitination of H2B. Moreover, we provide evidence for control of the fiber development by the ubiquitin-26S proteasome pathway.

The phosphatidylinositol synthase gene (GhPIS) contributes to longer, stronger, and finer fibers in cotton

Cotton fibers are the most important natural raw material used in textile industries world-wide. Fiber length, strength, and fineness are the three major traits which determine the quality and economic value of cotton. It is known that exogenous application of phosphatidylinositols (PtdIns), important structural phospholipids, can promote cotton fiber elongation. Here, we sought to increase the in planta production of PtdIns to improve fiber traits. Transgenic cotton plants were generated in which the expression of a cotton phosphatidylinositol synthase gene (i.e., GhPIS) was controlled by the fiber-specific SCFP promoter element, resulting in the specific up-regulation of GhPIS during cotton fiber development. We demonstrate that PtdIns content was significantly enhanced in transgenic cotton fibers and the elevated level of PtdIns stimulated the expression of genes involved in PtdIns phosphorylation as well as promoting lignin/lignin-like phenolic biosynthesis. Fiber length, strength and fineness were also improved in the transgenic plants as compared to the wild-type cotton, with no loss in overall fiber yield. Our data indicate that fiber-specific up-regulation of PtdIns synthesis is a promising strategy for cotton fiber quality improvement.

GhMYB1 regulates SCW stage-specific expression of the GhGDSL promoter in the fibres of Gossypium hirsutum L

Secondary cell wall (SCW) biosynthesis is an important stage of the cotton fibre development, and its transcriptional regulation is poorly understood. We selected the Gossypium hirsutum GDSL (GhGDSL) lipase/hydrolase gene (CotAD_74480), which is expressed during SCW biosynthesis (19 through to 25 days postanthesis; DPA), for study. T₁ -transgenic cotton lines expressing the β-glucuronidase (gus) reporter under the control of a 1026-bp promoter fragment of GhGDSL (P_GhGDSL ) showed 19 DPA stage-specific increase in GUS expression. 5' deletion indicated that the 194-bp fragment between -788 and -594 relative to the transcription start site was essential for this stage-specific expression. Site-directed mutagenesis of eight transcription factor binding sites within P_GhGDSL demonstrated that the MYB1AT motif (AAACCA) at -603/-598 was critical for the 19 DPA-specific reporter gene expressions. Yeast one-hybrid (Y1H) analysis identified nine proteins, including GhMYB1 (CotAD_64719) that bound to the P_GhGDSL promoter. Further, Y1H experiments using the 5' promoter deletions and individually mutated promoter motifs indicated that GhMYB1 interacted with P_GhGDSL at MYB1AT sequence. GhMYB1 was expressed specifically in fibre from 19 DPA, overlapping with the sharp rise in GhGDSL expression, indicating that it could regulate GhGDSL during fibre development. Analysis of genes co-expressed with GhMYB1 showed that it potentially regulates a number of other 19-25 DPA-specific genes in networks including those functioning in the cell wall and precursor synthesis, but not the major polysaccharide and protein components of the fibre SCW. GhGDSL and its promoter are therefore potential tools for the improvement of cotton fibre quality traits.

Tissue and cell-specific transcriptomes in cotton reveal the subtleties of gene regulation underlying the diversity of plant secondary cell walls

BACKGROUND

Knowledge of plant secondary cell wall (SCW) regulation and deposition is mainly based on the Arabidopsis model of a 'typical' lignocellulosic SCW. However, SCWs in other plants can vary from this. The SCW of mature cotton seed fibres is highly cellulosic and lacks lignification whereas xylem SCWs are lignocellulosic. We used cotton as a model to study different SCWs and the expression of the genes involved in their formation via RNA deep sequencing and chemical analysis of stem and seed fibre.

RESULTS

Transcriptome comparisons from cotton xylem and pith as well as from a developmental series of seed fibres revealed tissue-specific and developmentally regulated expression of several NAC transcription factors some of which are likely to be important as top tier regulators of SCW formation in xylem and/or seed fibre. A so far undescribed hierarchy was identified between the top tier NAC transcription factors SND1-like and NST1/2 in cotton. Key SCW MYB transcription factors, homologs of Arabidopsis MYB46/83, were practically absent in cotton stem xylem. Lack of expression of other lignin-specific MYBs in seed fibre relative to xylem could account for the lack of lignin deposition in seed fibre. Expression of a MYB103 homolog correlated with temporal expression of SCW CesAs and cellulose synthesis in seed fibres. FLAs were highly expressed and may be important structural components of seed fibre SCWs. Finally, we made the unexpected observation that cell walls in the pith of cotton stems contained lignin and had a higher S:G ratio than in xylem, despite that tissue's lacking many of the gene transcripts normally associated with lignin biosynthesis.

CONCLUSIONS

Our study in cotton confirmed some features of the currently accepted gene regulatory cascade for 'typical' plant SCWs, but also revealed substantial differences, especially with key downstream NACs and MYBs. The lignocellulosic SCW of cotton xylem appears to be achieved differently from that in Arabidopsis. Pith cell walls in cotton stems are compositionally very different from that reported for other plant species, including Arabidopsis. The current definition of a 'typical' primary or secondary cell wall might not be applicable to all cell types in all plant species.

Distribution of Silicified Microstructures, Regulation of Cinnamyl Alcohol Dehydrogenase and Lodging Resistance in Silicon and Paclobutrazol Mediated Oryza sativa

Lodging is a phenomenon that affects most of the cereal crops including rice, Oryza sativa. This is due to the fragile nature of herbaceous plants whose stems are non-woody, thus affecting its ability to grow upright. Silicon (Si), a beneficial nutrient is often used to toughen and protect plants from biotic and abiotic stresses. Deposition of Si in plant tissues enhances the rigidity and stiffness of the plant as a whole. Silicified cells provide the much needed strength to the culm to resist breaking. Lignin plays important roles in cell wall structural integrity, stem strength, transport, mechanical support, and plant pathogen defense. The aim of this study is to resolve effects of Si on formation of microstructure and regulation of cinnamyl alcohol dehydrogenase (CAD), a key gene responsible for lignin biosynthesis. Besides evaluating silicon, paclobutrazol (PBZ) a plant growth retartdant that reduces internode elongation is also incorporated in this study. Hardness, brittleness and stiffness were improved in presence of silicon thus reducing lodging. Scanning electron micrographs with the aid of energy dispersive x-ray (EDX) was used to map silicon distribution. Presence of trichomes, silica cells, and silica bodies were detected in silicon treated plants. Transcripts of CAD gene was also upregulated in these plants. Besides, phloroglucinol staining showed presence of lignified vascular bundles and sclerenchyma band. In conclusion, silicon treated rice plants showed an increase in lignin content, silicon content, and formation of silicified microstructures.

Genome-wide identification of multifunctional laccase gene family in cotton (Gossypium spp.); expression and biochemical analysis during fiber development

The single-celled cotton fibers, produced from seed coat epidermal cells are the largest natural source of textile fibers. The economic value of cotton fiber lies in its length and quality. The multifunctional laccase enzymes play important roles in cell elongation, lignification and pigmentation in plants and could play crucial role in cotton fiber quality. Genome-wide analysis of cultivated allotetraploid (G. hirsutum) and its progenitor diploid (G. arboreum and G. raimondii) cotton species identified 84, 44 and 46 laccase genes, respectively. Analysis of chromosomal location, phylogeny, conserved domain and physical properties showed highly conserved nature of laccases across three cotton species. Gene expression, enzymatic activity and biochemical analysis of developing cotton fibers was performed using G. arboreum species. Of the total 44, 40 laccases showed expression during different stages of fiber development. The higher enzymatic activity of laccases correlated with higher lignin content at 25 DPA (Days Post Anthesis). Further, analysis of cotton fiber phenolic compounds showed an overall decrease at 25 DPA indicating possible incorporation of these substrates into lignin polymer during secondary cell wall biosynthesis. Overall data indicate significant roles of laccases in cotton fiber development, and presents an excellent opportunity for manipulation of fiber development and quality.

A R2R3-MYB transcription factor that is specifically expressed in cotton (Gossypium hirsutum) fibers affects secondary cell wall biosynthesis and deposition in transgenic Arabidopsis

Secondary cell wall (SCW) is an important industrial raw material for pulping, papermaking, construction, lumbering, textiles and potentially for biofuel production. The process of SCW thickening of cotton fibers lays down the cellulose that will constitute the bulk (up to 96%) of the fiber at maturity. In this study, a gene encoding a MYB-domain protein was identified in cotton (Gossypium hirsutum) and designated as GhMYBL1. Quantitative real-time polymerase chain reaction (RT-PCR) analysis revealed that GhMYBL1 was specifically expressed in cotton fibers at the stage of secondary wall deposition. Further analysis indicated that this protein is a R2R3-MYB transcription factor, and is targeted to the cell nucleus. Overexpression of GhMYBL1 in Arabidopsis affected the formation of SCW in the stem xylem of the transgenic plants. The enhanced SCW thickening also occurred in the interfascicular fibers, xylary fibers and vessels of the GhMYBL1-overexpression transgenic plants. The expression of secondary wall-associated genes, such as CesA4, CesA7, CesA8, PAL1, F5H and 4CL1, were upregulated, and consequently, cellulose and lignin biosynthesis were enhanced in the GhMYBL1 transgenic plants. These data suggested that GhMYBL1 may participate in modulating the process of secondary wall biosynthesis and deposition of cotton fibers.

Metabolomic and transcriptomic insights into how cotton fiber transitions to secondary wall synthesis, represses lignification, and prolongs elongation

Background

The morphogenesis of single-celled cotton fiber includes extreme elongation and staged cell wall differentiation. Designing strategies for improving cotton fiber for textiles and other uses relies on uncovering the related regulatory mechanisms. In this research we compared the transcriptomes and metabolomes of two Gossypium genotypes, Gossypium barbadense cv Phytogen 800 and G. hirsutum cv Deltapine 90. When grown in parallel, the two types of fiber developed similarly except for prolonged fiber elongation in the G. barbadense cultivar. The data were collected from isolated fibers between 10 to 28 days post anthesis (DPA) representing: primary wall synthesis to support elongation; transitional cell wall remodeling; and secondary wall cellulose synthesis, which was accompanied by continuing elongation only in G. barbadense fiber.

Results

Of 206 identified fiber metabolites, 205 were held in common between the two genotypes. Approximately 38,000 transcripts were expressed in the fiber of each genotype, and these were mapped to the reference set and interpreted by homology to known genes. The developmental changes in the transcriptomes and the metabolomes were compared within and across genotypes with several novel implications. Transitional cell wall remodeling is a distinct stable developmental stage lasting at least four days (18 to 21 DPA). Expression of selected cell wall related transcripts was similar between genotypes, but cellulose synthase gene expression patterns were more complex than expected. Lignification was transcriptionally repressed in both genotypes. Oxidative stress was lower in the fiber of G. barbadense cv Phytogen 800 as compared to G. hirsutum cv Deltapine 90. Correspondingly, the G. barbadense cultivar had enhanced capacity for management of reactive oxygen species during its prolonged elongation period, as indicated by a 138-fold increase in ascorbate concentration at 28 DPA.

Conclusions

The parallel data on deep-sequencing transcriptomics and non-targeted metabolomics for two genotypes of single-celled cotton fiber showed that a discrete developmental stage of transitional cell wall remodeling occurs before secondary wall cellulose synthesis begins. The data showed how lignification can be transcriptionally repressed during secondary cell wall synthesis, and they implicated enhanced capacity to manage reactive oxygen species through the ascorbate-glutathione cycle as a positive contributor to fiber length.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1708-9) contains supplementary material, which is available to authorized users.

Proteomics profiling of fiber development and domestication in upland cotton (Gossypium hirsutum L.)

Comparative proteomic analyses were performed to detail the evolutionary consequences of strong directional selection for enhanced fiber traits in modern upland cotton (Gossypium hirsutum L.). Using two complementary proteomic approaches, 2-DE and iTRAQ LC-MS/MS, fiber proteomes were examined for four representative stages of fiber development. Approximately 1,000 protein features were characterized using each strategy, collectively resulting in the identification and functional categorization of 1,223 proteins. Unequal contributions of homoeologous proteins were detected for over a third of the fiber proteome, but overall expression was balanced with respect to the genome-of-origin in the allopolyploid G. hirsutum. About 30% of the proteins were differentially expressed during fiber development within wild and domesticated cotton. Notably, domestication was accompanied by a doubling of protein developmental dynamics for the period between 10 and 20 days following pollination. Expression levels of 240 iTRAQ proteins and 293 2-DE spots were altered by domestication, collectively representing multiple cellular and metabolic processes, including metabolism, energy, protein synthesis and destination, defense and stress response. Analyses of homoeolog-specific expression indicate that duplicated gene products in cotton fibers can be differently regulated in response to selection. These results demonstrate the power of proteomics for the analysis of crop domestication and phenotypic evolution.

Molecular cloning and functional analysis of nine cinnamyl alcohol dehydrogenase family members in Populus tomentosa

Nine CAD/CAD-like genes in P. tomentosa were classified into four classes based on expression patterns, phylogenetic analysis and biochemical properties with modification for the previous claim of SAD. Cinnamyl alcohol dehydrogenase (CAD) functions in monolignol biosynthesis and plays a critical role in wood development and defense. In this study, we isolated and cloned nine CAD/CAD-like genes in the Populus tomentosa genome. We investigated differential expression using microarray chips and found that PtoCAD1 was highly expressed in bud, root and vascular tissues (xylem and phloem) with the greatest expression in the root. Differential expression in tissues was demonstrated for PtoCAD3, PtoCAD6 and PtoCAD9. Biochemical analysis of purified PtoCADs in vitro indicated PtoCAD1, PtoCAD2 and PtoCAD8 had detectable activity against both coniferaldehyde and sinapaldehyde. PtoCAD1 used both substrates with high efficiency. PtoCAD2 showed no specific requirement for sinapaldehyde in spite of its high identity with so-called PtrSAD (sinapyl alcohol dehydrogenase). In addition, the enzymatic activity of PtoCAD1 and PtoCAD2 was affected by temperature. We classified these nine CAD/CAD-like genes into four classes: class I included PtoCAD1, which was a bone fide CAD with the highest activity; class II included PtoCAD2, -5, -7, -8, which might function in monolignol biosynthesis and defense; class III genes included PtoCAD3, -6, -9, which have a distinct expression pattern; class IV included PtoCAD12, which has a distinct structure. These data suggest divergence of the PtoCADs and its homologs, related to their functions. We propose genes in class II are a subset of CAD genes that evolved before angiosperms appeared. These results suggest CAD/CAD-like genes in classes I and II play a role in monolignol biosynthesis and contribute to our knowledge of lignin biosynthesis in P. tomentosa.

Fourier transform mid infrared spectroscopy applications for monitoring the structural plasticity of plant cell walls

Fourier transform mid-infrared (FT-MIR) spectroscopy has been extensively used as a potent, fast and non-destructive procedure for analyzing cell wall architectures, with the capacity to provide abundant information about their polymers, functional groups, and in muro entanglement. In conjunction with multivariate analyses, this method has proved to be a valuable tool for tracking alterations in cell walls. The present review examines recent progress in the use of FT-MIR spectroscopy to monitor cell wall changes occurring in muro as a result of various factors, such as growth and development processes, genetic modifications, exposition or habituation to cellulose biosynthesis inhibitors and responses to other abiotic or biotic stresses, as well as its biotechnological applications.

Proteomic profiling of developing cotton fibers from wild and domesticated Gossypium barbadense

Pima cotton (Gossypium barbadense) is widely cultivated because of its long, strong seed trichomes ('fibers') used for premium textiles. These agronomically advanced fibers were derived following domestication and thousands of years of human-mediated crop improvement. To gain an insight into fiber development and evolution, we conducted comparative proteomic and transcriptomic profiling of developing fiber from an elite cultivar and a wild accession. Analyses using isobaric tag for relative and absolute quantification (iTRAQ) LC-MS/MS technology identified 1317 proteins in fiber. Of these, 205 were differentially expressed across developmental stages, and 190 showed differential expression between wild and cultivated forms, 14.4% of the proteome sampled. Human selection may have shifted the timing of developmental modules, such that some occur earlier in domesticated than in wild cotton. A novel approach was used to detect possible biased expression of homoeologous copies of proteins. Results indicate a significant partitioning of duplicate gene expression at the protein level, but an approximately equal degree of bias for each of the two constituent genomes of allopolyploid cotton. Our results demonstrate the power of complementary transcriptomic and proteomic approaches for the study of the domestication process. They also provide a rich database for mining for functional analyses of cotton improvement or evolution.

Molecular cloning, functional analysis of three cinnamyl alcohol dehydrogenase (CAD) genes in the leaves of tea plant, Camellia sinensis

Cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.195) is considered to be a key enzyme in lignin biosynthesis, but little was known about CADs in tea plants (Camellia sinensis). A full-length cDNA sequence (CsCAD2) was isolated by suppressive subtractive hybridization (SSH) in Ectropis oblique feeding-induced tea plants, and another two full-length cDNA sequences (CsCAD1 and CsCAD3) were obtained from a transcriptome obtained by deep sequencing. However, they showed only 20-54% identities. Phylogenetic analysis revealed that they belonged to three different families. DNA gel blotting analysis revealed that two copies of CsCAD1 and CsCAD2 genes existed in tea genome, but CsCAD3 likely had only one copy. Recombinant proteins of these CsCADs were produced in Escherichia coli. The activity of purified recombinant CsCAD2 protein was up to 0.43 μmol min(-1) mg(-1). However, the other two recombinant proteins had lower activities, probably due to incomplete refolding. qRT-PCR analysis indicated that while CsCAD3 was strongly up-regulated in tea plants after E. oblique attack and mechanical damage, CsCAD1 and CsCAD2 showed only moderate or no changes in transcript levels. Treatment of defence-related hormones methyl jasmonate (MeJA) and salicylic acid (SA) elevated the expression of CsCAD1 and CsCAD2, but decreased the transcript abundance of CsCAD3. The transcript levels of CsCAD2 did not change after applying abscisic acid (ABA), whereas CsCAD1 and CsCAD3 were induced. These results suggested that these three CsCAD genes in tea plants may play a role in defense against insects and pathogens and adaptation to abiotic stresses and these genes likely have divergant functions.

Genetic mapping and characteristics of genes specifically or preferentially expressed during fiber development in cotton

Cotton fiber is an ideal model to study cell elongation and cell wall construction in plants. During fiber development, some genes and proteins have been reported to be specifically or preferentially expressed. Mapping of them will reveal the genomic distribution of these genes, and will facilitate selection in cotton breeding. Based on previous reports, we designed 331 gene primers and 164 protein primers, and used single-strand conformation polymorphism (SSCP) to map and integrate them into our interspecific BC₁ linkage map. This resulted in the mapping of 57 loci representing 51 genes or proteins on 22 chromosomes. For those three markers which were tightly linked with quantitative trait loci (QTLs), the QTL functions obtained in this study and gene functions reported in previous reports were consistent. Reverse transcription-polymerase chain reaction (RT-PCR) analysis of 52 polymorphic functional primers showed that 21 gene primers and 17 protein primers had differential expression between Emian22 (Gossypium hirsutum) and 3–79 (G. barbadense). Both RT-PCR and quantitative real-time PCR (qRT-PCR) analyses of the three markers tightly linked with QTLs were consistent with QTL analysis and field experiments. Gene Ontology (GO) categorization revealed that almost all 51 mapped genes belonged to multiple categories that contribute to fiber development, indicating that fiber development is a complex process regulated by various genes. These 51 genes were all specifically or preferentially expressed during fiber cell elongation and secondary wall biosynthesis. Therefore, these functional gene-related markers would be beneficial for the genetic improvement of cotton fiber length and strength.

Expression patterns of a cinnamyl alcohol dehydrogenase gene involved in lignin biosynthesis and environmental stress in Ginkgo biloba

The cinnamyl alcohol dehydrogenase (CAD) is a key enzyme in lignin biosynthesis as it catalyzes the final step in the synthesis of monolignols. A cDNA sequence encoding the CAD gene was isolated from the leaves of Ginkgo biloba L, designated as GbCAD1. The full-length cDNA of GbCAD1 was 1,494 bp containing a 1,074 bp open reading frame encoding a polypeptide of 357 amino acids with a calculated molecular mass of 38.7 kDa and an isoelectric point of 5.74. Comparative and bioinformatic analyses revealed that GbCAD1 showed extensive homology with CADs from other gymnosperm species. Southern blot analysis indicated that GbCAD1 belonged to a multi-gene family. Phylogenetic tree analysis revealed that GbCAD1 shared the same ancestor in evolution with other CADs and had a further relationship with other gymnosperm species. GbCAD1 was an enzyme being pH-dependent and temperature-sensitive, and showing a selected catalyzing. Tissue expression pattern analysis showed that GbCAD1 was constitutively expressed in stems and roots, especially in the parts of the pest and disease infection, with the lower expression being found in two- to four-year-old stem. Further analysis showed the change in lignin content had some linear correlation with the expression level of GbCAD1 mRNA in different tissues. The increased expression of GbCAD1 was detected when the seedling were treated with exogenous abscisic acid, salicylic acid, ethephon, ultraviolet and wounding. These results indicate that the GbCAD1 gene may play a role in the resistance mechanism to biotic and abiotic stresses as well as in tissue-specific developmental lignification.

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

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

Abstract Through exploring potential analogies between cotton seed trichomes (or cotton fiber) and arabidopsis shoot trichomes we discovered that CesAs from either the primary or secondary wall phylogenetic clades can support secondary wall thickening. CesA genes that typically support primary wall synthesis, AtCesA1,2,3,5, and 6, underpin expansion and secondary wall thickening of arabidopsis shoot trichomes. In contrast, apparent orthologs of CesA genes that support secondary wall synthesis in arabidopsis xylem, AtCesA4,7, and 8, are up-regulated for cotton fiber secondary wall deposition. These conclusions arose from: (a) analyzing the expression of CesA genes in arabidopsis shoot trichomes; (b) observing birefringent secondary walls in arabidopsis shoot trichomes with mutations in AtCesA4, 7, or 8; (c) assaying up-regulated genes during different stages of cotton fiber development; and (d) comparing genes that were co-expressed with primary or secondary wall CesAs in arabidopsis with genes up-regulated in arabidopsis trichomes, arabidopsis secondary xylem, or cotton fiber during primary or secondary wall deposition. Cumulatively, the data show that: (a) the xylem of arabidopsis provides the best model for secondary wall cellulose synthesis in cotton fiber; and (b) CesA genes within a "cell wall toolbox" are used in diverse ways for the construction of particular specialized cell walls.