Gossypium hirsutum L., is a widely cultivated cotton species around the world, but its production is seriously threatened by its susceptibility to chilling stress. Low temperature affects its germination, and the underlying molecular mechanisms are rarely known, particularly from a transcriptional perspective. In this study, transcriptomic profiles were analyzed and compared between two cotton varieties, the cold-tolerant variety KN27-3 and susceptible variety XLZ38. A total of 7535 differentially expressed genes (DEGs) were identified. Among them, the transcripts involved in energy metabolism were significantly enriched during germination based on analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, such as glycolysis/gluconeogenesis, tricarboxylic acid cycle (TCA cycle), and glyoxylate cycle (GAC). Results from further GO enrichment analysis show the earlier appearance of DNA integration, meristem growth, cotyledon morphogenesis, and other biological processes in KN27-3 compared with XLZ38 under chilling conditions. The synthesis of asparagine, GDP-mannose, and trehalose and the catabolic process of raffinose were activated. DEGs encoding antioxidants (spermidine) and antioxidase (CAT1, GPX4, DHAR2, and APX1) were much more up-regulated in embryos of KN27-3. The content of auxin (IAA), cis-zeatin riboside (cZR), and trans-zeatin riboside (tZR) in KN27-3 are higher than that in XLZ38 at five stages (from 12 h to 54 h). GA3 was expressed at a higher level in KN27-3 from 18 h to 54 h post imbibition compared to that in XLZ38. And abscisic acid (ABA) content of KN27-3 is lower than that in XLZ38 at five stages. Results from hormone content measurements and the related gene expression analysis indicated that IAA, CTK, and GA3 may promote germination of the cold-tolerant variety, while ABA inhibits it. These results expand the understanding of cottonseed germination and physiological regulations under chilling conditions by multiple pathways.

Background: Wall-associated kinases (WAK), one of the receptor-like kinases (RLK), function directly in the connection and communication between the plant cell wall and the cytoplasm. WAK genes are highly conserved and have been identified in plants, such as rice, but there is little research on the WAK gene family in cotton. Results: In the present study, we identified 29 GhWAK genes in Gossypium hirsutum. Phylogenetic analysis showed that cotton WAK proteins can be divided into five clades. The results of synteny and Ka/Ks analysis showed that the GhWAK genes mainly originated from whole genome duplication (WGD) and were then mainly under purifying selection. Transcriptome data and real-time PCR showed that 97% of GhWAK genes highly expressed in cotton fibers and ovules. β-glucuronidase (GUS) staining assays showed that GhWAK5 and GhWAK16 expressed in Arabidopsis leaf trichomes. Fourteen GhWAK genes were found to possess putative gibberellin (GA) response elements in the promoter regions, 13 of which were significantly induced by GA treatment. Ten GhWAK genes contained auxin (IAA) response elements and the expression level of nine GhWAKs significantly increased under auxin treatment. Conclusions: We provide a preliminary analysis of the WAK gene family in G. hirsutum, which sheds light on the potantial roles of GhWAK genes in cotton fiber cell development. Our data also provides a useful resource for future studies on the functional roles of GhWAK genes.

(1) Background: Upland cotton (Gossypium hirsutum L.) is the most important natural fiber worldwide, and it is extensively planted and plentifully used in the textile industry. Major cotton planting regions are frequently affected by abiotic stress, especially drought stress. Drought resistance is a complex, quantitative trait. A genome-wide association study (GWAS) constitutes an efficient method for dissecting the genetic architecture of complex traits. In this study, the drought resistance of a population of 316 upland cotton accessions was studied via GWAS. (2) Methods: GWAS methodology was employed to identify relationships between molecular markers or candidate genes and phenotypes of interest. (3) Results: A total of 8, 3, and 6 SNPs were associated with the euphylla wilting score (EWS), cotyledon wilting score (CWS), and leaf temperature (LT), respectively, based on a general linear model and a factored spectrally transformed linear mixed model. For these traits, 7 QTLs were found, of which 2 each were located on chromosomes A05, A11, and D03, and of which 1 was located on chromosome A01. Importantly, in the candidate regions WRKY70, GhCIPK6, SnRK2.6, and NET1A, which are involved in the response to abscisic acid (ABA), the mitogen-activated protein kinase (MAPK) signaling pathway and the calcium transduction pathway were identified in upland cotton at the seedling stage under drought stress according to annotation information and linkage disequilibrium (LD) block analysis. Moreover, RNA sequencing analysis showed that WRKY70, GhCIPK6, SnRK2.6, and NET1A were induced by drought stress, and the expression of these genes was significantly different between normal and drought stress conditions. (4) Conclusions: The present study should provide some genomic resources for drought resistance in upland cotton. Moreover, the germplasm of the different phenotypes, the detected SNPs and, the potential candidate genes will be helpful for molecular marker-assisted breeding studies about increased drought resistance in upland cotton.

Background: Appropriate plant architecture can improve the amount of cotton boll opening and allow increased planting density, thus increasing the level of cotton mechanical harvesting and cotton yields. The internodes of cotton fruiting branches are an important part of cotton plant architecture. Thus, studying the molecular mechanism of internode elongation in cotton fruiting branches is highly important.Results: In this study, we selected internodes of cotton fruiting branches at three different stages from two cultivars whose internode lengths differed significantly. A total of 76,331 genes were detected by transcriptome sequencing. By KEGG pathway analysis, we found that DEGs were significantly enriched in the plant hormone signal transduction pathway. The transcriptional data and qRT-PCR results showed that members of theGH3gene family, which are involved in auxin signal transduction, andCKXenzymes, which can reduce the level of CKs, were highly expressed in the cultivar XLZ77, which has relatively short internodes. Genes related to ethylene synthase (ACS),EIN2/3andERFin the ethylene signal transduction pathway and genes related toJAR1,COI1andMYC2in the JA signal transduction pathway were also highly expressed in XLZ77. Plant hormone determination results showed that the IAA and CK contents significantly decreased in cultivar XLZ77 compared with those in cultivar L28, while the ACC (the precursor of ethylene) and JA contents significantly increased. GO enrichment analysis revealed that the GO categories associated with promoting cell elongation, such as cell division, the cell cycle process and cell wall organization, were significantly enriched, and related genes were highly expressed in L28. However, genes related to the sphingolipid metabolic process and lignin biosynthetic process, whose expression can affect cell elongation, were highly expressed in XLZ77. In addition, 2067 TFs were differentially expressed. TheWRKY,ERFandbHLHTF families were the top three largest families whose members were active in the two varieties, and the expression levels of most of the genes encoding these TFs were upregulated in XLZ77.Conclusions: Auxin and CK are positive regulators of internode elongation in cotton branches. In contrast, ethylene and JA may act as negative regulators of internode elongation in cotton branches. Furthermore, theWRKY,ERFandbHLHTFs were identified as important inhibitors of internode elongation in cotton. In XLZ77(a short-internode variety), the mass synthesis of ethylene and amino acid conjugation of auxin led to the inhibition of plant cell elongation, while an increase in JA content and degradation of CKs led to a slow rate of cell division, which eventually resulted in a phenotype that presented relatively short internodes on the fruiting branches. The results of this study not only provide gene resources for the genetic improvement of cotton plant architecture but also lay a foundation for improved understanding of the molecular mechanism of the internode elongation of cotton branches.

Soil microbes are crucial in shaping the root-associated microbial communities. In this study, we analyzed the effect of the soil–root niche gradient on the diversity, composition, and assembly of the bacterial community and co-occurrence network of two cotton varieties. The results revealed that the bacterial communities in cotton soil–root compartment niches exhibited a skewed species abundance distribution, dominated by abundant taxa showing a strong spatial specificity. The assembly processes of the rhizosphere bacterial communities were mainly driven by stochastic processes, dominated by the enrichment pattern and supplemented by the depletion pattern to recruit bacteria from the bulk soil, resulting in a more stable bacterial community. The assembly processes of the endosphere bacterial communities were determined by processes dominated by the depletion pattern and supplemented by the enrichment pattern to recruit species from the rhizosphere, resulting in a decrease in the stability and complexity of the community co-occurrence network. The compartment niche shaped the diversity of the bacterial communities, and the cotton variety genotype was an important source of diversity in bacterial communities within the compartment niche. We suggest that the moderate taxa contribute to significantly more changes in the diversity of the bacterial community than the rare and abundant taxa during the succession of bacterial communities in the cotton root–soil continuum.