Genome-Wide Study of the GATL Gene Family in Gossypium hirsutum L. Reveals that GhGATL Genes Act on Pectin Synthesis to Regulate Plant Growth and Fiber Elongation

Daily glycome and transcriptome profiling reveals polysaccharide structures and correlated glycosyltransferases critical for cotton fiber growth

Cotton fiber is the most valuable naturally available material for the textile industry and the fiber length and strength are key determinants of its quality. Dynamic changes in the pectin, xyloglucan, xylan, and cellulose polysaccharide epitope content during fiber growth contribute to complex remodeling of fiber cell wall (CW) and quality. Detailed knowledge about polysaccharide compositional and structural alteration in the fiber during fiber elongation and strengthening is important to understand the molecular dynamics of fiber development and improve its quality. Here, large-scale glycome profiling coupled with fiber phenotype and transcriptome profiling was conducted on fiber collected daily covering the most critical window of fiber development. The profiling studies with high temporal resolution allowed us to identify specific polysaccharide epitopes associated with distinct fiber phenotypes that might contribute to fiber quality. This study revealed the critical role of highly branched RG-I pectin epitopes such as β-1,4-linked-galactans, β-1,6-linked-galactans, and arabinogalactans, in addition to earlier reported homogalacturonans and xyloglucans in the formation of cotton fiber middle lamella and contributing to fiber plasticity and elongation. We also propose the essential role of heteroxylans (Xyl-MeGlcA and Xyl-3Ar), as a guiding factor for secondary CW cellulose microfibril arrangement, thus contributing to fiber strength. Correlation analysis of profiles of polysaccharide epitopes from glycome data and expression profiles of glycosyltransferase-encoding genes from transcriptome data identified several key putative glycosyltransferases that are potentially involved in synthesizing the critical polysaccharide epitopes. The findings of this study provide a foundation to identify molecular factors that dictate important fiber traits.

Pan-genome analysis of GT64 gene family and expression response to Verticillium wilt in cotton

BACKGROUND

The GT64 subfamily, belonging to the glycosyltransferase family, plays a critical function in plant adaptation to stress conditions and the modulation of plant growth, development, and organogenesis processes. However, a comprehensive identification and systematic analysis of GT64 in cotton are still lacking.

RESULTS

This study used bioinformatics techniques to conduct a detailed investigation on the GT64 gene family members of eight cotton species for the first time. A total of 39 GT64 genes were detected, which could be classified into five subfamilies according to the phylogenetic tree. Among them, six genes were found in upland cotton. Furthermore, investigated the precise chromosomal positions of these genes and visually represented their gene structure details. Moreover, forecasted cis-regulatory elements in GhGT64s and ascertained the duplication type of the GT64 in the eight cotton species. Evaluation of the Ka/Ks ratio for similar gene pairs among the eight cotton species provided insights into the selective pressures acting on these homologous genes. Additionally, analyzed the expression profiles of the GT64 gene family. Overexpressing GhGT64_4 in tobacco improved its disease resistance. Subsequently, VIGS experiments conducted in cotton demonstrated reduced disease resistance upon silencing of the GhGT64_4, may indicate its involvement in affecting lignin and jasmonic acid biosynthesis pathways, thus impacting cotton resistance. Weighted Gene Co-expression Network Analysis (WGCNA) revealed an early immune response against Verticillium dahliae in G. barbadense compared to G. hirsutum. Quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR) analysis indicated that some GT64 genes might play a role under various biotic and abiotic stress conditions.

CONCLUSIONS

These discoveries enhance our knowledge of GT64 family members and lay the groundwork for future investigations into the disease resistance mechanisms of this gene in cotton.

Genome-wide identification and expression analysis of Ubiquitin-specific protease gene family in maize (Zea mays L.)

BACKGROUND

Ubiquitin-specific proteases (UBPs) are a large family of deubiquitinating enzymes (DUBs). They are widespread in plants and are critical for plant growth, development, and response to external stresses. However, there are few studies on the functional characteristics of the UBP gene family in the important staple crop, maize (Zea mays L.).

RESULTS

In this study, we performed a bioinformatic analysis of the entire maize genome and identified 45 UBP genes. Phylogenetic analysis indicated that 45 ZmUBP genes can be divided into 15 subfamilies. Analysis of evolutionary patterns and divergence levels indicated that ZmUBP genes were present before the isolation of dicotyledons, were highly conserved and subjected to purifying selection during evolution. Most ZmUBP genes exhibited different expression levels in different tissues and developmental stages. Based on transcriptome data and promoter element analysis, we selected eight ZmUBP genes whose promoters contained a large number of plant hormones and stress response elements and were up-regulated under different abiotic stresses for RT-qPCR analysis, results showed that these genes responded to abiotic stresses and phytohormones to varying degrees, indicating that they play important roles in plant growth and stress response.

CONCLUSIONS

In this study, the structure, location and evolutionary relationship of maize UBP gene family members were analyzed for the first time, and the ZmUBP genes that may be involved in stress response and plant growth were identified by combining promoter element analysis, transcriptome data and RT-qPCR analysis. This study informs research on the involvement of maize deubiquitination in stress response.

Physiological and molecular mechanisms of ZnO quantum dots mitigating cadmium stress in Salvia miltiorrhiza

This study delved into the physiological and molecular mechanisms underlying the mitigation of cadmium (Cd) stress in the model medicinal plant Salvia miltiorrhiza through the application of ZnO quantum dots (ZnO QDs, 3.84 nm). A pot experiment was conducted, wherein S. miltiorrhiza was subjected to Cd stress for six weeks with foliar application of 100 mg/L ZnO QDs. Physiological analyses demonstrated that compared to Cd stress alone, ZnO QDs improved biomass, reduced Cd accumulation, increased the content of photosynthetic pigments (chlorophyll and carotenoids), and enhanced the levels of essential nutrient elements (Ca, Mn, and Cu) under Cd stress. Furthermore, ZnO QDs significantly lowered Cd-induced reactive oxygen species (ROS) content, including H2O2, O2-, and MDA, while enhancing the activity of antioxidant enzymes (SOD, POD, APX, and GSH-PX). Additionally, ZnO QDs promoted the biosynthesis of primary and secondary metabolites, such as total protein, soluble sugars, terpenoids, and phenols, thereby mitigating Cd stress in S. miltiorrhiza. At the molecular level, ZnO QDs were found to activate the expression of stress signal transduction-related genes, subsequently regulating the expression of downstream target genes associated with metal transport, cell wall synthesis, and secondary metabolite synthesis via transcription factors. This activation mechanism contributed to enhancing Cd tolerance in S. miltiorrhiza. In summary, these findings shed light on the mechanisms underlying the mitigation of Cd stress by ZnO QDs, offering a potential nanomaterial-based strategy for enhancing Cd tolerance in medicinal plants.

Genome-Wide and Expression Pattern Analysis of the HIT4 Gene Family Uncovers the Involvement of GHHIT4_4 in Response to Verticillium Wilt in Gossypium hirsutum

Chromatin remodelers are essential for regulating plant growth, development, and responses to environmental stresses. HIT4 (HEAT-INTOLERANT 4) is a novel stress-induced chromatin remodeling factor that has been less studied in abiotic stress and stress resistance, particularly in cotton. In this study, we conducted a comprehensive analysis of the members of the HIT4 gene family in Gossypium hirsutum using bioinformatics methods, including phylogenetic relationships, gene organization, transcription profiles, phylogenetic connections, selection pressure, and stress response. A total of 18 HIT4 genes were identified in four cotton species, with six HIT4 gene members in upland cotton. Based on the evolutionary relationships shown in the phylogenetic tree, the 18 HIT4 protein sequences were classified into four distinct subgroups. Furthermore, we conducted chromosome mapping to determine the genomic locations of these genes and visually represented the structural characteristics of HIT4 in G. hirsutum. In addition, we predicted the regulatory elements in HIT4 in G. hirsutum and conducted an analysis of repetitive sequences and gene collinearity among HIT4 in four cotton species. Moreover, we calculated the Ka/Ks ratio for homologous genes to assess the selection pressure acting on HIT4. Using RNA-seq, we explored the expression patterns of HIT4 genes in G. hirsutum and Gossypium barbadense. Through weighted gene co-expression network analysis (WGCNA), we found that GHHIT4_4 belonged to the MEblue module, which was mainly enriched in pathways such as DNA replication, phagosome, pentose and glucuronate interconversions, steroid biosynthesis, and starch and sucrose metabolism. This module may regulate the mechanism of upland cotton resistance to Verticillium wilt through DNA replication, phagosome, and various metabolic pathways. In addition, we performed heterologous overexpression of GH_D11G0591 (GHHIT4_4) in tobacco, and the results showed a significant reduction in disease index compared to the wild type, with higher expression levels of disease resistance genes in the transgenic tobacco. After conducting a VIGS (virus-induced gene silencing) experiment in cotton, the results indicated that silencing GHHIT4_4 had a significant impact, the resistance to Verticillium wilt weakened, and the internode length of the plants significantly decreased by 30.7% while the number of true leaves increased by 41.5%. qRT-PCR analysis indicated that GHHIT4_4 mainly enhanced cotton resistance to Verticillium wilt by indirectly regulating the PAL, 4CL, and CHI genes. The subcellular localization results revealed that GHHIT4_4 was predominantly distributed in the mitochondria and nucleus. This study offers preliminary evidence for the involvement of the GHHIT4_4 in cotton resistance to Verticillium wilt and lays the foundation for further research on the disease resistance mechanism of this gene in cotton.

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.

Functional characterization of TBL genes revealed the role of GhTBL7 and GhTBL58 in cotton fiber elongation

TBL (Trichome Birefringence Like) gene family members are involved in trichome initiation and xylan acetylation in several plant species. In our research, we identified 102 TBLs from G. hirsutum. The phylogenetic tree classified TBL genes into five groups. Collinearity analysis of TBL genes indicated 136 paralogous gene pairs in G. hirsutum. Gene duplication indicated that WGD or segmental duplication contributed to the GhTBL gene family expansion. Promoter cis-elements of GhTBLs were related to growth and development, seed-specific regulation, light, and stress responses. GhTBL genes (GhTBL7, GhTBL15, GhTBL21, GhTBL25, GhTBL45, GhTBL54, GhTBL67, GhTBL72, and GhTBL77) exhibited upregulated response under exposure to cold, heat, NaCl, and PEG. GhTBL genes exhibited high expression during fiber development stages. Two GhTBL genes (GhTBL7 and GhTBL58) showed differential expression at 10 DPA fiber, as 10 DPA is a fast fiber elongation stage and fiber elongation is a very important stage of cotton fiber development. Subcellular localization of GhTBL7 and GhTBL58 revealed that these genes reside inside the cell membrane. Promoter GUS activity of GhTBL7 and GhTBL58 exhibited deep staining in roots. To further validate the role of these genes in cotton fiber elongation, we silenced these genes and observed a significant reduction in the fiber length at 10 DPA. In conclusion, the functional study of cell membrane-associated genes (GhTBL7 and GhTBL58) showed deep staining in root tissues and potential function during cotton fiber elongation at 10 DPA fiber.

Elucidation of the GAUT gene family in eight Rosaceae species and function analysis of PbrGAUT22 in pear pollen tube growth

The phylogenetic relationship and evolutionary history of the GAUT gene family were identified in 8 Rosaseae species. PbrGAUT22 was involved in controlling pollen tube growth by regulating the content of pectins. In plants, galacturonosyltransferases (GAUTs) were involved in homogalacturonan biosynthesis and functioned in maintaining pollen tube cell wall integrity. However, the feature and evolutionary history of the GAUT gene family in Rosaceae species and candidates in pear pollen tube growth remain unclear. Here, we identified 190 GAUT genes in 8 Rosaceae species, including Chinese white pear (Pyrus bretschneideri), European pear (Pyrus communis), apple (Malus × domestica), peach (Prunus persica), Japanese apricot (Prunus mume), sweet cherry (Prunus avium), woodland strawberry (Fragaria vesca) and black raspberry (Rubus occidentalis). Members in GAUT gene family were divided into 4 subfamilies according to the phylogenetic and structural analysis. Whole-genome duplication events and dispersed duplicates drove the expansion of the GAUT gene family. Among 23 pollen-expressed PbrGAUT genes in pear, PbrGAUT22 showed increased expression level during 1-6 h post-cultured pollen tubes. PbrGAUT22 was localized to the cytoplasm and plasma membrane. Knockdown of PbrGAUT22 expression in pollen tubes caused the decrease of pectin content and inhibited pear pollen tubes growth. Taken together, we investigated the identification and evolution of the GAUT gene family in Rosaceae species, and found that PbrGAUT22 played an essential role in the synthesis of pectin and the growth of pear pollen tubes.

Genomic identification of cotton SAC genes branded ovule and stress-related key genes in Gossypium hirsutum

SAC genes have been identified to play a variety of biological functions and responses to various stresses. Previously, SAC genes have been recognized in animals and Arabidopsis. For the very first time, we identified 157 SAC genes in eight cotton species including three diploids and five tetraploids with 23 SAC members in G. hirsutum. Evolutionary analysis classified all cotton SAC gene family members into five distinct groups. Cotton SAC genes showed conserved sequence logos and WGD or segmental duplication. Multiple synteny and collinearity analyses revealed gene family expansion and purifying selection pressure during evolution. G. hirsutum SAC genes showed uneven chromosomal distribution, multiple exons/introns, conserved protein motifs, and various growth and stress-related cis-elements. Expression pattern analysis revealed three GhSAC genes (GhSAC3, GhSAC14, and GhSAC20) preferentially expressed in flower, five genes (GhSAC1, GhSAC6, GhSAC9, GhSAC13, and GhSAC18) preferentially expressed in ovule and one gene (GhSAC5) preferentially expressed in fiber. Similarly, abiotic stress treatment verified that GhSAC5 was downregulated under all stresses, GhSAC6 and GhSAC9 were upregulated under NaCl treatment, and GhSAC9 and GhSAC18 were upregulated under PEG and heat treatment respectively. Overall, this study identified key genes related to flower, ovule, and fiber development and important genetic material for breeding cotton under abiotic stress conditions.

Genome-wide analysis elucidates the roles of GhHMA genes in different abiotic stresses and fiber development in upland cotton

The heavy metal-binding domain is involved in heavy metal transporting and plays a significant role in plant detoxification. However, the functions of HMAs are less well known in cotton. In this study, a total of 143 GhHMAs (heavy metal-binding domain) were detected by genome-wide identification in G. hirsutum L. All the GhHMAs were classified into four groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GhHMA genes was highly conserved. 212 paralogous GhHMA gene pairs were identified, and the segmental duplications were the main role to the expansion of GhHMAs. The Ka/Ks values suggested that the GhHMA gene family has undergone purifying selection during the long-term evolutionary process. GhHMA3 and GhHMA75 were located in the plasma membrane, while GhHMA26, GhHMA117 and GhHMA121 were located in the nucleus, respectively. Transcriptomic data and qRT-PCR showed that GhHMA26 exhibited different expression patterns in each tissue and during fiber development or under different abiotic stresses. Overexpressing GhHMA26 significantly promoted the elongation of leaf trichomes and also improved the tolerance to salt stress. Therefore, GhHMA26 may positively regulate fiber elongation and abiotic stress. Yeast two-hybrid assays indicated that GhHMA26 and GhHMA75 participated in multiple biological functions. Our results suggest some genes in the GhHMAs might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GhHMA genes in cotton.

Genome-Wide Identification and Expression Analyses of the Cotton AGO Genes and Their Potential Roles in Fiber Development and Stress Response

Argonaute proteins (AGOs) are indispensable components of RNA silencing. However, systematic characterization of the AGO genes have not been completed in cotton until now. In this study, cotton AGO genes were identified and analyzed with respect to their evolution and expression profile during biotic and abiotic stresses. We identified 14 GaAGO, 14 GrAGO, and 28 GhAGO genes in the genomes of Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum. Cotton AGO proteins were classified into four subgroups. Structural and functional conservation were observed in the same subgroups based on the analysis of the gene structure and conserved domains. Twenty-four duplicated gene pairs were identified in GhAGO genes, and all of them exhibited strong purifying selection during evolution. Moreover, RNA-seq analysis showed that most of the GhAGO genes exhibit high expression levels in the fiber initiation and elongation processes. Furthermore, the expression profiles of GhAGO genes tested by quantitative real-time polymerase chain reaction (qPCR) demonstrated that they were sensitive to Verticillium wilt infection and salt and drought stresses. Overall, our results will pave the way for further functional investigation of the cotton AGO gene family, which may be involved in fiber development and stress response.

Identification and expression analysis of the PtGATL genes under different nitrogen and carbon dioxide treatments in Populus trichocarpa

Pectin is one of the most important components of the plant cell wall. Galacturonosyltransferase-like (GATL) is an important enzyme involved in forming pectin in Arabidopsis thaliana. In this study, 12 PtGATL genes were identified and characterized based on the Populus trichocarpa genome using bioinformatics methods. The results showed that the PtGATLs contained four typical motifs, including DXD, LPPF, GLG, and HXXGXXKPW. According to phylogenetic analysis, PtGATLs were divided into six groups. Chromosome distribution and genome synteny analysis showed that there were 11 segmental-duplicated gene pairs with repeated fragments on chromosomes 2, 5, 7, 8, 10, and 14. Tissue-specific expression profiles indicated that these PtGATLs had different expression patterns. The transcription level of PtGATLs was regulated by different carbon dioxide and nitrogen concentrations. In conclusion, the identification and analysis of PtGATL genes in poplar provide important information on the gene function.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03129-y.

RAD gene family analysis in cotton provides some key genes for flowering and stress tolerance in upland cotton G. hirsutum

Background

RADIALIS (RAD), belongs to the MYB gene family and regulates a variety of functions including floral dorsoventral asymmetry in Antirrhinum majus and development of fruit proteins in Solanum lycopersicum. RAD genes contain an SNF2_N superfamily domain. Here, we comprehensively identified 68 RAD genes from six different species including Arabidopsis and five species of cotton.

Results

Phylogenetic analysis classified RAD genes into five groups. Gene structure, protein motifs and conserved amino acid residues indicated that GhRAD genes were highly conserved during the evolutionary process. Chromosomal location information showed that GhRAD genes were distributed unevenly on different chromosomes. Collinearity and selection pressure analysis indicated RAD gene family expansion in G. hirsutum and G. barbadense with purifying selection pressure. Further, various growth and stress related promotor cis-acting elements were observed. Tissue specific expression level indicated that most GhRAD genes were highly expressed in roots and flowers (GhRAD2, GhRAD3, GhRAD4 and GhRAD11). Next, GhRAD genes were regulated by phytohormonal stresses (JA, BL and IAA). Moreover, Ghi-miRN1496, Ghi-miR1440, Ghi-miR2111b, Ghi-miR2950a, Ghi-miR390a, Ghi-miR390b and Ghi-miR7495 were the miRNAs targeting most of GhRAD genes.

Conclusions

Our study revealed that RAD genes are evolutionary conserved and might be involved in different developmental processes and hormonal stress response. Data presented in our study could be used as the basis for future studies of RAD genes in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08248-z.

Genome-wide characterization and expression analysis of cystathionine β-synthase genes in plant development and abiotic stresses of cotton (Gossypium spp.)

Cystathionine β-synthase (CBS) domains containing proteins (CDCPs) form a large family and play roles in development via regulation of the thioredoxin system as well as abiotic and biotic stress responses of plant. However, the comprehensive study of CBS genes remained elusive in cotton. Here, we identified 237 CBS genes in 11 plant species and the phylogenetic analysis categorized CBS genes into four groups. Whole-genome or segmental with dispersed duplication events contributed to GhCBS gene family expansion. Moreover, orthologous/paralogous genes among three cotton species (G. hirsutum, G. arboreum, and G. raimondii) were detected from the syntenic map among eight plant species. Strong purifying selection for dicotyledonous and monocotyledonous CBS genes, and cis-elements related to plant growth and development, abiotic and hormonal response were observed. Transcriptomic data and qRT-PCR validation of 12 GhCBS genes indicated their critical role in ovule development as most of the genes showed high enrichment. Further, some of GhCBS (GhCBS5, GhCBS16, GhCBS17, GhCBS24, GhCBS25, GhCBS26, and GhCBS52) genes were regulated under various abiotic and hormonal treatments for different time points and involve in ovule and fiber development which provided key genes for future cotton breeding programs. In addition, transgenic tobacco plants overexpressing GhCBS4 transiently exhibited higher water and chlorophyll content indicating improved tolerance toward drought stress. Overall, this study provides the characterization of GhCBS genes for plant growth, abiotic and hormonal stresses, thereby, intimating their significance in cotton molecular breeding for resistant cultivars.

Identification, evolutionary analysis and functional diversification of RAV gene family in cotton (G. hirsutum L.)

Genome wide analysis, expression pattern analysis, and functional characterization of RAV genes highlight their roles in roots, stem development and hormonal response. RAV (Related to ABI3 and VP1) gene family members have been involved in tissues/organs growth and hormone signaling in various plant species. Here, we identified 247 RAVs from 12 different species with 33 RAV genes from G. hirsutum. Phylogenetic analysis classified RAV genes into four distinct groups. Analysis of gene structure showed that most GhRAVs lack introns. Motif distribution pattern and protein sequence logos indicated that GhRAV genes were highly conserved during the process of evolution. Promotor cis-acting elements revealed that promotor regions of GhRAV genes encode numerous elements related to plant growth, abiotic stresses and phytohormones. Chromosomal location information showed uneven distribution of 33 GhRAV genes on different chromosomes. Collinearity analysis identified 628 and 52 orthologous/ paralogous gene pairs in G. hirsutum and G. barbadense, respectively. Ka/Ks values indicated that GhRAV and GbRAV genes underwent strong purifying selection pressure. Selecton model and codon model selection revealed that GhRAV amino acids were under purifying selection and adaptive evolution exists among GhRAV proteins. Three dimensional structure of GhRAVs indicated the presence of numerous alpha helix and beta-barrels. Expression level revealed that some GhRAV genes exhibited high expression in roots (GhRAV3, GhRAV4, GhRAV11, GhRAV18, GhRAV20 and GhRAV30) and stem (GhRAV3 and GhRAV18), indicating their potential role in roots and stem development. GhRAV genes can be regulated by phytohormonal stresses (BL, JA and IAA). Our study provides a reference for future studies related to the functional analysis of GhRAVs in cotton.

Genome-wide analysis of ZAT gene family revealed GhZAT6 regulates salt stress tolerance in G. hirsutum

High salt environments can induce stress in different plants. The genes containing the ZAT domain constitute a family that belongs to a branch of the C₂H₂ family, which plays a vital role in responding to abiotic stresses. In this study, we identified 169 ZAT genes from seven plant species, including 44 ZAT genes from G. hirsutum. Phylogenetic tree analysis divided ZAT genes in six groups with conserved gene structure, protein motifs. Two C₂H₂ domains and an EAR domain and even chromosomal distribution on At and Dt sub-genome chromosomes of G. hirsutum was observed. GhZAT6 was primarily expressed in the root tissue and responded to NaCl and ABA treatments. Subcellular localization found that GhZAT6 was located in the nucleus and demonstrated transactivation activity during a transactivation activity assay. Arabidopsis transgenic lines overexpressing the GhZAT6 gene showed salt tolerance and grew more vigorously than WT on MS medium supplemented with 100 mmol NaCl. Additionally, the silencing of the GhZAT6 gene in cotton plants showed more obvious leaf wilting than the control plants, which were subjected to 400 mmol NaCl treatment. Next, the expressions of GhAPX1, GhFSD1, GhFSD2, and GhSOS3 were significantly lower in the GhZAT6-silenced plants treated with NaCl than the control. Based on these findings, GhZAT6 may be involved in the ABA pathway and mediate salt stress tolerance by regulating ROS-related gene expression.

Identification of BR biosynthesis genes in cotton reveals that GhCPD-3 restores BR biosynthesis and mediates plant growth and development

Brassinosteroid (BR) synthesis genes in different cotton species was comprehensively identified, and the participation of GhCPD-3 in the BR synthesis signaling pathway for regulating plant development was verified. Brassinosteroid is a natural steroidal phytohormone that plays fundamental roles in plant growth and development. In cotton, detailed characterization and functional validation of BR biosynthesis genes remain rare. Here, 16, 8 and 9 BR biosynthesis genes were identified in Gossypium hirsutum, Gossypium raimondii and Gossypium arboreum, respectively, and their phylogenetic relationships, gene structures, conserved motifs of the encoded proteins, chromosomal locations were determined and a synteny analysis was performed. Gossypium hirsutum and Arabidopsis BR biosynthesis genes closely clustered in the phylogenetic tree and fragment duplication was likely the primary cause promoting gene family expansion in G. hirsutum. Gene Ontology (GO) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis showed their relevance as BR biosynthesis genes. GhCPD-3 was highly expressed in roots and stems and the loci of single nucleotide polymorphisms (SNPs) were significantly associated with these traits.Ectopic overexpression of GhCPD-3 in the cpd91 Arabidopsis mutant rescued the mutant phenotype by increasing plant height and leaf size in comparison to those of cpd91 and WT plants. Moreover, overexpressed GhCPD-3 in cpd91 mutants showed greater hypocotyl and root lengths than those of cpd91 and WT plants under light and dark conditions, respectively, indicating that BR actively promotes hypocotyl and root growth. Similar to CPD (CONSTITUTIVE PHOTOMORPHOGENIC DWARF), GhCPD-3 restores BR biosynthesis thereby mediating plant growth and development.

Identification of GhLOG gene family revealed that GhLOG3 is involved in regulating salinity tolerance in cotton (Gossypium hirsutum L.)

Cytokinin (CK) is an important plant hormone that promotes plant cell division and differentiation, and participates in salt response under osmotic stress. LOGs (LONELY GUY) are CK-activating enzymes involved in CK synthesis. The LOG gene family has not been comprehensively characterized in cotton. In this study we identified 151 LOG genes from nine plant species, including 28 LOG genes in Gossypium hirsutum. Phylogenetic analysis divided LOG genes into three groups. Exon/intron structures and protein motifs of GhLOG genes were highly conserved. Synteny analysis revealed that several gene loci were highly conserved between the A and D sub-genomes of G. hirsutum with purifying selection pressure during evolution. Expression profiles showed that most LOG genes were constitutively expressed in eight different tissues. Furthermore, LOG genes can be regulated by abiotic stresses and phytohormone treatments. Moreover, subcellular localization revealed that GhLOG3_At resides inside the cell membrane. Overexpression of GhLOG3 enhanced salt tolerance in Arabidopsis. Virus-induced gene silencing (VIGS) of GhLOG3_At in cotton enhanced sensitivity of plants to salt stress with increased H₂O₂ contents and decreased chlorophyll and proline (PRO) activity. Our results suggested that GhLOG3_At induces salt stress tolerance in cotton, and provides a basis for the use of CK synthesis genes to regulate cotton growth and stress resistance.

Identification and Characterization of the ERF Subfamily B3 Group Revealed GhERF13.12 Improves Salt Tolerance in Upland Cotton

The APETALA2 (AP2)/ethylene response factor plays vital functions in response to environmental stimulus. The ethylene response factor (ERF) subfamily B3 group belongs to the AP2/ERF superfamily and contains a single AP2/ERF domain. Phylogenetic analysis of the ERF subfamily B3 group genes from Arabdiposis thaliana, Gossypium arboreum, Gossypium hirsutum, and Gossypium raimondii made it possible to divide them into three groups and showed that the ERF subfamily B3 group genes are conserved in cotton. Collinearity analysis identified172 orthologous/paralogous gene pairs between G. arboreum and G. hirsutum; 178 between G. hirsutum and G. raimondii; and 1,392 in G. hirsutum. The GhERF subfamily B3 group gene family experienced massive gene family expansion through either segmental or whole genome duplication events, with most genes showing signature compatible with the action of purifying selection during evolution. Most G. hirsutum ERF subfamily B3 group genes are responsive to salt stress. GhERF13.12 transgenic Arabidopsis showed enhanced salt stress tolerance and exhibited regulation of related biochemical parameters and enhanced expression of genes participating in ABA signaling, proline biosynthesis, and ROS scavenging. In addition, the silencing of the GhERF13.12 gene leads to increased sensitivity to salt stress in cotton. These results indicate that the ERF subfamily B3 group had remained conserved during evolution and that GhERF13.12 induces salt stress tolerance in Arabidopsis and cotton.

Genome-Wide Identification and Comparison of Cysteine Proteases in the Pollen Coat and Other Tissues in Maize

Cysteine proteases, belonging to the C1-papain family, play a major role in plant growth and development, senescence, and immunity. There is evidence to suggest that pollen cysteine protease (CP) (ZmCP03) is involved in regulating the anther development and pollen formation in maize. However, there is no report on the genome-wide identification and comparison of CPs in the pollen coat and other tissues in maize. In this study, a total of 38 homologous genes of ZmCP03 in maize were identified. Subsequently, protein motifs, conserved domains, gene structures, and duplication patterns of 39 CPs are analyzed to explore their evolutionary relationship and potential functions. The cis-elements were identified in the upstream sequence of 39 CPs, especially those that are related to regulating growth and development and responding to environmental stresses and hormones. The expression patterns of these genes displayed remarked difference at a tissue or organ level in maize based on the available transcriptome data in the public database. Quantitative reverse transcription PCR (RT-qPCR) analysis showed that ZmCP03 was preferably expressed at a high level in maize pollen. Analyses by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and immunoblot, immunofluorescence and immunogold electron microscopy all validated the cellular localization of ZmCP03 in both the pollen coat and pollen cytoplasm. In addition, 142 CP genes from Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa) and cotton (Gossypium hirsutum), together with 39 maize CPs, were retrieved to analyze their evolution by comparing with orthologous genes. The results suggested that ZmCP03 was relatively conservative and stable during evolution. This study may provide a referential evidence on the function of ZmCP03 in pollen development and germination in maize.