ACO is one of the rate-limiting enzymes in the biosynthesis of ethylene, and it plays a critical role in the regulation of plant growth and development. However, the function of ACO genes in cotton is not well studied. In this study, a total of 332 GhACOs, 187 GaACOs, and 181 GrACOs were identified in G. hirsutum, G. arboretum, and G. raimondii, respectively. Gene duplication analysis showed that whole-genome duplication (WGD) and tandem duplication were the major forces driving the generation of cotton ACO genes. In the promoters of GhACOs, there were cis-acting elements responding to stress, phytohormones, light, and circadian factors, indicating the possible involvement of GhACOs in these processes. Expression and co-expression analyses illustrated that most GhACOs were not only widely expressed in various tissues but also coexpressed with other genes in response to salt and drought stress. GhACO106_At overexpression in Arabidopsis promoted flowering and increased salt tolerance. These results provide a comprehensive overview of the ACO genes of cotton and lay the foundation for subsequent functional studies of these genes.

Background: Multiple C2 domains and transmembrane region proteins (MCTPs) may act as transport mediators of other regulators. Although increased number of MCTPs in higher plants implies their diverse and specific functions in plant growth and development, only a few plant MCTPs have been studied and no study on the MCTPs in cotton has been reported. Results: In this study, we identified 31 MCTPs in G. hirsutum, which were classified into five subfamilies according to the phylogenetic analysis. GhMCTPs from subfamily V exhibited isoelectric points (pIs) less than 7, whereas GhMCTPs from subfamily I, II, III and IV exhibited pIs more than 7.5, implying their distinct biological functions. In addition, GhMCTPs within subfamily III, IV and V exhibited more diverse physicochemical properties, domain architectures and expression patterns than GhMCTPs within subfamily I and II, suggesting that GhMCTPs within subfamily III, IV and V diverged to perform more diverse and specific functions. Analyses of conserved motifs and pIs indicated that the N-terminus was more divergent than the C-terminus and GhMCTPs’ functional divergence might be mainly contributed by the N-terminus. Furthermore, yeast two-hybrid assay indicated that the N-terminus was responsible to interact with target proteins. Phylogenetic analysis classified multiple N-terminal C2 domains into four subclades, suggesting that these C2 domains performed different molecular functions in mediating the transport of target proteins. Conclusions: Our systematic characterization of MCTPs in G. hirsutum will provide helpful information to further research GhMCTPs’ molecular roles in mediating other regulators’ transport to coordinate growth and development of various cotton tissues.

Background: The jasmonate ZIM domain (JAZ) protein is a key repressor of the jasmonate signal transduction pathway, which plays an important role in plant growth and development and defense responses. In this study, based on the published whole-genome data, we identified members of the JAZ gene family in Populus trichocarpa. Through a series of bioinformatic approaches, their expression patterns under various stress conditions have been analyzed to explore and excavate the endogenous resistance genes of poplar and provide a theoretical basis for breeding new varieties of poplar resistance. Results: A total of 13 PtJAZ genes have been identified in P. trichocarpa and designated as PtJAZ1–PtJAZ13. Those 13 PtJAZ genes were unevenly distributed on nine chromosomes, and they could be divided into four subfamilies. The gene structures and motif composition of the members derived from the same subfamily were similar. Collinearity analysis demonstrated that, compared with Arabidopsis thaliana and Oryza sativa, the most collinear pairs (13) were found in P. trichocarpa and Eucalyptus robusta. Cis-acting element analysis suggested that the promoter regions of PtJAZs contained a large number of hormones and stress response elements, of which abscisic acid (ABA) and methyl jasmonate (MeJA) hormone response elements were the most abundant. The PtJAZ genes not only had diverse expression patterns in different tissues, but they also responded to various abiotic and biotic stress conditions. The co-expression network and GO and KEGG analyses showed that JAZ genes were closely related to insect resistance. Conclusions: In this study, applying bioinformatic methods, 13 PtJAZ gene family members from P. trichocarpa were identified and comprehensively analyzed. By further studying the function of the poplar JAZ gene family, the aim is to select genes with better insect resistance and stress resistance so as to lay a solid foundation for the subsequent breeding of new poplar varieties.

Soil salinization conditions seriously restrict cotton yield and quality. Related studies have shown that the DUF4228 proteins are pivotal in plant resistance to abiotic stress. However, there has been no systematic identification and analysis of the DUF4228 gene family in cotton and their role in abiotic stress. In this study, a total of 308 DUF4228 genes were identified in four Gossypium species, which were divided into five subfamilies. Gene structure and protein motifs analysis showed that the GhDUF4228 proteins were conserved in each subfamily. In addition, whole genome duplication (WGD) events and allopolyploidization might play an essential role in the expansion of the DUF4228 genes. Besides, many stress-responsive (MYB, MYC) and hormone-responsive (ABA, MeJA) related cis-elements were detected in the promoters of the DUF4228 genes. The qRT-PCR results showed that GhDUF4228 genes might be involved in the response to abiotic stress. VIGS assays and the measurement of relative water content (RWC), Proline content, POD activity, and malondialdehyde (MDA) content indicated that GhDUF4228-67 might be a positive regulator of cotton response to salt stress. The results in this study systematically characterized the DUF4228s in Gossypium species and will provide helpful information to further research the role of DUF4228s in salt tolerance.

As one of the most important factors in alternative splicing (AS) events, serine/arginine-rich (SR) proteins not only participate in the growth and development of plants but also play pivotal roles in abiotic stresses. However, the research about SR proteins in cotton is still lacking. In this study, we performed an extensive comparative analysis of SR proteins and determined their phylogeny in the plant lineage. A total of 169 SR family members were identified from four Gossypium species, and these genes could be divided into eight distinct subfamilies. The domain, motif distribution and gene structure of cotton SR proteins are conserved within each subfamily. The expansion of SR genes is mainly contributed by WGD and allopolyploidization events in cotton. The selection pressure analysis showed that all the paralogous gene pairs were under purifying selection pressure. Many cis-elements responding to abiotic stress and phytohormones were identified in the upstream sequences of the GhSR genes. Expression profiling suggested that some GhSR genes may involve in the pathways of plant resistance to abiotic stresses. The WGCNA analysis showed that GhSCL-8 co-expressed with many abiotic responding related genes in a salt-responding network. The Y2H assays showed that GhSCL-8 could interact with GhSRs in other subfamilies. The subcellular location analysis showed that GhSCL-8 is expressed in the nucleus. The further VIGS assays showed that the silencing of GhSCL-8 could decrease salt tolerance in cotton. These results expand our knowledge of the evolution of the SR gene family in plants, and they will also contribute to the elucidation of the biological functions of SR genes in the future.

Abiotic stress, such as drought and salinity stress, seriously inhibit the growth and development of plants. Therefore, it is vital to understand the drought and salinity resistance mechanisms to enable cotton to provide more production under drought and salt conditions. In this study, we identified 8806 and 9108 differentially expressed genes (DEGs) through a comprehensive analysis of transcriptomic data related to the PEG-induced osmotic and salt stress in cotton. By performing weighted gene co-expression network analysis (WGCNA), we identified four co-expression modules in PEG treatment and five co-expression modules in salinity stress, which included 346 and 324 predicted transcription factors (TFs) in these modules, respectively. Correspondingly, whole genome duplication (WGD) events mainly contribute to the expansion of those TFs. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) analyses revealed those different modules were associated with stress resistance, including regulating macromolecule metabolic process, peptidase activity, transporter activity, lipid metabolic process, and responses to stimulus. Quantitative RT-PCR analysis was used to confirm the expression levels of 15 hub TFs in PEG6000 and salinity treatments. We found that the hub gene GhWRKY46 could alter salt and PEG-induced drought resistance in cotton through the virus-induced gene silencing (VIGS) method. Our results provide a preliminary framework for further investigation of the cotton response to salt and drought stress, which is significant to breeding salt- and drought-tolerant cotton varieties.