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

Cytochemical localization and synthesis mechanism of the glucomannan in pseudobulbs of Bletilla striata Reichb. f

The dried pseudobulbs of Bletilla striata, an important traditional Chinese medicine named BaiJi, have an extraordinary polysaccharide content and excellent prospects for medicinal effects. However, the distribution and molecular mechanism underlying biosynthesis are poorly understood. In this study, chemical and immunologic analyses were performed in representative tissues of B. striata, and the results showed that what are conventionally termed Bletilla striata polysaccharides (BSPs) are water-soluble polysaccharides deposited only in pseudobulbs. The structural component of BSPs is glucomannan, with a mannose:glucose mass ratio of ~3:2. BSPs are present in the parenchyma of the pseudobulbs in cells known as glucomannan idioblasts and distributed in the cytoplasm within cellular membranes, but are not contained in the vacuole. Comparative transcriptomics and bioinformatics analyses mapped the pathway from sucrose to BSP and identified BsGPI, BsmanA, and BsCSLAs as the key genes of BSP biosynthesis, suggesting that the functional differentiation of the cellulose synthase-like family A (CSLA) may be critical for the flow of glucomannan to the BSP or cell wall. Subsequently, virus-mediated gene silencing showed that silencing of two CSLAs (Bs03G11846 and Bs03G11849) led to a decrease in BSP content, and yeast two-hybrid and luciferase complementation experiments confirmed that four CSLAs (Bs03G11846, Bs03G11847, Bs03G11848, and Bs03G11849) can form homo- or heterodimers, suggesting that multiple CSLAs may form a large complex that functions in BSP synthesis. Our results provide cytological evidence of BSP and describe the isolation and characterization of candidate genes involved in BSP synthesis, laying a solid foundation for further research on its regulation mechanisms and the genetic engineering breeding of B. striata.

An aldehyde dehydrogenase gene, GhALDH7B4_A06, positively regulates fiber strength in upland cotton (Gossypium hirsutum L.)

High fiber strength (FS) premium cotton has significant market demand. Consequently, enhancing FS is a major objective in breeding quality cotton. However, there is a notable lack of known functionally applicable genes that can be targeted for breeding. To address this issue, our study used specific length-amplified fragment sequencing combined with bulk segregant analysis to study FS trait in an F2 population. Subsequently, we integrated these results with previous quantitative trait locus mapping results regarding fiber quality, which used simple sequence repeat markers in F2, F2:3, and recombinant inbred line populations. We identified a stable quantitative trait locus qFSA06 associated with FS located on chromosome A06 (90.74-90.83 Mb). Within this interval, we cloned a gene, GhALDH7B4_A06, which harbored a critical mutation site in coding sequences that is distinct in the two parents of the tested cotton line. In the paternal parent Ji228, the gene is normal and referred to as GhALDH7B4_A06O; however, there is a nonsense mutation in the maternal parent Ji567 that results in premature termination of protein translation, and this gene is designated as truncated GhALDH7B4_A06S. Validation using recombinant inbred lines and gene expression analysis revealed that this mutation site is correlated with cotton FS. Virus-induced gene silencing of GhALDH7B4 in cotton caused significant decreases in FS and fiber micronaire. Conversely, GhALDH7B4_A06O overexpression in Arabidopsis boosted cell wall component contents in the stem. The findings of our study provide a candidate gene for improving cotton fiber quality through molecular breeding.

GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, play a regulatory role in cotton fiber secondary cell wall biosynthesis

KEY MESSAGE

GhHB14_D10 and GhREV_D5 regulated secondary cell wall formation and played an important role in fiber development. Cotton serves as an important source of natural fiber, and the biosynthesis of the secondary cell wall plays a pivotal role in determining cotton fiber quality. Nevertheless, the intricacies of this mechanism in cotton fiber remain insufficiently elucidated. This study investigates the functional roles of GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, in secondary cell wall biosynthesis in cotton fibers. Both GhHB14_D10 and GhREV_D5 were found to be localized in the nucleus with transcriptional activation activity. Ectopic overexpression of GhHB14_D10 and GhREV_D5 in Arabidopsis resulted in changed xylem differentiation, secondary cell wall deposition, and expression of genes related to the secondary cell wall. Silencing of GhHB14_D10 and GhREV_D5 in cotton led to enhanced fiber length, reduced cell wall thickness, cellulose contents and expression of secondary cell wall-related genes. Moreover, GhHB14_D10's direct interaction with GhREV_D5, and transcriptional regulation of cellulose biosynthesis genes GhCesA4-4 and GhCesA7-2 revealed their collaborative roles in secondary cell wall during cotton fiber development. Overall, these results shed light on the roles of GhHB14_D10 and GhREV_D5 in secondary cell wall biosynthesis, offering a strategy for the genetic improvement of cotton fiber quality.

The transcription factor PbrMYB24 regulates lignin and cellulose biosynthesis in stone cells of pear fruits

Lignified stone cell content is a key factor used to evaluate fruit quality, influencing the economic value of pear (Pyrus pyrifolia) fruits. However, our understanding of the regulatory networks of stone cell formation is limited due to the complex secondary metabolic pathway. In this study, we used a combination of co-expression network analysis, gene expression profiles, and transcriptome analysis in different pear cultivars with varied stone cell content to identify a hub MYB gene, PbrMYB24. The relative expression of PbrMYB24 in fruit flesh was significantly correlated with the contents of stone cells, lignin, and cellulose. We then verified the function of PbrMYB24 in regulating lignin and cellulose formation via genetic transformation in homologous and heterologous systems. We constructed a high-efficiency verification system for lignin and cellulose biosynthesis genes in pear callus. PbrMYB24 transcriptionally activated multiple target genes involved in stone cell formation. On the one hand, PbrMYB24 activated the transcription of lignin and cellulose biosynthesis genes by binding to different cis-elements [AC-I (ACCTACC) element, AC-II (ACCAACC) element and MYB-binding sites (MBS)]. On the other hand, PbrMYB24 bound directly to the promoters of PbrMYB169 and NAC STONE CELL PROMOTING FACTOR (PbrNSC), activating the gene expression. Moreover, both PbrMYB169 and PbrNSC activated the promoter of PbrMYB24, enhancing gene expression. This study improves our understanding of lignin and cellulose synthesis regulation in pear fruits through identifying a regulator and establishing a regulatory network. This knowledge will be useful for reducing the stone cell content in pears via molecular breeding.

Genome-Wide Analysis of WRKY Transcription Factors Involved in Abiotic Stress and ABA Response in Caragana korshinskii

The WRKY transcription factor family plays a vital role in plant development and environmental response. However, the information of WRKY genes at the genome-wide level is rarely reported in Caragana korshinskii. In this study, we identified and renamed 86 CkWRKY genes, which were further classified into three groups through phylogenetic analysis. Most of these WRKY genes were clustered and distributed on eight chromosomes. Multiple sequence alignment revealed that the conserved domain (WRKYGQK) of the CkWRKYs was basically consistent, but there were also six variation types (WRKYGKK, GRKYGQK, WRMYGQK, WRKYGHK, WKKYEEK and RRKYGQK) that appeared. The motif composition of the CkWRKYs was quite conservative in each group. In general, the number of WRKY genes gradually increased from lower to higher plant species in the evolutionary analysis of 28 species, with some exceptions. Transcriptomics data and RT-qPCR analysis showed that the CkWRKYs in different groups were involved in abiotic stresses and ABA response. Our results provided a basis for the functional characterization of the CkWRKYs involved in stress resistance in C. korshinskii.

Strigolactones act downstream of gibberellins to regulate fiber cell elongation and cell wall thickness in cotton (Gossypium hirsutum)

Strigolactones (SLs) are a class of phytohormones that regulate plant shoot branching and adventitious root development. However, little is known regarding the role of SLs in controlling the behavior of the smallest unit of the organism, the single cell. Here, taking advantage of a classic single-cell model offered by the cotton (Gossypium hirsutum) fiber cell, we show that SLs, whose biosynthesis is fine-tuned by gibberellins (GAs), positively regulate cell elongation and cell wall thickness by promoting the biosynthesis of very long-chain fatty acids (VLCFAs) and cellulose, respectively. Furthermore, we identified two layers of transcription factors (TFs) involved in the hierarchical regulation of this GA-SL crosstalk. The top-layer TF GROWTH-REGULATING FACTOR 4 (GhGRF4) directly activates expression of the SL biosynthetic gene DWARF27 (D27) to increase SL accumulation in fiber cells and GAs induce GhGRF4 expression. SLs induce the expression of four second-layer TF genes (GhNAC100-2, GhBLH51, GhGT2, and GhB9SHZ1), which transmit SL signals downstream to two ketoacyl-CoA synthase genes (KCS) and three cellulose synthase (CesA) genes by directly activating their transcription. Finally, the KCS and CesA enzymes catalyze the biosynthesis of VLCFAs and cellulose, respectively, to regulate development of high-grade cotton fibers. In addition to providing a theoretical basis for cotton fiber improvement, our results shed light on SL signaling in plant development at the single-cell level.

Transcription factor FfMYB15 regulates the expression of cellulase gene FfCEL6B during mycelial growth of Flammulina filiformis

Background

Cellulose degradation can determine mycelial growth rate and affect yield during the growth of Flammulina filiformis. The degradation of cellulose requires the joint action of a variety of cellulases, and some cellulase-related genes have been detected in mushrooms. However, little is known about the transcriptional regulatory mechanisms of cellulose degradation.

Results

In this study, FfMYB15 that may regulate the expression of cellulase gene FfCEL6B in F. filiformis was identified. RNA interference (RNAi) showed that FfCEL6B positively regulated mycelial growth. Gene expression analyses indicated that the expression patterns of FfCEL6B and FfMYB15 in mycelia cultured on the 0.9% cellulose medium for different times were similar with a correlation coefficient of 0.953. Subcellular localization and transcriptional activity analyses implied that FfMYB15 was located in the nucleus and was a transcriptional activator. Electrophoretic mobility shift assay (EMSA) and dual-luciferase assays demonstrated that FfMYB15 could bind and activate FfCEL6B promoter by recognizing MYB cis-acting element.

Conclusions

This study indicated that FfCEL6B played an active role in mycelial growth of F. filiformis and was regulated by FfMYB15.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01932-z.

A transcription factor TaMYB5 modulates leaf rolling in wheat

Leaf rolling is an important agronomic trait in wheat (Triticum aestivum L.). Moderate leaf rolling keeps leaves upright and maintains the relatively normal photosynthesis of plants under drought stress. However, the molecular mechanism of wheat leaf rolling remains unclear. Here, we identified a candidate gene TaMYB5-3A that regulates leaf rolling by using a genome-wide association study (GWAS) in a panel of 323 wheat accessions. Phenotype analysis indicated that the leaves of tamyb5 mutants were flatter than that of the wild type under drought condition. A nucleotide variation in the TaMYB5-3A coding region resulted in a substitution of Thr to Lys, which corresponds to two alleles SNP-3A-1 and SNP-3A-2. The leaf rolling index (LRI) of the SNP-3A-1 genotype was significantly lower than that of the SNP-3A-2 genotype. In addition, TaMYB5-3A alleles were associated with canopy temperature (CT) in multiple environments. The CT of the SNP-3A-1 genotype was lower than that of the SNP-3A-2 genotype. Gene expression analysis showed that TaMYB5-3A was mainly expressed in leaves and down-regulated by PEG and ABA treatment. TaMYB5 induces TaNRL1 gene expression through the direct binding to the AC cis-acting element of the promoter of the target gene, which was validated by EMSA (electrophoretic mobility shift assay). Our results revealed a crucial molecular mechanism in wheat leaf rolling and provided the theoretical basis and a gene resource for crop breeding.

Molecular Mechanisms of Plant Trichome Development

Plant trichomes, protrusions formed from specialized aboveground epidermal cells, provide protection against various biotic and abiotic stresses. Trichomes can be unicellular, bicellular or multicellular, with multiple branches or no branches at all. Unicellular trichomes are generally not secretory, whereas multicellular trichomes include both secretory and non-secretory hairs. The secretory trichomes release secondary metabolites such as artemisinin, which is valuable as an antimalarial agent. Cotton trichomes, also known as cotton fibers, are an important natural product for the textile industry. In recent years, much progress has been made in unraveling the molecular mechanisms of trichome formation in Arabidopsis thaliana, Gossypium hirsutum, Oryza sativa, Cucumis sativus, Solanum lycopersicum, Nicotiana tabacum, and Artemisia annua. Here, we review current knowledge of the molecular mechanisms underlying fate determination and initiation, elongation, and maturation of unicellular, bicellular and multicellular trichomes in several representative plants. We emphasize the regulatory roles of plant hormones, transcription factors, the cell cycle and epigenetic modifications in different stages of trichome development. Finally, we identify the obstacles and key points for future research on plant trichome development, and speculated the development relationship between the salt glands of halophytes and the trichomes of non-halophytes, which provides a reference for future studying the development of plant epidermal cells.

Molecular studies of cellulose synthase supercomplex from cotton fiber reveal its unique biochemical properties

Cotton fiber is a highly elongated and thickened single cell that produces large quantities of cellulose, which is synthesized and assembled into cell wall microfibrils by the cellulose synthase complex (CSC). In this study, we report that in cotton (Gossypium hirsutum) fibers harvested during secondary cell wall (SCW) synthesis, GhCesA 4, 7, and 8 assembled into heteromers in a previously uncharacterized 36-mer-like cellulose synthase supercomplex (CSS). This super CSC was observed in samples prepared using cotton fiber cells harvested during the SCW synthesis period but not from cotton stem tissue or any samples obtained from Arabidopsis. Knock-out of any of GhCesA 4, 7, and 8 resulted in the disappearance of the CSS and the production of fiber cells with no SCW thickening. Cotton fiber CSS showed significantly higher enzyme activity than samples prepared from knock-out cotton lines. We found that the microfibrils from the SCW of wild-type cotton fibers may contain 72 glucan chains in a bundle, unlike other plant materials studied. GhCesA4, 7, and 8 restored both the dwarf and reduced vascular bundle phenotypes of their orthologous Arabidopsis mutants, potentially by reforming the CSC hexamers. Genetic complementation was not observed when non-orthologous CesA genes were used, indicating that each of the three subunits is indispensable for CSC formation and for full cellulose synthase function. Characterization of cotton CSS will increase our understanding of the regulation of SCW biosynthesis.

Cloning and prokaryotic expression of WRKY48 from Caragana intermedia

Caragana intermedia (C. intermedia) is a kind of drought-tolerant leguminous shrub. WRKY transcription factors are one of the largest family of transcription factors in plants and play critical regulatory roles in stress tolerance and the development of plants. In our study, CiWRKY48 was cloned from C. intermedia, analyzed using bioinformatics software, and expressed with a prokaryotic expression system. The results showed that the open reading frame (ORF) of CiWRKY48 was 1158bp, the molecular weight (MW) was 42 kDa, and its subcellular localization was in the nucleus. Additionally, fusion protein was obtained, and confirmed by western blotting. The stress resistance of the pET30a-His-MBP-CiWRKY48 transformed Escherichia coli expression strain was reduced under mannitol and salt treatment, compared with the control. Overall, our findings provided a foundation for uncovering the function of CiWRKY48.

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.

Characterization and Interaction Analysis of the Secondary Cell Wall Synthesis-Related Transcription Factor PmMYB7 in Pinus massoniana Lamb

In vascular plants, the importance of R2R3-myeloblastosis (R2R3-MYB) transcription factors (TFs) in the formation of secondary cell walls (SCWs) has long been a controversial topic due to the lack of empirical evidence of an association between TFs and downstream target genes. Here, we found that the transcription factor PmMYB7, which belongs to the R2R3-MYB subfamily, is involved in lignin biosynthesis in Pinus massoniana. PmMYB7 was highly expressed in lignified tissues and upon abiotic stress. As a bait carrier, the PmMYB7 protein had no toxicity or autoactivation in the nucleus. Forty-seven proteins were screened from the P. massoniana yeast library. These proteins were predicted to be mainly involved in resistance, abiotic stress, cell wall biosynthesis, and cell development. We found that the PmMYB7 protein interacted with caffeoyl CoA 3-O-methyltransferase-2 (PmCCoAOMT2)—which is involved in lignin biosynthesis—but not with beta-1, 2-xylosyltransferase (PmXYXT1) yeast two-hybrid (Y2H) studies. Our in vivo coimmunoprecipitation (Co-IP) assay further showed that the PmMYB7 and PmCCoAOMT2 proteins could interact. Therefore, we concluded that PmMYB7 is an upstream TF that can interact with PmCCoAOMT2 in plant cells. These findings lay a foundation for further research on the function of PmMYB7, lignin biosynthesis and molecular breeding in P. massoniana.

Novel insights into water-deficit-responsive mRNAs and lncRNAs during fiber development in Gossypium hirsutum

BACKGROUND

The fiber yield and quality of cotton are greatly and periodically affected by water deficit. However, the molecular mechanism of the water deficit response in cotton fiber cells has not been fully elucidated.

RESULTS

In this study, water deficit caused a significant reduction in fiber length, strength, and elongation rate but a dramatic increase in micronaire value. To explore genome-wide transcriptional changes, fibers from cotton plants subjected to water deficit (WD) and normal irrigation (NI) during fiber development were analyzed by transcriptome sequencing. Analysis showed that 3427 mRNAs and 1021 long noncoding RNAs (lncRNAs) from fibers were differentially expressed between WD and NI plants. The maximum number of differentially expressed genes (DEGs) and lncRNAs (DERs) was identified in fibers at the secondary cell wall biosynthesis stage, suggesting that this is a critical period in response to water deficit. Twelve genes in cotton fiber were differentially and persistently expressed at ≥ five time points, suggesting that these genes are involved in both fiber development and the water-deficit response and could potentially be used in breeding to improve cotton resistance to drought stress. A total of 540 DEGs were predicted to be potentially regulated by DERs by analysis of coexpression and genomic colocation, accounting for approximately 15.76% of all DEGs. Four DERs, potentially acting as target mimics for microRNAs (miRNAs), indirectly regulated their corresponding DEGs in response to water deficit.

CONCLUSIONS

This work provides a comprehensive transcriptome analysis of fiber cells and a set of protein-coding genes and lncRNAs implicated in the cotton response to water deficit, significantly affecting fiber quality during the fiber development stage.

Comparative transcriptome analysis of fiber and nonfiber tissues to identify the genes preferentially expressed in fiber development in Gossypium hirsutum

Cotton is an important natural fiber crop and economic crop worldwide. The quality of cotton fiber directly determines the quality of cotton textiles. Identifying cotton fiber development-related genes and exploring their biological functions will not only help to better understand the elongation and development mechanisms of cotton fibers but also provide a theoretical basis for the cultivation of new cotton varieties with excellent fiber quality. In this study, RNA sequencing technology was used to construct transcriptome databases for different nonfiber tissues (root, leaf, anther and stigma) and fiber developmental stages (7 days post-anthesis (DPA), 14 DPA, and 26 DPA) of upland cotton Coker 312. The sizes of the seven transcriptome databases constructed ranged from 4.43 to 5.20 Gb, corresponding to approximately twice the genome size of Gossypium hirsutum (2.5 Gb). Among the obtained clean reads, 83.32% to 88.22% could be compared to the upland cotton TM-1 reference genome. By analyzing the differential gene expression profiles of the transcriptome libraries of fiber and nonfiber tissues, we obtained 1205, 1135 and 937 genes with significantly upregulated expression at 7 DPA, 14 DPA and 26 DPA, respectively, and 124, 179 and 213 genes with significantly downregulated expression. Subsequently, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analyses were performed, which revealed that these genes were mainly involved in catalytic activity, carbohydrate metabolism, the cell membrane and organelles, signal transduction and other functions and metabolic pathways. Through gene annotation analysis, many transcription factors and genes related to fiber development were screened. Thirty-six genes were randomly selected from the significantly upregulated genes in fiber, and expression profile analysis was performed using qRT-PCR. The results were highly consistent with the gene expression profile analyzed by RNA-seq, and all of the genes were specifically or predominantly expressed in fiber. Therefore, our RNA sequencing-based comparative transcriptome analysis will lay a foundation for future research to provide new genetic resources for the genetic engineering of improved cotton fiber quality and for cultivating new transgenic cotton germplasms for fiber quality improvement.

Silencing GhJUB1L1 (JUB1-like 1) reduces cotton (Gossypium hirsutum) drought tolerance

Drought stress massively restricts plant growth and the yield of crops. Reducing the deleterious effects of drought is necessary for agricultural industry. The plant-specific NAC (NAM, ATAF1/2 and CUC2) transcription factors (TFs) are widely involved in the regulation of plant development and stress response. One of the NAC TF, JUNGBRUNNEN1 (JUB1), has been reported to involve in drought resistance in Arabidopsis. However, little is known of how the JUB1 gene respond to drought stress in cotton. In the present study, we cloned GhJUB1L1, a homologous gene of JUB1 in upland cotton. GhJUB1L1 is preferentially expressed in stem and leaf and could be induced by drought stress. GhJUB1L1 protein localizes to the cell nucleus, and the transcription activation region of which is located in the C-terminal region. Silencing GhJUB1L1 gene via VIGS () reduced cotton drought tolerance, and retarded secondary cell wall (SCW) development. Additionally, the expression of some drought stress-related genes and SCW synthesis-related genes were altered in the GhJUB1L1 silencing plants. Collectively, our findings indicate that GhJUB1L1 may act as a positive regulator in response to drought stress and SCW development in cotton. Our results enriched the roles of NAC TFs in cotton drought tolerance and laid a foundation for the cultivation of transgenic cotton with higher drought tolerance.

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.

The miR319-Targeted GhTCP4 Promotes the Transition from Cell Elongation to Wall Thickening in Cotton Fiber

Plant cell growth involves a complex interplay among cell-wall expansion, biosynthesis, and, in specific tissues, secondary cell wall (SCW) deposition, yet the coordination of these processes remains elusive. Cotton fiber cells are developmentally synchronous, highly elongated, and contain nearly pure cellulose when mature. Here, we report that the transcription factor GhTCP4 plays an important role in balancing cotton fiber cell elongation and wall synthesis. During fiber development the expression of miR319 declines while GhTCP4 transcript levels increase, with high levels of the latter promoting SCW deposition. GhTCP4 interacts with a homeobox-containing factor, GhHOX3, and repressing its transcriptional activity. GhTCP4 and GhHOX3 function antagonistically to regulate cell elongation, thereby establishing temporal control of fiber cell transition to the SCW stage. We found that overexpression of GhTCP4A upregulated and accelerated activation of the SCW biosynthetic pathway in fiber cells, as revealed by transcriptome and promoter activity analyses, resulting in shorter fibers with varied lengths and thicker walls. In contrast, GhTCP4 downregulation led to slightly longer fibers and thinner cell walls. The GhHOX3-GhTCP4 complex may represent a general mechanism of cellular development in plants since both are conserved factors in many species, thus providing us a potential molecular tool for the design of fiber traits.

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 expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

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.

Genome-wide analysis of BES1 genes in Gossypium revealed their evolutionary conserved roles in brassinosteroid signaling

Brassinosteroids (BRs), which are essential phytohormones for plant growth and development, are important for cotton fiber development. Additionally, BES1 transcription factors are critical for BR signal transduction. However, cotton BES1 family genes have not been comprehensively characterized. In this study, we identified 11 BES1 genes in G. arboreum, 11 in G. raimondii, 16 in G. barbadense, and 22 in G. hirsutum. The BES1 sequences were significantly conserved in the Arabidopsis thaliana, rice, and upland cotton genomes. A total of 94 BES1 genes from 10 different plant species were divided into three clades according to the neighbor-joining and minimum-evolution methods. Moreover, the exon/intron patterns and motif distributions were highly conserved among the A. thaliana and cotton BES1 genes. The collinearity among the orthologs from the At and Dt subgenomes was estimated. Segmental duplications in the At and Dt subgenomes were primarily responsible for the expansion of the cotton BES1 gene family. Of the GhBES1 genes, GhBES1.4_At/Dt exhibited BL-induced expression and was predominantly expressed in fibers. Furthermore, Col-0/mGhBES1.4_At plants produced curled leaves with long and bent petioles. These transgenic plants also exhibited decreased hypocotyl sensitivity to brassinazole and constitutive BR induced/repressed gene expression patterns. The constitutive BR responses of the plants overexpressing mGhBES1.4_At were similar to those of the bes1-D mutant.

GWAS Analysis and QTL Identification of Fiber Quality Traits and Yield Components in Upland Cotton Using Enriched High-Density SNP Markers

It is of great importance to identify quantitative trait loci (QTL) controlling fiber quality traits and yield components for future marker-assisted selection (MAS) and candidate gene function identifications. In this study, two kinds of traits in 231 F6:8 recombinant inbred lines (RILs), derived from an intraspecific cross between Xinluzao24, a cultivar with elite fiber quality, and Lumianyan28, a cultivar with wide adaptability and high yield potential, were measured in nine environments. This RIL population was genotyped by 122 SSR and 4729 SNP markers, which were also used to construct the genetic map. The map covered 2477.99 cM of hirsutum genome, with an average marker interval of 0.51 cM between adjacent markers. As a result, a total of 134 QTLs for fiber quality traits and 122 QTLs for yield components were detected, with 2.18-24.45 and 1.68-28.27% proportions of the phenotypic variance explained by each QTL, respectively. Among these QTLs, 57 were detected in at least two environments, named stable QTLs. A total of 209 and 139 quantitative trait nucleotides (QTNs) were associated with fiber quality traits and yield components by four multilocus genome-wide association studies methods, respectively. Among these QTNs, 74 were detected by at least two algorithms or in two environments. The candidate genes harbored by 57 stable QTLs were compared with the ones associated with QTN, and 35 common candidate genes were found. Among these common candidate genes, four were possibly "pleiotropic." This study provided important information for MAS and candidate gene functional studies.

The cotton (Gossypium hirsutum) NAC transcription factor (FSN1) as a positive regulator participates in controlling secondary cell wall biosynthesis and modification of fibers

Cotton (Gossypium hirsutum) fibers are the highly elongated and thickened single-cell trichomes on the seed epidermis. However, little is known about the molecular base of fiber cell wall thickening in detail. In this study, a cotton NAC transcription factor (GhFSN1) that is specifically expressed in secondary cell wall (SCW) thickening fibers was functionally characterized. The GhFSN1 transgenic cotton plants were generated to study how FSN1 regulates fiber SCW formation. Up-regulation of GhFSN1 expression in cotton resulted in an increase in SCW thickness of fibers but a decrease in fiber length. Transcriptomic analysis revealed that GhFSN1 activates or represses numerous downstream genes. GhFSN1 has the ability to form homodimers, binds to its promoter to activate itself, and might be degraded by the ubiquitin-mediated proteasome pathway. The direct targets of GhFSN1 include the fiber SCW-related GhDUF231L1, GhKNL1, GhMYBL1, GhGUT1 and GhIRX12 genes. GhFSN1 binds directly to a consensus sequence (GhNBS), (C/T)(C/G/T)TN(A/T)(G/T)(A/C/G)(A/G)(A/T/G)(A/T/G)AAG, which exists in the promoters of these SCW-related genes. Our data demonstrate that GhFSN1 acts as a positive regulator in controlling SCW formation of cotton fibers by activating its downstream SCW-related genes. Thus, these findings give us novel insights into comprehensive understanding of GhFSN1 function in fiber development.

Identification and Characterization of the Diverse Stress-Responsive R2R3-RMYB Transcription Factor from Hibiscus sabdariffa L

Various regulatory proteins play a fundamental role to manage the healthy plant growth under stress conditions. Differential display reverse transcriptase PCR and random amplification of cDNA ends (RACE) was used to explore the osmotic stress-responsive transcripts. We identified and characterized the salt stress-responsive R2R3 type RMYB transcription factor from Hibiscus sabdariffa which has an open reading frame of 690 bp, encoding 229 long chain amino acids. In silico analysis confirmed the conserved R2 and R3 domain as well as an NLS-1 localization site. The deduced amino acids of RMYB shared 83, 81, 80, 79, 72, 71, and 66% homology with Arabidopsis thaliana, Glycine max, Oryza sativa, Zea maize, Malus domestica, Populus tremula × Populus alba, and Medicago sativa specific MYB family, respectively. We observed the gene upregulation in stem, leaf, and root tissue in response to abiotic stress. Furthermore, RMYB gene was cloned into plant expression vector under CaMV35S promoter and transformed to Gossypium hirsutum: a local cotton cultivar. Overexpression of RMYB was observed in transgenic plants under abiotic stresses which further suggests its regulatory role in response to stressful conditions. The RMYB transcription factor-overexpressing in transgenic cotton plants may be used as potential agent for the development of stress tolerant crop cultivars.

Cotton functional genomics reveals global insight into genome evolution and fiber development

Due to the economic value of natural textile fiber, cotton has attracted much research attention, which has led to the publication of two diploid genomes and two tetraploid genomes. These big data facilitate functional genomic study in cotton, and allow researchers to investigate cotton genome structure, gene expression, and protein function on the global scale using high-throughput methods. In this review, we summarized recent studies of cotton genomes. Population genomic analyses revealed the domestication history of cultivated upland cotton and the roles of transposable elements in cotton genome evolution. Alternative splicing of cotton transcriptomes was evaluated genome-widely. Several important gene families like MYC, NAC, Sus and GhPLDα1 were systematically identified and classified based on genetic structure and biological function. High-throughput proteomics also unraveled the key functional proteins correlated with fiber development. Functional genomic studies have provided unprecedented insights into global-scale methods for cotton research.

Genome-wide identification and characterization of phospholipase C gene family in cotton (Gossypium spp.)

Phospholipase C (PLC) are important regulatory enzymes involved in several lipid and Ca²⁺-dependent signaling pathways. Previous studies have elucidated the versatile roles of PLC genes in growth, development and stress responses of many plants, however, the systematic analyses of PLC genes in the important fiber-producing plant, cotton, are still deficient. In this study, through genome-wide survey, we identified twelve phosphatidylinositol-specific PLC (PI-PLC) and nine non-specific PLC (NPC) genes in the allotetraploid upland cotton Gossypium hirsutum and nine PI-PLC and six NPC genes in two diploid cotton G. arboretum and G.raimondii, respectively. The PI-PLC and NPC genes of G. hirsutum showed close phylogenetic relationship with their homologous genes in the diploid cottons and Arabidopsis. Segmental and tandem duplication contributed greatly to the formation of the gene family. Expression profiling indicated that few of the PLC genes are constitutely expressed, whereas most of the PLC genes are preferentially expressed in specific tissues and abiotic stress conditions. Promoter analyses further implied that the expression of these PLC genes might be regulated by MYB transcription factors and different phytohormones. These results not only suggest an important role of phospholipase C members in cotton plant development and abiotic stress response but also provide good candidate targets for future molecular breeding of superior cotton cultivars.

Fine-mapping qFS07.1 controlling fiber strength in upland cotton (Gossypium hirsutum L.)

KEY MESSAGE: qFS07.1 controlling fiber strength was fine-mapped to a 62.6-kb region containing four annotated genes. RT-qPCR and sequence of candidate genes identified an LRR RLK gene as the most likely candidate. Fiber strength is an important component of cotton fiber quality and is associated with other properties, such as fiber maturity, fineness, and length. Stable QTL qFS07.1, controlling fiber strength, had been identified on chromosome 7 in an upland cotton recombinant inbred line (RIL) population from a cross (CCRI35 × Yumian1) described in our previous studies. To fine-map qFS07.1, an F₂ population with 2484 individual plants from a cross between recombinant line RIL014 and CCRI35 was established. A total of 1518 SSR primer pairs, including 1062, designed from chromosome 1 of the Gossypium raimondii genome and 456 from chromosome 1 of the G. arboreum genome (corresponding to the QTL region) were used to fine-map qFS07.1, and qFS07.1 was mapped into a 62.6-kb genome region which contained four annotated genes on chromosome A07 of G. hirsutum. RT-qPCR and comparative analysis of candidate genes revealed a leucine-rich repeat protein kinase (LRR RLK) family protein to be a promising candidate gene for qFS07.1. Fine mapping and identification of the candidate gene for qFS07.1 will play a vital role in marker-assisted selection (MAS) and the study of mechanism of cotton fiber development.

Molecular characterization of GhPLDα1 and its relationship with secondary cell wall thickening in cotton fibers

Phospholipase D (PLD) hydrolyzes phospholipids to generate a free polar head group (e.g., choline) and a second messenger phosphatidic acid and plays diverse roles in plant growth and development, including seed germination, leaf senescence, root hair growth, and hypocotyl elongation. However, the function of PLD in cotton remains largely unexplored. Here, the comprehensive molecular characterization of GhPLDα1 was explored with its role in upland cotton (Gossypium hirsutum) fiber development. The GhPLDα1 gene was cloned successfully, and a sequence alignment showed that GhPLDα1 contains one C2 domain and two HKD (HxKxxxxD) domains. Quantitative reverse transcriptase-polymerase chain reaction measured the expression of GhPLDα1 in various cotton tissues with the highest level in fibers at 20 days post anthesis (d.p.a.). Fluorescent microscopy and immunoblotting in tobacco epidermis showed the GhPLDα1 distribution in both cell membranes and the cytoplasm. An activity assay indicated changes in PLDα enzyme activity in developing fiber cells with a peak level at 20 d.p.a., coinciding with the onset of cellulose accumulation and the increased H₂O₂ content during fiber development. Furthermore, the inhibition of PLDα activity obviously decreased the cellulose and H₂O₂ contents of in vitro-cultured cotton fibers. These results provide important evidence explaining the relationship of GhPLDα1 with secondary cell wall thickening in cotton fibers in that GhPLDα1 may correlate with the increased H₂O₂ content at the onset of secondary cell wall thickening, ultimately promoting cellulose biosynthesis.