Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Strigolactones act downstream of gibberellins to regulate fiber cell elongation and cell wall thickness in cotton (Gossypium hirsutum)

Strigolactones (SLs) are a class of phytohormones that regulate plant shoot branching and adventitious root development. However, little is known regarding the role of SLs in controlling the behavior of the smallest unit of the organism, the single cell. Here, taking advantage of a classic single-cell model offered by the cotton (Gossypium hirsutum) fiber cell, we show that SLs, whose biosynthesis is fine-tuned by gibberellins (GAs), positively regulate cell elongation and cell wall thickness by promoting the biosynthesis of very long-chain fatty acids (VLCFAs) and cellulose, respectively. Furthermore, we identified two layers of transcription factors (TFs) involved in the hierarchical regulation of this GA-SL crosstalk. The top-layer TF GROWTH-REGULATING FACTOR 4 (GhGRF4) directly activates expression of the SL biosynthetic gene DWARF27 (D27) to increase SL accumulation in fiber cells and GAs induce GhGRF4 expression. SLs induce the expression of four second-layer TF genes (GhNAC100-2, GhBLH51, GhGT2, and GhB9SHZ1), which transmit SL signals downstream to two ketoacyl-CoA synthase genes (KCS) and three cellulose synthase (CesA) genes by directly activating their transcription. Finally, the KCS and CesA enzymes catalyze the biosynthesis of VLCFAs and cellulose, respectively, to regulate development of high-grade cotton fibers. In addition to providing a theoretical basis for cotton fiber improvement, our results shed light on SL signaling in plant development at the single-cell level.

Identification and Functional Analysis of bZIP Genes in Cotton Response to Drought Stress

The basic leucine zipper (bZIP) transcription factors, which harbor a conserved bZIP domain composed of two regions, a DNA-binding basic region and a Leu Zipper region, operate as important switches of transcription networks in eukaryotes. However, this gene family has not been systematically characterized in cotton (Gossypium hirsutum). Here, we identified 197 bZIP family members in cotton. The chromosome distribution pattern indicates that the GhbZIP genes have undergone 53 genome-wide segmental and 7 tandem duplication events which contribute to the expansion of the cotton bZIP family. Phylogenetic analysis showed that cotton GhbZIP proteins cluster into 13 subfamilies, and homologous protein pairs showed similar characteristics. Inspection of the DNA-binding basic region and leucine repeat heptads within the bZIP domains indicated different DNA-binding site specificities as well as dimerization properties among different groups. Comprehensive expression analysis indicated the most highly and differentially expressed genes in root and leaf that might play significant roles in cotton response to drought stress. GhABF3D was identified as a highly and differentially expressed bZIP family gene in cotton leaf and root under drought stress treatments that likely controls drought stress responses in cotton. These data provide useful information for further functional analysis of the GhbZIP gene family and its potential application in crop improvement.

Systematic analysis of Histidine photosphoto transfer gene family in cotton and functional characterization in response to salt and around tolerance

BACKGROUND

Phosphorylation regulated by the two-component system (TCS) is a very important approach signal transduction in most of living organisms. Histidine phosphotransfer (HP) is one of the important members of the TCS system. Members of the HP gene family have implications in plant stresses tolerance and have been deeply studied in several crops. However, upland cotton is still lacking with complete systematic examination of the HP gene family.

RESULTS

A total of 103 HP gene family members were identified. Multiple sequence alignment and phylogeny of HPs distributed them into 7 clades that contain the highly conserved amino acid residue "XHQXKGSSXS", similar to the Arabidopsis HP protein. Gene duplication relationship showed the expansion of HP gene family being subjected with whole-genome duplication (WGD) in cotton. Varying expression profiles of HPs illustrates their multiple roles under altering environments particularly the abiotic stresses. Analysis is of transcriptome data signifies the important roles played by HP genes against abiotic stresses. Moreover, protein regulatory network analysis and VIGS mediated functional approaches of two HP genes (GhHP23 and GhHP27) supports their predictor roles in salt and drought stress tolerance.

CONCLUSIONS

This study provides new bases for systematic examination of HP genes in upland cotton, which formulated the genetic makeup for their future survey and examination of their potential use in cotton production.

Identification, Classification and Characterization Analysis of FBXL Gene in Cotton

F-box/LR (FBXL), Leucine-rich repeats in F-box proteins, belongs to the Skp1-Cullin1-F-box protein (SCF) E3 ligase family. FBXL genes play important roles in plant growth, such as plant hormones, responses to environmental stress, and floral organ development. Here, a total of 518 FBXL genes were identified and analyzed in six plant species. Phylogenetic analysis showed that AtFBXLs, VvFBXLs, and GrFBXLs were clustered into three subfamilies (Ⅰ-Ⅲ). Based on the composition of the F-box domain and carboxyl-terminal amino acid sequence, FBXL proteins were classified into three types (Type-A/-B/-C). Whole-genome duplication (WGD) along with tandem duplications and segmental contributed to the expansion of this gene family. The result indicates that four cotton species are also divided into three subfamilies. FBXLs in cotton were classified into three clades by phylogenetic and structural analyses. Furthermore, expression analyses indicated that the expression patterns of GhFBXLs in different cotton tissues were different. The highly expressed of GH_A07G2363 in 5-8 mm anthers, indicates that this gene might play a role in the reproductive process, providing candidate genes for future studies on cotton fertility materials. This study provides an original functional opinion and a useful interpretation of the FBXL protein family in 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.

Systematical Characterization of the Cotton Di19 Gene Family and the Role of GhDi19-3 and GhDi19-4 as Two Negative Regulators in Response to Salt Stress

Drought-induced 19 (Di19) protein is a Cys2/His2 (C2H2) type zinc-finger protein, which plays a crucial role in plant development and in response to abiotic stress. This study systematically investigated the characteristics of the GhDi19 gene family, including the member number, gene structure, chromosomal distribution, promoter cis-elements, and expression profiles. Transcriptomic analysis indicated that some GhDi19s were up-regulated under heat and salt stress. Particularly, two nuclear localized proteins, GhDi19-3 and GhDi19-4, were identified as being in potential salt stress responsive roles. GhDi19-3 and GhDi19-4 decreased sensitivity under salt stress through virus-induced gene silencing (VIGS), and showed significantly lower levels of H2O2, malondialdehyde (MDA), and peroxidase (POD) as well as significantly increased superoxide dismutase (SOD) activity. This suggested that their abilities were improved to effectively reduce the reactive oxygen species (ROS) damage. Furthermore, certain calcium signaling and abscisic acid (ABA)-responsive gene expression levels showed up- and down-regulation changes in target gene-silenced plants, suggesting that GhDi19-3 and GhDi19-4 were involved in calcium signaling and ABA signaling pathways in response to salt stress. In conclusion, GhDi19-3 and GhDi19-4, two negative transcription factors, were found to be responsive to salt stress through calcium signaling and ABA signaling pathways.

A Comprehensive Analysis of the DUF4228 Gene Family in Gossypium Reveals the Role of GhDUF4228-67 in Salt Tolerance

Soil salinization conditions seriously restrict cotton yield and quality. Related studies have shown that the DUF4228 proteins are pivotal in plant resistance to abiotic stress. However, there has been no systematic identification and analysis of the DUF4228 gene family in cotton and their role in abiotic stress. In this study, a total of 308 DUF4228 genes were identified in four Gossypium species, which were divided into five subfamilies. Gene structure and protein motifs analysis showed that the GhDUF4228 proteins were conserved in each subfamily. In addition, whole genome duplication (WGD) events and allopolyploidization might play an essential role in the expansion of the DUF4228 genes. Besides, many stress-responsive (MYB, MYC) and hormone-responsive (ABA, MeJA) related cis-elements were detected in the promoters of the DUF4228 genes. The qRT-PCR results showed that GhDUF4228 genes might be involved in the response to abiotic stress. VIGS assays and the measurement of relative water content (RWC), Proline content, POD activity, and malondialdehyde (MDA) content indicated that GhDUF4228-67 might be a positive regulator of cotton response to salt stress. The results in this study systematically characterized the DUF4228s in Gossypium species and will provide helpful information to further research the role of DUF4228s in salt tolerance.

Identification of elite fiber quality loci in upland cotton based on the genotyping-by-target-sequencing technology

Genome-wide association studies (GWAS) of fiber quality traits of upland cotton were conducted to identify the single-nucleotide polymorphic (SNP) loci associated with cotton fiber quality, which lays the foundation for the mining of elite] cotton fiber gene resources and its application in molecular breeding. A total of 612 upland cotton accessions were genotyped using the ZJU Cotton Chip No. 1 40K chip array via the liquid-phase probe hybridization-based genotyping-by-target-sequencing (GBTS) technology. In the present study, five fiber quality traits, namely fiber length, fiber strength, micronaire, uniformity and elongation, showed different degrees of variation in different environments. The average coefficient of variation of fiber strength was the greatest, whereas the average coefficient of variation of uniformity was the least. Significant or extremely significant correlations existed among the five fiber quality traits, especially fiber length, strength, uniformity and elongation all being significantly negative correlated with micronaire. Population cluster analysis divided the 612 accessions into four groups: 73 assigned to group I, 226 to group II, 220 to group III and 93 to group IV. Genome-wide association studies of five fiber quality traits in five environments was performed and a total of 42 SNP loci associated with target traits was detected, distributed on 19 chromosomes, with eight loci associated with fiber length, five loci associated with fiber strength, four loci associated with micronaire, twelve loci associated with fiber uniformity and thirteen loci associated with fiber elongation. Of them, seven loci were detected in more than two environments. Nine SNP loci related to fiber length, fiber strength, uniformity and elongation were found on chromosome A07, seven loci related to fiber length, fiber strength, micronaire and elongation were detected on chromosome D01, and five loci associated with fiber length, uniformity and micronaire were detected on chromosome D11. The results from this study could provide more precise molecular markers and genetic resources for cotton breeding for better fiber quality in the future.

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

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

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Status and prospects of genome-wide association studies in cotton

Over the last two decades, the use of high-density SNP arrays and DNA sequencing have allowed scientists to uncover the majority of the genotypic space for various crops, including cotton. Genome-wide association study (GWAS) links the dots between a phenotype and its underlying genetics across the genomes of populations. It was first developed and applied in the field of human disease genetics. Many areas of crop research have incorporated GWAS in plants and considerable literature has been published in the recent decade. Here we will provide a comprehensive review of GWAS studies in cotton crop, which includes case studies on biotic resistance, abiotic tolerance, fiber yield and quality traits, current status, prospects, bottlenecks of GWAS and finally, thought-provoking question. This review will serve as a catalog of GWAS in cotton and suggest new frontiers of the cotton crop to be studied with this important tool.

Identification and functional analysis of PIN family genes in Gossypium barbadense

Background

PIN proteins are an important class of auxin polar transport proteins that play an important regulatory role in plant growth and development. However, their characteristics and functions have not been identified in Gossypium barbadense.

Methods

PIN family genes were identified in the cotton species G. barbadense, Gossypium hirsutum, Gossypium raimondii, and Gossypium arboreum, and detailed bioinformatics analyses were conducted to explore the roles of these genes in G. barbadense using transcriptome data and quantitative reverse-transcription polymerase chain reaction (qRT-PCR) technology. Functional verification of the genes was performed using virus-induced gene silencing (VIGS) technology.

Results

A total of 138 PIN family genes were identified in the four cotton species; the genes were divided into seven subgroups. GbPIN gene family members were widely distributed on 20 different chromosomes, and most had repeated duplication events. Transcriptome analysis showed that some genes had differential expression patterns in different stages of fiber development. According to 'PimaS-7' and '5917' transcript component association analysis, the transcription of five genes was directly related to endogenous auxin content in cotton fibers. qRT-PCR analysis showed that the GbPIN7 gene was routinely expressed during fiber development, and there were significant differences among materials. Transient silencing of the GbPIN7 gene by VIGS led to significantly higher cotton plant growth rates and significantly lower endogenous auxin content in leaves and stems. This study provides comprehensive analyses of the roles of PIN family genes in G. barbadense and their expression during cotton fiber development. Our results will form a basis for further PIN auxin transporter research.

CottonMD: a multi-omics database for cotton biological study

Cotton is an important economic crop, and many loci for important traits have been identified, but it remains challenging and time-consuming to identify candidate or causal genes/variants and clarify their roles in phenotype formation and regulation. Here, we first collected and integrated the multi-omics datasets including 25 genomes, transcriptomes in 76 tissue samples, epigenome data of five species and metabolome data of 768 metabolites from four tissues, and genetic variation, trait and transcriptome datasets from 4180 cotton accessions. Then, a cotton multi-omics database (CottonMD, http://yanglab.hzau.edu.cn/CottonMD/) was constructed. In CottonMD, multiple statistical methods were applied to identify the associations between variations and phenotypes, and many easy-to-use analysis tools were provided to help researchers quickly acquire the related omics information and perform multi-omics data analysis. Two case studies demonstrated the power of CottonMD for identifying and analyzing the candidate genes, as well as the great potential of integrating multi-omics data for cotton genetic breeding and functional genomics research.

Sequence Characteristics and Expression Analysis of GhCIPK23 Gene in Upland Cotton (Gossypium hirsutum L.)

CIPK (calcineurin B-like-interacting protein kinase) is a kind of serine/threonine protein kinase widely existing in plants, and it plays an important role in plant growth and development and stress response. To better understand the biological functions of the GhCIPK23 gene in upland cotton, the coding sequence (CDS) of the GhCIPK23 gene was cloned in upland cotton, and its protein sequence, evolutionary relationship, subcellular localization, expression pattern and cis-acting elements in the promoter region were analyzed. Our results showed that the full-length CDS of GhCIPK23 was 1368 bp, encoding a protein with 455 amino acids. The molecular weight and isoelectric point of this protein were 50.83 KDa and 8.94, respectively. The GhCIPK23 protein contained a conserved N-terminal protein kinase domain and C-terminal regulatory domain of the CIPK gene family member. Phylogenetic tree analysis demonstrated that GhCIPK23 had a close relationship with AtCIPK23, followed by OsCIPK23, and belonged to Group A with AtCIPK23 and OsCIPK23. The subcellular localization experiment indicated that GhCIPK23 was located in the plasma membrane. Tissue expression analysis showed that GhCIPK23 had the highest expression in petals, followed by sepals, and the lowest in fibers. Stress expression analysis showed that the expression of the GhCIPK23 gene was in response to drought, salt, low-temperature and exogenous abscisic acid (ABA) treatment, and had different expression patterns under different stress conditions. Further cis-acting elements analysis showed that the GhCIPK23 promoter region had cis-acting elements in response to abiotic stress, phytohormones and light. These results established a foundation for understanding the function of GhCIPK23 and breeding varieties with high-stress tolerance in cotton.

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

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

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

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Natural variation of GhSI7 increases seed index in cotton

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

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons

Wild relatives of domesticated plants provide a rich resource for crop improvement and a valuable comparative perspective for understanding genomic, physiological, and agricultural traits. Here, we provide high-quality reference genomes of one early domesticated form of the economically most important cotton species, Gossypium hirsutum, and two other wild species, to clarify evolutionary relationships and understand the genomic changes that characterize these species and their close relatives. We document abundant gene resources involved in adaptation to environmental challenges, highlighting the potential for introgression of favorable genes into domesticated cotton and for increasing resilience to climate variability. Our study complements other recent genomic analyses in the cotton genus and provides a valuable foundation for breeding improved cotton varieties.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)₆, Ge] and Gossypium stephensii [(AD)₇, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)₁, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates—including phenotypic differentiation, genetic isolation, and genetic convergence—that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

ATP-citrate lyase B (ACLB) negatively affects cell death and resistance to Verticillium wilt

BACKGROUND

ATP-citrate lyase (ACL) plays a pivotal role in histone acetylation and aerobic glycolysis. In plant, ACL is a heteromeric enzyme composed of ACLA (45 kD) and ACLB (65 kD). So far, the function of ACL genes in cotton still remains unknown.

RESULTS

Here, we identified three ACLA homologous sequences and two ACLB homologous in each genome/sub-genome of cotton species. Silencing ACLB in cotton led to cell death at newly-grown leaves and stem apexes. Simultaneously, in ACLB-silenced plants, transcription factors related to senescence including SGR, WRKY23 and Osl57 were observed to be activated. Further investigation showed that excessive H₂O₂ was accumulated, salicylic acid-dependent defense response and pathogenesis-related gene expressions were evidently enhanced in ACLB-silenced plants, implying that knockdown of ACLB genes leads to hypersensitive response-like cell death in cotton seedlings. However, as noted, serious cell death happened in newly-grown leaves and stem apexes in ACLB-silenced plants, which led to the failure of subsequent fungal pathogenicity assays. To confirm the role of ACLB gene in regulating plant immune response, the dicotyledonous model plant Arabidopsis was selected for functional verification of ACLB gene. Our results indicate the resistance to Verticillium dahliae infection in the Arabidopsis mutant aclb-2 were enhanced without causing strong cell death. Ectopic expression of GausACLB-2 in Arabidopsis weakened its resistance to V. dahliae either in Col-0 or in aclb-2 background, in which the expression level of ACLB is negatively correlated with the resistance to V. dahliae.

CONCLUSIONS

These results indicate that ACLB has a new function in negatively affecting the induction of plant defense response and cell death in cotton, which provides theoretical guidance for developing cotton varieties with resistance against Verticillium wilt.

Deficiency of a peroxisomal NADP-isocitrate dehydrogenase leads to dwarf plant and defect seed in upland cotton

The NADP-isocitrate dehydrogenase-encoded gene GH_D13G1452 with a C-terminus tripeptide Proline-Lysine-Leucine was localized in the peroxisome. It was highly expressed in stems and ovules of 15 days post-anthesis and responded to multiple external stimuli in upland cotton. An upland cotton mutant (Ghpericdh) was identified by flanking sequence amplification and genome variation detection that exogenous sequence was inserted in the middle of the 12th intron of GH_D13G1452, resulting in the deficiency of gene expression. The Ghpericdh mutant displayed a dwarf plant phenotype when grown under field or greenhouse conditions, and GH_D13G1452 functioned as an incomplete dominance on plant height. The germination rate of mutant seed from greenhouse-grown plants was dramatically lower than that from field-grown plants, which indicated that GhperICDH plays a critical role in seed maturation and germination. Therefore, GH_D13G1452 is indispensable in the development of stems and seeds and functions in the adaptability of cotton to the environment. The Ghpericdh mutant provides insight into the function of peroxisomal ICDH and may contribute to the genetic improvement in cotton.

Genome-wide identification of GAD family genes suggests GhGAD6 functionally respond to Cd2+ stress in cotton

Glutamate decarboxylase (GAD) mainly regulated the biosynthesis of γ-aminobutyric acid (GABA) and played an important role in plant growth and stress resistance. To explore the potential function of GAD in cotton growth, the genome-wide identification, structure, and expression analysis of GAD genes were performed in this study. There were 10, 9, 5, and 5 GAD genes identified in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. GAD was divided into four clades according to the protein motif composition, gene structure, and phylogenetic relationship. The segmental duplication was the main way of the GAD gene family evolution. Most GhGADs respond to abiotic stress. Clade Ⅲ GAD was induced by Cd²⁺ stress, especially GhGAD6, and silencing GhGAD6 would lead to more serious Cd²⁺ poisoning in cotton. The oxidative damage caused by Cd²⁺ stress was relieved by increasing the GABA content. It was speculated that the decreased expression of GhGAD6 reduced the content of GABA in vivo and caused the accumulation of ROS. This study will further expand our understanding of the relationship between the evolution and function of the GhGAD gene family and provide new genetic resources for cotton breeding under environmental stress and phytoremediation.

Genome-wide identification, evolution and function analysis of UGTs superfamily in cotton

Glycosyltransferases mainly catalyse the glycosylation reaction in living organisms and widely exists in plants. UGTs have been identified from G. raimondii, G. arboreum and G. hirsutum. However, Genome-wide systematic analysis of UGTs superfamily have not been studied in G. barbadense. 752 UGTs were identified from four cotton species and grouped into 18 clades, of which R was newly discovered clades. Most UGTs were clustered at both ends of the chromosome and showed a heterogeneous distribution. UGT proteins were widely distributed in cells, with the highest distribution in chloroplasts. UGTs of the same clade shared similar intron/exon structural features. During evolution, the gene family has undergone strong selection for purification. UGTs were significantly enriched in “transcriptional activity (GO:0016758)” and “metabolic processes (GO:0008152)”. Genes from the same clade differed in function under various abiotic stresses. The analysis of cis-acting element and qRT–PCR may indicate that GHUGTs play important roles in plant growth, development and abiotic stress. We further found that GHUGT74-2 plays an important role under submergence. The study broadens the understanding of UGTs in terms of gene characteristics, evolutionary processes, and gene function in cotton and provides a new way to systematically and globally understand the structure–function relationship of multigene families in the evolutionary process.

COTTONOMICS: a comprehensive cotton multi-omics database

The rapid advancement of sequencing technology, including next-generation sequencing (NGS), has greatly improved sequencing efficiency and decreased cost. Consequently, huge amounts of genomic, transcriptomic and epigenetic data concerning cotton species have been generated and released. These large-scale data provide immense opportunities for the study of cotton genomic structure and evolution, population genetic diversity and genome-wide mining of excellent genes for important traits. However, the complexity of NGS data also causes distress, as it cannot be utilized easily. Here, we presented the cotton omics data platform COTTONOMICS (http://cotton.zju.edu.cn/), an easily accessible web database that integrates 32.5 TB of omics data including seven assembled genomes, resequencing data from 1180 allotetraploid cotton accessions and RNA-sequencing (RNA-seq), small RNA-sequencing (smRNA-seq), Chromatin Immunoprecipitation sequencing (ChIP-seq), DNase hypersensitive sites sequencing (DNase-seq) and Bisulfite sequencing (BS-seq). COTTONOMICS allows users to employ various search scenarios and retrieve information concerning the cotton genomes, genomic variation (Single nucleotide polymorphisms (SNPs) and Insertion and Deletion (InDels)), gene expression, smRNA expression, epigenetic regulation and quantitative trait locus (QTLs). The user-friendly web interface offers a variety of modules for storing, retrieving, analyzing and visualizing cotton multi-omics data to diverse ends, thereby enabling users to decipher cotton population genetics and identify potential novel genes that influence agronomically beneficial traits.

Database URL: http://cotton.zju.edu.cn

Genome sequence of Gossypium anomalum facilitates interspecific introgression breeding

Crop wild relatives are an important reservoir of natural biodiversity. However, incorporating wild genetic diversity into breeding programs is often hampered by reproductive barriers and a lack of accurate genomic information. We assembled a high-quality, accurately centromere-anchored genome of Gossypium anomalum, a stress-tolerant wild cotton species. We provided a strategy to discover and transfer agronomically valuable genes from wild diploid species to tetraploid cotton cultivars. With a (Gossypium hirsutum × G. anomalum)2 hexaploid as a bridge parent, we developed a set of 74 diploid chromosome segment substitution lines (CSSLs) of the wild cotton species G. anomalum in the G. hirsutum background. This set of CSSLs included 70 homozygous substitutions and four heterozygous substitutions, and it collectively contained about 72.22% of the G. anomalum genome. Twenty-four quantitative trait loci associated with plant height, yield, and fiber qualities were detected on 15 substitution segments. Integrating the reference genome with agronomic trait evaluation of the CSSLs enabled location and cloning of two G. anomalum genes that encode peroxiredoxin and putative callose synthase 8, respectively, conferring drought tolerance and improving fiber strength. We have demonstrated the power of a high-quality wild-species reference genome for identifying agronomically valuable alleles to facilitate interspecific introgression breeding in crops.

Evaluation of genetic gains of some quantitative characters in Egyptian cotton cross (Giza 86 × Menoufi) under water deficit stress

This work was carried out to select cotton genotypes adapted to semi-arid climate conditions cultivated under irrigation for high yields and the standards of the fiber quality properties required by the textile industry. Also to determine the predicted and realized gains from different selection indices to improve some economic characters under water stress conditions. Except for lint percentage and Pressley index, F4 generation reduced PCV and GCV values for all studied characters due to reduction in genetic variability and heterozygosity due to different selection procedures that exhausted a significant part of variability. Except for fiber length and micronaire reading, mean performance in the F4 generation was revealed to be higher than those in the F3 generation for all studied characters. However, micronaire reading was lower (desirable) in F4 than F3 generation. Generally, genotypic correlations were higher than phenotypic correlations. Direct selection for lint index (Ped.3) was the most efficient in improving lint cotton yield/plant and bolls/plant. However, the multiplicative index involving all studied characters (I.5) exhibited the highest values for boll weight. Also, the Ped.2 index (direct selection for lint percentage) proved to be the most efficient in improving seed and lint indexes. Direct selection for lint cotton yield/plant (Ped.1) could produce the highest desirable values for lint percentage and seed per boll with a relatively reasonable yield. A selection index involving yield and its components (I.3) is recommended in improving uniformity index, fiber strength, and micronaire reading. The superior five families released from these indices in F4 generation exceeded the better parent for lint cotton yield/plant, bolls/plant, boll weight, seeds/boll, lint index, and reasonable fiber traits. These families could be continued to further generations as breeding material for developing water deficit tolerant genotypes.

Retrieving a disrupted gene encoding phospholipase A for fibre enhancement in allotetraploid cultivated cotton

After polyploidization originated from one interspecific hybridization event in Gossypium, Gossypium barbadense evolved to produce extra-long staple fibres than Gossypium hirsutum (Upland cotton), which produces a higher fibre yield. The genomic diversity between G. barbadense and G. hirsutum thus provides a genetic basis for fibre trait variation. Recently, rapid accumulation of gene disruption or deleterious mutation was reported in allotetraploid cotton genomes, with unknown impacts on fibre traits. Here, we identified gene disruptions in allotetraploid G. hirsutum (18.14%) and G. barbadense (17.38%) through comparison with their presumed diploid progenitors. Relative to conserved genes, these disrupted genes exhibited faster evolution rate, lower expression level and altered gene co-expression networks. Within a module regulating fibre elongation, a hub gene experienced gene disruption in G. hirsutum after polyploidization, with a 2-bp deletion in the coding region of GhNPLA1D introducing early termination of translation. This deletion was observed in all of the 34 G. hirsutum landraces and 36 G. hirsutum cultivars, but not in 96% of 57 G. barbadense accessions. Retrieving the disrupted gene GhNPLA1D using its homoeolog GhNPLA1A achieved longer fibre length in G. hirsutum. Further enzyme activity and lipids analysis confirmed that GhNPLA1A encodes a typical phospholipase A and promotes cotton fibre elongation via elevating intracellular levels of linolenic acid and 34:3 phosphatidylinositol. Our work opens a strategy for identifying disrupted genes and retrieving their functions in ways that can provide valuable resources for accelerating fibre trait enhancement in cotton breeding.

Genome-wide analysis of serine carboxypeptidase-like protein (SCPL) family and functional validation of Gh_SCPL42 unchromosome conferring cotton Verticillium der Verticillium wilt stress in Gossypium hirsutum

BACKGROUND

Serine carboxypeptidase-like protein (SCPL) plays an important role in response to stress in plant. However, our knowledge of the function of the SCPL gene family is limited.

RESULTS

In this study, a comprehensive and systematic analysis of SCPL gene family was conducted to explore the phylogeny and evolution of the SCPL gene in Gossypium hirsutum. The phenotype and molecular mechanism of silencing of the Gh_SCPL42 under Verticillium wilt stress was also studied. Our results showed that 96 SCPL genes were observed in genome of G. hirsutum, which distributed on 25 chromosomes and most of them were located in the nucleus. The phylogenetic tree analysis showed that members of SCPL gene family can be divided into three subgroups in G. hirsutum, which are relatively conservative in evolution. SCPL gene has a wide range of tissue expression types in G. hirsutum. Promoter analysis showed that the most cis-acting elements related to MeJA and ABA were contained. Through RNA-seq combined with genotyping, it was found that 11 GhSCPL genes not only had significant expression changes during Verticillium wilt stress but also had differential SNPs in the upstream, downstream, exonic or intronic regions. The expression of these 11 genes in the resistant (Zhongzhimian 2) and susceptible (Junmian 1) materials was further analyzed by qRT-PCR, it was found that 6 genes showed significant expression differences in the two materials. Among them, Gh_SCPL42 has the most obvious expression change. Furthermore, virus-induced gene silencing (VIGS) showed necrosis and yellowing of leaves and significantly higher disease severity index (DSI) and disease severity rate (DSR) values in VIGS plants than in control silenced Gh_SCPL42 plants. Moreover, the expression levels of genes related to the SA and JA pathways were significantly downregulated. These results show that Gh_SCPL42 might improve resistance to Verticillium wilt through the SA and JA pathways in G. hirsutum.

CONCLUSION

In conclusion, our findings indicated that Gh_SCPL42 gene plays an important role in resistance to Verticillium wilt in cotton. It was provided an important theoretical basis for further research on the function of SCPL gene family and the molecular mechanism of resistance to Verticillium wilt in cotton.

Functional divergence of GLP genes between G. barbadense and G. hirsutum in response to Verticillium dahliae infection

Germin-like proteins (GLPs) play important roles in plant disease resistance but are rarely reported in cotton. We compared the expression of GLPs in Verticillium dahliae inoculate G. hirsutum (susceptible) and G. barbadense (resistant) and enriched 11 differentially expressed GLPs. 2741 GLP proteins identified from 53 species determined that GLP probably originated from algae and could be classified into 7 clades according to phylogenetic analysis, among which Clade I is likely the most ancient. Cotton GLP (two allopolyploids and two diploids) genes within a shared clade were highly conserved. Intriguingly, clade VII genes were mainly located in gene clusters that derived from the expansion of LTR transposons. Clade VII members expressed mainly in root which is the first battle against Verticillium dahlia and could be induced more intensely in G. barbadense than G. hirsutum. The GLP genes are resistant to Verticillium dahliae, which can be further investigated against Verticillium wilt.

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.

Genome-Wide Identification and Expression Analyses of the Cotton AGO Genes and Their Potential Roles in Fiber Development and Stress Response

Argonaute proteins (AGOs) are indispensable components of RNA silencing. However, systematic characterization of the AGO genes have not been completed in cotton until now. In this study, cotton AGO genes were identified and analyzed with respect to their evolution and expression profile during biotic and abiotic stresses. We identified 14 GaAGO, 14 GrAGO, and 28 GhAGO genes in the genomes of Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum. Cotton AGO proteins were classified into four subgroups. Structural and functional conservation were observed in the same subgroups based on the analysis of the gene structure and conserved domains. Twenty-four duplicated gene pairs were identified in GhAGO genes, and all of them exhibited strong purifying selection during evolution. Moreover, RNA-seq analysis showed that most of the GhAGO genes exhibit high expression levels in the fiber initiation and elongation processes. Furthermore, the expression profiles of GhAGO genes tested by quantitative real-time polymerase chain reaction (qPCR) demonstrated that they were sensitive to Verticillium wilt infection and salt and drought stresses. Overall, our results will pave the way for further functional investigation of the cotton AGO gene family, which may be involved in fiber development and stress response.

Impact of water deficiency on leaf cuticle lipids and gene expression networks in cotton (Gossypium hirsutum L.)

BACKGROUND

Water deficit (WD) has serious effect on the productivity of crops. Formation of cuticular layer with increased content of wax and cutin on leaf surfaces is closely related to drought tolerance. Identification of drought tolerance associated wax components and cutin monomers and the genes responsible for their biosynthesis is essential for understanding the physiological and genetic mechanisms underlying drought tolerance and improving crop drought resistance.

RESULT

In this study, we conducted comparative phenotypic and transcriptomic analyses of two Gossypium hirsutum varieties that are tolerant (XL22) or sensitive (XL17) to drought stress. XL17 consumed more water than XL22, particularly under the WD conditions. WD significantly induced accumulation of most major wax components (C29 and C31 alkanes) and cutin monomers (palmitic acid and stearic acid) in leaves of both XL22 and XL17, although accumulation of the major cutin monomers, i.e., polyunsaturated linolenic acid (C18:3n-3) and linoleic acid (C18:2n-6), were significantly repressed by WD in both XL22 and XL17. According to the results of transcriptome analysis, although many genes and their related pathways were commonly induced or repressed by WD in both XL22 and XL17, WD-induced differentially expressed genes specific to XL22 or XL17 were also evident. Among the genes that were commonly induced by WD were the GhCER1 genes involved in biosynthesis of alkanes, consistent with the observation of enhanced accumulation of alkanes in cotton leaves under the WD conditions. Interestingly, under the WD conditions, several GhCYP86 genes, which encode enzymes catalyzing the omega-hydroxylation of fatty acids and were identified to be the hub genes of one of the co-expression gene modules, showed a different expression pattern between XL22 and XL17 that was in agreement with the WD-induced changes of the content of hydroxyacids or fatty alcohols in these two varieties.

CONCLUSION

The results contribute to our comprehending the physiological and genetic mechanisms underlying drought tolerance and provide possible solutions for the difference of drought resistance of different cotton varieties.

Structure and character analysis of cotton response regulator genes family reveals that GhRR7 responses to draught stress

Background

Cytokinin signal transduction is mediated by a two-component system (TCS). Two-component systems are utilized in plant responses to hormones as well as to biotic and abiotic environmental stimuli. In plants, response regulatory genes (RRs) are one of the main members of the two-component system (TCS).

Method

From the aspects of gene structure, evolution mode, expression type, regulatory network and gene function, the evolution process and role of RR genes in the evolution of the cotton genome were analyzed.

Result

A total of 284 RR genes in four cotton species were identified. Including 1049 orthologous/paralogous gene pairs were identified, most of which were whole genome duplication (WGD). The RR genes promoter elements contain phytohormone responses and abiotic or biotic stress-related cis-elements. Expression analysis showed that RR genes family may be negatively regulate and involved in salt stress and drought stress in plants. Protein regulatory network analysis showed that RR family proteins are involved in regulating the DNA-binding transcription factor activity (COG5641) pathway and HP kinase pathways. VIGS analysis showed that the GhRR7 gene may be in the same regulatory pathway as GhAHP5 and GhPHYB, ultimately negatively regulating cotton drought stress by regulating POD, SOD, CAT, H₂O₂ and other reactive oxygen removal systems.

Conclusion

This study is the first to gain insight into RR gene members in cotton. Our research lays the foundation for discovering the genes related to drought and salt tolerance and creating new cotton germplasm materials for drought and salt tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00394-2.

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

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

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.

Advanced genes expression pattern greatly contributes to divergence in Verticillium wilt resistance between Gossypium barbadense and Gossupium hirsutum

Verticillium, representing one of the world's major pathogens, causes Verticillium wilt in important woody species, ornamentals, agricultural, etc., consequently resulting in a serious decline in production and quality, especially in cotton. Gossupium hirutum and Gossypium barbadense are two kinds of widely cultivated cotton species that suffer from Verticillium wilt, while G. barbadense has much higher resistance toward it than G. hirsutum. However, the molecular mechanism regarding their divergence in Verticillium wilt resistance remains largely unknown. In the current study, G. barbadense cv. Hai7124 and G. hirsutum acc. TM-1 were compared at 0, 12, 24, 48, 72, 96, 120, and 144 h post-inoculation (hpi) utilizing high throughput RNA-Sequencing. As a result, a total of 3,549 and 4,725 differentially expressed genes (DEGs) were identified, respectively. In particular, the resistant type Hai7124 displayed an earlier and faster detection and signaling response to the Verticillium dahliae infection and demonstrated higher expression levels of defense-related genes over TM-1 with respect to transcription factors, plant hormone signal transduction, plant-pathogen interaction, and nucleotide-binding leucine-rich repeat (NLR) genes. This study provides new insights into the molecular mechanisms of divergence in Verticillium wilt resistance between G. barbadense and G. hirsutum and important candidate genes for breeding V. dahliae resistant cotton cultivars.

GhBZR3 suppresses cotton fiber elongation by inhibiting very-long-chain fatty acid biosynthesis

The BRASSINAZOLE-RESISTANT (BZR) transcription factor is a core component of brassinosteroid (BR) signaling and is involved in the development of many plant species. BR is essential for the initiation and elongation of cotton fibers. However, the mechanism of BR-regulating fiber development and the function of BZR is poorly understood in Gossypium hirsutum L. (cotton). Here, we identified a BZR family transcription factor protein referred to as GhBZR3 in cotton. Overexpression of GhBZR3 in Arabidopsis caused shorter root hair length, hypocotyl length, and hypocotyl cell length, indicating that GhBZR3 negatively regulates cell elongation. Pathway enrichment analysis from VIGS-GhBZR3 cotton plants found that fatty acid metabolism and degradation might be the regulatory pathway that is primarily controlled by GhBZR3. Silencing GhBZR3 expression in cotton resulted in taller plant height as well as longer fibers. The very-long-chain fatty acid (VLCFA) content was also significantly increased in silenced GhBZR3 plants compared with the wild type. The GhKCS13 promoter, a key gene for VLCFA biosynthesis, contains two GhBZR3 binding sites. The results of yeast one-hybrid, electrophoretic mobility shift, and luciferase assays revealed that GhBZR3 directly interacted with the GhKCS13 promoter to suppress gene expression. Taken together, these results indicate that GhBZR3 negatively regulates cotton fiber development by reducing VLCFA biosynthesis. This study not only deepens our understanding of GhBZR3 function in cotton fiber development, but also highlights the potential of improving cotton fiber length and plant growth using GhBZR3 and its related genes in future cotton breeding programs.

Gb_ANR-47 Enhances the Resistance of Gossypium barbadense to Fusarium oxysporum f. sp. vasinfectum (FOV) by Regulating the Content of Proanthocyanidins

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

A genome-wide comparative evolutionary analysis of zinc finger-BED transcription factor genes in land plants

Zinc finger (Zf)-BED proteins are a novel superfamily of transcription factors that controls numerous activities in plants including growth, development, and cellular responses to biotic and abiotic stresses. Despite their important roles in gene regulation, little is known about the specific functions of Zf-BEDs in land plants. The current study identified a total of 750 Zf-BED-encoding genes in 35 land plant species including mosses, bryophytes, lycophytes, gymnosperms, and angiosperms. The gene family size was somewhat proportional to genome size. All identified genes were categorized into 22 classes based on their specific domain architectures. Of these, class I (Zf-BED_DUF-domain_Dimer_Tnp_hAT) was the most common in the majority of the land plants. However, some classes were family-specific, while the others were species-specific, demonstrating diversity at different classification levels. In addition, several novel functional domains were also predicated including WRKY and nucleotide-binding site (NBS). Comparative genomics, transcriptomics, and proteomics provided insights into the evolutionary history, duplication, divergence, gene gain and loss, species relationship, expression profiling, and structural diversity of Zf-BEDs in land plants. The comprehensive study of Zf-BEDs in Gossypium sp., (cotton) also demonstrated a clear footprint of polyploidization. Overall, this comprehensive evolutionary study of Zf-BEDs in land plants highlighted significant diversity among plant species.

Equilibrated evolution of the mixed auto-/allopolyploid haplotype-resolved genome of the invasive hexaploid Prussian carp

Understanding genome evolution of polyploids requires dissection of their often highly similar subgenomes and haplotypes. Polyploid animal genome assemblies so far restricted homologous chromosomes to a 'collapsed' representation. Here, we sequenced the genome of the asexual Prussian carp, which is a close relative of the goldfish, and present a haplotype-resolved chromosome-scale assembly of a hexaploid animal. Genome-wide comparisons of the 150 chromosomes with those of two ancestral diploid cyprinids and the allotetraploid goldfish and common carp revealed the genomic structure, phylogeny and genome duplication history of its genome. It consists of 25 syntenic, homeologous chromosome groups and evolved by a recent autoploid addition to an allotetraploid ancestor. We show that de-polyploidization of the alloploid subgenomes on the individual gene level occurred in an equilibrated fashion. Analysis of the highly conserved actinopterygian gene set uncovered a subgenome dominance in duplicate gene loss of one ancestral chromosome set.

Characterization of a Novel Insect-Induced Sesquiterpene Synthase GbTPS1 Based on the Transcriptome of Gossypium barbadense Feeding by Cotton Bollworm

When attacked by insect herbivores, plants initiate sophisticated defenses mediated by complex signaling networks and usually release a blend of functional volatiles such as terpenes against infestation. The extra-long staple cotton Gossypium barbadense cultivated worldwide as natural textile fiber crop is frequently exposed to a variety of herbivores, such as cotton bollworm Helicoverpa armigera. However, little is known about insect-induced transcriptional changes and molecular mechanisms underlying subsequent defense responses in G. barbadense. In the current study, transcriptome changes in G. barbadense infested with chewing H. armigera larvae were investigated, and we identified 5,629 differentially expressed genes (DEGs) in the infested cotton leaves compared with non-infested controls. H. armigera feeding triggered complex signaling networks in which almost all (88 out of 90) DEGs associated with the jasmonic acid (JA) pathway were upregulated, highlighting a central role for JA in the defense responses of G. barbadense against target insects. All DEGs involved in growth-related photosynthesis were downregulated, whereas most DEGs associated with defense-related transcript factors and volatile secondary metabolism were upregulated. It was noteworthy that a terpene synthase gene in the transcriptome data, GbTPS1, was strongly expressed in H. armigera-infested G. barbadense leaves. The upregulation of GbTPS1 in qPCR analysis also suggested an important role for GbTPS1 in herbivore-induced cotton defense. In vitro assays showed that recombinant GbTPS1 catalyzed farnesyl pyrophosphate and neryl diphosphate to produce three sesquiterpenes (selinene, α-gurjunene, and β-elemene) and one monoterpene (limonene), respectively. Moreover, these catalytic products of GbTPS1 were significantly elevated in G. barbadense leaves after H. armigera infestation, and elemene and limonene had repellent effects on H. armigera larvae in a dual-choice bioassay and increased larval mortality in a no-choice bioassay. These findings provide a valuable insight into understanding the transcriptional changes reprogramming herbivore-induced sesquiterpene biosynthesis in G. barbadense infested by H. armigera, which help elucidate the molecular mechanisms underlying plant defense against insect pests.

Genome-wide association study reveals novel quantitative trait loci and candidate genes of lint percentage in upland cotton based on the CottonSNP80K array

Thirty-four SNPs corresponding with 22 QTLs for lint percentage, including 13 novel QTLs, was detected via GWAS. Two candidate genes underlying this trait were also identified. Cotton (Gossypium spp.) is an important natural textile fiber and oilseed crop cultivated worldwide. Lint percentage (LP, %) is one of the important yield components, and increasing LP is a core goal of cotton breeding improvement. However, the genetic and molecular mechanisms underlying LP in upland cotton remain unclear. Here, we performed a genome-wide association study (GWAS) for LP based on 254 upland cotton accessions in four environments as well as the best linear unbiased predictors using the high-density CottonSNP80K array. In total, 41,413 high-quality single-nucleotide polymorphisms (SNPs) were screened, and 34 SNPs within 22 quantitative trait loci (QTLs) were significantly associated with LP. In total, 175 candidate genes were identified from two major genomic loci (GR1 and GR2), and 50 hub genes were identified through GO enrichment and weighted gene co-expression network analysis. Two candidate genes (Gh_D01G0162 and Gh_D07G0463), which may participate in early fiber development to affect the number of fiber protrusions and LP, were also identified. Their genetic variation and expression were verified by linkage disequilibrium blocks, haplotypes, and quantitative real-time polymerase chain reaction, respectively. The weighted gene interaction network analysis showed that the expression of Gh_D07G0463 was significantly correlated with that of Gh_D01G0162. These identified SNPs, QTLs and candidate genes provide important insights into the genetic and molecular mechanisms underlying variations in LP and serve as a foundation for LP improvement via marker-assisted breeding.

Expression Patterns Divergence of Reciprocal F1 Hybrids Between Gossypium hirsutum and Gossypium barbadense Reveals Overdominance Mediating Interspecific Biomass Heterosis

Hybrid breeding has provided an impetus to the process and achievement of a higher yield and quality of crops. Interspecific hybridization is critical for resolving parental genetic diversity bottleneck problems. The reciprocal interspecific hybrids and their parents (Gossypium hirsutum and Gossypium barbadense) have been applied in this study to elucidate the transcription regulatory mechanism of early biomass heterosis. Phenotypically, the seed biomass, plant height over parent heterosis, leaf area over parent heterosis, and fresh and dry biomass were found to be significantly higher in hybrids than in parents. Analysis of leaf areas revealed that the one-leaf stage exhibits the most significant performance in initial vegetative growth vigor and larger leaves in hybrids, increasing the synthesis of photosynthesis compounds and enhancing photosynthesis compound synthesis. Comparative transcriptome analysis showed that transgressive down-regulation (TDR) is the main gene expression pattern in the hybrids (G. hirsutum × G. barbadense, HB), and it was found that the genes of photosystem I and Adenosine triphosphate (ATP)-binding may promote early growth vigor. Transgressive up-regulation (TUR) is the major primary gene expression pattern in the hybrids (G. barbadense × G. hirsutum, BH), and photosystem II-related genes mediated the performance of early biomass heterosis. The above results demonstrated that overdominance mediates biomass heterosis in interspecific hybrid cotton and the supervisory mechanism divergence of hybrids with different females. Photosynthesis and other metabolic process are jointly involved in controlling early biomass heterosis in interspecific hybrid cotton. The expression pattern data of transcriptome sequencing were supported using the qRT-PCR analysis. Our findings could be useful in theoretical and practical studies of early interspecific biomass heterosis, and the results provide potential resources for the theoretical and applied research on early interspecific biomass heterosis.

A GWAS identified a major QTL for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a MAGIC population of Upland cotton and a meta-analysis of QTLs for Fusarium wilt resistance

A major QTL conferring resistance to Fusarium wilt race 4 in a narrow region of chromosome D02 was identified in a MAGIC population of 550 RILs of Upland cotton. Numerous studies have been conducted to investigate the genetic basis of Fusarium wilt (FW, caused by Fusarium oxysporum f. sp. vasinfectum, FOV) resistance using bi-parental and association mapping populations in cotton. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs), together with their 11 Upland cotton (Gossypium hirsutum) parents, was used to identify QTLs for FOV race 4 (FOV4) resistance. Among the parents, Acala Ultima, M-240 RNR, and Stoneville 474 were the most resistant, while Deltapine Acala 90, Coker 315, and Stoneville 825 were the most susceptible. Twenty-two MAGIC lines were consistently resistant to FOV4. Through a genome-wide association study (GWAS) based on 473,516 polymorphic SNPs, a major FOV4 resistance QTL within a narrow region on chromosomes D02 was detected, allowing identification of 14 candidate genes. Additionally, a meta-analysis of 133 published FW resistance QTLs showed a D subgenome and individual chromosome bias and no correlation between homeologous chromosome pairs. This study represents the first GWAS study using a largest genetic population and the most comprehensive meta-analysis for FW resistance in cotton. The results illustrated that 550 lines were not enough for high resolution mapping to pinpoint a candidate gene, and experimental errors in phenotyping cotton for FW resistance further compromised the accuracy and precision in QTL localization and identification of candidate genes. This study identified FOV4-resistant parents and MAGIC lines, and the first major QTL for FOV4 resistance in Upland cotton, providing useful information for developing FOV4-resistant cultivars and further genomic studies towards identification of causal genes for FOV4 resistance in cotton.

Consequences of Drought Stress Encountered During Seedling Stage on Physiology and Yield of Cultivated Cotton

Survival of living organisms depends on the availability of water resources required for agriculture. In the current scenario of limited water resources, it is our priority to maximise the yield potential of crops with a minimum supply of available water. In this study, we evaluated seven cultivated varieties of Gossypium hirsutum (FH-114, FH-152, FH-326, FH-492, FH-942, VH-327 and FH-NOOR) for their tolerance, yield potential and fibre quality under water shortages. We also studied the effect of drought stress on osmoregulation, chlorophyll content, antioxidant (peroxidase and catalase) activity, lipid peroxidation and secondary metabolite accumulation in the varieties under study. It was revealed that three varieties (FH-114, FH-152 and VH-327) exhibited a lower stress susceptibility index and more tolerance to drought stress. All the varieties demonstrated enhanced proline and malondialdehyde content, but no significant change in chlorophyll content was observed under limited water supply. Antioxidant activity offered by catalase and phenolic content was enhanced in FH-492 whilst peroxidase activity increased in FH-114 and FH-326. Phenolic content was highest in FH-942 and decreased significantly in the remaining varieties. Ginning outturn of the cotton varieties increased in VH-327 (19.8%) and FH-326 (3.7%), was not affected in FH-114 and FH-492 and was reduced in FH-152, FH-942 and FH-NOOR. All cotton varieties tested showed an increase in micronaire thickness when exposed to drought stress as early as the seedling stage. This study highlights the evaluation and screening of cotton varieties for their response to drought stress in terms of yield and fibre quality when exposed to water shortages during plant development and can help in devising irrigation plans.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

The Current Progresses in the Genes and Networks Regulating Cotton Plant Architecture

Cotton is the most important source of natural fiber in the world as well as a key source of edible oil. The plant architecture and flowering time in cotton are crucial factors affecting cotton yield and the efficiency of mechanized harvest. In the model plant arabidopsis, the functions of genes related to plant height, inflorescence structure, and flowering time have been well studied. In the model crops, such as tomato and rice, the similar genetic explorations have greatly strengthened the economic benefits of these crops. Plants of the Gossypium genus have the characteristics of perennials with indeterminate growth and the cultivated allotetraploid cottons, G. hirsutum (Upland cotton), and G. barbadense (Sea-island cotton), have complex branching patterns. In this paper, we review the current progresses in the identification of genes affecting cotton architecture and flowering time in the cotton genome and the elucidation of their functional mechanisms associated with branching patterns, branching angle, fruit branch length, and plant height. This review focuses on the following aspects: (i) plant hormone signal transduction pathway; (ii) identification of cotton plant architecture QTLs and PEBP gene family members; (iii) functions of FT/SFT and SP genes; (iv) florigen and anti-florigen systems. We highlight areas that require further research, and should lay the groundwork for the targeted bioengineering of improved cotton cultivars with flowering times, plant architecture, growth habits and yields better suited for modern, mechanized cultivation.

Dynamic Quantitative Trait Loci Mapping for Plant Height in Recombinant Inbred Line Population of Upland Cotton

Plant height (PH) is a key plant architecture trait for improving the biological productivity of cotton. Ideal PH of cotton is conducive to lodging resistance and mechanized harvesting. To detect quantitative trait loci (QTL) and candidate genes of PH in cotton, a genetic map was constructed with a recombinant inbred line (RIL) population of upland cotton. PH phenotype data under nine environments and three best linear unbiased predictions (BLUPs) were used for QTL analyses. Based on restriction-site-associated DNA sequence (RAD-seq), the genetic map contained 5,850 single-nucleotide polymorphism (SNP) markers, covering 2,747.12 cM with an average genetic distance of 0.47 cM. Thirty-seven unconditional QTL explaining 1.03-12.50% of phenotypic variance, including four major QTL and seven stable QTL, were identified. Twenty-eight conditional QTL explaining 3.27-28.87% of phenotypic variance, including 1 major QTL, were identified. Importantly, five QTL, including 4 stable QTL, were both unconditional and conditional QTL. Among the 60 PH QTL (including 39 newly identified), none of them were involved in the whole period of PH growth, indicating that QTL related to cotton PH development have dynamic expression characteristics. Based on the functional annotation of Arabidopsis homologous genes and transcriptome data of upland cotton TM-1, 14 candidate genes were predicted within 10 QTL. Our research provides valuable information for understanding the genetic mechanism of PH development, which also increases the economic production of cotton.

Potential Roles of 1-Aminocyclopropane-1-carboxylic Acid Synthase Genes in the Response of Gossypium Species to Abiotic Stress by Genome-Wide Identification and Expression Analysis

Ethylene plays a pivotal role in plant stress resistance and 1-aminocyclopropane-1-carboxylic acid synthase (ACS) is the rate-limiting enzyme in ethylene biosynthesis. Upland cotton (Gossypium hirsutum L.) is the most important natural fiber crop, but the function of ACS in response to abiotic stress has rarely been reported in this plant. We identified 18 GaACS, 18 GrACS, and 35 GhACS genes in Gossypiumarboreum, Gossypium raimondii and Gossypiumhirsutum, respectively, that were classified as types I, II, III, or IV. Collinearity analysis showed that the GhACS genes were expanded from diploid cotton by the whole-genome-duplication. Multiple alignments showed that the C-terminals of the GhACS proteins were conserved, whereas the N-terminals of GhACS10 and GhACS12 were different from the N-terminals of AtACS10 and AtACS12, probably diverging during evolution. Most type II ACS genes were hardly expressed, whereas GhACS10/GhACS12 were expressed in many tissues and in response to abiotic stress; for example, they were highly and hardly expressed at the early stages of cold and heat exposure, respectively. The GhACS genes showed different expression profiles in response to cold, heat, drought, and salt stress by quantitative PCR analysis, which indicate the potential roles of them when encountering the various adverse conditions, and provide insights into GhACS functions in cotton’s adaptation to abiotic stress.

Characterization and Functional Analysis of GhNAC82, A NAM Domain Gene, Coordinates the Leaf Senescence in Upland Cotton (Gossypium hirsutum L.)

In the process of growth and development, cotton exhibits premature senescence under various abiotic stresses, impairing yield and fiber quality. NAC (NAM, ATAF1,2, and CUC2) protein widely distributed in land plants, play the critical role in responding to abiotic stress and regulating leaf senescence. We have identified and functional analyzed a NAM domain gene GhNAC82 in upland cotton, it was located on the A11 chromosome 4,921,702 to 4,922,748 bp, only containing one exon. The spatio-temporal expression pattern analysis revealed that it was highly expressed in root, torus, ovule and fiber development stage. The results of qRT-PCR validated that GhNAC82 negatively regulated by salt stress, drought stress, H₂O₂ stress, IAA treatment, and ethylene treatment, positively regulated by the ABA and MeJA treatment. Moreover, heterologous overexpression of GhNAC82 results in leaf premature senescence and delays root system development in Arabidopsis thaliana. The phenotype of delayed-senescence was performed after silencing GhNAC82 by VIGS in premature cotton. Taken together, GhNAC82 was involved in different abiotic stress pathways and play important roles in negatively regulating leaf premature senescence.