Live-cell imaging of the cytoskeleton in elongating cotton fibres

The transcription factor GhMYB4 represses lipid transfer and sucrose transporter genes and inhibits fiber cell elongation in cotton

Cotton (Gossypium hirsutum) fiber is a highly elongated single cell with a thickened cell wall. MYB transcription factors are important regulators of plant cell elongation; however, the molecular mechanism involved in regulating fiber elongation remains to be explored. Here, we present evidence that the R2R3-MYB transcription factor GhMYB4 negatively regulates cotton fiber cell elongation by suppressing the expression of 2 crucial genes previously reported to affect fiber development: lipid transfer protein 4 (GhLTP4) and sucrose transporter 12 (GhSWEET12). GhMYB4 is preferentially expressed in elongating fiber cells. Knockdown of GhMYB4 in cotton results in longer fiber cells, whereas overexpression of GhMYB4 in Arabidopsis leads to reduced plant height and root length. Transcriptomic and lipidomic analyses revealed that GhMYB4 is involved in coordinating 3 interconnected biological processes, namely lipid content regulation, auxin signaling, and sugar metabolism. Additionally, we showed that GhMYB4 inhibits the expression of GhLTP4 and GhSWEET12 by binding to the MYB cis-element (TTTAGTG) in their respective promoters. Interestingly, basic helix-loop-helix transcription factor 105 (GhbHLH105) and MYB transcription factor 212 (GhMYB212) counteract the inhibitory effects of GhMYB4 on the expression of GhLTP4 and GhSWEET12, respectively. These findings provide insights into the complex molecular mechanisms regulating plant cell elongation.

Sphingolipid inhibitor response gene GhMYB86 controls fiber elongation by regulating microtubule arrangement

Although the cell membrane and cytoskeleton play essential roles in cellular morphogenesis, the interaction between the membrane and cytoskeleton is poorly understood. Cotton fibers are extremely elongated single cells, which makes them an ideal model for studying cell development. Here, we used the sphingolipid biosynthesis inhibitor, fumonisin B1 (FB1), and found that it effectively suppressed the myeloblastosis (MYB) transcription factor GhMYB86, thereby negatively affecting fiber elongation. A direct target of GhMYB86 is GhTUB7, which encodes the tubulin protein, the major component of the microtubule cytoskeleton. Interestingly, both the overexpression of GhMYB86 and GhTUB7 caused an ectopic microtubule arrangement at the fiber tips, and then leading to shortened fibers. Moreover, we found that GhMBE2 interacted with GhMYB86 and that FB1 and reactive oxygen species induced its transport into the nucleus, thereby enhancing the promotion of GhTUB7 by GhMYB86. Overall, we established a GhMBE2-GhMYB86-GhTUB7 regulation module for fiber elongation and revealed that membrane sphingolipids affect fiber elongation by altering microtubule arrangement.

Cotton Pectate Lyase GhPEL48_Dt Promotes Fiber Initiation Mediated by Histone Acetylation

GhPEL48_Dt, a Pectate lyase (PEL, EC4.2.2.2), is a crucial enzyme involved in cell-wall modification and pectin degradation. Studies have shown that the GhPEL48_Dt also plays a significant role in cotton-fiber development; however, the specific function and regulatory mechanism of GhPEL48_Dt in cotton-fiber development are still not fully understood. Here, we found that the histone deacetylase inhibitor-Trichostatin A significantly reduces the transcript levels of GhPEL48_Dt and its enzyme activity. Further, silencing of GhPEL48_Dt significantly inhibits the initiation and elongation of cotton fibers by promoting pectin degradation, and the heterologous expression of GhPEL48_Dt promotes the development of trichomes and root hairs in Arabidopsis, which suggests that GhPEL48_Dt plays a positive and conserved role in single cell i.e., fiber, root hair, and leaf trichome development. Collectively, this paper provides a comprehensive analysis of the fundamental characteristics and functions of GhPEL48_Dt in fiber development, including the regulatory role of histone acetylation on GhPEL48_Dt, which contributes to the understanding of pectin degradation pathways and establishes a theoretical foundation for elucidating its regulatory mechanism.

Deciphering the dynamic expression network of fiber elongation and the functional role of the GhTUB5 gene for fiber length in cotton based on an introgression population of upland cotton

INTRODUCTION

Interspecific introgression between Gossypium hirsutum and G. barbadense allows breeding cotton with outstanding fiber length (FL). However, the dynamic gene regulatory network of FL-related genes has not been characterized, and the functional mechanism through which the hub gene GhTUB5 mediates fiber elongation has yet to be determined.

METHODS

Coexpression analyses of 277 developing fiber transcriptomes integrated with QTL mapping using 250 introgression lines of different FL phenotypes were conducted to identify genes related to fiber elongation. The function of GhTUB5 was determined by ectopic expression of two TUB5 alleles in Arabidopsis and knockout of GhTUB5 in upland cotton. Yeast two-hybrid, split-luciferase and pull-down assays were conducted to screen for interacting proteins, and upstream genes were identified by yeast one-hybrid, dual-LUC and electrophoretic mobility shift assays.

RESULTS

The 32,612, 30,837 and 30,277 genes expressed at 5, 10 and 15 days postanthesis (dpa) were grouped into 19 distinct coexpression modules, and 988 genes in the MEblack module were enriched in the cell wall process and exhibited significant associations with FL. A total of 20 FL-QTLs were identified, each explaining 3.34-16.04 % of the phenotypic variance in the FL. Furthermore, several FL-QTLs contained 15 genes that were differentially expressed in the MEblack module including the tubulin beta gene (TUB5). Compared with the wild type, the overexpression of GhTUB5 and GbTUB5 in Arabidopsis suppressed root cell length but promoted cellulose synthesis. Knockout of GhTUB5 resulted in longer fiber lines. Protein-based experiments revealed that GhTUB5 interacts with GhZFP6. Additionally, GhTUB5 was directly activated by GhHD-ZIP7, a homeobox-leucine zipper transcription factor, and its paralogous gene was previously reported to mediate fiber elongation.

CONCLUSION

This study opens a new avenue to dissect functional mechanism of cotton fiber elongation. Our findings provide some molecular details on how GhTUB5 mediates the FL phenotype in cotton.

Unraveling the lncRNA-miRNA-mRNA Regulatory Network Involved in Poplar Coma Development through High-Throughput Sequencing

Poplar coma, the fluff-like appendages of seeds originating from the differentiated surface cells of the placenta and funicle, aids in the long-distance dispersal of seeds in the spring. However, it also poses hazards to human safety and causes pollution in the surrounding environment. Unraveling the regulatory mechanisms governing the initiation and development of coma is essential for addressing this issue comprehensively. In this study, strand-specific RNA-seq was conducted at three distinct stages of coma development, revealing 1888 lncRNAs and 52,810 mRNAs. The expression profiles of lncRNAs and mRNAs during coma development were analyzed. Subsequently, potential target genes of lncRNAs were predicted through co-localization and co-expression analyses. Integrating various types of sequencing data, lncRNA-miRNA-TF regulatory networks related to the initiation of coma were constructed. Utilizing identified differentially expressed genes encoding kinesin and actin, lncRNA-miRNA-mRNA regulatory networks associated with the construction and arrangement of the coma cytoskeleton were established. Additionally, relying on differentially expressed genes encoding cellulose synthase, sucrose synthase, and expansin, lncRNA-miRNA-mRNA regulatory networks related to coma cell wall synthesis and remodeling were developed. This study not only enhances the comprehension of lncRNA but also provides novel insights into the molecular mechanisms governing the initiation and development of poplar coma.

Transmembrane formins as active cargoes of membrane trafficking

Formins are a large, evolutionarily old family of cytoskeletal regulators whose roles include actin capping and nucleation, as well as modulation of microtubule dynamics. The plant class I formin clade is characterized by a unique domain organization, as most of its members are transmembrane proteins with possible cell wall-binding motifs exposed to the extracytoplasmic space-a structure that appears to be a synapomorphy of the plant kingdom. While such transmembrane formins are traditionally considered mainly as plasmalemma-localized proteins contributing to the organization of the cell cortex, we review, from a cell biology perspective, the growing evidence that they can also, at least temporarily, reside (and in some cases also function) in endomembranes including secretory and endocytotic pathway compartments, the endoplasmic reticulum, the nuclear envelope, and the tonoplast. Based on this evidence, we propose that class I formins may thus serve as 'active cargoes' of membrane trafficking-membrane-embedded proteins that modulate the fate of endo- or exocytotic compartments while being transported by them.

The synergistic effect of multiple organic macromolecules on the formation of calcium oxalate raphides of Musa spp

Needle-like calcium oxalate crystals called raphides are unique structures in the plant kingdom. Multiple biomacromolecules work together in the regulatory and transportation pathways to form raphides; however, the mechanism by which this occurs remains unknown. Using banana (Musa spp.), this study combined in vivo methods including confocal microscopy, transmission electron microscopy, and Q Exactive mass spectrometry to identify the main biomolecules, such as vesicles, together with the compositions of lipids and proteins in the crystal chamber, which is the membrane compartment that surrounds each raphide during its formation. Simulations of the vesicle transportation process and the synthesis of elongated calcium oxalate crystals in vitro were then conducted, and the results suggested that the vesicles carrying amorphous calcium oxalate and proteins embedded in raphides are transported along actin filaments. These vesicles subsequently fuse with the crystal chamber, utilizing the proteins embedded in the raphides as a template for the final formation of the structure. Our findings contribute to the fundamental understanding of the regulation of the diverse biomacromolecules that are crucial for raphide formation. Moreover, the implications of these findings extend to other fields such as materials science, and particularly the synthesis of functionalized materials.

The Microtubule End Binding Protein Mal3 Is Essential for the Dynamic Assembly of Microtubules during Magnaporthe oryzae Growth and Pathogenesis

Cytoskeletal microtubules (MTs) play crucial roles in many aspects of life processes in eukaryotic organisms. They dynamically assemble physiologically important MT arrays under different cell conditions. Currently, aspects of MT assembly underlying the development and pathogenesis of the model plant pathogenic fungus Magnaporthe oryzae (M. oryzae) are unclear. In this study, we characterized the MT plus end binding protein MoMal3 in M. oryzae. We found that knockout of MoMal3 results in defects in hyphal polar growth, appressorium-mediated host penetration and nucleus division. Using high-resolution live-cell imaging, we further found that the MoMal3 mutant assembled a rigid MT in parallel with the MT during hyphal polar growth, the cage-like network in the appressorium and the stick-like spindle in nuclear division. These aberrant MT organization patterns in the MoMal3 mutant impaired actin-based cell growth and host infection. Taken together, these findings showed that M. oryzae relies on MoMal3 to assemble elaborate MT arrays for growth and infection. The results also revealed the assembly mode of MTs in M. oryzae, indicating that MTs are pivotal for M. oryzae growth and host infection and may be new targets for devastating fungus control.

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.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Brassinosteroids regulate cotton fiber elongation by modulating very-long-chain fatty acid biosynthesis

Brassinosteroid (BR), a growth-promoting phytohormone, regulates many plant growth processes including cell development. However, the mechanism by which BR regulates fiber growth is poorly understood. Cotton (Gossypium hirsutum) fibers are an ideal single-cell model in which to study cell elongation due to their length. Here we report that BR controls cotton fiber elongation by modulating very-long-chain fatty acid (VLCFA) biosynthesis. BR deficiency reduces the expression of 3-ketoacyl-CoA synthases (GhKCSs), the rate-limiting enzymes involved in VLCFA biosynthesis, leading to lower saturated VLCFA contents in pagoda1 (pag1) mutant fibers. In vitro ovule culture experiments show that BR acts upstream of VLCFAs. Silencing of BRI1-EMS-SUPPRESOR 1.4 (GhBES1.4), encoding a master transcription factor of the BR signaling pathway, significantly reduces fiber length, whereas GhBES1.4 overexpression produces longer fibers. GhBES1.4 regulates endogenous VLCFA contents and directly binds to BR RESPONSE ELEMENTS (BRREs) in the GhKCS10_At promoter region, which in turn regulates GhKCS10_At expression to increase endogenous VLCFA contents. GhKCS10_At overexpression promotes cotton fiber elongation, whereas GhKCS10_At silencing inhibits cotton fiber growth, supporting a positive regulatory role for GhKCS10_At in fiber elongation. Overall, these results uncover a mechanism of fiber elongation through crosstalk between BR and VLCFAs at the single-cell level.

Transcriptome, Ectopic Expression and Genetic Population Analysis Identify Candidate Genes for Fiber Quality Improvement in Cotton

Comparative transcriptome analysis of fiber tissues between Gossypium barbadense and Gossypium hirsutum could reveal the molecular mechanisms underlying high-quality fiber formation and identify candidate genes for fiber quality improvement. In this study, 759 genes were found to be strongly upregulated at the elongation stage in G. barbadense, which showed four distinct expression patterns (I-IV). Among them, the 346 genes of group IV stood out in terms of the potential to promote fiber elongation, in which we finally identified 42 elongation-related candidate genes by comparative transcriptome analysis between G. barbadense and G. hirsutum. Subsequently, we overexpressed GbAAR3 and GbTWS1, two of the 42 candidate genes, in Arabidopsis plants and validated their roles in promoting cell elongation. At the secondary cell wall (SCW) biosynthesis stage, 2275 genes were upregulated and exhibited five different expression profiles (I-V) in G. barbadense. We highlighted the critical roles of the 647 genes of group IV in SCW biosynthesis and further picked out 48 SCW biosynthesis-related candidate genes by comparative transcriptome analysis. SNP molecular markers were then successfully developed to distinguish the SCW biosynthesis-related candidate genes from their G. hirsutum orthologs, and the genotyping and phenotyping of a BC3F5 population proved their potential in improving fiber strength and micronaire. Our results contribute to the better understanding of the fiber quality differences between G. barbadense and G. hirsutum and provide novel alternative genes for fiber quality improvement.

VILLIN2 regulates cotton defense against Verticillium dahliae by modulating actin cytoskeleton remodeling

The active structural change of actin cytoskeleton is a general host response upon pathogen attack. This study characterized the function of the cotton (Gossypium hirsutum) actin-binding protein VILLIN2 (GhVLN2) in host defense against the soilborne fungus Verticillium dahliae. Biochemical analysis demonstrated that GhVLN2 possessed actin-binding, -bundling, and -severing activities. A low concentration of GhVLN2 could shift its activity from actin bundling to actin severing in the presence of Ca2+. Knockdown of GhVLN2 expression by virus-induced gene silencing reduced the extent of actin filament bundling and interfered with the growth of cotton plants, resulting in the formation of twisted organs and brittle stems with a decreased cellulose content of the cell wall. Upon V. dahliae infection, the expression of GhVLN2 was downregulated in root cells, and silencing of GhVLN2 enhanced the disease tolerance of cotton plants. The actin bundles were less abundant in root cells of GhVLN2-silenced plants than in control plants. However, upon infection by V. dahliae, the number of actin filaments and bundles in the cells of GhVLN2-silenced plants was raised to a comparable level as those in control plants, with the dynamic remodeling of the actin cytoskeleton appearing several hours in advance. GhVLN2-silenced plants exhibited a higher incidence of actin filament cleavage in the presence of Ca2+, suggesting that pathogen-responsive downregulation of GhVLN2 could activate its actin-severing activity. These data indicate that the regulated expression and functional shift of GhVLN2 contribute to modulating the dynamic remodeling of the actin cytoskeleton in host immune responses against V. dahliae.

Studying the Organization of the Actin Cytoskeleton in the Multicellular Trichomes of Tomato

Trichomes are unique polarized cells of the plant epidermis that provide an excellent model system for studying the plant cytoskeleton. Unlike Arabidopsis trichomes that are unicellular with a typical three-branch shape, the trichomes in tomato (Solanum lycopersicum) are multicellular with additional morphology and function diversity. Technically, it is hard to image tomato trichomes at the subcellular level because of their size and because they can be easily damaged. Here, we describe the methods we have used for the visualization and quantification of cytoskeletal arrangements in tomato trichomes which are at different developmental stages, using both live-cell imaging and phalloidin staining after fixation.

Visualization of Cytoskeleton Organization and Dynamics in Elongating Cotton Fibers by Live-Cell Imaging

Cotton fibers are extremely elongated single cells and have long been regarded as an ideal model to investigate polarized plant cell elongation. Actin filaments (F-actin), as well as the cortical microtubules (CMTs), play vital roles in polarized cell growth and morphogenesis. We have generated stable transgenic cotton plants expressing fluorescent markers for the actin and microtubule cytoskeletons. Further live-cell imaging identified dynamic features of the F-actin and cortical microtubule (CMT) architectures and discovered that cotton fibers elongate in a unique tip-biased diffuse growth mode. Here, we describe methods for preparing growing cotton fiber samples, as well as the visualization of cytoskeletal organization and dynamics by live-cell imaging. Combined with comprehensive image analyses, these methods can be used to identify how cytoskeleton organization and dynamics determine cell morphogenesis in highly polarized cotton fibers.

Cotton microtubule-associated protein GhMAP20L5 mediates fiber elongation through the interaction with the tubulin GhTUB13

Targeting proteins for Xklp2 (TPX2s) comprise a class of MAPs that are essential for plant growth and development by regulating the dynamic changes of microtubules (MTs) and proper formation of cytoskeleton. However, the function of TPX2 proteins in cotton fiber development remains poorly understood. Here, we identified the function of a fiber elongation-specific TPX2 protein, GhMAP20L5, in cotton. Suppressed GhMAP20L5 gene expression in cotton (GhMAP20L5i) significantly reduced fiber elongation rate, fiber length and lint percentage. GhMAP20L5i fibers had thinner and looser secondary cell walls (SCW), and incompact helix twists. GhMAP20L5 specifically interacted with the tubulin GhTUB13 on the cytoskeleton. Gene coexpression analysis showed that GhMAP20L5 involved in multiple pathways related to cytoskeleton establishment and fiber cell wall formation and affected cellulase genes expressions. In summary, our results revealed that GhMAP20L5 is important for fiber development by regulating cytoskeleton establishment and the cellulose deposition in cotton.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes for Fuzz Development in Gossypium hirsutum

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

Single-cell RNA-seq reveals fate determination control of an individual fibre cell initiation in cotton (Gossypium hirsutum)

Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single-cell RNA sequencing (scRNA-seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non-fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA-seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25-like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip-biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

GhRabA4c coordinates cell elongation via regulating actin filament–dependent vesicle transport

GhRabA4c is required for cotton fiber cell elongation via functioning in actin filament assembly and bundling, vesicle transport, and deposition of multiple cell wall components.

Plant cell expands via a tip growth or diffuse growth mode. In plants, RabA is the largest group of Rab GTPases that regulate vesicle trafficking. The functions of RabA protein in modulating polarized expansion in tip growth cells have been demonstrated. However, whether and how RabA protein functions in diffuse growth plant cells have never been explored. Here, we addressed this question by examining the role of GhRabA4c in cotton fibers. GhRabA4c was preferentially expressed in elongating fibers with its protein localized to endoplasmic reticulum and Golgi apparatus. Over- and down-expression of GhRabA4c in cotton lead to longer and shorter fibers, respectively. GhRabA4c interacted with GhACT4 to promote the assembly of actin filament to facilitate vesicle transport for cell wall synthesis. Consistently, GhRabA4c-overexpressed fibers exhibited increased content of wall components and the transcript levels of the genes responsible for the synthesis of cell wall materials. We further identified two MYB proteins that directly regulate the transcription of GhRabA4c. Collectively, our data showed that GhRabA4c promotes diffused cell expansion by supporting vesicle trafficking and cell wall synthesis.

Fluorescent cytoskeletal markers reveal associations between the actin and microtubule cytoskeleton in rice cells

Rice (Oryza sativa) is one of our main food crops, feeding ∼3.5 billion people worldwide. An increasing number of studies note the importance of the cytoskeleton, including actin filaments and microtubules, in rice development and environmental responses. Yet, reliable in vivo cytoskeleton markers are lacking in rice, which limits our knowledge of cytoskeletal functions in living cells. Therefore, we generated bright fluorescent marker lines of the actin and microtubule cytoskeletons in rice, suitable for live-cell imaging in a wide variety of rice tissues. Using these lines, we show that actin bundles and microtubules engage and co-function during pollen grain development, how the cytoskeletal components are coordinated during root cell development, and that the actin cytoskeleton is robust and facilitates microtubule responses during salt stress. Hence, we conclude that our cytoskeletal marker lines, highlighted by our findings of cytoskeletal associations and dynamics, will substantially further future investigations in rice biology.

Microtubules play a crucial role in regulating actin organization and cell initiation in cotton fibers

Dynamic organization of actin and microtubule cytoskeletons directs a distinct expansion behavior of cotton fiber initiation from cell elongation. Cotton fibers are highly elongated single cells derived from the ovule epidermis. Although actin and microtubule (MT) cytoskeletons have been implicated in cell elongation and secondary wall deposition, their roles in fiber initiation is poorly understood. Here, we used fluorescent probes and pharmacological approaches to study the roles of these cytoskeletal components during cotton fiber initiation. Both cytoskeletons align along the growth axis in initiating fibers. The dorsal view of ovules shows that unlike the fine actin filaments (AFs) in nonfiber cells, the AFs in fiber cells are dense and bundled. MTs are randomized in fiber cells and well-ordered in nonfiber cells. The pharmacological experiments revealed that the depolymerization of AFs and MTs assisted fiber initiation. Both AF stabilization and depolymerization inhibited fiber elongation. In contrast, the proper depolymerization of MTs promoted cell elongation, although the MT-stabilizing drug consistently resulted in a negative effect. Notably, we found that the organization of AFs was correlated with MT dynamics. Stabilizing the MTs by taxol treatment promoted the formation of AF bundles (in fiber initials) and transversely aligned AFs (in elongating fibers), whereas depolymerizing the MTs by oryzalin treatment promoted the fragmentation of AFs. Collectively, our data indicates that MTs plays a crucial role in regulating AF organization and early development of cotton fibers.

Dynamic roles and intricate mechanisms of ethylene in epidermal hair development in Arabidopsis and cotton

Ethylene affects many aspects of plant growth and development, including root hairs and trichomes growth in Arabidopsis, as well as fiber development in cotton, though the underlying mechanism is unclear. In this article, we update the research progress associated with the main genes in ethylene biosynthesis and signaling pathway, and we propose a clear ethylene pathway based on genome-wide identification of homologues in cotton. Expression pattern analysis using transcriptome data revealed that some candidate genes may contribute to cotton fiber development through the ethylene pathway. Moreover, we systematically summarized the effects of ethylene on the development of epidermal hair and the underlying regulatory mechanisms in Arabidopsis. Based on the knowledge of ethylene-promoted cell differentiation, elongation, and development in different tissues or plants, we advised a possible regulatory network for cotton fiber development with ethylene as the hub. Importantly, we emphasized the roles of ethylene as an important node in regulating cotton vegetative growth, and stress resistance, and suggested utilizing multiple methods to subtly modify ethylene synthesis or signaling in a tissue or spatiotemporal-specific manner to clarify its exact effect on architecture, adaptability of the plant, and fiber development, paving the way for basic research and genetic improvement of the cotton crop.

Cell twisting during desiccation reveals axial asymmetry in wall organization

Plant cell size and shape are tuned to their function and specified primarily by cellulose microfibril (CMF) patterning of the cell wall. Arabidopsis thaliana leaf trichomes are unicellular structures that act as a physical defense to deter insect feeding. This highly polarized cell type employs a strongly anisotropic cellulose wall to extend and taper, generating sharply pointed branches. During elongation, the mechanisms by which shifts in fiber orientation generate cells with predictable sizes and shapes are unknown. Specifically, the axisymmetric growth of trichome branches is often thought to result from axisymmetric CMF patterning. Here, we analyzed the direction and degree of twist of branches after desiccation to reveal the presence of an asymmetric cell wall organization with a left-hand bias. CMF organization, quantified using computational modeling, suggests a limited reorientation of microfibrils during growth and a maximum branch length limited by the wall axial stiffness. The model provides a mechanism for CMF asymmetry, which occurs after the branch bending stiffness becomes low enough that ambient bending affects the principal stresses. After this stage, the CMF synthesis results in a constant bending stiffness for longer branches. The bending vibration natural frequencies of branches with respect to their length are also discussed.

Optimizing the Protein Fluorescence Reporting System for Somatic Embryogenesis Regeneration Screening and Visual Labeling of Functional Genes in Cotton

Protein fluorescence reporting systems are of crucial importance to in-depth life science research, providing systematic labeling tools for visualization of microscopic biological activities in vivo and revolutionizing basic research. Cotton somatic cell regeneration efficiency is low, causing difficulty in cotton transformation. It is conducive to screening transgenic somatic embryo using the fluorescence reporting system. However, available fluorescence labeling systems in cotton are currently limited. To optimize the fluorescence reporting system of cotton with an expanded range of available fluorescent proteins, we selected 11 fluorescent proteins covering red, green, yellow, and cyan fluorescence colors and expressed them in cotton. Besides mRuby2 and G3GFP, the other nine fluorescent proteins (mCherry, tdTomato, sfGFP, Clover, EYFP, YPet, mVenus, mCerulean, and ECFP) were stably and intensely expressed in transgenic callus and embryo, and inherited in different cotton organs derive from the screened embryo. In addition, transgenic cotton expressing tdTomato appears pink under white light, not only for callus and embryo tissues but also various organs of mature plants, providing a visual marker in the cotton genetic transformation process, accelerating the evaluation of transgenic events. Further, we constructed transgenic cotton expressing mCherry-labeled organelle markers in vivo that cover seven specific subcellular compartments: plasma membrane, endoplasmic reticulum, tonoplast, mitochondrion, plastid, Golgi apparatus, and peroxisome. We also provide a simple and highly efficient strategy to quickly determine the subcellular localization of uncharacterized proteins in cotton cells using organelle markers. Lastly, we built the first cotton stomatal fluorescence reporting system using stomata-specific expression promoters (ProKST1, ProGbSLSP, and ProGC1) to drive Clover expression. The optimized fluorescence labeling system for transgenic somatic embryo screening and functional gene labeling in this study offers the potential to accelerating somatic cell regeneration efficiency and the in vivo monitoring of diverse cellular processes in cotton.

GhADF6-mediated actin reorganization is associated with defence against Verticillium dahliae infection in cotton

Several studies have revealed that actin depolymerizing factors (ADFs) participate in plant defence responses; however, the functional mechanisms appear intricate and need further exploration. In this study, we identified an ADF6 gene in upland cotton (designated as GhADF6) that is evidently involved in cotton's response to the fungal pathogen Verticillium dahliae. GhADF6 binds to actin filaments and possesses actin severing and depolymerizing activities in vitro and in vivo. When cotton root (the site of the fungus invasion) was inoculated with the pathogen, the expression of GhADF6 was markedly down-regulated in the epidermal cells. By virus-induced gene silencing analysis, the down-regulation of GhADF6 expression rendered the cotton plants tolerant to V. dahliae infection. Accordingly, the abundance of actin filaments and bundles in the root cells was significantly higher than that in the control plant, which phenocopied that of the V. dahliae-challenged wild-type cotton plant. Altogether, our results provide evidence that an increase in filament actin (F-actin) abundance as well as dynamic actin remodelling are required for plant defence against the invading pathogen, which are likely to be fulfilled by the coordinated expressional regulation of the actin-binding proteins, including ADF.

Elongation of wood fibers combines features of diffuse and tip growth

Xylem fibers are highly elongated cells that are key constituents of wood, play major physiological roles in plants, comprise an important terrestrial carbon reservoir, and thus have enormous ecological and economic importance. As they develop, from fusiform initials, their bodies remain the same length while their tips elongate and intrude into intercellular spaces. To elucidate mechanisms of tip elongation, we studied the cell wall along the length of isolated, elongating aspen xylem fibers and used computer simulations to predict the forces driving the intercellular space formation required for their growth. We found pectin matrix epitopes (JIM5, LM7) concentrated at the tips where cellulose microfibrils have transverse orientation, and xyloglucan epitopes (CCRC-M89, CCRC-M58) in fiber bodies where microfibrils are disordered. These features are accompanied by changes in cell wall thickness, indicating that while the cell wall elongates strictly at the tips, it is deposited all over fibers. Computer modeling revealed that the intercellular space formation needed for intrusive growth may only require targeted release of cell adhesion, which allows turgor pressure in neighboring fiber cells to 'round' the cells creating spaces. These characteristics show that xylem fibers' elongation involves a distinct mechanism that combines features of both diffuse and tip growth.

Recent Advances and Future Perspectives in Cotton Research

Cotton is not only the world's most important natural fiber crop, but it is also an ideal system in which to study genome evolution, polyploidization, and cell elongation. With the assembly of five different cotton genomes, a cotton-specific whole-genome duplication with an allopolyploidization process that combined the A- and D-genomes became evident. All existing A-genomes seemed to originate from the A₀-genome as a common ancestor, and several transposable element bursts contributed to A-genome size expansion and speciation. The ethylene production pathway is shown to regulate fiber elongation. A tip-biased diffuse growth mode and several regulatory mechanisms, including plant hormones, transcription factors, and epigenetic modifications, are involved in fiber development. Finally, we describe the involvement of the gossypol biosynthetic pathway in the manipulation of herbivorous insects, the role of GoPGF in gland formation, and host-induced gene silencing for pest and disease control. These new genes, modules, and pathways will accelerate the genetic improvement of cotton.

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.

HOMEODOMAIN PROTEIN8 mediates jasmonate-triggered trichome elongation in tomato

Plant hormones can adjust the physiology and development of plants to enhance their adaptation to biotic and abiotic challenges. Jasmonic acid (JA), one of the immunity hormones in plants, triggers genome-wide transcriptional changes in response to insect attack and wounding. Although JA is known to affect the development of trichomes, epidermal appendages that form a protective barrier against various stresses, it remains unclear how JA interacts with developmental programs that regulate trichome development. In this study, we show that JA affects trichome length in tomato by releasing the transcriptional repression mediated by Jasmonate ZIM (JAZ) proteins. We identified SlJAZ4, a negative regulator preferentially expressed in trichomes, as the critical component in JA signaling in tomato trichomes. We also identified a homeodomain-leucine zipper gene, SlHD8, as the downstream regulator of JA signaling that promotes trichome elongation. SlHD8 is also highly expressed in trichomes and physically interacts with SlJAZ4. Loss-of-function mutations in SlHD8 caused shorter trichomes, a phenotype that was only partially rescued by methyl jasmonate treatment. Our dual-luciferase and chromatin immunoprecipitation-quantitative PCR assays revealed that SlHD8 regulates trichome elongation by directly binding to the promoters of a set of cell-wall-loosening protein genes and activating their transcription. Together, our findings define SlHD8-SlJAZ4 as a key module mediating JA-induced trichome elongation in tomato.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development

BACKGROUND

Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored.

RESULTS

Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes.

CONCLUSION

Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.

A Modified Actin (Gly65Val Substitution) Expressed in Cotton Disrupts Polymerization of Actin Filaments Leading to the Phenotype of Ligon Lintless-1 (Li1) Mutant

Dynamic remodeling of the actin cytoskeleton plays a central role in the elongation of cotton fibers, which are the most important natural fibers in the global textile industry. Here, a high-resolution mapping approach combined with comparative sequencing and a transgenic method revealed that a G65V substitution in the cotton actin Gh_D04G0865 (GhACT17D in the wild-type) is responsible for the Gossypium hirsutum Ligon lintless-1 (Li₁) mutant (GhACT17DM). In the mutant GhACT17DM from Li₁ plant, Gly65 is substituted with valine on the lip of the nucleotide-binding domain of GhACT17D, which probably affects the polymerization of F-actin. Over-expression of GhACT17DM, but not GhACT17D, driven by either a CaMV35 promoter or a fiber-specific promoter in cotton produced a Li₁-like phenotype. Compared with the wild-type control, actin filaments in Li₁ fibers showed higher growth and shrinkage rates, decreased filament skewness and parallelness, and increased filament density. Taken together, our results indicate that the incorporation of GhACT17DM during actin polymerization disrupts the establishment and dynamics of the actin cytoskeleton, resulting in defective fiber elongation and the overall dwarf and twisted phenotype of the Li₁ mutant.

Microtubules exert early, partial, and variable control of cotton fiber diameter

Variable cotton fiber diameter is set early in anisotropic elongation by cell-type-specific processes involving the temporal and spatial regulation of microtubules in the apical region. Cotton fibers are single cells that originate from the seed epidermis of Gossypium species. Then, they undergo extreme anisotropic elongation and limited diametric expansion. The details of cellular morphogenesis determine the quality traits that affect fiber uses and value, such as length, strength, and diameter. Lower and more consistent diameter would increase the competitiveness of cotton fiber with synthetic fiber, but we do not know how this trait is controlled. The complexity of the question is indicated by the existence of fibers in two major width classes in the major commercial species: broad and narrow fibers exist in commonly grown G. hirsutum, whereas G. barbadense produces only narrow fiber. To further understand how fiber diameter is controlled, we used ovule cultures, morphology measurements, and microtubule immunofluorescence to observe the effects of microtubule antagonists on fiber morphology, including shape and diameter within 80 µm of the apex. The treatments were applied at either one or two days post-anthesis during different stages of fiber morphogenesis. The results showed that inhibiting the presence and/or dynamic activity of microtubules caused larger diameter tips to form, with greater effects often observed with earlier treatment. The presence and geometry of a microtubule-depleted-zone below the apex were transiently correlated with the apical diameter of the narrow tip types. Similarly, the microtubule antagonists had somewhat different effects between tip types. Overall, the results demonstrate cell-type-specific mechanisms regulating fiber expansion within 80 µm of the apex, with variation in the impact of microtubules between tip types and over developmental time.

Genome-Wide Association Analysis Reveals Loci and Candidate Genes Involved in Fiber Quality Traits Under Multiple Field Environments in Cotton (Gossypium hirsutum)

Fiber length, fiber strength, and fiber micronaire are the main fiber quality parameters in cotton. Thus, mining the elite and stable loci/alleles related to fiber quality traits and elucidating the relationship between the two may accelerate genetic improvement of fiber quality in cotton. Here, genome-wide association analysis (GWAS) was performed for fiber quality parameters based on phenotypic data, and 56,010 high-quality single nucleotide polymorphisms (SNPs) using 242 upland cotton accessions under 12 field environments were obtained. Phenotypic analysis exhibited that fiber length (FL) had a positive correlation with fiber strength (FS) and had a negative correlation with fiber micronaire (Mic). Genetic analysis also indicated that FL, FS, and Mic had high heritability of more than 80%. A total of 67 stable quantitative trait loci (QTLs) were identified through GWAS analysis, including 31 for FL, 21 for FS, and 22 for Mic. Of them, three pairs homologous QTLs were detected between A and D subgenomes, and seven co-located QTLs with two fiber quality parameters were found. Compared with the reported QTLs, 34 co-located with previous studies, and 33 were newly revealed. Integrated with transcriptome analysis, we selected 256, 244, and 149 candidate genes for FL, FS, and Mic, respectively. Gene Ontology (GO) analysis showed that most of the genes located in QTLs interval of the three fiber quality traits were involved in sugar biosynthesis, sugar metabolism, microtubule, and cytoskeleton organization, which played crucial roles in fiber development. Through correlation analysis between haplotypes and phenotypes, three genes (GH_A05G1494, GH_D11G3097, and GH_A05G1082) predominately expressed in fiber development stages were indicated to be potentially responsible for FL, FS, and Mic, respectively. The GH_A05G1494 encoded a protein containing SGS-domain, which is related to tubulin-binding and ubiquitin-protein ligase binding. The GH_D11G3097 encoded 20S proteasome beta subunit G1, and was involved in the ubiquitin-dependent protein catabolic process. The GH_A05G1082 encoded RAN binding protein 1 with a molecular function of GTPase activator activity. These results provide new insights and candidate loci/genes for the improvement of fiber quality in cotton.

Genome-Wide Identification and Comparative Analysis of Myosin Gene Family in Four Major Cotton Species

Myosin protein as a molecular motor, binding with Actin, plays a significant role in various physiological activities such as cell division, movement, migration, and morphology; however, there are only a few studies on plant Myosin gene family, particularly in cotton. A total of 114 Myosin genes were found in Gossypium hirsutum, Gossypium barbadense, Gossypium raimondii, and Gossypium arboreum. All Myosins could be grouped into six groups, and for each group of these genes, similar gene structures are found. Study of evolution suggested that the whole genome duplications event occurring about 13-20 MYA (millions of years ago) is the key explanation for Myosins expanse in cotton. Cis-element and qPCR analysis revealed that plant hormones such as abscisic acid, methyl jasmonate, and salicylic acid can control the expression of Myosins. This research provides useful information on the function of Myosin genes in regulating plant growth, production, and fiber elongation for further studies.

Gossypium Genomics: Trends, Scope, and Utilization for Cotton Improvement

Cotton (Gossypium spp.) is the most important natural fiber crop worldwide. The diversity of Gossypium species also provides an ideal model for investigating evolution and domestication of polyploids. However, the huge and complex cotton genome hinders genomic research. Technical advances in high-throughput sequencing and bioinformatics analysis have now largely overcome these obstacles, bringing about a new era of cotton genomics. Here, we review recent progress in Gossypium genomics based on whole genome sequencing, resequencing, and comparative genomics, which have provided insights about the genomic basis of fiber biogenesis and the landscape of cotton functional genomics. We address current challenges and present multidisciplinary genomics-enabled breeding strategies covering the breadth of high fiber yield, quality, and environmental resilience for future cotton breeding programs.

CytoSeg 2.0: automated extraction of actin filaments

MOTIVATION

Actin filaments (AFs) are dynamic structures that substantially change their organization over time. The dynamic behavior and the relatively low signal-to-noise ratio during live-cell imaging have rendered the quantification of the actin organization a difficult task.

RESULTS

We developed an automated image-based framework that extracts AFs from fluorescence microscopy images and represents them as networks, which are automatically analyzed to identify and compare biologically relevant features. Although the source code is freely available, we have now implemented the framework into a graphical user interface that can be installed as a Fiji plugin, thus enabling easy access by the research community.

AVAILABILITY AND IMPLEMENTATION

CytoSeg 2.0 is open-source software under the GPL and is available on Github: https://github.com/jnowak90/CytoSeg2.0.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Magnaporthe oryzae fimbrin organizes actin networks in the hyphal tip during polar growth and pathogenesis

Magnaporthe oryzae causes rice blast disease, but little is known about the dynamic restructuring of the actin cytoskeleton during its polarized tip growth and pathogenesis. Here, we used super-resolution live-cell imaging to investigate the dynamic organization of the actin cytoskeleton in M. oryzae during hyphal tip growth and pathogenesis. We observed a dense actin network at the apical region of the hyphae and actin filaments originating from the Spitzenkörper (Spk, the organizing center for hyphal growth and development) that formed branched actin bundles radiating to the cell membrane. The actin cross-linking protein Fimbrin (MoFim1) helps organize this actin distribution. MoFim1 localizes to the actin at the subapical collar, the actin bundles, and actin at the Spk. Knockout of MoFim1 resulted in impaired Spk maintenance and reduced actin bundle formation, preventing polar growth, vesicle transport, and the expansion of hyphae in plant cells. Finally, transgenic rice (Oryza sativa) expressing RNA hairpins targeting MoFim1 exhibited improved resistance to M. oryzae infection, indicating that MoFim1 represents an excellent candidate for M. oryzae control. These results reveal the dynamics of actin assembly in M. oryzae during hyphal tip development and pathogenesis, and they suggest a mechanism in which MoFim1 organizes such actin networks.

Cortical Microtubule Organization during Petal Morphogenesis in Arabidopsis

Cortical microtubules guide the direction and deposition of cellulose microfibrils to build the cell wall, which in turn influences cell expansion and plant morphogenesis. In the model plant Arabidopsis thaliana (Arabidopsis), petal is a relatively simple organ that contains distinct epidermal cells, such as specialized conical cells in the adaxial epidermis and relatively flat cells with several lobes in the abaxial epidermis. In the past two decades, the Arabidopsis petal has become a model experimental system for studying cell expansion and organ morphogenesis, because petals are dispensable for plant growth and reproduction. Recent advances have expanded the role of microtubule organization in modulating petal anisotropic shape formation and conical cell shaping during petal morphogenesis. Here, we summarize recent studies showing that in Arabidopsis, several genes, such as SPIKE1, Rho of plant (ROP) GTPases, and IPGA1, play critical roles in microtubule organization and cell expansion in the abaxial epidermis during petal morphogenesis. Moreover, we summarize the live-confocal imaging studies of Arabidopsis conical cells in the adaxial epidermis, which have emerged as a new cellular model. We discuss the microtubule organization pattern during conical cell shaping. Finally, we propose future directions regarding the study of petal morphogenesis and conical cell shaping.

Spatiotemporal cytoskeleton organizations determine morphogenesis of multicellular trichomes in tomato

Plant trichomes originate from epidermal cell, forming protective structure from abiotic and biotic stresses. Different from the unicellular trichome in Arabidopsis, tomato trichomes are multicellular structure and can be classified into seven different types based on cell number, shape and the presence of glandular cells. Despite the importance of tomato trichomes in insect resistance, our understanding of the tomato trichome morphogenesis remains elusive. In this study, we quantitatively analyzed morphological traits of trichomes in tomato and further performed live imaging of cytoskeletons in stably transformed lines with actin and microtubule markers. At different developmental stages, two types of cytoskeletons exhibited distinct patterns in different trichome cells, ranging from transverse, spiral to longitudinal. This gradual transition of actin filament angle from basal to top cells could correlate with the spatial expansion mode in different cells. Further genetic screen for aberrant trichome morphology led to the discovery of a number of independent mutations in SCAR/WAVE and ARP2/3 complex, which resulted in actin bundling and distorted trichomes. Disruption of microtubules caused isotropic expansion while abolished actin filaments entirely inhibited axial extension of trichomes, indicating that microtubules and actin filaments may control distinct aspects of trichome cell expansion. Our results shed light on the roles of cytoskeletons in the formation of multicellular structure of tomato trichomes.