Genome-wide identification of GhAAI genes reveals that GhAAI66 triggers a phase transition to induce early flowering

Transcriptional and translational landscape fine-tune genome annotation and explores translation control in cotton

INTRODUCTION

The unavailability of intergenic region annotation in whole genome sequencing and pan-genomics hinders efforts to enhance crop improvement.

OBJECTIVES

Despite advances in research, the impact of post-transcriptional regulation on fiber development and translatome profiling at different stages of fiber growth in cotton (G. hirsutum) remains unexplored.

METHODS

We utilized a combination of reference-guided de novo transcriptome assembly and ribosome profiling techniques to uncover the hidden mechanisms of translational control in eight distinct tissues of upland cotton.

RESULTS

Our study identified P-site distribution at three-nucleotide periodicity and dominant ribosome footprint at 27 nucleotides. Specifically, we have detected 1,589 small open reading frames (sORFs), including 1,376 upstream ORFs (uORFs) and 213 downstream ORFs (dORFs), as well as 552 long non-coding RNAs (lncRNAs) with potential coding functions, which fine-tune the annotation of the cotton genome. Further, we have identified novel genes and lncRNAs with strong translation efficiency (TE), while sORFs were found to affect mRNA transcription levels during fiber elongation. The reliability of these findings was confirmed by the high consistency in correlation and synergetic fold change between RNA-sequencing (RNA-seq) and Ribosome-sequencing (Ribo-seq) analyses. Additionally, integrated omics analysis of the normal fiber ZM24 and short fiber pag1 cotton mutant revealed several differentially expressed genes (DEGs), and fiber-specific expressed (high/low) genes associated with sORFs (uORFs and dORFs). These findings were further supported by the overexpression and knockdown of GhKCS6, a gene associated with sORFs in cotton, and demonstrated the potential regulation of the mechanism governing fiber elongation on both the transcriptional and post-transcriptional levels.

CONCLUSION

Reference-guided transcriptome assembly and the identification of novel transcripts fine-tune the annotation of the cotton genome and predicted the landscape of fiber development. Our approach provided a high-throughput method, based on multi-omics, for discovering unannotated ORFs, hidden translational control, and complex regulatory mechanisms in crop plants.

GhCOL2 Positively Regulates Flowering by Activating the Transcription of GhHD3A in Upland Cotton (Gossypium hirsutum L.)

Plants have evolved sophisticated signaling networks to adjust flowering time, ensuring successful reproduction. Two crucial flowering regulators, FLOWERING LOCUS T (FT) and CONSTANS (CO), play pivotal roles in regulating flowering across various species. Previous studies have indicated that suppressing Gossypium hirsutum CONSTANS-LIKE 2 (GhCOL2), a homolog of Arabidopsis CO, leads to delayed flowering in cultivated cotton. However, the underlying regulatory mechanisms remain unknown. In this study, a yeast one-hybrid and dual-LUC expression assays were used to elucidate the molecular mechanism through which GhCOL2 regulates the transcription of GhHD3A. RT-qPCR was used to examine the expression of GhCOL2 and GhHD3A. Our findings reveal that GhCOL2 directly binds to CCACA cis-elements and atypical CORE (TGTGTATG) cis-elements in the promoter regions of HEADING DATE 3 A (HD3A), thereby activating GhHD3A transcription. Notably, GhCOL2 and GhHD3A exhibited high expression levels in the adult stage and low levels in the juvenile stage. Interestingly, the expression of GhCOL2 and GhHD3A varied significant between the two cotton varieties (Tx2094 and Maxxa). In summary, our study enhances the understanding of the molecular mechanism by which cotton GhCOL2-GhHD3A regulates flowering at the molecular level. Furthermore, it contributes to a broader comprehension of the GhCOL2-GhHD3A model in G. hirsutum.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

OsAAI1 Increases Rice Yield and Drought Tolerance Dependent on ABA-Mediated Regulatory and ROS Scavenging Pathway

In this study, we investigated the function of OsAAI1 in yield and drought tolerance by constructing overexpression line OE-OsAAI1 and mutant line osaai1. Bioinformatics analysis showed that the AAI gene-OsAAI1- belongs to the HPS_like subfamily of the AAI_LTSS superfamily, and OsAAI1 was localized in the nucleus. The expression of OsAAI1 was significantly induced by ABA and drought stress. OsAAI1 overexpression (OE19) significantly increased, and gene mutant (osaai1-1) repressed plant height, primary root length, lateral root number, grain size and yield in rice. Moreover, physiological and biochemical analyses showed that osaai1 was sensitive to drought stress, while OE19 enhanced the drought tolerance in rice. DAB and NBT staining revealed that under drought treatment, osaai1 accumulated a large amount of ROS compared with the wild type, while OE19 accumulated the least, and CAT, APX, GPX, GR activities were higher in OE19 and lower in osaai1, suggesting that OE19 improves rice tolerance to drought stress by enhancing ROS scavenging ability. OE19 also induce the expression of ABA-mediated regulatory pathway genes and enhance accumulation of ABA content in rice seedling. Predictably, OE19 displayed enhanced sensitivity to ABA, and ROS accumulation was significantly higher than in wild type and osaai1 under 3 µM ABA treatment. Thus, these results suggest that OsAAI1 is a positive regulator of rice yield and drought tolerance dependent on the ABA-mediated regulatory and ROS scavenging pathway.

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.

Genome-wide identification of the geranylgeranyl pyrophosphate synthase (GGPS) gene family involved in chlorophyll synthesis in cotton

BACKGROUND

Geranylgeranyl pyrophosphate synthase (GGPS) is a structural enzyme of the terpene biosynthesis pathway that is involved in regulating plant photosynthesis, growth and development, but this gene family has not been systematically studied in cotton.

RESULTS

In the current research, genome-wide identification was performed, and a total of 75 GGPS family members were found in four cotton species, Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii. The GGPS genes were divided into three subgroups by evolutionary analysis. Subcellular localization prediction showed that they were mainly located in chloroplasts and plastids. The closely related GGPS contains a similar gene structure and conserved motif, but some genes are quite different, resulting in functional differentiation. Chromosome location analysis, collinearity and selection pressure analysis showed that many fragment duplication events occurred in GGPS genes. Three-dimensional structure analysis and conservative sequence analysis showed that the members of the GGPS family contained a large number of α-helices and random crimps, and all contained two aspartic acid-rich domains, DDxxxxD and DDxxD (x is an arbitrary amino acid), suggesting its key role in function. Cis-regulatory element analysis showed that cotton GGPS may be involved in light response, abiotic stress and other processes. A GGPS gene was silenced successfully by virus-induced gene silencing (VIGS), and it was found that the chlorophyll content in cotton leaves decreased significantly, suggesting that the gene plays an important role in plant photosynthesis.

CONCLUSIONS

In total, 75 genes were identified in four Gossypium species by a series of bioinformatics analysis. Gene silencing from GGPS members of G. hirsutum revealed that GGPS plays an important regulatory role in photosynthesis. This study provides a theoretical basis for the biological function of GGPS in cotton growth and development.

Genome-wide identification of NHX (Na+/H+ antiporter) gene family in Cucurbita L. and functional analysis of CmoNHX1 under salt stress

Soil salinization, which is the accumulation of salt in soil, can have a negative impact on crop growth and development by creating an osmotic stress that can reduce water uptake and cause ion toxicity. The NHX gene family plays an important role in plant response to salt stress by encoding for Na⁺/H⁺ antiporters that help regulate the transport of sodium ions across cellular membranes. In this study, we identified 26 NHX genes in three cultivars of Cucurbita L., including 9 Cucurbita moschata NHXs (CmoNHX1-CmoNHX9), 9 Cucurbita maxima NHXs (CmaNHX1-CmaNHX9) and 8 Cucurbita pepo NHXs (CpNHX1-CpNHX8). The evolutionary tree splits the 21 NHX genes into three subfamilies: the endosome (Endo) subfamily, the plasma membrane (PM) subfamily, and the vacuole (Vac) subfamily. All the NHX genes were irregularly distributed throughout the 21 chromosomes. 26 NHXs were examined for conserved motifs and intron-exon organization. These findings suggested that the genes in the same subfamily may have similar functions while genes in other subfamilies may have functional diversity. The circular phylogenetic tree and collinearity analysis of multi-species revealed that Cucurbita L. had a substantially greater homology relationship than Populus trichocarpa and Arabidopsis thaliana in terms of NHX gene homology. We initially examined the cis-acting elements of the 26 NHXs in order to investigate how they responded to salt stress. We discovered that the CmoNHX1, CmaNHX1, CpNHX1, CmoNHX5, CmaNHX5, and CpNHX5 all had numerous ABRE and G-box cis-acting elements that were important to salt stress. Previous transcriptome data showed that in the mesophyll and veins of leaves, many CmoNHXs and CmaNHXs, such as CmoNHX1, responded significantly to salt stress. In addition, we heterologously expressed in A. thaliana plants in order to further confirm the response of CmoNHX1 to salt stress. The findings demonstrated that during salt stress, A. thaliana that had CmoNHX1 heterologously expression was found to have decreased salt tolerance. This study offers important details that will aid in further elucidating the molecular mechanism of NHX under salt stress.

Characterization and expression analysis of bHLH transcription factors reveal their putative regulatory effects on nectar spur development in Aquilegia species

Nectar spur is a hollow extension of certain flower parts and shows strikingly diverse size and shape in Aquilegia. Nectar spur development is involved in cell division and expansion processes. The basic helix-loop-helix (bHLH) transcription factors (TFs) control a diversity of organ morphogenesis, including cell division and cell expansion processes. However, the role of bHLH genes in nectar spur development in Aquilegia is mainly unknown. We conducted a genome-wide identification of the bHLH gene family in Aquilegia to determine structural characteristics and phylogenetic relationships, and to analyze expression profiles of these genes during the development of nectar spur in spurless and spurred species. A total of 120 AqbHLH genes were identified from the Aquilegia coerulea genome. The phylogenetic tree showed that AqbHLH proteins were divided into 15 subfamilies, among which S7 and S8 subfamilies occurred marked expansion. The AqbHLH genes in the same clade had similar motif composition and gene structure characteristics. Conserved residue analysis indicated nineteen residues with conservation of more than 50% were found in the four conserved regions. In the upstream sequence of AqbHLH genes, the light-responsive element was the most abundant cis-acting element. Eighteen AqbHLH genes showed syntenic relationships, and eight genes from four syntenic pairs underwent tandem duplications. According to the expression profiling analysis by public RNA-Seq data and qRT-PCR results, five AqbHLH genes, including AqbHLH027, AqbHLH046, AqbHLH082, AqbHLH083 and AqbHLH092, were differentially expressed between different tissues in A. coerulea at early developmental stages, as well as between spurless and spurred Aquilegia species. Of them, AqbHLH046 was not only highly expressed in spur compared with blade, but also showed higher expression levels in spurred species than spurless specie, suggesting it plays an essential role in the development of spur by regulating cell division. This study lays a foundation to investigate the function of AqbHLH genes family in nectar spur development, and has potential implications for speciation and genetic breeding in the genus Aquilegia.

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.

Recent advances and future perspectives in early-maturing cotton research

Cotton's fundamental requirements for long periods of growth and specific seasonal temperatures limit the global arable areas that can be utilized to cultivate cotton. This constraint can be alleviated by breeding for early-maturing varieties. By delaying the sowing dates without impacting the boll-opening time, early-maturing varieties not only mitigate the yield losses brought on by unfavorable weathers in early spring and late autumn but also help reducing the competition between cotton and other crops for arable land, thereby optimizing the cropping system. This review presents studies and breeding efforts for early-maturing cotton, which efficiently pyramid early maturity, high-quality, multiresistance traits, and suitable plant architecture by leveraging pleiotropic genes. Attempts are also made to summarize our current understanding of the molecular mechanisms underlying early maturation, which involves many pathways such as epigenetic, circadian clock, and hormone signaling pathways. Moreover, new avenues and effective measures are proposed for fine-scale breeding of early-maturing crops to ensure the healthy development of the agricultural industry.

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.

Updated role of ABA in seed maturation, dormancy, and germination

•

Functional ABA biosynthesis genes show specific roles for ABA accumulation at different stages of seed development and seedling establishment.

•

De novo ABA biosynthesis during embryogenesis is required for late seed development, maturation, and induction of primary dormancy.

•

ABA plays multiple roles with the key LAFL hub to regulate various downstream signaling genes in seed and seedling development.

•

Key ABA signaling genes ABI3, ABI4, and ABI5 play important multiple functions with various cofactors during seed development such as de-greening, desiccation tolerance, maturation, dormancy, and seed vigor.

•

The crosstalk between ABA and other phytohormones are complicated and important for seed development and seedling establishment.

Background

Seed is vital for plant survival and dispersion, however, its development and germination are influenced by various internal and external factors. Abscisic acid (ABA) is one of the most important phytohormones that influence seed development and germination. Until now, impressive progresses in ABA metabolism and signaling pathways during seed development and germination have been achieved. At the molecular level, ABA biosynthesis, degradation, and signaling genes were identified to play important roles in seed development and germination. Additionally, the crosstalk between ABA and other hormones such as gibberellins (GA), ethylene (ET), Brassinolide (BR), and auxin also play critical roles. Although these studies explored some actions and mechanisms by which ABA-related factors regulate seed morphogenesis, dormancy, and germination, the complete network of ABA in seed traits is still unclear.

Aim of review

Presently, seed faces challenges in survival and viability. Due to the vital positive roles in dormancy induction and maintenance, as well as a vibrant negative role in the seed germination of ABA, there is a need to understand the mechanisms of various ABA regulators that are involved in seed dormancy and germination with the updated knowledge and draw a better network for the underlying mechanisms of the ABA, which would advance the understanding and artificial modification of the seed vigor and longevity regulation.

Key scientific concept of review

Here, we review functions and mechanisms of ABA in different seed development stages and seed germination, discuss the current progresses especially on the crosstalk between ABA and other hormones and signaling molecules, address novel points and key challenges (e.g., exploring more regulators, more cofactors involved in the crosstalk between ABA and other phytohormones, and visualization of active ABA in the plant), and outline future perspectives for ABA regulating seed associated traits.

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.

The emerging role of jasmonate in the control of flowering time

Plants dynamically synchronize their flowering time with changes in the internal and external environments through a variety of signaling pathways to maximize fitness. In the last two decades, the major pathways associated with flowering, including the photoperiod, vernalization, age, autonomous, gibberellin, and ambient temperature pathways, have been extensively analyzed. In recent years, an increasing number of signals, such as sugar, thermosensory, stress, and certain hormones, have been shown to be involved in fine-tuning flowering time. Among these signals, the jasmonate signaling pathway has a function in the determination of flowering time that has not been systematically summarized. In this review, we present an overview of current knowledge of jasmonate control of flowering and discuss jasmonate crosstalk with other signals (such as gibberellin, defense, and touch) during floral transition.

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.

Genome-wide Identification and Evolution of the PP2C Gene Family in Eight Rosaceae Species and Expression Analysis Under Stress in Pyrus bretschneideri

Type 2C protein phosphatase (PP2C) plays an essential role in abscisic acid (ABA) signaling transduction processes. In the current study, we identify 719 putative PP2C genes in eight Rosaceae species, including 118 in Chinese white pear, 110 in European pear, 73 in Japanese apricot, 128 in apple, 74 in peach, 65 in strawberry, 78 in sweet cherry, and 73 in black raspberry. Further, the phylogenetic analysis categorized PbrPP2C genes of Chinese white pear into twelve subgroups based on the phylogenic analysis. We observed that whole-genome duplication (WGD) and dispersed gene duplication (DSD) have expanded the Rosaceae PP2C family despite simultaneous purifying selection. Expression analysis finds that PbrPP2C genes have organ-specific functions. QRT-PCR validation of nine PbrPP2C genes of subgroup A indicates a role in ABA-mediated response to abiotic stress. Finally, we find that five PbrPP2C genes of subgroup A function in the nucleus. In summary, our research suggests that the PP2C family functions to modulate ABA signals and responds to abiotic stress.

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 and Analysis of GhEXO Gene Family Indicated That GhEXO7_At Promotes Plant Growth and Development Through Brassinosteroid Signaling in Cotton (Gossypium hirsutum L.)

Brassinosteroids (BRs), an efficient plant endogenous hormone, significantly promotes plant nutrient growth adapting to biological and abiotic adversities. BRs mainly promote plant cell elongation by regulating gene expression patterns. EXORDIUM (EXO) genes have been characterized as the indicators of BR response genes. Cotton, an ancient crop, is of great economic value and its fibers can be made into all kinds of fabrics. However, EXO gene family genes have not been full identified in cotton. 175 EXO genes were identified in nine plant species, of which 39 GhEXO genes in Gossypium hirsutum in our study. A phylogenetic analysis grouped all of the proteins encoded by the EXO genes into five major clades. Sequence identification of conserved amino acid residues among monocotyledonous and dicotyledonous species showed a high level of conservation across the N and C terminal regions. Only 25% the GhEXO genes contain introns besides conserved gene structure and protein motifs distribution. The 39 GhEXO genes were unevenly distributed on the 18 At and Dt sub-genome chromosomes. Most of the GhEXO genes were derived from gene duplication events, while only three genes showed evidence of tandem duplication. Homologous locus relationships showed that 15 GhEXO genes are located on collinear blocks and that all orthologous/paralogous gene pairs had Ka > Ks values, indicating purifying selection pressure. The GhEXO genes showed ubiquitous expression in all eight tested cotton tissues and following exposure to three phytohormones, IAA, GA, and BL. Furthermore, GhEXO7_At was mainly expressed in response to BL treatment, and was predominantly expressed in the fibers. GhEXO7_At was found to be a plasma membrane protein, and its ectopic expression in Arabidopsis mediated BR-regulated plant growth and development with altered expression of DWF4, CPD, KCS1, and EXP5. Additionally, the functions of GhEXO7_At were confirmed by virus-induced gene silencing (VIGS) in cotton. This study will provide important genetic resources for future cotton breeding programs.

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 Evolutionary Analysis of Gossypium Tubby-Like Protein (TLP) Gene Family and Expression Analyses During Salt and Drought Stress

Tubby-like proteins (TLPs) possess a highly conserved closed β barrel tubby domain at C-terminal and N-terminal F-box. The role of TLP gene family members has been widely discussed in numerous organisms; however, the detailed genome-wide study of this gene family in Gossypium species has not been reported till date. Here, we systematically identified 105 TLP gene family members in cotton (Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense) genomes and classified them into eight phylogenetic groups. Cotton TLP12 gene family members clustered into two groups, 4 and 8. They experienced higher evolutionary pressure in comparison to others, indicating the faster evolution in both diploid as well as in tetraploid cotton. Cotton TLP gene family members expanded mainly due to segmental duplication, while only one pair of tandem duplication was found in cotton TLPs paralogous gene pairs. Subsequent qRT-PCR validation of seven putative key candidate genes of GhTLPs indicated that GhTLP11A and GhTLP12A.1 genes were highly sensitive to salt and drought stress. The co-expression network, pathways, and cis-regulatory elements of GhTLP11A and GhTLP12A.1 genes confirmed their functional importance in salt and drought stress responses. This study proposes the significance of GhTLP11A and GhTLP12A.1 genes in exerting control over salt and drought stress responses in G. hirsutum and also provides a reference for future research, elaborating the biological roles of G. hirsutum TLPs in both stress responses.

Genome-wide identification of the BASS gene family in four Gossypium species and functional characterization of GhBASSs against salt stress

Bile acid sodium symporter (BASS) family proteins encode a class of sodium/solute symporters. Even though the sodium transporting property of BASSs in mammals was well studied, their sodium transportability and functional roles in plant salt tolerance remained largely unknown. Here, BASS family members from 4 cotton species, as well as 30 other species were identified. Then, they were designated as members of BASS1 to BASS5 subfamilies according to their sequence similarity and phylogenetic relationships. There were 8, 11, 16 and 18 putative BASS genes in four cotton species. While whole-genome duplications (WGD) and segmental duplications rendered the expansion of the BASS gene family in cotton, BASS gene losses occurred in the tetraploid cotton during the evolution from diploids to allotetraploids. Concerning functional characterizations, the transcript profiling of GhBASSs revealed that they not only preferred tissue-specific expression but also were differently induced by various stressors and phytohormones. Gene silencing and overexpression experiments showed that GhBASS1 and GhBASS3 positively regulated, whereas GhBASS2, GhBASS4 and GhBASS5 negatively regulated plant salt tolerance. Taken together, BASS family genes have evolved before the divergence from the common ancestor of prokaryotes and eukaryotes, and GhBASSs are plastidial sodium-dependent metabolite co-transporters that can influence plant salt tolerance.

High-resolution temporal dynamic transcriptome landscape reveals a GhCAL-mediated flowering regulatory pathway in cotton (Gossypium hirsutum L.)

The transition from vegetative to reproductive growth is very important for early maturity in cotton. However, the genetic control of this highly dynamic and complex developmental process remains unclear. A high-resolution tissue- and stage-specific transcriptome profile was generated from six developmental stages using 72 samples of two early-maturing and two late-maturing cotton varieties. The results of histological analysis of paraffin sections showed that flower bud differentiation occurred at the third true leaf stage (3TLS) in early-maturing varieties, but at the fifth true leaf stage (5TLS) in late-maturing varieties. Using pairwise comparison and weighted gene co-expression network analysis, 5312 differentially expressed genes were obtained, which were divided into 10 gene co-expression modules. In the MElightcyan module, 46 candidate genes regulating cotton flower bud differentiation were identified and expressed at the flower bud differentiation stage. A novel key regulatory gene related to flower bud differentiation, GhCAL, was identified in the MElightcyan module. Anti-GhCAL transgenic cotton plants exhibited late flower bud differentiation and flowering time. GhCAL formed heterodimers with GhAP1-A04/GhAGL6-D09 and regulated the expression of GhAP1-A04 and GhAGL6-D09. GhAP1-A04- and GhAGL6-D09-silenced plants also showed significant late flowering. Finally, we propose a new flowering regulatory pathway mediated by GhCAL. This study elucidated the molecular mechanism of cotton flowering regulation and provides good genetic resources for cotton early-maturing breeding.

Genome-Wide Identification and Expression Analysis of the NHX (Sodium/Hydrogen Antiporter) Gene Family in Cotton

The sodium/hydrogen antiporter (NHX) gene family with the Na⁺/H⁺ exchange protein domain is a transporter of sodium and hydrogen ions and plays an important role in the response of plants to salt stress. Studying the response of cotton to salt stress through comprehensive identification and analysis of NHX genes in several species and their roles in salt tolerance mechanisms is of great significance. In this study, 23, 24, 12, and 12 NHX genes were identified from Gossypium hirsutum (Gh), G. barbadense, G. arboreum and G. raimondii, respectively. Phylogenetic analysis showed that these genes were mainly divided into three clades with significant subcellular localization, namely, endosome (Endo-class), plasma membrane (PM-class) and vacuole (Vac-class). By analyzing the structure of NHX genes and proteins, each branch of the NHX gene family was found to be structurally conserved, and collinearity analysis showed that NHX genes were mainly expressed through whole genome and segmental duplication. The non-synonymous (Ka)/synonymous (Ks) values showed that the NHX gene family experienced strong purifying selection during long-term evolution. Cis-acting element analysis showed that the NHX gene family may be related to the regulation of abscisic acid (ABA) and methyl jasmonate (MeJA) hormones. Additionally, transcriptomic data analysis and qRT-PCR showed that GhNHXs exhibited different expression patterns in each tissue and under different salinities. These results provide an important reference for us to further understand and analyze the molecular regulation mechanism of cotton NHX genes.

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

Background

GGPP (geranylgeranyl diphosphate) is produced in the isoprenoid pathway and mediates the function of various plant metabolites, which is synthesized by GGPPS (GGPP synthases) in plants. GGPPS characterization has not been performed in any plant species except Arabidopsis thaliana. Here, we performed a complete computational and bioinformatics analysis of GGPPS and detected their transcription expression pattern in Gossypium hirsutum for the first time so that to explore their evolutionary relationship and potential functions. Finally, we unravelled evolutionary relationship, conserved sequence logos, gene duplication and potential involvement in plant development and abiotic stresses tolerance of GGPPS genes in G. hirsutum and other plant species.

Results

A total of 159 GGPPS genes from 18 plant species were identified and evolutionary analysis divided these GGPPS genes into five groups to indicate their divergence from a common ancestor. Further, GGPPS family genes were conserved during evolution and underwent segmental duplication. The identified 25 GhGGPPS genes showed diverse expression pattern particularly in ovule and fiber development indicating their vital and divers roles in the fiber development. Additionally, GhGGPPS genes exhibited wide range of responses when subjected to abiotic (heat, cold, NaCl and PEG) stresses and hormonal (BL, GA, IAA, SA and MeJA) treatments, indicating their potential roles in various biotic and abiotic stresses tolerance.

Conclusions

The GGPPS genes are evolutionary conserved and might be involve in different developmental stages and stress response. Some potential key genes (e.g. GhGGPP4, GhGGPP9, and GhGGPP15) were suggested for further study and provided valuable source for cotton breeding to improve fiber quality and resistant to various stresses.

Gossypium Genomics: Trends, Scope, and Utilization for Cotton Improvement

Cotton (Gossypium spp.) is the most important natural fiber crop worldwide. The diversity of Gossypium species also provides an ideal model for investigating evolution and domestication of polyploids. However, the huge and complex cotton genome hinders genomic research. Technical advances in high-throughput sequencing and bioinformatics analysis have now largely overcome these obstacles, bringing about a new era of cotton genomics. Here, we review recent progress in Gossypium genomics based on whole genome sequencing, resequencing, and comparative genomics, which have provided insights about the genomic basis of fiber biogenesis and the landscape of cotton functional genomics. We address current challenges and present multidisciplinary genomics-enabled breeding strategies covering the breadth of high fiber yield, quality, and environmental resilience for future cotton breeding programs.