Genome-wide identification of the GhARF gene family reveals that GhARF2 and GhARF18 are involved in cotton fibre cell initiation

Strigolactone-induced degradation of SMXL7 and SMXL8 contributes to gibberellin- and auxin-mediated fiber cell elongation in cotton

Cotton (Gossypium) fiber length, a key trait determining fiber yield and quality, is highly regulated by a class of recently identified phytohormones, strigolactones (SLs). However, the underlying molecular mechanisms of SL signaling involved in fiber cell development are largely unknown. Here, we show that the SL signaling repressors MORE AXILLARY GROWTH2-LIKE7 (GhSMXL7) and GhSMXL8 negatively regulate cotton fiber elongation. Specifically, GhSMXL7 and GhSMXL8 inhibit the polyubiquitination and degradation of the gibberellin (GA)-triggered DELLA protein (GhSLR1). Biochemical analysis revealed that GhSMXL7 and GhSMXL8 physically interact with GhSLR1, which interferes with the association of GhSLR1 with the E3 ligase GA INSENSITIVE2 (GhGID2), leading to the repression of GA signal transduction. GhSMXL7 also interacts with the transcription factor GhHOX3, preventing its binding to the promoters of essential fiber elongation regulatory genes. Moreover, both GhSMXL7 and GhSMXL8 directly bind to the promoter regions of the AUXIN RESPONSE FACTOR (ARF) genes GhARF18-10A, GhARF18-10D, and GhARF19-7D to suppress their expression. Cotton plants in which GhARF18-10A, GhARF18-10D, and GhARF19-7D transcript levels had been reduced by virus-induced gene silencing (VIGS) displayed reduced fiber length compared with control plants. Collectively, our findings reveal a mechanism illustrating how SL integrates GA and auxin signaling to coordinately regulate plant cell elongation at the single-cell level.

Exploring Evolutionary Pathways and Abiotic Stress Responses through Genome-Wide Identification and Analysis of the Alternative Oxidase (AOX) Gene Family in Common Oat (Avena sativa)

The alternative oxidase (AOX), a common terminal oxidase in the electron transfer chain (ETC) of plants, plays a crucial role in stress resilience and plant growth and development. Oat (Avena sativa), an important crop with high nutritional value, has not been comprehensively studied regarding the AsAOX gene family. Therefore, this study explored the responses and potential functions of the AsAOX gene family to various abiotic stresses and their potential evolutionary pathways. Additionally, we conducted a genome-wide analysis to explore the evolutionary conservation and divergence of AOX gene families among three Avena species (Avena sativa, Avena insularis, Avena longiglumis) and four Poaceae species (Avena sativa, Oryza sativa, Triticum aestivum, and Brachypodium distachyon). We identified 12 AsAOX, 9 AiAOX, and 4 AlAOX gene family members. Phylogenetic, motif, domain, gene structure, and selective pressure analyses revealed that most AsAOXs, AiAOXs, and AlAOXs are evolutionarily conserved. We also identified 16 AsAOX segmental duplication pairs, suggesting that segmental duplication may have contributed to the expansion of the AsAOX gene family, potentially preserving these genes through subfunctionalization. Chromosome polyploidization, gene structural variations, and gene fragment recombination likely contributed to the evolution and expansion of the AsAOX gene family as well. Additionally, we hypothesize that AsAOX2 may have potential function in resisting wounding and heat stresses, while AsAOX4 could be specifically involved in mitigating wounding stress. AsAOX11 might contribute to resistance against chromium and waterlogging stresses. AsAOX8 may have potential fuction in mitigating ABA-mediated stress. AsAOX12 and AsAOX5 are most likely to have potential function in mitigating salt and drought stresses, respectively. This study elucidates the potential evolutionary pathways of the AsAOXs gene family, explores their responses and potential functions to various abiotic stresses, identifies potential candidate genes for future functional studies, and facilitates molecular breeding applications in A. sativa.

Genetic Regulatory Perturbation of Gene Expression Impacted by Genomic Introgression in Fiber Development of Allotetraploid Cotton

Interspecific genomic introgression is an important evolutionary process with respect to the generation of novel phenotypic diversity and adaptation. A key question is how gene flow perturbs gene expression networks and regulatory interactions. Here, an introgression population of two species of allopolyploid cotton (Gossypium) to delineate the regulatory perturbations of gene expression regarding fiber development accompanying fiber quality change is utilized. De novo assembly of the recipient parent (G. hirsutum Emian22) genome allowed the identification of genomic variation and introgression segments (ISs) in 323 introgression lines (ILs) from the donor parent (G. barbadense 3-79). It documented gene expression dynamics by sequencing 1,284 transcriptomes of developing fibers and characterized genetic regulatory perturbations mediated by genomic introgression using a multi-locus model. Introgression of individual homoeologous genes exhibiting extreme low or high expression bias can lead to a parallel expression bias in their non-introgressed duplicates, implying a shared yet divergent regulatory fate of duplicated genes following allopolyploidy. Additionally, the IL N182 with improved fiber quality is characterized, and the candidate gene GhFLAP1 related to fiber length is validated. This study outlines a framework for understanding introgression-mediated regulatory perturbations in polyploids, and provides insights for targeted breeding of superior upland cotton fiber.

GhMAX2 Contributes to Auxin-Mediated Fiber Elongation in Cotton (Gossypium hirsutum)

Strigolactones (SLs) represent a new group of phytohormones that play a pivotal role in the regulation of plant shoot branching and the development of adventitious roots. In cotton (Gossypium hirsutum, Gh), SLs play a crucial role in the regulation of fiber cell elongation and secondary cell wall thickness. However, the underlying molecular mechanisms of SL signaling involved in fiber cell development are largely unknown. In this study, we report two SL-signaling genes, GhMAX2-3 and GhMAX2-6, which positively regulate cotton fiber elongation. Further protein-protein interaction and degradation assays showed that the repressor of the auxin cascade GhIAA17 serves as a substrate for the F-box E3 ligase GhMAX2. The in vivo ubiquitination assay suggested that GhMAX2-3 and GhMAX2-6 ubiquitinate GhIAA17 and coordinately degrade GhIAA17 with GhTIR1. The findings of this investigation offer valuable insights into the roles of GhMAX2-mediated SL signaling in cotton and establish a solid foundation for future endeavors aimed at optimizing cotton plant cultivation.

Auxin resistant 2 and short hypocotyl 2 regulate cotton fiber initiation and elongation

Auxin, a pivotal regulator of diverse plant growth processes, remains central to development. The auxin-responsive genes auxin/indole-3-acetic acids (AUX/IAAs) are indispensable for auxin signal transduction, which is achieved through intricate interactions with auxin response factors (ARFs). Despite this, the potential of AUX/IAAs to govern the development of the most fundamental biological unit, the single cell, remains unclear. In this study, we harnessed cotton (Gossypium hirsutum) fiber, a classic model for plant single-cell investigation, to determine the complexities of AUX/IAAs. Our research identified 2 pivotal AUX/IAAs, auxin resistant 2 (GhAXR2) and short hypocotyl 2 (GhSHY2), which exhibit opposite control over fiber development. Notably, suppressing GhAXR2 reduced fiber elongation, while silencing GhSHY2 fostered enhanced fiber elongation. Investigating the mechanistic intricacies, we identified specific interactions between GhAXR2 and GhSHY2 with distinct ARFs. GhAXR2's interaction with GhARF6-1 and GhARF23-2 promoted fiber cell development through direct binding to the AuxRE cis-element in the constitutive triple response 1 promoter, resulting in transcriptional inhibition. In contrast, the interaction of GhSHY2 with GhARF7-1 and GhARF19-1 exerted a negative regulatory effect, inhibiting fiber cell growth by activating the transcription of xyloglucan endotransglucosylase/hydrolase 9 and cinnamate-4-hydroxylase. Thus, our study reveals the intricate regulatory networks surrounding GhAXR2 and GhSHY2, elucidating the complex interplay of multiple ARFs in AUX/IAA-mediated fiber cell growth. This work enhances our understanding of single-cell development and has potential implications for advancing plant growth strategies and agricultural enhancements.

GhWER controls fiber initiation and early elongation by regulating ethylene signaling pathway in cotton (Gossypium hirsutum)

Cotton fibers are specialized single-cell trichomes derived from epidermal cells, similar to root hairs and trichomes in Arabidopsis. While the MYB-bHLH-WD40 (MBW) complex has been shown to regulate initiation of both root hairs and trichomes in Arabidopsis, the role of their homologous gene in cotton fiber initiation remains unknown. In this study, we identified a R2R3 MYB transcription factor (TF), GhWER, which exhibited a significant increase in expression within the outer integument of ovule at -1.5 DPA (days post anthesis). Its expression peaked at -1 DPA and then gradually decreased. Knockout of GhWER using CRISPR technology inhibited the initiation and early elongation of fiber initials, resulting in the shorter mature fiber length. Additionally, GhWER interacted with two bHLH TF, GhDEL65 and GhbHLH121, suggesting a potential regulatory complex for fiber development. RNA-seq analysis of the outer integument of the ovule at -1.5 DPA revealed that the signal transduction pathways of ethylene, auxin and gibberellin were affected in the GhWER knockout lines. Further examination demonstrated that GhWER directly activated ethylene signaling genes, including ACS1 and ETR2. These findings highlighted the biological function of GhWER in regulating cotton fiber initiation and early elongation, which has practical significance for improving fiber quality and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01477-6.

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.

The transcription factor ERF108 interacts with AUXIN RESPONSE FACTORs to mediate cotton fiber secondary cell wall biosynthesis

Phytohormones play indispensable roles in plant growth and development. However, the molecular mechanisms underlying phytohormone-mediated regulation of fiber secondary cell wall (SCW) formation in cotton (Gossypium hirsutum) remain largely underexplored. Here, we provide mechanistic evidence for functional interplay between the APETALA2/ethylene response factor (AP2/ERF) transcription factor GhERF108 and auxin response factors GhARF7-1 and GhARF7-2 in dictating the ethylene-auxin signaling crosstalk that regulates fiber SCW biosynthesis. Specifically, in vitro cotton ovule culture revealed that ethylene and auxin promote fiber SCW deposition. GhERF108 RNA interference (RNAi) cotton displayed remarkably reduced cell wall thickness compared with controls. GhERF108 interacted with GhARF7-1 and GhARF7-2 to enhance the activation of the MYB transcription factor gene GhMYBL1 (MYB domain-like protein 1) in fibers. GhARF7-1 and GhARF7-2 respond to auxin signals that promote fiber SCW thickening. GhMYBL1 RNAi and GhARF7-1 and GhARF7-2 virus-induced gene silencing (VIGS) cotton displayed similar defects in fiber SCW formation as GhERF108 RNAi cotton. Moreover, the ethylene and auxin responses were reduced in GhMYBL1 RNAi plants. GhMYBL1 directly binds to the promoters of GhCesA4-1, GhCesA4-2, and GhCesA8-1 and activates their expression to promote cellulose biosynthesis, thereby boosting fiber SCW formation. Collectively, our findings demonstrate that the collaboration between GhERF108 and GhARF7-1 or GhARF7-2 establishes ethylene-auxin signaling crosstalk to activate GhMYBL1, ultimately leading to the activation of fiber SCW biosynthesis.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture

BACKGROUND

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production.

GENOMICS APPROACHES IN COTTON

In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research.

CONCLUSIONS

This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

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

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

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.

Genome-Wide Identification and Expression Analysis of Auxin Response Factor Gene Family in Linum usitatissimum

Auxin response factors (ARFs) are critical components of the auxin signaling pathway, and are involved in diverse plant biological processes. However, ARF genes have not been investigated in flax (Linum usitatissimum L.), an important oilseed and fiber crop. In this study, we comprehensively analyzed the ARF gene family and identified 33 LuARF genes unevenly distributed on the 13 chromosomes of Longya-10, an oil-use flax variety. Detailed analysis revealed wide variation among the ARF family members and predicted nuclear localization for all proteins. Nineteen LuARFs contained a complete ARF structure, including DBD, MR, and CTD, whereas the other fourteen lacked the CTD. Phylogenetic analysis grouped the LuARFs into four (I-V) clades. Combined with sequence analysis, the LuARFs from the same clade showed structural conservation, implying functional redundancy. Duplication analysis identified twenty-seven whole-genome-duplicated LuARF genes and four tandem-duplicated LuARF genes. These duplicated gene pairs' K_a/K_s ratios suggested a strong purifying selection pressure on the LuARF genes. Collinearity analysis revealed that about half of the LuARF genes had homologs in other species, indicating a relatively conserved nature of the ARFs. The promoter analysis identified numerous hormone- and stress-related elements, and the qRT-PCR experiment revealed that all LuARF genes were responsive to phytohormone (IAA, GA3, and NAA) and stress (PEG, NaCl, cold, and heat) treatments. Finally, expression profiling of LuARF genes in different tissues by qRT-PCR indicated their specific functions in stem or capsule growth. Thus, our findings suggest the potential functions of LuARFs in flax growth and response to an exogenous stimulus, providing a basis for further functional studies on these genes.

Genome-wide identification of the geranylgeranyl pyrophosphate synthase (GGPS) gene family involved in chlorophyll synthesis in cotton

BACKGROUND

Geranylgeranyl pyrophosphate synthase (GGPS) is a structural enzyme of the terpene biosynthesis pathway that is involved in regulating plant photosynthesis, growth and development, but this gene family has not been systematically studied in cotton.

RESULTS

In the current research, genome-wide identification was performed, and a total of 75 GGPS family members were found in four cotton species, Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii. The GGPS genes were divided into three subgroups by evolutionary analysis. Subcellular localization prediction showed that they were mainly located in chloroplasts and plastids. The closely related GGPS contains a similar gene structure and conserved motif, but some genes are quite different, resulting in functional differentiation. Chromosome location analysis, collinearity and selection pressure analysis showed that many fragment duplication events occurred in GGPS genes. Three-dimensional structure analysis and conservative sequence analysis showed that the members of the GGPS family contained a large number of α-helices and random crimps, and all contained two aspartic acid-rich domains, DDxxxxD and DDxxD (x is an arbitrary amino acid), suggesting its key role in function. Cis-regulatory element analysis showed that cotton GGPS may be involved in light response, abiotic stress and other processes. A GGPS gene was silenced successfully by virus-induced gene silencing (VIGS), and it was found that the chlorophyll content in cotton leaves decreased significantly, suggesting that the gene plays an important role in plant photosynthesis.

CONCLUSIONS

In total, 75 genes were identified in four Gossypium species by a series of bioinformatics analysis. Gene silencing from GGPS members of G. hirsutum revealed that GGPS plays an important regulatory role in photosynthesis. This study provides a theoretical basis for the biological function of GGPS in cotton growth and development.

LkARF7 and LkARF19 overexpression promote adventitious root formation in a heterologous poplar model by positively regulating LkBBM1

Cuttage propagation involves adventitious root formation induced by auxin. In our previous study, Larix kaempferi BABY BOOM 1 (LkBBM1), which is known to regulate adventitious root formation, was affected by auxin. However, the relationship between LkBBM1 and auxin remains unclear. Auxin response factors (ARFs) are a class of important transcription factors in the auxin signaling pathway and modulate the expression of early auxin-responsive genes by binding to auxin response elements. In the present study, we identified 14 L. kaempferi ARFs (LkARFs), and found LkARF7 and LkARF19 bound to LkBBM1 promoter and enhanced its transcription using yeast one-hybrid, ChIP-qPCR, and dual-luciferase assays. In addition, the treatment with naphthalene acetic acid promoted the expression of LkARF7 and LkARF19. We also found that overexpression of these two genes in poplar promoted adventitious root formation. Furthermore, LkARF19 interacted with the DEAD-box ATP-dependent RNA helicase 53-like protein to form a heterodimer to regulate adventitious root formation. Altogether, our results reveal an additional regulatory mechanism underlying the control of adventitious root formation by auxin.

Poplar coma morphogenesis and miRNA regulatory networks by combining ovary tissue sectioning and deep sequencing

Poplar coma, commonly referred to as "seed hairs", is a tuft of trichomes attached to the seed coat that helps seed dispersal. However, they can also trigger health impacts for humans, including sneezing, shortness of breath, and skin irritation. Despite efforts to study the regulatory mechanism of herbaceous trichome formation, poplar coma remains poorly understood. In this study, we showed that the epidermal cells of the funiculus and placenta are the origin of poplar coma based on observations of paraffin sections. Small RNA (sRNA) and degradome libraries were also constructed at three stages of poplar coma development, including initiation and elongation stages. Based on 7,904 miRNA-target pairs identified by small RNA and degradome sequencing, we constructed a miRNA-transcript factor and a stage-specific miRNA regulatory network. By combining paraffin section observation and deep sequencing, our research will provide greater insight into the molecular mechanisms of poplar coma development.

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.

Genome-Scale Investigation of GARP Family Genes Reveals Their Pivotal Roles in Nutrient Stress Resistance in Allotetraploid Rapeseed

The GARP genes are plant-specific transcription factors (TFs) and play key roles in regulating plant development and abiotic stress resistance. However, few systematic analyses of GARPs have been reported in allotetraploid rapeseed (Brassica napus L.) yet. In the present study, a total of 146 BnaGARP members were identified from the rapeseed genome based on the sequence signature. The BnaGARP TFs were divided into five subfamilies: ARR, GLK, NIGT1/HRS1/HHO, KAN, and PHL subfamilies, and the members within the same subfamilies shared similar exon-intron structures and conserved motif configuration. Analyses of the Ka/Ks ratios indicated that the GARP family principally underwent purifying selection. Several cis-acting regulatory elements, essential for plant growth and diverse biotic and abiotic stresses, were identified in the promoter regions of BnaGARPs. Further, 29 putative miRNAs were identified to be targeting BnaGARPs. Differential expression of BnaGARPs under low nitrate, ammonium toxicity, limited phosphate, deficient boron, salt stress, and cadmium toxicity conditions indicated their potential involvement in diverse nutrient stress responses. Notably, BnaA9.HHO1 and BnaA1.HHO5 were simultaneously transcriptionally responsive to these nutrient stresses in both hoots and roots, which indicated that BnaA9.HHO1 and BnaA1.HHO5 might play a core role in regulating rapeseed resistance to nutrient stresses. Therefore, this study would enrich our understanding of molecular characteristics of the rapeseed GARPs and will provide valuable candidate genes for further in-depth study of the GARP-mediated nutrient stress resistance in rapeseed.

Roles of Auxin in the Growth, Development, and Stress Tolerance of Horticultural Plants

Auxin, a plant hormone, regulates virtually every aspect of plant growth and development. Many current studies on auxin focus on the model plant Arabidopsis thaliana, or on field crops, such as rice and wheat. There are relatively few studies on what role auxin plays in various physiological processes of a range of horticultural plants. In this paper, recent studies on the role of auxin in horticultural plant growth, development, and stress response are reviewed to provide novel insights for horticultural researchers and cultivators to improve the quality and application of horticultural crops.

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.

Genome-wide identification and expression profile of GhGRF gene family in Gossypium hirsutum L

Background

Cotton is the primary source of renewable natural fiber in the textile industry and an important biodiesel crop. Growth regulating factors (GRFs) are involved in regulating plant growth and development.

Methods

Using genome-wide analysis, we identified 35 GRF genes in Gossypium hirsutum.

Results

Chromosomal location information revealed an uneven distribution of GhGRF genes, with maximum genes on chromosomes A02, A05, and A12 from the At sub-genome and their corresponding D05 and D12 from the Dt sub-genome. In the phylogenetic tree, 35 GRF genes were divided into five groups, including G1, G2, G3, G4, and G5. The majority of GhGRF genes have two to three introns and three to four exons, and their deduced proteins contained conserved QLQ and WRC domains in the N-terminal end of GRFs in Arabidopsis and rice. Sequence logos revealed that GRF genes were highly conserved during the long-term evolutionary process. The CDS of the GhGRF gene can complement MiRNA396a. Moreover, most GhGRF genes transcripts developed high levels of ovules and fibers. Analyses of promoter cis-elements and expression patterns indicated that GhGRF genes play an essential role in regulating plant growth and development by coordinating the internal and external environment and multiple hormone signaling pathways. Our analysis indicated that GhGRFs are ideal target genes with significant potential for improving the molecular structure of cotton.

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.

An Overview of Cotton Gland Development and Its Transcriptional Regulation

Cotton refers to species in the genus Gossypium that bear spinnable seed coat fibers. A total of 50 species in the genus Gossypium have been described to date. Of these, only four species, viz. Gossypium, hirsutum, G. barbadense, G. arboretum, and G. herbaceum are cultivated; the rest are wild. The black dot-like structures on the surfaces of cotton organs or tissues, such as the leaves, stem, calyx, bracts, and boll surface, are called gossypol glands or pigment glands, which store terpenoid aldehydes, including gossypol. The cotton (Gossypium hirsutum) pigment gland is a distinctive structure that stores gossypol and its derivatives. It provides an ideal system for studying cell differentiation and organogenesis. However, only a few genes involved in the process of gland formation have been identified to date, and the molecular mechanisms underlying gland initiation remain unclear. The terpenoid aldehydes in the lysigenous glands of Gossypium species are important secondary phytoalexins (with gossypol being the most important) and one of the main defenses of plants against pests and diseases. Here, we review recent research on the development of gossypol glands in Gossypium species, the regulation of the terpenoid aldehyde biosynthesis pathway, discoveries from genetic engineering studies, and future research directions.

Genome-wide analysis and expression profiling of Cation/H+ exchanger (CAX) family genes reveal likely functions in cadmium stress responses in poplar

Cadmium, a toxic heavy metal, seriously affects human health and ecological security. The cation/H⁺ exchanger (CAX) family is a unique metal transporter that plays a crucial role in Cd acquisition, transfer, and remission in plants. Although there are many studies related to the genome-wide analysis of Populus trichocarpa, little research has been done on the CAX family genes, especially concerning Cd stress. In this study, genome-wide analysis of the Populus CAX family identified seven stress-related CAX genes. The evolutionary tree indicated that the CaCA family genes were grouped into four clusters. Moreover, seven pairs of genes were derived by segmental duplication in poplars. Cis-acting element analysis identified numerous stress-related elements in the promoters of diverse PtrCAXs. Furthermore, some PtrCAXs were up-regulated by drought, beetle, and mechanical damage, indicating their possible function in regulating stress response. Under cadmium stress, all CAX genes in the roots were up-regulated. Our findings suggest that plants may regulate their response to Cd stress through the TF-CAXs module. Comprehensively investigating the CAX family provides a scientific basis for the phytoremediation of heavy metal pollution by Populus.

Genome-Wide Investigation of the PtrCHLP Family Reveals That PtrCHLP3 Actively Mediates Poplar Growth and Development by Regulating Photosynthesis

Chlorophyll (Chl) plays a crucial role in plant photosynthesis. The geranylgeraniol reductase gene (CHLP) participates in the terminal hydrogenation of chlorophyll biosynthesis. Although there are many studies related to the genome-wide analysis of Populus trichocarpa, little research has been conducted on CHLP family genes, especially those concerning growth and photosynthesis. In this study, three CHLP genes were identified in Populus. The evolutionary tree indicated that the CHLP family genes were divided into six groups. Moreover, one pair of genes was derived from segmental duplications in Populus. Many elements related to growth were detected by cis-acting element analysis of the promoters of diverse PtrCHLPs. Furthermore, PtrCHLPs exhibit different tissue expression patterns. In addition, PtrCHLP3 is preferentially expressed in the leaves and plays an important role in regulating chlorophyll biosynthesis. Silencing of PtrCHLP3 in poplar resulted in a decrease in chlorophyll synthesis in plants, thus blocking electron transport during photosynthesis. Furthermore, inhibition of PtrCHLP3 expression in poplar can inhibit plant growth through the downregulation of photosynthesis. Ultimately, PtrCHLP3 formed a co-expression network with photosynthesis and chlorophyll biosynthesis-related genes, which synergistically affected the growth and photosynthesis of poplars. Thus, this study provides genetic resources for the improved breeding of fast-growing tree traits.

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.

GhARF16-1 modulates leaf development by transcriptionally regulating the GhKNOX2-1 gene in cotton

The leaf is a crucial organ evolved with remarkable morphological diversity to maximize plant photosynthesis. The leaf shape is a key trait that affects photosynthesis, flowering rates, disease resistance and yield. Although many genes regulating leaf development have been identified in the past years, the precise regulatory architecture underlying the generation of diverse leaf shapes remains to be elucidated. We used cotton as a reference model to probe the genetic framework underlying divergent leaf forms. Comparative transcriptome analysis revealed that the GhARF16-1 and GhKNOX2-1 genes might be potential regulators of leaf shape. We functionally characterized the auxin-responsive factor ARF16-1 acting upstream of GhKNOX2-1 to determine leaf morphology in cotton. The transcription of GhARF16-1 was significantly higher in lobed-leaved cotton than in smooth-leaved cotton. Furthermore, the overexpression of GhARF16-1 led to the up-regulation of GhKNOX2-1 and resulted in more and deeper serrations in cotton leaves, similar to the leaf shape of cotton plants overexpressing GhKNOX2-1. We found that GhARF16-1 specifically bound to the promoter of GhKNOX2-1 to induce its expression. The heterologous expression of GhARF16-1 and GhKNOX2-1 in Arabidopsis led to lobed and curly leaves, and a genetic analysis revealed that GhKNOX2-1 is epistatic to GhARF16-1 in Arabidopsis, suggesting that the GhARF16-1 and GhKNOX2-1 interaction paradigm also functions to regulate leaf shape in Arabidopsis. To our knowledge, our results uncover a novel mechanism by which auxin, through the key component ARF16-1 and its downstream-activated gene KNOX2-1, determines leaf morphology in eudicots.

Genome-Wide Analysis of Ribosomal Protein GhRPS6 and Its Role in Cotton Verticillium Wilt Resistance

Verticillium wilt is threatening the world’s cotton production. The pathogenic fungus Verticillium dahliae can survive in the soil in the form of microsclerotia for a long time, colonize through the root of cotton, and invade into vascular bundles, causing yellowing and wilting of cotton leaves, and in serious cases, leading to plant death. Breeding resistant varieties is the most economical and effective method to control Verticillium wilt. In previous studies, proteomic analysis was carried out on different cotton varieties inoculated with V. dahliae strain Vd080. It was found that GhRPS6 was phosphorylated after inoculation, and the phosphorylation level in resistant cultivars was 1.5 times than that in susceptible cultivars. In this study, knockdown of GhRPS6 expression results in the reduction of SA and JA content, and suppresses a series of defensive response, enhancing cotton plants susceptibility to V. dahliae. Overexpression in Arabidopsis thaliana transgenic plants was found to be more resistant to V. dahliae. Further, serines at 237 and 240 were mutated to phenylalanine, respectively and jointly. The transgenic Arabidopsis plants demonstrated that seri-237 compromised the plant resistance to V. dahliae. Subcellular localization in Nicotiana benthamiana showed that GhRPS6 was localized in the nucleus. Additionally, the pathogen inoculation and phosphorylation site mutation did not change its localization. These results indicate that GhRPS6 is a potential molecular target for improving resistance to Verticillium wilt in cotton. This lays a foundation for breeding disease-resistant varieties.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-Wide Identification of Wheat WRKY Gene Family Reveals That TaWRKY75-A Is Referred to Drought and Salt Resistances

It is well known that WRKY transcription factors play essential roles in plants' response to diverse stress responses, especially to drought and salt stresses. However, a full comprehensive analysis of this family in wheat is still missing. Here we used in silico analysis and identified 124 WRKY genes, including 294 homeologous copies from a high-quality reference genome of wheat (Triticum aestivum). We also found that the TaWRKY gene family did not undergo gene duplication rather than gene loss during the evolutionary process. The TaWRKY family members displayed different expression profiles under several abiotic stresses, indicating their unique functions in the mediation of particular responses. Furthermore, TaWRKY75-A was highly induced after polyethylene glycol and salt treatments. The ectopic expression of TaWRKY75-A in Arabidopsis enhanced drought and salt tolerance. A comparative transcriptome analysis demonstrated that TaWRKY75-A integrated jasmonic acid biosynthetic pathway and other potential metabolic pathways to increase drought and salt resistances in transgenic Arabidopsis. Our study provides valuable insights into the WRKY family in wheat and will generate a useful genetic resource for improving wheat breeding.

Genomic signatures and candidate genes of lint yield and fibre quality improvement in Upland cotton in Xinjiang

Xinjiang has been the largest and highest yield cotton production region not only in China, but also in the world. Improvements in Upland cotton cultivars in Xinjiang have occurred via pedigree selection and/or crossing of elite alleles from the former Soviet Union and other cotton producing regions of China. But it is unclear how genomic constitutions from foundation parents have been selected and inherited. Here, we deep-sequenced seven historic foundation parents, comprising four cultivars introduced from the former Soviet Union (108Ф, C1470, 611Б and KK1543) and three from United States and Africa (DPL15, STV2B and UGDM), and re-sequenced sixty-nine Xinjiang modern cultivars. Phylogenetic analysis of more than 2 million high-quality single nucleotide polymorphisms allowed their classification two groups, suggesting that Xinjiang Upland cotton cultivars were not only spawned from 108Ф, C1470, 611Б and KK1543, but also had a close kinship with DPL15, STV2B and UGDM. Notably, identity-by-descent (IBD) tracking demonstrated that the former Soviet Union cultivars have made a huge contribution to modern cultivar improvement in Xinjiang. A total of 156 selective sweeps were identified. Among them, apoptosis-antagonizing transcription factor gene (GhAATF1) and mitochondrial transcription termination factor family protein gene (GhmTERF1) were highly involved in the determination of lint percentage. Additionally, the auxin response factor gene (GhARF3) located in inherited IBD segments from 108Ф and 611Б was highly correlated with fibre quality. These results provide an insight into the genomics of artificial selection for improving cotton production and facilitate next-generation precision breeding of cotton and other crops.

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

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

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

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

Genome-Wide Analysis of Cotton Auxin Early Response Gene Families and Their Roles in Somatic Embryogenesis

Auxin is well known to regulate growth and development processes. Auxin early response genes serve as a critical component of auxin signaling and mediate auxin regulation of diverse physiological processes. In the present study, a genome-wide identification and comprehensive analysis of auxin early response genes were conducted in upland cotton. A total of 71 auxin response factor (ARF), 86 Auxin/Indole-3-Acetic Acid (Aux/IAA), 63 Gretchen Hagen3 (GH3), and 194 small auxin upregulated RNA (SAUR) genes were identified in upland cotton, respectively. Phylogenetic analysis revealed that the ARF, GH3, and SAUR families were likely subject to extensive evolutionary divergence between Arabidopsis and upland cotton, while the Aux/IAA family was evolutionary conserved. Expression profiles showed that the ARF, Aux/IAA, GH3, and SAUR family genes were extensively involved in embryogenic competence acquisition of upland cotton callus. The Aux/IAA family genes generally showed a higher expression level in the non-embryogenic callus (NEC) of highly embryogenic cultivar CCRI24 than that of recalcitrant cultivar CCRI12, which may be conducive to initializing the embryogenic transformation. Auxin early response genes were tightly co-expressed with most of the known somatic embryogenesis (SE) related genes, indicating that these genes may regulate upland cotton SE by interacting with auxin early response genes.

Insight into the B3Transcription Factor Superfamily and Expression Profiling of B3 Genes in Axillary Buds after Topping in Tobacco(Nicotiana tabacum L.)

Members of the plant-specific B3 transcription factor superfamily play important roles in various growth and developmental processes in plants. Even though there are many valuable studies on B3 genes in other species, little is known about the B3 superfamily in tobacco. We identified 114 B3 proteins from tobacco using comparative genome analysis. These proteins were classified into four subfamilies based on their phylogenetic relationships, and include the ARF, RAV, LAV, and REM subfamilies. The chromosomal locations, gene structures, conserved protein motifs, and sub-cellular localizations of the tobacco B3 proteins were analyzed. The patterns of exon-intron numbers and arrangement and the protein structures of the tobacco B3 proteins were in general agreement with their phylogenetic relationships. The expression patterns of 114 B3 genes revealed that many B3 genes show tissue-specific expression. The expression levels of B3 genes in axillary buds after topping showed that the REM genes are mainly up-regulated in response to topping, while the ARF genes are down-regulated after topping.

A Pivotal Role of Hormones in Regulating Cotton Fiber Development

Cotton is the main source of renewable fiber in the world and is primarily used for textile production. Cotton fibers are single cells differentiated from the ovule epidermis and are an excellent model system for studying cell elongation, polyploidization, and cell wall biosynthesis. Plant hormones, which are present in relatively low concentrations, play important roles in various developmental processes, and recently, multiple reports have revealed the pivotal roles of hormones in regulating cotton fiber development. For example, exogenous application of hormones has been shown to promote the initiation and growth of fiber cells. However, a comprehensive understanding about phytohormone regulating fiber development is still unknown. Here, we focus on recent advances in elucidating the roles of multiple phytohormones in the control of fiber development, namely auxin, gibberellin, brassinosteroid, ethylene, cytokinin, abscisic acid, and strigolactones. We not only review the identification of genes involved in hormone biosynthetic and signaling pathways but also discuss the mechanisms of these phytohormones in regulating the initiation and elongation of fiber cells in cotton. Auxin, gibberellin, brassinosteroid, ethylene, jasmonic acid, and strigolactones play positive roles in fiber development, whereas cytokinin and abscisic acid inhibit fiber growth. Our aim is to provide a comprehensive review of the role of phytohormones in cotton fiber development that will serve as the basis for further elucidation of the mechanisms by which plant hormones regulate fiber growth.