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

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

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.

Genome-Wide and Expression Pattern Analysis of the DVL Gene Family Reveals GhM_A05G1032 Is Involved in Fuzz Development in G. hirsutum

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

Applications of Virus-Induced Gene Silencing in Cotton

Virus-induced gene silencing (VIGS) is an RNA-mediated reverse genetics technique that has become an effective tool to investigate gene function in plants. Cotton is one of the most important economic crops globally. In the past decade, VIGS has been successfully applied in cotton functional genomic studies, including those examining abiotic and biotic stress responses and vegetative and reproductive development. This article summarizes the traditional vectors used in the cotton VIGS system, the visible markers used for endogenous gene silencing, the applications of VIGS in cotton functional genomics, and the limitations of VIGS and how they can be addressed in cotton.

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.

Activation of Gossypium hirsutum ACS6 Facilitates Fiber Development by Improving Sucrose Metabolism and Transport

Cotton fiber yield depends on the density of fiber cell initials that form on the ovule epidermis. Fiber initiation is triggered by MYB-MIXTA-like transcription factors (GhMMLs) and requires a sucrose supply. Ethylene or its precursor ACC (1-aminocyclopropane-1-carboxylic acid) is suggested to affect fiber yield. The Gossypium hirsutum (L.) genome contains 35 ACS genes (GhACS) encoding ACC synthases. Here, we explored the role of a GhACS family member in the regulation of fiber initiation. Expression analyses showed that the GhACS6.3 gene pair was specifically expressed in the ovules during fiber initiation (3 days before anthesis to 5 days post anthesis, -3 to 5 DPA), especially at -3 DPA, whereas other GhACS genes were expressed at very low or undetectable levels. The expression profile of GhACS6.3 during fiber initial development was confirmed by qRT-PCR analysis. Transgenic lines overexpressing GhACS6.3 (GhACS6.3-OE) showed increased ACC accumulation in ovules, which promoted the formation of fiber initials and fiber yield components. This was accompanied by increased transcript levels of GhMML3 and increased transcript levels of genes encoding sucrose transporters and sucrose synthase. These findings imply that GhACS6.3 activation is required for fiber initial development. Our results lay the foundation for further research on increasing cotton fiber production.

Advances in understanding epigenetic regulation of plant trichome development: a comprehensive review

Plant growth and development are controlled by a complex gene regulatory network, which is currently a focal point of research. It has been established that epigenetic factors play a crucial role in plant growth. Trichomes, specialized appendages that arise from epidermal cells, are of great significance in plant growth and development. As a model system for studying plant development, trichomes possess both commercial and research value. Epigenetic regulation has only recently been implicated in the development of trichomes in a limited number of studies, and microRNA-mediated post-transcriptional regulation appears to dominate in this context. In light of this, we have conducted a review that explores the interplay between epigenetic regulations and the formation of plant trichomes, building upon existing knowledge of hormones and transcription factors in trichome development. Through this review, we aim to deepen our understanding of the regulatory mechanisms underlying trichome formation and shed light on future avenues of research in the field of epigenetics as it pertains to epidermal hair growth.

Global identification of natural antisense transcripts in Gossypium hirsutum and Gossypium barbadense under chilling stress

Natural antisense transcripts (NATs) in model plants have been recognized as important regulators of gene expression under abiotic stresses. However, the functional roles of NATs in crops under low temperature are still unclear. Here, we identified 815 and 689 NATs from leaves of Gossypium hirsutum and G. barbadense under chilling stress. Among those, 224 NATs were identified as interspecific homologs between the two species. The correlation coefficients for expression of NATs and their cognate sense genes (CSG) were 0.43 and 0.37 in G. hirsutum and G. barbadense, respectively. Furthermore, expression of interspecific NATs and CSGs alike was highly consistent under chilling stress with correlation coefficients of 0.90-0.91. Four cold-associated NATs were selected for functional validation using virus-induced gene silencing (VIGS). Our results suggest that CAN1 engage in the molecular regulation of chilling stress by regulating SnRK2.8 expression. This highly conserved NAT have valuable potential for applications in breeding cold-tolerant cotton.

An insight into the roles of regulatory ncRNAs in plants: An abiotic stress and developmental perspective

Different environmental cues lead to changes in physiology, biochemistry and molecular status of plant's growth. Till date, various genes have been accounted for their role in regulating plant development and response to abiotic stress. Excluding genes that code for a functional protein in a cell, a large chunk of the eukaryotic transcriptome consists of non-coding RNAs (ncRNAs) which lack protein coding capacity but are still functional. Recent advancements in Next Generation Sequencing (NGS) technology have led to the unearthing of different types of small and large non-coding RNAs in plants. Non-coding RNAs are broadly categorised into housekeeping ncRNAs and regulatory ncRNAs which work at transcriptional, post-transcriptional and epigenetic levels. Diverse ncRNAs play different regulatory roles in nearly all biological processes including growth, development and response to changing environments. This response can be perceived and counteracted by plants using diverse evolutionarily conserved ncRNAs like miRNAs, siRNAs and lncRNAs to participate in complex molecular regimes by activating gene-ncRNA-mRNA regulatory modules to perform the downstream function. Here, we review the current understanding with a focus on recent advancements in the functional studies of the regulatory ncRNAs at the nexus of abiotic stresses and development. Also, the potential roles of ncRNAs in imparting abiotic stress tolerance and yield improvement in crop plants are also discussed with their future prospects.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Overexpression of cotton Trihelix transcription factor GhGT-3b_A04 enhances resistance to Verticillium dahliae and affects plant growth in Arabidopsis thaliana

Verticillium wilt is a soil-borne fungal disease that severely affects cotton fiber yield and quality. Herein, a cotton Trihelix family gene, GhGT-3b_A04, was strongly induced by the fungal pathogen Verticillium dahliae. Overexpression of the gene in Arabidopsis thaliana enhanced the plant's resistance to Verticillium wilt but inhibited the growth of rosette leaves. In addition, the primary root length, root hair number, and root hair length increased in GhGT-3b_A04-overexpressing plants. The density and length of trichomes on the rosette leaves also increased. GhGT-3b_A04 localized to the nucleus, and transcriptome analysis revealed that it induced gene expression for salicylic acid synthesis and signal transduction and activated gene expression for disease resistance. The gene expression for auxin signal transduction and trichome development was reduced in GhGT-3b_A04-overexpressing plants. Our results highlight important regulatory genes for Verticillium wilt resistance and cotton fiber quality improvement. The identification of GhGT-3b_A04 and other important regulatory genes can provide crucial reference information for future research on transgenic cotton breeding.

Cell-specific clock-controlled gene expression program regulates rhythmic fiber cell growth in cotton

BACKGROUND

The epidermis of cotton ovule produces fibers, the most important natural cellulose source for the global textile industry. However, the molecular mechanism of fiber cell growth is still poorly understood.

RESULTS

Here, we develop an optimized protoplasting method, and integrate single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) to systematically characterize the cells of the outer integument of ovules from wild type and fuzzless/lintless (fl) cotton (Gossypium hirsutum). By jointly analyzing the scRNA-seq data from wildtype and fl, we identify five cell populations including the fiber cell type and construct the development trajectory for fiber lineage cells. Interestingly, by time-course diurnal transcriptomic analysis, we demonstrate that the primary growth of fiber cells is a highly regulated circadian rhythmic process. Moreover, we identify a small peptide GhRALF1 that circadian rhythmically controls fiber growth possibly through oscillating auxin signaling and proton pump activity in the plasma membrane. Combining with scATAC-seq, we further identify two cardinal cis-regulatory elements (CREs, TCP motif, and TCP-like motif) which are bound by the trans factors GhTCP14s to modulate the circadian rhythmic metabolism of mitochondria and protein translation through regulating approximately one third of genes that are highly expressed in fiber cells.

CONCLUSIONS

We uncover a fiber-specific circadian clock-controlled gene expression program in regulating fiber growth. This study unprecedentedly reveals a new route to improve fiber traits by engineering the circadian clock of fiber cells.

A glutathione S-transferase GhTT19 determines flower petal pigmentation via regulating anthocyanin accumulation in cotton

Anthocyanin accumulations in the flowers can improve seed production of hybrid lines, and produce higher commodity value in cotton fibre. However, the genetic mechanism underlying the anthocyanin pigmentation in cotton petals is poorly understood. Here, we showed that the red petal phenotype was introgressed from Gossypium bickii through recombination with the segment containing the R₃ ^bic region in the A07 chromosome of Gossypium hirsutum variety LR compared with the near-isogenic line of LW with white flower petals. The cyanidin-3-O-glucoside (Cy3G) was the major anthocyanin in red petals of cotton. A GhTT19 encoding a TT19-like GST was mapped to the R₃ ^bic site associated with red petals via map-based cloning, but GhTT19 homologue gene from the D genome was not expressed in G. hirsutum. Intriguingly, allelic variations in the promoters between GhTT19^LW and GhTT19^LR , rather than genic regions, were found as genetic causal of petal colour variations. GhTT19-GFP was found localized in both the endoplasmic reticulum and tonoplast for facilitating anthocyanin transport. An additional MYB binding element found only in the promoter of GhTT19^LR , but not in that of GhTT19^LW , enhanced its transactivation by the MYB activator GhPAP1. The transgenic analysis confirmed the function of GhTT19 in regulating the red flower phenotype in cotton. The essential light signalling component GhHY5 bonded to and activated the promoter of GhPAP1, and the GhHY5-GhPAP1 module together regulated GhTT19 expression to mediate the light-activation of petal anthocyanin pigmentation in cotton. This study provides new insights into the molecular mechanisms for anthocyanin accumulation and may lay a foundation for faster genetic improvement of cotton.

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.

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

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

Duplicate mutations of GhCYP450 lead to the production of ms5m6 male sterile line in cotton

The duplicated male sterile genes ms₅m₆ in cotton were map-based cloned and validated by the virus-induced gene silencing assays. Duplicate mutations of the GhCYP450 gene encoding a cytochrome P450 protein are responsible for the male sterility in cotton. The utilization of male sterility in cotton plays a vital role in improving yield and fiber quality. A complete male sterile line (ms₅ms₆) has been extensively used to develop hybrid cotton worldwide. Using Zhongkang-A (ZK-A) developed by transferring Bt and ms₅ms₆ genes into the commercial cultivar Zhongmiansuo 12, the duplicate genes were map-based cloned and confirmed via the virus-induced gene silencing (VIGS) assays. The duplicate mutations of GhCYP450 genes encoding a cytochrome P450 protein were responsible for producing male sterility in ms₅ms₆ in cotton. Sequence alignment showed that GhCYP450-Dt in ZK-A differed in two critical aspects from the fertile wild-type TM-1: GhCYP450-Dt has three amino acid (D98E, E168K, G198R) changes in the coding region and a 7-bp (GGAAAAA) insertion in the promoter domain; GhCYP450-At appears to be premature termination of GhCYP450 translation. Further morphological observation and cytological examination of GhCYP450-silenced plants induced by VIGS exhibited shorter filaments and no mature pollen grains. These results indicate that GhCYP450 is essential for pollen exine formation and pollen development for male fertility. Investigating the mechanisms of ms₅ms₆ male sterility will deepen our understanding of the development and utilization of heterosis.

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.

Carpel-specific down-regulation of GhCKXs in cotton significantly enhances seed and fiber yield

Cytokinin is considered to be an important driver of seed yield. To increase the yield of cotton while avoiding the negative consequences caused by constitutive overproduction of cytokinin, we down-regulated specifically the carpel genes for cytokinin oxidase/dehydrogenase (CKX), a key negative regulator of cytokinin levels, in transgenic cotton. The carpel-specific down-regulation of CKXs significantly enhanced cytokinin levels in the carpels. The elevated cytokinin promoted the expression of carpel- and ovule-development-associated genes, GhSTK2, GhAG1, and GhSHP, boosting ovule formation and thus producing more seeds in the ovary. Field experiments showed that the carpel-specific increase of cytokinin significantly increased both seed yield and fiber yield of cotton, without resulting in detrimental phenotypes. Our study details the regulatory mechanism of cytokinin signaling for seed development, and provides an effective and feasible strategy for yield improvement of seed crops.

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.

A deletion/duplication in the Ligon lintless-2 locus induces siRNAs that inhibit cotton fiber cell elongation

Most cultivated cotton (Gossypium hirsutum L.) varieties have two types of seed fibers: short fuzz fiber strongly adhered to the seed coat, and long lint fiber used in the textile industry. The Ligon lintless-2 (Li2) cotton mutant has a normal vegetative phenotype but produces very short lint fiber on the seeds. The Li2 mutation is controlled by a single dominant gene. We discovered a large structural rearrangement at the end of chromosome D13 in the Li2 mutant based on whole-genome sequencing and genetic mapping of segregating populations. The rearrangement contains a 177-kb deletion and a 221-kb duplication positioned as a tandem inverted repeat. The gene Gh_D13G2437 is located at the junction of the inverted repeat in the duplicated region. During transcription such structure spontaneously forms self-complementary hairpin RNA of Gh_D13G2437 followed by production of small interfering RNA (siRNA). Gh_D13G2437 encodes a Ran-Binding Protein 1 (RanBP1) that preferentially expresses during cotton fiber elongation. The abundance of siRNA produced from Gh_D13G2437 reciprocally corresponds with the abundance of highly homologous (68%-98% amino acid sequence identity) RanBP1 family transcripts during fiber elongation, resulting in a shorter fiber phenotype in the Li2. Overexpression of Gh_D13G2437 in the Li2 mutant recovered the long lint fiber phenotype. Taken together, our findings revealed that siRNA-induced silencing of a family of RanBP1s inhibit elongation of cotton fiber cells in the Li2 mutant.

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

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

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

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

Natural variation of GhSI7 increases seed index in cotton

qSI07.1, a major QTL for seed index in cotton, was fine-mapped to a 17.45-kb region, and the candidate gene GhSI7 was verified in transgenic plants. Improving production to meet human needs is a vital objective in cotton breeding. The yield-related trait seed index is a complex quantitative trait, but few candidate genes for seed index have been characterized. Here, a major QTL for seed index qSI07.1 was fine-mapped to a 17.45-kb region by linkage analysis and substitutional mapping. Only GhSI7, encoding the transcriptional regulator STERILE APETALA, was contained in the candidate region. Association test and genetic analysis indicated that an 845-bp-deletion in its intron was responsible for the seed index variation. Origin analysis revealed that this variation was unique in Gossypium hirsutum and originated from race accessions. Overexpression of GhSI7 (haplotype 2) significantly increased the seed index and organ size in cotton plants. Our findings provided a diagnostic marker for breeding and selecting cotton varieties with high seed index, and laid a foundation for further studies to understand the molecular mechanism of cotton seed morphogenesis.

Visible gland constantly traces virus-induced gene silencing in cotton

A virus-induced gene silencing (VIGS) system was established to induce endogenous target gene silencing by post-transcriptional gene silencing (PTGS), which is a powerful tool for gene function analysis in plants. Compared with stable transgenic plant via Agrobacterium-mediated gene transformation, phenotypes after gene knockdown can be obtained rapidly, effectively, and high-throughput through VIGS system. This approach has been successfully applied to explore unknown gene functions involved in plant growth and development, physiological metabolism, and biotic and abiotic stresses in various plants. In this system, GhCLA1 was used as a general control, however, silencing of this gene leads to leaf albino, wilting, and plant death ultimately. As such, it cannot indicate the efficiency of target gene silencing throughout the whole plant growth period. To address this question, in this study, we developed a novel marker gene, Gossypium PIGMENT GLAND FORMATION GENE (GoPGF), as the control to trace the efficiency of gene silencing in the infected tissues. GoPGF has been proved a key gene in gland forming. Suppression of GoPGF does not affect the normal growth and development of cotton. The number of gland altered related to the expression level of GoPGF gene. So it is a good marker that be used to trace the whole growth stages of plant. Moreover, we further developed a method of friction inoculation to enhance and extend the efficiency of VIGS, which facilitates the analysis of gene function in both the vegetative stage and reproductive stage. This improved VIGS technology will be a powerful tool for the rapid functional identification of unknown genes in genomes.

Cotton genes GhMML1 and GhMML2 control trichome branching when ectopically expressed in tobacco

Trichomes exhibit extraordinary diversity in shape, ultrastructure, distribution, secretion capability, biological functions, and morphological differences, which are strongly associated with their multifunction. Previous researches showed MIXTA-like transcription factors involved in regulating trichome initiation and patterning via forming MYB-bHLH-WD40 transcriptional activator complex to induce the expression of downstream genes. Here, we report the characteristics and role of GhMML1 and GhMML2, members of subgroup 9 of the R2R3-type MYB TFs. GhMML1 and GhMML2 were preferentially targeted to the nucleus and prominently expressed in the early stage during fiber development. Ectopic expression of GhMML1 and GhMML2 respectively in the transgenic tobacco plants changed the morphological characteristics of leaf trichomes; that is, the unbranched trichomes turned into multiple branched, and in the meantime, the density of trichomes was reduced on the surface of the leaf. Y2H and LCI assay revealed that both GhMML1 and GhMML2 could physically interact with a bZIP transcription factor family protein (GhbZIP) in vivo and in vitro. It has been reported that GhbZIP's homolog TAG3 in Arabidopsis is involved in the asymmetric growth of leaves and flowers via direct interaction with BOP1. Taken together, our results demonstrated that two MYB MIXTA-like proteins, GhMML1 and GhMML2, together with GhbZIP might form a multimeric complex to involve in trichome development. This study highlights the importance of MIXTA-like genes from TF subgroup 9 and will help to uncover the molecular mechanism underlying differential trichomes and their development.

Broad-spectrum fungal resistance in sorghum is conferred through the complex regulation of an immune receptor gene embedded in a natural antisense transcript

Sorghum (Sorghum bicolor), the fifth most widely grown cereal crop globally, provides food security for millions of people. Anthracnose caused by the fungus Colletotrichum sublineola is a major disease of sorghum worldwide. We discovered a major fungal resistance locus in sorghum composed of the nucleotide-binding leucine-rich repeat receptor gene ANTHRACNOSE RESISTANCE GENE1 (ARG1) that is completely nested in an intron of a cis-natural antisense transcript (NAT) gene designated CARRIER OF ARG1 (CARG). Susceptible genotypes express CARG and two alternatively spliced ARG1 transcripts encoding truncated proteins lacking the leucine-rich repeat domains. In resistant genotypes, elevated expression of an intact allele of ARG1, attributed to the loss of CARG transcription and the presence of miniature inverted-repeat transposable element sequences, resulted in broad-spectrum resistance to fungal pathogens with distinct virulence strategies. Increased ARG1 expression in resistant genotypes is also associated with higher histone H3K4 and H3K36 methylation. In susceptible genotypes, lower ARG1 expression is associated with reduced H3K4 and H3K36 methylation and increased expression of NATs of CARG. The repressive chromatin state associated with H3K9me2 is low in CARG-expressing genotypes within the CARG exon and higher in genotypes with low CARG expression. Thus, ARG1 is regulated by multiple mechanisms and confers broad-spectrum, strong resistance to fungal pathogens.

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.

The GhMYB36 transcription factor confers resistance to biotic and abiotic stress by enhancing PR1 gene expression in plants

Drought and Verticillium wilt disease are two main factors that limit cotton production, which necessitates the identification of key molecular switch to simultaneously improve cotton resistance to Verticillium dahliae and tolerance to drought stress. R2R3-type MYB proteins could play such a role because of their conserved functions in plant development, growth, and metabolism regulation, however, till date a MYB gene conferring the desired resistance to both biotic and abiotic stresses has not been found in cotton. Here, we describe the identification of GhMYB36, a gene encoding a R2R3-type MYB protein in Gossypium hirsutum, which confers drought tolerance and Verticilium wilt resistance in both Arabidopsis and cotton. GhMYB36 was highly induced by PEG-simulated drought stress in G. hirsutum. GhMYB36-silenced cotton plants were more sensitive to both drought stress and Verticillium wilt. GhMYB36 overexpression in transgenic Arabidopsis and cotton plants gave rise to improved drought tolerance and Verticillium wilt resistance. Transient expression of fused GhMYB36-GFP in tobacco cells was able to localize GhMYB36 in the cell nucleus. In addition, RNA-seq analysis together with qRT-PCR validation in transgenic Arabidopsis overexpressing GhMYB36 revealed significantly enhanced PR1 expression. Luciferase interaction assays indicated that GhMYB36 are probably bound to the promoter of PR1 to activate its expression and the interaction, which was further verified by Yeast one hybrid assay. Taken together, our results suggest that GhMYB36 functions as a transcription factor that is involved in drought tolerance and Verticillium wilt resistance in Arabidopsis and cotton by enhancing PR1 expression.

Genome-wide identification and expression analysis of GL2-interacting-repressor (GIR) genes during cotton fiber and fuzz development

GL2-interacting-repressor (GIR) family members may contribute to fiber/fuzz formation via a newly discovered unique pathway in Gossypium arboreum. There are similarities between cotton fiber development and the formation of trichomes and root hairs. The GL2-interacting-repressors (GIRs) are crucial regulators of root hair and trichome formation. The GaFzl gene, annotated as GaGIR1, is negatively associated with trichome development and fuzz initiation. However, there is relatively little available information regarding the other GIR genes in cotton, especially regarding their effects on cotton fiber development. In this study, 21 GIR family genes were identified in the diploid cotton species Gossypium arboreum; these genes were divided into three groups. The GIR genes were characterized in terms of their phylogenetic relationships, structures, chromosomal distribution and evolutionary dynamics. These GIR genes were revealed to be unequally distributed on 12 chromosomes in the diploid cotton genome, with no GIR gene detected on Ga06. The cis-acting elements in the promoter regions were predicted to be responsive to light, phytohormones, defense activities and stress. The transcriptomic data and qRT-PCR results revealed that most GIR genes were not differentially expressed between the wild-type control and the fuzzless mutant line. Moreover, 14 of 21 family genes were expressed at high levels, indicating these genes may play important roles during fiber development and fuzz formation. Furthermore, Ga01G0231 was predominantly expressed in root samples, suggestive of a role in root hair formation rather than in fuzz initiation and development. The results of this study have enhanced our understanding of the GIR genes and their potential utility for improving cotton fiber through breeding.

Identification, evolutionary analysis and functional diversification of RAV gene family in cotton (G. hirsutum L.)

Genome wide analysis, expression pattern analysis, and functional characterization of RAV genes highlight their roles in roots, stem development and hormonal response. RAV (Related to ABI3 and VP1) gene family members have been involved in tissues/organs growth and hormone signaling in various plant species. Here, we identified 247 RAVs from 12 different species with 33 RAV genes from G. hirsutum. Phylogenetic analysis classified RAV genes into four distinct groups. Analysis of gene structure showed that most GhRAVs lack introns. Motif distribution pattern and protein sequence logos indicated that GhRAV genes were highly conserved during the process of evolution. Promotor cis-acting elements revealed that promotor regions of GhRAV genes encode numerous elements related to plant growth, abiotic stresses and phytohormones. Chromosomal location information showed uneven distribution of 33 GhRAV genes on different chromosomes. Collinearity analysis identified 628 and 52 orthologous/ paralogous gene pairs in G. hirsutum and G. barbadense, respectively. Ka/Ks values indicated that GhRAV and GbRAV genes underwent strong purifying selection pressure. Selecton model and codon model selection revealed that GhRAV amino acids were under purifying selection and adaptive evolution exists among GhRAV proteins. Three dimensional structure of GhRAVs indicated the presence of numerous alpha helix and beta-barrels. Expression level revealed that some GhRAV genes exhibited high expression in roots (GhRAV3, GhRAV4, GhRAV11, GhRAV18, GhRAV20 and GhRAV30) and stem (GhRAV3 and GhRAV18), indicating their potential role in roots and stem development. GhRAV genes can be regulated by phytohormonal stresses (BL, JA and IAA). Our study provides a reference for future studies related to the functional analysis of GhRAVs in cotton.

GhMYB7 promotes secondary wall cellulose deposition in cotton fibres by regulating GhCesA gene expression through three distinct cis-elements

Cotton fibre is the most important source for natural textiles. The secondary cell walls (SCWs) of mature cotton fibres contain the highest proportion of cellulose content (> 90%) in any plant. The onset and progression of SCW cellulose synthesis need to be tightly controlled to balance fibre elongation and cell wall deposition. However, regulatory mechanisms that control cellulose synthesis during cotton fibre growth remain elusive. Here, we conducted genetic and functional analyses demonstrating that the R2R3-MYB GhMYB7 controls cotton fibre cellulose synthesis. Overexpression of GhMYB7 in cotton sped up SCW cellulose biosynthesis in fibre cells, and led to shorter fibres with thicker walls. By contrast, RNA interference (RNAi) silencing of GhMYB7 delayed fibre SCW cellulose synthesis and resulted in elongated fibres with thinner walls. Furthermore, we demonstrated that GhMYB7 regulated cotton fibre SCW cellulose synthases by directly binding to three distinct cis-elements in the respective GhCesA4, GhCesA7 and GhCesA8 promoters. We found that this regulatory mechanism of cellulose synthesis was 'hi-jacked' also by other GhMYBs. Together, our findings uncover a hitherto-unknown mechanism that cotton fibre employs to regulate SCW cellulose synthesis. Our results also provide a strategy for genetic improvement of SCW thickness of cotton fibre.

Making sense of the natural antisense transcript puzzle

In plants, thousands of genes are associated with antisense transcription, which often produces noncoding RNAs. Although widespread, sense-antisense pairs have been implicated in a limited variety of functions in plants and are often thought to form extensive dsRNA stretches triggering gene silencing. In this opinion, we show that evidence does not support gene silencing as a major role for antisense transcription. In fact, it is more likely that antisense transcripts play diverse functions in gene regulation. We propose a general framework for the initial functional dissection of antisense transcripts, suggesting testable hypotheses relying on an experiment-based decision tree. By moving beyond the gene silencing paradigm, we argue that a broad and diverse role for natural antisense transcription will emerge.

Function analysis of a cotton R2R3 MYB transcription factor GhMYB3 in regulating plant trichome development

Cotton is an important fibre-producing crop. Cotton fibres consist of highly elongated trichomes derived from the ovule. To improve the quality of cotton, it is necessary to identify the genes regulating fibre development. GhMYB3 was identified through bioinfomatic analysis and introduced to Arabidopsis and cotton to observe the phenotype. Protein inteaction and promoter bingding assays were conducted to explore the role of GhMYB3 in trichome fibre growth. Cotton fibre development might share a similar regulatory mechanism to Arabidopsis leaf trichomes, which is determined by the essential regulatory complex, MYB-bHLH-WD40. The GL1-like R2R3 MYB transcription factor GhMYB3 interacts with the AtGL3 protein involved in Arabidopsis trichome development. Ectopic expression of GhMYB3 could rescue the glabrous phenotype of the Arabidopsis gl1 mutant and produced more ectopic trichomes on inflorescence stems and floral organs, confirming its orthologous function in plant trichome development. The expression of GhMYB3 increased in response to exogenous gibberellin (GA₃ ), auxin (IAA) and methyl jasmonate (MeJA). Overexpression of this gene in cotton leads to a slight increase in fibre length and lint percentage, possibly by activating the transcription of its downstream gene GhRDL1 or other fibre-related genes. The results increase our understanding of the key role of GhMYB3 in positively controlling plant trichome development, and this gene could be a potential target for molecular breeding in cotton.

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.

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.

Genome-Wide Introgression and Quantitative Trait Locus Mapping Reveals the Potential of Asian Cotton (Gossypium arboreum) in Improving Upland Cotton (Gossypium hirsutum)

Gossypium arboreum (2n=2x=26, A₂), the putative progenitor of the A_t-subgenome of Gossypium hirsutum (2n=4x=52, AD), is a repository of genes of interesting that have been eliminated during evolution/domestication of G. hirsutum. However, its valuable genes remain untapped so far due to species isolation. Here, using a synthetic amphiploid (AADDA₂A₂) previously reported, we developed a set of 289 G. arboreum chromosome segment introgression lines (ILs) in G. hirsutum by expanding the backcrossing population and through precise marker-assisted selection (MAS) although complex chromosomal structural variations existed between parents which severely hindered introgression. Our results showed the total coverage length of introgressed segments was 1,116.29 Mb, representing 78.48% of the A_t-subgenome in the G. hirsutum background, with an average segment-length of 8.69 Mb. A total of 81 co- quantitative trait loci (QTLs) for yield and fiber quality were identified by both the RSTEP-ADD-based QTL mapping and the genome-wide association study (GWAS) analysis, with 1.01-24.78% of the phenotypic variance explained. Most QTLs for boll traits showed negative additive effects, but G. arboreum still has the potential to improve boll-number traits in G. hirsutum. Most QTLs for fiber quality showed negative additive effects, implying these QTLs were domesticated in G. hirsutum compared with G. arboreum and, a small quantity of fiber quality QTLs showing positive additive effects, conversely; however, indicates that G. arboreum has the underlying genes of enhancing fiber quality of G. hirsutum. This study provides new insights into the breeding genetic potential of G. arboreum, lays the foundation for further mining favorable genes of interest, and provides guidance for inter-ploidy gene transference from relatives into cultivated crops.

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.

Mapping-by-sequencing the locus of EMS-induced mutation responsible for tufted-fuzzless seed phenotype in cotton

Cotton fiber mutants are valuable resources for studying functions of altered genes and their roles in fiber development. The n₄^t is a recessive tufted-fuzzless seed mutant created through chemical mutagenesis with ethyl methanesulfonate. Genetic analysis indicated that the tufted-fuzzless phenotype is controlled by a single recessive locus. In this study, we developed an F₂ population of 602 progeny plants and sequenced the genomes of the parents and two DNA bulks from F₂ progenies showing the mutant phenotype. We identified DNA sequence variants between the tufted-fuzzless mutant and wild type by aligning the sequence reads to the reference TM-1 genome and designed subgenome-specific SNP markers. We mapped the n₄^t locus on chromosome D04 within a genomic interval of about 411 kb. In this region, seven genes showed significant differential expression between the tufted-fuzzless mutant and wild type. Possible candidate genes are discussed in this study. The utilization of the n₄^t mutant along with other fiber mutants will facilitate our understanding of the molecular mechanisms of cotton fiber cell growth and development.

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.

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.

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.

Genetic mapping and transcriptomic characterization of a new fuzzless-tufted cottonseed mutant

Fiber mutants are unique and valuable resources for understanding the genetic and molecular mechanisms controlling initiation and development of cotton fibers that are extremely elongated single epidermal cells protruding from the seed coat of cottonseeds. In this study, we reported a new fuzzless-tufted cotton mutant (Gossypium hirsutum) and showed that fuzzless-tufted near-isogenic lines (NILs) had similar agronomic traits and a higher ginning efficiency compared to their recurrent parents with normal fuzzy seeds. Genetic analysis revealed that the mutant phenotype is determined by a single incomplete dominant locus, designated N5. The mutation was fine mapped to an approximately 250-kb interval containing 33 annotated genes using a combination of bulked segregant sequencing, SNP chip genotyping, and fine mapping. Comparative transcriptomic analysis using 0-6 days post-anthesis (dpa) ovules from NILs segregating for the phenotypes of fuzzless-tufted (mutant) and normal fuzzy cottonseeds (wild-type) uncovered candidate genes responsible for the mutant phenotype. It also revealed that the flanking region of the N5 locus is enriched with differentially expressed genes (DEGs) between the mutant and wild-type. Several of those DEGs are members of the gene families with demonstrated roles in cell initiation and elongation, such as calcium-dependent protein kinase and expansin. The transcriptome landscape of the mutant was significantly reprogrammed in the 6 dpa ovules and, to a less extent, in the 0 dpa ovules, but not in the 2 and 4 dpa ovules. At both 0 and 6 dpa, the reprogrammed mutant transcriptome was mainly associated with cell wall modifications and transmembrane transportation, while transcription factor activity was significantly altered in the 6 dpa mutant ovules. These results imply a similar molecular basis for initiation of lint and fuzz fibers despite certain differences.

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.

MIXTAs and phytohormones orchestrate cotton fiber development

Cotton is the largest source of natural fiber for textile industry in the world. Cotton fibers are seed trichomes that make cotton unique among plants. Cotton fibers originate from ovule epidermal cells and serve as an excellent model to study the process of cell differentiation in plants. Characterization of factors contributing to fiber development will help to reveal general mechanisms of cell differentiation in plants. Transcription factors (TFs), especially MYB-MIXTA-like (MML) factors, appear to have evolved unique roles in fiber development. In addition, phytohormones including brassinosteroids, jasmonic acid, GA and auxin also play an important role in regulating fiber development. Here, we summarize the mechanisms of MIXTAs and phytohormones orchestrating cotton fiber development. The progress in understanding molecular basis of fiber development will facilitate future genetic engineering and breeding to improve cotton fiber quality and yield.

Response of phytohormone mediated plant homeodomain (PHD) family to abiotic stress in upland cotton (Gossypium hirsutum spp.)

BACKGROUND

The sequencing and annotations of cotton genomes provide powerful theoretical support to unravel more physiological and functional information. Plant homeodomain (PHD) protein family has been reported to be involved in regulating various biological processes in plants. However, their functional studies have not yet been carried out in cotton.

RESULTS

In this study, 108, 55, and 52 PHD genes were identified in G. hirsutum, G. raimondii, and G. arboreum, respectively. A total of 297 PHD genes from three cotton species, Arabidopsis, and rice were divided into five groups. We performed chromosomal location, phylogenetic relationship, gene structure, and conserved domain analysis for GhPHD genes. GhPHD genes were unevenly distributed on each chromosome. However, more GhPHD genes were distributed on At_05, Dt_05, and At_07 chromosomes. GhPHD proteins depicted conserved domains, and GhPHD genes exhibiting similar gene structure were clustered together. Further, whole genome duplication (WGD) analysis indicated that purification selection greatly contributed to the functional maintenance of GhPHD gene family. Expression pattern analysis based on RNA-seq data showed that most GhPHD genes showed clear tissue-specific spatiotemporal expression patterns elucidating the multiple functions of GhPHDs in plant growth and development. Moreover, analysis of cis-acting elements revealed that GhPHDs may respond to a variety of abiotic and phytohormonal stresses. In this regard, some GhPHD genes showed good response against abiotic and phytohormonal stresses. Additionally, co-expression network analysis indicated that GhPHDs are essential for plant growth and development, while GhPHD genes response against abiotic and phytohormonal stresses may help to improve plant tolerance in adverse environmental conditions.

CONCLUSION

This study will provide useful information to facilitate further research related to the vital roles of GhPHD gene family in plant growth and development.