Two cotton fiber-associated glycosyltransferases, GhGT43A1 and GhGT43C1, function in hemicellulose glucuronoxylan biosynthesis during plant development

Metabolite Profiling and Biological Activity Assessment of Paeonia ostii Anthers and Pollen Using UPLC-QTOF-MS

Paeonia ostii is an important economic oil and medicinal crop. Its anthers are often used to make tea in China with beneficial effects on human health. However, the metabolite profiles, as well as potential biological activities of P. ostii anthers and the pollen within anthers have not been systematically analyzed, which hinders the improvement of P. ostii utilization. With comprehensive untargeted metabolomic analysis using UPLC-QTOF-MS, we identified a total of 105 metabolites in anthers and pollen, mainly including phenylpropanoids, polyketides, organic acids, benzenoids, lipids, and organic oxygen compounds. Multivariate statistical analysis revealed the metabolite differences between anthers and pollen, with higher carbohydrates and flavonoids content in pollen and higher phenolic content in anthers. Meanwhile, both anthers and pollen extracts exhibited antioxidant activity, antibacterial activity, α-glucosidase and α-amylase inhibitory activity. In general, the anther stage of S4 showed the highest biological activity among all samples. This study illuminated the metabolites and biological activities of anthers and pollen of P. ostii, which supports the further utilization of them.

A wheat heat shock transcription factor gene, TaHsf-7A, regulates seed dormancy and germination

Heat shock transcription factors (Hsfs) play multifaceted roles in plant growth, development, and responses to environmental factors. However, their involvement in seed dormancy and germination processes has remained elusive. In this study, we identified a wheat class B Hsf gene, TaHsf-7A, with higher expression in strong-dormancy varieties compared to weak-dormancy varieties during seed imbibition. Specifically, TaHsf-7A expression increased during seed dormancy establishment and subsequently declined during dormancy release. Through the identification of a 1-bp insertion (ins)/deletion (del) variation in the coding region of TaHsf-7A among wheat varieties with different dormancy levels, we developed a CAPS marker, Hsf-7A-1319, resulting in two allelic variations: Hsf-7A-1319-ins and Hsf-7A-1319-del. Notably, the allele Hsf-7A-1319-ins correlated with a reduced seed germination rate and elevated dormancy levels, while Hsf-7A-1319-del exhibited the opposite trend across 175 wheat varieties. The association of TaHsf-7A allelic status with seed dormancy and germination levels was confirmed in various genetically modified species, including Arabidopsis, rice, and wheat. Results from the dual luciferase assay demonstrated notable variations in transcriptional activity among transformants harboring distinct TaHsf-7A alleles. Furthermore, the levels of abscisic acid (ABA) and gibberellin (GA), along with the expression levels of ABA and GA biosynthesis genes, showed significant differences between transgenic rice lines carrying different alleles of TaHsf-7A. These findings represent a significant step towards a comprehensive understanding of TaHsf-7A's involvement in the dormancy and germination processes of wheat seeds.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Transcriptome and proteome profiling revealed the key genes and pathways involved in the fiber quality formation in brown cotton

Colored cotton is also called eco-cotton because of its natural color fiber. It is inferior in yield and quality than white cotton. The underlying regulatory genes involved in fiber quality and pigment synthesis are not well understood. This study aimed to investigate the transcriptomic and proteomic changes during fiber development in a brown cotton cultivar (Z161) and a white cotton cultivar. The differential proteins with the same expression trend as genes were significantly and positively correlated with corresponding fold changes in expression. Enrichment analysis revealed that Z161, enriched in fiber elongation genes related to flavonoid biosynthesis, phenylalanine metabolism, glutathione metabolism, and many more genes (proteins) are up-regulated. Moreover, 164 glycosyltransferases genes, 15 MYB-bHLH-WD40 genes, and other transcription factors such as C2H2 (12), ERF (11), and NAC (7) were preferentially expressed in Z161. Weighted correlation network analysis identified fatty acid synthesis and energy metabolism as the principal metabolic pathways in both cotton genotypes during fiber development. Identified 15 hub genes will provide important insights for genetic manipulation of fiber quality and pigment deposition balance in brown cotton fibers.

Multifaceted analyses reveal carbohydrate metabolism mainly affecting the quality of postharvest bamboo shoots

Bamboo shoot is one of nutritious vegetables in China. However, the edible quality of fresh bamboo shoots deteriorates easily after harvest. Here, morphological, physiological, transcriptomic and microRNA sequencing analyses were conducted to investigate the postharvest characteristics of moso bamboo (Phyllostachys edulis) shoots. Rapid decreases of soluble sugars, structural polysaccharides and hydrolyzed tannins, and increases of lignin and condensed tannins were observed in the postharvest bamboo shoots. Differentially expressed genes (DEGs) and miRNAs with opposite trends were mainly enriched in structural polysaccharide metabolism, starch and sucrose metabolism and glycolysis pathways, which were consistent with the changes of carbohydrates. A co-expression network of carbohydrate metabolism was constructed, which was verified by qPCR and yeast one-hybrid (Y1H) assay. Furthermore, the function of one hub glycosyltransferase gene was validated in Arabidopsis, which confirmed that it was involved in xylan biosynthesis. These results are of great significance for revealing the carbohydrate metabolism mechanisms of postharvest bamboo shoots and provide a potential candidate gene for molecular breeding related to xylan in the future.

The DUF288 domain containing proteins GhSTLs participate in cotton fiber cellulose synthesis and impact on fiber elongation

Cotton is one of the most important economic crops in the world, with over 90 % cellulose in the mature fiber. However, the cellulose synthesis mechanism in cotton fibers is poorly understood. Here, we identified four DUF288 domain containing proteins, which we designated GhSTL1-4. These four GhSTL genes are highly expressed in 6 days post anthesis (dpa) and 20 dpa cotton fibers. They are localized to the Golgi apparatus, and can rescue the growth defects in primary cell wall (PCW) and secondary cell wall (SCW) of cellulose synthesis of the Arabidopsis stl1stl2 double mutant at varying degrees. Silencing of GhSTLs resulted in reduced cellulose content and shorter fibers. In addition, split-ubiquitin membrane yeast two-hybrid analysis showed that GhSTL1 and GhSTL4 can interact with PCW-related GhCesA6-1/6-3 and SCW-associated GhCesA7-1/7-2. GhSTL3 can interact with SCW-related GhCesA4-3. These interactions are further confirmed by firefly luciferase complementation imaging assay. Together, we demonstrate that GhSTLs can selectively interact with both the PCW and SCW-associated GhCesAs and impact on cellulose synthesis and fiber development. Our findings provide insights into the mechanism underlying cellulose biosynthesis in cotton fibers, and offer potential candidate genes to coordinate PCW and SCW cellulose synthesis of cotton fibers for developing elite cotton varieties with enhanced fiber quality.

GhMYB7 promotes secondary wall cellulose deposition in cotton fibres by regulating GhCesA gene expression through three distinct cis-elements

Cotton fibre is the most important source for natural textiles. The secondary cell walls (SCWs) of mature cotton fibres contain the highest proportion of cellulose content (> 90%) in any plant. The onset and progression of SCW cellulose synthesis need to be tightly controlled to balance fibre elongation and cell wall deposition. However, regulatory mechanisms that control cellulose synthesis during cotton fibre growth remain elusive. Here, we conducted genetic and functional analyses demonstrating that the R2R3-MYB GhMYB7 controls cotton fibre cellulose synthesis. Overexpression of GhMYB7 in cotton sped up SCW cellulose biosynthesis in fibre cells, and led to shorter fibres with thicker walls. By contrast, RNA interference (RNAi) silencing of GhMYB7 delayed fibre SCW cellulose synthesis and resulted in elongated fibres with thinner walls. Furthermore, we demonstrated that GhMYB7 regulated cotton fibre SCW cellulose synthases by directly binding to three distinct cis-elements in the respective GhCesA4, GhCesA7 and GhCesA8 promoters. We found that this regulatory mechanism of cellulose synthesis was 'hi-jacked' also by other GhMYBs. Together, our findings uncover a hitherto-unknown mechanism that cotton fibre employs to regulate SCW cellulose synthesis. Our results also provide a strategy for genetic improvement of SCW thickness of cotton fibre.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes Contributing to Fuzz Development in Gossypium arboreum

Fuzzless mutants are ideal materials to decipher the regulatory network and mechanism underlying fuzz initiation and formation. In this study, we utilized two Gossypium arboreum accessions differing in fuzz characteristics to explore expression pattern differences and discriminate genes involved in fuzz development using RNA sequencing. Gene ontology (GO) analysis was conducted and found that DEGs were mainly enriched in the regulation of transcription, metabolic processes and oxidation–reduction-related processes. Weighted gene co-expression network analysis discerned the MEmagenta module highly associated with a fuzz/fuzzless trait, which included a total of 50 hub genes differentially expressed between two materials. GaFZ, which negatively regulates trichome and fuzz formation, was found involved in MEmagenta cluster1. In addition, twenty-eight hub genes in MEmagenta cluster1 were significantly up-regulated and expressed in fuzzless mutant DPL972. It is noteworthy that Ga04G1219 and Ga04G1240, which, respectively, encode Fasciclin-like arabinogalactan protein 18(FLA18) and transport protein, showed remarkable differences of expression level and implied that they may be involved in protein glycosylation to regulate fuzz formation and development. This module and hub genes identified in this study will provide new insights on fiber and fuzz formation and be useful for the molecular design breeding of cotton genetic improvement.

A cotton NAC domain transcription factor, GhFSN5, negatively regulates secondary cell wall biosynthesis and anther development in transgenic Arabidopsis

NAC domain transcription factors (TFs) are plant-specific transcriptional regulators, some of which play crucial roles in secondary cell wall (SCW) biosynthesis in plants. Cotton is one of the most important natural fiber producing crops, whose mature fiber SCW contains more than 90% cellulose with very small amounts of xylan and lignin, but little is known about the molecular mechanism underlying fiber SCW formation. We previously identified seven fiber preferentially expressed NAC members, GhFSN1-7. One, GhFSN1, was demonstrated to positively regulate fiber SCW thickening, but the functions of other GhFSN members remain unknown. In this study, roles of GhFSN5 were dissected. qRT-PCR analysis showed that GhFSN5 was predominantly transcribed during the fiber SCW thickening stage. In addition, a large number of fiber SCW biosynthetic genes and SCW-related TFs were co-expressed with GhFSN5. Heterologous expression of GhFSN5 in Arabidopsis resulted in plants with smaller siliques and severe sterility. Anther dehiscence in transgenic lines was not substantially affected, but most pollen was collapsed and nonviable. Furthermore, cellulose and lignin contents in inflorescence stems as well as roots were reduced in transgenic lines, compared with the wild type. Moreover, a set of SCW biosynthetic genes for cellulose, xylan and lignin and several transcription factors involved in regulation of SCW formation were down-regulated in transgenic plants. Our findings indicate that GhFSN5 acts as a negative regulator of SCW formation and anther development and expands our understanding of transcriptional regulation of SCW biosynthesis.

Genome-Wide Identification and Characterization of Glycosyltransferase Family 47 in Cotton

The glycosyltransferase (GT) 47 family is involved in the biosynthesis of xylose, pectin and xyloglucan and plays a significant role in maintaining the normal morphology of the plant cell wall. However, the functions of GT47s are less well known in cotton. In the present study, a total of 53, 53, 105 and 109 GT47 genes were detected by genome-wide identification in Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. All the GT47s were classified into six major groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GT47 genes was highly conserved. Collinearity analysis showed that GT47 gene family expansion occurred in Gossypium spp. mainly through whole-genome duplication and that segmental duplication mainly promoted GT47 gene expansion within the A and D subgenomes. The Ka/Ks values suggested that the GT47 gene family has undergone purifying selection during the long-term evolutionary process. Transcriptomic data and qRT-PCR showed that GhGT47 genes exhibited different expression patterns in each tissue and during fiber development. Our results suggest that some genes in the GhGT47 family might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GT47 genes in cotton.

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.

Xylan in the Middle: Understanding Xylan Biosynthesis and Its Metabolic Dependencies Toward Improving Wood Fiber for Industrial Processing

Lignocellulosic biomass, encompassing cellulose, lignin and hemicellulose in plant secondary cell walls (SCWs), is the most abundant source of renewable materials on earth. Currently, fast-growing woody dicots such as Eucalyptus and Populus trees are major lignocellulosic (wood fiber) feedstocks for bioproducts such as pulp, paper, cellulose, textiles, bioplastics and other biomaterials. Processing wood for these products entails separating the biomass into its three main components as efficiently as possible without compromising yield. Glucuronoxylan (xylan), the main hemicellulose present in the SCWs of hardwood trees carries chemical modifications that are associated with SCW composition and ultrastructure, and affect the recalcitrance of woody biomass to industrial processing. In this review we highlight the importance of xylan properties for industrial wood fiber processing and how gaining a greater understanding of xylan biosynthesis, specifically xylan modification, could yield novel biotechnology approaches to reduce recalcitrance or introduce novel processing traits. Altering xylan modification patterns has recently become a focus of plant SCW studies due to early findings that altered modification patterns can yield beneficial biomass processing traits. Additionally, it has been noted that plants with altered xylan composition display metabolic differences linked to changes in precursor usage. We explore the possibility of using systems biology and systems genetics approaches to gain insight into the coordination of SCW formation with other interdependent biological processes. Acetyl-CoA, s-adenosylmethionine and nucleotide sugars are precursors needed for xylan modification, however, the pathways which produce metabolic pools during different stages of fiber cell wall formation still have to be identified and their co-regulation during SCW formation elucidated. The crucial dependence on precursor metabolism provides an opportunity to alter xylan modification patterns through metabolic engineering of one or more of these interdependent pathways. The complexity of xylan biosynthesis and modification is currently a stumbling point, but it may provide new avenues for woody biomass engineering that are not possible for other biopolymers.

Downregulating aspen xylan biosynthetic GT43 genes in developing wood stimulates growth via reprograming of the transcriptome

Xylan is one of the main compounds determining wood properties in hardwood species. The xylan backbone is thought to be synthesized by a synthase complex comprising two members of the GT43 family. We downregulated all GT43 genes in hybrid aspen (Populus tremula × tremuloides) to understand their involvement in xylan biosynthesis. All three clades of the GT43 family were targeted for downregulation using RNA interference individually or in different combinations, either constitutively or specifically in developing wood. Simultaneous downregulation in developing wood of the B (IRX9) and C (IRX14) clades resulted in reduced xylan Xyl content relative to reducing end sequence, supporting their role in xylan backbone biosynthesis. This was accompanied by a higher lignocellulose saccharification efficiency. Unexpectedly, GT43 suppression in developing wood led to an overall growth stimulation, xylem cell wall thinning and a shift in cellulose orientation. Transcriptome profiling of these transgenic lines indicated that cell cycling was stimulated and secondary wall biosynthesis was repressed. We suggest that the reduced xylan elongation is sensed by the cell wall integrity surveying mechanism in developing wood. Our results show that wood-specific suppression of xylan-biosynthetic GT43 genes activates signaling responses, leading to increased growth and improved lignocellulose saccharification.

ccNET: Database of co-expression networks with functional modules for diploid and polyploid Gossypium

Plant genera with both diploid and polyploid species are a common evolutionary occurrence. Polyploids, especially allopolyploids such as cotton and wheat, are a great model system for heterosis research. Here, we have integrated genome sequences and transcriptome data of Gossypium species to construct co-expression networks and identified functional modules from different cotton species, including 1155 and 1884 modules in G. arboreum and G. hirsutum, respectively. We overlayed the gene expression results onto the co-expression network. We further provided network comparison analysis for orthologous genes across the diploid and allotetraploid Gossypium. We also constructed miRNA-target networks and predicted PPI networks for both cotton species. Furthermore, we integrated in-house ChIP-seq data of histone modification (H3K4me3) together with cis-element analysis and gene sets enrichment analysis tools for studying possible gene regulatory mechanism in Gossypium species. Finally, we have constructed an online ccNET database (http://structuralbiology.cau.edu.cn/gossypium) for comparative gene functional analyses at a multi-dimensional network and epigenomic level across diploid and polyploid Gossypium species. The ccNET database will be beneficial for community to yield novel insights into gene/module functions during cotton development and stress response, and might be useful for studying conservation and diversity in other polyploid plants, such as T. aestivum and Brassica napus.

Functional conservation and divergence of Miscanthus lutarioriparius GT43 gene family in xylan biosynthesis

BACKGROUND

Xylan is the most abundant un-cellulosic polysaccharides of plant cell walls. Much progress in xylan biosynthesis has been gained in the model plant species Arabidopsis. Two homologous pairs Irregular Xylem 9 (IRX9)/9L and IRX14/14L from glycosyltransferase (GT) family 43 have been proved to play crucial roles in xylan backbone biosynthesis. However, xylan biosynthesis in grass such as Miscanthus remains poorly understood.

RESULTS

We characterized seven GT43 members in M. lutarioriparius, a promising bioenergy crop. Quantitative real-time RT-PCR (qRT-PCR) analysis revealed that the expression of MlGT43 genes was ubiquitously detected in the tissues examined. In-situ hybridization demonstrated that MlGT43A-B and MlGT43F-G were specifically expressed in sclerenchyma, while MlGT43C-E were expressed in both sclerenchyma and parenchyma. All seven MlGT43 proteins were localized to Golgi apparatus. Overexpression of MlGT43A-E but not MlGT43F and MlGT43G in Arabidopsis irx9 fully or partially rescued the mutant defects, including morphological changes, collapsed xylem and increased xylan contents, whereas overexpression of MlGT43F and MlGT43G but not MlGT43A-E complemented the defects of irx14, indicating that MlGT43A-E are functional orthologues of IRX9, while MlGT43F and MlGT43G are functional orthologues of IRX14. However, overexpression of all seven MlGT43 genes could not rescue the mucilage defects of irx14 seeds. Furthermore, transient transactivation analyses of MlGT43A-E reporters demonstrated that MlGT43A and MlGT43B but not MlGT43C-E were differentially activated by MlSND1, MlMYB46 or MlVND7.

CONCLUSION

The results demonstrated that all seven MlGT43s are functionally conserved in xylan biosynthesis during secondary cell wall formation but diversify in seed coat mucilage xylan biosynthesis. The results obtained provide deeper insight into xylan biosynthesis in grass, which lay the foundation for genetic modification of grass cell wall components and structure to better suit for next-generation biofuel production.

Cotton GhMYB7 is predominantly expressed in developing fibers and regulates secondary cell wall biosynthesis in transgenic Arabidopsis

The secondary cell wall in mature cotton fibers contains over 90% cellulose with low quantities of xylan and lignin. However, little is known regarding the regulation of secondary cell wall biosynthesis in cotton fibers. In this study, we characterized an R2R3-MYB transcription factor, GhMYB7, in cotton. GhMYB7 is expressed at a high level in developing fibers and encodes a MYB protein that is targeted to the cell nucleus and has transcriptional activation activity. Ectopic expression of GhMYB7 in Arabidopsis resulted in small, curled, dark green leaves and also led to shorter inflorescence stems. A cross-sectional assay of basal stems revealed that cell wall thickness of vessels and interfascicular fibers was higher in transgenic lines overexpressing GhMYB7 than in the wild type. Constitutive expression of GhMYB7 in Arabidopsis activated the expression of a suite of secondary cell wall biosynthesis-related genes (including some secondary cell wall-associated transcription factors), leading to the ectopic deposition of cellulose and lignin. The ectopic deposition of secondary cell walls may have been initiated before the cessation of cell expansion. Moreover, GhMYB7 was capable of binding to the promoter regions of AtSND1 and AtCesA4, suggesting that GhMYB7 may function upstream of NAC transcription factors. Collectively, these findings suggest that GhMYB7 is a potential transcriptional activator, which may participate in regulating secondary cell wall biosynthesis of cotton fibers.