Transcriptome profiling, molecular biological, and physiological studies reveal a major role for ethylene in cotton fiber cell elongation

Potential Roles of 1-Aminocyclopropane-1-carboxylic Acid Synthase Genes in the Response of Gossypium Species to Abiotic Stress by Genome-Wide Identification and Expression Analysis

Ethylene plays a pivotal role in plant stress resistance and 1-aminocyclopropane-1-carboxylic acid synthase (ACS) is the rate-limiting enzyme in ethylene biosynthesis. Upland cotton (Gossypium hirsutum L.) is the most important natural fiber crop, but the function of ACS in response to abiotic stress has rarely been reported in this plant. We identified 18 GaACS, 18 GrACS, and 35 GhACS genes in Gossypiumarboreum, Gossypium raimondii and Gossypiumhirsutum, respectively, that were classified as types I, II, III, or IV. Collinearity analysis showed that the GhACS genes were expanded from diploid cotton by the whole-genome-duplication. Multiple alignments showed that the C-terminals of the GhACS proteins were conserved, whereas the N-terminals of GhACS10 and GhACS12 were different from the N-terminals of AtACS10 and AtACS12, probably diverging during evolution. Most type II ACS genes were hardly expressed, whereas GhACS10/GhACS12 were expressed in many tissues and in response to abiotic stress; for example, they were highly and hardly expressed at the early stages of cold and heat exposure, respectively. The GhACS genes showed different expression profiles in response to cold, heat, drought, and salt stress by quantitative PCR analysis, which indicate the potential roles of them when encountering the various adverse conditions, and provide insights into GhACS functions in cotton’s adaptation to abiotic stress.

The bHLH/HLH transcription factors GhFP2 and GhACE1 antagonistically regulate fiber elongation in cotton

Basic helix-loop-helix/helix-loop-helix (bHLH/HLH) transcription factors play important roles in cell elongation in plants. However, little is known about how bHLH/HLH transcription factors antagonistically regulate fiber elongation in cotton (Gossypium hirsutum). In this study, we report that two bHLH/HLH transcription factors, fiber-related protein 2 (GhFP2) and ACTIVATOR FOR CELL ELONGATION 1 (GhACE1), function in fiber development of cotton. GhFP2 is an atypical bHLH protein without the basic region, and its expression is regulated by brassinosteroid (BR)-related BRASSINAZOLE RESISTANT 1 (GhBZR1). Overexpression of GhFP2 in cotton hindered fiber elongation, resulting in shortened fiber length. In contrast, suppression of GhFP2 expression in cotton promoted fiber development, leading to longer fibers compared with the wild-type. GhFP2 neither contains a DNA-binding domain nor has transcriptional activation activity. Furthermore, we identified GhACE1, a bHLH protein that interacts with GhFP2 and positively regulates fiber elongation. GhACE1 could bind to promoters of plasma membrane intrinsic protein 2;7 (GhPIP2;7) and expansions 8 (GhEXP8) for directly activating their expression, but the interaction between GhFP2 and GhACE1 suppressed transcriptional activation of these target genes by GhACE1. Taken together, our results indicate that GhACE1 promotes fiber elongation by activating expressions of GhPIP2;7 and GhEXP8, but its transcription activation on downstream genes may be obstructed by BR-modulated GhFP2. Thus, our data reveal a key mechanism for fiber cell elongation through a pair of antagonizing HLH/bHLH transcription factors in cotton.

In silico analysis and expression profiling of Expansin A4, BURP domain protein RD22-like and E6-like genes associated with fiber quality in cotton

BACKGROUND

To supply high-quality cotton fibre for the textile industry, the development of long, strong and fine fibre cotton varieties is imperative. An interlinked approach was used to comprehend the role of fibre genes by analyzing interspecific progenies of cotton species. Wild Gossypium species and races are rich source of genetic polymorphism due to environmental dispersal and continuous natural selection. These genetic resources hold mass of outclass genes that can be used in cotton improvement breeding programs to exploit possible traits such as fibre quality, abiotic stress tolerance, and disease and insect resistance. Therefore, use of new molecular techniques such as genomics, transcriptomics and bioinformatics is very important to utilize the genetic potential of wild species in cotton improvement programs.

METHODS

Interspecific lines and Gossypium species used in the study were grown at Central Cotton Research Institute (CCRI), Multan. After retrieving DNA sequence of the genes from NCBI, the primers for gene expression and full-length gene sequence were designed. Expression profiling of Expansin A4, BURP Domain protein RD22-like and E6-like fibre genes was performed through Real Time PCR. BLAST and DNA sequence alignment was conducted for sequence comparison of interspecific lines and Gossypium species. Different in silico analysis were used for characterization of fibre genes and identification of cis acting promoter elements in promoter region.

RESULTS

Variable expression of genes related to fibre development was observed at different stages. BLAST and DNA sequence alignment demonstrated resemblance of interspecific lines with G. hirsutum. In silico analysis on the sequence data also confirmed the role of Expansin A4, BURP Domain protein RD22-like and E6-like fibre genes in fibre development. Genetic engineering is also recommended by transferring E6-like, Expansin A4 and BURP Domain RD22-like genes in local cotton cultivars. Similarly, several stress tolerant and light responsive cis acting elements were identified through promotor analysis, which may contribute for fibre development in the breeding programs.

CONCLUSION

Expansin A4, BURP Domain RD22-like and E6-like have positive role in fibre development with variable expression at fiber length and strength associated stages.

Molecular Regulation of Cotton Fiber Development: A Review

Cotton (Gossypium spp.) is an economically important natural fiber crop. The quality of cotton fiber has a substantial effect on the quality of cotton textiles. The identification of cotton fiber development-related genes and exploration of their biological functions will not only enhance our understanding of the elongation and developmental mechanisms of cotton fibers but also provide insights that could aid the cultivation of new cotton varieties with improved fiber quality. Cotton fibers are single cells that have been differentiated from the ovule epidermis and serve as a model system for research on single-cell differentiation, growth, and fiber production. Genes and fiber formation mechanisms are examined in this review to shed new light on how important phytohormones, transcription factors, proteins, and genes linked to fiber development work together. Plant hormones, which occur in low quantities, play a critically important role in regulating cotton fiber development. Here, we review recent research that has greatly contributed to our understanding of the roles of different phytohormones in fiber development and regulation. We discuss the mechanisms by which phytohormones regulate the initiation and elongation of fiber cells in cotton, as well as the identification of genes involved in hormone biosynthetic and signaling pathways that regulate the initiation, elongation, and development of cotton fibers.

Transcriptome Time-Course Analysis in the Whole Period of Cotton Fiber Development

Gossypium hirsutum and Gossypium barbadense are the widely cultivated tetraploid cottons around the world, which evolved great differences in the fiber yield and quality due to the independent domestication process. To reveal the genetic basis of the difference, we integrated 90 samples from ten time points during the fiber developmental period for investigating the dynamics of gene expression changes associated with fiber in G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124 and acc. 3-79. Globally, 44,484 genes expressed in all three cultivars account for 61.14% of the total genes. About 61.39% (N = 3,412) of the cotton transcription factors were involved in fiber development, which consisted of 58 cotton TF families. The differential analysis of intra- and interspecies showed that 3 DPA had more expression changes. To discover the genes with temporally changed expression profiles during the whole fiber development, 1,850 genes predominantly expressed in G. hirsutum and 1,050 in G. barbadense were identified, respectively. Based on the weighted gene co-expression network and time-course analysis, several candidate genes, mainly involved in the secondary cell wall synthesis and phytohormones, were identified in this study, underlying possibly the transcriptional regulation and molecular mechanisms of the fiber quality differences between G. barbadense and G. hirsutum. The quantitative real-time PCR validation of the candidate genes was consistent with the RNA-seq data. Our study provides a strong rationale for the analysis of gene function and breeding of high-quality cotton.

Identification of key proteins related to high-quality fiber in Upland cotton via proteomics analysis

The dynamics of cotton fiber elongation and microfibirl deposition orientation were delineated; advanced ethylene synthesis and redox reaction homeostasis may be crucial for high-quality fiber formation. Fiber length, strength, and fineness determine the use and commercial value of cotton fiber, but their underlying molecular mechanisms remain obscure. We compared the dynamic change trajectories of length, diameter and microfibril orientation angle of the fibers produced by an introgression line SY6167 which generates high-quality fibers even better than Sea island cotton with those of the common-quality fibers from TM-1 across 5 to 30 days post anthesis (DPA). The proteomes were profiled and compared at six representative time points using 2-DE and MS/MS. 14 proteins differentially expressed inside each of cotton line temporally and significantly different tween the two lines were identified. The dynamic change trajectories of fiber length and microfibril angle are close to "s" and reverse "s" growth curves, respectively. SY6167 and TM-1 fibers entered the logarithmic elongation phase simultaneously at 10 DPA, and SY6167 kept elongating logarithmically for 2 more days than TM-1. In parallel to logarithmic elongation, microfibril orientation angles dived sharply, and SY6167 declined faster for a shorter duration than TM-1. 53% of the identified proteins are related to redox homeostasis, and most of them are expressed at higher levels in SY6167 during logarithmic elongation. 1-Aminocyclopropane-1-Carboxylic Acid Oxidase (ACO) started to accumulate at 16 DPA in SY6167, and its encoding genes were highly expressed at this stage, with a much higher level than TM-1. These findings suggest high-quality fibers are associated with high expression of the proteins related to stress and redox homeostasis, the continuously elevated expression of ethylene synthesis ACO gene may play an essential role.

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.

The Pivotal Role of Major Chromosomes of Sub-Genomes A and D in Fiber Quality Traits of Cotton

Lack of precise information about the candidate genes involved in a complex quantitative trait is a major obstacle in the cotton fiber quality improvement, and thus, overall genetic gain in conventional phenotypic selection is low. Recent molecular interventions and advancements in genome sequencing have led to the development of high-throughput molecular markers, quantitative trait locus (QTL) fine mapping, and single nucleotide polymorphisms (SNPs). These advanced tools have resolved the existing bottlenecks in trait-specific breeding. This review demonstrates the significance of chromosomes 3, 7, 9, 11, and 12 of sub-genomes A and D carrying candidate genes for fiber quality. However, chromosome 7 carrying SNPs for stable and potent QTLs related to fiber quality provides great insights for fiber quality-targeted research. This information can be validated by marker-assisted selection (MAS) and transgene in Arabidopsis and subsequently in cotton.

Development and Utilization of Functional Kompetitive Allele-Specific PCR Markers for Key Genes Underpinning Fiber Length and Strength in Gossypium hirsutum L

Fiber length (FL) and fiber strength (FS) are the important indicators of fiber quality in cotton. Longer and stronger fibers are preferred for manufacturing finer yarns in the textile industry. Functional markers (FMs) designed from polymorphic sites within gene sequences attributing to phenotypic variation are highly efficient when used for marker-assisted selection (MAS) in breeding superior varieties with longer FL and higher FS. The aims of this study were to develop FMs via kompetitive allele-specific PCR (KASP) assays and to validate the efficacy of the FMs for allele discrimination and the potential value in practice application. We used four single-nucleotide polymorphism markers and 360 cotton accessions and found that two FMs, namely, D11_24030087 and A07_72204443, could effectively differentiate accessions of different genotypes with higher consistency to phenotype. The appeared frequencies of varieties harbored Hap2 (elite alleles G and T) with longer FL (> the mean of accessions with non-elite allele, 28.50 mm) and higher FS (> the mean of accessions with non-elite allele, 29.06 cN•tex-1) were 100 and 72.7%, respectively, which was higher than that of varieties harbored only on a single elite allele (G or T, 77.9 or 61.9%), suggesting a favorable haplotype for selecting varieties with superior FL and FS. These FMs could be valuable for the high-throughput selection of superior materials by providing genotypic information in cotton breeding programs.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

The hexokinase Gene Family in Cotton: Genome-Wide Characterization and Bioinformatics Analysis

Hexokinase (HXK) is involved in hexose phosphorylation, sugar sensing, and signal transduction, all of which regulate plant growth and adaptation to stresses. Gossypium hirsutum L. is one of the most important fiber crops in the world, however, little is known about the HXKs gene family in G. hirsutum L. We identified 17 GhHXKs from the allotetraploid G. hirsutum L. genome (AADD). G. raimondii (DD) and G. arboreum (AA) are the diploid progenitors of G. hirsutum L. and contributed equally to the At_genome and Dt_genome GhHXKs genes. The chromosomal locations and exon-intron structures of GhHXK genes among cotton species are conservative. Phylogenetic analysis grouped the HXK proteins into four and three groups based on whether they were monocotyledons and dicotyledons, respectively. Duplication event analysis demonstrated that HXKs in G. hirsutum L. primarily originated from segmental duplication, which prior to diploid hybridization. Experiments of qRT-PCR, transcriptome and promoter cis-elements demonstrated that GhHXKs' promoters have auxin and GA responsive elements that are highly expressed in the fiber initiation and elongation stages, while the promoters contain ABA-, MeJA-, and SA-responsive elements that are highly expressed during the synthesis of the secondary cell wall. We performed a comprehensive analysis of the GhHXK gene family is a vital fiber crop, which lays the foundation for future studies assessing its role in fiber development.

Subgenome Bias and Temporal Postponement of Gene Expression Contributes to the Distinctions of Fiber Quality in Gossypium Species

As two cultivated widely allotetraploid cotton species, although Gossypium hirsutum and Gossypium barbadense evolved from the same ancestor, they differ in fiber quality; the molecular mechanism of that difference should be deeply studied. Here, we performed RNA-seq of fiber samples from four G. hirsutum and three G. barbadense cultivars to compare their gene expression patterns on multiple dimensions. We found that 15.90-37.96% of differentially expressed genes showed biased expression toward the A or D subgenome. In particular, interspecific biased expression was exhibited by a total of 330 and 486 gene pairs at 10 days post-anthesis (DPA) and 20 DPA, respectively. Moreover, 6791 genes demonstrated temporal differences in expression, including 346 genes predominantly expressed at 10 DPA in G. hirsutum (TM-1) but postponed to 20 DPA in G. barbadense (Hai7124), and 367 genes predominantly expressed at 20 DPA in TM-1 but postponed to 25 DPA in Hai7124. These postponed genes mainly participated in carbohydrate metabolism, lipid metabolism, plant hormone signal transduction, and starch and sucrose metabolism. In addition, most of the co-expression network and hub genes involved in fiber development showed asymmetric expression between TM-1 and Hai7124, like three hub genes detected at 10 DPA in TM-1 but not until 25 DPA in Hai7124. Our study provides new insights into interspecific expression bias and postponed expression of genes associated with fiber quality, which are mainly tied to asymmetric hub gene network. This work will facilitate further research aimed at understanding the mechanisms underlying cotton fiber improvement.

GSK3s: nodes of multilayer regulation of plant development and stress responses

Glycogen synthase kinase 3 (GSK3) family members are highly conserved serine/threonine protein kinases in eukaryotes. Unlike animals, plants have evolved with multiple homologs of GSK3s involved in a diverse array of biological processes. Emerging evidence suggests that GSK3s act as signaling hubs for integrating perception and transduction of diverse signals required for plant development and responses to abiotic and biotic cues. Here we review recent advances in understanding the molecular interactions between GSK3s and an expanding spectrum of their upstream regulators and downstream substrates in plants. We further discuss how GSK3s act as key signaling nodes of multilayer regulation of plant development and stress response through either being regulated at the post-translational level or regulating their substrates via phosphorylation.

Comparative transcriptome analysis uncovers cell wall reorganization and repressed cell division during cotton fiber initiation

Background

Tetraploid cotton plants serve as prime natural fiber source for the textile industry. Although various omics studies have revealed molecular basis for fiber development, a better understanding of transcriptional regulation mechanism regulating lint fiber initiation is necessary to meet global natural fiber demand.

Results

Here, we aimed to perform transcriptome sequencing to identify DEGs (differentially expressed genes) in ovules of the cotton variety Xu142 and its fibreless mutant Xu142fl during early lint fiber initiation period. Totally, 5516 DEGs including 1840 upregulated and 3676 downregulated were identified. GO enrichment analysis revealed that the downregulated DEGs were mainly associated with biological processes such as transcription related biosynthesis and metabolism, organic cyclic compound biosynthesis and metabolism, photosynthesis, and plant cell wall organization, with molecular functions involving transcription related binding, organic cyclic compound binding, and dioxygenase activity, while the upregulated DEGs were associated with DNA replication and phospholipid biosynthetic related processes. Among the 490 DEGs annotated as transcription factor genes, 86.5% were downregulated in the mutant including the Malvaceae-specific MMLs, expression patterns of which were confirmed during the central period of lint fiber initiation. Investigation of the 16 genes enriched in the cell wall organization revealed that 15 were EXPA coding genes.

Conclusions

Overall, our data indicate that lint fiber initiation is a complicated process involving cooperation of multiple transcription factor families, which might ultimately lead to the reorganization of the cell wall and terminated cell division of the differentiating fiber initials.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12861-021-00247-3.

L2, a chloroplast metalloproteinase, regulates fruit ripening by participating in ethylene autocatalysis under the control of ethylene response factors

Although autocatalytic ethylene biosynthesis plays an important role in the ripening of climacteric fruits, our knowledge of the network that promotes it remains limited. We identified white fruit (wf), a tomato mutant that produces immature fruit that are white and that ripen slowly. We found that an inversion on chromosome 10 disrupts the LUTESCENT2 (L2) gene, and that white fruit is allelic to lutescent2. Using CRISPR/Cas9 technology we knocked out L2 in wild type tomato and found that the l2-cr mutants produced phenotypes that were very similar to white fruit (lutescent2). In the l2-cr fruit, chloroplast development was impaired and the accumulation of carotenoids and lycopene occurred more slowly than in wild type. During fruit ripening in l2-cr mutants, the peak of ethylene release was delayed, less ethylene was produced, and the expression of ACO genes was significantly suppressed. We also found that exogenous ethylene induces the expression of L2 and that ERF.B3, an ethylene response factor, binds to the promoter of the L2 gene and activates its transcription. Thus, the expression of L2 is regulated by exogenous ethylene. Taken together, our results indicate that ethylene may affect the expression of L2 gene and that L2 participates in autocatalytic ethylene biosynthesis during tomato fruit ripening.

The regulation of ACC synthase protein turnover: a rapid route for modulating plant development and stress responses

The phytohormone ethylene regulates plant growth, development, and stress responses. The strict fine-tuning of the regulation of ethylene biosynthesis contributes to the diverse roles of ethylene in plants. Pyridoxal 5'-phosphate-dependent 1-aminocyclopropane-1-carboxylic acid synthase, a rate-limiting enzyme in ethylene biosynthesis, is central and often rate-limiting to regulate ethylene concentration in plants. The post-translational regulation of ACS is a major pathway controlling ethylene biosynthesis in response to various stimuli. We conclude that the regulation of ACS turnover may serve as a central hub for the rapid integration of developmental, environmental, and hormonal signals, all of which influence plant growth and stress responses.

Plant Cell Identity in the Era of Single-Cell Transcriptomics

High-throughput single-cell transcriptomic approaches have revolutionized our view of gene expression at the level of individual cells, providing new insights into their heterogeneity, identities, and functions. Recently, technical challenges to the application of single-cell transcriptomics to plants have been overcome, and many plant organs and tissues have now been subjected to analyses at single-cell resolution. In this review, we describe these studies and their impact on our understanding of the diversity, differentiation, and activities of plant cells. We particularly highlight their impact on plant cell identity, including unprecedented views of cell transitions and definitions of rare and novel cell types. We also point out current challenges and future opportunities for the application and analyses of single-cell transcriptomics in plants. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Phosphorylation of WRKY16 by MPK3-1 is essential for its transcriptional activity during fiber initiation and elongation in cotton (Gossypium hirsutum)

Cotton, one of the most important crops in the world, produces natural fiber materials for the textile industry. WRKY transcription factors play important roles in plant development and stress responses. However, little is known about whether and how WRKY transcription factors regulate fiber development of cotton so far. In this study, we show that a fiber-preferential WRKY transcription factor, GhWRKY16, positively regulates fiber initiation and elongation. GhWRKY16-silenced transgenic cotton displayed a remarkably reduced number of fiber protrusions on the ovule and shorter fibers compared to the wild-type. During early fiber development, GhWRKY16 directly binds to the promoters of GhHOX3, GhMYB109, GhCesA6D-D11, and GhMYB25 to induce their expression, thereby promoting fiber initiation and elongation. Moreover, GhWRKY16 is phosphorylated by the mitogen-activated protein kinase GhMPK3-1 at residues T-130 and S-260. Phosphorylated GhWRKY16 directly activates the transcription of GhMYB25, GhHOX3, GhMYB109, and GhCesA6D-D11 for early fiber development. Thus, our data demonstrate that GhWRKY16 plays a crucial role in fiber initiation and elongation, and that GhWRKY16 phosphorylation by GhMPK3-1 is essential for the transcriptional activation on downstream genes during the fiber development of cotton.

VvMYB114 mediated by miR828 negatively regulates trichome development of Arabidopsis

Trichome is a specialized structure differentiated during the morphogenesis of plant leaf epidermal cells. In recent years, with the continuous researches on trichome development of Arabidopsis and other plants, more and more genes related to trichome morphogenesis have been discovered, including R2R3-type MYB genes. In this study, we cloned a R2R3-type MYB family gene from grape, VvMYB114, a target gene of vvi-miR828. qRT-PCR showed that VvMYB114 mRNA accumulated during grape fruit ripening, and VvMYB114 protein had transcriptional activation activity. Heterologous overexpression of VvMYB114 in Arabidopsis reduced the number of trichome on leaves and stems. Mutating the miR828-binding site in VvMYB114 without altering amino-acid sequence had no effect on trichome development in Arabidopsis. The results showed a different role of the regulation of miR828 to VvMYB114 in Arabidopsis from in grape, which indicated the functional divergence of miRNA targeting homoeologous genes in different species played an important roles in evolution and useful trait selection.

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.

Large-fragment insertion activates gene GaFZ (Ga08G0121) and is associated with the fuzz and trichome reduction in cotton (Gossypium arboreum)

Cotton seeds are typically covered by lint and fuzz fibres. Natural 'fuzzless' mutants are an ideal model system for identifying genes that regulate cell initiation and elongation. Here, using a genome-wide association study (GWAS), we identified a ~ 6.2 kb insertion, larINDEL_FZ , located at the end of chromosome 8, composed of a ~ 5.0 kb repetitive sequence and a ~ 1.2 kb fragment translocated from chromosome 12 in fuzzless Gossypium arboreum. The presence of larINDEL_FZ was associated with a fuzzless seed and reduced trichome phenotypes in G. arboreum. This distant insertion was predicted to be an enhancer, located ~ 18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). Ectopic overexpression of GaFZ in Arabidopsis thaliana and G. hirsutum suggested that GaFZ negatively modulates fuzz and trichome development. Co-expression and interaction analyses demonstrated that GaFZ might impact fuzz fibre/trichome development by repressing the expression of genes in the very-long-chain fatty acid elongation pathway. Thus, we identified a novel regulator of fibre/trichome development while providing insights into the importance of noncoding sequences in cotton.

LCM and RNA-seq analyses revealed roles of cell cycle and translational regulation and homoeolog expression bias in cotton fiber cell initiation

Background

Cotton fibers provide a powerful model for studying cell differentiation and elongation. Each cotton fiber is a singular and elongated cell derived from epidermal-layer cells of a cotton seed. Efforts to understand this dramatic developmental shift have been impeded by the difficulty of separation between fiber and epidermal cells.

Results

Here we employed laser-capture microdissection (LCM) to separate these cell types. RNA-seq analysis revealed transitional differences between fiber and epidermal-layer cells at 0 or 2 days post anthesis. Specifically, down-regulation of putative cell cycle genes was coupled with upregulation of ribosome biosynthesis and translation-related genes, which may suggest their respective roles in fiber cell initiation. Indeed, the amount of fibers in cultured ovules was increased by cell cycle progression inhibitor, Roscovitine, and decreased by ribosome biosynthesis inhibitor, Rbin-1. Moreover, subfunctionalization of homoeologs was pervasive in fiber and epidermal cells, with expression bias towards 10% more D than A homoeologs of cell cycle related genes and 40–50% more D than A homoeologs of ribosomal protein subunit genes. Key cell cycle regulators were predicted to be epialleles in allotetraploid cotton. MYB-transcription factor genes displayed expression divergence between fibers and ovules. Notably, many phytohormone-related genes were upregulated in ovules and down-regulated in fibers, suggesting spatial-temporal effects on fiber cell development.

Conclusions

Fiber cell initiation is accompanied by cell cycle arrest coupled with active ribosome biosynthesis, spatial-temporal regulation of phytohormones and MYB transcription factors, and homoeolog expression bias of cell cycle and ribosome biosynthesis genes. These valuable genomic resources and molecular insights will help develop breeding and biotechnological tools to improve cotton fiber production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07579-1.

Cotton pan-genome retrieves the lost sequences and genes during domestication and selection

BACKGROUND

Millennia of directional human selection has reshaped the genomic architecture of cultivated cotton relative to wild counterparts, but we have limited understanding of the selective retention and fractionation of genomic components.

RESULTS

We construct a comprehensive genomic variome based on 1961 cottons and identify 456 Mb and 357 Mb of sequence with domestication and improvement selection signals and 162 loci, 84 of which are novel, including 47 loci associated with 16 agronomic traits. Using pan-genome analyses, we identify 32,569 and 8851 non-reference genes lost from Gossypium hirsutum and Gossypium barbadense reference genomes respectively, of which 38.2% (39,278) and 14.2% (11,359) of genes exhibit presence/absence variation (PAV). We document the landscape of PAV selection accompanied by asymmetric gene gain and loss and identify 124 PAVs linked to favorable fiber quality and yield loci.

CONCLUSIONS

This variation repertoire points to genomic divergence during cotton domestication and improvement, which informs the characterization of favorable gene alleles for improved breeding practice using a pan-genome-based approach.

Identify of Fast-Growing Related Genes Especially in Height Growth by Combining QTL Analysis and Transcriptome in Salix matsudana (Koidz)

The study on the fast-growing traits of trees, mainly valued by tree height (TH) and diameter at breast height (DBH), is of great significance to promote the development of the forest industry. Quantitative trait locus (QTL) mapping based on high-density genetic maps is an efficient approach to identify genetic regions for fast-growing traits. In our study, a high-density genetic map for the F₁ population was constructed. The genetic map had a total size of 5,484.07 centimorgan (cM), containing 5,956 single nucleotide polymorphisms (SNPs) based on Specific Length Amplified Fragment sequencing. Six fast-growing related stable QTL were identified on six chromosomes, and five stable QTL were identified by a principal component analysis (PCA). By combining the RNA-seq analysis for the two parents and two progenies with the qRT-PCR analysis, four candidate genes, annotated as DnaJ, 1-aminocyclopropane-1-carboxylate oxidase 1 (ACO1), Caffeic acid 3-O-methyltransferase 1 (COMT1), and Dirigent protein 6 (DIR6), that may regulate height growth were identified. Several lignin biosynthesis-related genes that may take part in height growth were detected. In addition, 21 hotspots in this population were found. The results of this study will provide an important foundation for further studies on the molecular and genetic regulation of TH and DBH.

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.

GhPIPLC2D promotes cotton fiber elongation by enhancing ethylene biosynthesis

Inositol-1,4,5-trisphosphate (IP₃) is an important second messenger and one of the products of phosphoinositide-specific phospholipase C (PIPLC)-mediated phosphatidylinositol (4,5) bisphosphate (PIP₂) hydrolysis. However, the function of IP₃ in cotton is unknown. Here, we characterized the function of GhPIPLC2D in cotton fiber elongation. GhPIPLC2D was preferentially expressed in elongating fibers. Suppression of GhPIPLC2D transcripts resulted in shorter fibers and decreased IP₃ accumulation and ethylene biosynthesis. Exogenous application of linolenic acid (C18:3) and phosphatidylinositol (PI), the precursor of IP₃, improved IP₃ and myo-inositol-1,2,3,4,5,6-hexakisphosphate (IP₆) accumulation, as well as ethylene biosynthesis. Moreover, fiber length in GhPIPLC2D-silenced plant was reduced after exogenous application of IP₆ and ethylene. These results indicate that GhPIPLC2D positively regulates fiber elongation and IP₃ promotes fiber elongation by enhancing ethylene biosynthesis. Our study broadens our understanding of the function of IP₃ in cotton fiber elongation and highlights the possibility of cultivating better cotton varieties by manipulating GhPIPLC2D in the future.

•

GhPIPLC2D positively regulates cotton fiber elongation

•

GhPIPLC2D cleaves PIP₂ into IP₃, which could be phosphorylated to IP₆

•

IP₆ enhances fiber elongation via improving ethylene biosynthesis

Selection and Validation of Reference Genes for qRT-PCR Analysis in the Oil-Rich Tuber Crop Tiger Nut (Cyperus esculentus) Based on Transcriptome Data

Tiger nut (Cyperus esculentus), a perennial C₄ plant of the Cyperaceae family, is an unconventional crop that is distinguished by its oil-rich tubers, which also possesses the advantages of strong resistance, wide adaptability, short life periods, and large biomass. To facilitate studies on gene expression in this species, we identified and validated a series of reference genes (RGs) based on transcriptome data, which can be employed as internal controls for qRT-PCR analysis in tiger nut. Fourteen putative candidate RGs were identified and evaluated across nine different tissues of two cultivars, and the RGs were analyzed using three different algorithms (geNorm, NormFinder, and BestKeeper). The stability rankings of the candidate RGs were merged into consensus lists with RankAggreg. For the below-ground storage organ of tiger nut, the optimal RGs were TUB4 and UCE2 in different developmental stages of tubers. UCE2 and UBL5 were the most stably expressed RGs among all tissues, while Rubisco and PGK exhibited the lowest expression stability. UCE2, UBL5 and Rubisco were compared to normalize the expression levels of the caleosin (CLO) and diacylglycerol acyltransferase 2-2 (DGAT2-2) genes across the same tissues. Our results showed that the RGs identified in this study, which exhibit more uniform expression patterns, may be utilized for the normalization of qRT-PCR results, promoting further research on gene expression in various tissues of tiger nut.

Analysis and review of trichomes in plants

BACKGROUND

Trichomes play a key role in the development of plants and exist in a wide variety of species.

RESULTS

In this paper, it was reviewed that the structure and morphology characteristics of trichomes, alongside the biological functions and classical regulatory mechanisms of trichome development in plants. The environment factors, hormones, transcription factor, non-coding RNA, etc., play important roles in regulating the initialization, branching, growth, and development of trichomes. In addition, it was further investigated the atypical regulation mechanism in a non-model plant, found that regulating the growth and development of tea (Camellia sinensis) trichome is mainly affected by hormones and the novel regulation factors.

CONCLUSIONS

This review further displayed the complex and differential regulatory networks in trichome initiation and development, provided a reference for basic and applied research on trichomes in plants.

QTL controlling fiber quality traits under salt stress in upland cotton (Gossypium hirsutum L.)

QTL for fiber quality traits under salt stress discerned candidate genes controlling fatty acid metabolism. Salinity stress seriously affects plant growth and limits agricultural productivity of crop plants. To dissect the genetic basis of response to salinity stress, a recombinant inbred line population was developed to compare fiber quality in upland cotton (Gossypium hirsutum L.) under salt stress and normal conditions. Based on three datasets of (1) salt stress, (2) normal growth, and (3) the difference value between salt stress and normal conditions, 51, 70, and 53 QTL were mapped, respectively. Three QTL for fiber length (FL) (qFL-Chr1-1, qFL-Chr5-5, and qFL-Chr24-4) were detected under both salt and normal conditions and explained 4.26%, 9.38%, and 3.87% of average phenotypic variation, respectively. Seven genes within intervals of two stable QTL (qFL-Chr1-1 and qFL-Chr5-5) were highly expressed in lines with extreme long fiber. A total of 35 QTL clusters comprised of 107 QTL were located on 18 chromosomes and exhibited pleiotropic effects. Thereinto, two clusters were responsible for improving five fiber quality traits, and 6 influenced FL and fiber strength (FS). The QTL with positive effect for fiber length exhibited active effects on fatty acid synthesis and elongation, but the ones with negative effect played passive roles on fatty acid degradation under salt stress.

Genome-wide characterization of the WAK gene family and expression analysis under plant hormone treatment in cotton

BACKGROUND

Wall-associated kinases (WAK), one of the receptor-like kinases (RLK), function directly in the connection and communication between the plant cell wall and the cytoplasm. WAK genes are highly conserved and have been identified in plants, such as rice, but there is little research on the WAK gene family in cotton.

RESULTS

In the present study, we identified 29 GhWAK genes in Gossypium hirsutum. Phylogenetic analysis showed that cotton WAK proteins can be divided into five clades. The results of synteny and Ka/Ks analysis showed that the GhWAK genes mainly originated from whole genome duplication (WGD) and were then mainly under purifying selection. Transcriptome data and real-time PCR showed that 97% of GhWAK genes highly expressed in cotton fibers and ovules. β-glucuronidase (GUS) staining assays showed that GhWAK5 and GhWAK16 expressed in Arabidopsis leaf trichomes. Fourteen GhWAK genes were found to possess putative gibberellin (GA) response elements in the promoter regions, 13 of which were significantly induced by GA treatment. Ten GhWAK genes contained auxin (IAA) response elements and the expression level of nine GhWAKs significantly increased under auxin treatment.

CONCLUSIONS

We provide a preliminary analysis of the WAK gene family in G. hirsutum, which sheds light on the potantial roles of GhWAK genes in cotton fiber cell development. Our data also provides a useful resource for future studies on the functional roles of GhWAK genes.

Genome-Wide Analysis of MDHAR Gene Family in Four Cotton Species Provides Insights into Fiber Development via Regulating AsA Redox Homeostasis

Monodehydroasorbate reductase (MDHAR) (EC1.6.5.4), a key enzyme in ascorbate-glutathione recycling, plays important roles in cell growth, plant development and physiological response to environmental stress via control of ascorbic acid (AsA)-mediated reduction/oxidation (redox) regulation. Until now, information regarding MDHAR function and regulatory mechanism in Gossypium have been limited. Herein, a genome-wide identification and comprehensive bioinformatic analysis of 36 MDHAR family genes in four Gossypium species, Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, were performed, indicating their close evolutionary relationship. Expression analysis of GhMDHARs in different cotton tissues and under abiotic stress and phytohormone treatment revealed diverse expression features. Fiber-specific expression analysis showed that GhMDHAR1A/D, 3A/D and 4A/D were preferentially expressed in fiber fast elongating stages to reach peak values in 15-DPA fibers, with corresponding coincident observances of MDHAR enzyme activity, AsA content and ascorbic acid/dehydroascorbic acid (AsA/DHA) ratio. Meanwhile, there was a close positive correlation between the increase of AsA content and AsA/DHA ratio catalyzed by MDHAR and fiber elongation development in different fiber-length cotton cultivars, suggesting the potential important function of MDHAR for fiber growth. Following H₂O₂ stimulation, GhMDHAR demonstrated immediate responses at the levels of mRNA, enzyme, the product of AsA and corresponding AsA/DHA value, and antioxidative activity. These results for the first time provide a comprehensive systemic analysis of the MDHAR gene family in plants and the four cotton species and demonstrate the contribution of MDHAR to fiber elongation development by controlling AsA-recycling-mediated cellular redox homeostasis.

UDP-Glucose Dehydrogenases: Identification, Expression, and Function Analyses in Upland Cotton (Gossypium hirsutum)

UDP-glucose dehydrogenase (UGD; EC1.1.1.22) is a NAD⁺-dependent enzyme that catalyzes the two-fold oxidation of UDP-glucose (UDP-Glc) to produce UDP-glucuronic acid and plays an important role in plant cell wall synthesis. A total of 42 UGD genes from four Gossypium genomes including G. hirsutum, G. arboretum, G. barbadense, and G. raimondii were identified and found that the UGD gene family has conservative evolution patterns in gene structure and protein domain. The growth of fibers can be effectively promoted after adding the UDP-Glc to the medium, and the GhUGD gene expression enhanced. In addition, the transgenic Arabidopsis lines over-expressing GH_D12G1806 had longer root lengths and higher gene expression level than the wild-type plants of Columbia-0. These results indicated that UGD may play important roles in cotton fiber development and has a guiding significance for dissecting fiber development mechanism.

GrTCP11, a Cotton TCP Transcription Factor, Inhibits Root Hair Elongation by Down-Regulating Jasmonic Acid Pathway in Arabidopsis thaliana

TCP transcription factors play important roles in diverse aspects of plant development as transcriptional activators or repressors. However, the functional mechanisms of TCPs are not well understood, especially in cotton fibers. Here, we identified a total of 37 non-redundant TCP proteins from the diploid cotton (Gossypium raimondii), which showed great diversity in the expression profile. GrTCP11, an ortholog of AtTCP11, was preferentially expressed in cotton anthers and during fiber initiation and secondary cell wall synthesis stages. Overexpression of GrTCP11 in Arabidopsis thaliana reduced root hair length and delayed flowering. It was found that GrTCP11 negatively regulated genes involved in jasmonic acid (JA) biosynthesis and response, such as AtLOX4, AtAOS, AtAOC1, AtAOC3, AtJAZ1, AtJAZ2, AtMYC2, and AtERF1, which resulted in a decrease in JA concentration in the overexpressed transgenic lines. As with the JA-deficient mutant dde2-2, the transgenic line 4-1 was insensitive to 50 μM methyl jasmonate, compared with the wild-type plants. The results suggest that GrTCP11 may be an important transcription factor for cotton fiber development, by negatively regulating JA biosynthesis and response.

Auxin Directly Upregulates GhRAC13 Expression to Promote the Onset of Secondary Cell Wall Deposition in Cotton Fibers

Cotton fibers are single cells that show a relatively independent developmental process of cell differentiation, elongation, and secondary wall deposition. Auxin promotes fiber cell protrusion from the surface of the ovule. However, the role of auxin at other stages of cotton fiber development remains largely unknown. To gain a deeper insight into this aspect, we measured indoleacetic acid (IAA) content in developing fibers. Results showed an increase in IAA content at the transition stage from elongation to secondary cell wall deposition. Subsequently, we investigated the differences between two transgenic cottons that show upregulated and downregulated fiber auxin levels, respectively. In planta analysis revealed that, in addition to promoting cell elongation, auxin regulated the time of initiation of reactive oxygen species (ROS) production and secondary wall deposition in cotton fibers. This was closely correlated with the upregulated expression of GhRAC13, which regulates ROS-triggered cellulose synthesis. We found multiple putative auxin-responsive elements existed within the promoter region of GhRAC13, and IAA could induce proGhRAC13 activity. The dual-luciferase reporter assay further proved the activation of proGhRAC13 by GhARF5, an auxin-signaling activator. Altogether, our results suggest a role of auxin in promoting the onset of secondary growth by directly upregulating GhRAC13 expression in cotton fibers.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

Improved lint yield under field conditions in cotton over-expressing transcription factors regulating fibre initiation

Only a few transcription factors (TFs) regulating which cells of the ovule epidermis differentiate into lint fibres have been identified in cotton (Gossypium hirsutum L.). In this study, the effect on lint yield and fibre quality of over-expressing three TFs in cotton, GhHD-1, GhMYB25 and GhMYB25Like, and their double and triple combinations, were evaluated in field experiments over two seasons. The expression of single or stacked TFs were all driven either by an ovule-specific promoter, FBP 7, or a constitutive promoter, Stunt 7, in a Coker 315 background. TF type, either singly or in combination, was found to be the most significant factor affecting lint yield. Among 64 transgenic lines tested, seven were higher yielding than null segregant lines in one or both seasons and were all from the sets with single and double over-expressed TF combinations. A reduced yield was associated with the set of triple combinations. The two most stable high yielding lines across the seasons recorded 12-22% higher yields than the nulls, although were not competitive to locally adapted commercial controls. Over-expression of TFs singly or in combination did not significantly alter fibre length and strength, but sometimes increased fibre micronaire. There were positive relationships between lint yield and lint percentage and lint yield and fibre density amongst the transgenic lines. Our preliminary results suggest that manipulating TF expression, either singly or in pairs, can increase the density of fibres initiated on developing seeds and fibre yields under field conditions while maintaining overall fibre quality.

Sequencing Multiple Cotton Genomes Reveals Complex Structures and Lays Foundation for Breeding

Cotton is a major fiber plant, which provides raw materials for clothing, protecting humans from the harsh environment of cold or hot weathers, enriching the culture and custom of human societies. Due to its importance, the diploid and tetraploid genomes of different cotton plants have been repeatedly sequenced to obtain their complete and fine genome sequences. These valuable genome data sets revealed the evolutionary past of the cotton plants, which were recursively affected by polyploidization, with a decaploidization contributing to the formation of the genus Gossypium, and a neo-tetraploidization contributing to the formation of nowadays widely cultivated cotton plants. Post-polyploidization genome instability resulted in numerous structural changes of the genomes, such as gene loss, DNA inversion and translocation, illegitimate recombination, and accumulation of repetitive sequences, and functional innovation accompanied by elevated evolutionary rates of genes. Many these changes have been asymmetric between subgnomes of the tetraploid cottons, rendering their divergent profiles of biological regulation and function. The availability of whole-genome sequences has now paved the way to identify and clone functional genes, e.g., those relating to fiber development, and to enhance breeding efforts to cultivate cottons to produce high-yield and high-quality fibers, and to resist environmental and biological stress.

Transcriptome Analysis Reveals a Gene Expression Pattern Associated with Fuzz Fiber Initiation Induced by High Temperature in Gossypium barbadense

Gossypium barbadense is an important source of natural textile fibers, as is Gossypium hirsutum. Cotton fiber development is often affected by various environmental factors, such as abnormal temperature. However, little is known about the underlying mechanisms of temperature regulating the fuzz fiber initiation. In this study, we reveal that high temperatures (HT) accelerate fiber development, improve fiber quality, and induced fuzz initiation of a thermo-sensitive G. barbadense variety L7009. It was proved that fuzz initiation was inhibited by low temperature (LT), and 4 dpa was the stage most susceptible to temperature stress during the fuzz initiation period. A total of 43,826 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Of these, 9667 were involved in fiber development and temperature response with 901 transcription factor genes and 189 genes related to plant hormone signal transduction. Further analysis of gene expression patterns revealed that 240 genes were potentially involved in fuzz initiation induced by high temperature. Functional annotation revealed that the candidate genes related to fuzz initiation were significantly involved in the asparagine biosynthetic process, cell wall biosynthesis, and stress response. The expression trends of sixteen genes randomly selected from the RNA-seq data were almost consistent with the results of qRT-PCR. Our study revealed several potential candidate genes and pathways related to fuzz initiation induced by high temperature. This provides a new view of temperature-induced tissue and organ development in Gossypium barbadense.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

Understanding the role of phytohormones in cotton fiber development through omic approaches; recent advances and future directions

Cotton is among the most important fiber crops for the textile-based industry, thanks to its cellulose-rich mature fibers. The fiber initiation and elongation are one of the best models for deciphering mechanisms of single-cell differentiation and growth, that also target of fiber development programs. During the last couple of decades, high yielding omics approaches (genomics, transcriptomics, and proteomics), have helped in the identification of several genes and gene products involved in fiber development along with functional relationship to phytohormones. For example, MYB transcription factor family and Sus gene family have been evidenced by controlling cotton fiber initiation. Most importantly, the biosynthesis, responses, and transporting of phytohormones is documented to participate in the initiation of cotton fibers. Herein, in this review, the reliable genetic evidence by manipulating the above genes in cotton have been summarized to describe the relationships among key phytohormones, transcription factors, proteins, and downstream fiber growth-related genes such as Sus. The effect of other important factors such as ROS, fatty acid metabolism, and actin (globular multi-functional proteins) over fiber development has also been discussed. The challenges and deficiencies in the research of cotton fiber development have been mentioned along with a future perspective to discover new crucial genes using multiple omics analysis.

The Effect of Ethylene on the Color Change and Resistance to Botrytis cinerea Infection in 'Kyoho' Grape Fruits

The formation of grape quality and the mechanism of resistance against foreign pathogens affect the storage stability of fruits during post-harvest handling. Ethylene plays a crucial role in regulating the ripeness of fruits and can be used as an exogenous regulator to resist exogenous pathogens. In this study, we used different concentrations of ethephon for treatment of grape fruits before veraison, analyzed the anthocyanin content, soluble solids, titratable acid, and determined fruit firmness and cell wall metabolism-related enzymes during fruit development. Results showed that exogenous ethephon promoted the early coloration of grape fruits and increased the coloring-related genes myeloblastosis A1(MYBA1), myeloblastosis A2(MYBA2), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3’-hydroxylase gene (F3’H), flavonoid 3’, 5’hydroxylase (F3’5’H), 3-O-flavonoid glucosyltransferase (UFGT), and glutathione S-transferase (GST), softening related genes Polygalacturonase(PG), pectinate lyases(PL) and Pectin methylesterase( PME, as well as ethylene metabolism pathway-related genes 1-aminocyclopropane-1-carboxylic acid synthase 1(ACS1), 1-aminocyclopropane-1-carboxylic acid oxidase 2 (ACO2), ethylene receptor gene(ETR2), and ethylene-insensitive 3 (EIN3). Ethephon treatment also increased soluble solids and decreased titratable acid in grape fruit. Fruits pretreated with ethephon were inoculated with Botrytis cinerea, which led to resistance in grape fruit through activation of the antioxidant system. The expression levels of disease resistance-related genes including VvPAD4, VvPIP1, VvNAC26, VvDREB, VvAPX, Vvpgip, VvWRKY70, VvMYC2, VvNPR1 also increased in inoculated fruit with pathogen following ethephon pretreatment. Furthermore, we monitored ethylene response factor 1(ERF1) transcription factor, which could interact with protein EIN3 during ethylene signal transduction and mediate fruit resistance against B. cinerea infection. Meanwhile, overexpression of VvERF1 vectorin strawberry fruits reduced the susceptibility to B. cinerea infection. We suggest that ethylene can induce resistance in ripened fruits after B. cinerea infection and provide adequate postharvest care.

Comparative phosphoproteomic analysis of BR-defective mutant reveals a key role of GhSK13 in regulating cotton fiber development

Brassinosteroid (BR), a steroid phytohormone, whose signaling transduction pathways include a series of phosphorylation and dephosphorylation events, and GSK3s are the main negative regulator kinases. BRs have been shown to play vital roles in cotton fiber elongation. However, the underlying mechanism is still elusive. In this study, fibers of a BR-defective mutant Pagoda 1 (pag1), and its corresponding wild-type (ZM24) were selected for a comparative global phosphoproteome analysis at critical developmental time points: fast-growing stage (10 days after pollination (DPA)) and secondary cell wall synthesis stage (20 DPA). Based on the substrate characteristics of GSK3, 900 potential substrates were identified. Their GO and KEGG annotation results suggest that BR functions in fiber development by regulating GhSKs (GSK3s of Gossypium hirsutum L.) involved microtubule cytoskeleton organization, and pathways of glucose, sucrose and lipid metabolism. Further experimental results revealed that among the GhSK members identified, GhSK13 not only plays a role in BR signaling pathway, but also functions in developing fiber by respectively interacting with an AP2-like ethylene-responsive factor GhAP2L, a nuclear transcription factor Gh_DNF_YB19, and a homeodomain zipper member GhHDZ5. Overall, our phosphoproteomic research advances the understanding of fiber development controlled by BR signal pathways especially through GhSKs, and also offers numbers of target proteins for improving cotton fiber quality.

Genome-wide identification of the mitogen-activated protein kinase (MAPK) family in cotton (Gossypium hirsutum) reveals GhMPK6 involved in fiber elongation

Mitogen-activated protein kinases (MAPKs) are important in regulating plant development as well as stress response. In this study, we genome-widely identified 56 MAPK genes in upland cotton. These MAPK genes unequally distribute on 22 chromosomes of cotton genome, but no MAPK gene is located on At_Chr6, Dt_Chr6, At_Chr13 and Dt_Chr13. The exons and introns in GhMAPK gene family vary widely at the position, number and length. Furthermore, GhMAPK family can be divided into 4 groups (A, B, C and D), and the TEY type of T-loop exists in three groups (A, B and C), but the TDY type of T-loop is only in group D. Further study revealed that some GhMAPK genes (including GhMPK6) are preferentially expressed in elongating fibers. GhMPK6 maintains a high phosphorylation level in elongating fibers, and its phosphorylation was enhanced in fibers by phytohormones brassinosteroid (BR), ethylene and indole-3-acetic acid (IAA). Additionally, GhMPK6 could interact with GhMKK2-2 and GhMKK4, suggesting that GhMKK2-2/4-GhMPK6 module may be involved in phosphorylation of its downstream proteins for regulating fiber elongation of cotton.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li ₂ ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li ₂ , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li ₂ genotype) SSR7 and SSR18 showed complete linkage with Li ₂ Nonrandom distribution of markers with null alleles suggests that the Li ₂ phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li ₂ fiber, and virus-induced silencing of these genes in the wild type induced Li ₂ -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li ₂ mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

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.

Nitrogen and mepiquat chloride can affect fiber quality and cotton yield

ABSTRACT The use of mepiquat chloride (MC) in cotton (Gossypium hirsutum L.) cultivation has increased significantly in recent years. The use of MC, a growth regulator, results in higher precocity and lower height and consequently increases yield and reduces crop costs. The objective of this study was to evaluate the effects of nitrogen (N) and MC doses on fiber quality and cotton yield. This study was conducted during the 2012/13 and 2013/14 harvests in the experimental area of the Fundação de Apoio à Pesquisa Agropecuária de Chapadão, located in the municipality of Chapadão do Sul, MS, Brazil, at 18° 48’ 45.9” S, 52° 36’ 3” W, having an altitude of 820 m. The experimental design was a randomized blocks in a 5 × 4 factorial scheme comprised of five N doses (0, 40, 80, 120, 160 kg ha-1) and four MC doses (0, 50, 100, 150 mL ha-1) with four repetitions. Analyses of fiber quality included: length, uniformity, elongation, strength, micronaire index, reflectance, yellowing degrees and short fiber content. Regarding cotton yield, green color index (GCI), plant height at harvest time, and cotton seed yield were determined. The best results for fiber quality and productivity occurr with the application of 76.8 mL ha-1 of MC to the cotton crops when the N dose is greater than 80 kg ha-1. Increasing N doses causes an increase in cotton yield and micronaire index, but there is a decrease in fiber uniformity.

A basic helix-loop-helix protein (GhFP1) promotes fibre elongation of cotton (Gossypium hirsutum) by modulating brassinosteroid biosynthesis and signalling

Basic helix-loop-helix (bHLH) proteins are involved in transcriptional networks controlling a number of biological processes in plants. However, little information is known on the roles of bHLH proteins in cotton fibre development so far. Here, we show that a cotton bHLH protein (GhFP1) positively regulates fibre elongation. GhFP1 transgenic cotton and Arabidopsis plants were generated to study how GhFP1 regulates fibre cell elongation. Fibre length of the transgenic cotton overexpressing GhFP1 was significantly longer than that of wild-type, whereas suppression of GhFP1 expression hindered fibre elongation. Furthermore, overexpression of GhFP1 in Arabidopsis promoted trichome development. Expression of the brassinosteroid (BR)-related genes was markedly upregulated in fibres of GhFP1 overexpression cotton, but downregulated in GhFP1-silenced fibres. BR content in the transgenic fibres was significantly altered, relative to that in wild-type. Moreover, GhFP1 protein could directly bind to the promoters of GhDWF4 and GhCPD to activate expression of these BR-related genes. Therefore, our data suggest that GhFP1 as a positive regulator participates in controlling fibre elongation by activating BR biosynthesis and signalling. Additionally, homodimerisation of GhFP1 may be essential for its function, and interaction between GhFP1 and other cotton bHLH proteins may interfere with its DNA-binding activity.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.