Gossypium barbadense genome sequence provides insight into the evolution of extra-long staple fiber and specialized metabolites

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.

Analysis of the histidine kinase gene family and the role of GhHK8 in response to drought tolerance in cotton

As an important member of the two-component system (TCS), histidine kinases (HKs) play important roles in various plant developmental processes and signal transduction in response to a wide range of biotic and abiotic stresses. So far, the HK gene family has not been investigated in Gossypium. In this study, a total of 177 HK gene family members were identified in cotton. They were further divided into seven groups, and the protein characteristics, genetic relationship, gene structure, chromosome location, collinearity, and cis-elements identification were comprehensively analyzed. Whole genome duplication (WGD) / segmental duplication may be the reason why the number of HK genes doubled in tetraploid Gossypium species. Expression analysis revealed that most cotton HK genes were mainly expressed in the reproductive organs and the fiber at initial stage. Gene expression analysis revealed that HK family genes are involved in cotton abiotic stress, especially drought stress and salt stress. In addition, gene interaction networks showed that HKs were involved in the regulation of cotton abiotic stress, especially drought stress. VIGS experiments have shown that GhHK8 is a negative regulatory factor in response to drought stress. Our systematic analysis provided insights into the characteristics of the HK genes in cotton and laid a foundation for further exploring their potential in drought stress resistance in 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.

Identification of salt stress-tolerant candidate genes in the BC2F2 population at the seedling stages of G. hirsutum and G. darwinii using NGS-based bulked segregant analysis

Salinity is a major threat to the yield and productivity of cotton seedlings. In the present study, we developed a BC₂F₂ population of cotton plants from Gossypium darwinii (5-7) and Gossypium hirsutum (CCRI 12-4) salt-susceptible parents to identify salt-resistant candidate genes. The Illumina HiSeq™ strategy was used with bulked segregant analysis. Salt-resistant and salt-susceptible DNA bulks were pooled by using 30 plants from a BC₂F₂ population. Next-generation sequencing (NGS) technology was used for the sequencing of parents and both bulks. Four significant genomic regions were identified: the first genomic region was located on chromosome 18 (1.86 Mb), the second and third genomic regions were on chromosome 25 (1.06 Mb and 1.94 Mb, respectively), and the fourth was on chromosome 8 (1.41 Mb). The reads of bulk1 and bulk2 were aligned to the G. darwinii and G. hirsutum genomes, respectively, leading to the identification of 20,664,007 single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels). After the screening, 6,573 polymorphic markers were obtained after filtration of the candidate regions. The SNP indices in resistant and susceptible bulks and Δ(SNP-index) values of resistant and susceptible bulks were measured. Based on the higher Δ(SNP-index) value, six effective polymorphic SNPs were selected in a different chromosome. Six effective SNPs were linked to five candidate genes in four genomic regions. Further validation of these five candidate genes was carried out using reverse transcription-quantitative polymerase chain reaction (RT-qPCR), resulting in an expression profile that showed two highly upregulated genes in the salt-tolerant species G. darwinii, i.e., Gohir.D05G367800 and Gohir.D12G239100; however, the opposite was shown in G. hirsutum, for which all genes, except one, showed partial expression. The results indicated that Gohir.D05G367800 and Gohir.D12G239100 may be salt-tolerant genes. We are confident that this study could be helpful for the cloning, transformation, and development of salt-resistant cotton varieties.

Genome-wide identification of terpenoid synthase family genes in Gossypium hirsutum and functional dissection of its subfamily cadinene synthase A in gossypol synthesis

Plant terpenoid synthase (TPS) family genes participate in metabolite synthesis, hormones, gossypol, etc. Here, we genome-widely identified TPS family genes in 12 land plant species. Four hundred and thirty TPS-related genes were divided into seven subfamilies. The TPS-c in Bryophytes was suggested to be the earliest subfamily, followed by the TPS-e/f and TPS-h presence in ferns. TPS-a, the largest number of genes, was derived from monocotyledonous and dicotyledonous plants. Collinearity analysis showed that 38 out of the 76 TPS genes in G. hirsutum were collinear within G. arboreum and G. raimondii. Twenty-one GhTPS-a genes belong to the cadinene synthase (GhCDN) subfamily and were divided into five groups, A, B, C, D, and E. The special cis-elements in the promoters of 12 GhCDN-A genes suggested that the JA and ethylene signaling pathways may be involved in their expression regulation. When 12 GhCDN-A genes were simultaneously silenced through virus-induced gene silencing, the glandular color of GhCDN-A-silenced plants was lighter than that of the control, supported by a gossypol content decrease based on HPLC testing, suggesting that GhCDN-A subgroup genes participate in gossypol synthesis. According to RNA-seq analysis, gossypol synthesis-related genes and disease-resistant genes in the glandular variety exhibited upregulated expression compared to the glandless variety, whereas hormone signaling-related genes were downregulated. All in all, these results revealed plant TPS gene evolution rules and dissected the TPS subfamily, GhCDN-A, function in gossypol synthesis in cotton.

Structural variation (SV)-based pan-genome and GWAS reveal the impacts of SVs on the speciation and diversification of allotetraploid cottons

Structural variations (SVs) have long been described as being involved in the origin, adaption, and domestication of species. However, the underlying genetic and genomic mechanisms are poorly understood. Here, we report a high-quality genome assembly of Gossypium barbadense acc. Tanguis, a landrace that is closely related to formation of extra-long-staple (ELS) cultivated cotton. An SV-based pan-genome (Pan-SV) was then constructed using a total of 182 593 non-redundant SVs, including 2236 inversions, 97 398 insertions, and 82 959 deletions from 11 assembled genomes of allopolyploid cotton. The utility of this Pan-SV was then demonstrated through population structure analysis and genome-wide association studies (GWASs). Using segregation mapping populations produced through crossing ELS cotton and the landrace along with an SV-based GWAS, certain SVs responsible for speciation, domestication, and improvement in tetraploid cottons were identified. Importantly, some of the SVs presently identified as associated with the yield and fiber quality improvement had not been identified in previous SNP-based GWAS. In particular, a 9-bp insertion or deletion was found to associate with elimination of the interspecific reproductive isolation between Gossypium hirsutum and G. barbadense. Collectively, this study provides new insights into genome-wide, gene-scale SVs linked to important agronomic traits in a major crop species and highlights the importance of SVs during the speciation, domestication, and improvement of cultivated crop species.

Genome-wide identification and expression analysis of the HD2 protein family and its response to drought and salt stress in Gossypium species

Histone deacetylase 2 (HD2) proteins play an important role in the regulation of gene expression. This helps with the growth and development of plants and also plays a crucial role in responses to biotic and abiotic stress es. HD2s comprise a C₂H₂-type Zn²⁺ finger at their C-terminal and an HD2 label, deacetylation and phosphorylation sites, and NLS motifs at their N-terminal. In this study, a total of 27 HD2 members were identified, using Hidden Markov model profiles, in two diploid cotton genomes (Gossypium raimondii and Gossypium arboretum) and two tetraploid cotton genomes (Gossypium hirsutum and Gossypium barbadense). These cotton HD2 members were classified into 10 major phylogenetic groups (I-X), of which group III was found to be the largest with 13 cotton HD2 members. An evolutionary investigation showed that the expansion of HD2 members primarily occurred as a result of segmental duplication in paralogous gene pairs. Further qRT-PCR validation of nine putative genes using RNA-Seq data suggested that GhHDT3D.2 exhibits significantly higher levels of expression at 12h, 24h, 48h, and 72h of exposure to both drought and salt stress conditions compared to a control measure at 0h. Furthermore, gene ontology, pathways, and co-expression network study of GhHDT3D.2 gene affirmed their significance in drought and salt stress responses.

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.

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.

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.

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.

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.

An overview of salinity stress, mechanism of salinity tolerance and strategies for its management in cotton

Salinity stress is one of the primary threats to agricultural crops resulting in impaired crop growth and development. Although cotton is considered as reasonably salt tolerant, it is sensitive to salt stress at some critical stages like germination, flowering, boll formation, resulting in reduced biomass and fiber production. The mechanism of partial ion exclusion (exclusion of Na+ and/or Cl-) in cotton appears to be responsible for the pattern of uptake and accumulation of harmful ions (Na+ and Cl) in tissues of plants exposed to saline conditions. Maintaining high tissue K+/Na+ and Ca2+/Na+ ratios has been proposed as a key selection factor for salt tolerance in cotton. The key adaptation mechanism in cotton under salt stress is excessive sodium exclusion or compartmentation. Among the cultivated species of cotton, Egyptian cotton (Gossypium barbadense L.) exhibit better salt tolerance with good fiber quality traits as compared to most cultivated cotton and it can be used to improve five quality traits and transfer salt tolerance into Upland or American cotton (Gossypium hirsutum L.) by interspecific introgression. Cotton genetic studies on salt tolerance revealed that the majority of growth, yield, and fiber traits are genetically determined, and controlled by quantitative trait loci (QTLs). Molecular markers linked to genes or QTLs affecting key traits have been identified, and they could be utilized as an indirect selection criterion to enhance breeding efficiency through marker-assisted selection (MAS). Transfer of genes for compatible solute, which are an important aspect of ion compartmentation, into salt-sensitive species is, theoretically, a simple strategy to improve tolerance. The expression of particular stress-related genes is involved in plant adaptation to environmental stressors. As a result, enhancing tolerance to salt stress can be achieved by marker assisted selection added with modern gene editing tools can boost the breeding strategies that defend and uphold the structure and function of cellular components. The intent of this review was to recapitulate the advancements in salt screening methods, tolerant germplasm sources and their inheritance, biochemical, morpho-physiological, and molecular characteristics, transgenic approaches, and QTLs for salt tolerance in cotton.

Transcriptome analysis of the 2,4-dichlorophenoxyacetic acid (2,4-D)-tolerant cotton chromosome substitution line CS-B15sh and its susceptible parental lines G. hirsutum L. cv. Texas Marker-1 and G. barbadense L. cv. Pima 379

The cotton chromosome substitution line, CS-B15sh, exhibits 41% lower injury from 2,4-D when applied at the field recommended rate of 1.12 kg ae ha-1 (1×) than does Texas Marker-1 (TM-1). CS-B15sh was developed in the genetic background of Gossypium hirsutum L. cv TM-1 and has chromosome introgression on the short arm of chromosome 15 from Gossypium barbadense L. cv. Pima 379. In a previous experiment, we observed reduced translocation of [14C]2,4-D outside the treated leaf tissue in CS-B15sh, which contrasted with an increased translocation of the herbicide in the tissues above and below the treated leaf in TM-1. Our results indicate a potential 2,4-D tolerance mechanism in CS-B15sh involving altered movement of 2,4-D. Here, we used RNA sequencing (RNA-seq) to determine the differential expression of genes between 2,4-D-challenged and control plants of the tolerant (CS-B15sh) and susceptible lines (TM-1 and Pima 379). Several components of the 2,4-D/auxin-response pathway-including ubiquitin E3 ligase, PB1|AUX/IAA, ARF transcription factors, and F-box proteins of the SCFTIR1/AFB complex-were upregulated with at least threefold higher expression in TM-1 compared with CS-B15sh, while both Pima 379 and TM-1 showed the same fold change expression for PB1|AUX/IAA mRNA. Some genes associated with herbicide metabolism, including flavin monooxygenase (Gohir.A01G174100) and FAD-linked oxidase (Gohir.D06G002600), exhibited at least a twofold increase in CS-B15sh than in TM-1 (the gene was not expressed in Pima 379), suggesting a potential relationship between the gene's expression and 2,4-D tolerance. It is interesting to note that glutathione S-transferase was differentially expressed in both CS-B15sh and Pima 379 but not in TM-1, while cytochrome P450 and other genes involved in the oxidation-reduction process were significantly expressed only in CS-B15sh in response to 2,4-D. Gene set enrichment analysis on the union DEGs of the three cotton genotypes revealed the depletion of transcripts involved in photosynthesis and enrichment of transcripts involved in ABA response and signaling.

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.

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.

Differential seedling growth and tolerance indices reflect drought tolerance in cotton

BACKGROUND

Cotton production is adversely effected by drought stress. It is exposed to drought stress at various critical growth stages grown under a water scarcity environment. Roots are the sensors of plants; they detect osmotic stress under drought stress and play an important role in plant drought tolerance mechanisms. The seedling stage is very sensitive to drought stress, and it needed to explore the methods and plant characteristics that contribute to drought tolerance in cotton.

RESULTS

Initially, seedlings of 18 genotypes from three Gossypium species: G. hirsutum, G. barbadense, and G. arboreum, were evaluated for various seedling traits under control (NS) and drought stress (DS). Afterward, six genotypes, including two of each species, one tolerant and one susceptible, were identified based on the cumulative drought sensitivity response index (CDSRI). Finally, growth rates (GR) were examined for shoot and root growth parameters under control and DS in experimental hydroponic conditions. A significant variation of drought stress responses was observed across tested genotypes and species. CDSRI allowed here to identify the drought-sensitive and drought-resistant cultivar of each investigated species. Association among root and shoots growth traits disclosed influential effects of enduring the growth under DS. The traits including root length, volume, and root number were the best indicators with significantly higher differential responses in the tolerant genotypes. These root growth traits, coupled with the accumulation of photosynthates and proline, were also the key indicators of the resistance to drought stress.

CONCLUSION

Tolerant genotypes have advanced growth rates and the capacity to cop with drought stress by encouraging characteristics, including root differential growth traits coupled with physiological traits such as chlorophyll and proline contents. Tolerant and elite genotypes of G. hirsutum were more tolerant of drought stress than obsolete genotypes of G. barbadense and G. arboreum. Identified genotypes have a strong genetic basis of drought tolerance, which can be used in cotton breeding programs.

Genome-wide identification and functional analysis of class III peroxidases in Gossypium hirsutum

Class III peroxidase (PRX) genes play essential roles in various processes, such as auxin catabolism, removal of H₂O₂, crosslinking cell wall components, and response to biotic and abiotic stresses. In this study, we identified 166, 78 and 89 PRX genes from G. hirsutum, G. arboretum and G. raimondii, respectively. These PRX genes were classified into seven subfamilies based on phylogenetic tree analysis and the classification of PRX genes in Arabidopsis. Segmental duplication and purifying selection were the major factors driving the evolution of GhPRXs. GO and KEGG enrichment analysis revealed that GhPRX genes were mainly associated with responding to oxidative stresses, peroxidase activities and phenylpropanoid biosynthesis pathways. Transcriptome data analysis showed that GhPRX genes expression were significantly different in microspore development between the sterility line-JinA and the maintainer line MB177. We confirmed the up-regulation of GhPRX107 and down-regulation of GhPRX128 in the sterile line compared to its maintainer line using qRT-PCR, suggesting their roles in pollen fertility. In addition, silencing GhPRX107 in cotton showed a significant decrease of the reactive oxygen species (ROS) levels of microsporocyte stage anthers compared to control. Overexpressing GhPRX107 in Arabidopsis significantly increased the ROS levels of anthers compared to wild type. In conclusion, we identified GhPRX107 as a determinant of ROS levels in anther. This work sets a foundation for PRX studies in pollen development.

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.

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.

Construction of a high-density genetic map and identification of QTLs related to agronomic and physiological traits in an interspecific (Gossypium hirsutum × Gossypium barbadense) F2 population

Background

Advances in genome sequencing technology, particularly restriction-site associated DNA sequence (RAD-seq) and whole-genome resequencing, have greatly aided the construction of cotton interspecific genetic maps based on single nucleotide polymorphism (SNPs), Indels, and other types of markers. High-density genetic maps can improve accuracy of quantitative trait locus (QTL) mapping, narrow down location intervals, and facilitate identification of the candidate genes.

Result

In this study, 249 individuals from an interspecific F2 population (TM-1 and Hai7124) were re-sequenced, yielding 6303 high-confidence bin markers spanning 5057.13 cM across 26 cotton chromosomes. A total of 3380 recombination hot regions RHRs were identified which unevenly distributed on the 26 chromosomes. Based on this map, 112 QTLs relating to agronomic and physiological traits from seedling to boll opening stage were identified, including 15 loci associated with 14 traits that contained genes harboring nonsynonymous SNPs. We analyzed the sequence and expression of these ten candidate genes and discovered that GhRHD3 (GH_D10G0500) may affect fiber yield while GhGPAT6 (GH_D04G1426) may affect photosynthesis efficiency.

Conclusion

Our research illustrates the efficiency of constructing a genetic map using binmap and QTL mapping on the basis of a certain size of the early-generation population. High-density genetic map features high recombination exchanges in number and distribution. The QTLs and the candidate genes identified based on this high-density genetic map may provide important gene resources for the genetic improvement of 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.

GhCDNC and GhCYP706B1 genes mediate gossypol biosynthesis in upland cotton

BACKGROUND

In cotton the identification and characterization of natural defense is a cost-effective, sustainable, and environment-friendly strategy to combat cotton pests. The secondary metabolites traits in cotton plant i.e., toxic gossypol glands play significant role for development and self-defense mechanism. To utilize gossypol in breeding implements, the understanding of gossypol initiation biosynthesis genes has vital importance at reproductive organ development stages.

METHODS

Cotton germplasm of 100 genotypes screened visually based on gossypol glandedness and a core set of ten genotypes was developed. Further three genotypes FH-330 (high glanding), F-280 (low glanding) and IRMA-197 (glandless) were used for determining the transcript abundance of twelve gossypol biosynthesis genes.

RESULTS

Out of 100, germplasm categorized as (76) high glanding, (22) medium glanding), one genotype for each (low glanding) and (glandless) category. Real-time qPCR analysis revealed varied expression patterns among selected three genotypes. Out of twelve, three genes CYP706B1, CDNC and 2ODD-1 had strong expression levels in all tested tissues in high glanded genotype, while, slight or no expression of these genes was recorded in low glanding and glandless genotype, respectively. The shell of developing boll (10, 20, 30 DPA), and developing embryo (20, 30 DPA) showed substantially medium to maximum expression, respectively while high to medium expression was recorded in sepals and leaf tissue.

CONCLUSIONS

Our study demonstrated that CYP706B1, CDNC and 2ODD-1 are the most promising genes involved in gossypol biosynthesis. Developing boll shell, developing embryo, leaf and sepal also have significant ability to synthesize gossypol. This will provide scientists a way to manipulate gossypol contents in economically important organs of cotton plant for targeted breeding.

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.

Genome-wide characterization and expression analysis of Erf gene family in cotton

BACKGROUND

AP2/ERF transcription factors are important in a variety of biological activities, including plant growth, development, and responses to biotic and abiotic stressors. However, little study has been done on cotton's AP2/ERF genes, although cotton is an essential fibre crop. We were able to examine the tissue and expression patterns of AP2/ERF genes in cotton on a genome-wide basis because of the recently published whole genome sequence of cotton. Genome-wide analysis of ERF gene family within two diploid species (G. arboreum & G. raimondii) and two tetraploid species (G. barbadense, G. hirsutum) was performed.

RESULTS

A total of 118, 120, 213, 220 genes containing the sequence of single AP2 domain were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum respectively. The identified genes were unevenly distributed across 13/26 chromosomes of A and D genomes of cotton. Synteny and collinearity analysis revealed that segmental duplications may have played crucial roles in the expansion of the cotton ERF gene family, as well as tandem duplications played a minor role. Cis-acting elements of the promoter sites of Ghi-ERFs genes predict the involvement in multiple hormone responses and abiotic stresses. Transcriptome and qRT-PCR analysis revealed that Ghi-ERF-2D.6, Ghi-ERF-12D.13, Ghi-ERF-6D.1, Ghi-ERF-7A.6 and Ghi-ERF-11D.5 are candidate genes against salinity tolerance in upland cotton.

CONCLUSION

Overwhelmingly, the present study paves the way to better understand the evolution of cotton ERF genes and lays a foundation for future investigation of ERF genes in improving salinity stress tolerance in cotton.

Molecular traits and functional analysis of Rapid Alkalinization Factors (RALFs) in four Gossypium species

Rapid Alkalinization Factors (RALFs) are plant-secreted, cysteine-rich polypeptides which are known to play essential roles in plant developmental processes and in several defense mechanisms. So far, RALF polypeptides have not been investigated in the Gossypium genus. In this study, 42, 38, 104 and 120 RALFs were identified from diploid G. arboreum and G. raimondi and tetraploid G. hirsutum and G. barbadense, respectively. These were further divided into four groups. Protein characteristics, sequence alignment, gene structure, conserved motifs, chromosomal location and cis-element identification were comprehensively analyzed. Whole genome duplication (WGD) /segmental duplication may be the reason why the number of RALF genes doubled in tetraploid Gossypium species. Expression patterns analysis showed that GhRALFs had different transcript accumulation patterns in the tested tissues and were differentially expressed in response to various abiotic stresses. Furthermore, GhRALF41-3 over-expressing (OE) plants showed reduction in root length and developed later with short stems and small rosettes than that of the wild type. The GhRALF14-8 and GhRALF27-8 OE plants, especially the latter, showed increase in seed abortion. Both transgenic Arabidopsis and VIGS cotton demonstrate that three GhRALFs are negative regulators in response to salt stress. Our systematic analyses provided insights into the characterization of RALF genes in Gossypium, which forms genetic basis for further exploration in their potential applications in cotton production.

Core Collections: Is There Any Value for Cotton Breeding?

Global plant breeding activities are reliant on the available genetic variation held in extant varieties and germplasm collections. Throughout the mid- to late 1900s, germplasm collecting efforts were prioritized for breeding programs to archive precious material before it disappeared and led to the development of the numerous large germplasm resources now available in different countries. In recent decades, however, the maintenance and particularly the expansion of these germplasm resources have come under threat, and there has been a significant decline in investment in further collecting expeditions, an increase in global biosecurity restrictions, and restrictions placed on the open exchange of some commercial germplasm between breeders. The large size of most genebank collections, as well as constraints surrounding the availability and reliability of accurate germplasm passport data and physical or genetic characterization of the accessions in collections, limits germplasm utilization by plant breeders. To overcome these constraints, core collections, defined as a representative subset of the total germplasm collection, have gained popularity. Core collections aim to increase germplasm utilization by containing highly characterized germplasm that attempts to capture the majority of the variation in a whole collection. With the recent availability of many new genetic tools, the potential to unlock the value of these resources can now be realized. The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program supplies 100% of the cotton cultivars grown in Australia. The program is reliant on the use of plant genetic resources for the development of improved cotton varieties to address emerging challenges in pest and disease resistance as well as the global changes occurring in the climate. Currently, the CSIRO germplasm collection is actively maintained but underutilized by plant breeders. This review presents an overview of the Australian cotton germplasm resources and discusses the appropriateness of a core collection for 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.

The Evolution and Expression Profiles of EC1 Gene Family during Development in Cotton

Fertilization is essential to sexual reproduction of flowering plants. EC1 (EGG CELL 1) proteins have a conserved cysteine spacer characteristic and play a crucial role in double fertilization process in many plant species. However, to date, the role of EC1 gene family in cotton is fully unknown. Hence, detailed bioinformatics analysis was explored to elucidate the biological mechanisms of EC1 gene family in cotton. In this study, we identified 66 genes in 10 plant species in which a total of 39 EC1 genes were detected from cotton genome. Phylogenetic analysis clustered the identified EC1 genes into three families (I-III) and all of them contain Prolamin-like domains. A good collinearity was observed in the synteny analysis of the orthologs from cotton genomes. Whole-genome duplication was determined to be one of the major impetuses for the expansion of the EC1 gene family during the process of evolution. qRT-PCR analysis showed that EC1 genes were highly expressed in reproductive tissues under multiple stresses, signifying their potential role in enhancing stress tolerance or responses. Additionally, gene interaction networks showed that EC1 genes may be involved in cell stress and response transcriptional regulator in the synergid cells and activate the expression of genes required for pollen tube guidance. Our results provide novel functional insights into the evolution and functional elucidation of EC1 gene family in cotton.

Genome-wide analysis of the CalS gene family in cotton reveals their potential roles in fiber development and responses to stress

Callose deposition occurs during plant growth and development, as well as when plants are under biotic and abiotic stress. Callose synthase is a key enzyme for the synthesis of callose. In this study, 27, 28, 16, and 15 callose synthase family members were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium raimondii, and Gossypium arboreum using the sequence of Arabidopsis callose synthase. The CalSs were divided into five groups by phylogenetic, gene structure, and conservative motif analysis. The conserved motifs and gene structures of CalSs in each group were highly similar. Based on the analysis of cis-acting elements, it is inferred that GhCalSs were regulated by abiotic stress. WGD/Segmental duplication promoted the amplification of the CalS gene in cotton, and purification selection had an important function in the CalS family. The transcriptome data and qRT-PCR under cold, heat, salt, and PEG treatments showed that GhCalSs were involved in abiotic stress. The expression patterns of GhCalSs were different in various tissues. We predicted that GhCalS4, which was highly expressed in fibers, had an important effect on fiber elongation. Hence, these results help us understand the role of GhCalSs in fiber development and stress response.

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 Identification and Characterization of GhCOMT Gene Family during Fiber Development and Verticillium Wilt Resistance in Cotton

Caffeic acid O-methyltransferases (COMTs) play an essential role in lignin synthesis procession, especially in the plant’s phenylalanine metabolic pathway. The content of COMT genes in cotton and the relationship between their expression patterns have not been studied clearly in cotton. In this study, we have identified 190 COMT genes in cotton, which were classified into three groups (I, II and III), and mapped on the cotton chromosomes. In addition, we found that 135 of the 190 COMT genes result from dispersed duplication (DSD) and whole-genome duplication (WGD), indicating that DSD and WGD were the main forces driving COMT gene expansion. The Ka/Ks analysis showed that GhCOMT43 and GhCOMT41 evolved from GaCOMT27 and GrCOMT14 through positive selection. The results of qRT-PCR showed that GhCOMT13, GhCOMT28, GhCOMT39 and GhCOMT55 were related to lignin content during the cotton fiber development. GhCOMT28, GhCOMT39, GhCOMT55, GhCOMT56 and GhCOMT57 responded to Verticillium Wilt (VW) and maybe related to VW resistance through lignin synthesis. Conclusively, this study found that GhCOMTs were highly expressed in the secondary wall thickening stage and VW. These results provide a clue for studying the functions of GhCOMTs in the development of cotton fiber and VW resistance and could lay a foundation for breeding cotton cultivates with higher quantity and high resistance to VW.

Genome-wide identification and analysis of Oleosin gene family in four cotton species and its involvement in oil accumulation and germination

BACKGROUND

Cotton is not only a major textile fiber crop but also a vital oilseed, industrial, and forage crop. Oleosins are the structural proteins of oil bodies, influencing their size and the oil content in seeds. In addition, the degradation of oleosins is involved in the mobilization of lipid and oil bodies during seed germination. However, comprehensive identification and the systematic analysis of the Oleosin gene (OLEOs) family have not been conducted in cotton.

RESULTS

An in-depth analysis has enabled us to identify 25 and 24 OLEOs in tetraploid cotton species G. hirsutum and G. barbadense, respectively, while 12 and 13 OLEOs were identified in diploid species G. arboreum and G. raimondii, respectively. The 74 OLEOs were further clustered into three lineages according to the phylogenetic tree. Synteny analysis revealed that most of the OLEOs were conserved and that WGD or segmental duplications might drive their expansion. The transmembrane helices in GhOLEO proteins were predicted, and three transmembrane models were summarized, in which two were newly proposed. A total of 24 candidate miRNAs targeting GhOLEOs were predicted. Three highly expressed oil-related OLEOs, GH_A07G0501 (SL), GH_D10G0941 (SH), and GH_D01G1686 (U), were cloned, and their subcellular localization and function were analyzed. Their overexpression in Arabidopsis increased seed oil content and decreased seed germination rates.

CONCLUSION

We identified OLEO gene family in four cotton species and performed comparative analyses of their relationships, conserved structure, synteny, and gene duplication. The subcellular localization and function of three highly expressed oil-related OLEOs were detected. These results lay the foundation for further functional characterization of OLEOs and improving seed oil content.

Design and application of hydrophilic bimetallic metal-organic framework magnetic nanoparticles for rapid capture of exosomes

Functional materials with good biocompatibility have been widely used in the study of genomics, proteomics and disease diagnosis, which has improved the progress of life science. In this paper, the material not only exhibited a strong affinity to the phosphate groups on the exosomal membrane due to the coexistence of Zr-O clusters and Ti⁴⁺, but also owned great hydrophilicity to reduce non-specific adsorption of contaminated proteins, achieving the separation and purification of exosomes from complex biosamples. The model exosomes extracted by ultracentrifugation (UC) were used to evaluate the feasibility of Fe₃O₄@UiO-66-NH₂@PA-Ti⁴⁺ capturing exosomes. The process of Fe₃O₄@UiO-66-NH₂@PA-Ti⁴⁺ capturing exosomes was simple to operate with a high recovery rate (97.3%) within a short time (5 min). Then Fe₃O₄@UiO-66-NH₂@PA-Ti⁴⁺ was further applied to capture exosomes in media and urine followed by the downstream proteomics analysis. 348 and 284 exosomal proteins were identified for cell medium and urine, respectively. This work shows great potential of the material for subsequent function research of disease-related exosomes by separating exosomes rapidly and efficiently.

Comparative Dynamic Transcriptome Reveals the Delayed Secondary-Cell-Wall Thickening Results in Altered Lint Percentage and Fiber Elongation in a Chromosomal Segment Substitution Line of Cotton (Gossypium hirsutum L.)

Lint percentage (LP) is an important yield component in cotton that is usually affected by initial fiber number and cell wall thickness. To explore how fiber cell wall development affects LP, phenotypic identification and dynamic transcriptome analysis were conducted using a single segment substitution line of chromosome 15 (SL15) that harbors a major quantitative trait locus (QTL) for LP. Compared to its recurrent parent LMY22, SL15 did not differ in initial fiber number, but the fiber cell wall thickness and single-fiber weight decreased significantly, altering LP. The comparative transcriptome profiles revealed that the secondary cell wall (SCW) development phase of SL15 was relatively delayed. Meanwhile, the expression of genes related to cell expansion decreased more slightly in SL15 with fiber development, resulting in relatively higher expression at SL15_25D than at LMY22_25D. SCW development-related genes, such as GhNACs and GhMYBs, in the putative NAC-MYB-CESA network differentially expressed at SL15_25D, along with the lower expression of CESA6, CSLC12, and CSLA2. The substituted chromosomal interval was further investigated, and found 6 of 146 candidate genes were differentially expressed in all four cell development periods including 10, 15, 20 and 25 DPA. Genetic variation and co-expression analysis showed that GH_D01G0052, GH_D01G0099, GH_D01G0100, and GH_D01G0140 may be important candidate genes associated with qLP-C15-1. Our results provide novel insights into cell wall development and its relationship with LP, which is beneficial for lint yield and fiber quality improvement.

High-quality genome assembly and resequencing of modern cotton cultivars provide resources for crop improvement

Cotton produces natural fiber for the textile industry. The genetic effects of genomic structural variations underlying agronomic traits remain unclear. Here, we generate two high-quality genomes of Gossypium hirsutum cv. NDM8 and Gossypium barbadense acc. Pima90, and identify large-scale structural variations in the two species and 1,081 G. hirsutum accessions. The density of structural variations is higher in the D-subgenome than in the A-subgenome, indicating that the D-subgenome undergoes stronger selection during species formation and variety development. Many structural variations in genes and/or regulatory regions potentially influencing agronomic traits were discovered. Of 446 significantly associated structural variations, those for fiber quality and Verticillium wilt resistance are located mainly in the D-subgenome and those for yield mainly in the A-subgenome. Our research provides insight into the role of structural variations in genotype-to-phenotype relationships and their potential utility in crop improvement.

High-quality genomes of two cultivated tetraploid cottons Gossypium hirsutum cv. NDM8 and Gossypium barbadense acc. Pima90 and resequencing of 1,081 G. hirsutum accessions provide insights into the role of structural variations.

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.

Genome-Wide Identification and Evolutionary Analysis of Gossypium Tubby-Like Protein (TLP) Gene Family and Expression Analyses During Salt and Drought Stress

Tubby-like proteins (TLPs) possess a highly conserved closed β barrel tubby domain at C-terminal and N-terminal F-box. The role of TLP gene family members has been widely discussed in numerous organisms; however, the detailed genome-wide study of this gene family in Gossypium species has not been reported till date. Here, we systematically identified 105 TLP gene family members in cotton (Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense) genomes and classified them into eight phylogenetic groups. Cotton TLP12 gene family members clustered into two groups, 4 and 8. They experienced higher evolutionary pressure in comparison to others, indicating the faster evolution in both diploid as well as in tetraploid cotton. Cotton TLP gene family members expanded mainly due to segmental duplication, while only one pair of tandem duplication was found in cotton TLPs paralogous gene pairs. Subsequent qRT-PCR validation of seven putative key candidate genes of GhTLPs indicated that GhTLP11A and GhTLP12A.1 genes were highly sensitive to salt and drought stress. The co-expression network, pathways, and cis-regulatory elements of GhTLP11A and GhTLP12A.1 genes confirmed their functional importance in salt and drought stress responses. This study proposes the significance of GhTLP11A and GhTLP12A.1 genes in exerting control over salt and drought stress responses in G. hirsutum and also provides a reference for future research, elaborating the biological roles of G. hirsutum TLPs in both stress responses.

Population structure and genetic diversity of the Pee Dee cotton breeding program

Accelerated marker-assisted selection and genomic selection breeding systems require genotyping data to select the best parents for combining beneficial traits. Since 1935, the Pee Dee (PD) cotton germplasm enhancement program has developed an important genetic resource for upland cotton (Gossypium hirsutum L.), contributing alleles for improved fiber quality, agronomic performance, and genetic diversity. To date, a detailed genetic survey of the program's eight historical breeding cycles has yet to be undertaken. The objectives of this study were to evaluate genetic diversity across and within-breeding groups, examine population structure, and contextualize these findings relative to the global upland cotton gene pool. The CottonSNP63K array was used to identify 17,441 polymorphic markers in a panel of 114 diverse PD genotypes. A subset of 4597 markers was selected to decrease marker density bias. Identity-by-state pairwise distance varied substantially, ranging from 0.55 to 0.97. Pedigree-based estimates of relatedness were not very predictive of observed genetic similarities. Few rare alleles were present, with 99.1% of SNP alleles appearing within the first four breeding cycles. Population structure analysis with principal component analysis, discriminant analysis of principal components, fastSTRUCTURE, and a phylogenetic approach revealed an admixed population with moderate substructure. A small core collection (n < 20) captured 99% of the program's allelic diversity. Allele frequency analysis indicated potential selection signatures associated with stress resistance and fiber cell growth. The results of this study will steer future utilization of the program's germplasm resources and aid in combining program-specific beneficial alleles and maintaining genetic diversity.

Transcriptome analysis reveals genes potentially related to high fiber strength in a Gossypium hirsutum line IL9 with Gossypium mustelinum introgression

Cotton (Gossypium L.) is the most important fiber crop worldwide. Here, transcriptome analysis was conducted on developing fibers of a G. mustelinum introgression line, IL9, and its recurrent parent, PD94042, at 17 and 21 days post-anthesis (dpa). Differentially expressed genes (DEGs) of PD94042 and IL9 were identified. Gene Ontology (GO) enrichment analysis showed that the annotated DEGs were rich in two main biological processes and two main molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis likewise showed that the annotated DEGs were mainly enriched in metabolic pathways and biosynthesis of secondary metabolites. In total, 52 DEGs were selected as candidate genes based on comparison of the DEGs and GO function annotation information. Quantitative real-time PCR (RT-qPCR) analysis results for 12 randomly selected DEGs were consistent with transcriptome analysis. SNP identification based on G. mustelinum chromatin segment introgression showed that 394 SNPs were identified in 268 DEGs, and two genes with known functions were identified within fiber strength quantitative trait loci (QTL) regions or near the confidence intervals. We identified 52 key genes potentially related to high fiber strength in a G. mustelinum introgression line and provided significant insights into the study of cotton fiber quality improvement.

GhKWL1 Upregulates GhERF105 but Its Function Is Impaired by Binding with VdISC1, a Pathogenic Effector of Verticillium dahliae

Verticillium wilt, caused by Verticillium dahliae, is a devastating disease for many important crops, including cotton. Kiwellins (KWLs), a group of cysteine-rich proteins synthesized in many plants, have been shown to be involved in response to various phytopathogens. To evaluate genes for their function in resistance to Verticillium wilt, we investigated KWL homologs in cotton. Thirty-five KWL genes (GhKWLs) were identified from the genome of upland cotton (Gossypium hirsutum). Among them, GhKWL1 was shown to be localized in nucleus and cytosol, and its gene expression is induced by the infection of V. dahliae. We revealed that GhKWL1 was a positive regulator of GhERF105. Silencing of GhKWL1 resulted in a decrease, whereas overexpression led to an increase in resistance of transgenic plants to Verticillium wilt. Interestingly, through binding to GhKWL1, the pathogenic effector protein VdISC1 produced by V. dahliae could impair the defense response mediated by GhKWL1. Therefore, our study suggests there is a GhKWL1-mediated defense response in cotton, which can be hijacked by V. dahliae through the interaction of VdISC1 with GhKWL1.

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 Comparative Analysis of Heat Shock Transcription Factors Provides Novel Insights for Evolutionary History and Expression Characterization in Cotton Diploid and Tetraploid Genomes

Heat shock transcription factors (HSFs) are involved in environmental stress response and plant development, such as heat stress and flowering development. According to the structural characteristics of the HSF gene family, HSF genes were classified into three major types (HSFA, HSFB, and HSFC) in plants. Using conserved domains of HSF genes, we identified 621 HSF genes among 13 cotton genomes, consisting of eight diploid and five tetraploid genomes. Phylogenetic analysis indicated that HSF genes among 13 cotton genomes were grouped into two different clusters: one cluster contained all HSF genes of HSFA and HSFC, and the other cluster contained all HSF genes of HSFB. Comparative analysis of HSF genes in Arabidopsis thaliana, Gossypium herbaceum (A1), Gossypium arboreum (A2), Gossypium raimondii (D5), and Gossypium hirsutum (AD1) genomes demonstrated that four HSF genes were inherited from a common ancestor, A0, of all existing cotton A genomes. Members of the HSF gene family in G. herbaceum (A1) genome indicated a significant loss compared with those in G. arboretum (A2) and G. hirsutum (AD1) A genomes. However, HSF genes in G. raimondii (D5) showed relative loss compared with those in G. hirsutum (AD1) D genome. Analysis of tandem duplication (TD) events of HSF genes revealed that protein-coding genes among different cotton genomes have experienced TD events, but only the two-gene tandem array was detected in Gossypium thurberi (D1) genome. The expression analysis of HSF genes in G. hirsutum (AD1) and Gossypium barbadense (AD2) genomes indicated that the expressed HSF genes were divided into two different groups, respectively, and the expressed HSF orthologous genes between the two genomes showed totally different expression patterns despite the implementation of the same abiotic stresses. This work will provide novel insights for the study of evolutionary history and expression characterization of HSF genes in different cotton genomes and a widespread application model for the study of HSF gene families in plants.

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.

Parallel and Intertwining Threads of Domestication in Allopolyploid Cotton

The two cultivated allopolyploid cottons, Gossypium hirsutum and Gossypium barbadense, represent a remarkable example of parallel independent domestication, both involving dramatic morphological transformations under selection from wild perennial plants to annualized row crops. Deep resequencing of 643 newly sampled accessions spanning the wild-to-domesticated continuum of both species, and their allopolyploid relatives, are combined with existing data to resolve species relationships and elucidate multiple aspects of their parallel domestication. It is confirmed that wild G. hirsutum and G. barbadense were initially domesticated in the Yucatan Peninsula and NW South America, respectively, and subsequently spread under domestication over 4000-8000 years to encompass most of the American tropics. A robust phylogenomic analysis of infraspecific relationships in each species is presented, quantify genetic diversity in both, and describe genetic bottlenecks associated with domestication and subsequent diffusion. As these species became sympatric over the last several millennia, pervasive genome-wide bidirectional introgression occurred, often with striking asymmetries involving the two co-resident genomes of these allopolyploids. Diversity scans revealed genomic regions and genes unknowingly targeted during domestication and additional subgenomic asymmetries. These analyses provide a comprehensive depiction of the origin, divergence, and adaptation of cotton, and serve as a rich resource for cotton improvement.

Genome-wide identification and expression analysis of the plant specific LIM genes in Gossypium arboreum under phytohormone, salt and pathogen stress

Asiatic cotton (Gossypium arboreum) cultivated as ‘desi cotton’ in India, is renowned for its climate resilience and robustness against biotic and abiotic stresses. The genome of G. arboreum is therefore, considered as a valued reserve of information for discovering novel genes or gene functions for trait improvements in the present context of cotton cultivation world-wide. In the present study, we carried out genome-wide analysis of LIM gene family in desi cotton and identified twenty LIM domain proteins (GaLIMs) which include sixteen animals CRP-like GaLIMs and four plant specific GaLIMs with presence (GaDA1) or absence (GaDAR) of UIM (Ubiquitin Interacting Motifs). Among the sixteen CRP-like GaLIMs, eleven had two conventional LIM domains while, five had single LIM domain which was not reported in LIM gene family of the plant species studied, except in Brassica rapa. Phylogenetic analysis of these twenty GaLIM proteins in comparison with LIMs of Arabidopsis, chickpea and poplar categorized them into distinct αLIM1, βLIM1, γLIM2, δLIM2 groups in CRP-like LIMs, and GaDA1 and GaDAR in plant specific LIMs group. Domain analysis had revealed consensus [(C-X₂-C-X₁₇-H-X₂-C)-X₂-(C-X₂-C-X₁₇-C-X₂-H)] and [(C-X₂-C-X₁₇-H-X₂-C)-X₂-(C-X₄-C-X₁₅-C-X₂-H)] being conserved as first and/or second LIM domains of animal CRP-like GaLIMs, respectively. Interestingly, single LIM domain containing GaLIM15 was found to contain unique consensus with longer inter-zinc-motif spacer but shorter second zinc finger motif. All twenty GaLIMs showed variable spatio-temporal expression patterns and accordingly further categorized into distinct groups of αLIM1, βLIM1, γLIM2 δLIM2 and plant specific LIM (DA1/DAR). For the first time, response of GaDA1/DAR under the influence of biotic and abiotic stresses were studied in cotton, involving treatments with phytohormones (Jasmonic acid and Abscisic acid), salt (NaCl) and wilt causing pathogen (Fusarium oxysporum). Expressions patterns of GaDA1/DAR showed variable response and identified GaDA2 as a probable candidate gene for stress tolerance in G. arboreum.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development

BACKGROUND

Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored.

RESULTS

Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes.

CONCLUSION

Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.