Genomics-enabled analysis of the emergent disease cotton bacterial blight

Comparative Genomics Identifies Conserved and Variable TAL Effectors in African Strains of the Cotton Pathogen Xanthomonas citri pv. malvacearum

Strains of Xanthomonas citri pv. malvacearum cause bacterial blight of cotton, a potentially serious threat to cotton production worldwide, including in sub-Saharan countries. Development of disease symptoms, such as water soaking, has been linked to the activity of a class of type 3 effectors, called transcription activator-like (TAL) effectors, which induce susceptibility genes in the host's cells. To gain further insight into the global diversity of the pathogen, to elucidate their repertoires of TAL effector genes, and to better understand the evolution of these genes in the cotton-pathogenic xanthomonads, we sequenced the genomes of three African strains of X. citri pv. malvacearum using nanopore technology. We show that the cotton-pathogenic pathovar of X. citri is a monophyletic lineage containing at least three distinct genetic subclades, which appear to be mirrored by their repertoires of TAL effectors. We observed an atypical level of TAL effector gene pseudogenization, which might be related to resistance genes that are deployed to control the disease. Our work thus contributes to a better understanding of the conservation and importance of TAL effectors in the interaction with the host plant, which can inform strategies for improving resistance against bacterial blight in cotton.

Gene B5 in Cotton Confers High and Broad Resistance to Bacterial Blight and Conditions High Amounts of Sesquiterpenoid Phytoalexins

Bacterial blight resistance gene B₅ has received little attention since it was first described in 1950. A near-isogenic line (NIL) of Gossypium hirsutum cotton, AcB₅, was generated in an otherwise bacterial-blight-susceptible 'Acala 44' background. The introgressed locus B₅ in AcB₅ conferred strong and broad-spectrum resistance to bacterial blight. Segregation patterns of test crosses under Oklahoma field conditions indicated that AcB₅ is likely homozygous for resistance at two loci with partial dominance gene action. In controlled-environment conditions, two of the four copies of B₅ were required for effective resistance. Contrary to expectations of gene-for-gene theory, AcB₅ conferred high resistance toward isogenic strains of Xanthomonas citri subsp. malvacearum carrying cloned avirulence genes avrB4, avrb7, avrBIn, avrB101, and avrB102, respectively, and weaker resistance toward the strain carrying cloned avrb6. The hypothesis that each B gene, in the absence of a polygenic complex, triggers sesquiterpenoid phytoalexin production was tested by measurement of cadalene and lacinilene phytoalexins during resistant responses in five NILs carrying different B genes, four other lines carrying multiple resistance genes, as well as susceptible Ac44E. Phytoalexin production was an obvious, but variable, response in all nine resistant lines. AcB₅ accumulated an order of magnitude more of all four phytoalexins than any of the other resistant NILs. Its total levels were comparable to those detected in OK1.2, a highly resistant line that possesses several B genes in a polygenic background.

Improving cassava bacterial blight resistance by editing the epigenome

Pathogens rely on expression of host susceptibility (S) genes to promote infection and disease. As DNA methylation is an epigenetic modification that affects gene expression, blocking access to S genes through targeted methylation could increase disease resistance. Xanthomonas phaseoli pv. manihotis, the causal agent of cassava bacterial blight (CBB), uses transcription activator-like20 (TAL20) to induce expression of the S gene MeSWEET10a. In this work, we direct methylation to the TAL20 effector binding element within the MeSWEET10a promoter using a synthetic zinc-finger DNA binding domain fused to a component of the RNA-directed DNA methylation pathway. We demonstrate that this methylation prevents TAL20 binding, blocks transcriptional activation of MeSWEET10a in vivo and that these plants display decreased CBB symptoms while maintaining normal growth and development. This work therefore presents an epigenome editing approach useful for crop improvement.

A comparison of ImageJ and machine learning based image analysis methods to measure cassava bacterial blight disease severity

BACKGROUND

Methods to accurately quantify disease severity are fundamental to plant pathogen interaction studies. Commonly used methods include visual scoring of disease symptoms, tracking pathogen growth in planta over time, and various assays that detect plant defense responses. Several image-based methods for phenotyping of plant disease symptoms have also been developed. Each of these methods has different advantages and limitations which should be carefully considered when choosing an approach and interpreting the results.

RESULTS

In this paper, we developed two image analysis methods and tested their ability to quantify different aspects of disease lesions in the cassava-Xanthomonas pathosystem. The first method uses ImageJ, an open-source platform widely used in the biological sciences. The second method is a few-shot support vector machine learning tool that uses a classifier file trained with five representative infected leaf images for lesion recognition. Cassava leaves were syringe infiltrated with wildtype Xanthomonas, a Xanthomonas mutant with decreased virulence, and mock treatments. Digital images of infected leaves were captured overtime using a Raspberry Pi camera. The image analysis methods were analyzed and compared for the ability to segment the lesion from the background and accurately capture and measure differences between the treatment types.

CONCLUSIONS

Both image analysis methods presented in this paper allow for accurate segmentation of disease lesions from the non-infected plant. Specifically, at 4-, 6-, and 9-days post inoculation (DPI), both methods provided quantitative differences in disease symptoms between different treatment types. Thus, either method could be applied to extract information about disease severity. Strengths and weaknesses of each approach are discussed.

Cotton miR393-TIR1 Module Regulates Plant Defense Against Verticillium dahliae via Auxin Perception and Signaling

Plant auxin is essential in plant growth and development. However, the molecular mechanisms of auxin involvement in plant immunity are unclear. Here, we addressed the function of the cotton (Gossypium hirsutum) miR393-TIR1 module in plant defense against Verticillium dahliae infection via auxin perception and signaling. GhTIR1 was directedly cleaved by ghr-miR393 according to mRNA degradome data, 5'-RACE analysis, and a GUS reporter assay. Ghr-miR393 knockdown significantly increased plant susceptibility to V. dahliae compared to the control, while ghr-miR393 overexpression and GhTIR1 knockdown significantly increased plant resistance. External indole-3-acetic acid (IAA) application significantly enhanced susceptibility to V. dahliae in ghr-miR393 knockdown and control plants compared to mock treatment, and only slightly increased susceptibility in overexpressing ghr-miR393 and GhTIR1-silenced plants. Application of external PEO-IAA (an auxin antagonist) had a contrary trend with IAA application. Based on yeast two-hybrid and bimolecular fluorescence complementation assays, GhTIR1 interacted with GhIAA14 in the nucleus, and GhIAA14 knockdown reduced plant resistance to V. dahliae infection. The results suggested that the ghr-miR393-GhTIR1 module regulates plant defense via auxin perception and signaling. Additionally, simultaneous knockdown of GhTIR1 and GhICS1 significantly increased plant susceptibility to V. dahliae compared to the control, indicating that salicylic acid (SA) accumulation is vital for the ghr-miR393-GhTIR1 module to regulates plant resistance. Transcriptome data also demonstrated that GhTIR1 knockdown significantly downregulated expression of auxin-related genes and upregulated expression of SA-related genes. Overall, the ghr-miR393-GhTIR1 module participates in plant response to V. dahliae infection via IAA perception and signaling partially depending on the SA defense pathway.

Patterns of Seed-to-Seedling Transmission of Xanthomonas citri pv. malvacearum, the Causal Agent of Cotton Bacterial Blight

Cotton bacterial blight (CBB), caused by Xanthomonas citri pv. malvacearum, was a major disease of cotton in the United States in the early part of the twentieth century. The reemergence of CBB revealed many gaps in our understanding of this important disease. In this study, we employed a wild-type (WT) field isolate of X. citri pv. malvacearum from Georgia (U.S.A.) to generate a nonpathogenic hrcV mutant lacking a functional type-III secretion system (T3SS-). We tagged the WT and T3SS- strains with an auto-bioluminescent Tn7 reporter and compared colonization patterns of CBB-susceptible and CBB-resistant cotton seedlings using macroscopic image analysis and bacterial load enumeration. WT and T3SS- X. citri pv. malvacearum strains colonized cotton cotyledons of CBB-resistant and CBB-susceptible cotton cultivars. However, X. citri pv. malvacearum populations were significantly higher in CBB-susceptible seedlings inoculated with the WT strain. Additionally, WT and T3SS- X. citri pv. malvacearum strains systemically colonized true leaves, although at different rates. Finally, we observed that seed-to-seedling transmission of X. citri pv. malvacearum may involve systemic spread through the vascular tissue of cotton plants. These findings yield novel insights into potential X. citri pv. malvacearum reservoirs for CBB outbreaks.

SWEET genes and TAL effectors for disease resistance in plants: Present status and future prospects

SWEET genes encode sugar transporter proteins and often function as susceptibility (S) genes. Consequently, the recessive alleles of these SWEET genes provide resistance. This review summarizes the available literature on the molecular basis of the role of SWEET genes (as S genes) in the host and corresponding transcription activator-like effectors (TALEs) secreted by the pathogen. The review has four major sections, which follow a brief introduction: The first part gives some details about the occurrence and evolution of SWEET genes in approximately 30 plant species; the second part gives some details about systems where (a) SWEET genes with and without TALEs and (b) TALEs without SWEET genes cause different diseases; the third part summarizes the available information about TALEs along with interfering/truncated TALEs secreted by the pathogens; this section also summarizes the available information on effector-binding elements (EBEs) available in the promoters of either the SWEET genes or the Executor R genes; the code that is used for binding of TALEs to EBEs is also described in this section; the fourth part gives some details about the available approaches that are being used or can be used in the future for exploiting SWEET genes for developing disease-resistant cultivars. The review concludes with a section giving conclusions and future possibilities of using SWEET genes for developing disease-resistant cultivars using different approaches, including conventional breeding and genome editing.

The GhSWEET42 Glucose Transporter Participates in Verticillium dahliae Infection in Cotton

The SWEET (sugars will eventually be exported transporter) proteins, a family of sugar transporters, mediate sugar diffusion across cell membranes. Pathogenic fungi can acquire sugars from plant cells to satisfy their nutritional demands for growth and infection by exploiting plant SWEET sugar transporters. However, the mechanism underlying the sugar allocation in cotton plants infected by Verticillium dahliae, the causative agent of Verticillium wilt, remains unclear. In this study, observations of the colonization of cotton roots by V. dahliae revealed that a large number of conidia had germinated at 48-hour post-inoculation (hpi) and massive hyphae had appeared at 96 hpi. The glucose content in the infected roots was significantly increased at 48 hpi. On the basis of an evolutionary analysis, an association analysis, and qRT-PCR assays, GhSWEET42 was found to be closely associated with V. dahliae infection in cotton. Furthermore, GhSWEET42 was shown to encode a glucose transporter localized to the plasma membrane. The overexpression of GhSWEET42 in Arabidopsis thaliana plants led to increased glucose content, and compromised their resistance to V. dahliae. In contrast, knockdown of GhSWEET42 expression in cotton plants by virus-induced gene silencing (VIGS) led to a decrease in glucose content, and enhanced their resistance to V. dahliae. Together, these results suggest that GhSWEET42 plays a key role in V. dahliae infection in cotton through glucose translocation, and that manipulation of GhSWEET42 expression to control the glucose level at the infected site is a useful method for inhibiting V. dahliae infection.

Gene tagging via CRISPR-mediated homology-directed repair in cassava

Research on a few model plant-pathogen systems has benefitted from years of tool and resource development. This is not the case for the vast majority of economically and nutritionally important plants, creating a crop improvement bottleneck. Cassava bacterial blight (CBB), caused by Xanthomonas axonopodis pv. manihotis (Xam), is an important disease in all regions where cassava (Manihot esculenta Crantz) is grown. Here, we describe the development of cassava that can be used to visualize one of the initial steps of CBB infection in vivo. Using CRISPR-mediated homology-directed repair (HDR), we generated plants containing scarless insertion of GFP at the 3' end of CBB susceptibility (S) gene MeSWEET10a. Activation of MeSWEET10a-GFP by the transcription activator-like (TAL) effector TAL20 was subsequently visualized at transcriptional and translational levels. To our knowledge, this is the first such demonstration of HDR via gene editing in cassava.

Evaluation and genome-wide association study of resistance to bacterial blight race 18 in U.S. Upland cotton germplasm

Bacterial blight (BB), caused by Xanthomonas citri pv. malvacearum (Xcm), is a destructive disease to cotton production in many countries. In the U.S., Xcm race 18 is the most virulent and widespread race and can cause serious yield losses. Planting BB-resistant cotton cultivars is the most effective method of controlling this disease. In this study, 335 U.S. Upland cotton accessions were evaluated for resistance to race 18 using artificial inoculations by scratching cotyledons on an individual plant basis in a greenhouse. The analysis of variance detected significant genotypic variation in disease incidence, and 50 accessions were resistant including 38 lines with no symptoms on either cotyledons or true leaves. Many of the resistant lines were developed in the MAR (multi-adversity resistance) breeding program at Texas A&M University, whereas others were developed before race 18 was first reported in the U.S. in 1973, suggesting a broad base of resistance to race 18. A genome-wide association study (GWAS) based on 26,301 single nucleotide polymorphic (SNP) markers detected 11 quantitative trait loci (QTL) anchored by 79 SNPs, including three QTL on each of the three chromosomes A01, A05 and D02, and one QTL on each of D08 and D10. This study has identified a set of obsolete Upland germplasm with resistance to race 18 and specific chromosomal regions delineated by SNPs for resistance. The results will assist in breeding cotton for BB resistance and facilitate further genomic studies in fine mapping resistance genes to enhance the understanding of the genetic basis of BB resistance in cotton.

Escalation in the host-pathogen arms race: A host resistance response corresponds to a heightened bacterial virulence response

The zig-zag model of host-pathogen interaction describes the relative strength of defense response across a spectrum of pathogen-induced plant phenotypes. A stronger defense response results in increased resistance. Here, we investigate the strength of pathogen virulence during disease and place these findings in the context of the zig-zag model. Xanthomonas vasicola pv. holcicola (Xvh) causes sorghum bacterial leaf streak. Despite being widespread, this disease has not been described in detail at the molecular level. We divided diverse sorghum genotypes into three groups based on disease symptoms: water-soaked lesions, red lesions, and resistance. Bacterial growth assays confirmed that these three phenotypes represent a range of resistance and susceptibility. To simultaneously reveal defense and virulence responses across the spectrum of disease phenotypes, we performed dual RNA-seq on Xvh-infected sorghum. Consistent with the zig-zag model, the expression of plant defense-related genes was strongest in the resistance interaction. Surprisingly, bacterial virulence genes related to the type III secretion system (T3SS) and type III effectors (T3Es) were also most highly expressed in the resistance interaction. This expression pattern was observed at multiple time points within the sorghum-Xvh pathosystem. Further, a similar expression pattern was observed in Arabidopsis infected with Pseudomonas syringae for effector-triggered immunity via AvrRps4 but not AvrRpt2. Specific metabolites were able to repress the Xvh virulence response in vitro and in planta suggesting a possible signaling mechanism. Taken together, these findings reveal multiple permutations of the continually evolving host-pathogen arms race from the perspective of host defense and pathogen virulence responses.

The arms race between plants and pathogens is a complex process. To dissect the plant defense and pathogen virulence responses simultaneously, we used sorghum and Xanthomonas vasicola pv. holcicola, as a model pathosystem. We performed dual RNA-seq on infected sorghum with a range of disease phenotypes. Our characterization of this pathosystem demonstrates that genes related to the plant defense and pathogen virulence responses are most highly induced during a resistance interaction. We observed a similar pattern of escalation in Arabidopsis infected with Pseudomonas syringae. These observations support a conceptual model of fluidity between the different stages of plant immunity and pathogen virulence.

Mechanistic insights into host adaptation, virulence and epidemiology of the phytopathogen Xanthomonas

Xanthomonas is a well-studied genus of bacterial plant pathogens whose members cause a variety of diseases in economically important crops worldwide. Genomic and functional studies of these phytopathogens have provided significant understanding of microbial-host interactions, bacterial virulence and host adaptation mechanisms including microbial ecology and epidemiology. In addition, several strains of Xanthomonas are important as producers of the extracellular polysaccharide, xanthan, used in the food and pharmaceutical industries. This polymer has also been implicated in several phases of the bacterial disease cycle. In this review, we summarise the current knowledge on the infection strategies and regulatory networks controlling virulence and adaptation mechanisms from Xanthomonas species and discuss the novel opportunities that this body of work has provided for disease control and plant health.

Unravelling Cotton Nonexpressor of Pathogenesis-Related 1(NPR1)-Like Genes Family: Evolutionary Analysis and Putative Role in Fiber Development and Defense Pathway

The nonexpressor of pathogenesis-related 1 (NPR1) family plays diverse roles in gene regulation in the defense and development signaling pathways in plants. Less evidence is available regarding the significance of the NPR1-like gene family in cotton (Gossypium species). Therefore, to address the importance of the cotton NPR1-like gene family in the defense pathway, four Gossypium species were studied: two tetraploid species, G.hirsutum and G. barbadense, and their two potential ancestral diploids, G. raimondii and G. arboreum. In this study, 12 NPR1-like family genes in G. hirsutum were recognized, including six genes in the A-subgenome and six genes in the D-subgenome. Based on the phylogenetic analysis, gene and protein structural features, cotton NPR-like proteins were grouped into three different clades. Our analysis suggests the significance of cis-regulatory elements in the upstream region of cotton NPR1-like genes in hormonal signaling, biotic stress conditions, and developmental processes. The quantitative expression analysis for different developmental tissues and fiber stages (0 to 25 days post-anthesis), as well as salicylic acid induction, confirmed the distinct function of different cotton NPR genes in defense and fiber development. Altogether, this study presents specifications of conservation in the cotton NPR1-like gene family and their functional divergence for development of fiber and defense properties.

Genomic diversifications of five Gossypium allopolyploid species and their impact on cotton improvement

Polyploidy is an evolutionary innovation for many animals and all flowering plants, but its impact on selection and domestication remains elusive. Here we analyze genome evolution and diversification for all five allopolyploid cotton species, including economically important Upland and Pima cottons. Although these polyploid genomes are conserved in gene content and synteny, they have diversified by subgenomic transposon exchanges that equilibrate genome size, evolutionary rate heterogeneities and positive selection between homoeologs within and among lineages. These differential evolutionary trajectories are accompanied by gene-family diversification and homoeolog expression divergence among polyploid lineages. Selection and domestication drive parallel gene expression similarities in fibers of two cultivated cottons, involving coexpression networks and N⁶-methyladenosine RNA modifications. Furthermore, polyploidy induces recombination suppression, which correlates with altered epigenetic landscapes and can be overcome by wild introgression. These genomic insights will empower efforts to manipulate genetic recombination and modify epigenetic landscapes and target genes for crop improvement.

Sequencing and genomic diversification of five allopolyploid cotton species provide insights into polyploid genome evolution and epigenetic landscapes for cotton improvement.

Identification of a virulence tal gene in the cotton pathogen, Xanthomonas citri pv. malvacearum strain Xss-V2-18

Background

Bacterial blight of cotton (BBC), which is caused by the bacterium Xanthomonas citri pv. malvacearum (Xcm), is a destructive disease in cotton. Transcription activator-like effectors (TALEs), encoded by tal-genes, play critical roles in the pathogenesis of xanthomonads. Characterized strains of cotton pathogenic Xcm harbor 8–12 different tal genes and only one of them is functionally decoded. Further identification of novel tal genes in Xcm strains with virulence contributions are prerequisite to decipher the Xcm-cotton interactions.

Results

In this study, we identified six tal genes in Xss-V₂–18, a highly-virulent strain of Xcm from China, and assessed their role in BBC. RFLP-based Southern hybridization assays indicated that Xss-V₂–18 harbors the six tal genes on a plasmid. The plasmid-encoded tal genes were isolated by cloning BamHI fragments and screening clones by colony hybridization. The tal genes were sequenced by inserting a Tn5 transposon in the DNA encoding the central repeat region (CRR) of each tal gene. Xcm TALome evolutionary relationship based on TALEs CRR revealed relatedness of Xss-V₂–18 to MSCT1 and MS14003 from the United States. However, Tal2 of Xss-V₂–18 differs at two repeat variable diresidues (RVDs) from Tal6 and Tal26 in MSCT1 and MS14003, respectively, inferred functional dissimilarity. The suicide vector pKMS1 was then used to construct tal deletion mutants in Xcm Xss-V₂–18. The mutants were evaluated for pathogenicity in cotton based on symptomology and growth in planta. Four mutants showed attenuated virulence and all contained mutations in tal2. One tal2 mutant designated M2 was further investigated in complementation assays. When tal2 was introduced into Xcm M2 and expressed in trans, the mutant was complemented for both symptoms and growth in planta, thus indicating that tal2 functions as a virulence factor in Xcm Xss-V₂–18.

Conclusions

Overall, the results demonstrated that Tal2 is a major pathogenicity factor in Xcm strain Xss-V₂–18 that contributes significantly in BBC. This study provides a foundation for future efforts aimed at identifying susceptibility genes in cotton that are targeted by Tal2.

Engineering Disease-Resistant Cassava

Manihot esculenta Crantz (cassava) is a food crop originating from South America grown primarily for its starchy storage roots. Today, cassava is grown in the tropics of South America, Africa, and Asia with an estimated 800 million people relying on it as a staple source of calories. In parts of sub-Saharan Africa, cassava is particularly crucial for food security. Cassava root starch also has use in the pharmaceutical, textile, paper, and biofuel industries. Cassava has seen strong demand since 2000 and production has increased consistently year-over-year, but potential yields are hampered by susceptibility to biotic and abiotic stresses. In particular, bacterial and viral diseases can cause severe yield losses. Of note are cassava bacterial blight (CBB), cassava mosaic disease (CMD), and cassava brown streak disease (CBSD), all of which can cause catastrophic losses for growers. In this article, we provide an overview of the major microbial diseases of cassava, discuss current and potential future efforts to engineer new sources of resistance, and conclude with a discussion of the regulatory hurdles that face biotechnology.

Return of old foes - recurrence of bacterial blight and Fusarium wilt of cotton

Bacterial blight of cotton, caused by Xanthomonas citri subsp. malvacearum, and Fusarium wilt of cotton, caused by Fusarium oxysporum f. sp. vasinfectum, contribute cotton losses worldwide. Resurgences of these diseases in the United States were reported in recent years. There is a pressing need to understand pathogenicity and host responses to the pathogens and develop effective strategies for disease prevention and management. Here, we discuss the current status of bacterial blight and Fusarium wilt of cotton in the field as well as the knowledge of cotton resistance and susceptibility to these pathogens. In addition, we aim to provide insights into how these diseases are recurring and possible methods to use current technologies for biological control of these pathogens.

Development of a qPCR Protocol to Detect the Cotton Bacterial Blight Pathogen, Xanthomonas citri pv. malvacearum, from Cotton Leaves and Seeds

Bacterial blight, historically a seed-borne disease of cotton (Gossypium hirsutum) is caused by Xanthomonas citri pv. malvacearum, resulted in significant economic losses prior to development of resistant varieties and implementation of acid-delinting of planting seed. Periodic outbreaks have been associated with seed since the early twentieth century; of note, the disease has experienced a resurgence since 2011. Effective management of bacterial blight is dependent on accurate diagnosis and detection of the pathogen. Currently, detection of X. citri pv. malvacearum is performed by time-consuming microbiological methods. In this study, a novel and sensitive TaqMan-based qPCR protocol was developed to test for X. citri pv. malvacearum in cotton plant tissue. The primers developed are specific to five races of X. citri pv. malvacearum, but not to other Xanthomonas species or cotton-associated nonpathogenic bacteria. The efficiency of this assay was evaluated on artificially inoculated cotton leaves and seed, on naturally infected cotton leaves, and on bolls and seed originating from bacterial blight symptomatic bolls. The protocol's efficiency from artificially inoculated plant tissue was 10² copies g^-1 and 37 copies from 1 g seed for leaves and seed, respectively. In addition, X. citri pv. malvacearum was detected from 94% of the seed samples originating from blight symptomatic bolls. The qPCR protocol provides a rapid and accurate method for diagnosis and detection of bacterial blight and offers a tool for monitoring X. citri pv. malvacearum and potentially reducing its spread in seed.

Evolution and Stress Responses of Gossypium hirsutum SWEET Genes

The SWEET (sugars will eventually be exported transporters) proteins are sugar efflux transporters containing the MtN3_saliva domain, which affects plant development as well as responses to biotic and abiotic stresses. These proteins have not been functionally characterized in the tetraploid cotton, Gossypium hirsutum, which is a widely cultivated cotton species. In this study, we comprehensively analyzed the cotton SWEET gene family. A total of 55 putative G. hirsutumSWEET genes were identified. The GhSWEET genes were classified into four clades based on a phylogenetic analysis and on the examination of gene structural features. Moreover, chromosomal localization and an analysis of homologous genes in Gossypium arboreum, Gossypium raimondii, and G. hirsutum suggested that a whole-genome duplication, several tandem duplications, and a polyploidy event contributed to the expansion of the cotton SWEET gene family, especially in Clade III and IV. Analyses of cis-acting regulatory elements in the promoter regions, expression profiles, and artificial selection revealed that the GhSWEET genes were likely involved in cotton developmental processes and responses to diverse stresses. These findings may clarify the evolution of G. hirsutum SWEET gene family and may provide a foundation for future functional studies of SWEET proteins regarding cotton development and responses to abiotic stresses.