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

The genetic basis of leaf hair development in upland cotton (Gossypium hirsutum)

Trichomes, which originate from the epidermal cell of aerial organs, provide plants with defense and secretion functions. Although numerous genes have been implicated in trichome development, the molecular mechanisms underlying trichome cell formation in plants remain incompletely understood. Here, we using genome-wide association study (GWAS) across 1037 diverse accessions in upland cotton (Gossypium hirsutum) to identify three loci associated with leaf pubescence (hair) amount, located on chromosome A06 (LPA1), A08 (LPA2) and A11 (LPA3), respectively. GhHD1, a previously characterized candidate gene, was identified on LPA1 and encodes an HD-Zip transcription factor. For LPA2 and LPA3, we identified two candidate genes, GhGIR1 and GhGIR2, both encoding proteins with WD40 and RING domains that act as inhibitors of leaf hair formation. Expression analysis revealed that GhHD1 was predominantly expressed in hairy accessions, whereas GhGIR1 and GhGIR2 were expressed in hairless accessions. Silencing GhHD1 or overexpressing GhGIR1 in hairy accessions induced in a hairless phenotype, whereas silencing GhGIR2 in hairless accessions resulted in a hairy phenotype. We also demonstrated that GhHD1 interact with both GhGIR1 and GhGIR2, and GhGIR1 can interact with GhGIR2. Further investigation indicated that GhHD1 functions as a transcriptional activator, binding to the promoters of the GhGIR1 and GhGIR2 to active their expression, whereas GhGIR1 and GhGIR2 can suppress the transcriptional activation of GhHD1. Our findings shed light on the intricate regulatory network involving GhHD1, GhGIR1 and GhGIR2 in the initiation and development of plant epidermal hairs in cotton.

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Identification and characterization of candidate genes for primary root length in Asiatic cotton (Gossypium arboreum L.)

KEY MESSAGE

One major gene controlling primary root length (PRL) in Gossypium arboreum is identified and this research provides a theoretical basis for root development for cotton. Primary root elongation is an essential process in plant root system structure. Here, we investigated the primary root length (PRL) of 215 diploid cotton (G. arboreum) accessions at 5, 8, 10, 15 days after sowing. A Genome-wide association study was performed for the PRL, resulting in 49 significant SNPs associated with 32 putative candidate genes. The SNP with the strongest signal (Chr07_8047530) could clearly distinguish the PRLs between accessions with two haplotypes. GamurG is the only gene that showed higher relative expression in the long PRL genotypes than the short PRL genotypes, which indicated it was the most likely candidate gene for regulating PRL. Moreover, the GamurG-silenced cotton seedlings showed a shorter PRL, while the GamurG-overexpressed Arabidopsis exhibited a significantly longer PRL. Our findings provide insight into the regulation mechanism of cotton root growth and will facilitate future breeding programs to optimize the root system structure in cotton.

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.

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

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

GoSTR, a negative modulator of stem trichome formation in cotton

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

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.

A HD-ZIP transcription factor specifies fates of multicellular trichomes via dosage-dependent mechanisms in tomato

Hair-like structures are shared by most living organisms. The hairs on plant surfaces, commonly referred to as trichomes, form diverse types to sense and protect against various stresses. However, it is unclear how trichomes differentiate into highly variable forms. Here, we show that a homeodomain leucine zipper (HD-ZIP) transcription factor named Woolly controls the fates of distinct trichomes in tomato via a dosage-dependent mechanism. The autocatalytic reinforcement of Woolly is counteracted by an autoregulatory negative feedback loop, creating a circuit with a high or low Woolly level. This biases the transcriptional activation of separate antagonistic cascades that lead to different trichome types. Our results identify the developmental switch of trichome formation and provide mechanistic insights into the progressive fate specification in plants, as well as a path to enhancing plant stress resistance and the production of beneficial chemicals.

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.

Integrating Genome-wide association and whole transcriptome analysis to reveal genetic control of leaf traits in Gossypium arboreum L

Leaves are important organs for crop photosynthesis and transpiration, and their morphological characteristics can directly reflect the growth state of plants. Accurate measurement of leaf traits and mining molecular markers are of great significance to the study of cotton growth. Here, we performed a Genome-wide association study on 7 leaf traits in 213 Asian cotton accessions. 32 significant SNPs and 44 genes were identified. A field experiment showed significant difference in leaf hair and leaf area between DPL971 and its natural mutant DPL972. We also compared the leaf transcriptome difference between DPL971 and DPL972, and found a batch of differentially expressed genes and non-coding RNAs (including lncRNAs, microRNAs, and circRNAs). After integrating the GWAS and transcriptome results, we finally selected two coding genes (Ga03G2383 and Ga05G3412) and two microRNAs (hbr-miR156, unconservative_Chr03_contig343_2364) as the candidate for leaf traits. Those findings will provide important genomic resources for cotton leaf improvement breeding.

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.

Genome-Wide Analysis of SRNF Genes in Gossypium hirsutum Reveals the Role of GhSRNF18 in Primary Root Growth

Root systems are instrumental for water and nutrient uptake and the anchorage of plants in the soil. Root regulating GL2-interacting repressors (GIRs) contain a Short RING-like Zinc-Finger (SRNF) domain, but there has been no comprehensive characterization about this gene family in any plant species. Here, we renamed the GIR-like proteins as SRNF proteins due to their conserved domain and identified 140 SRNF genes from 16 plant species including 24 GhSRNF genes in Gossypium hirsutum. Phylogenetic analysis of the SRNFs revealed both similarities and divergences between five subfamilies. Notably, synteny analysis revealed that polyploidization and whole-genome duplication contribute to the expansion of the GhSRNF gene family. Various cis-acting regulatory elements were shown to be pertinent to light, phytohormone, defense responsive, and meristem regulation. Furthermore, GhSRNF2/15 were predominantly expressed in root, whereas the expression of GhSRNF18 is positively correlated with the primary root (PR) length in G. hirsutum, quantified by quantitative real-time PCR (qRT-PCR). Over-expression of GhSRNF18 in Arabidopsis and virus-induced gene silencing (VIGS) of GhSRNF18 in G. hirsutum has revealed the role of GhSRNF18 in PR growth. The over-expression of GhSRNF18 in Arabidopsis resulted in an increase of meristematic activities and auxin accumulations in PRs, which were consistent with the transcriptomic data. Our results suggested that GhSRNF18 positively regulates PR growth. This study increased our understanding of the SRNF gene family in plants and provided a novel rationale for the further investigation of cotton root morphogenesis regulated by the GhSRNFs.

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.

Genome-wide association and transcriptome analysis of root color-related genes in Gossypium arboreum L

The significant number loci and candidate genes of root color in Gossypium arboreum are identified and provide a theoretical basis of root color for cotton. A stimulating phenomenon was observed on the 4th day of sowing in the root color of some G. arboreum accessions that turned red. To disclose the genetic mechanisms of root color formation via genome and transcript levels, we identified the significant number of SNPs and candidate genes that are related to root color through genome-wide association study (GWAS) and RNAseq analysis in G. arboreum. Initially, 215 no. of G. arboreum accessions was collected, and the colors of root on the 4th, 6th and 9th day of germination were recorded. The GWAS demonstrated that 225 significant SNPs and 47 candidate genes have been identified totally. The strongest signal SNP A04_91824 could greatly distinguish the root color with most "C" allele accessions have displayed white and "T" allele accessions displayed red. RNAseq was performed on accessions having the white and red root, and results revealed that 12 and 138 DEGs were detected on 2nd and 4th day, respectively. ACD6, UFGT, and LYM2 were the most related genes of root color, later, verified by qRT-PCR. The mature zone of red and the white roots was observed by the histological section method, and results shown that cells were more closely arranged in the white root, and both average cell length and cell width were longer in the red root. This study will be helpful to cotton breeders for utilization of several elite genes and related SNPs related to root color, in addition to find linkage with economically important traits of interests.