The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres

NAC transcription factor GbNTL9 modifies the accumulation and organization of cellulose microfibrils to enhance cotton fiber strength

INTRODUCTION

Fiber strength is a critical determinant of fiber quality, with stronger fibers being highly preferred in the cotton textile industry. However, the genetic basis and the specific regulatory mechanism underlying the formation of cotton fiber strength remain largely unknown.

OBJECTIVES

To explore fiber strength-related genes, QTL mapping, map-based cloning, and gene function verification were conducted in a backcross inbred line BS41 derived from interspecific hybridization between upland cotton and sea-island cotton.

METHODS

Upland cotton Emian22 (E22) and an interspecific backcross inbred line (BIL) BS41 were used as parents to construct secondary segregation populations for BSA and QTL mapping of fiber strength. The candidate gene GbNTL9 was identified through map-based cloning and expression analysis. The function of NTL9 was determined through transgenic experiments and cytological observations. The regulatory mechanisms of NTL9 were explored using RNA-seq, RT-qPCR, yeast two-hybrid, bimolecular fluorescence complementation, and yeast one-hybrid.

RESULTS

A major QTL for fiber strength, qFS-A11-1, was mapped to a 14.6-kb genomic region using segregating populations from E22 × BS41. GbNTL9, which encodes a NAC transcription factor, was identified as the candidate gene. Overexpression of both upland cotton genotype NTL9E22 and sea-island genotype NTL9BS41 in upland cotton enhanced fiber strength by facilitating the dense accumulation and orderly organization of cellulose microfibrils within the cell wall. Transcriptomic analysis revealed that NTL9 inhibited the expression of genes involved in secondary wall synthesis, such as CESA4, CESA7, and CESA8, thereby delaying cell wall cellulose deposition and altering the microfibril deposition pattern. NTL9 interacted with MYB6 and functioned as a downstream gene in the ethylene signaling pathway. Additionally, an effective gene marker NTL9-24 was developed to distinguish haplotypes from G. barbadense and G. hirsutum for fiber quality breeding program.

CONCLUSION

Our findings demonstrate that GbNTL9 positively regulates fiber strength through altering the microfibril deposition pattern, and provide a new insight into the molecular mechanism underlying fiber strength.

Cotton2035: From genomics research to optimized breeding

Cotton is the world's most important natural fiber crop and serves as an ideal model for studying plant genome evolution, cell differentiation, elongation, and cell wall biosynthesis. The first draft genome assembly for Gossypium raimondii, completed in 2012, marked the beginning of global efforts in studying cotton genomics. Over the past decade, the cotton research community has continued to assemble and refine the genomes for both wild and cultivated Gossypium species. With the accumulation of de novo genome assemblies and resequencing data across virous cotton populations, significant progress has been made in uncovering the genetic basis of key agronomic traits. Achieving the goal of cotton genomics-to-breeding (G2B) will require a deeper understanding of the spatiotemporal regulatory mechanisms involved in genome information storage and expression. We advocate for a cotton ENCODE project to systematically decode the functional elements and regulatory networks within the cotton genome. Technological advances, particularly on single-cell sequencing and high-resolution spatiotemporal omics, will be essential for elucidating these regulatory mechanisms. By integrating multi-omics data, genome editing tools, and artificial intelligence, these efforts will empower the genomics-driven strategies needed for future cotton G2B breeding.

Unveiling the Role of GhP5CS1 in Cotton Salt Stress Tolerance: A Comprehensive Genomic and Functional Analysis of P5CS Genes

Proline, a critical osmoregulatory compound, is integral to various plant stress responses. The P5CS gene, which encodes the rate-limiting enzyme in proline biosynthesis, known as ∆1-pyrroline-5-carboxylate synthetase, is fundamental to these stress response pathways. While the functions of P5CS genes in plants have been extensively documented, their specific roles in cotton remain inadequately characterized. In this study, we identified 40 P5CS genes across four cotton species with diverse sequence lengths and molecular weights. Phylogenetic analysis of 100 P5CS genes from nine species revealed three subgroups, with Gossypium hirsutum closely related to Gossypium barbadense. Collinearity analysis highlighted significant differences in collinear gene pairs, indicating evolutionary divergence among P5CS genes in tetraploid and diploid cotton. Exon-intron structures and conserved motifs correlated with phylogenetic relationships, suggesting functional differentiation. Stress-responsive elements in P5CS promoters suggest involvement in abiotic stress. Expression analysis under salt stress revealed differential expressions of GhP5CS genes, with GhP5CS1 emerging as a potential key regulator. Virus-induced gene silencing confirmed the pivotal role of GhP5CS1 in cotton's salt stress response, as evidenced by increased salt sensitivity in the silenced plants. This study enhances our understanding of the functional diversity and roles of P5CS genes in cotton under stress conditions.

A CC-NB-ARC-LRR Gene Regulates Bract Morphology in Cotton

Bracts are leaf-like structures in flowering plants. They serve multiple functions such as attracting pollinators, aiding tolerance of abiotic stressors, and conducting photosynthesis. While previous studies extensively examine bract function, the molecular mechanisms underlying bract growth remain unknown. Here, the map-based isolation and characterization of a crucial factor responsible for cotton bract development, identified from a mutant known as frego bract (fg), discovered by Frego in 1945 are presented. This gene, named Ghfg, encodes a CC-NB-ARC-LRR (CNL) family protein. Through analysis of bract form in plants with virus-induced gene silencing (VIGS) and transgenic plants, this gene is confirmed to be the causal gene under the fg locus. Furthermore, high-resolution single-cell transcriptomic landscape of cotton bracts is generated, which reveals differences related to auxin in proliferating cells from TM-1 and T582; differences in auxin distribution and ROS accumulation are experimentally verified. These findings suggest that GhFG is in a self-activated state in the fg mutant, and its activity leads to ROS accumulation that impacts auxin distribution and transport. Finally, an island cotton variety with the frego bract trait is developed, demonstrating a novel solution for reducing the high impurity rate caused by bract remnants.

Advances in genome sequencing and artificially induced mutation provides new avenues for cotton breeding

Cotton production faces challenges in fluctuating environmental conditions due to limited genetic variation in cultivated cotton species. To enhance the genetic diversity crucial for this primary fiber crop, it is essential to augment current germplasm resources. High-throughput sequencing has significantly impacted cotton functional genomics, enabling the creation of diverse mutant libraries and the identification of mutant functional genes and new germplasm resources. Artificial mutation, established through physical or chemical methods, stands as a highly efficient strategy to enrich cotton germplasm resources, yielding stable and high-quality raw materials. In this paper, we discuss the good foundation laid by high-throughput sequencing of cotton genome for mutant identification and functional genome, and focus on the construction methods of mutant libraries and diverse sequencing strategies based on mutants. In addition, the important functional genes identified by the cotton mutant library have greatly enriched the germplasm resources and promoted the development of functional genomes. Finally, an innovative strategy for constructing a cotton CRISPR mutant library was proposed, and the possibility of high-throughput screening of cotton mutants based on a UAV phenotyping platform was discussed. The aim of this review was to expand cotton germplasm resources, mine functional genes, and develop adaptable materials in a variety of complex environments.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we resequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions was clustered into 3 groups: G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of 5 fiber quality traits across multiple field environments identified a total of 512 qtls (main-effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs colocated with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from Gossypium hirsutum, and the recombination of domesticated elite allelic variation were 3 major contributors to improve the fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variations associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced to G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H is localized in the plasma membrane, while GB_D03G0092B is in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

GhiPLATZ17 and GhiPLATZ22, zinc-dependent DNA-binding transcription factors, promote salt tolerance in upland cotton

KEY MESSAGE

The utilization of transcriptome analysis, functional validation, VIGS, and DAB techniques have provided evidence that GhiPLATZ17 and GhiPLATZ22 play a pivotal role in improving the salt tolerance of upland cotton. PLATZ (Plant AT-rich sequences and zinc-binding proteins) are known to be key regulators in plant growth, development, and response to salt stress. In this study, we comprehensively analyzed the PLATZ family in ten cotton species in response to salinity stress. Gossypium herbaceum boasts 25 distinct PLATZ genes, paralleled by 24 in G. raimondii, 25 in G. arboreum, 46 in G. hirsutum, 48 in G. barbadense, 43 in G. tomentosum, 67 in G. mustelinum, 60 in G. darwinii, 46 in G. ekmanianum, and a total of 53 PLATZ genes attributed to G. stephensii. The PLATZ gene family shed light on the hybridization and allopolyploidy events that occurred during the evolutionary history of allotetraploid cotton. Ka/Ks analysis suggested that the PLATZ gene family underwent intense purifying selection during cotton evolution. Analysis of synteny and gene collinearity revealed a complex pattern of segmental and dispersed duplication events to expand PLATZ genes in cotton. Cis-acting elements and gene expressions revealed that GhiPLATZ exhibited salt stress resistance. Transcriptome analysis, functional validation, virus-induced gene silencing (VIGS), and diaminobenzidine staining (DAB) demonstrated that GhiPLATZ17 and GhiPLATZ22 enhance salt tolerance in upland cotton. The study can potentially advance our understanding of identifying salt-resistant genes in cotton.

High-quality Gossypium hirsutum and Gossypium barbadense genome assemblies reveal the landscape and evolution of centromeres

Centromere positioning and organization are crucial for genome evolution; however, research on centromere biology is largely influenced by the quality of available genome assemblies. Here, we combined Oxford Nanopore and Pacific Biosciences technologies to de novo assemble two high-quality reference genomes for Gossypium hirsutum (TM-1) and Gossypium barbadense (3-79). Compared with previously published reference genomes, our assemblies show substantial improvements, with the contig N50 improved by 4.6-fold and 5.6-fold, respectively, and thus represent the most complete cotton genomes to date. These high-quality reference genomes enable us to characterize 14 and 5 complete centromeric regions for G. hirsutum and G. barbadense, respectively. Our data revealed that the centromeres of allotetraploid cotton are occupied by members of the centromeric repeat for maize (CRM) and Tekay long terminal repeat families, and the CRM family reshapes the centromere structure of the At subgenome after polyploidization. These two intertwined families have driven the convergent evolution of centromeres between the two subgenomes, ensuring centromere function and genome stability. In addition, the repositioning and high sequence divergence of centromeres between G. hirsutum and G. barbadense have contributed to speciation and centromere diversity. This study sheds light on centromere evolution in a significant crop and provides an alternative approach for exploring the evolution of polyploid plants.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

Repression of GhTUBB1 Reduces Plant Height in Gossypium hirsutum

The original 'Green Revolution' genes are associated with gibberellin deficiency. However, in some species, mutations in these genes cause pleiotropic phenotypes, preventing their application in dwarf breeding. The development of novel genotypes with reduced plant height will resolve this problem. In a previous study, we obtained two dwarf lines, L28 and L30, by introducing the Ammopiptanthus mongolicus (Maxim. ex Kom.) Cheng f. C-repeat-binding factor 1 (AmCBF1) into the upland cotton variety R15. We found that Gossypium hirsutum Tubulin beta-1 (GhTUBB1) was downregulated in L28 and L30, which suggested that this gene may have contributed to the dwarf phenotype of L28 and L30. Here, we tested this hypothesis by silencing GhTUBB1 expression in R15 and found that decreased expression resulted in a dwarf phenotype. Interestingly, we found that repressing AmCBF1 expression in L28 and L30 partly recovered the expression of GhTUBB1. Thus, AmCBF1 expression presented a negative relationship with GhTUBB1 expression in L28 and L30. Moreover, yeast one-hybrid and dual-luciferase assays suggest that AmCBF1 negatively regulates GhTUBB1 expression by directly binding to C-repeat/dehydration-responsive (CRT/DRE) elements in the GhTUBB1 promoter, potentially explaining the dwarf phenotypes of L28 and L30. This study elucidates the regulation of GhTUBB1 expression by AmCBF1 and suggests that GhTUBB1 may be a new target gene for breeding dwarf and compact cultivars.

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.

GhIMP10D, an inositol monophosphates family gene, enhances ascorbic acid and antioxidant enzyme activities to confer alkaline tolerance in Gossypium hirsutum L

BACKGROUND

Inositol monophosphates (IMP) are key enzymes in the ascorbic acid (AsA) synthesis pathways, which play vital roles in regulating plant growth and development and stresses tolerance. To date, no comprehensive analysis of the expression profile of IMP genes and their functions under abiotic stress in cotton has been reported.

RESULTS

In this study, the genetic characteristics, phylogenetic evolution, cis-acting elements and expression patterns of IMP gene family in cotton were systematically analyzed. A total of 28, 27, 13 and 13 IMP genes were identified in Gossypium hirsutum (G. hirsutum), Gossypium barbadense (G. barbadense), Gossypium arboreum (G. arboreum), and Gossypium raimondii (G. raimondii), respectively. Phylogenetic analysis showed that IMP family genes could cluster into 3 clades. Structure analysis of genes showed that GhIMP genes from the same subgroup had similar genetic structure and exon number. And most GhIMP family members contained hormone-related elements (abscisic acid response element, MeJA response element, gibberellin response element) and stress-related elements (low temperature response element, defense and stress response element, wound response element). After exogenous application of abscisic acid (ABA), some GhIMP genes containing ABA response elements positively responded to alkaline stress, indicating that ABA response elements played an important role in response to alkaline stress. qRT-PCR showed that most of GhIMP genes responded positively to alkaline stress, and GhIMP10D significantly upregulated under alkaline stress, with the highest up-regulated expression level. Virus-induced gene silencing (VIGS) experiment showed that compared with 156 plants, MDA content of pYL156:GhIMP10D plants increased significantly, while POD, SOD, chlorophyII and AsA content decreased significantly.

CONCLUSIONS

This study provides a thorough overview of the IMP gene family and presents a new perspective on the evolution of this gene family. In particular, some IMP genes may be involved in alkaline stress tolerance regulation, and GhIMP10D showed high expression levels in leaves, stems and roots under alkaline stress, and preliminary functional verification of GhIMP10D gene suggested that it may regulate tolerance to alkaline stress by regulating the activity of antioxidant enzymes and the content of AsA. This study contributes to the subsequent broader discussion of the structure and alkaline resistance of IMP genes in cotton.

Genetic dissection of lint percentage in short-season cotton using combined QTL mapping and RNA-seq

KEY MESSAGE

In total, 17 QTLs for lint percentage in short-season cotton, including three stable QTLs, were detected. Twenty-eight differentially expressed genes located within the stable QTLs were identified, and two genes were validated by qRT-PCR. The breeding and use of short-season cotton have significant values in addressing the question of occupying farmlands with either cotton or cereals. However, the fiber yields of short-season cotton varieties are significantly lower than those of middle- and late-maturing varieties. How to effectively improve the fiber yield of short-season cotton has become a focus of cotton research. Here, a high-density genetic map was constructed using genome resequencing and an RIL population generated from the hybridization of two short-season cotton accessions, Dong3 and Dong4. The map contained 4960 bin markers across the 26 cotton chromosomes and spanned 3971.08 cM, with an average distance of 0.80 cM between adjacent markers. Based on the genetic map, quantitative trait locus (QTL) mapping for lint percentage (LP, %), an important yield component trait, was performed. In total, 17 QTLs for LP, including three stable QTLs, qLP-A02, qLP-D04, and qLP-D12, were detected. Three out of 11 non-redundant QTLs overlapped with previously reported QTLs, whereas the other eight were novel QTLs. A total of 28 differentially expressed genes associated with the three stable QTLs were identified using RNA-seq of ovules and fibers at different seed developmental stages from the parental materials. The two genes, Ghir_A02G017640 and Ghir_A02G018500, may be related to LP as determined by further qRT-PCR validation. This study provides useful information for the genetic dissection of LP and promotes the molecular breeding of short-season cotton.

Genome-Wide Identification and Expression Analysis of RLCK-VII Subfamily Genes Reveal Their Roles in Stress Responses of Upland Cotton

Receptor-like cytoplasmic kinase VII (RLCK-VII) subfamily members are vital players in plant innate immunity and are also involved in plant development and abiotic stress tolerance. As a widely cultivated textile crop, upland cotton (Gossypium hirsutum) attaches great importance to the cotton industry worldwide. To obtain details of the composition, phylogeny, and putative function of RLCK-VII genes in upland cotton, genome-wide identification, evolutionary event analysis, and expression pattern examination of RLCK-VII subfamily genes in G. hirsutum were performed. There are 129 RLCK-VII members in upland cotton (GhRLCKs) and they were divided into nine groups based on their phylogenetic relationships. The gene structure and sequence features are relatively conserved within each group, which were divided based on their phylogenetic relationships, and consistent with those in Arabidopsis. The phylogenetic analysis results showed that RLCK-VII subfamily genes evolved in plants before the speciation of Arabidopsis and cotton, and segmental duplication was the major factor that caused the expansion of GhRLCKs. The diverse expression patterns of GhRLCKs in response to abiotic stresses (temperature, salt, and drought) and V. dahliae infection were observed. The candidates that may be involved in cotton's response to these stresses are highlighted. GhRLCK7 (GhRLCK7A and D), which is notably induced by V. dahliae infection, was demonstrated to positively regulate cotton defense against V. dahliae by the loss-of-function assay in cotton. These findings shed light on the details of the RLCK-VII subfamily in cotton and provide a scaffold for the further function elucidation and application of GhRLCKs for the germplasm innovation of cotton.

UAV-based time-series phenotyping reveals the genetic basis of plant height in upland cotton

Plant height (PH) is an important agronomic trait affecting crop architecture, biomass, resistance to lodging and mechanical harvesting. Elucidating the genetic governance of plant height is crucial because of the global demand for high crop yields. However, during the rapid growth period of plants the PH changes a lot on a daily basis, which makes it difficult to accurately phenotype the trait by hand on a large scale. In this study, an unmanned aerial vehicle (UAV)-based remote-sensing phenotyping platform was applied to obtain time-series PHs of 320 upland cotton accessions in three different field trials. The results showed that the PHs obtained from UAV images were significantly correlated with ground-based manual measurements, for three trials (R2  = 0.96, 0.95 and 0.96). Two genetic loci on chromosomes A01 and A11 associated with PH were identified by genome-wide association studies (GWAS). GhUBP15 and GhCUL1 were identified to influence PH in further analysis. We obtained a time series of PH values for three field conditions based on remote sensing with UAV. The key genes identified in this study are of great value for the breeding of ideal plant architecture in cotton.

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 identification and functional analysis of ICE genes reveal that Gossypium thurberi "GthICE2" is responsible for cold and drought stress tolerance

Cold stress has been found to have a negative impact on cotton growth and annual production. To address this issue, the utilization of cold-tolerant gene resources from wild species of Gossypium is crucial for genetic improvements in cultivated cotton. ICE (inducer of CBF expression) are the key regulators of cold tolerance in plants, however, there is relatively little information on ICE genes in cotton. Herein, we performed comprehensive bioinformatics analyses of the ICE gene family in eight cotton species. Phylogenetic analysis showed that 52 ICE genes were clustered into four subgroups. Cis-regulatory elements analysis suggests that the expression of ICE genes might be regulated by light, plant hormones, and various environment stresses. Higher expression of GthICE2 was observed in leaves as compared to roots and stems, in response to cold, drought, and exogenous hormone ABA. Furthermore, overexpression of GthICE2 in A. thaliana led to higher germination and survival rates, longer root length, lower ion leakage, and induction under cold and drought stress. Histochemical staining showed that oxidative damage in transgenic lines was much lower compared to wild-type plants. Lower MDA contents and higher SOD and POD activities were observed in overexpressed plants. Y1H and LUC assays revealed that GthICE2 might activate the expression of GthCBF4, a cold-responsive gene, by connecting with the MYC cis-element present in the promoter of GthCBF4. GthICE2 confers cold and drought stress tolerance in cotton. Our findings add significantly to the existing knowledge regarding cold stress tolerance and helps to elucidate cold response mechanisms in cotton.

Identification and analysis of isoflavone reductase gene family in Gossypium hirsutum L

Isoflavone reductase (IFR) is a key enzyme controlling isoflavone synthesis and widely involved in response to various stresses. In this study, the IFR genes in four Gossypium species and other 7 species were identified and analyzed in the whole genome, and the physicochemical properties, gene structures, cis-acting elements, chromosomal locations, collinearity relationships and expression patterns of IFR genes were systematically analyzed. 28, 28, 14 and 15 IFR genes were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii, respectively, which were divided into five clades according to the evolutionary tree and gene structure. Collinear analysis showed that segmental duplication and whole genome duplication were the main driving forces in the process of evolution, and most genes underwent pure selection. Gene structure analysis showed that IFR gene family was relatively conserved. Cis-element analysis of promoter showed that most GhIFR genes contain cis-elements related to abiotic stresses and plant hormones. Analysis of GhIFR gene expression under different stresses showed that GhIFR genes were involved in the response to drought, salt, heat and cold stresses through corresponding network mechanisms, especially GhIFR9A. Phenotypic analysis after silencing GhIFR9A gene by VIGS was shown that GhIFR9A gene was involved in the response to salt stress. This study laid a foundation for the subsequent functional study of cotton IFR genes.

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.

Assembly and annotation of the Gossypium barbadense L. 'Pima-S6' genome raise questions about the chromosome structure and gene content of Gossypium barbadense genomes

BACKGROUND

Gossypium barbadense L. Pima cotton is known for its resistance to Fusarium wilt and for producing fibers of superior quality highly prized in the textile market. We report a high-quality genome assembly and annotation of Pima-S6 cotton and its comparison at the chromosome and protein level to other ten Gossypium published genome assemblies.

RESULTS

Synteny and orthogroup analyses revealed important differences on chromosome structure and annotated proteins content between our Pima-S6 and other publicly available G. barbadense assemblies, and across Gossypium assemblies in general. Detailed synteny analyses revealed chromosomal rearrangements between Pima-S6 and other Pima genomes on several chromosomes, with three major inversions in chromosomes A09, A13 and D05, raising questions about the true chromosome structure of Gossypium barbadense genomes.

CONCLUSION

Analyses of the re-assembled and re-annotated genome of the close relative G. barbadense Pima 3-79 using our Pima-S6 assembly suggest that contig placement of some recent G. barbadense assemblies might have been unduly influenced by the use of the G. hirsutum TM-1 genome as the anchoring reference. The Pima-S6 reference genome provides a valuable genomic resource and offers new insights on genomic structure, and can serve as G. barbadense genome reference for future assemblies and further support FOV4-related studies and breeding efforts.

A DEK domain-containing protein GhDEK2D mediated Gossypium hirsutum enhanced resistance to Verticillium dahliae

DEK is associated with DNA replication and break repair, mRNA splicing, and transcriptional regulation, which had been studied in humans and mammals. The function of DEK in plants was poorly understood. In this study, GhDEK2D was identified in Gossypium hirsutum by genome-wide and post-translational modifications. GhDEK2D had been phosphorylated, acetylated and ubiquitylated under Verticillium dahliae (Vd) challenge. The GhDEK2D-silenced cotton decreased resistance against Vd. In GhDEK2D-silenced cotton plants, the reactive oxygen species was activated, the callose, xylogen, hypersensitive reaction (HR) and expression levels of defense-related genes were reduced. Homozygous overexpressing-GhDEK2D transgenic Arabidopsis lines were more resistant to Verticillium wilt (Vw). We propose that GhDEK2D was a potential molecular target for improving resistance to Vw in cotton.

Highly efficient Agrobacterium rhizogenes-mediated hairy root transformation for gene editing analysis in cotton

CRIPSR/Cas9 gene editing system is an effective tool for genome modification in plants. Multiple target sites are usually designed and the effective target sites are selected for editing. Upland cotton (Gossypium hirsutum L., hereafter cotton) is allotetraploid and is commonly considered as difficult and inefficient to transform, it is important to select the effective target sites that could result in the ideal transgenic plants with the CRISPR-induced mutations. In this study, Agrobacterium rhizogenes-mediated hairy root method was optimized to detect the feasibility of the target sites designed in cotton phytoene desaturase (GhPDS) gene. A. rhizogenes showed the highest hairy root induction (30%) when the bacteria were cultured until OD₆₀₀ reached to 0.8. This procedure was successfully applied to induce hairy roots in the other three cultivars (TM-1, Lumian-21, Zhongmian-49) and the mutations were detected in GhPDS induced by CRISPR/Cas9 system. Different degrees of base deletions at two sgRNAs (sgRNA5 and sgRNA10) designed in GhPDS were detected in R15 hairy roots. Furthermore, we obtained an albino transgenic cotton seeding containing CRISPR/Cas9-induced gene editing mutations in sgRNA10. The hairy root transformation system established in this study is sufficient for selecting sgRNAs in cotton, providing a technical basis for functional genomics research of cotton.

Comparative Genomics and Functional Studies of Putative m6A Methyltransferase (METTL) Genes in Cotton

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Analysis of transcriptome data and quantitative trait loci enables the identification of candidate genes responsible for fiber strength in Gossypium barbadense

Gossypium barbadense possesses a superior fiber quality because of its fiber length and strength. An in-depth analysis of the underlying genetic mechanism could aid in filling the gap in research regarding fiber strength and could provide helpful information for Gossypium barbadense breeding. Three quantitative trait loci related to fiber strength were identified from a Gossypium barbadense recombinant inbred line (PimaS-7 × 5917) for further analysis. RNA sequencing was performed in the fiber tissues of PimaS-7 × 5917 0–35 days postanthesis. Four specific modules closely related to the secondary wall-thickening stage were obtained using the weighted gene coexpression network analysis. In total, 55 genes were identified as differentially expressed from 4 specific modules. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for enrichment analysis, and Gbar_D11G032910, Gbar_D08G020540, Gbar_D08G013370, Gbar_D11G033670, and Gbar_D11G029020 were found to regulate fiber strength by playing a role in the composition of structural constituents of cytoskeleton and microtubules during fiber development. Quantitative real-time PCR results confirmed the accuracy of the transcriptome data. This study provides a quick strategy for exploring candidate genes and provides new insights for improving fiber strength in cotton.

Systematic analysis of CNGCs in cotton and the positive role of GhCNGC32 and GhCNGC35 in salt tolerance

BACKGROUND

Cyclic nucleotide-gated ion channels (CNGCs) are calcium-permeable channels that participate in a variety of biological functions, such as signaling pathways, plant development, and environmental stress and stimulus responses. Nevertheless, there have been few studies on CNGC gene family in cotton.

RESULTS

In this study, a total of 114 CNGC genes were identified from the genomes of 4 cotton species. These genes clustered into 5 main groups: I, II, III, IVa, and IVb. Gene structure and protein motif analysis showed that CNGCs on the same branch were highly conserved. In addition, collinearity analysis showed that the CNGC gene family had expanded mainly by whole-genome duplication (WGD). Promoter analysis of the GhCNGCs showed that there were a large number of cis-acting elements related to abscisic acid (ABA). Combination of transcriptome data and the results of quantitative RT-PCR (qRT-PCR) analysis revealed that some GhCNGC genes were induced in response to salt and drought stress and to exogenous ABA. Virus-induced gene silencing (VIGS) experiments showed that the silencing of the GhCNGC32 and GhCNGC35 genes decreased the salt tolerance of cotton plants (TRV:00). Specifically, physiological indexes showed that the malondialdehyde (MDA) content in gene-silenced plants (TRV:GhCNGC32 and TRV:GhCNGC35) increased significantly under salt stress but that the peroxidase (POD) activity decreased. After salt stress, the expression level of ABA-related genes increased significantly, indicating that salt stress can trigger the ABA signal regulatory mechanism.

CONCLUSIONS

we comprehensively analyzed CNGC genes in four cotton species, and found that GhCNGC32 and GhCNGC35 genes play an important role in cotton salt tolerance. These results laid a foundation for the subsequent study of the involvement of cotton CNGC genes in salt tolerance.

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.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

lncRNA7 and lncRNA2 modulate cell wall defense genes to regulate cotton resistance to Verticillium wilt

In plants, long noncoding RNAs (lncRNAs) regulate disease resistance against fungi and other pathogens. However, the specific mechanism behind this regulation remains unclear. In this study, we identified disease resistance-related lncRNAs as well as their regulating genes and assessed their functions by infection of cotton (Gossypium) chromosome segment substitution lines with Verticillium dahliae. Our results demonstrated that lncRNA7 and its regulating gene Pectin methylesterase inhibitor 13 (GbPMEI13) positively regulated disease resistance via the silencing approach, while ectopic overexpression of GbPMEI13 in Arabidopsis (Arabidopsis thaliana) promoted growth and enhanced resistance to V. dahliae. In contrast, lncRNA2 and its regulating gene Polygalacturonase 12 (GbPG12) negatively regulated resistance to V. dahliae. We further found that fungal disease-related agents, including the pectin-derived oligogalacturonide (OG), could downregulate the expression of lncRNA2 and GbPG12, leading to pectin accumulation. Conversely, OG upregulated the expression of lncRNA7, which encodes a plant peptide phytosulfokine (PSK-α), which was confirmed by lncRNA7 overexpression and Ultra Performance Liquid Chromatography Tandem Mass Spectrometry (UPLC-MS) experiments. We showed that PSK-α promoted 3-Indoleacetic acid (IAA) accumulation and activated GbPMEI13 expression through Auxin Response Factor 5. Since it is an inhibitor of pectin methylesterase (PME), GbPMEI13 promotes pectin methylation and therefore increases the resistance to V. dahliae. Consistently, we also demonstrated that GbPMEI13 inhibits the mycelial growth and spore germination of V. dahliae in vitro. In this study, we demonstrated that lncRNA7, lncRNA2, and their regulating genes modulate cell wall defense against V. dahliae via auxin-mediated signaling, providing a strategy for cotton breeding.

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

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

Genomic and GWAS analyses demonstrate phylogenomic relationships of Gossypium barbadense in China and selection for fibre length, lint percentage and Fusarium wilt resistance

Sea Island cotton (Gossypium barbadense) is the source of the world's finest fibre quality cotton, yet relatively little is understood about genetic variations among diverse germplasms, genes underlying important traits and the effects of pedigree selection. Here, we resequenced 336 G. barbadense accessions and identified 16 million SNPs. Phylogenetic and population structure analyses revealed two major gene pools and a third admixed subgroup derived from geographical dissemination and interbreeding. We conducted a genome-wide association study (GWAS) of 15 traits including fibre quality, yield, disease resistance, maturity and plant architecture. The highest number of associated loci was for fibre quality, followed by disease resistance and yield. Using gene expression analyses and VIGS transgenic experiments, we confirmed the roles of five candidate genes regulating four key traits, that is disease resistance, fibre length, fibre strength and lint percentage. Geographical and temporal considerations demonstrated selection for the superior fibre quality (fibre length and fibre strength), and high lint percentage in improving G. barbadense in China. Pedigree selection breeding increased Fusarium wilt disease resistance and separately improved fibre quality and yield. Our work provides a foundation for understanding genomic variation and selective breeding of Sea Island cotton.

Inheritance, QTLs, and Candidate Genes of Lint Percentage in Upland Cotton

Cotton (Gossypium spp.) is an important natural fiber plant. Lint percentage (LP) is one of the most important determinants of cotton yield and is a typical quantitative trait with high variation and heritability. Many cotton LP genetic linkages and association maps have been reported. This work summarizes the inheritance, quantitative trait loci (QTLs), and candidate genes of LP to facilitate LP genetic study and molecular breeding. More than 1439 QTLs controlling LP have been reported. Excluding replicate QTLs, 417 unique QTLs have been identified on 26 chromosomes, including 243 QTLs identified at LOD >3. More than 60 are stable, major effective QTLs that can be used in marker-assisted selection (MAS). More than 90 candidate genes for LP have been reported. These genes encode MYB, HOX, NET, and other proteins, and most are preferentially expressed during fiber initiation and elongation. A putative molecular regulatory model of LP was constructed and provides the foundation for the genetic study and molecular breeding of LP.

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BACKGROUND

The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.

RESULTS

In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.

CONCLUSION

We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

Development of an efficient and precise adenine base editor (ABE) with expanded target range in allotetraploid cotton (Gossypium hirsutum)

BACKGROUND

Base editors (BEs) display diverse applications in a variety of plant species such as Arabidopsis, rice, wheat, maize, soybean, and cotton, where they have been used to mediate precise base pair conversions without the collateral generation of undesirable double-stranded breaks (DSB). Studies of single-nucleotide polymorphisms (SNPs) underpinning plant traits are still challenging, particularly in polyploidy species where such SNPs are present in multiple copies, and simultaneous modification of all alleles would be required for functional analysis. Allotetraploid cotton has a number of homoeologous gene pairs located in the A and D sub-genomes with considerable SNPs, and it is desirable to develop adenine base editors (ABEs) for efficient and precise A-to-G single-base editing without DSB in such complex genome.

RESULTS

We established various ABE vectors based on different engineered adenosine deaminase (TadA) proteins fused to Cas9 variants (dCas9, nCas9), enabling efficient A to G editing up to 64% efficiency on-target sites of the allotetraploid cotton genome. Comprehensive analysis showed that GhABE7.10n exhibited the highest editing efficiency, with the main editing sites specifically located at the position A5 (counting the PAM as positions 21-23). Furthermore, DNA and RNA off-target analysis of cotton plants edited with GhABE7.10n and GhABE7.10d by whole genome and whole-transcriptome sequencing revealed no DNA off-target mutations, while very low-level RNA off-target mutations were detected. A new base editor, namely GhABE7.10dCpf1 (7.10TadA + dCpf1), that recognizes a T-rich PAM, was developed for the first time. Targeted A-to-G substitutions generated a single amino acid change in the cotton phosphatidyl ethanolamine-binding protein (GhPEBP), leading to a compact cotton plant architecture, an ideotype for mechanized harvesting of modern cotton production.

CONCLUSIONS

Our data illustrate the robustness of adenine base editing in plant species with complex genomes, which provides efficient and precise toolkit for cotton functional genomics and precise molecular breeding.

A Calmodulin-Like Gene (GbCML7) for Fiber Strength and Yield Improvement Identified by Resequencing Core Accessions of a Pedigree in Gossypium barbadense

Sea Island cotton (Gossypium barbadense) is world-renowned for its superior natural fiber. Although fiber strength is one of the most important fiber quality traits, genes contributing to fiber strength are poorly understood. Production of sea island cotton also is inextricably linked to improving its relatively low yield, thus enhancing the importance of joint improvement of both fiber quality and yield. We used genomic variation to uncover the genetic evidence of trait improvement resulting from pedigree breeding of Sea Island cotton. This pedigree was aimed at improving fiber strength and yielded an elite cultivar, XH35. Using a combination of genome-wide association study (GWAS) and selection screens, we detected 82 putative fiber-strength-related genes. Expression analysis confirmed a calmodulin-like gene, GbCML7, which enhanced fiber strength in a specific haplotype. This gene is a major-effect gene, which interacts with a minor-effect gene, GbTUA3, facilitating the enhancement of fiber strength in a synergistic fashion. Moreover, GbCML7 participates in the cooperative improvement of fiber strength, fiber length, and fiber uniformity, though a slight compromise exists between the first two of these traits and the latter. Importantly, GbCML7 is shown to boost yield in some backgrounds by increasing multiple yield components to varying degrees, especially boll number. Our work provides valuable genomic evidence and a key genetic factor for the joint improvement of fiber quality and yield in Sea Island cotton.

Dynamic 3D genome architecture of cotton fiber reveals subgenome-coordinated chromatin topology for 4-staged single-cell differentiation

BACKGROUND

Despite remarkable advances in our knowledge of epigenetically mediated transcriptional programming of cell differentiation in plants, little is known about chromatin topology and its functional implications in this process.

RESULTS

To interrogate its significance, we establish the dynamic three-dimensional (3D) genome architecture of the allotetraploid cotton fiber, representing a typical single cell undergoing staged development in plants. We show that the subgenome-relayed switching of the chromatin compartment from active to inactive is coupled with the silencing of developmentally repressed genes, pinpointing subgenome-coordinated contribution to fiber development. We identify 10,571 topologically associating domain-like (TAD-like) structures, of which 25.6% are specifically organized in different stages and 75.23% are subject to partition or fusion between two subgenomes. Notably, dissolution of intricate TAD-like structure cliques showing long-range interactions represents a prominent characteristic at the later developmental stage. Dynamic chromatin loops are found to mediate the rewiring of gene regulatory networks that exhibit a significant difference between the two subgenomes, implicating expression bias of homologous genes.

CONCLUSIONS

This study sheds light on the spatial-temporal asymmetric chromatin structures of two subgenomes in the cotton fiber and offers a new insight into the regulatory orchestration of cell differentiation in plants.

Analyzing Quantitative Trait Loci for Fiber Quality and Yield-Related Traits From a Recombinant Inbred Line Population With Gossypium hirsutum Race palmeri as One Parent

Fiber quality and yield-related traits are important agronomic traits in cotton breeding. To detect the genetic basis of fiber quality and yield related traits, a recombinant inbred line (RIL) population consisting of 182 lines was established from a cross between Gossypium hirsutum cultivar CCRI35 and G. hirsutum race palmeri accession TX-832. The RIL population was deeply genotyped using SLAF-seq and was phenotyped in six environments. A high-density genetic linkage map with 15,765 SNP markers and 153 SSR markers was constructed, with an average distance of 0.30 cM between adjacent markers. A total of 210 fiber quality quantitative trait loci (QTLs) and 73 yield-related QTLs were identified. Of the detected QTLs, 62 fiber quality QTLs and 10 yield-related QTLs were stable across multiple environments. Twelve and twenty QTL clusters were detected on the At and Dt subgenome, respectively. Twenty-three major QTL clusters were further validated through associated analysis and five candidate genes of four stable fiber quality QTLs were identified. This study revealed elite loci influencing fiber quality and yield and significant phenotypic selection regions during G. hirsutum domestication, and set a stage for future utilization of molecular marker assisted breeding in cotton breeding programs.

CottonGVD: A Comprehensive Genomic Variation Database for Cultivated Cottons

Cultivated cottons are the most important economic crop, which produce natural fiber for the textile industry. In recent years, the genetic basis of several essential traits for cultivated cottons has been gradually elucidated by decoding their genomic variations. Although an abundance of resequencing data is available in public, there is still a lack of a comprehensive tool to exhibit the results of genomic variations and genome-wide association study (GWAS). To assist cotton researchers in utilizing these data efficiently and conveniently, we constructed the cotton genomic variation database (CottonGVD; http://120.78.174.209/ or http://db.cngb.org/cottonGVD). This database contains the published genomic information of three cultivated cotton species, the corresponding population variations (SNP and InDel markers), and the visualized results of GWAS for major traits. Various built-in genomic tools help users retrieve, browse, and query the variations conveniently. The database also provides interactive maps (e.g., Manhattan map, scatter plot, heatmap, and linkage disequilibrium block) to exhibit GWAS and expression GWAS results. Cotton researchers could easily focus on phenotype-associated loci visualization, and they are interested in and screen for candidate genes. Moreover, CottonGVD will continue to update by adding more data and functions.

CottonGen: The Community Database for Cotton Genomics, Genetics, and Breeding Research

Over the last eight years, the volume of whole genome, gene expression, SNP genotyping, and phenotype data generated by the cotton research community has exponentially increased. The efficient utilization/re-utilization of these complex and large datasets for knowledge discovery, translation, and application in crop improvement requires them to be curated, integrated with other types of data, and made available for access and analysis through efficient online search tools. Initiated in 2012, CottonGen is an online community database providing access to integrated peer-reviewed cotton genomic, genetic, and breeding data, and analysis tools. Used by cotton researchers worldwide, and managed by experts with crop-specific knowledge, it continuous to be the logical choice to integrate new data and provide necessary interfaces for information retrieval. The repository in CottonGen contains colleague, gene, genome, genotype, germplasm, map, marker, metabolite, phenotype, publication, QTL, species, transcriptome, and trait data curated by the CottonGen team. The number of data entries housed in CottonGen has increased dramatically, for example, since 2014 there has been an 18-fold increase in genes/mRNAs, a 23-fold increase in whole genomes, and a 372-fold increase in genotype data. New tools include a genetic map viewer, a genome browser, a synteny viewer, a metabolite pathways browser, sequence retrieval, BLAST, and a breeding information management system (BIMS), as well as various search pages for new data types. CottonGen serves as the home to the International Cotton Genome Initiative, managing its elections and serving as a communication and coordination hub for the community. With its extensive curation and integration of data and online tools, CottonGen will continue to facilitate utilization of its critical resources to empower research for cotton crop improvement.

Quantitative Trait Loci and Transcriptome Analysis Reveal Genetic Basis of Fiber Quality Traits in CCRI70 RIL Population of Gossypium hirsutum

Upland cotton (Gossypium hirsutum) is widely planted around the world for its natural fiber, and producing high-quality fiber is essential for the textile industry. CCRI70 is a hybrid cotton plant harboring superior yield and fiber quality, whose recombinant inbred line (RIL) population was developed from two upland cotton varieties (sGK156 and 901-001) and were used here to investigate the source of high-quality related alleles. Based on the material of the whole population, a high-density genetic map was constructed using specific locus-amplified fragment sequencing (SLAF-seq). It contained 24,425 single nucleotide polymorphism (SNP) markers, spanning a distance of 4,850.47 centimorgans (cM) over 26 chromosomes with an average marker interval of 0.20 cM. In evaluating three fiber quality traits in nine environments to detect multiple environments stable quantitative trait loci (QTLs), we found 289 QTLs, of which 36 of them were stable QTLs and 18 were novel. Based on the transcriptome analysis for two parents and two RILs, 24,941 unique differentially expressed genes (DEGs) were identified, 473 of which were promising genes. For the fiber strength (FS) QTLs, 320 DEGs were identified, suggesting that pectin synthesis, phenylpropanoid biosynthesis, and plant hormone signaling pathways could influence FS, and several transcription factors may regulate fiber development, such as GAE6, C4H, OMT1, AFR18, EIN3, bZIP44, and GAI. Notably, the marker D13_56413025 in qFS-chr18-4 provides a potential basis for enhancing fiber quality of upland cotton via marker-assisted breeding and gene cloning of important fiber quality traits.

Genome-wide expression analysis of phospholipase A1 (PLA1) gene family suggests phospholipase A1-32 gene responding to abiotic stresses in cotton

Cotton is the most important crop for the production of natural fibres used in the textile industry. High salinity, drought, cold and high temperature represent serious abiotic stresses, which seriously threaten cotton production. Phospholipase A_S has an irreplaceable role in lipid signal transmission, growth and development and stress events. Phospholipase A can be divided into three families: PLA1, PLA2 and pPLA. Among them, the PLA1 family is rarely studied in plants. In order to study the potential functions of the PLA1 family in cotton, the bioinformatics analysis of the PLA1 family was correlated with cotton adversity, and tissue-specific analysis was performed. Explore the structure-function relationship of PLA1 members. It is found that the expression of GbPLA1-32 gene is affected by a variety of environmental stimuli, indicating that it plays a very important role in stress and hormone response, and closely associates the cotton adversity with this family. Through further functional verification, we found that virus-induced GbPLA1-32 gene silencing (VIGS) caused Gossypium barbadense to be sensitive to salt stress. This research provides an important basis for further research on the molecular mechanism of cotton resistance to abiotic stress.

Identification of Raf-Like Kinases B Subfamily Genes in Gossypium Species Revealed GhRAF42 Enhanced Salt Tolerance in Cotton

Salinity is a critical abiotic factor that significantly reduces agricultural production. Cotton is an important fiber crop and a pioneer on saline soil, hence genetic architecture that underpins salt tolerance should be thoroughly investigated. The Raf-like kinase B-subfamily (RAF) genes were discovered to regulate the salt stress response in cotton plants. However, understanding the RAFs in cotton, such as Enhanced Disease Resistance 1 and Constitutive Triple Response 1 kinase, remains a mystery. This study obtained 29, 28, 56, and 54 RAF genes from G. arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively. The RAF gene family described allopolyploidy and hybridization events in allotetraploid cotton evolutionary connections. Ka/Ks analysis advocates that cotton evolution was subjected to an intense purifying selection of the RAF gene family. Interestingly, integrated analysis of synteny and gene collinearity suggested dispersed and segmental duplication events involved in the extension of RAFs in cotton. Transcriptome studies, functional validation, and virus-induced gene silencing on salt treatments revealed that GhRAF42 is engaged in salt tolerance in upland cotton. This research might lead to a better understanding of the role of RAFs in plants and the identification of suitable candidate salt-tolerant genes for cotton breeding.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

The Gossypium anomalum genome as a resource for cotton improvement and evolutionary analysis of hybrid incompatibility

Cotton is an important crop that has been the beneficiary of multiple genome sequencing efforts, including diverse representatives of wild species for germplasm development. Gossypium anomalum is a wild African diploid species that harbors stress-resistance and fiber-related traits with potential application to modern breeding efforts. In addition, this species is a natural source of cytoplasmic male sterility and a resource for understanding hybrid lethality in the genus. Here, we report a high-quality de novo genome assembly for G. anomalum and characterize this genome relative to existing genome sequences in cotton. In addition, we use the synthetic allopolyploids 2(A2D1) and 2(A2D3) to discover regions in the G. anomalum genome potentially involved in hybrid lethality, a possibility enabled by introgression of regions homologous to the D3 (Gossypium davidsonii) lethality loci into the synthetic 2(A2D3) allopolyploid.

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log₁₀(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

Structural insights into homotrimeric assembly of cellulose synthase CesA7 from Gossypium hirsutum

Cellulose is one of the most abundant organic polymers in nature. It contains multiple β-1,4-glucan chains synthesized by cellulose synthases (CesAs) on the plasma membrane of higher plants. CesA subunits assemble into a pseudo-sixfold symmetric cellulose synthase complex (CSC), known as a 'rosette complex'. The structure of CesA remains enigmatic. Here, we report the cryo-EM structure of the homotrimeric CesA7 from Gossypium hirsutum at 3.5-angstrom resolution. The GhCesA7 homotrimer shows a C3 symmetrical assembly. Each protomer contains seven transmembrane helices (TMs) which form a channel potentially facilitating the release of newly synthesized glucans. The cytoplasmic glycosyltransferase domain (GT domain) of GhCesA7 protrudes from the membrane, and its catalytic pocket is directed towards the TM pore. The homotrimer GhCesA7 is stabilized by the transmembrane helix 7 (TM7) and the plant-conserved region (PCR) domains. It represents the building block of CSCs and facilitates microfibril formation. This structure provides insight into how eukaryotic cellulose synthase assembles and provides a mechanistic basis for the improvement of cotton fibre quality in the future.

Genome-Wide Comparative Analysis of Flowering-Time Genes; Insights on the Gene Family Expansion and Evolutionary Perspective

In polyploids, whole genome duplication (WGD) played a significant role in genome expansion, evolution and diversification. Many gene families are expanded following polyploidization, with the duplicated genes functionally diversified by neofunctionalization or subfunctionalization. These mechanisms may support adaptation and have likely contributed plant survival during evolution. Flowering time is an important trait in plants, which affects critical features, such as crop yields. The flowering-time gene family is one of the largest expanded gene families in plants, with its members playing various roles in plant development. Here, we performed genome-wide identification and comparative analysis of flowering-time genes in three palnt families i.e., Malvaceae, Brassicaceae, and Solanaceae, which indicate these genes were expanded following the event/s of polyploidization. Duplicated genes have been retained during evolution, although genome reorganization occurred in their flanking regions. Further investigation of sequence conservation and similarity network analyses provide evidence for functional diversification of duplicated genes during evolution. These functionally diversified genes play important roles in plant development and provide advantages to plants for adaptation and survival in response to environmental changes encountered during evolution. Collectively, we show that flowering-time genes were expanded following polyploidization and retained as large gene family by providing advantages from functional diversification during evolution.

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.

Genome-wide identification of the BASS gene family in four Gossypium species and functional characterization of GhBASSs against salt stress

Bile acid sodium symporter (BASS) family proteins encode a class of sodium/solute symporters. Even though the sodium transporting property of BASSs in mammals was well studied, their sodium transportability and functional roles in plant salt tolerance remained largely unknown. Here, BASS family members from 4 cotton species, as well as 30 other species were identified. Then, they were designated as members of BASS1 to BASS5 subfamilies according to their sequence similarity and phylogenetic relationships. There were 8, 11, 16 and 18 putative BASS genes in four cotton species. While whole-genome duplications (WGD) and segmental duplications rendered the expansion of the BASS gene family in cotton, BASS gene losses occurred in the tetraploid cotton during the evolution from diploids to allotetraploids. Concerning functional characterizations, the transcript profiling of GhBASSs revealed that they not only preferred tissue-specific expression but also were differently induced by various stressors and phytohormones. Gene silencing and overexpression experiments showed that GhBASS1 and GhBASS3 positively regulated, whereas GhBASS2, GhBASS4 and GhBASS5 negatively regulated plant salt tolerance. Taken together, BASS family genes have evolved before the divergence from the common ancestor of prokaryotes and eukaryotes, and GhBASSs are plastidial sodium-dependent metabolite co-transporters that can influence plant salt tolerance.

Genome-wide association analysis reveals quantitative trait loci and candidate genes involved in yield components under multiple field environments in cotton (Gossypium hirsutum)

BACKGROUND

Numerous quantitative trait loci (QTLs) and candidate genes associated with yield-related traits have been identified in cotton by genome-wide association study (GWAS) analysis. However, most of the phenotypic data were from a single or few environments, and the stable loci remained to be validated under multiple field environments.

RESULTS

Here, 242 upland cotton accessions collected from different origins were continuously investigated for phenotypic data of four main yield components, including boll weight (BW) and lint percentage (LP) under 13 field environments, and boll number per plant (BN) and seed index (SI) under 11 environments. Correlation analysis revealed a positive correlation between BN and LP, BW and SI, while SI had a negative correlation with LP and BN. Genetic analysis indicated that LP had the highest heritability estimates of 94.97%, followed by 92.08% for SI, 86.09% for BW, and 72.92% for BN, indicating LP and SI were more suitable traits for genetic improvement. Based on 56,010 high-quality single nucleotide polymorphisms (SNPs) and GWAS analysis, a total of 95 non-redundant QTLs were identified, including 12 of BN, 23 of BW, 45 of LP, and 33 of SI, respectively. Of them, 10 pairs of homologous QTLs were detected between A and D sub-genomes. We also found that 15 co-located QTLs with more than two traits and 12 high-confidence QTLs were detected under more than six environments, respectively. Further, two NET genes (GH_A08G0716 and GH_A08G0783), located in a novel QTL hotspot (qtl24, qtl25 and qlt26) were predominately expressed in early fiber development stages, exhibited significant correlation with LP and SI. The GH_A07G1389 in the stable qtl19 region encoded a tetratricopeptide repeat (TPR)-like superfamily protein and was a homologous gene involved in short fiber mutant ligon lintless-y (Li_y), implying important roles in cotton yield.

CONCLUSIONS

The present study provides a foundation for understanding the regulatory mechanisms of yield components and may enhance yield improvement through molecular breeding in cotton.

Novel Structural Variation and Evolutionary Characteristics of Chloroplast tRNA in Gossypium Plants

Cotton is one of the most important fiber and oil crops in the world. Chloroplast genomes harbor their own genetic materials and are considered to be highly conserved. Transfer RNAs (tRNAs) act as "bridges" in protein synthesis by carrying amino acids. Currently, the variation and evolutionary characteristics of tRNAs in the cotton chloroplast genome are poorly understood. Here, we analyzed the structural variation and evolution of chloroplast tRNA (cp tRNA) based on eight diploid and two allotetraploid cotton species. We also investigated the nucleotide evolution of chloroplast genomes in cotton species. We found that cp tRNAs in cotton encoded 36 or 37 tRNAs, and 28 or 29 anti-codon types with lengths ranging from 60 to 93 nucleotides. Cotton chloroplast tRNA sequences possessed specific conservation and, in particular, the Ψ-loop contained the conserved U-U-C-X3-U. The cp tRNAs of Gossypium L. contained introns, and cp tRNA^Ile contained the anti-codon (C-A-U), which was generally the anti-codon of tRNA^Met. The transition and transversion analyses showed that cp tRNAs in cotton species were iso-acceptor specific and had undergone unequal rates of evolution. The intergenic region was more variable than coding regions, and non-synonymous mutations have been fixed in cotton cp genomes. On the other hand, phylogeny analyses indicated that cp tRNAs of cotton were derived from several inferred ancestors with greater gene duplications. This study provides new insights into the structural variation and evolution of chloroplast tRNAs in cotton plants. Our findings could contribute to understanding the detailed characteristics and evolutionary variation of the tRNA family.