Epidermal cell differentiation in cotton mediated by the homeodomain leucine zipper gene, GhHD-1

GhHDZ76, a cotton HD-Zip transcription factor, involved in regulating the initiation and early elongation of cotton fiber development in G. hirsutum

In this study, the whole HD-Zip family members of G. hirsutum were identified, and GhHDZ76 was classified into the HD-Zip IV subgroup. GhHDZ76 was predominantly expressed in the 0-5 DPA of fiber development stage and localized in the nucleus. Overexpression of GhHDZ76 significantly increased the length and density of trichomes in Arabidopsis thaliana. The fiber length of GhHDZ76 knockout lines by CRISPR/Cas9 was significantly shorter than WT at the early elongation and mature stage, indicating that GhHDZ76 positively regulate the fiber elongation. Scanning electron microscopy showed that the number of ovule surface protrusion of 0 DPA of GhHDZ76 knockout lines was significantly lower than WT, suggesting that GhHDZ76 can also promote the initiation of fiber development. The transcript level of GhWRKY16, GhRDL1, GhEXPA1 and GhMYB25 genes related to fiber initiation and elongation in GhHDZ76 knockout lines were significantly decreased. Yeast two-hybrid and Luciferase complementation imaging (LCI) assays showed that GhHDZ76 can interact with GhWRKY16 directly. As a transcription factor, GhHDZ76 has transcriptional activation activity, which could bind to L1-box elements of the promoters of GhRDL1 and GhEXPA1. Double luciferase reporter assay showed that the GhWRKY16 could enhance the transcriptional activity of GhHDZ76 to pGhRDL1, but it did not promote the transcriptional activity of GhHDZ76 to pGhEXPA1. GhHDZ76 protein may also promote the transcriptional activity of GhWRKY16 to the downstream target gene GhMYB25. Our results provided a new gene resource for fiber development and a theoretical basis for the genetic improvement of cotton fiber quality.

Phosphatidic acid interacts with an HD-ZIP transcription factor GhHOX4 to influence its function in fiber elongation of cotton (Gossypium hirsutum)

Upland cotton, the mainly cultivated cotton species in the world, provides over 90% of natural raw materials (fibers) for the textile industry. The development of cotton fibers that are unicellular and highly elongated trichomes on seeds is a delicate and complex process. However, the regulatory mechanism of fiber development is still largely unclear in detail. In this study, we report that a homeodomain-leucine zipper (HD-ZIP) IV transcription factor, GhHOX4, plays an important role in fiber elongation. Overexpression of GhHOX4 in cotton resulted in longer fibers, while GhHOX4-silenced transgenic cotton displayed a "shorter fiber" phenotype compared with wild type. GhHOX4 directly activates two target genes, GhEXLB1D and GhXTH2D, for promoting fiber elongation. On the other hand, phosphatidic acid (PA), which is associated with cell signaling and metabolism, interacts with GhHOX4 to hinder fiber elongation. The basic amino acids KR-R-R in START domain of GhHOX4 protein are essential for its binding to PA that could alter the nuclear localization of GhHOX4 protein, thereby suppressing the transcriptional regulation of GhHOX4 to downstream genes in the transition from fiber elongation to secondary cell wall (SCW) thickening during fiber development. Thus, our data revealed that GhHOX4 positively regulates fiber elongation, while PA may function in the phase transition from fiber elongation to SCW formation by negatively modulating GhHOX4 in cotton.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum

KEY MESSAGE

GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

Dissecting the major genetic components underlying cotton lint development

Numerous genetic loci and several functionally characterized genes have been linked to determination of lint percentage (lint%), one of the most important cotton yield components, but we still know little about the major genetic components underlying lint%. Here, we first linked the genetic loci containing MYB25-like_At and HD1_At to the fiberless seed trait of 'SL1-7-1' and found that MYB25-like_At and HD1_At were very lowly expressed in 'SL1-7-1' ovules during fiber initiation. We then dissected the genetic components involved in determination of lint% using segregating populations derived from crosses of fuzzless mutants and intermediate segregants with different lint%, which not only confirmed the HD1_At locus but identified the HD1_Dt locus as being the major genetic components contributing to fiber initiation and lint%. The segregating populations also allowed us to evaluate the relative contributions of MYB25-like_At, MYB25-like_Dt, HD1_At, and HD1_Dt to lint%. Haplotype analysis of an Upland cotton (Gossypium hirsutum) population with 723 accessions (including 81 fuzzless seed accessions) showed that lint% of the accessions with the LP allele (higher lint%) at MYB25-like_At, MYB25-like_Dt, or HD1_At was significantly higher than that with the lp allele (lower lint%). The lint% of the Upland cotton accessions with 3 or 4 LP alleles at MYB25-like and HD1 was significantly higher than that with 2 LP alleles. The results prompted us to propose a strategy for breeding high-yielding cotton varieties, i.e. pyramiding the LP alleles of MYB25-like and HD1 with new lint% LP alleles without negative impact on seed size and fiber quality.

Genetic basis of nectar guide trichome variation between bumblebee- and self-pollinated monkeyflowers (Mimulus): role of the MIXTA-like gene GUIDELESS

Nectar guide trichomes play crucial ecological roles in bee-pollinated flowers, as they serve as footholds and guides for foraging bees to access the floral rewards. However, the genetic basis of natural variation in nectar guide trichomes among species remains poorly understood. In this study, we performed genetic analysis of nectar guide trichome variation between two closely related monkeyflower (Mimulus) species, the bumblebee-pollinated Mimulus lewisii and self-pollinated M. parishii. We demonstrate that a MIXTA-like R2R3-MYB gene, GUIDELESS, is a major contributor to the nectar guide trichome length variation between the two species. The short-haired M. parishii carries a recessive allele due to non-synonymous substitutions in a highly conserved motif among MIXTA-like MYB proteins. Furthermore, our results suggest that besides GUIDELESS, additional loci encoding repressors of trichome elongation also contribute to the transition from bumblebee-pollination to selfing. Taken together, these results suggest that during a pollination syndrome switch, changes in seemingly complex traits such as nectar guide trichomes could have a relatively simple genetic basis, involving just a few genes of large effects.

Beyond skin-deep: targeting the plant surface for crop improvement

The above-ground plant surface is a well-adapted tissue layer that acts as an interface between the plant and its surrounding environment. As such, its primary role is to protect against desiccation and maintain the gaseous exchange required for photosynthesis. Further, this surface layer provides a barrier against pathogens and herbivory, while attracting pollinators and agents of seed dispersal. In the context of agriculture, the plant surface is strongly linked to post-harvest crop quality and yield. The epidermal layer contains several unique cell types adapted for these functions, while the non-lignified above-ground plant organs are covered by a hydrophobic cuticular membrane. This review aims to provide an overview of the latest understanding of the molecular mechanisms underlying crop cuticle and epidermal cell formation, with focus placed on genetic elements contributing towards quality, yield, drought tolerance, herbivory defence, pathogen resistance, pollinator attraction, and sterility, while highlighting the inter-relatedness of plant surface development and traits. Potential crop improvement strategies utilizing this knowledge are outlined in the context of the recent development of new breeding techniques.

Systematically and Comprehensively Understanding the Regulation of Cotton Fiber Initiation: A Review

Cotton fibers provide an important source of raw materials for the textile industry worldwide. Cotton fiber is a kind of single cell that differentiates from the epidermis of the ovule and provides a perfect research model for the differentiation and elongation of plant cells. Cotton fiber initiation is the first stage throughout the entire developmental process. The number of fiber cell initials on the seed ovule epidermis decides the final fiber yield. Thus, it is of great significance to clarify the mechanism underlying cotton fiber initiation. Fiber cell initiation is controlled by complex and interrelated regulatory networks. Plant phytohormones, transcription factors, sugar signals, small signal molecules, functional genes, non-coding RNAs, and histone modification play important roles during this process. Here, we not only summarize the different kinds of factors involved in fiber cell initiation but also discuss the mechanisms of these factors that act together to regulate cotton fiber initiation. Our aim is to synthesize a systematic and comprehensive review of different factors during fiber initiation that will provide the basics for further illustrating these mechanisms and offer theoretical guidance for improving fiber yield in future molecular breeding work.

GoSTR, a negative modulator of stem trichome formation in cotton

Trichomes, the outward projection of plant epidermal tissue, provide an effective defense against stress and insect pests. Although numerous genes have been identified to be involved in trichome development, the molecular mechanism for trichome cell fate determination is not well enunciated. Here, we reported GoSTR functions as a master repressor for stem trichome formation, which was isolated by map-based cloning based on a large F2 segregating population derived from a cross between TM-1 (pubescent stem) and J220 (smooth stem). Sequence alignment revealed a critical G-to-T point mutation in GoSTR's coding region that converted codon 2 from GCA (Alanine) to TCA (Serine). This mutation occurred between the majority of Gossypium hirsutum with pubescent stem (GG-haplotype) and G. barbadense with glabrous stem (TT-haplotype). Silencing of GoSTR in J220 and Hai7124 via virus-induced gene silencing resulted in the pubescent stems but no visible change in leaf trichomes, suggesting stem trichomes and leaf trichomes are genetically distinct. Yeast two-hybrid assay and luciferase complementation imaging assay showed GoSTR interacts with GoHD1 and GoHOX3, two key regulators of trichome development. Comparative transcriptomic analysis further indicated that many transcription factors such as GhMYB109, GhTTG1, and GhMYC1/GhDEL65 which function as positive regulators of trichomes were significantly upregulated in the stem from the GoSTR-silencing plant. Taken together, these results indicate that GoSTR functions as an essential negative modulator of stem trichomes and its transcripts will greatly repress trichome cell differentiation and growth. This study provided valuable insights for plant epidermal hair initiation and differentiation research.

Role of Actin Dynamics and GhACTIN1 Gene in Cotton Fiber Development: A Prototypical Cell for Study

Cotton crop is considered valuable for its fiber and seed oil. Cotton fiber is a single-celled outgrowth from the ovule epidermis, and it is a very dynamic cell for study. It has four distinct but overlapping developmental stages: initiation, elongation, secondary cell wall synthesis, and maturation. Among the various qualitative characteristics of cotton fiber, the important ones are the cotton fiber staple length, tensile strength, micronaire values, and fiber maturity. Actin dynamics are known to play an important role in fiber elongation and maturation. The current review gives an insight into the cotton fiber developmental stages, the qualitative traits associated with cotton fiber, and the set of genes involved in regulating these developmental stages and fiber traits. This review also highlights some prospects for how biotechnological approaches can improve cotton fiber quality.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton industry.

A novel regulatory complex mediated by Lanata (Ln) controls multicellular trichome formation in tomato

Trichomes that originate from plant aerial epidermis act as mechanical and chemical barriers against herbivores. Although several regulators have recently been identified, the regulatory pathway underlying multicellular trichome formation remains largely unknown in tomato. Here, we report a novel HD-ZIP IV transcription factor, Lanata (Ln), a missense mutation which caused the hairy phenotype. Biochemical analyses demonstrate that Ln separately interacts with two trichome regulators, Woolly (Wo) and Hair (H). Genetic and molecular evidence demonstrates that Ln directly regulates the expression of H. The interaction between Ln and Wo can increase trichome density by enhancing the expression of SlCycB2 and SlCycB3, which we previously showed are involved in tomato trichome formation. Furthermore, SlCycB2 represses the transactivation of the SlCycB3 gene by Ln and vice versa. Our findings provide new insights into the novel regulatory network controlling multicellular trichome formation in tomato.

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

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

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Genome-wide identification, expression and evolution analysis of OVATE family proteins in cotton (Gossypium spp.)

OVATE family proteins (OFPs) are plant-specific transcription factors with a conserved OVATE domain. Although OFPs have been reported to regulate many aspects of plant growth and development, little is known about their evolution, structure, and function in fiber development in cotton. In this study, 174 OFPs were identified from four species of Gossypium namely, G. hirsutum, G. barbadense, G. arboreum, and G. raimondii. These OFPs were grouped into 6 sub-families by using phylogenetic analysis, and members within the same sub-family had similar conserved motifs. Chromosomal localization revealed that OFPs are distributed in cotton genome unevenly. Gene structure analysis showed that most of OFPs were intronless. Moreover, Ka/Ks analysis exhibited that OFPs were gone through purifying selection processes during evolution. Multiple cis-acting elements were observed in promoter region of OFPs, which are responsive to light, phytohormone, biotic stresses, growth and developmental related cis-acting elements. In addition, OFPs play important role in fiber and ovule development. In conclusion, this study provides a systematic analysis of cotton OFPs and provides the foundation for further studies on biological functioning of cotton OFPs.

Tomato ARPC1 regulates trichome morphology and density and terpene biosynthesis

Based on transcriptomic analysis of wild-type and mutant tomato plants, ARPC1 was found to be important for trichome formation and development and it plays a key role in terpene synthesis. Trichomes are protruding epidermal cells in plant species. They function as the first defense layer against biotic and abiotic stresses. Despite the essential role of tomato trichomes in defense against herbivores, the understanding of their development is still incomplete. Therefore, the aim of this study was to identify genes involved in trichome formation and morphology and terpene synthesis, using transcriptomic techniques. To achieve this, we examined leaf morphology and compared the expression levels of some putative genes involved in trichome formation between wild-type (WT) and hairless-3 (hl-3) tomato mutant. The hl-3 plants displayed swollen and distorted trichomes and reduced trichome density (type I and IV) and terpene synthesis compared with that of the WT plants. Gene expression analysis showed that Actin-Related Protein Component1 (ARPC1) was expressed more highly in the WT than in the hl-3 mutant, indicating its critical role in trichome morphology and density. Additionally, the expression of MYC1 and several terpene synthase genes (TPS9, 12, 20), which are involved in type VI trichome initiation and terpene synthesis, was lower in the hl-3 mutant than in the WT plants. Moreover, transformation of the hl-3 mutant with WT ARPC1 restored normal trichome structure and density, and terpene synthesis. Structural and amino acid sequence analysis showed that there was a missplicing mutation in the hl-3 mutant, which was responsible for the abnormal trichome structure and density, and impaired terpene synthesis. Overall, the findings of this study demonstrated that ARPC1 is involved in regulating trichome structure and terpene synthesis in tomato.

Genome-Wide Identification and Characterization of CPR5 Genes in Gossypium Reveals Their Potential Role in Trichome Development

Trichomes protect plants against insects, microbes, herbivores, and abiotic damages and assist seed dispersal. The function of CPR5 genes have been found to be involved in the trichome development but the research on the underlying genetic and molecular mechanisms are extremely limited. Herein, genome wide identification and characterization of CPR5 genes was performed. In total, 26 CPR5 family members were identified in Gossypium species. Phylogenetic analysis, structural characteristics, and synteny analysis of CPR5s showed the conserved evolution relationships of CPR5. The promoter analysis of CPR5 genes revealed hormone, stress, and development-related cis-elements. Gene ontology (GO) enrichment analysis showed that the CPR5 genes were largely related to biological regulation, developmental process, multicellular organismal process. Protein-protein interaction analysis predicted several trichome development related proteins (SIM, LGO, and GRL) directly interacting with CPR5 genes. Further, nine putative Gossypium-miRNAs were also identified, targeting Gossypium CPR5 genes. RNA-Seq data of G. arboreum (with trichomes) and G. herbaceum (with no trichomes) was used to perform the co-expression network analysis. GheCPR5.1 was identified as a hub gene in a co-expression network analysis. RT-qPCR of GheCPR5.1 gene in different tissues suggests that this gene has higher expressions in the petiole and might be a key candidate involved in the trichome development. Virus induced gene silencing of GheCPR5.1 (Ghe02G17590) confirms its role in trichome development and elongation. Current results provide proofs of the possible role of CPR5 genes and provide preliminary information for further studies of GheCPR5.1 functions in trichome development.

Molecular Mechanisms of Plant Trichome Development

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

Cotton Breeding in Australia: Meeting the Challenges of the 21st Century

The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program is the sole breeding effort for cotton in Australia, developing high performing cultivars for the local industry which is worth∼AU$3 billion per annum. The program is supported by Cotton Breeding Australia, a Joint Venture between CSIRO and the program's commercial partner, Cotton Seed Distributors Ltd. (CSD). While the Australian industry is the focus, CSIRO cultivars have global impact in North America, South America, and Europe. The program is unique compared with many other public and commercial breeding programs because it focuses on diverse and integrated research with commercial outcomes. It represents the full research pipeline, supporting extensive long-term fundamental molecular research; native and genetically modified (GM) trait development; germplasm enhancement focused on yield and fiber quality improvements; integration of third-party GM traits; all culminating in the release of new commercial cultivars. This review presents evidence of past breeding successes and outlines current breeding efforts, in the areas of yield and fiber quality improvement, as well as the development of germplasm that is resistant to pests, diseases and abiotic stressors. The success of the program is based on the development of superior germplasm largely through field phenotyping, together with strong commercial partnerships with CSD and Bayer CropScience. These relationships assist in having a shared focus and ensuring commercial impact is maintained, while also providing access to markets, traits, and technology. The historical successes, current foci and future requirements of the CSIRO cotton breeding program have been used to develop a framework designed to augment our breeding system for the future. This will focus on utilizing emerging technologies from the genome to phenome, as well as a panomics approach with data management and integration to develop, test and incorporate new technologies into a breeding program. In addition to streamlining the breeding pipeline for increased genetic gain, this technology will increase the speed of trait and marker identification for use in genome editing, genomic selection and molecular assisted breeding, ultimately producing novel germplasm that will meet the coming challenges of the 21st Century.

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.

Multi-Dimensional Molecular Regulation of Trichome Development in Arabidopsis and Cotton

Plant trichomes are specialized epidermal cells that are widely distributed on plant aerial tissues. The initiation and progression of trichomes are controlled in a coordinated sequence of multiple molecular events. During the past decade, major breakthroughs in the molecular understanding of trichome development were achieved through the characterization of various trichomes defective mutants and trichome-associated genes, which revealed a highly complex molecular regulatory network underlying plant trichome development. This review focuses on the recent millstone in plant trichomes research obtained using genetic and molecular studies, as well as 'omics' analyses in model plant Arabidopsis and fiber crop cotton. In particular, we discuss the latest understanding and insights into the underlying molecular mechanisms of trichomes formation at multiple dimensions, including at the chromatin, transcriptional, post-transcriptional, and post-translational levels. We summarize that the integration of multi-dimensional trichome-associated genes will enable us to systematically understand the molecular regulation network that landscapes the development of the plant trichomes. These advances will enable us to address the unresolved questions regarding the molecular crosstalk that coordinate concurrent and ordered the changes in cotton fiber initiation and progression, together with their possible implications for genetic improvement of cotton fiber.

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.

Two R2R3-MYB genes cooperatively control trichome development and cuticular wax biosynthesis in Prunus persica

The fruit surface has an enormous impact on the external appearance and postharvest shelf-life of fruit. Here, we report two functionally redundant genes, PpMYB25 and PpMYB26, involved in regulation of fruit skin texture in peach. PpMYB25 can activate transcription of PpMYB26 and they both induce trichome development and cuticular wax accumulation, resulting in peach fruit with a fuzzy and dull appearance. By contrast, nonfunctional mutation of PpMYB25 caused by an insertional retrotransposon in the last exon in nectarine fails to activate transcription of PpMYB26, resulting in nectarine fruit with a smooth and shiny appearance due to loss of trichome initiation and decreased cuticular wax accumulation. Secondary cell wall biosynthesis in peach fruit pubescence is controlled by a transcriptional regulatory network, including the master regulator PpNAC43 and its downstream MYB transcription factors such as PpMYB42, PpMYB46 and PpMYB83. Our results show that PpMYB25 and PpMYB26 coordinately regulate fruit pubescence and cuticular wax accumulation and their simultaneous perturbation results in the origin of nectarine, which is botanically classified as a subspecies of peach.

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.

RNA-Seq with a novel glabrous-ZM24fl reveals some key lncRNAs and the associated targets in fiber initiation of cotton

BACKGROUND

Cotton fiber is an important natural resource for textile industry and an excellent model for cell biology study. Application of glabrous mutant cotton and high-throughput sequencing facilitates the identification of key genes and pathways for fiber development and cell differentiation and elongation. LncRNA is a type of ncRNA with more than 200 nt in length and functions in the ways of chromatin modification, transcriptional and post-transcriptional modification, and so on. However, the detailed lncRNA and associated mechanisms for fiber initiation are still unclear in cotton.

RESULTS

In this study, we used a novel glabrous mutant ZM24fl, which is endowed with higher somatic embryogenesis, and functions as an ideal receptor for cotton genetic transformation. Combined with the high-throughput sequencing, fatty acid pathway and some transcription factors such as MYB, ERF and bHLH families were identified the important roles in fiber initiation; furthermore, 3,288 lncRNAs were identified, and some differentially expressed lncRNAs were also analyzed. From the comparisons of ZM24_0 DPA vs ZM24_-2 DPA and fl_0 DPA vs ZM24_0 DPA, one common lncRNA MSTRG 2723.1 was found that function upstream of fatty acid metabolism, MBY25-mediating pathway, and pectin metabolism to regulate fiber initiation. In addition, other lncRNAs MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1 were also showed potential important roles in fiber development; and the co-expression analysis between lncRNAs and targets showed the distinct models of different lncRNAs and complicated interaction between lncRNAs in fiber development of cotton.

CONCLUSIONS

From the above results, a key lncRNA MSTRG 2723.1 was identified that might mediate some key genes transcription of fatty acid metabolism, MYB25-mediating pathway, and pectin metabolism to regulate fiber initiation of ZM24 cultivar. Co-expression analysis implied that some other important lncRNAs (e.g., MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1) were also showed the different regulatory model and interaction between them, which proposes some valuable clues for the lncRNAs associated mechanisms in fiber development.

An Unexpected Regulatory Sequence from Rho-Related GTPase6 Confers Fiber-Specific Expression in Upland Cotton

Cotton fibers, single seed trichomes derived from ovule epidermal cells, are the major source of global textile fibers. Fiber-specific promoters are desirable to study gene function and to modify fiber properties during fiber development. Here, we revealed that Rho-related GTPase6 (GhROP6) was expressed preferentially in developing fibers. A 1240 bp regulatory region of GhROP6, which contains a short upstream regulatory sequence, the first exon, and the partial first intron, was unexpectedly isolated and introduced into transgenic cotton for analyzing promoter activity. The promoter of GhROP6 (proChROP6) conferred a specific expression in ovule surface, but not in the other floral organs and vegetative tissues. Reverse transcription PCR analysis indicated that proGhROP6 directed full-length transcription of the fused ß-glucuronidase (GUS) gene. Further investigation of GUS staining showed that proChROP6 regulated gene expression in fibers and ovule epidermis from fiber initiation to cell elongation stages. The preferential activity was enriched in fiber cells after anthesis and reached to peak on flowering days. By comparison, proGhROP6 was a mild promoter with approximately one-twenty-fifth of the strength of the constitutive promoter CaMV35S. The promoter responded to high-dosage treatments of auxin, gibberellin and salicylic acid and slightly reduced GUS activity under the in vitro treatment. Collectively, our data suggest that the GhROP6 promoter has excellent activity in initiating fibers and has potential for bioengineering of cotton fibers.

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.

An HD-ZIP-MYB complex regulates glandular secretory trichome initiation in Artemisia annua

Plant glandular secretory trichomes (GSTs) produce various specialized metabolites. Increasing GST density represents a strategy to enhance the yield of these chemicals; however, the gene regulatory network that controls GST initiation remains unclear. In a previous study of Artemisia annua L., we found that a HD-ZIP IV transcription factor, AaHD1, promotes GST initiation by directly regulating AaGSW2. Here, we identified two AaHD1-interacting transcription factors, namely AaMIXTA-like 2 (AaMYB16) and AaMYB5. Through the generation and characterization of transgenic plants, we found that AaMYB16 is a positive regulator of GST initiation, whereas AaMYB5 has the opposite effect. Notably, neither of them regulates GST formation independently. Rather, they act competitively, by interacting and modulating AaHD1 promoter binding activity. Additionally, the phytohormone jasmonic acid (JA) was shown to be associated with the AaHD1-AaMYB16/AaMYB5 regulatory network through transcriptional regulation via a JASMONATE-ZIM DOMAIN (JAZ) protein repressor. These results bring new insights into the mechanism of GST initiation through regulatory complexes, which appear to have similar functions in a range of vascular plant taxa.

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.

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

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

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.

Looking into 'hair tonics' for cotton fiber initiation

Cotton fiber is the most important source of cellulose for the global textile industry. These hair-like single-celled trichomes develop from ovule epidermis. They are classified into long spinnable lint and short fuzz. A key objective in the cotton industry is to breed elite cultivars with fuzzless seeds carrying high lint yield. Molecular basis underlying lint and fuzz initiation remains obscure. Recent studies indicate fiber initiation is under the control of MYB-bHLH-WDR (MBW) transcription factor complex. Based on molecular genetic studies and gene expression patterns linking fiber phenotypes, we propose that specific but different sets of MBW genes are required to precisely regulate the initiation of the lint and fuzz fibers. Emerging evidence further points to sugar signaling as a 'hair-tonic' to boost fiber initiation through interaction with MBW complex and auxin signaling. An integrative model is provided as a conceptual framework for future studies to dissect the molecular network responsible for cotton fiber initiation.

Elucidation of sequence polymorphism in fuzzless-seed cotton lines

Most commercially produced cotton cultivars have two types of fibers on the seed coat, short fuzz and long lint. Lint fiber is used in the textile industry, while fuzz is considered an undesirable trait. Both types of fibers are believed to be controlled by the same regulators; however, their mechanisms of actions are still obscure. Cotton fiber mutants provide an excellent system to study the genes that regulate fiber development. Here we described four uncharacterized and three previously reported cotton mutants with fuzzless seed phenotypes. To evaluate whether or not the genes previously associated with fuzzless seed phenotypes have mutations we sequenced whole genomic DNA of seven mutants and wild type varieties. We identified multiple polymorphic changes among the tested genes. Non-synonymous SNPs in the coding region of the MML3-A gene was common in the six mutant lines tested in this study, showing both dominant and recessive fuzzless phenotypes. We have mapped the locus of the causative mutation for one of the uncharacterized fuzzless lines using an F₂ population that originated from a cross between the dominant fuzzless mutant and a wild type. Further, we have clarified the current knowledge about the causative n₂ mutations by analyzing the sequence data and previously reported mapping data. The key genes and possible mechanisms of fiber differentiation are discussed in this study.

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.

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.

Genetic Identification and Transcriptome Analysis of Lintless and Fuzzless Traits in Gossypium arboreum L

Cotton fibres, as single cells arising from the seed coat, can be classified as lint and fuzz according to their final length. Gossypium arboreum is a cultivated diploid cotton species and a potential donor of the A subgenome of the more widely grown tetraploid cottons. In this study, we performed genetic studies on one lintless and seven fuzzless G. arboreum accessions. Through association and genetic linkage analyses, a recessive locus on Chr06 containing GaHD-1 was found to be the likely gene underlying the lintless trait. GaHD-1 carried a mutation at a splicing acceptor site that resulted in alternative splicing and a deletion of 247 amino acid from the protein. The regions containing GaGIR1 and GaMYB25-like were found to be associated with fuzz development in G. arboreum, with the former being the major contributor. Comparative transcriptome analyses using 0-5 days post-anthesis (dpa) ovules from lintless, fuzzless, and normal fuzzy seed G. arboreum accessions revealed gene modules and hub genes potentially important for lint and fuzz initiation and development. Three significant modules and 26 hub genes associated with lint fibre initiation were detected by weighted gene co-expression network analysis. Similar analyses identified three vital modules and 10 hub genes to be associated with fuzz development. The findings in this study contribute to understanding the complex molecular mechanism(s) regulating fibre initiation and development and indicate that G. arboreum may have fibre developmental pathways different from tetraploid cotton. It also provides candidate genes for further investigation into modifying fibre development in G. arboreum.

Updates on molecular mechanisms in the development of branched trichome in Arabidopsis and nonbranched in cotton

Trichomes are specialized epidermal cells and a vital plant organ that protect plants from various harms and provide valuable resources for plant development and use. Some key genes related to trichomes have been identified in the model plant Arabidopsis thaliana through glabrous mutants and gene cloning, and the hub MYB-bHLH-WD40, consisting of several factors including GLABRA1 (GL1), GL3, TRANSPARENT TESTA GLABRA1 (TTG1), and ENHANCER OF GLABRA3 (EGL3), has been established. Subsequently, some upstream transcription factors, phytohormones and epigenetic modification factors have also been studied in depth. In cotton, a very important fibre and oil crop globally, in addition to the key MYB-like factors, more important regulators and potential molecular mechanisms (e.g. epigenetic modifiers, distinct metabolic pathways) are being exploited during different fibre developmental stages. This occurs due to increased cotton research, resulting in the discovery of more complex regulation mechanisms from the allotetraploid genome of cotton. In addition, some conservative as well as specific mediators are involved in trichome development in other species. This study summarizes molecular mechanisms in trichome development and provides a detailed comparison of the similarities and differences between Arabidopsis and cotton, analyses the possible reasons for the discrepancy in identification of regulators, and raises future questions and foci for understanding trichome development in more detail.

Preferential insertion of a Ty1 LTR-retrotransposon into the A sub-genome's HD1 gene significantly correlated with the reduction in stem trichomes of tetraploid cotton

Stem trichomes and seed fibers originate from epidermal cells and partially share a regulatory pathway at the molecular level. In Gossypium barbadense, two insertions of a Ty1 long-terminal repeat-retrotransposon [transposable element TE1 and TE2] in a homeodomain-leucine zipper gene (HD1) result in glabrous stems. The primers used to identify the TE insertions in G. barbadense were applied to screen for the same events in 81 modern G. hirsutum varieties and 31 wild races. Three wild races were found carrying the same TEs as G. barbadense. However, the TE insertions in two of these wild races occurred at different sites (4th exon), therefore, named TE3, while the TE in the other wild race occurred at the same site as TE2. An RNA sequencing and qRT-PCR analysis indicated that the loss of HD1 function was caused by the TE insertion. Genetic mapping revealed a strong association between glabrous stems and TE3 insertions, confirming that HD1 is a critical gene for stem trichome initiation in G. hirsutum, as in G. barbadense. Using the long-terminal repeat sequence as a query to search against the Texas Marker-1 reference genome sequence, we found that the TE occurred after tetraploid cotton formation and evolved at different rates in G. hirsutum and G. barbadense. Interestingly, at least three independent insertion events of the same retrotransposon occurred preferentially in the A sub-genome's HD1 gene, but not in the D sub-genome of G. hirsutum or G. barbadense, suggesting that an unknown TE insertion mechanism and resultant gene function changes may have hastened cotton speciation.

Transcriptome analysis reveals differences in the mechanisms of fiber initiation and elongation between long- and short-fiber cotton (Gossypium hirsutum L.) lines

Background

Improving the yield and fiber quality of upland cotton is a goal of plant breeders. However, increasing the yield and quality of cotton fibers is becoming more urgent. While the growing human population needs more cotton fiber, climate change is reducing the amount of land on which cotton can be planted, or making it difficult to ensure that water and other resources will be available in optimal quantities. The most logical means of improving yield and quality is understanding and manipulating the genes involved. Here, we used comparative transcriptomics to explore differences in gene expression between long- and short-fiber cotton lines to identify candidate genes useful for cotton improvement.

Results

Light and electron microscopy revealed that the initial fiber density was significantly greater in our short-fiber group (SFG) than in our long-fiber group (LFG). Compared with the SFG fibers, the LFG fibers were longer at all developmental stages. Comparison of the LFG and SFG transcriptomes revealed a total of 3538 differentially expressed genes (DEGs). Notably, at all three developmental stages examined, two expression patterns, consistently downregulated (profile 0) and consistently upregulated (profile 7), were identified, and both were significantly enriched in the SFG and LFG. Twenty-two DEGs known to be involved in fiber initiation were detected in profile 0, while 31 DEGs involved in fiber elongation were detected in profile 7. Functional annotation suggested that these DEGs, which included ERF1, TUA2, TUB1, and PER64, affect fiber elongation by participating in the ethylene response, microtubule synthesis, and/or the peroxidase (POD) catalytic pathway. qRT-PCR was used to confirm the RNA sequencing results for select genes.

Conclusions

A comparison of SFG and LFG transcription profiles revealed modest but important differences in gene expression between the groups. Notably, our results confirm those of previous studies suggesting that genes involved in ethylene, tubulin, and POD pathways play important roles in fiber development. The 22 consistently downregulated DEGs involved in fiber initiation and the 31 consistently upregulated genes involved in fiber elongation are seemingly good candidate genes for improving fiber initiation and elongation in cotton.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5986-5) contains supplementary material, which is available to authorized users.

Cotton fiber elongation requires the transcription factor GhMYB212 to regulate sucrose transportation into expanding fibers

Cotton is white gold across the globe and composed of fiber cells derived from the outer integument of cotton ovules. Fiber elongation uses sucrose as a direct carbon source. The molecular mechanism transcriptionally controlling sucrose transport from ovules into the elongating fibers remains elusive. In this study the involvement of GhMYB212 in the regulation of sucrose transportion into expanding fibers was investigated. GhMYB212 RNAi plants (GhMYB212i) accumulated less sucrose and glucose in developing fibers, and had shorter fibers and a lower lint index. RNA-seq and protein-DNA binding assays revealed that GhMYB212 was closely linked to the pathways of sucrose and starch transportation and metabolism, directly controling the expression of a sucrose transporter gene GhSWEET12. GhSWEET12 RNAi plants (GhSWEET12i) possessed similar fiber phenotypes to those of GhMYB212i. Exogenous sucrose supplementation in ovule cultures did not rescue the shorter fiber phenotype of GhMYB212i and GhSWEET12i. This finding supported the idea that the attenuated rate of sucrose transport from the outer seed coat into the fibers is responsible for the retardation of fiber elongation. Current investigations support the idea that GhMYB212 functions as the main regulator of fiber elongation by controlling the expression of GhSWEET12, and therefore it is important to study cell expansion and sugar transportation during seed development.

Role of GhHDA5 in H3K9 deacetylation and fiber initiation in Gossypium hirsutum

Cotton fibers are single-celled trichomes that initiate from the epidermal cells of the ovules at or before anthesis. Here, we identified that the histone deacetylase (HDAC) activity is essential for proper cotton fiber initiation. We further identified 15 HDACs homoeologs in each of the A- and D-subgenomes of Gossypium hirsutum. Few of these HDAC homoeologs expressed preferentially during the early stages of fiber development [-1, 0 and 6 days post-anthesis (DPA)]. Among them, GhHDA5 expressed significantly at the time of fiber initiation (-1 and 0 DPA). The in vitro assay for HDAC activity indicated that GhHDA5 primarily deacetylates H3K9 acetylation marks. Moreover, the reduced expression of GhHDA5 also suppresses fiber initiation and lint yield in the RNA interference (RNAi) lines. The 0 DPA ovules of GhHDA5RNAi lines also showed alterations in reactive oxygen species homeostasis and elevated autophagic cell death in the developing fibers. The differentially expressed genes (DEGs) identified through RNA-seq of RNAi line (DEP12) and their pathway analysis showed that GhHDA5 modulates expression of many stress and development-related genes involved in fiber development. The reduced expression of GhHDA5 in the RNAi lines also resulted in H3K9 hyper-acetylation on the promoter region of few DEGs assessed by chromatin immunoprecipitation assay. The positively co-expressed genes with GhHDA5 showed cumulative higher expression during fiber initiation, and gene ontology annotation suggests their involvement in fiber development. Furthermore, the predicted protein interaction network in the positively co-expressed genes indicates HDA5 modulates fiber initiation-specific gene expression through a complex involving reported repressors.

Genome-wide classification, evolutionary analysis and gene expression patterns of the kinome in Gossypium

The protein kinase (PK, kinome) family is one of the largest families in plants and regulates almost all aspects of plant processes, including plant development and stress responses. Despite their important functions, comprehensive functional classification, evolutionary analysis and expression patterns of the cotton PK gene family has yet to be performed on PK genes. In this study, we identified the cotton kinomes in the Gossypium raimondii, Gossypium arboretum, Gossypium hirsutum and Gossypium barbadense genomes and classified them into 7 groups and 122-24 subfamilies using software HMMER v3.0 scanning and neighbor-joining (NJ) phylogenetic analysis. Some conserved exon-intron structures were identified not only in cotton species but also in primitive plants, ferns and moss, suggesting the significant function and ancient origination of these PK genes. Collinearity analysis revealed that 16.6 million years ago (Mya) cotton-specific whole genome duplication (WGD) events may have played a partial role in the expansion of the cotton kinomes, whereas tandem duplication (TD) events mainly contributed to the expansion of the cotton RLK group. Synteny analysis revealed that tetraploidization of G. hirsutum and G. barbadense contributed to the expansion of G. hirsutum and G. barbadense PKs. Global expression analysis of cotton PKs revealed stress-specific and fiber development-related expression patterns, suggesting that many cotton PKs might be involved in the regulation of the stress response and fiber development processes. This study provides foundational information for further studies on the evolution and molecular function of cotton PKs.

SSR-based association mapping of fiber quality in upland cotton using an eight-way MAGIC population

The quality of fiber is significant in the upland cotton industry. As complex quantitative traits, fiber quality traits are worth studying at a genetic level. To investigate the genetic architecture of fiber quality traits, we conducted an association analysis using a multi-parent advanced generation inter-cross (MAGIC) population developed from eight parents and comprised of 960 lines. The reliable phenotypic data for six major fiber traits of the MAGIC population were collected from five environments in three locations. Phenotypic analysis showed that the MAGIC lines have a wider variation amplitude and coefficient than the founders. A total of 284 polymorphic SSR markers among eight parents screened from a high-density genetic map were used to genotype the MAGIC population. The MAGIC population showed abundant genetic variation and fast linkage disequilibrium (LD) decay (0.76 cM, r² > 0.1), which revealed the advantages of high efficiency and power in QTL exploration. Association mapping via a mixed linear model identified 52 significant loci associated with six fiber quality traits; 14 of them were mapped in reported QTL regions with fiber-related or other agronomic traits. Nine markers demonstrated the pleiotropism that controls more than two fiber traits. Furthermore, two SSR markers, BNL1231 and BNL3452, were authenticated as hotspots that were mapped with multi-traits. In addition, we provided candidate regions and screened six candidate genes for identified loci according to the LD decay distance. Our results provide valuable QTL for further genetic mapping and will facilitate marker-based breeding for fiber quality in cotton.

Genome-wide identification and characterization of TCP transcription factor genes in upland cotton (Gossypium hirsutum)

TCP proteins are plant-specific transcription factors (TFs), and perform a variety of physiological functions in plant growth and development. In this study, 74 non-redundant TCP genes were identified in upland cotton (Gossypium hirsutum L.) genome. Cotton TCP family can be classified into two classes (class I and class II) that can be further divided into 11 types (groups) based on their motif composition. Quantitative RT-PCR analysis indicated that GhTCPs display different expression patterns in cotton tissues. The majority of these genes are preferentially or specifically expressed in cotton leaves, while some GhTCP genes are highly expressed in initiating fibers and/or elongating fibers of cotton. Yeast two-hybrid results indicated that GhTCPs can interact with each other to form homodimers or heterodimers. In addition, GhTCP14a and GhTCP22 can interact with some transcription factors which are involved in fiber development. These results lay solid foundation for further study on the functions of TCP genes during cotton fiber development.

The cotton β-galactosyltransferase 1 (GalT1) that galactosylates arabinogalactan proteins participates in controlling fiber development

Arabinogalactan proteins (AGPs) are highly glycosylated proteins that play pivotal roles in diverse developmental processes in plants. Type-II AG glycans, mostly O-linked to the hydroxyproline residues of the protein backbone, account for up to 95% w/w of the AGP, but their functions are still largely unclear. Cotton fibers are extremely elongated single-cell trichomes on the seed epidermis; however, little is known of the molecular basis governing the regulation of fiber cell development. Here, we characterized the role of a CAZy glycosyltransferase 31 (GT31) family member, GhGalT1, in cotton fiber development. The fiber length of the transgenic cotton overexpressing GhGalT1 was shorter than that of the wild type, whereas in the GhGalT1-silenced lines there was a notable increase in fiber length compared with wild type. The carbohydrate moieties of AGPs were altered in fibers of GhGalT1 transgenic cotton. The galactose: arabinose ratio of AG glycans was higher in GhGalT1 overexpression fibers, but was lower in GhGalT1-silenced lines, compared with that in the wild type. Overexpression of GhGalT1 upregulates transcript levels of a broad range of cell wall-related genes, especially the fasciclin-like AGP (FLA) backbone genes. An enzyme activity assay demonstrated that GhGalT1 is a β-1,3-galactosyltransferase (β-1,3-GalT) involved in biosynthesis of the β-1,3-galactan backbone of the type-II AG glycans of AGPs. We also show that GhGalT1 can form homo- and heterodimers with other cotton GT31 family members to facilitate AG glycan assembly of AGPs. Thus, our data demonstrate that GhGalT1 influences cotton fiber development via controlling the glycosylation of AGPs, especially FLAs.

Genome-wide analysis of the HD-ZIP IV transcription factor family in Gossypium arboreum and GaHDG11 involved in osmotic tolerance in transgenic Arabidopsis

HD-ZIP IV proteins belong to the homeodomain-leucine zipper (HD-ZIP) transcription factor family and are involved in trichome development and drought stress in plants. Although some functions of the HD-ZIP IV group are well understood in Arabidopsis, little is known about their function in cotton. In this study, HD-ZIP genes were identified from three Gossypium species (G. arboreum, G. raimondii and G. hirsutum) and clustered into four families (HD-ZIP I, II, III and IV) to separate HD-ZIP IV from the other three families. Systematic analyses of phylogeny, gene structure, conserved domains, and expression profiles in different plant tissues and the expression patterns under osmotic stress in leaves were further conducted in G. arboreum. More importantly, ectopic overexpression of GaHDG11, a representative of the HD-ZIP IV family, confers enhanced osmotic tolerance in transgenic Arabidopsis plants, possibly due to elongated primary root length, lower water loss rates, high osmoprotectant proline levels, significant levels of antioxidants CAT, and/or SOD enzyme activity with reduced levels of MDA. Taken together, these observations may lay the foundation for future functional analysis of cotton HD-ZIP IV genes to unravel their biological roles in cotton.

HOMEODOMAIN PROTEIN 1 is required for jasmonate-mediated glandular trichome initiation in Artemisia annua

Glandular trichomes are generally considered biofactories that produce valuable chemicals. Increasing glandular trichome density is a very suitable way to improve the productivity of these valuable metabolites, but little is known about the regulation of glandular trichome formation. Phytohormone jasmonate (JA) promotes glandular trichome initiation in various plants, but its mechanism is also unknown. By searching transcription factors regulated by JA in Artemisia annua, we identified a novel homeodomain-leucine zipper transcription factor, HOMEODOMAIN PROTEIN 1 (AaHD1), which positively controls both glandular and nonglandular trichome initiations. Overexpression of AaHD1 in A. annua significantly increased glandular trichome density without harming plant growth. Consequently, the artemisinin content was improved. AaHD1 interacts with A. annua jasmonate ZIM-domain 8 (AaJAZ8), which is a repressor of JA, thereby resulting in decreased transcriptional activity. AaHD1 knockdown lines show decreased sensitivity to JA on glandular trichome initiation, which indicates that AaHD1 plays an important role in JA-mediated glandular trichome initiation. We identified a new transcription factor that promotes A. annua glandular trichome initiation and revealed a novel molecular mechanism by which a homeodomain protein transduces JA signal to promote glandular trichome initiation. Our results also suggested a connection between glandular and nonglandular trichome formations.

Small interfering RNAs from bidirectional transcripts of GhMML3_A12 regulate cotton fiber development

Natural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development. Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development. The extremely low expression of GhMML3_A12 in N1 is associated with NAT production, driven by its 3' antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton.

Fine mapping and RNA-Seq unravels candidate genes for a major QTL controlling multiple fiber quality traits at the T1 region in upland cotton

BACKGROUND

Improving fiber quality is a major challenge in cotton breeding, since the molecular basis of fiber quality traits is poorly understood. Fine mapping and candidate gene prediction of quantitative trait loci (QTL) controlling cotton fiber quality traits can help to elucidate the molecular basis of fiber quality. In our previous studies, one major QTL controlling multiple fiber quality traits was identified near the T1 locus on chromosome 6 in Upland cotton.

RESULTS

To finely map this major QTL, the F2 population with 6975 individuals was established from a cross between Yumian 1 and a recombinant inbred line (RIL118) selected from a recombinant inbred line population (T586 × Yumian 1). The QTL was mapped to a 0.28-cM interval between markers HAU2119 and SWU2302. The QTL explained 54.7 % (LOD = 222.3), 40.5 % (LOD = 145.0), 50.0 % (LOD = 194.3) and 30.1 % (LOD = 100.4) of phenotypic variation with additive effects of 2.78, -0.43, 2.92 and 1.90 units for fiber length, micronaire, strength and uniformity, respectively. The QTL region corresponded to a 2.7-Mb interval on chromosome 10 in the G. raimondii genome sequence and a 5.3-Mb interval on chromosome A06 in G. hirsutum. The fiber of Yumian 1 was much longer than that of RIL118 from 3 DPA to 7 DPA. RNA-Seq of ovules at 0 DPA and fibers at 5 DPA from Yumian 1 and RIL118 showed four genes in the QTL region of the G. raimondii genome to be extremely differentially expressed. RT-PCR analysis showed three genes in the QTL region of the G. hirsutum genome to behave similarly.

CONCLUSIONS

This study mapped a major QTL influencing four fiber quality traits to a 0.28-cM interval and identified three candidate genes by RNA-Seq and RT-PCR analysis. Integration of fine mapping and RNA-Seq is a powerful strategy to uncover candidates for QTL in large genomes.

RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.)

RNA interference (RNAi), is a powerful new technology in the discovery of genetic sequence functions, and has become a valuable tool for functional genomics of cotton (Gossypium sp.). The rapid adoption of RNAi has replaced previous antisense technology. RNAi has aided in the discovery of function and biological roles of many key cotton genes involved in fiber development, fertility and somatic embryogenesis, resistance to important biotic and abiotic stresses, and oil and seed quality improvements as well as the key agronomic traits including yield and maturity. Here, we have comparatively reviewed seminal research efforts in previously used antisense approaches and currently applied breakthrough RNAi studies in cotton, analyzing developed RNAi methodologies, achievements, limitations, and future needs in functional characterizations of cotton genes. We also highlighted needed efforts in the development of RNAi-based cotton cultivars, and their safety and risk assessment, small and large-scale field trials, and commercialization.