Genome-wide identification and analysis of the evolution and expression patterns of the cellulose synthase gene superfamily in Gossypium species

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

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.

Identification and bioinformatic analysis of the CaCesA/Csls family members and the expression of the CaCslD1 in the flower buds of CMS/Rf system in pepper

The cellulose synthase gene superfamily contains cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, which synthesize cellulose and hemicellulose in plant cell walls and play a crucial role in plant growth and development. However, the CesA/Csl gene family has not been reported in pepper. Therefore, the genome-wide research of the CaCesA/CaCsl gene family was conducted in pepper. In this study, a total of 39 CaCesA/CaCsls genes (10 CesAs genes and 29 Csls genes) were identified in pepper and unevenly distributed on 11 chromosomes. These CaCesA/Csls were divided into seven subfamilies (CesAs, CslAs, CslBs, CslCs, CslDs, CslEs, CslGs), and most of CaCesA/Csls genes are closely related to AtCesA/Csls genes. The cis-acting elements in the promoters of CaCesA/Csls genes are mainly related to hormone response and stress response. There are ten collinear gene pairs between the CesA/Csls gene family of pepper and Arabidopsis, and four fragment duplication gene pairs of the CaCesA/Csls genes were discovered. RNA-seq analysis shows that the majority of CaCesA/Csls are expressed in a variety of plant tissues, indicating that most CaCesA/Csls gene expression patterns are not organ-specific, and CaCslD1/D4 have the highest expression in anthers, followed by petal, ovary, and F9. RNA-seq analysis shows that most CaCesA/Csls are responsive to five hormones (IAA, GA₃, ABA, SA, and MeJA). The tissue-specific expression analysis of the CaCslD1 gene shows that the CaCslD1 gene is expressed specifically in flowers. In the flower buds IV of cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system, CaCslD1 reach the highest expression respectively. However, the relative expression level of CaCslD1 in the fertile accessions is extremely significantly higher than in the sterile accessions. This study shows an overall understanding of the CaCesA/Csls gene family and provides a new insight for understanding the function of CaCslD1 in pollen development and exploring the fertility restoration of CMS in pepper.

Species-Specific Gene Expansion of the Cellulose synthase Gene Superfamily in the Orchidaceae Family and Functional Divergence of Mannan Synthesis-Related Genes in Dendrobium officinale

Plant Cellulose synthase genes constitute a supergene family that includes the Cellulose synthase (CesA) family and nine Cellulose synthase-like (Csl) families, the members of which are widely involved in the biosynthesis of cellulose and hemicellulose. However, little is known about the Cellulose synthase superfamily in the family Orchidaceae, one of the largest families of angiosperms. In the present study, we identified and systematically analyzed the CesA/Csl family members in three fully sequenced Orchidaceae species, i.e., Dendrobium officinale, Phalaenopsis equestris, and Apostasia shenzhenica. A total of 125 Cellulose synthase superfamily genes were identified in the three orchid species and classified into one CesA family and six Csl families: CslA, CslC, CslD, CslE, CslG, and CslH according to phylogenetic analysis involving nine representative plant species. We found species-specific expansion of certain gene families, such as the CslAs in D. officinale (19 members). The CesA/Csl families exhibited sequence divergence and conservation in terms of gene structure, phylogeny, and deduced protein sequence, indicating multiple origins via different evolutionary processes. The distribution of the DofCesA/DofCsl genes was investigated, and 14 tandemly duplicated genes were detected, implying that the expansion of DofCesA/DofCsl genes may have originated via gene duplication. Furthermore, the expression profiles of the DofCesA/DofCsl genes were investigated using transcriptome sequencing and quantitative Real-time PCR (qRT-PCR) analysis, which revealed functional divergence in different tissues and during different developmental stages of D. officinale. Three DofCesAs were highly expressed in the flower, whereas DofCslD and DofCslC family genes exhibited low expression levels in all tissues and at all developmental stages. The 19 DofCslAs were differentially expressed in the D. officinale stems at different developmental stages, among which six DofCslAs were expressed at low levels or not at all. Notably, two DofCslAs (DofCslA14 and DofCslA15) showed significantly high expression in the stems of D. officinale, indicating a vital role in mannan synthesis. These results indicate the functional redundancy and specialization of DofCslAs with respect to polysaccharide accumulation. In conclusion, our results provide insights into the evolution, structure, and expression patterns of CesA/Csl genes and provide a foundation for further gene functional analysis in Orchidaceae and other plant species.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

Cellulose synthase-like protein OsCSLD4 plays an important role in the response of rice to salt stress by mediating abscisic acid biosynthesis to regulate osmotic stress tolerance

Cell wall polysaccharide biosynthesis enzymes play important roles in plant growth, development and stress responses. The functions of cell wall polysaccharide synthesis enzymes in plant growth and development have been well studied. In contrast, their roles in plant responses to environmental stress are poorly understood. Previous studies have demonstrated that the rice cell wall cellulose synthase-like D4 protein (OsCSLD4) is involved in cell wall polysaccharide synthesis and is important for rice growth and development. This study demonstrated that the OsCSLD4 function-disrupted mutant nd1 was sensitive to salt stress, but insensitive to abscisic acid (ABA). The expression of some ABA synthesis and response genes was repressed in nd1 under both normal and salt stress conditions. Exogenous ABA can restore nd1-impaired salt stress tolerance. Moreover, overexpression of OsCSLD4 can enhance rice ABA synthesis gene expression, increase ABA content and improve rice salt tolerance, thus implying that OsCSLD4-regulated rice salt stress tolerance is mediated by ABA synthesis. Additionally, nd1 decreased rice tolerance to osmotic stress, but not ion toxic tolerance. The results from the transcriptome analysis showed that more osmotic stress-responsive genes were impaired in nd1 than salt stress-responsive genes, thus indicating that OsCSLD4 is involved in rice salt stress response through an ABA-induced osmotic response pathway. Intriguingly, the disruption of OsCSLD4 function decreased grain width and weight, while overexpression of OsCSLD4 increased grain width and weight. Taken together, this study demonstrates a novel plant salt stress adaptation mechanism by which crops can coordinate salt stress tolerance and yield.

Genome-wide bioinformatics analysis of Cellulose Synthase gene family in common bean (Phaseolus vulgaris L.) and the expression in the pod development

Background

CesA and Csl gene families, which belong to the cellulose synthase gene superfamily, plays an important role in the biosynthesis of the plant cell wall. Although researchers have investigated this gene superfamily in several model plants, to date, no comprehensive analysis has been conducted in the common bean.

Results

In this study, we identified 39 putative cellulose synthase genes from the common bean genome sequence. Then, we performed a bioinformatics analysis of this gene family involving sequence alignment, phylogenetic analysis, gene structure, collinearity analysis and chromosome location. We found all members possess a cellulose_synt domain. Phylogenetic analysis revealed that these cellulose synthase genes may be classified into five subfamilies, and that members in the same subfamily share conserved exon-intron distribution and motif compositions. Abundant and distinct cis-acting elements in the 2 k basepairs upstream regulatory regions indicate that the cellulose synthase gene family may plays a vital role in the growth and development of common bean. Moreover, the 39 cellulose synthase genes are distributed on 10 of the 11 chromosomes. Additionally expression analysis shows that all CesA/Csl genes selected are constitutively expressed in the pod development.

Conclusions

This research reveals both the putative biochemical and physiological functions of cellulose synthase genes in common bean and implies the importance of studying non-model plants to understand the breadth and diversity of cellulose synthase genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01026-0.

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

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

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log₁₀(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

Phenotypic effects of changes in the FTVTxK region of an Arabidopsis secondary wall cellulose synthase compared with results from analogous mutations in other isoforms

Understanding protein structure and function relationships in cellulose synthase (CesA), including divergent isomers, is an important goal. Here, we report results from mutant complementation assays that tested the ability of sequence variants of AtCesA7, a secondary wall CesA of Arabidopsis thaliana, to rescue the collapsed vessels, short stems, and low cellulose content of the irx3-1 AtCesA7 null mutant. We tested a catalytic null mutation and seven missense or small domain changes in and near the AtCesA7 FTVTSK motif, which lies near the catalytic domain and may, analogously to bacterial CesA, exist within a substrate "gating loop." A low-to-high gradient of rescue occurred, and even inactive AtCesA7 had a small positive effect on stem cellulose content but not stem elongation. Overall, secondary wall cellulose content and stem length were moderately correlated, but the results were consistent with threshold amounts of cellulose supporting particular developmental processes. Vibrational sum frequency generation microscopy allowed tissue-specific analysis of cellulose content in stem xylem and interfascicular fibers, revealing subtle differences between selected genotypes that correlated with the extent of rescue of the collapsing xylem phenotype. Similar tests on PpCesA5 from the moss Physcomitrium (formerly Physcomitrella) patens helped us to synergize the AtCesA7 results with prior results on AtCesA1 and PpCesA5. The cumulative results show that the FTVTxK region is important for the function of an angiosperm secondary wall CesA as well as widely divergent primary wall CesAs, while differences in complementation results between isomers may reflect functional differences that can be explored in further work.

Advances about the Roles of Membranes in Cotton Fiber Development

Cotton fiber is an extremely elongated single cell derived from the ovule epidermis and is an ideal model for studying cell development. The plasma membrane is tremendously expanded and accompanied by the coordination of various physiological and biochemical activities on the membrane, one of the three major systems of a eukaryotic cell. This review compiles the recent progress and advances for the roles of the membrane in cotton fiber development: the functions of membrane lipids, especially the fatty acids, sphingolipids, and phytosterols; membrane channels, including aquaporins, the ATP-binding cassette (ABC) transporters, vacuolar invertase, and plasmodesmata; and the regulation mechanism of membrane proteins, such as membrane binding enzymes, annexins, and receptor-like kinases.

The cellulose synthase (CesA) gene family in four Gossypium species: phylogenetics, sequence variation and gene expression in relation to fiber quality in Upland cotton

Cellulose synthases (CesAs) are multi-subunit enzymes found on the plasma membrane of plant cells and play a pivotal role in cellulose production. The cotton fiber is mainly composed of cellulose, and the genetic relationships between CesA genes and cotton fiber yield and quality are not fully understood. Through a phylogenetic analysis, the CesA gene family in diploid Gossypium arboreum and Gossypium raimondii, as well as tetraploid Gossypium hirsutum ('TM-1') and Gossypium barbadense ('Hai-7124' and '3-79'), was divided into 6 groups and 15 sub-groups, with each group containing two to five homologous genes. Most CesA genes in the four species are highly collinear. Among the five cotton genomes, 440 and 1929 single nucleotide polymorphisms (SNPs) in the CesA gene family were identified in exons and introns, respectively, including 174 SNPs resulting in amino acid changes. In total, 484 homeologous SNPs between the A and D genomes were identified in diploids, while 142 SNPs were detected between the two tetraploids, with 32 and 82 SNPs existing within G. hirsutum and G. barbadense, respectively. Additionally, 74 quantitative trait loci near 18 GhCesA genes were associated with fiber quality. One to four GhCesA genes were differentially expressed (DE) in ovules at 0 and 3 days post anthesis (DPA) between two backcross inbred lines having different fiber lengths, but no DE genes were identified between these lines in developing fibers at 10 DPA. Twenty-seven SNPs in above DE CesA genes were detected among seven cotton lines, including one SNP in Ghi_A08G03061 that was detected in four G. hirsutum genotypes. This study provides the first comprehensive characterization of the cotton CesA gene family, which may play important roles in determining cotton fiber quality.

Genome-wide characterization of the cellulose synthase gene superfamily in Pyrus bretschneideri and reveal its potential role in stone cell formation

Members of the cellulose synthase (CesA) and cellulose synthase-like (Csl) families from the cellulose synthase gene superfamily participate in cellulose and hemicellulose synthesis in the plasma membrane. The members of this superfamily are vital for cell wall construction during plant growth and development. However, little is known about their function in pear fruit, a model for Rosaceae species and for fleshy fruit development. In our research, a total of 36 CesA/Csl family members were identified from the pear and were grouped into six subfamilies (CesA, CslB, CslC, CslD, CslE, and CslG) according to phylogenetic relationships. We performed a protein sequence physicochemical analysis, phylogenetic tree construction, a gene structure, a conserved domain, and chromosomal localization analysis. The results indicated that most of the CesA/Csl genes from pear are closely related to genes in Arabidopsis, but these families have unique characteristics in terms of their gene structure, chromosomal localization, phylogeny, and deduced protein sequences, suggesting that they have evolved through different processes. Tissue expression analysis results showed that most of the CesA/Csl genes were constitutively expressed at different levels in different organs. Furthermore, the expression levels of four genes (Pbr032894.2, Pbr016107.1, Pbr00518.1, and Pbr034218.1) tended to first increase and then decrease during fruit development, implying that these four genes may be involved in the development of stone cells of pear fruit. Our results may help elucidate the evolutionary history and functional differences of the CesA/Csl genes in pear and lay a foundation for further investigation of the CesA/Csl genes in pear and other Rosaceae species.

Insight into the relationship between S-lignin and fiber quality based on multiple research methods

Cotton (Gossypium hirsutum) is an important cash crop, providing people with high quality natural fiber. Lignin is the main component of cotton fiber, second only to cellulose. As a main substance filled in the cellulose framework during the secondary wall thickening process, lignin plays a key role in the formation of cotton fiber quality. However, the mechanism behind it is still unclear. In this research, we screened candidate genes involved in lignin biosynthesis based on analysis of cotton genome and transcriptome sequence data. The authenticity of the transcriptome data was verified by qRT-PCR assay. Total 62 genes were identified from nine gene families. In the process, we found the key gene GhCAD7 that affects the biosynthesis of S-lignin and the ratio of syringyl/guaiacyl (S/G). In addition, in combination with the metabolites and transcriptome profiles of the line 0-153 with high fiber quality and the line sGK9708 with low fiber quality during cotton fiber development, we speculate that the ratio of syringyl/guaiacyl (S/G) is inseparable from the quality of cotton fiber. Finally, the S-type lignin synthesis branch may play a more important role in the formation of high-quality fiber. This work provides insights into the synthesis of lignin in cotton and lays the foundation for future research into improving fiber quality.

Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population

Cotton is widely cultivated globally because it provides natural fibre for the textile industry and human use. To identify quantitative trait loci (QTLs)/genes associated with fibre quality and yield, a recombinant inbred line (RIL) population was developed in upland cotton. A consensus map covering the whole genome was constructed with three types of markers (8295 markers, 5197.17 centimorgans (cM)). Six fibre yield and quality traits were evaluated in 17 environments, and 983 QTLs were identified, 198 of which were stable and mainly distributed on chromosomes 4, 6, 7, 13, 21 and 25. Thirty-seven QTL clusters were identified, in which 92.8% of paired traits with significant medium or high positive correlations had the same QTL additive effect directions, and all of the paired traits with significant medium or high negative correlations had opposite additive effect directions. In total, 1297 genes were discovered in the QTL clusters, 414 of which were expressed in two RNA-Seq data sets. Many genes were discovered, 23 of which were promising candidates. Six important QTL clusters that included both fibre quality and yield traits were identified with opposite additive effect directions, and those on chromosome 13 (qClu-chr13-2) could increase fibre quality but reduce yield; this result was validated in a natural population using three markers. These data could provide information about the genetic basis of cotton fibre quality and yield and help cotton breeders to improve fibre quality and yield simultaneously.

Cotton CSLD3 restores cell elongation and cell wall integrity mainly by enhancing primary cellulose production in the Arabidopsis cesa6 mutant

Overexpression of cotton cellulose synthase like D3 (GhCSLD3) gene partially rescued growth defect of atcesa6 mutant with restored cell elongation and cell wall integrity mainly by enhancing primary cellulose production. Among cellulose synthase like (CSL) family proteins, CSLDs share the highest sequence similarity to cellulose synthase (CESA) proteins. Although CSLD proteins have been implicated to participate in the synthesis of carbohydrate-based polymers (cellulose, pectins and hemicelluloses), and therefore plant cell wall formation, the exact biochemical function of CSLD proteins remains controversial and the function of the remaining CSLD genes in other species have not been determined. In this study, we attempted to illustrate the function of CSLD proteins by overexpressing Arabidopsis AtCSLD2, -3, -5 and cotton GhCSLD3 genes in the atcesa6 mutant, which has a background that is defective for primary cell wall cellulose synthesis in Arabidopsis. We found that GhCSLD3 overexpression partially rescued the growth defect of the atcesa6 mutant during early vegetative growth. Despite the atceas6 mutant having significantly reduced cellulose contents, the defected cell walls and lower dry mass, GhCSLD3 overexpression largely restored cell wall integrity (CWI) and improved the biomass yield. Our result suggests that overexpression of the GhCSLD protein enhances primary cell wall synthesis and compensates for the loss of CESAs, which is required for cellulose production, therefore rescuing defects in cell elongation and CWI.

Identification of TPX2 Gene Family in Upland Cotton and Its Functional Analysis in Cotton Fiber Development

Microtubules (MTs) are of importance to fiber development. The Xklp2 (TPX2) proteins as a class of microtubule-associated proteins (MAPs) play a key role in plant growth and development by regulating the dynamic changes of microtubules (MTs). However, the mechanism underlying this is unknown. The interactions between TPX2 proteins and tubulin protein, which are the main structural components, have not been studied in fiber development of upland cotton. Therefore, a genome-wide analysis of the TPX2 family was firstly performed in Gossypiumhirsutum L. This study identified 41 GhTPX2 sequences in the assembled G. hirsutum genome by a series of bioinformatic methods. Generally, this gene family is phylogenetically grouped into six subfamilies, and 41 G. hirsutum TPX2 genes (GhTPX2s) are distributed across 21 chromosomes. A heatmap of the TPX2 gene family showed that homologous GhTPX2 genes, GhWDLA2/7 and GhWDLA4/9, have large differences in expression levels between two upland cotton recombinant inbred lines (69307 and 69362) that are different in fiber quality at 15 and 20 days post anthesis. The relative data indicate that these four genes are down-regulated under oryzalin, which causes microtubule depolymerization, as determined via qRT-PCR. A subcellular localization experiment suggested that GhWDLA2 and GhWDLA7 are localized to the microtubule cytoskeleton, and GhWDLA4 and GhWDLA9 are only localized to the nucleus. However, only GhWDLA7 between GhWDLA2 and GhWDLA7 interacted with GhTUA2 in the yeast two-hybrid assay. These results lay the foundation for further function study of the TPX2 gene family.

Genome-wide characterization of the cellulose synthase gene superfamily in Solanum lycopersicum

The cellulose synthase gene superfamily, which includes the cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, plays a vital role in the biosynthesis of cellulose and hemicellulose in plants. However, these genes have not been extensively studied in tomato (Solanum lycopersicum), a model for Solanaceae plants and for fleshy fruit development. Here, we identified and systematically analyzed 38 CesA/Csl family members that contained cellulose synthase domain regions, and categorized their encoded proteins into 6 subfamilies (CesA, CslA, CslB, CslD, CslE, and CslG) based on phylogenetic analysis. Most CesA/Csl genes from tomato are closely related to those from Arabidopsis, but the families have distinct features regarding gene structure, chromosome distribution and localization, phylogeny, and deduced protein sequence, indicating that they arose via different evolutionary process. Furthermore, expression analysis of CesA/Csl genes in different tissues at various developmental stages showed that most CesAs were constitutively expressed with differential expression levels in various organs; three CslD genes were expressed specifically in flowers, and four CesA and five Csl putative genes were preferentially expressed in fruits. Our results provide insight into the general characteristics of the CesA/Csl genes in tomato, and lay the foundation for further functional studies of CesA/Csl genes in tomato and other Solanaceae species.