Revealing the Complete Bispecific Phosphatase Genes (DUSPs) across the Genome and Investigating the Expression Patterns of GH_A11G3500 Resistance against Verticillium wilt.

DUSPs, a diverse group of protein phosphatases, play a pivotal role in orchestrating cellular growth and development through intricate signaling pathways. Notably, they actively participate in the MAPK pathway, which governs crucial aspects of plant physiology, including growth regulation, disease resistance, pest resistance, and stress response. DUSP is a key enzyme, and it is the enzyme that limits the rate of cell metabolism. At present, complete understanding of the DUSP gene family in cotton and its specific roles in resistance to Verticillium wilt (VW) remains elusive. To address this knowledge gap, we conducted a comprehensive identification and analysis of four key cotton species: Gossypium arboreum, Gossypium barbadense, Gossypium hirsutum, and Gossypium raimondii. The results revealed the identification of a total of 120 DUSP genes in the four cotton varieties, which were categorized into six subgroups and randomly distributed at both ends of 26 chromosomes, predominantly localized within the nucleus. Our analysis demonstrated that closely related DUSP genes exhibited similarities in terms of the conserved motif composition and gene structure. A promoter analysis performed on the GhDUSP gene promoter revealed the presence of several cis-acting elements, which are associated with abiotic and biotic stress responses, as well as hormone signaling. A tissue expression pattern analysis demonstrated significant variations in GhDUSP gene expression under different stress conditions, with roots exhibiting the highest levels, followed by stems and leaves. In terms of tissue-specific detection, petals, leaves, stems, stamens, and receptacles exhibited higher expression levels of the GhDUSP gene. The gene expression analysis results for GhDUSPs under stress suggest that DUSP genes may have a crucial role in the cotton response to stress in cotton. Through Virus-Induced Gene Silencing (VIGS) experiments, the silencing of the target gene significantly reduced the resistance efficiency of disease-resistant varieties against Verticillium wilt (VW). Consequently, we conclude that GH_A11G3500-mediated bispecific phosphorylated genes may serve as key regulators in the resistance of G. hirsutum to Verticillium wilt (VW). This study presents a comprehensive structure designed to provide an in-depth understanding of the potential biological functions of cotton, providing a strong foundation for further research into molecular breeding and resistance to plant pathogens.

Rapidly mining candidate cotton drought resistance genes based on key indicators of drought resistance.

BACKGROUND: Focusing on key indicators of drought resistance is highly important for quickly mining candidate genes related to drought resistance in cotton. RESULTS: In the present study, drought resistance was identified in drought resistance-related RIL populations during the flowering and boll stages, and multiple traits were evaluated; these traits included three key indicators: plant height (PH), single boll weight (SBW) and transpiration rate (Tr). Based on these three key indicators, three groups of extreme mixing pools were constructed for BSA-seq. Based on the mapping interval of each trait, a total of 6.27 Mb QTL intervals were selected on chromosomes A13 (3.2 Mb), A10 (2.45 Mb) and A07 (0.62 Mb) as the focus of this study. Based on the annotation information and qRT‒PCR analysis, three key genes that may be involved in the drought stress response of cotton were screened: GhF6'H1, Gh3AT1 and GhPER55. qRT‒PCR analysis of parental and extreme germplasm materials revealed that the expression of these genes changed significantly under drought stress. Cotton VIGS experiments verified the important impact of key genes on cotton drought resistance. CONCLUSIONS: This study focused on the key indicators of drought resistance, laying the foundation for the rapid mining of drought-resistant candidate genes in cotton and providing genetic resources for directed molecular breeding of drought resistance in cotton.

Genome-Wide and Expression Pattern Analysis of the DVL Gene Family Reveals GhM_A05G1032 Is Involved in Fuzz Development in G. hirsutum.

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

Glyceraldehyde-3-phosphate dehydrogenase Gh_GAPDH9 is associated with drought resistance in Gossypium hirsutum.

Background: Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is the central enzyme of glycolysis and plays important regulatory roles in plant growth and development and responses to adverse stress conditions. However, studies on the characteristics and functions of cotton GAPDH family genes are still lacking. Methods: In this study, genome-wide identification of the cotton GAPDH gene family was performed, and the phylogeny, gene structures, promoter progenitors and expression profiles of upland cotton GAPDH gene family members were explored by bioinformatics analysis to highlight potential functions. The functions of GhGAPDH9 in response to drought stress were initially validated based on RNA-seq, qRT‒PCR, VIGS techniques and overexpression laying a foundation for further studies on the functions of GAPDH genes. Results: This study is the first systematic analysis of the cotton GAPDH gene family, which contains a total of 84 GAPDH genes, among which upland cotton contains 27 members. Quantitative, phylogenetic and covariance analyses of the genes revealed that the GAPDH gene family has been conserved during the evolution of cotton. Promoter analysis revealed that most cis-acting elements were related to MeJA and ABA. Based on the identified promoter cis-acting elements and RNA-seq data, it was hypothesized that Gh_GAPDH9, Gh_GAPDH11, Gh_GAPDH19 and Gh_GAPDH21 are involved in the response of cotton to abiotic stress. The expression levels of the Gh_GAPDH9 gene in two drought-resistant and two drought-sensitive materials were analyzed by qRT‒PCR and found to be high early in the treatment period in the drought-resistant material. The silencing of Gh_GAPDH9 based on virus-induced gene silencing (VIGS) technology resulted in significant leaf wilting or whole-plant dieback in silenced plants after drought stress compared to the control. The content of-malondialdehyde (MDA) in cotton leaves was significantly increased, and the content of proline (Pro) and chlorophyll (Chl) was reduced. In addition, the leaf wilting and dryness of transgenic lines under drought stress were lower than those of wild-type Arabidopsis, indicating that Gh_GAPDH9 is a positive regulator of drought resistance. In conclusion, our results demonstrate that GAPDH genes play an important role in the response of cotton to abiotic stresses and provide preliminary validation of the function of the Gh_GAPDH9 gene under drought stress. These findings provide an important theoretical basis for further studies on the function of the Gh_GAPDH9 gene and the molecular mechanism of the drought response in cotton.

Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L.

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

GBSOT4 Enhances the Resistance of Gossypium barbadense to Fusarium oxysporum f. sp. vasinfectum (FOV) by Regulating the Content of Flavonoid.

Sulfotransferases (SOTs) (EC 2.8.2.-) are sulfate regulatory proteins in a variety of organisms that have been previously shown to be involved in regulating a variety of physiological and biological processes, such as growth, development, adaptation to land, stomatal closure, drought tolerance, and response to pathogen infection. However, there is a lack of comprehensive identification and systematic analysis of SOT in cotton, especially in G. barbadense. In this study, we used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structure, expression patterns, evolutionary relationships, selection pressure and stress response of SOT gene family members in G. barbadense. In this study, a total of 241 SOT genes were identified in four cotton species, among which 74 SOT gene members were found in G. barbadense. According to the phylogenetic tree, 241 SOT protein sequences were divided into five distinct subfamilies. We also mapped the physical locations of these genes on chromosomes and visualized the structural information of SOT genes in G. barbadense. We also predicted the cis-acting elements of the SOT gene in G. barbadense, and we identified the repetitive types and collinearity analysis of SOT genes in four cotton species. We calculated the Ka/Ks ratio between homologous gene pairs to elucidate the selective pressure between SOT genes. Transcriptome data were used to explore the expression patterns of SOT genes, and then qRT-PCR was used to detect the expression patterns of GBSOT4, GBSOT17 and GBSOT33 under FOV stress. WGCNA (weighted gene co-expression network analysis) showed that GB_A01G0479 (GBSOT4) belonged to the MEblue module, which may regulate the resistance mechanism of G. barbadense to FOV through plant hormones, signal transduction and glutathione metabolism. In addition, we conducted a VIGS (virus-induced gene silencing) experiment on GBSOT4, and the results showed that after FOV inoculation, the plants with a silenced target gene had more serious leaf wilting, drying and cracking than the control group, and the disease index of the plants with the silenced target gene was significantly higher than that of the control group. This suggests that GBSOT4 may be involved in protecting the production of G. barbadense from FOV infection. Subsequent metabolomics analysis showed that some flavonoid metabolites, such as Eupatorin-5-methylether (3'-hydroxy-5,6,7,4'-tetramethoxyflavone, were accumulated in cotton plants in response to FOV infection.

Using association mapping and local interval haplotype association analysis to improve the cotton drought stress response.

Drought stress has a serious impact on the growth and development of cotton. To explore the relevant molecular mechanism of the drought stress response in cotton, gene mapping based on the QTL interval mapped by simplified genome BSA-seq of the drought-resistance-related RIL population was performed. A QTL region spanning 2.02 Mb on chromosome D07 was selected, and 201 resource materials were genotyped using 9 KASP markers in the interval. After local interval haplotype association analysis, the overlap of the 110 kb peak region confirmed the reliability of this region, and at the same time, the role of GhGF14-30, the only gene in the overlapping region, was modeled in the response of cotton to drought stress. qRTPCR analysis of the materials and population parents proved that this gene plays a role in the drought stress response in cotton. Virus-induced gene silencing proved the importance of this gene in drought-sensitive materials, and drought-resistance-related marker genes also proved that the GhGF14-30 gene may play an important role in the ABA and SOS signaling pathways. This study provides a basis for mining drought stress response functional genes in cotton and lays the foundation for the molecular mechanism of the GhGF14-30 gene in response to drought stress in cotton.

Development and identification of molecular markers of GhHSP70-26 related to heat tolerance in cotton.

Heat stress significantly affect plant growth and development, which is an important factor contributing to crop yield loss. However, heat shock proteins (HSPs) in plants can effectively alleviate cell damage caused by heat stress. In order to rapidly and accurately cultivate heat-tolerant cotton varieties, this study conducted correlation analysis between heat tolerance index and insertion/deletion (In/Del) sites of GhHSP70-26 promoter in 39 cotton materials, so as to find markers related to heat tolerance function of cotton, which can be used in molecular marker-assisted breeding. The results showed the natural variation allele (Del22 bp) type at -1590 bp upstream of GhHSP70-26 promoter (haplotype2, Hap2) in cotton (Gossypium spp.) promoted GhHSP70-26 expression under heat stress. The relative expression level of GhHSP70-26 of M-1590-Del22 cotton materials were significantly higher than that of M-1590-In type cotton materials under heat stress (40 ℃). Also, M-1590-Del22 material had lower conductivity and less cell damage after heat stress, indicating that it is a heat resistant cotton material. The Hap1 (M-1590-In) promoter was mutated into Hap1del22, and Hap1 and Hap1del22 were fused with GUS to transform Arabidopsis thaliana. Furthermore, Hap1del22 promoter had higher induction activity than Hap1 under heat stress and abscisic acid (ABA) treatment in transgenic Arabidopsis thaliana. Further analysis confirmed that M-1590-Del22 was the dominant heat-resistant allele. In summary, these results identify a key and previously unknown natural variation in GhHSP70-26 with respect to heat tolerance, providing a valuable functional molecular marker for genetic breeding of cotton and other crops with heat tolerance.

Yield-based drought tolerance index evaluates the drought tolerance of cotton germplasm lines in the interaction of genotype-by-environment.

Background: Cotton is an economically important crop in China, and drought has seriously affected cotton production. Understanding genetic variation, genotype ×environment interactions, and the associations between these traits is critical for developing improved cotton varieties with high drought tolerance. Methods: To screen ideal drought-resistant cotton germplasm lines and excellent genotypes, the yield traits of 103 cotton germplasm lines were analyzed. Cotton resource material was planted under normal watering and water deficit conditions for three consecutive years. The yield traits under normal irrigation and water stress conditions were measured, and then five screening indicators were calculated based on the cotton yield per plant under the two water treatments to determine the ideal genotype and most accurate identification indicators. Results: The results of correlation analysis and principal component analysis showed that the geometric mean productivity (GMP), mean productivity (MP), and stress tolerance index (STI) were significantly positively correlated with yield under water stress and could be used to distinguish genotypes with high drought tolerance. Among the experimental germplasm lines, some had higher STI and GMP values, indicating their higher drought tolerance. This result indicates that best linear unbiased prediction (BLUP) analysis of the STI and GMP under drought stress can effectively improve screening for drought tolerance in cotton germplasm lines. The results from the screening index, three-dimensional map, and genotype ×environment (GGE) biplots were consistent with the above results. We determined that CQJ-5, Xin lu zao 45, Bellsno, Zhong R 2016 and ND 359-5 are drought-tolerant genotypes that can be used to breed drought-tolerant germplasm lines that produce high and stable yields.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes for Fuzz Development in Gossypium hirsutum.

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

Quantitative trait locus mapping and identification of candidate genes for resistance to Verticillium wilt in four recombinant inbred line populations of Gossypium hirsutum.

Improving resistance to Verticillium wilt is of great significance for achieving high and stable yields of Upland cotton (Gossypium hirsutum). To deeply understand the genetic basis of cotton resistance to Verticillium wilt, Verticillium wilt-resistant Upland Lumianyan 28 and four Verticillium wilt-susceptible Acala cotton cultivars were used to create four recombinant inbred line (RIL) populations of 469 families through nested hybridization. Phenotypic data collected in five stressful environments were used to select resistant and sensitive lines and create a mixed pool of extreme phenotypes for BSA-seq. A total of 8 QTLs associated with Verticillium wilt resistance were identified on 4 chromosomes, of which qVW-A12-5 was detected simultaneously in the RIL populations and in one of the RIL populations and was identified for the first time. According to the sequence comparison and transcriptome analysis of candidate genes in the QTL interval between parents and pools, 4 genes were identified in the qVW-A12-5 interval. qRT-PCR of parental and phenotypically extreme lines revealed that Gh_CPR30 was induced by and may be a candidate gene for resistance to Verticillium wilt in G. hirsutum. Furthermore, VIGS technology revealed that the disease severity index (DSI) of the Gh_CPR30-silenced plants was significantly higher than that of the control. These results indicate that the Gh_CPR30 gene plays an important role in the resistance of G. hirsutum to Verticillium wilt, and the study provides a molecular basis for analyzing the molecular mechanism underlying G. hirsutum resistance to Verticillium wilt.

Efficient virus-mediated genome editing in cotton using the CRISPR/Cas9 system.

Plant virus-mediated sgRNA delivery and expression have great advantages; sgRNA expression can rapidly expand and accumulate along with virus replication and movement, resulting in efficient gene editing efficiency. In this study, a VIGE system based on cotton leaf crumple virus (CLCrV) was established using cotton overexpressing Cas9 (Cas9-OE) as the VIGE receptor. CLCrV-mediated VIGE could not only target and knock out the GhMAPKKK2, GhCLA1 and GhPDS genes subgroup A and D genome sequences but also achieve double mutation of GhCLA1 and GhPDS genes at the same time. These results verified the effectiveness and efficiency of this system. In addition, the off-target effect assay demonstrated that the CLCrV-mediated VIGE system not only has high gene editing efficiency but also high gene editing specificity in cotton. We further explored whether the FT-sgRNA strategy could transport sgRNA to cotton apical meristem (SAM) over long distances to avoid using tissue culture to obtain stable genetic mutants. The results showed that the sgRNA fused with FT mRNA at the 5' end could also efficiently achieve targeted editing of endogenous genes in cotton, but it was difficult to detect heritable mutant progeny. The above results showed that the CLCrV-mediated VIGE system provided an accurate and rapid validation tool for screening effective sgRNAs in cotton.

Protein S-Acyl Transferase GhPAT27 Was Associated with Verticillium wilt Resistance in Cotton.

Protein palmitoylation is an ability of the frame of the cell marker protein is one of the most notable reversible changes after translation. However, studies on protein palmitoylation in cotton have not yet been performed. In our current research, the PAT gene family was systematically identified and bioinformatically analyzed in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, and 211 PAT genes were authenticated and classified into six subfamilies. Sixty-nine PAT genes were identified in upland cotton, mainly at the ends of its the 26 chromosomes of upland cotton. The majority of these genes are located in the nucleus of the plant. Gene structure analysis revealed that each member encodes a protein that which contains at least one DHHC structural domain. Cis-acting element analysis indicated that GhPATs genes are mainly involved in hormone production, light response and stress response. Gene expression pattern analysis indicated that most GhPATs genes were differentially expressed upon induction by pathogenic bacteria, drought, salt, hot and cold stresses, and some GhPATs could be induced by multiple abiotic stresses simultaneously. GhPATs genes showed different expression patterns in tissue-specific assays and were found to be preferentially expressed in roots, followed by expression in stems and leaves. Virus-induced gene silencing (VIGS) experiments showed that cotton was significantly less resistant to Verticillium dahliae when GhPAT27 was silenced. We conclude that the GhPAT27 gene, which mediates S-palmitoylation acetylation, may be involved in the regulation of upland cotton resistance to Verticillium wilt (VW). Overall, this work has provided a fundamental framework for understanding the latent capabilities of GhPATs and a solid foundation for molecular breeding and plant pathogen resistance in cotton.

Identification and functional analysis of PIN family genes in Gossypium barbadense.

Background: PIN proteins are an important class of auxin polar transport proteins that play an important regulatory role in plant growth and development. However, their characteristics and functions have not been identified in Gossypium barbadense. Methods: PIN family genes were identified in the cotton species G. barbadense, Gossypium hirsutum, Gossypium raimondii, and Gossypium arboreum, and detailed bioinformatics analyses were conducted to explore the roles of these genes in G. barbadense using transcriptome data and quantitative reverse-transcription polymerase chain reaction (qRT-PCR) technology. Functional verification of the genes was performed using virus-induced gene silencing (VIGS) technology. Results: A total of 138 PIN family genes were identified in the four cotton species; the genes were divided into seven subgroups. GbPIN gene family members were widely distributed on 20 different chromosomes, and most had repeated duplication events. Transcriptome analysis showed that some genes had differential expression patterns in different stages of fiber development. According to 'PimaS-7' and '5917' transcript component association analysis, the transcription of five genes was directly related to endogenous auxin content in cotton fibers. qRT-PCR analysis showed that the GbPIN7 gene was routinely expressed during fiber development, and there were significant differences among materials. Transient silencing of the GbPIN7 gene by VIGS led to significantly higher cotton plant growth rates and significantly lower endogenous auxin content in leaves and stems. This study provides comprehensive analyses of the roles of PIN family genes in G. barbadense and their expression during cotton fiber development. Our results will form a basis for further PIN auxin transporter research.

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

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

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

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

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

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

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton.

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

Quantitative Trait Locus Mapping and Identification of Candidate Genes for Resistance to Fusarium Wilt Race 7 Using a Resequencing-Based High Density Genetic Bin Map in a Recombinant Inbred Line Population of Gossypium barbadense.

Fusarium wilt caused by Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive diseases in cotton (Gossypium spp.) production, and use of resistant cultivars is the most cost-effective method managing the disease. To understand the genetic basis of cotton resistance to FOV race 7 (FOV7), this study evaluated a recombinant inbred line (RIL) population of 110 lines of G. barbadense from a cross between susceptible Xinhai 14 and resistant 06-146 in eight tests and constructed a high-density genetic linkage map with resequencing-based 933,845 single-nucleotide polymorphism (SNP) markers covering a total genetic distance of 2483.17 cM. Nine quantitative trait loci (QTLs) for FOV7 resistance were identified, including qFOV7-D03-1 on chromosome D03 in two tests. Through a comparative analysis of gene expression and DNA sequence for predicted genes within the QTL region between the two parents and selected lines inoculated with FOV7, GB_D03G0217 encoding for a calmodulin (CaM)-like (CML) protein was identified as a candidate gene. A further analysis confirmed that the expression of GB_D03G0217 was suppressed, leading to increased disease severity in plants of the resistant parent with virus induced gene silencing (VIGS).

Nanopore-Based Comparative Transcriptome Analysis Reveals the Potential Mechanism of High-Temperature Tolerance in Cotton (Gossypium hirsutum L.).

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant ('Xinluzao36') and heat-sensitive ('Che61-72') cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.