Verticillium wilt (VW), also known as “cotton cancer,” is one of the most destructive diseases in global cotton production that seriously impacts fiber yield and quality. Despite numerous attempts, little significant progress has been made in improving the VW resistance of upland cotton. The development of chromosome segment substitution lines (CSSLs) from Gossypium hirsutum × G. barbadense has emerged as a means of simultaneously developing new cotton varieties with high-yield, superior fiber, and resistance to VW. In this study, VW-resistant investigations were first conducted in an artificial greenhouse, a natural field, and diseased nursery conditions, resulting in the identification of one stably VW-resistant CSSL, MBI8255, and one VW-susceptible G. hirsutum, CCRI36, which were subsequently subjected to biochemical tests and transcriptome sequencing during V991 infection (0, 1, and 2 days after inoculation). Eighteen root samples with three replications were collected to perform multiple comparisons of enzyme activity and biochemical substance contents. The findings indicated that VW resistance was positively correlated with peroxidase and polyphenol oxidase activity, but negatively correlated with malondialdehyde content. Additionally, RNA sequencing was used for the same root samples, resulting in a total of 77,412 genes, of which 23,180 differentially expressed genes were identified from multiple comparisons between samples. After Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the expression profiles identified using Short Time-series Expression Miner, we found that the metabolic process in the biological process, as well as the pathways of phenylpropanoid biosynthesis and plant hormone signal transduction, participated significantly in the response to VW. Gene functional annotation and expression quantity analysis indicated the important roles of the phenylpropanoid metabolic pathway and oxidation-reduction process in response to VW, which also provided plenty of candidate genes related to plant resistance. This study concentrates on the preliminary response to V991 infection by comparing the VW-resistant CSSL and its VW-susceptible recurrent parent. Not only do our findings facilitate the culturing of new resistant varieties with high yield and superior performance, but they also broaden our understanding of the mechanisms of cotton resistance to VW.

Increasing soil salinization has led to severe reductions in plant yield and quality, and investigating the molecular mechanism of salt stress response is therefore an urgent priority. In this study, we systematically analyzed the response of cotton roots to salt stress using single-cell transcriptomics technology; 56 281 high-quality cells were obtained from 5-day-old lateral root tips of Gossypium arboreum under natural growth conditions and different salt treatments. Ten cell types with an array of novel marker genes were identified and confirmed by in situ RNA hybridization, and pseudotime analysis of some specific cell types revealed their potential differentiation trajectories. Prominent changes in cell numbers under salt stress were observed for outer epidermal and inner endodermal cells, which were significantly enriched in response to stress, amide biosynthetic process, glutathione metabolism, and glycolysis/gluconeogenesis. Analysis of differentially expressed genes identified in multiple comparisons revealed other functional aggregations concentrated on plant-type primary cell wall biogenesis, defense response, phenylpropanoid biosynthesis, and metabolic pathways. Some candidate differentially expressed genes encoding transcription factors or associated with plant hormones also responsive to salt stress were identified, and the function of Ga03G2153, annotated as auxin-responsive GH3.6, was confirmed by virus-induced gene silencing. The GaGH3.6-silenced plants showed a severe stress-susceptible phenotype, and physiological and biochemical measurements indicated that they suffered more significant oxidative damage. These results suggest that GaGH3.6 might participate in cotton salt tolerance by regulating redox processes. We thus construct a transcriptional atlas of salt-stressed cotton roots at single-cell resolution, enabling us to explore cellular heterogeneity and differentiation trajectories and providing valuable insights into the molecular mechanisms that underlie plant stress tolerance.

How to develop new cotton varieties possessing high yield traits of Upland cotton and superior fiber quality traits of Sea Island cotton remains a key task for cotton breeders and researchers. While multiple attempts bring in little significant progresses, the development of Chromosome Segment Substitution Lines (CSSLs) from Gossypium barbadense in G. hirsutum background provided ideal materials for aforementioned breeding purposes in upland cotton improvement. Based on the excellent fiber performance and relatively clear chromosome substitution segments information identified by Simple Sequence Repeat (SSR) markers, two CSSLs, MBI9915 and MBI9749, together with the recurrent parent CCRI36 were chosen to conduct transcriptome sequencing during the development stages of fiber elongation and Secondary Cell Wall (SCW) synthesis (from 10DPA and 28DPA), aiming at revealing the mechanism of fiber development and the potential contribution of chromosome substitution segments from Sea Island cotton to fiber development of Upland cotton. In total, 15 RNA-seq libraries were constructed and sequenced separately, generating 705.433 million clean reads with mean GC content of 45.13% and average Q30 of 90.26%. Through multiple comparisons between libraries, 1801 differentially expressed genes (DEGs) were identified, of which the 902 up-regulated DEGs were mainly involved in cell wall organization and response to oxidative stress and auxin, while the 898 down-regulated ones participated in translation, regulation of transcription, DNA-templated and cytoplasmic translation based on GO annotation and KEGG enrichment analysis. Subsequently, STEM software was performed to explicate the temporal expression pattern of DEGs. Two peroxidases and four flavonoid pathway-related genes were identified in the “oxidation-reduction process”, which could play a role in fiber development and quality formation. Finally, the reliability of RNA-seq data was validated by quantitative real-time PCR of randomly selected 20 genes. The present report focuses on the similarities and differences of transcriptome profiles between the two CSSLs and the recurrent parent CCRI36 and provides novel insights into the molecular mechanism of fiber development, and into further exploration of the feasible contribution of G. barbadense substitution segments to fiber quality formation, which will lay solid foundation for simultaneously improving fiber yield and quality of upland cotton through CSSLs.

Pectin methylesterase (PME, EC 3.1.1.11) is a hydrolytic enzyme that utilizes pectin as substrates, and plays a significant role in regulating pectin reconstruction thereby regulating plant growth. Pectin is one of the important components of the plant cell wall, which forms the main structural material of cotton fiber. In this research, cotton genome information was used to identify PMEs. We identified 80 (GaPME01-GaPME80) PME genes from diploid G. arboreum (A genome), 78 (GrPME01-GrPME78) PME genes from G. raimondii (D genome), and 135 (GhPME001-GhPME135) PME genes from tetraploid cotton G. hirsutum (AD genome). We further analyzed their gene structure, conserved domain, gene expression, and systematic evolution to lay the foundation for deeper research on the function of PMEs. Phylogenetic data indicated that members from the same species demonstrated relatively high sequence identities and genetic similarities. Analysis of gene structures showed that most of the PMEs genes had 2–3 exons, with a few having a variable number of exons from 4 to 6. There are nearly no differences in the gene structure of PMEs among the three (two diploid and one tetraploid) cotton species. Selective pressure analysis showed that the Ka/Ks value for each of the three cotton species PME families was less than one. Conserved domain analysis showed that PMEs members had a relatively conserved C-terminal pectinesterase domain (PME) while the N-terminus was less conserved. Moreover, some of the family members contained a pectin methylesterase inhibitor (PMEI) domain. The Ka/Ks ratios suggested that the duplicated PMEs underwent purifying selection after the duplication events. This study provided an important basis for further research on the functions of cotton PMEs. Results from qRT-PCR indicated that the expression level of different PMEs at various fiber developmental stages was different. Moreover, some of the PMEs showed fiber predominant expression in secondary wall thickening indicating tissue-specific expression patterns.

Cotton (Gossypium spp. L.) is a major origin of natural fiber, and is projected at 117 million bales worldwide for 2021/22. A variety of biotic and abiotic stresses have considerable negative impacts on cotton. The significantly decreased applications of chemical insecticidal sprays in the agro-ecosystem have greatly affected the biodiversity and dynamics of primary and secondary insects. Various control measures were taken around the globe to increase production costs. Temperature, drought, and salinity, and biotic stresses such as bacteria, viruses, fungi, nematodes, insects, and mites cause substantial losses to cotton crops. Here, we summarize a number of biotic and abiotic stresses upsetting Bt cotton crop with present and future biotechnology solution strategies that include a refuge strategy, multi-gene pyramiding, the release of sterile insects, seed mixing, RNAi, CRISPR/Cas9, biotic signaling, and the use of bioagents. Surveillance of insect resistance, monitoring of grower compliance, and implementation of remedial actions can lead to the sustainable use of cotton across the globe.