The Chalcone Isomerase Family in Cotton: Whole-Genome Bioinformatic and Expression Analyses of the Gossypium barbadense L. Response to Fusarium Wilt Infection

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.

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

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

Cloning, Identification, and Functional Analysis of the Chalcone Isomerase Gene from Astragalus sinicus

Astragalus sinicus is an important winter-growing cover crop. It is widely utilized, not only as a cover crop for its benefits in fertilizing the soil but also as a landscape ground cover plant. Anthocyanins are involved in the pigmentation of plants in leaves and flowers, which is a crucial characteristic trait for A. sinicus. The formation of anthocyanins depends significantly on the enzyme chalcone isomerase (CHI). However, research on the CHI gene of A. sinicus remains unexplored. The rapid amplification of cDNA ends (RACE) approach was used in this research to clone the CHI sequence from A. sinicus (AsiCHI). The expression profiles of the AsiCHI gene in multiple tissues of A. sinicus were subsequently examined by qRT-PCR (Quantitative Real-Time PCR). Furthermore, the function of the AsiCHI was identified by the performance of ectopic expression in Arabidopsis (Arabidopsis thaliana). The outcomes revealed that the full-length cDNA of the AsiCHI gene (GeneBank: OQ870547) measured 972 bp in length and included an open reading frame of 660 bp. The encoded protein contains 219 amino acids with a molecular weight of 24.14 kDa and a theoretical isoelectric point of 5.11. In addition, the remarkable similarity between the AsiCHI protein and the CHI proteins of other Astragalus species was demonstrated by the sequence alignment and phylogenetic analysis. Moreover, the highest expression level of AsiCHI was observed in leaves and showed a positive correlation with anthocyanin content. The functional analysis further revealed that the overexpression of AsiCHI enhanced the anthocyanidin accumulation in the transgenic lines. This study provided a better understanding of AsiCHI and elucidated its role in anthocyanin production.

Integrated transcriptome and microRNA sequencing analyses reveal gene responses in poplar leaves infected by the novel pathogen bean common mosaic virus (BCMV)

Recently, a novel poplar mosaic disease caused by bean common mosaic virus (BCMV) was investigated in Populus alba var. pyramidalis in China. Symptom characteristics, physiological performance of the host, histopathology, genome sequences and vectors, and gene regulation at the transcriptional and posttranscriptional levels were analyzed and RT-qPCR (quantitative reverse transcription PCR) validation of expression was performed in our experiments. In this work, the mechanisms by which the BCMV pathogen impacts physiological performance and the molecular mechanisms of the poplar response to viral infection were reported. The results showed that BCMV infection decreased the chlorophyll content, inhibited the net photosynthesis rate (Pn) and stomatal conductance (Gs), and significantly changed chlorophyll fluorescence parameters in diseased leaves. Transcriptome analysis revealed that the expression of the majority of DEGs (differentially expressed genes) involved in the flavonoid biosynthesis pathway was promoted, but the expression of all or almost all DEGs associated with photosynthesis-antenna proteins and the photosynthesis pathway was inhibited in poplar leaves, suggesting that BCMV infection increased the accumulation of flavonoids but decreased photosynthesis in hosts. Gene set enrichment analysis (GSEA) illustrated that viral infection promoted the expression of genes involved in the defense response or plant-pathogen interaction. MicroRNA-seq analysis illustrated that 10 miRNA families were upregulated while 6 families were downregulated in diseased poplar leaves; moreover, miR156, the largest family with the most miRNA members and target genes, was only differentially upregulated in long-period disease (LD) poplar leaves. Integrated transcriptome and miRNA-seq analyses revealed 29 and 145 candidate miRNA-target gene pairs; however, only 17 and 76 pairs, accounting for 2.2% and 3.2% of all DEGs, were authentically negatively regulated in short-period disease (SD) and LD leaves, respectively. Interestingly, 4 miR156/SPL (squamosa promoter-binding-like protein) miRNA-target gene pairs were identified in LD leaves: the miR156 molecules were upregulated, but SPL genes were downregulated. In conclusion, BCMV infection significantly changed transcriptional and posttranscriptional gene expression in poplar leaves, inhibited photosynthesis, increased the accumulation of flavonoids, induced systematic mosaic symptoms, and decreased physiological performance in diseased poplar leaves. This study elucidated the fine-tuned regulation of poplar gene expression by BCMV; moreover, the results also suggested that miR156/SPL modules played important roles in the virus response and development of viral systematic symptoms in plant virus disease.

Integrative analysis of transcriptome and metabolome reveals the effect of DNA methylation of chalcone isomerase gene in promoter region on Lithocarpus polystachyus Rehd flavonoids

Metabolite biosynthesis is regulated by gene expression, which is altered by DNA methylation in the promoter region. Chalcone isomerase (CHI) gene encodes a key enzyme in the Lithocarpus polystachyus Rehd flavonoid pathway, and the expression of L. polystachyus CHI (LpCHI) is closely related to the synthesis of flavonoid metabolites. In this study, we analyzed the DNA methylation site of the LpCHI promoter and its effect on gene expression and metabolite accumulation. The proportions of three types of LpCHI promoter DNA methylation are 7.5%, 68.75%, 18.75%, determined by bisulfite sequencing. Transcriptome sequencing shows that LpCHI is strongly up-regulated in LpCHI promoter methylation Type A but down-regulated in LpCHI promoter methylation Type B and Type C. The expression of LpCHI shows no significant difference between Type B and Type C. Moreover, nine kinds of differentially expressed transcription factors (DETFs) bind to seven CpG-sites of the LpCHI promoter region to regulate LpCHI expression. The results of metabolomics show that differentially accumulated flavonoids are higher in LpCHI promoter methylation Type A than in LpCHI promoter methylation Type B and Type C. Additionally, a positive correlation was found between the LpCHI expression and flavonoids accumulation. These results show that the effect of CpG site-specificity on gene transcription is great than that of overall promoter DNA methylation on gene transcription. The mechanisms of flavonoid genes regulating metabolite accumulation are further revealed.

Genome-Wide Classification and Evolutionary Analysis Reveal Diverged Patterns of Chalcone Isomerase in Plants

Flavonoids as a class of important secondary metabolites are widely present in land plants, and chalcone isomerase (CHI) is the key rate-limiting enzyme that participates in catalyzing the stereospecific isomerization of chalcones to yield their corresponding flavanones. However, the phylogenetic dynamics and functional divergence of CHI family genes during the evolutionary path of green plants remains poorly understood. Here, a total of 122 CHI genes were identified by performing a genome-wide survey of 15 representative green plants from the most ancestral basal plant chlorophyte algae to higher angiosperm plants. Phylogenetic, orthologous groups (OG) classification, and genome structure analysis showed that the CHI family genes have evolved into four distinct types (types I–IV) containing eight OGs after gene duplication, and further studies indicated type III CHIs consist of three subfamilies (FAP1, FAP2, and FAP3). The phylogeny showed FAP3 CHIs as an ancestral out-group positioned on the outer layers of the main branch, followed by type IV CHIs, which are placed in an evolutionary intermediate between FAP3 CHIs and bona fide CHIs (including type I and type II). The results imply a potential intrinsic evolutionary connection between CHIs existing in the green plants. The amino acid substitutions occurring in several residues have potentially affected the functional divergence between CHI proteins. This is supported by the analysis of transcriptional divergence and cis-acting element analysis. Evolutionary dynamics analyses revealed that the differences in the total number of CHI family genes in each plant are primarily attributed to the lineage-specific expansion by natural selective forces. The current studies provide a deeper understanding of the phylogenetic relationships and functional diversification of CHI family genes in green plants, which will guide further investigation on molecular characteristics and biological functions of CHIs.

Genome-Wide Identification of the SAMS Gene Family in Upland Cotton (Gossypium hirsutum L.) and Expression Analysis in Drought Stress Treatments

Cotton is an important commercial crop whose growth and yield are severely affected by drought. S-adenosylmethionine (SAM) is widely involved in the plant stress response and growth regulation; however, the role of the S-adenosylmethionine synthase (SAMS) gene family in this process is poorly understood. Here, we systematically analyzed the expression of SAMS genes in Upland Cotton (Gossypium hirsutum L.). A total of 16 SAMS genes were identified, each with a similar predicted structure. A large number of cis-acting elements involved in the response to abiotic stress were predicted based on promoter analysis, indicating a likely important role in abiotic stress responses. The results of qRT-PCR validation showed that GhSAMS genes had different expression patterns after drought stress and in response to drought stress. Analysis of a selected subset of GhSAMS genes showed increased expression in cultivar Xinluzhong 39 (drought resistant) when compared to cultivar Xinluzao 26 (drought-sensitive) upland cotton. This study provides important relevant information for further study of SAMS genes in drought resistance research of upland cotton, which is helpful for drought-resistance improvement of upland cotton.

Genome-Wide Identification and Expression Profiles of 13 Key Structural Gene Families Involved in the Biosynthesis of Rice Flavonoid Scaffolds

Flavonoids are a class of key polyphenolic secondary metabolites with broad functions in plants, including stress defense, growth, development and reproduction. Oryza sativa L. (rice) is a well-known model plant for monocots, with a wide range of flavonoids, but the key flavonoid biosynthesis-related genes and their molecular features in rice have not been comprehensively and systematically characterized. Here, we identified 85 key structural gene candidates associated with flavonoid biosynthesis in the rice genome. They belong to 13 families potentially encoding chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonol synthase (FLS), leucoanthocyanidin dioxygenase (LDOX), anthocyanidin synthase (ANS), flavone synthase II (FNSII), flavanone 2-hydroxylase (F2H), flavonoid 3′-hydroxylase (F3′H), flavonoid 3′,5′-hydroxylase (F3′5′H), dihydroflavonol 4-reductase (DFR), anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR). Through structural features, motif analyses and phylogenetic relationships, these gene families were further grouped into five distinct lineages and were examined for conservation and divergence. Subsequently, 22 duplication events were identified out of a total of 85 genes, among which seven pairs were derived from segmental duplication events and 15 pairs were from tandem duplications, demonstrating that segmental and tandem duplication events play important roles in the expansion of key flavonoid biosynthesis-related genes in rice. Furthermore, these 85 genes showed spatial and temporal regulation in a tissue-specific manner and differentially responded to abiotic stress (including six hormones and cold and salt treatments). RNA-Seq, microarray analysis and qRT-PCR indicated that these genes might be involved in abiotic stress response, plant growth and development. Our results provide a valuable basis for further functional analysis of the genes involved in the flavonoid biosynthesis pathway in rice.

Genome-wide identification and expression analysis of the 2OG-Fe(II) oxygenase gene family in upland cotton (Gossypium hirsutum L.)

The 2OG-Fe(II) oxygenase (RF) family of enzyme proteins can affect bulliform cells and cause leaf curling. However, there are few studies related to this family in cotton, and there has been no systematic analysis of RF genes. Here, we determined 25 RF genes in the complete genome sequence of upland cotton (Gossypium hirsutum L.) and 11 RF genes in the complete genome sequence of Arabidopsis thaliana. Cotton RF proteins can be divided into three categories. Whole genome/fragment and scattered replication events played an important role in the expansion of the RF gene family. qRT-PCR analysis results showed that RF genes respond to drought stress Pairwise comparison results showed that the expression of RF genes in Shi yuan 321 was higher than that in Kui 85-174. Overall, genome-wide identification approach was used to further analyze the related functions of the RF gene family, which may include the response to drought stress, in cotton.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01065-4.

Molecular characterization, expression and interaction of MAPK, MAPKK and MAPKKK genes in upland cotton

Mitogen-activated protein kinase (MAPK) signaling cascades, consisting of three types of sequentially phosphorylated kinases (MAPKKK, MAPKK, and MAPK), play vital roles in various processes including plant development and stress response. In this study, 52 GhMAPKs, 23 GhMAPKKs, and 166 GhMAPKKKs were identified in upland cotton. Chromosomal locations, gene duplication and structure, motifs, cis-regulatory elements, and protein subcellular localization were further analyzed. With the identified MAPK cascade genes in G. arboretum and G. raimondii, a syntenic diagram of three cotton species was constructed. The interactions of seven GhMAPK cascade genes were investigated. Two complete signaling modules were defined: The GhMEKK24/GhMEKK31-GhMAPKK9-GhMAPK10 and GhMEKK3/GhMEKK24/GhMEKK31-GhMAPKK16-GhMAPK10/GhMAPK11 cascades. Moreover, interaction networks and the interaction pairs were combined with their expression patterns and demonstrated that the network mediated by the MAPK signaling cascade participates in abiotic stress signaling. Our research provides a foundation for studying the molecular mechanism of the MAPK signaling pathway under abiotic stress.

Adaptive Evolution of Chalcone Isomerase Superfamily in Fagaceae

Chalcone Isomerase (CHI) catalyzes the biosynthesis of flavonoids and secondary metabolism in plants. Currently, there is no systematic analysis of CHIs gene family in Fagaceae which is available. In this study, twenty-two CHI proteins were identified in five species of the Fagaceae family. The CHI superfamily in Fagaceae can be classified into three subfamilies and five groups using phylogenetic analysis, analysis of physicochemical properties, and structural prediction. Results indicated that serine (Ser) and isoleucine (Ile) residues determine the substrate preferred by active Type I Fagaceae CHI, and the chalcone isomerase-like (CHIL) of Fagaceae had active site residues. Adaptive analysis of CHIs showed that CHIs are subject to selection pressure. The active CHI gene of Fagaceae was located in the cytoplasm, and it had the typical gene structure of CHI and contains four exons. All the twenty-two identified CHIs had the conserved domain motif 3, and the different groups had their own structural characteristics. In the process of fatty acid binding protein (FAP) evolution to CHIL and CHI, the physical and chemical properties of proteins also had significant differences in addition to changes in protein functions.