Next Article in Journal
Integrated Bioinformatics Analysis Confirms the Diagnostic Value of Nourin-Dependent miR-137 and miR-106b in Unstable Angina Patients
Next Article in Special Issue
Special Issue: Cotton Molecular Genomics and Genetics 2.0
Previous Article in Journal
Systemic Effects of Homoarginine Supplementation on Arginine Metabolizing Enzymes in Rats with Heart Failure with Preserved Ejection Fraction
Previous Article in Special Issue
Comprehensive Evaluation and Transcriptome Analysis Reveal the Salt Tolerance Mechanism in Semi-Wild Cotton (Gossypium purpurascens)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
College of Agronomy, Xinjiang Agricultural University, 311 Nongda East Road, Urumqi 830052, China
2
The State Key Laboratory of Genetic Improvement and Germplasm Innovation of Crop Resistance in Arid Desert Regions (Preparation), No. 403, Nanchang Road, Urumqi 830052, China
3
College of Life Science, Xinjiang Agricultural University, 311 Nongda East Road, Urumqi 830052, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(19), 14775; https://doi.org/10.3390/ijms241914775
Submission received: 15 August 2023 / Revised: 22 September 2023 / Accepted: 25 September 2023 / Published: 30 September 2023
(This article belongs to the Special Issue Cotton Molecular Genomics and Genetics 2.0)

Abstract

:
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.

1. Introduction

Cotton is an economically important crop in China, and cotton production accounts for the largest proportion of natural textile fiber production [1]. Gossypium barbadense plants are highly susceptible to Fusarium wilt, the occurrence of which leads to a sharp drop in cotton yield, the deterioration of fiber quality, and severe production failure. Therefore, cotton Fusarium wilt is referred to as cotton “cancer” [2]. Fusarium wilt is caused by Fusarium oxysporum f. sp. vasinfectum (Fov), a global soil-borne pathogenic fungus that can infect more than 100 kinds of plants [3].
To respond to infection by various pathogens, plants have evolved two main innate immune mechanisms [4]: pathogen-associated molecular pattern (PAMP)-triggered immunity (PT1) and effector-triggered immunity (ET1) [5,6]. Pathogens can complete invasion and colonization processes in the host by secreting effector molecules to inhibit the defense signaling pathway [7] or by directly entering the host [8]. In addition, terrestrial plants have evolved many specific metabolites, including a wide variety of flavonoids, to respond to various biotic and abiotic stresses [9]. Relevant studies have shown that the chalcone isomerase (CHI) gene plays an important role in plant disease resistance. In soybean, the overexpression of the CHI gene enhanced resistance to Phytophthora sojae [10]. In cotton, CHI genes have been predicted to be regulated by miRNAs in response to pathogen challenge [11]. In addition, the GbCHI protein inhibits spore germination and mycelial growth of Verticillium dahliae [12].
Studies have shown that flavonoids are secondary metabolites involved in plant development and defense. They have received extensive attention due to their antioxidant, free radical scavenging, and antibacterial functions, which are of great significance to human health [13]. CHI is ubiquitous in higher plants and is the second most important enzyme in the flavonoid synthesis pathway. Chalcone in plants is converted to dihydroflavones through the action of CHI, and then further converted to various flavanones and flavonoids. The CHI gene family can be divided into four categories [14]. Type I CHI proteins can use only 6-hydroxychalcone (naringenin chalcone) as a substrate in vascular plants to generate 5-hydroxyflavanone (2s-naringenin); type II CHI proteins not only have the catalytic function of type I proteins, but can also use 6′-deoxychalcone as a substrate to generate 5-deoxyflavanone (liquiritigenin). Moreover, the type II proteins are found only in legumes [15]. Type III CHI proteins localize to plastids and are the sites of new fatty acid (FA) biosynthesis in plant cells. In vivo, the expression profile of type III CHI proteins is associated with core FA biosynthetic enzymes. In vitro, these proteins are FA-binding proteins (FAPs) [16]. Notably, FAPs represent a specific period in plant evolution, that is, the adaptive evolution of a nonenzymatic ancestor to a stereospecific and catalytically perfected enzyme [17]. Type IV CHI proteins are found only in terrestrial plants and lack the conventional CHI activity. Previously, it was found that the CHI-like protein EFP is a type IV CHI and is involved in the early stage of the flavonoid biosynthesis pathway of Japanese morning glory, playing a role in ensuring the production of flavonoids. Its function is conserved in different terrestrial plant species, and it is a catalyst for the synthesis of flavonoids and anthocyanins [18,19,20]. Therefore, the classic type I protein and the legume-specific type II CHI are truly catalytically active CHIs that can efficiently produce naringenin.
Gene duplication is one of the core causes for the existence of multiple gene copies. Gene duplication can lead to the formation of families with new functions that may contribute to speciation [21,22,23,24], and the types of gene duplication include whole-genome duplication (WGD), tandem and segmental duplication, transposon-mediated duplication, and reverse transcription-based replication [23,25]. The most common evolutionary model is that one of the repeats retains the function of the original gene, while the other repeat is a pseudogene [25,26]. Alternatively, both copies coexist and slowly evolve over a long period, which provides functional novelty in genetic resources for the evolution of species, that is, neofunctionalization [22].
To date, the relationship between the functional evolution of CHI family genes and Fusarium wilt in G. barbadense has rarely been reported. According to the identification of the cotton CHI gene family, there are 10 CHI genes in G. barbadense, as shown in our previous studies [27]. Furthermore, some homologous genes of the cotton CHI gene family may have been involved in the Fusarium wilt resistance pathway of G. barbadense during evolution. Over the course of evolution, the cotton CHI genes may have assumed the role of disease resistance genes in G. barbadense. Our research further clarifies the molecular specificity and function of cotton CHI genes and their potential application value in cotton genetic improvement. In addition, the relationships of several CHI homolog genes with the same function in the process of genome evolution and the basis for genetic redundancy and phenotype robustness were further clarified.
At present, there are few reports on the synergistic relationship and possible neofunctionalization differences of paralogous CHI family paralogous genes involved in Fusarium wilt resistance in cotton during evolution. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs; two resistant offspring and two susceptible offspring) and the parents after Fov infection showed that CHI genes have important functions in cotton [27]. In addition, four paralogous genes in the cotton CHI gene family were found to be evolutionarily involved in the resistance to Fusarium wilt in G. barbadense. Over the course of evolution, the four paralogs may have assumed the role of disease resistance genes in G. barbadense [27]. In this study, we identified the molecular specificity and function of the four paralogous genes in G. barbadense and found that these genes function synergistically. This study reveals that homologous genes are synergistically regulated in the homeostasis of secondary metabolites due to gene duplication, and the antibacterial effect of flavonoids can be exerted through the “one major, three secondary” pattern of the GbCHI gene family (one major gene, GbCHI05, and three minor genes GbCHI01, GbCHI06 and GbCHI09). This study provides evidence for genetic redundancy and phenotypic robustness. In addition, this study offers a strategy for ensuring plant stability to maintain normal growth and development in the face of genetic and environmental changes.

2. Results

2.1. Cloning and Bioinformatics Analysis of the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 Genes

The open reading frame (ORF) sequences of GbCHI01, GbCHI05, GbCHI06 and GbCHI09, with sizes of 618 bp, 684 bp, 618 bp and 684 bp, respectively, were cloned (Figure S1). The proteins encoded by GbCHI01, GbCHI05, GbCHI06 and GbCHI09 were analyzed by the biological software program ProtParam on the ExPASy website. The four proteins were all hydrophilic proteins. The physicochemical properties of paralogs on the same chromosome were similar, while those of paralogs on different chromosomes were completely different (Table S4). The tertiary structures of the proteins encoded by the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes of G. barbadense were analyzed, and the α-helix and β-sheet structures were folded to determine the correctly folded protein structure that could perform the biological functions of the genes (Figure S2). The overall structure of the GbCHI01 and GbCHI06 proteins resembled an upside-down bouquet, while the overall structure of the GbCHI05 and GbCHI09 proteins resembled an upright bouquet. The two bouquets were similar to each other, differing by only a 180-degree horizontal rotation. Motif Scan software (https://swissmodel.expasy.org/ (accessed on 14 August 2023)) was used to analyze the motifs and domains of the GbCHI01, GbCHI05, GbCHI06, and GbCHI09 proteins. The four proteins had casein kinase II phosphorylation sites, N-myristoylation sites, protein kinase C phosphorylation sites and other functional motifs. GbCHI01 and GbCHI06 had cAMP- and cGMP-dependent protein kinase phosphorylation sites, tyrosine kinase phosphorylation sites and BRCA2 repeat motif functions, and GbCHI05 and GbCHI09 had N-glycosylation sites, an ATP/GTP-binding site motif A (P-loop) and a Big-1 (bacterial Ig-like domain 1) domain profile with primitive functions (Table S5). The GbCHI01, GbCHI05, GbCHI06 and GbCHI09 protein sequences were aligned using multiple sequence alignment by CLUSTALW and ESPript3.0 online software, and online analysis with Swiss Model software delineated the secondary structure of the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 proteins. The GbCHI01 and GbCHI06 proteins had eight α-helices and 11 β-sheets. The GbCHI05 and GbCHI09 proteins had nine α-helices and 11 β-sheets (Figure S3).

2.2. Subcellular Localization of the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 Genes

The GbCHI01, GbCHI05, GbCHI06 and GbCHI09 proteins may be localized to the nucleus and cell membrane (Figure 1). Based on the structural and functional analysis of CHI proteins described above, it is speculated that the nucleus and cell membrane are the main sites where CHI proteins function. The paralogous genes have the same subcellular localization and may perform the same function.

2.3. Silencing of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 Reduces Resistance to Fusarium Wilt in Cotton

The GbCHI01, GbCHI05, GbCHI06, and GbCHI09 genes of G. barbadense 06-146 were silenced by virus-induced gene silencing (VIGS) technology, and cloroplastos alterados 1 (CLA1) was silenced as a positive control. Positive control (pTRV2::CLA1) cotton plants exhibited an albino phenotype after injection (Figure 2). The silencing efficiency measurements show that the expression levels of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 were significantly decreased (Figure 2), indicating that GbCHI01, GbCHI05, GbCHI06 and GbCHI09 were successfully silenced. The experimental group with silenced GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) and the control group (pTRV2::00) were observed 28 days after infection with Fusarium wilt.
The plants of the experimental group and the control group exhibited leaf yellowing, wilting and defoliation, and the plants of the experimental group exhibited more severe yellowing and wilting than those of the control group (Figure 2). Statistical analysis of the disease index (DI) showed that the DI of the experimental group was significantly higher than that of the control group. In conclusion, silencing GbCHI01, GbCHI05, GbCHI06 and GbCHI09 in G. barbadense 06-146 significantly reduced the resistance of cotton seedlings to Fusarium wilt (Figure 2).
The results of the recovery test of Fusarium wilt in stem segments show that in the same parts of cotton stem segments, the control plants developed fewer Fusarium wilt stem segments than the gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) (Figure 3). The gene-silenced plants displayed much higher Fov race 7 biomass. These results indicate that the gene-silenced plants had a greater extent of cotton Fov infection, and the ability of the gene-silenced plants to resist cotton Fusarium wilt was weaker than that of the control plants. The ability of cotton seedlings to resist cotton Fusarium wilt decreased with the silencing of the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes (Figure 3).

2.4. Overexpression of GbCHI05 Improves Resistance to Fusarium Wilt in Arabidopsis thaliana and Cotton

After the overexpression of the GbCHI05 gene in T3 Arabidopsis thaliana seeds, the wild type (WT) and chi05 mutant were cultured in a greenhouse. The wild type (WT) and chi05 mutant began to show yellow margins and curled leaves 14 days after Fov infection, while the transgenic A. thaliana plants showed good leaf growth (Figure 4). At 21 days of infection, WT A. thaliana and the chi05 mutant were found to have more severe disease than the GbCHI05-overexpressing strains (Figure 4). Statistical analysis of the DI showed that the DIs of the WT and chi05 mutant were significantly higher than those of the GbCHI05-overexpressing strains. The results show that the GbCHI05 gene improved the resistance of A. thaliana to Fusarium wilt.
The results of the Fusarium wilt recovery assay show that the control plants developed fewer Fusarium wilt-infected stem segments than the GbCHI05 mutant A. thaliana plants (Figure 4). The GbCHI05 mutant A. thaliana plants displayed much higher Fov race 7 biomass. The GbCHI05-overexpressing A. thaliana plants developed fewer Fusarium wilt-infected stem segments than the control plants (Figure 4). The GbCHI05-overexpressing A. thaliana plants displayed much higher Fov race 7 biomass. These results indicate that the GbCHI05 mutant A. thaliana had a greater extent of Fov infection, and that the GbCHI05 mutant A. thaliana had a weaker ability to resist Fusarium wilt than the control plants (Figure 4). The GbCHI05-overexpressing A. thaliana plants were less susceptible to Fusarium wilt than the control plants, exhibiting a greater ability to resist cotton Fusarium wilt (Figure 4).
After the overexpression of the GbCHI05 gene in T3 06-146 seeds, the 06-146 (WT) seeds were cultured in the greenhouse. The expression level of the GbCHI05 gene in GbCHI05-overexpressing 06-146 was significantly higher than that in the recipient materials, indicating that the GbCHI05 gene was overexpressed in 06-146 (Figure 5). After 28 days of infection, the disease in control plants was more severe than that in GbCHI05-overexpressing strains (Figure 5). The results indicate that cotton seedlings overexpressing the GbCHI05 gene were resistant to cotton Fusarium wilt. The DI was statistically analyzed, and the DI of the control group was found to be higher than that of the experimental group. The GbCHI05-overexpressing cotton plants developed fewer Fusarium wilt-infected stem segments than the control plants (Figure 5). The GbCHI05-overexpressing cotton plants displayed much higher Fov race 7 biomass. The overexpression of the GbCHI05 gene in 06-146 enhanced the resistance of the cotton seedlings to Fusarium wilt.

2.5. Detection of Genes in Flavonoid Metabolic Pathways and Genes Related to Fusarium Wilt

To study the effects of silencing the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes on the genes in the flavonoid metabolic pathway, the expression levels of synthetic genes in the flavonoid metabolic pathway were determined. At the same time, the expression levels of the previously reported GbERF-like and Gbar_D03G002290 genes were also measured and analyzed. The study found that when GbCHI01 was silenced, the expression levels of GbCHI05, GbCHI06, GbCHI09, GbDFR, GbF3H, GbFLS, GbANR and GbANS were upregulated to varying degrees, the expression of GbC4H was decreased, the expression of the GbERF-like gene was upregulated, and the expression of Gbar_D03G002290 was downregulated (Figure 6). When GbCHI05 was silenced, the expression levels of GbCHI01, GbCHI06, GbCHI09, GbDFR, GbF3H, GbFLS, and GbANS were all downregulated to varying degrees (Figure 6). The expression levels of GbANR and GbC4H first decreased, then increased and then decreased again, and GbERF-like expression was upregulated. Gbar_D03G002290 expression was downregulated first and then upregulated. When GbCHI06 was silenced, the expression levels of GbCHI06, GbCHI09, GbDFR, GbF3H and GbFLS were upregulated to different degrees, and the expression levels of GbCHI01, GbANR, GbANS, GbC4H, GbERF-like and Gbar_D03G002290 were all downregulated to different degrees (Figure 6). When the GbCHI09 gene was silenced, the expression levels of GbCHI01, GbCHI05, GbCHI06, GbDFR, GbANS, GbC4H, GbERF-like and Gbar_D03G002290 were all downregulated to varying degrees, and the expression levels of GbF3H, GbFLS and GbANR were upregulated to varying degrees (Figure 6). The results show that the four homologous genes had similar regulatory genes under the conditions of Fov infection, but the differences in regulation were unique to each gene involved in the regulation of flavonoid metabolic flux.

2.6. Flavonoid Content after Silencing of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 and GbCHI05 Overexpression

The following linear regression equation was established according to the standard curve of rutin (Figure S4): y = 0.1322x + 0.0968, R² = 0.9925. The flavonoid content was calculated as follows: (A508 − 0.0968)/0.1322 × 100. Statistical analysis of the flavonoid content showed that the flavonoid content in control plants was higher than that in experimental plants without Fov infection, and the flavonoid content in control plants was lower than that in experimental plants without Fov infection (Figure 7). The study results show that the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes are involved in the synthesis of flavonoids, but they are not the only genes that control the synthesis of flavonoids (Figure 7). The synthesis of flavonoids is regulated by multiple genes. After GbCHI05 gene overexpression, the levels of flavanones in GbCHI05-overexpressing cotton seedlings were higher than those in the recipient material 06-146 (Figure S5).

2.7. Levels of Endogenous Hormones and ROS in GbCHI01, GbCHI05, GbCHI06 and GbCHI09 Gene-Silenced Plants after Fov Infection

The levels of endogenous MeJA and SA in the samples were determined by ELISA, and it was found that with the prolongation of infection time, the levels of MeJA and SA increased in both the control group and experimental group (Figure 8). In addition, in the GbCHI01, GbCHI06 and GbCHI09 gene-silenced plants, the endogenous MeJA content in the samples of the experimental group was lower than that in the samples of the control group, while the endogenous SA content in the samples of the experimental group was lower than that in the samples of the control group (Figure 8). In GbCHI05 gene-silenced plants, the endogenous MeJA content in the experimental group samples was lower than that in the control group samples. However, the amount of endogenous SA in the experimental group samples was lower than that in the control group samples, except at 2 h and 4 h (Figure 8).
The levels of ROS in the samples were determined by ELISA. The ROS content in the test group was lower than that in the control group in GbCHI01 and GbCHI05 gene-silenced plants. In GbCHI06 and GbCHI09 gene-silenced plants, the ROS content in the experimental group was higher than that in the control group (Figure 9). In GbCHI01 and GbCHI06 gene-silenced plants, the ROS content in the samples before infection was higher than that in the samples after infection. In GbCHI05 and GbCHI09 gene-silenced plants, the ROS content in the samples before infection was lower than that in the samples after infection (Figure 9).

3. Discussion

During plant development and evolution, genes control biological processes through precise spatial and temporal expression. Genes regulate biological processes and respond to environmental changes through changes in expression. Therefore, evolutionary adaptation and changes in genes are of great significance for the ability of individual plants to adapt to environmental changes over tens of millions of years [28,29]. The copy number changes caused by gene duplication and evolution and the differences in copy structure and function are important contributors to the construction and improvement of gene regulatory networks [30,31,32,33,34]. In addition, gene duplication leading to redundancy in gene segments may help species acquire adaptation and genetic resistance to changes in environmental conditions [35].

3.1. Effect of Gene Structural Variations on the Neofunctionalization Differentiation of Homologous Genes and the Synergistic Functions of Homologous Genes

Different genomic structural variations (insertions or deletions) can lead to the reorganization of topologically associated domains (TADs). Cotton species-specific TAD boundary formation occurs more frequently in inactive chromatin regions, which is related to cotton species-specific transposon insertion. At the same time, the chromatin interactions between genes are more conserved than those between noncoding regulatory regions and genes among different cotton species. It was found that structural variation, transposon insertion and changes in chromatin openness drove the evolution of regulatory networks between noncoding regulatory elements and homologous genes, providing new ideas for the study of structural variation and function in cotton [36].
We investigated the effect of gene structural variation on the neofunctionalization of homologous genes and the synergistic function of homologous genes during GbCHI gene evolution, and found that the sequences of paralogous CHI genes in different chromosomal subgroups were highly conserved, with only a few base changes. In addition, the three-dimensional structure prediction showed that the proteins encoded by the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes were similar in spatial structure (Figure S2). The protein sequences of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 have different amino acids at individual sites (Table S5), which may lead to differences in their secondary structures, leading to differences in activation or enzymatic reaction rates (Figure S3).
As the main executors of gene functions, proteins need to be located in the correct subcellular compartments to function normally, maintain orderly and efficient progression of various complex biochemical processes in the organism and ensure the normal life activities of the organism [37]. Therefore, protein subcellular localization is closely related to protein function, which is the basis for studying the biological function of proteins and has important significance for proteomics [38]. The function of a gene and the metabolic pathways in which it participates are largely related to the location of the protein encoded by the gene in the organelle. Studying the subcellular localization of a gene can enable identification of the specific location of protein expression in the cell, which allows speculation around its function and the metabolic pathways involved [39]. In this study, we found that the nucleus and cell membrane might be the main functional sites of the CHI protein, which was similar to the localization of CHI proteins in Ipomoea batatas (Linn.) Lam. [40] and Vitis vinifera L. [41], indicating that CHI may play essential roles as an enzyme and in secondary metabolite synthesis during cell reproduction (Figure 1). Research on the subcellular localization of the protein encoded by the CHI gene provides a theoretical basis for further analysis of the function of the gene and screening of its interacting proteins.

3.2. Relationship between the Neofunctionalization of Homologous Genes and Flavonoid Metabolic Pathway Genes

Both whole-genome and small-scale duplications are referred to as gene duplications and are ubiquitous in plant genomes [25,42,43,44,45]. Many redundant repeat segments are generated after gene duplication, which makes genetic variation abundantly selective [25,43]. Genetic variation can be accomplished by pseudogenization, subfunctionalization or neofunctionalization [25,43,46,47]. In addition, paralog evolution by “active compensation”, in which one or more paralogs in the genome are upregulated by transcription to replace the impaired activity of other paralogs, manifests in the genome in the form of redundancy [48,49,50]. This approach offers the possibility for species to maintain stability during genetic changes and under changes in the external environment, and allows genes to remain stably inherited under selection pressure [51,52]. However, although duplication can provide many redundant fragments that allow species to cope with impaired gene activity, it also promotes genetic variation [51,53,54]. This genetic variation can occur simultaneously in the coding and regulatory region sequences of genes, and is the source of the diversity of paralogous genes. Species produce duplications for short periods, and the variation caused by this short period of duplication-induced diversity in the function of homologous genes is unknown. To date, functional studies of homologous genes have been conducted only within a single system or in a few species, and the compensatory relationships and functional changes between genes have not been systematically studied [44,48,55]. Therefore, studying the compensatory relationships of paralogous genes in Fusarium wilt resistance and their differences after neofunctionalization in G. barbadense will provide evidence for the evolutionary model after lineage gene duplication.
Our study showed that silencing GbCHI01, GbCHI05, GbCHI06 and GbCHI09 in G. barbadense 06-146 reduced the resistance of cotton seedlings to Fusarium wilt (Figure 2 and Figure 3). The overexpression of the GbCHI05 gene in 06-146 and A. thaliana enhanced the resistance of cotton and A. thaliana seedlings to Fusarium wilt (Figure 4, Figure 5 and Figure 6). Therefore, the roles of these paralogous genes are not equivalent; instead, these genes work together to achieve robust genetic channelization. To study the effect of silencing the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes on the genes in the flavonoid metabolic pathway, the expression levels of synthetic genes in the metabolic pathway were measured (Figure 7).
The experiment showed that under the conditions of Fov infection, the silencing of the homologous genes led to similar changes in the expression levels of the same genes, but the expression levels had opposite trends, and some genes exhibited different expression levels due to gene silencing. The expression patterns of GbDFR, GbF3H and GbFLS were the same when GbCHI01, GbCHI06 and GbCHI09 were silenced, respectively. The expression patterns of GbANS were different when GbCHI05 was silenced compared to when other genes were silenced. In GbCHI01 gene-silenced plants, the expression pattern of GbANR was different from that in other gene-silenced plants. This indicates that the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes have similar regulatory networks but may also have their own unique regulatory pathways. These genes may cooperate with and restrict each other, and jointly affect the synthesis and functions of secondary metabolites. However, the differences in their regulation may be the key to the divergence of the neofunctionalization of genes and the maintenance of gene dosage balance.

3.3. Relationship between Neofunctionalization of Differentiation of Homologous Genes and Flavonoid Content

Flavonoids are the largest class of polyphenols, and the initial substrates for plant flavonoid synthesis are coumaroyl-CoA and malonyl-CoA, derived from the phenylpropane metabolic pathway [56,57]. First, phenylalanine produces colorless naringenin under the action of phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumaric acid-CoA ligase (4CL), chalcone synthase (CHS) and CHI. Naringenin, as a major metabolite, enters the synthetic pathways of other flavonoids. Since the synthetic pathways of flavonoids are conserved in plants, many enzymes can change the basic skeleton of flavonoids under different external conditions, thereby generating different kinds of flavonoids [58,59].
The CHI gene was first isolated from French pea (Pisum sativum), and it has subsequently been isolated and cloned from various species [60,61,62,63,64,65,66,67]. Whether CHI is expressed in plants, and the level of expression, affect the metabolism of flavonoids in plants. Relevant studies have shown that when CHI is inactive or mutated, the chalcone content in plants increases, and the flavonoid content is significantly reduced compared with that in the presence of wild type CHI. Decreased CHI activity in carnation (Dianthus caryophyllus Linn.) [68] and cyclamen (Cyclamen persicum Mill.) [69] resulted in a change in the color of the flowers, leading to the production of yellow flowers. In addition, after the expression of the CHI gene in carnation (D. caryophyllus Linn.) was reduced, chalcone accumulated in large quantities, and yellow flowers were produced. Similarly, the mutation of CHI changed the flower color of carnations to yellow [70]. The mutation of the CHI gene in petunia (Petunia hybrida Vilmorin.) led to the formation of yellow or green pollen due to the accumulation of chalcone [60]. After the RNAi-mediated inhibition of CHI gene expression in tobacco (Nicotiana tabacum L.), the anthocyanin content in the petals decreased, the chalcone content increased, and the petals turned yellow [71]. Relevant studies have shown that the flavonoid content in plants can be significantly increased by increasing the activity of the CHI enzyme or overexpressing the associated gene. For example, the petunia CHIA gene was introduced into tomato peel, and the flavonol content was found to increase 79-fold [72]. The jellyfish snow lotus CHI gene was transferred into Xinjiang snow lotus (Saussurea involucrata Kar.) and overexpressed. The apigenin content in transgenic Xinjiang snow lotus with the CHI gene was 12 times higher than that in the control group, and the total flavonoid content increased 4-fold [73]. The total flavonoid content in transgenic tobacco was 5 times higher than that in the wild type when the jellyfish yalus CHI gene was transferred into tobacco [74]. Studies have shown that overexpression of the CHI gene can also improve the antioxidant capacity of potato (Solanum tuberosum L.) [75]. Flavonoids not only have antibacterial effects, but also are associated with a variety of other biological processes. The results of previous studies on the main synthases in the flavonoid metabolism pathway are relatively clear, but the regulatory network of the flavonoid metabolism pathway has a certain complexity. This network plays an important role in various biological processes, and it is very important to study the key rate-limiting enzyme-encoding genes and homologous gene regulatory networks to elucidate the role of metabolites in the biological processes.
Our experimental results show that the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes are all involved in the synthesis of flavonoids, and they are also involved in the defense provided by flavonoids via antibacterial activity in G. barbadense. The synthesis of flavonoid antibacterial substances is catalyzed and regulated by the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes (Figure 7). After the overexpression of the GbCHI05 gene in 06-146, the flavanones content in the transgenic strains was higher than that in the recipient materials (Figure S5). Flavanones, as intermediate metabolites [76], inhibit Fusarium wilt by feedback regulation and the synthesis of other flavonoid compounds downstream. These results suggest that GbCHI05 gene silencing might rapidly stabilize the metabolic flux of flavonoids through feedback regulation or the rapid activation of other homologous genes to confer defense against Fusarium wilt in G. barbadense.
Theoretical models of gene functional differentiation after whole-genome duplication involve gene dosage balance [25], subfunctionalization [25,77,78,79] and neofunctionalization [25,45]. Among them, neofunctionalization can occur via numerous different biological processes, such as the acquisition of different transcriptional regulation pathways [80], the subcellular relocalization of proteins [81,82,83] and the acquisition of different protein targets [84,85,86]. However, when the original function of the gene is indispensable, how does neofunctionalization occur? This poorly understood and curious phenomenon is called Ohno’s dilemma [87].
Given that the four genes have the same subcellular localization, the whole-genome duplication-based evolution of the four homologous genes may be explained by the gene dosage balance theory. The four genes cooperate to maintain the stability of flavonoid metabolic flux and synthesize downstream secondary metabolites in an orderly manner, according to the gene dosage balance theory, thus achieving Fusarium wilt resistance. The four homologous genes play an important role in jointly maintaining the dynamic balance of flavonoid synthesis (Figure 8). Through the synergistic regulation of the four genes, the synthetic small-molecular compounds associated with disease resistance are released one at a time, avoiding the one-time release of antibacterial substances, which would affect the response of plants to other stresses and lead to excessive immune damage in the plants. During the process of selection and evolution in plants, the four genes of the CHI gene family have specific functions, acting as a cooperative group.

3.4. Effect of Endogenous Hormones and ROS on the Neofunctionalization of Homologous Genes

Three phytohormones, MeJA, SA and ethylene (ET), play important regulatory roles in the plant’s defense against pathogens and pests and in the response to abiotic stresses. JA and ET are usually involved in the defense response against necrotizing pathogens. Generally, the MeJA and SA defense pathways are antagonistic to each other, and they act synergistically with the ET defense pathway [88]. Studies have shown that the concentration of JA increases after pathogen infection or tissue damage, and the exogenous application of JA induces the expression of defense-related genes [89]. SA is commonly associated with the defense against biotrophic and semibiotrophic pathogens, as well as the establishment of systemic acquired resistance [90].
Recent studies have shown that JA signaling plays an important role in systemic immunity in A. thaliana L. by participating in long-distance information transmission [91]. The rapid accumulation of and increase in JA levels in leaf phloem exudates, and systemwide in leaves, suggests that JA may play a role as a mobile signal involved in pathogen immunity in A. thaliana L. [91]. Research on the SA content in Verticillium wilt-resistant G. barbadense varieties and Verticillium wilt-susceptible upland cotton varieties have has that the SA content in Verticillium wilt-resistant G. barbadense varieties before and after inoculation was significantly higher than that in Verticillium wilt-susceptible upland cotton varieties [92]. Recent studies have shown that methyl salicylate production is induced during pathogen infection, and this compound acts as a flow inducer in systemic acquired resistance in tobacco [93].
This study showed that silencing the GbCHI01, GbCHI06 and GbCHI09 genes reduced the accumulation of endogenous MeJA and increased the accumulation of endogenous SA (Figure 9). However, after the GbCHI05 gene was silenced, the accumulation of endogenous MeJA and endogenous SA was reduced. These four genes might regulate the metabolic flux of flavonoids via the regulation of the levels of endogenous MeJA and SA and simultaneously trigger the immune mechanism of G. barbadense against Fusarium wilt (Figure 9). After the GbCHI05 gene was silenced, the changes in the accumulation of endogenous SA were different from those after the silencing of the other genes. These results indicate that the change in the SA content in plants after GbCHI05 gene silencing may be the key signal or factor that can stabilize the metabolic flux of flavonoid synthesis in plants. In the early stage of resistance to Fusarium wilt, the flavonoid metabolic pathway may be mainly mediated by MeJA, and in the later stage, it may be mainly mediated by SA. In addition, MeJA and SA do not act independently in the resistance of G. barbadense to Fusarium wilt. The regulatory network associated with MeJA and SA regulation in Fusarium wilt resistance in G. barbadense and the underlying mechanisms remain to be studied. However, MeJA and SA, as signals, play an important role in the mechanisms by which flavonoid metabolic pathway structural genes regulate Fusarium wilt resistance in G. barbadense (Figure 9).
ROS are key signaling molecules that enable cells to respond quickly to different stimuli. In plants, ROS play a crucial role in the perception of abiotic stress and biological stress, the integration of different environmental signals, and the activation of stress response networks. In addition, ROS play a key role in the activation and training of plant defense mechanisms [94]. ROS metabolism-related genes also share a certain relationship with fungal pathogenicity. ROS accumulation negatively regulates the pathogenicity of Fusarium oxysporum f. sp. momordicae [95]. Research data show that QSOX1 in plants is a REDOX sensor that restricts the production of ROS by connecting ROS and active nitrogen signals, and negatively regulates plant immunity [96]. ROS can directly inhibit the growth of pathogens and participate in the disease resistance process as signaling molecules. The rapid production of ROS is an important sign of the activation of the plant’s defense system. After TaNOX10 gene knockout in Triticum aestivum L., the extracellular ROS concentration in T. aestivum L. was significantly reduced, which confirmed that TaNOX10 participated in T. aestivum L. resistance to Puccinia striiformis by promoting ROS production. Genetic experiments have confirmed that TaWRKY19 can specifically bind to W-box elements in the TaNOX10 promoter and transcriptionally inhibit the expression of the TaNOX10 gene [97]. When cotton plants sense infection by pathogens such as Verticillium dahliae, GhCaM7 binds to Ca2+ as a secondary messenger to enhance the defense response of cotton by activating the JA, ROS and MAPK pathways and changing the cell osmotic potential [98]. In this study, silencing the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes directly affected ROS accumulation, indicating that the four CHI family genes may directly confer resistance to Fusarium wilt and activate other defense systems of cotton by regulating ROS metabolism (Figure 10).
Gene and genome duplications occur very frequently during plant evolution, and paralogous genes generated by gene duplication during genome evolution provide the basis for genetic redundancy and phenotypic robustness [99]. Gene duplication is a hallmark of plant genome evolution and underlies the genetic interactions that shape phenotypic diversity. Whole-genome duplication events duplicate an entire pathway or network, and the genes in the duplicated pathway or network can evolve new functions through the neofunctionalization and subfunctionalization of genes [100]. Compensation is a major form of homologous gene interaction, with many genes acting as “backup copies” of the original genes. However, how the compensatory relationship changes as allelic variations accumulate is unclear. This makes the result of gene editing less predictable, adding to the difficulty of molecular breeding-based crop improvement efforts [101]. In this study, four genes of the CHI gene family related to Fusarium wilt resistance in G. barbadense were identified, and how they might act cooperatively to confer resistance to Fusarium wilt was explored. The results of this paper reveal a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous SA and MeJA by paralogous genes, thereby achieving disease resistance.
Due to the coregulation of paralogous genes in the homeostasis of secondary metabolites due to gene duplication, the antibacterial effects of flavonoids can be exerted through the “one major, three secondary” mode in the CHI gene family (one main gene, GbCHI05, and three minor genes, GbCHI01, GbCHI06 and GbCHI09) (Figure 10). This study demonstrates the ability of plants to maintain normal growth and development in the face of genetic and environmental changes, and provides a theoretical basis for studying the evolutionary patterns of plant homologous genes and for using homologous genes for molecular breeding.

4. Materials and Methods

4.1. Cloning of GbCHI01, GbCHI05, GbCHI06, and GbCHI09 and Bioinformatics Analysis

The G. barbadense cultivar used in the test was a resistant cultivar (06-146). Total RNA from cotton was extracted with the Plant Polysaccharide Polyphenol RNA Extraction Kit (TianGen, Beijing, China), and cotton cDNA was synthesized with the Reverse Transcription Kit (Applied Biological Materials, Vancouver, CA, USA). According to the gene sequences of GB_A05G4332 (GbCHI01), GB_A13G0219 (GbCHI05), GB_D04G0136 (GbCHI06) and GB_D13G0209 (GbCHI09) [27], Premier 5 software was used to design specific amplification primers for the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 genes (Table S1). The amplification program was as follows: 98 °C for 5 min; 35 cycles of 98 °C for 30 s, 57 °C for 30 s, and 72 °C for 1 min for 40 s; and 72 °C for 10 min. The fragments of G. barbadense GB_A05G4332 (GbCHI01), GB_A13G0219 (GbCHI05), GB_D04G0136 (GbCHI06) and GB_D13G0209 (GbCHI09) were amplified and inserted in the pLB-Simple Vector according to the instructions of the Lethal Based Simple Fast Cloning Kit (Tiangen, Beijing, China).
Properties such as the instability index and overall average hydrophilicity were predicted using the online tool ProtParam. PHYRE2 was used to predict the tertiary structure of the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 proteins. Motif Scan software was used to analyze the motifs and domains of GbCHI family proteins. Using MEGA 7.0 and ESPript 2.2 software, the amino acid sequences of GbCHI family proteins were aligned, and the secondary structure of the CHI protein was described according to the tertiary structure model predicted by the Swiss Model.

4.2. Subcellular Localization of the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 Genes

The GbCHI01, GbCHI05, GbCHI06, and GbCHI09 genes were inserted into subcellular localization vectors with pEarlryGate101 primers as described in Table S2, and the recombinant plasmids were transformed into Agrobacterium tumefaciens GV3101 and transiently expressed in 4-week-old tobacco for 48 h [102]. YFP fluorescence signals were examined for transient expression analysis.

4.3. Construction of the VIGS Vector and Fov Infection of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 Gene-Silenced Plants

The pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09 vectors were designed and constructed using infusion technology. The vector primers are shown in Table S2. Each vector was transformed into A. tumefaciens GV3101. The G. barbadense cultivar used in the test was a resistant cultivar (06-146). The gene silencing vectors induced by tobacco mosaic virus were pTRV1 and pTRV2, and the positive control vector for silencing was pTRV2::CLA1.
Cotton seedlings (06-146) were planted in a cotton culture room with a photoperiod of 16 h (light)/8 h (dark), a temperature of 25 °C, and a relative humidity of 60–70% for 8 days. Cotton seedlings with fully expanded cotyledons were injected with the spore solution, and then the injected cotton seedlings were cultured in the dark for 24 h in a cotton culture room [103]. After this cultivation period, the cotton seedlings were planted under suitable conditions. After the above two weeks, injected cotton seedlings were taken from the experimental group and the control group to measure the silencing efficiency. At the three-leaf stage of cotton seedlings with GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene silencing, the experimental group and the control group were infected with Fusarium wilt [27]. The Fusarium strain used was Fov race 7 (a spore solution of approximately 107 spores/mL) [27]. The experiment was repeated 3 times, and at least 30 cotton seedlings were used in each experiment. The statistical analysis of disease symptoms was performed according to the five-level standard [104]. Control plants (pTRV2::00) and GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) were infected with Fusarium wilt. Afterward, the stem segment 3 mm above the cotyledon node was taken for the Fusarium wilt recovery assay. The stems were surface-disinfected with 0.5% sodium hypochlorite (NaClO), 70% ethanol and sterile distilled water and then cultivated on solid potato dextrose agar (PDA) solid for 3 days at 28 °C [105].

4.4. Construction of the Overexpression Vector and Fov Infection of GbCHI05-Overexpressing Plants

The pCAMBIA3301::GbCHI05 vector was designed and constructed using infusion technology. The primers are shown in Table S2. The pCAMBIA3301::GbCHI05 vector was transformed into A. tumefaciens GV3101. The GbCHI05 gene was overexpressed in Columbia type 0 A. thaliana and cotton (06-146) by the floral dip method and pollen tube channel method. Positive A. thaliana seedlings were screened with herbicides at 200 mL/L until T3 seeds were obtained. The GbCHI05 mutants were purchased from the AraShare Technology Service Center. Positive cotton seedlings were screened with herbicides at 90 parts per million (ppm) and qRT–PCR until T3 seeds were obtained.
The seeds of wild type A. thaliana, GbCHI05 mutants and three GbCHI05 transgenic strains were selected and sterilized with 70% alcohol and 0.7% hypochlorous acid. The seeds were planted in the soil after 3 days of vernalization at 4 °C and then cultured in a greenhouse for 7 days with a photoperiod of 16 h (light)/8 h (dark), a temperature of 20 °C, and a relative humidity of 60–70%. The wild type A. thaliana, GbCHI05 mutants and three GbCHI05 transgenic strains were transplanted into nutrient pots. After 4 weeks, the wild type A. thaliana, GbCHI05 mutants and three GbCHI05 transgenic strains were infected with Fusarium wilt [27], and the phenotype was observed. After infection with Fusarium wilt, control (WT), mutant (chi05) and the GbCHI05-overexpressing A. thaliana plants sampled below the first 3 mm of the stem were used for the Fusarium wilt recovery assay [105]. The details of the protocols for identification of Fusarium wilt resistance and the Fusarium wilt recovery assay of cotton overexpression plants are described in the methods above.

4.5. Genes Associated with Metabolic Pathways of Flavonoids and Fusarium Wilt

After infection with Fusarium wilt, samples were taken at 0, 2, 4, 8, 12, 24, 48, and 72 h. The samples of the experimental group and the control group were quickly frozen in liquid nitrogen and stored at −80 °C. Metabolic pathway-related genes (GbC4H, GbCHS, GbCHI01, GbCHI05, GbCHI06, GbCHI09, GbDFR, GbF3H, GbANR, GbFLS and GbANS) and Fusarium wilt-related genes (GbERF-like and Gbar_D03G002290) were detected. Among them, the GbERF-like gene is involved in the resistance to Fusarium wilt in G. barbadense mediated by SA [30], Gbar_D03G002290 (Gh_D03G0209) is the key gene for resistance to Fusarium wilt [3], and the Gbar_D03G002290 gene is homologous to the Gh_D03G0209 gene in upland cotton. Total RNA of cotton was extracted with a plant polysaccharide and polyphenol RNA extraction kit (Tiangen, Beijing, China), and cotton cDNA was synthesized with a reverse transcription kit (Applied Biological Materials, Vancouver, CA, USA). PCR amplification was performed with a real-time quantitative PCR instrument (Applied Biosystems, Foster City, CA, USA) and the SYBR Green dye method (Applied Biological Materials, Vancouver, CA, USA). The reaction conditions were as follows: 3 min at 95 °C; 40 cycles of 15 s at 95 °C and 1 min at 60 °C. The cotton ubiquitin 7 (UBQ7) gene was used as an internal reference, and the PCR primer sequences are shown in Table S3. The biological experiments were repeated 3 times, and the relative gene expression level was analyzed by the 2−ΔΔCt calculation method [106].

4.6. Extraction and Quantification of Total Flavonoids

An appropriate amount of rutin was weighed and dissolved in 60% ethanol to prepare the rutin standard solution (0.2 mg/mL). Then, 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 mL of the rutin standard solution (0.2 mg/mL) was placed in 10 mL volumetric flasks, and 60% ethanol was added to a constant volume to prepare the rutin standard reaction mixture. The flavonoid content was determined by an ultraviolet–visible spectrophotometer by measuring the absorption peak at 508 nm. Taking the rutin concentration (mg/mL) as the abscissa and the absorbance at 508 nm as the ordinate, a standard curve was drawn, and the linear regression equation (y = ax + b) was obtained. The absorbance at 508 nm of the reaction mixture of each prepared sample was measured, and the total flavonoid content of the unit sample was calculated according to the linear regression equation.
After infection with Fusarium wilt, the stems and leaves above the base of the cotyledons of cotton seedlings in the experimental and control groups were taken as experimental materials at 0 h and 96 h, respectively. The experimental materials were dried in a constant-temperature electric heating box at 60 °C for 2 h, pulverized in a mortar and passed through a 40-mesh sample sieve, and the powder was collected in a sampling Ziploc bag and marked. A total of 0.05 g of the experimental material was accurately weighed and placed into a 10 mL centrifuge tube. After adding 3 mL of 60% ethanol to the centrifuge tube, the centrifuge tube was heated in a 70 °C constant-temperature water bath for 4 h, and the total flavonoids in the experimental material were fully dissolved in 60% ethanol after heating. After centrifugation, the supernatant solution in the centrifuge tube was placed in a 10 mL volumetric flask [107]. In total 2 mL of the extract solution was placed in a 10 mL volumetric flask, and 0.5 mL of 5% sodium nitrite solution and 0.5 mL of 10% aluminum nitrate solution were added in turn. Then, the mixture was evenly mixed and allowed to stand for 6 min. Then, 5 mL of 4% sodium hydroxide solution and 2 mL of 60% ethanol were added and mixed for 15 min to obtain the sample reaction mixture. Each treatment was repeated three times.

4.7. Determination of the Levels of Endogenous Hormones, Flavanone and ROS in Plants

After the experimental group and the control group were infected with Fusarium wilt, samples were taken at 0, 4, 8, 12, 24, 48 and 72 h. The levels of endogenous MeJA and SA in the samples were determined by the enzyme-linked immunosorbent assay (ELISA) in GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants. The assays were performed according to the instructions for the Plant Methyl Jasmonate ELISA Kit (Jiangsu Jingmei Biotechnology, Yancheng, China) and the Plant Salicylic Acid ELISA Kit (Jiangsu Jingmei Biotechnology, Yancheng, China).
After the experimental group and the control group were infected with Fusarium wilt, samples were taken at 0 and 96 h. The levels of ROS in the samples were determined by enzyme-linked immunosorbent assay (ELISA) in GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants. The assays were performed according to the instructions for the Plant Reactive Oxygen Species ELISA Kit (Jiangsu Jingmei Biotechnology, Yancheng, China). The levels of flavanone in the samples were determined by enzyme-linked immunosorbent assay (ELISA) in GbCHI05-overexpressing cotton plants. The assays were performed according to the instructions for the Flavanone ELISA Kit (Jiangsu Jingmei Biotechnology, Yancheng, China).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241914775/s1.

Author Contributions

Conceptualization, X.D., X.L. and Q.C. (Quanjia Chen); methodology, Q.Z., X.D., C.W. and Y.L.; software, X.C., Y.C. and K.Z.; validation, Q.Z. and Q.C. (Qin Chen); formal analysis, Q.Z. and Y.Q.; investigation, X.D., Y.C., Y.L. and K.Z.; resources, Y.Q., X.C., Q.C. (Qin Chen) and X.L.; data curation, Q.Z., X.D. and C.W.; writing—original draft, Q.Z.; writing—review and editing, Q.Z. and Q.C. (Quanjia Chen); supervision, Y.Q.; funding acquisition, Q.C. (Quanjia Chen) and X.D. contributed equally to the paper and were listed as the first authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Major Science and Technology Projects of Xinjiang Uygur Autonomous Region (2021A02001-3).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data generated or analyzed during this study are included in this article and its supplementary data files.

Acknowledgments

The author is especially grateful for the support of the Key Laboratory of Crop Genetic Improvement and Germplasm Innovation of Xinjiang Agricultural University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, J.; Fang, H.; Zhou, H.; Sanogo, S.; Ma, Z. Genetics, breeding, and marker-assisted selection for Verticillium wilt resistance in cotton. Crop. Sci. 2014, 54, 1289–1303. [Google Scholar] [CrossRef]
  2. Cox, K.L.; Babilonia, K.; Wheeler, T.; He, P.; Shan, L. Return of old foes—Recurrence of bacterial blight and Fusarium wilt of cotton. Curr. Opin. Plant Biol. 2019, 50, 95–103. [Google Scholar] [CrossRef]
  3. Liu, S.; Zhang, X.; Xiao, S.; Ma, J.; Shi, W.; Qin, T.; Xi, H.; Nie, X.; You, C.; Xu, Z.; et al. A single-nucleotide mutation in a GLUTAMATE RECEPTOR-LIKE gene confers resistance to Fusarium wilt in Gossypium hirsutum. Adv. Sci. 2021, 8, 2002723. [Google Scholar] [CrossRef] [PubMed]
  4. Jones, D.A.; Takemoto, D. Plant innate immunity—Direct and indirect recognition of general and specific pathogen-associated molecules. Curr. Opin. Immunol. 2004, 16, 48–62. [Google Scholar] [CrossRef]
  5. Jones, J.D.G.; Dangl, J.L. The plant immune system. Nature 2006, 444, 323–329. [Google Scholar] [CrossRef] [PubMed]
  6. Medzhitov, R.; Janeway, C.A. Innate immunity: The virtues of a nonclonal system of recognition. Cell 1997, 91, 295–298. [Google Scholar] [CrossRef]
  7. Schechter, L.M.; Vencato, M.; Jordan, K.L.; Schneider, S.E.; Schneider, D.J.; Collmer, A. Multiple approaches to a complete inventory of Pseudomonas syringae pv. tomato DC3000 type III secretion system effector proteins. Mol. Plant-Microbe Interact. 2006, 19, 1180–1192. [Google Scholar] [CrossRef]
  8. Yu, X.; Tang, J.; Wang, Q.; Ye, W.; Tao, K.; Duan, S.; Lu, C.; Yang, X.; Dong, S.; Zheng, X.; et al. The RxLR effector Avh241 from Phytophthora sojae requires plasma membrane localization to induce plant cell death. New Phytol. 2012, 196, 247–260. [Google Scholar] [CrossRef]
  9. Nejat, N.; Mantri, N. Plant immune system: Crosstalk between responses to biotic and abiotic stresses the missing link in understanding plant defence. Curr. Issues Mol. Biol. 2017, 23, 1–16. [Google Scholar] [CrossRef]
  10. Zhou, Y.; Huang, J.L.; Zhang, X.L.; Zhu, L.M.; Wang, X.F.; Guo, N.; Zhao, J.M.; Xing, H. Overexpression of chalcone isomerase (CHI) increases resistance against Phytophthora sojae in Soybean. J. Plant Biol. 2018, 61, 309–319. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Wang, W.; Chen, J.; Liu, J.; Xia, M.; Shen, F. Identification of miRNAs and their targets in cotton inoculated with Verticillium dahliae by high-throughput sequencing and degradome analysis. Int. J. Mol. Sci. 2015, 16, 14749–14768. [Google Scholar] [CrossRef]
  12. Song, C.; Xia, S.; Wang, X.; Zhang, J.; Qin, H.; Zhang, Y.; Feng, C.; Bie, S. Cloning and expression analysis of chalcone synthase and chalcone isomerase encoding genes in Gossypium hirsutum. Agric. Biotechnol. 2018, 7, 15–26. [Google Scholar]
  13. Zhang, X.M.; Qi, X.G.; Chen, Y.B.; Yao, H.Y.; Zheng, Y.J. Research advances of flavonoids in strawberry. North. Hortic. 2020, 1, 128–133. (In Chinese) [Google Scholar]
  14. Ralston, L.; Subramanian, S.; Matsuno, M.; Yu, O. Partial reconstruction of flavonoid and isoflavonoid biosynthesis in yeast using soybean type I and type II chalcone isomerases. Plant Physiol. 2005, 137, 1375–1388. [Google Scholar] [CrossRef]
  15. Kimura, Y.; Aoki, T.; Ayabe, S.-i. Chalcone isomerase isozymes with different substrate specificities towards 6′-Hydroxy- and 6′-Deoxychalcones in cultured cells of glycyrrhiza echinata, a leguminous plant pro-ducing 5-deoxyflavonoids. Plant Cell Physiol. 2001, 42, 1169–1173. [Google Scholar] [CrossRef] [PubMed]
  16. Ngaki, M.N.; Louie, G.V.; Philippe, R.N.; Manning, G.; Pojer, F.; Bowman, M.E.; Li, L.; Larsen, E.; Wurtele, E.S.; Noel, J.P. Evolution of the chalcone-isomerase fold from fatty-acid binding to stereospecific catalysis. Nature 2012, 485, 530–533. [Google Scholar] [CrossRef]
  17. Hartmann, T. From waste products to ecochemicals: Fifty years research of plant secondary metabolism. Phytochemistry 2007, 68, 2831–2846. [Google Scholar] [CrossRef]
  18. Kaltenbach, M.; Burke, J.R.; Dindo, M.; Pabis, A.; Munsberg, F.S.; Rabin, A.; Kamerlin, S.C.L.; Noel, J.P.; Tawfik, D.S. Evolution of chalcone isomerase from a noncatalytic ancestor. Nat. Chem. Biol. 2018, 14, 548–555. [Google Scholar] [CrossRef]
  19. Ban, Z.; Qin, H.; Mitchell, A.J.; Liu, B.; Zhang, F.; Weng, J.K.; Dixon, R.A.; Wang, G. Noncatalytic chalcone isomerase-fold proteins in Humulus lupulus are auxiliary components in prenylated flavonoid biosynthesis. Proc. Natl. Acad. Sci. USA 2018, 115, E5223–E5232. [Google Scholar] [CrossRef]
  20. Cheng, A.X.; Zhang, X.; Han, X.J.; Zhang, Y.Y.; Gao, S.; Liu, C.J.; Lou, H.X. Identification of chalcone isomerase in the basal land plants reveals an ancient evolution of enzymatic cyclization activity for synthesis of flavonoids. New Phytol. 2017, 217, 909–924. [Google Scholar] [CrossRef] [PubMed]
  21. Ohno, S. Evolution by Gene Duplication; Springer: New York, NY, USA, 1970. [Google Scholar]
  22. Flagel, L.E.; Wendel, J.F. Gene duplication and evolutionary novelty in plants. New Phytol. 2009, 183, 557–564. [Google Scholar] [CrossRef]
  23. Magadum, S.; Banerjee, U.; Murugan, P.; Gangapur, D.; Ravikesavan, R. Gene duplication as a major force in evolution. J. Genet. 2013, 92, 155–161. [Google Scholar] [CrossRef]
  24. Birchler, J.A.; Yang, H. The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation. Plant Cell 2022, 34, 2466–2474. [Google Scholar] [CrossRef] [PubMed]
  25. Panchy, N.; Lehti-Shiu, M.; Shiu, S.H. Evolution of gene duplication in plants. Plant Physiol. 2016, 171, 2294–2316. [Google Scholar] [CrossRef] [PubMed]
  26. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef] [PubMed]
  27. Zu, Q.L.; Qu, Y.Y.; Ni, Z.Y.; Zheng, K.; Chen, Q.; Chen, Q.J. The chalcone isomerase family in cotton: Whole-genome bioinformatic and expression analyses of the Gossypium barbadense L. response to Fusarium wilt infection. Genes 2019, 10, 1006. [Google Scholar] [CrossRef]
  28. Guo, W.; Jin, L.; Miao, Y.; He, X.; Hu, Q.; Guo, K.; Zhu, L.; Zhang, X. An ethylene response-related factor, GbERF1-like, from Gossypium barbadense improves resistance to Verticillium dahliae via activating lignin synthesis. Plant Mol. Biol. 2016, 91, 305–318. [Google Scholar] [CrossRef]
  29. Gu, Z.; Rifkin, S.A.; White, K.P.; Li, W.-H. Duplicate genes increase gene expression diversity within and between species. Nat. Genet. 2004, 36, 577–579. [Google Scholar] [CrossRef]
  30. Teichmann, S.A.; Babu, M.M. Gene regulatory network growth by duplication. Nat. Genet. 2004, 36, 492–496. [Google Scholar] [CrossRef]
  31. Voordeckers, K.; Pougach, K.; Verstrepen, K.J. How do regulatory networks evolve and expand throughout evolution? Curr. Opin. Biotechnol. 2015, 34, 180–188. [Google Scholar] [CrossRef]
  32. Chen, Z.J. Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu. Rev. Plant Biol. 2007, 58, 377–406. [Google Scholar] [CrossRef]
  33. Strygina, K.V.; Khlestkina, E.K. Myc-like transcriptional factors in wheat: Structural and functional organization of the subfamily I members. BMC Plant Biol. 2019, 19, 50. [Google Scholar] [CrossRef] [PubMed]
  34. Strygina, K.V.; Börner, A.; Khlestkina, E.K. Identification and characterization of regulatory network components for anthocyanin synthesis in barley aleurone. BMC Plant Biol. 2017, 17, 184. [Google Scholar] [CrossRef]
  35. Strygina, K.V.; Khlestkina, E.K. Structural and functional divergence of the Mpc1 genes in wheat and barley. BMC Evol. Biol. 2019, 19, 45. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, Z.J.; Birchler, J.A. Polyploid and Hybrid Genomics; John Wiley & Sons, Inc.: Oxford, UK, 2013. [Google Scholar]
  37. Wang, M.; Li, J.; Qi, Z.; Long, Y.; Pei, L.; Huang, X.; Grover, C.E.; Du, X.; Xia, C.; Wang, P.; et al. Genomic innovation and regulatory rewiring during evolution of the cotton genus Gossypium. Nat. Genet. 2022, 54, 1959–1971. [Google Scholar] [CrossRef]
  38. Lunn, J.E. Compartmentation in plant metabolism. J. Exp. Bot. 2006, 58, 35–47. [Google Scholar] [CrossRef]
  39. Tanz, S.K.; Castleden, I.; Small, I.D.; Millar, A.H. Fluorescent protein tagging as a tool to define the subcellular distribution of proteins in plants. Front. Plant Sci. 2013, 4, 214. [Google Scholar] [CrossRef]
  40. Rao, L.P.; Su, W.J.; Liu, Y.; Song, T.X.; Soviguidi, D.R.J.; Zhang, W.Y.; Yang, X.M. Cloning and expression characteristics of chalcone isomerase gene IbCHI from sweet potato. J. Plant Resour. Environ. 2020, 20, 1–10. (In Chinese) [Google Scholar]
  41. Wang, H.; Wang, W.; Zhan, J.; Huang, W.; Xu, H. Tissue-specific accumulation and subcellular localization of chalcone isomerase (CHI) in grapevine. Plant Cell Tissue Organ Cult. 2019, 137, 125–137. [Google Scholar] [CrossRef]
  42. Xing, H.R.; Liu, L.J.; Liu, G.Z. Advancement of protein subcellular localization in plants. Acta Agric. Boreali Sin. 2006, 21, 1–16. [Google Scholar]
  43. Kuzmin, E.; Taylor, J.S.; Boone, C. Retention of duplicated genes in evolution. Trends Genet. 2022, 38, 59–72. [Google Scholar] [CrossRef]
  44. Jiao, W.B.; Patel, V.; Klasen, J.; Liu, F.; Pecinkova, P.; Ferrand, M.; Gy, I.; Camilleri, C.; Effgen, S.; Koornneef, M.; et al. The evolutionary dynamics of genetic incompatibilities introduced by duplicated genes in Arabidopsis thaliana. Mol. Biol. Evol. 2021, 38, 1225–1240. [Google Scholar] [CrossRef]
  45. Innan, H.; Kondrashov, F. The evolution of gene duplications: Classifying and distinguishing between models. Nat. Rev. Genet. 2010, 11, 97–108. [Google Scholar] [CrossRef]
  46. Conrad, B.; Antonarakis, S.E. Gene duplication: A drive for phenotypic diversity and cause of human disease. Annu. Rev. Genom. Hum. Genet. 2007, 8, 17–35. [Google Scholar] [CrossRef] [PubMed]
  47. Albalat, R.; Cañestro, C. Evolution by gene loss. Nat. Rev. Genet. 2016, 17, 379–391. [Google Scholar] [CrossRef]
  48. Diss, G.; Ascencio, D.; DeLuna, A.; Landry, C.R. Molecular mechanisms of paralogous compensation and the robustness of cellular networks. J. Exp. Zool. B Mol. Dev. Evol. 2013, 322, 488–499. [Google Scholar] [CrossRef] [PubMed]
  49. El-Brolosy, M.A.; Stainier, D.Y.R. Genetic compensation: A phenomenon in search of mechanisms. PLoS Genet. 2017, 13, e1006780. [Google Scholar] [CrossRef]
  50. Hanada, K.; Sawada, Y.; Kuromori, T.; Klausnitzer, R.; Saito, K.; Toyoda, T.; Shinozaki, K.; Li, W.-H.; Hirai, M.Y. Functional compensation of primary and secondary metabolites by duplicate genes in Arabidopsis thaliana. Mol. Biol. Evol. 2011, 28, 377–382. [Google Scholar] [CrossRef]
  51. Kafri, R.; Bar-Even, A.; Pilpel, Y. Transcription control reprogramming in genetic backup circuits. Nat. Genet. 2005, 37, 295–299. [Google Scholar] [CrossRef]
  52. DeLuna, A.; Vetsigian, K.; Shoresh, N.; Hegreness, M.; Colón-González, M.; Chao, S.; Kishony, R. Exposing the fitness contribution of duplicated genes. Nat. Genet. 2008, 40, 676–681. [Google Scholar] [CrossRef] [PubMed]
  53. Freeling, M. Bias in plant gene content following different sorts of duplication: Tandem, whole-genome, segmental, or by transposition. Annu. Rev. Plant Biol. 2009, 60, 433–453. [Google Scholar] [CrossRef] [PubMed]
  54. Vavouri, T.; Semple, J.I.; Lehner, B. Widespread conservation of genetic redundancy during a billion years of eukaryotic evolution. Trends Genet. 2008, 24, 485–488. [Google Scholar] [CrossRef]
  55. Rodriguez-Leal, D.; Xu, C.; Kwon, C.T.; Soyars, C.; Demesa-Arevalo, E.; Man, J.; Liu, L.; Lemmon, Z.H.; Jones, D.S.; Van Eck, J.; et al. Evolution of buffering in a genetic circuit controlling plant stem cell proliferation. Nat. Genet. 2019, 51, 786–792. [Google Scholar] [CrossRef] [PubMed]
  56. Qiao, X.Y.; Ma, C.L.; Chen, L. Plant flavonoid biosynthesis pathway and regulation of its important genes. Nat. Prod. Res. Dev. 2009, 21, 354–360. [Google Scholar]
  57. Xing, W.; Jin, X.L. Recent advances of MYB transcription factors involved in the regulation of flavonoid biosynthesis. Mol. Plant Breed. 2015, 13, 689–696. [Google Scholar]
  58. Guo, X.W.; Huang, C.L.; Wu, Z.Y.; Zhang, X.H.; Luo, C.; Cheng, X. Molecular regulation of plant flavonoid biosynthesis pathhway. North. Hortic. 2011, 1, 204–207. (In Chinese) [Google Scholar]
  59. Pandey, A.; Agrawal, S.B. Ultraviolet-b radiation: A potent regulator of flavonoids biosynthesis, accumulation and functions in plants. Curr. Sci. 2020, 119, 176–185. [Google Scholar] [CrossRef]
  60. Bednar, R.A.; Hadcock, J.R. Purification and characterization of chalcone isomerase from soybeans. J. Biol. Chem. 1988, 263, 9582–9588. [Google Scholar] [CrossRef]
  61. Grotewold, E.; Peterson, T. Isolation and characterization of a maize gene encoding chalcone flavonone isomerase. Mol. Gen. Genet. 1994, 242, 1–8. [Google Scholar] [CrossRef]
  62. Shoeva, O.Y.; Khlestkina, E.K.; Berges, H.; Salina, E.A. The homoeologous genes encoding chalcone–flavanone isomerase in Triticum aestivum L.: Structural characterization and expression in different parts of wheat plant. Gene 2014, 538, 334–341. [Google Scholar] [CrossRef]
  63. McKhann, H.I.; Hirsch, A.M. Isolation of chalcone synthase and chalcone isomerase cDNAs from alfalfa (Medicago sativa L.): Highest transcript levels occur in young roots and root tips. Plant Mol. Biol. 1994, 24, 767–777. [Google Scholar] [CrossRef]
  64. Geng, L.; Duan, X.; Liang, C.; Shu, C.; Song, F.; Zhang, J. Mining tissue-specific contigs from peanut (Arachis hypogaea L.) for promoter cloning by deep transcriptome sequencing. Plant Cell Physiol. 2014, 55, 1793–1801. [Google Scholar] [CrossRef] [PubMed]
  65. Qin, J.C.; Zhu, L.; Gao, M.J.; Wu, X.; Pan, H.Y.; Zhang, Y.S.; Li, X. Cloning and functional characterization of a chalcone isomerase from Trigonella foenum-graecum L. Planta Med. 2010, 77, 765–770. [Google Scholar] [CrossRef]
  66. Forkmann, G.; Dangelmayr, B. Genetic control of chalcone isomerase activity in flowers of Dianthus caryophyllus. Biochem. Genet. 1980, 18, 519–527. [Google Scholar] [CrossRef] [PubMed]
  67. Cheng, H.; Li, L.; Cheng, S.; Cao, F.; Wang, Y.; Yuan, H. Molecular cloning and function assay of a chalcone isomerase gene (GbCHI) from Ginkgo biloba. Plant Cell Rep. 2010, 30, 49–62. [Google Scholar] [CrossRef]
  68. Kuhn, B.; Forkmann, G.; Seyffert, W. Genetic control of chalcone-flavanone isomerase activity in Callistephus chinensis. Planta 1978, 138, 199–203. [Google Scholar] [CrossRef]
  69. Takamura, T.; Tomihama, T.; Miyajima, I. Inheritance of yellow-flowered characteristic and yellow pigments in diploid cyclamen (Cyclamen persicum Mill.) cultivars. Sci. Hortic. 1995, 64, 55–63. [Google Scholar] [CrossRef]
  70. Itoh, Y.; Higeta, D.; Suzuki, A.; Yoshida, H.; Ozeki, Y. Excision of transposable elements from the chalcone isomerase and dihydroflavonol 4-reductase genes may contribute to the variegation of the yellow-flowered carnation (Dianthus caryophyllus). Plant Cell Physiol. 2002, 43, 578–585. [Google Scholar] [CrossRef]
  71. Nishihara, M.; Nakatsuka, T.; Yamamura, S. Flavonoid components and flower color change in transgenic tobacco plants by suppression of chalcone isomerase gene. FEBS Lett. 2005, 579, 6074–6078. [Google Scholar] [CrossRef]
  72. White, L.R.; Petrovitch, H.; Ross, G.W.; Masaki, K.; Hardman, J.; Nelson, J.; Davis, D.; Markesbery, W. Brain aging and midlife Tofu consumption. J. Am. Coll. Nutr. 2000, 19, 242–255. [Google Scholar] [CrossRef]
  73. Li, F.X.; Jin, Z.P.; Zhao, D.X.; Cheng, L.Q.; Fu, C.X.; Ma, F. Overexpression of the Saussurea medusa chalcone isomerase gene in S. involucrata hairy root cultures enhances their biosynthesis of apigenin. Phytochemistry 2006, 67, 553–560. [Google Scholar] [CrossRef]
  74. Li, F.; Jin, Z.; Qu, W.; Zhao, D.; Ma, F. Cloning of a cDNA encoding the Saussurea medusa chalcone isomerase and its expression in transgenic tobacco. Plant Physiol. Biochem. 2006, 44, 455–461. [Google Scholar] [CrossRef]
  75. Lukaszewicz, M.; Matysiak-Kata, I.; Skala, J.; Fecka, I.; Cisowski, W.; Szopa, J. Antioxidant capacity manipulation in transgenic potato tuber by changes in phenolic compounds content. J. Agric. Food Chem. 2004, 52, 1526–1533. [Google Scholar] [CrossRef]
  76. Shi, L.; Feng, X.E.; Cui, J.R.; Fang, L.H.; Du, G.H.; Li, Q.S. Synthesis and biological activity of flavanone derivatives. Bioorg. Med. Chem. Lett. 2010, 20, 5466–5468. [Google Scholar] [CrossRef] [PubMed]
  77. Hughes, A.L. The evolution of functionally novel proteins after gene duplication. Proc. R. Soc. Lond. B Biol. Sci. 1994, 256, 119–124. [Google Scholar] [CrossRef]
  78. Force, A.; Lynch, M.; Pickett, F.B.; Amores, A.; Yan, Y.L.; Postlethwait, J. Preservation of duplicate genes by complementary, degenerative mutations. Genetics 1999, 151, 1531–1545. [Google Scholar] [CrossRef]
  79. Rastogi, S.; Liberles, D.A. Subfunctionalization of duplicated genes as a transition state to neofunctionalization. BMC Evol. Biol. 2005, 5, 28. [Google Scholar] [CrossRef] [PubMed]
  80. Hallin, J.; Landry, C.R. Regulation plays a multifaceted role in the retention of gene duplicates. PLoS Biol. 2019, 17, e3000519. [Google Scholar] [CrossRef] [PubMed]
  81. Ren, L.L.; Liu, Y.J.; Liu, H.J.; Qian, T.T.; Qi, L.W.; Wang, X.R.; Zeng, Q.Y. Subcellular relocalization and positive selection play key roles in the retention of duplicate genes of populus class III peroxidase family. Plant Cell 2014, 26, 2404–2419. [Google Scholar] [CrossRef] [PubMed]
  82. Qiu, Y.; Tay, Y.V.; Ruan, Y.; Adams, K.L. Divergence of duplicated genes by repeated partitioning of splice forms and subcellular localization. New Phytol. 2019, 225, 1011–1022. [Google Scholar] [CrossRef]
  83. von der Dunk, S.H.A.; Snel, B. Recurrent sequence evolution after independent gene duplication. BMC Evol. Biol. 2020, 20, 98. [Google Scholar] [CrossRef]
  84. He, X.; Zhang, J. Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics 2005, 169, 1157–1164. [Google Scholar] [CrossRef]
  85. Gibson, T.A.; Goldberg, D.S. Questioning the ubiquity of neofunctionalization. PLoS Comput. Biol. 2009, 5, e1000252. [Google Scholar] [CrossRef] [PubMed]
  86. Diss, G.; Gagnon-Arsenault, I.; Dion-Coté, A.M.; Vignaud, H.; Ascencio, D.I.; Berger, C.M.; Landry, C.R. Gene duplication can impart fragility, not robustness, in the yeast protein interaction network. Science 2017, 355, 630–634. [Google Scholar] [CrossRef]
  87. Bergthorsson, U.; Andersson, D.I.; Roth, J.R. Ohno’s dilemma: Evolution of new genes under continuous selection. Proc. Natl. Acad. Sci. USA 2007, 104, 17004–17009. [Google Scholar] [CrossRef]
  88. Bari, R.; Jones, J.D.G. Role of plant hormones in plant defence responses. Plant Mol. Biol. 2008, 69, 473–488. [Google Scholar] [CrossRef] [PubMed]
  89. Wasternack, C. Jasmonates: An update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann. Bot. 2007, 100, 681–697. [Google Scholar] [CrossRef]
  90. Grant, M.; Lamb, C. Systemic immunity. Curr. Opin. Plant Biol. 2006, 9, 414–420. [Google Scholar] [CrossRef]
  91. Truman, W.; Bennett, M.H.; Kubigsteltig, I.; Turnbull, C.; Grant, M. Arabidopsis systemic immunity uses conserved defense signaling pathways and is mediated by jasmonates. Proc. Natl. Acad. Sci. USA 2007, 104, 1075–1080. [Google Scholar] [CrossRef]
  92. Zhang, Y.; Wang, X.; Rong, W.; Yang, J.; Li, Z.; Wu, L.; Zhang, G.; Ma, Z. Histochemical analyses reveal that stronger intrinsic defenses in Gossypium barbadense than in G. hirsutum are associated with resistance to Verticillium dahliae. Mol. Plant-Microbe Interact. 2017, 30, 984–996. [Google Scholar] [CrossRef]
  93. Park, S.W.; Kaimoyo, E.; Kumar, D.; Mosher, S.; Klessig, D.F. Methyl salicylate is a critical mobile signal for plant systemic acquired resistance. Science 2007, 318, 113–116. [Google Scholar] [CrossRef] [PubMed]
  94. Mittler, R.; Zandalinas, S.I.; Fichman, Y.; Van Breusegem, F. Reactive oxygen species signalling in plant stress responses. Nat. Rev. Mol. Cell Biol. 2022, 23, 663–679. [Google Scholar] [CrossRef] [PubMed]
  95. Wei, C.; Wen, C.; Zhang, Y.; Du, H.; Zhong, R.; Guan, Z.; Wang, M.; Qin, Y.; Wang, F.; Song, L.; et al. The FomYjeF protein influences the sporulation and virulence of Fusarium oxysporum f. sp. momordicae. Int. J. Mol. Sci. 2023, 24, 7260. [Google Scholar] [CrossRef]
  96. Chae, H.B.; Kim, M.G.; Kang, C.H.; Park, J.H.; Lee, E.S.; Lee, S.U.; Chi, Y.H.; Paeng, S.K.; Bae, S.B.; Wi, S.D.; et al. Redox sensor QSOX1 regulates plant immunity by targeting GSNOR to modulate ROS generation. Mol. Plant 2021, 14, 1312–1327. [Google Scholar] [CrossRef]
  97. Wang, N.; Fan, X.; He, M.; Hu, Z.; Tang, C.; Zhang, S.; Lin, D.; Gan, P.; Wang, J.; Huang, X.; et al. Transcriptional repression of TaNOX10 by TaWRKY19 compromises ROS generation and enhances wheat susceptibility to stripe rust. Plant Cell 2022, 34, 1784–1803. [Google Scholar] [CrossRef]
  98. Zhang, L.; Wu, Y.; Yu, Y.; Zhang, Y.; Wei, F.; Zhu, Q.H.; Zhou, J.; Zhao, L.; Zhang, Y.; Feng, Z.; et al. Acetylation of GhCaM7 enhances cotton resistance to Verticillium dahliae. Plant J. 2023, 114, 1405–1424. [Google Scholar] [CrossRef] [PubMed]
  99. De Smet, R.; Sabaghian, E.; Li, Z.; Saeys, Y.; Van de Peer, Y. Coordinated functional divergence of genes after genome duplication in Arabidopsis thaliana. Plant Cell 2017, 29, 2786–2800. [Google Scholar] [CrossRef]
  100. Kwon, C.T.; Tang, L.; Wang, X.; Gentile, I.; Hendelman, A.; Robitaille, G.; Van Eck, J.; Xu, C.; Lippman, Z.B. Dynamic evolution of small signalling peptide compensation in plant stem cell control. Nat. Plants 2022, 8, 346–355. [Google Scholar] [CrossRef]
  101. Wei, Q.; Liu, J.; Guo, F.; Wang, Z.; Zhang, X.; Yuan, L.; Ali, K.; Qiang, F.; Wen, Y.; Li, W.; et al. Kinase regulators evolved into two families by gain and loss of ability to bind plant steroid receptors. Plant Physiol. 2023, 191, 1167–1185. [Google Scholar] [CrossRef]
  102. Qiao, Y.; Shi, J.; Zhai, Y.; Hou, Y.; Ma, W. Phytophthora effector targets a novel component of small RNA pathway in plants to promote infection. Proc. Natl. Acad. Sci. USA 2015, 112, 5850–5855. [Google Scholar] [CrossRef]
  103. Gao, X.; Britt, R.C.; Shan, L.; He, P. Agrobacterium-mediated virus-induced gene silencing assay in cotton. J. Vis. Exp. 2011, 54, e2938. [Google Scholar] [CrossRef]
  104. Nazeer, W.; Tipu, A.L.; Ahmad, S.; Mahmood, K.; Mahmood, A.; Zhou, B. Evaluation of cotton leaf curl virus resistance in BC1, BC2, and BC3 progenies from an interspecific cross between Gossypium arboreum and Gossypium hirsutum. PLoS ONE 2014, 9, e111861. [Google Scholar] [CrossRef] [PubMed]
  105. Zu, Q.; Qu, Y.; Su, X.; Zheng, K.; Chen, Q.; Deng, X.; Gao, W.; Zhao, J.; Chen, Q. GbC4H regulates the metabolic flow of flavonoids and inhibits the occurrence of Fusarium wilt in sea island cotton. Plant Growth Regul. 2023, 101, 87–97. [Google Scholar] [CrossRef]
  106. Du, L.L.; Wu, P.H.; Qu, Y.Y.; Li, X.Y.; Chen, Q.J. Analysis and evluation of Fusarium wilt resistance for diferent populations in Sea island cotton (Gossypium barbadense). Mol. Plant Breed. 2015, 13, 1517–1525. [Google Scholar]
  107. Wang, H.; Wang, P.; Qin, J.; Zhou, H. Research profile of the extraction, purification and determination of flavonoids in plants. J. Jilin Univ. Chem. Technol. 2018, 35, 44–47. (In Chinese) [Google Scholar]
Figure 1. Results of subcellular localization analysis of the GbCHI gene from cotton in epidermal cells of Nicotiana benthamiana L. Bars, 25 µm.
Figure 1. Results of subcellular localization analysis of the GbCHI gene from cotton in epidermal cells of Nicotiana benthamiana L. Bars, 25 µm.
Ijms 24 14775 g001
Figure 2. Resistance phenotype of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced cotton plants to Fusarium wilt. (a) When Cloroplastos alterados 1 (CLA1) was used as a positive control, the positive control (pTRV2::CLA1) cotton plants showed an albino phenotype. (b) Resistance phenotypes of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09). (c) The GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene silencing efficiency test. In terms of statistical significance, “**” indicates p < 0.01. (d) The disease index of the phenotypes. a, b, and c above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Figure 2. Resistance phenotype of GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced cotton plants to Fusarium wilt. (a) When Cloroplastos alterados 1 (CLA1) was used as a positive control, the positive control (pTRV2::CLA1) cotton plants showed an albino phenotype. (b) Resistance phenotypes of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09). (c) The GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene silencing efficiency test. In terms of statistical significance, “**” indicates p < 0.01. (d) The disease index of the phenotypes. a, b, and c above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Ijms 24 14775 g002
Figure 3. Fusarium wilt recovery assay in silenced plants. The stems 3 mm above the cotyledon node of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) after infection with Fusarium wilt were taken. (a) The stems 3 mm above the cotyledon node of the control plants and the GbCHI01 gene-silenced plants. (b) The stems 3 mm above the cotyledon node of the control plants and the GbCHI05 gene-silenced plants after infection with Fusarium wilt were taken. (c) The stems 3 mm above the cotyledon node of the control plants and the GbCHI06 after infection with Fusarium wilt were taken. (d) The stems 3 mm above the cotyledon node of the control plants and the GbCHI09 gene-silenced plants after infection with Fusarium wilt were taken.
Figure 3. Fusarium wilt recovery assay in silenced plants. The stems 3 mm above the cotyledon node of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) after infection with Fusarium wilt were taken. (a) The stems 3 mm above the cotyledon node of the control plants and the GbCHI01 gene-silenced plants. (b) The stems 3 mm above the cotyledon node of the control plants and the GbCHI05 gene-silenced plants after infection with Fusarium wilt were taken. (c) The stems 3 mm above the cotyledon node of the control plants and the GbCHI06 after infection with Fusarium wilt were taken. (d) The stems 3 mm above the cotyledon node of the control plants and the GbCHI09 gene-silenced plants after infection with Fusarium wilt were taken.
Ijms 24 14775 g003
Figure 4. (a) Phenotypic identification of GbCHI05 gene overexpression in Arabidopsis thaliana for resistance to Fusarium wilt. WT: wild type. chi05: a GbCHI05 gene mutant. OE-1, OE-2 and OE-3: three GbCHI05 transgenic lines. (b) Fusarium wilt recovery assay in A. thaliana GbCHI05-overexpressing plants. After infection with Fusarium wilt, control (WT), mutant (chi05) and GbCHI05-overexpressing A. thaliana plants were grown, and the 3 mm stem segment of the hypocotyledonary axis was examined. (c) The disease index of the phenotypes. a, b, and c above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Figure 4. (a) Phenotypic identification of GbCHI05 gene overexpression in Arabidopsis thaliana for resistance to Fusarium wilt. WT: wild type. chi05: a GbCHI05 gene mutant. OE-1, OE-2 and OE-3: three GbCHI05 transgenic lines. (b) Fusarium wilt recovery assay in A. thaliana GbCHI05-overexpressing plants. After infection with Fusarium wilt, control (WT), mutant (chi05) and GbCHI05-overexpressing A. thaliana plants were grown, and the 3 mm stem segment of the hypocotyledonary axis was examined. (c) The disease index of the phenotypes. a, b, and c above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Ijms 24 14775 g004
Figure 5. Phenotypic identification of GbCHI05 gene overexpression in cotton against Fusarium wilt and Fusarium wilt recovery assay. (a) Resistance phenotype analysis under GbCHI05 gene overexpression in cotton. (b) Determination of GbCHI05 expression level in cotton plants overexpressing this gene. (c) Fusarium wilt recovery assay. WT: 06-146 (Gossypium barbadense L.). OE-1: transgenic lines of the GbCHI05 gene in T3 seedlings. The receptor of the transgenic material is 06-146. In terms of statistical significance, “**” indicates p < 0.01. (d) The disease index of the phenotypes. In terms of statistical significance, “*” indicates p < 0.05.
Figure 5. Phenotypic identification of GbCHI05 gene overexpression in cotton against Fusarium wilt and Fusarium wilt recovery assay. (a) Resistance phenotype analysis under GbCHI05 gene overexpression in cotton. (b) Determination of GbCHI05 expression level in cotton plants overexpressing this gene. (c) Fusarium wilt recovery assay. WT: 06-146 (Gossypium barbadense L.). OE-1: transgenic lines of the GbCHI05 gene in T3 seedlings. The receptor of the transgenic material is 06-146. In terms of statistical significance, “**” indicates p < 0.01. (d) The disease index of the phenotypes. In terms of statistical significance, “*” indicates p < 0.05.
Ijms 24 14775 g005
Figure 6. Detection of genes involved in the flavonoid pathway and genes related to Fusarium wilt in control plants (pTRV2::00) and GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09). Samples were taken at 0, 2, 4, 8, 12, 24, 48 and 72 h; (a) GbCHI01 gene-silenced plants; (b) GbCHI05 gene-silenced plants; (c) GbCHI06 gene-silenced plants; (d) GbCHI09 gene-silenced plants. In terms of statistical significance, “*” is used to indicate p < 0.05, and “**” indicates p < 0.01.
Figure 6. Detection of genes involved in the flavonoid pathway and genes related to Fusarium wilt in control plants (pTRV2::00) and GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09). Samples were taken at 0, 2, 4, 8, 12, 24, 48 and 72 h; (a) GbCHI01 gene-silenced plants; (b) GbCHI05 gene-silenced plants; (c) GbCHI06 gene-silenced plants; (d) GbCHI09 gene-silenced plants. In terms of statistical significance, “*” is used to indicate p < 0.05, and “**” indicates p < 0.01.
Ijms 24 14775 g006
Figure 7. Establishment of a rutin standard curve for measuring flavonoid content in cotton. Flavonoid content of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) before and after infection with the causal agent of Fusarium wilt. a, b, c, d, e, f and g above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Figure 7. Establishment of a rutin standard curve for measuring flavonoid content in cotton. Flavonoid content of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) before and after infection with the causal agent of Fusarium wilt. a, b, c, d, e, f and g above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Ijms 24 14775 g007
Figure 8. Endogenous hormone content (MeJA and SA) after Fov infection. pTRV2::00: control plants. pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09: GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants. Samples were taken at 0, 2, 4, 8, 12, 24, 48 and 72 h; “*” indicates p < 0.05, and “**” indicates p < 0.01.
Figure 8. Endogenous hormone content (MeJA and SA) after Fov infection. pTRV2::00: control plants. pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09: GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants. Samples were taken at 0, 2, 4, 8, 12, 24, 48 and 72 h; “*” indicates p < 0.05, and “**” indicates p < 0.01.
Ijms 24 14775 g008
Figure 9. ROS contents after Fov infection. ROS contents of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) before and after infection with the causal agent of Fusarium wilt. a, b, c, d, e, f, g, h and i above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Figure 9. ROS contents after Fov infection. ROS contents of the control plants (pTRV2::00) and the GbCHI01, GbCHI05, GbCHI06 and GbCHI09 gene-silenced plants (pTRV2::GbCHI01, pTRV2::GbCHI05, pTRV2::GbCHI06 and pTRV2::GbCHI09) before and after infection with the causal agent of Fusarium wilt. a, b, c, d, e, f, g, h and i above the columns indicate significant differences (p < 0.05) according to one-way ANOVA.
Ijms 24 14775 g009
Figure 10. The “one major, three secondary” mode in the CHI gene family.
Figure 10. The “one major, three secondary” mode in the CHI gene family.
Ijms 24 14775 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zu, Q.; Deng, X.; Qu, Y.; Chen, X.; Cai, Y.; Wang, C.; Li, Y.; Chen, Q.; Zheng, K.; Liu, X.; et al. Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L. Int. J. Mol. Sci. 2023, 24, 14775. https://doi.org/10.3390/ijms241914775

AMA Style

Zu Q, Deng X, Qu Y, Chen X, Cai Y, Wang C, Li Y, Chen Q, Zheng K, Liu X, et al. Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L. International Journal of Molecular Sciences. 2023; 24(19):14775. https://doi.org/10.3390/ijms241914775

Chicago/Turabian Style

Zu, Qianli, Xiaojuan Deng, Yanying Qu, Xunji Chen, Yongsheng Cai, Caoyue Wang, Ying Li, Qin Chen, Kai Zheng, Xiaodong Liu, and et al. 2023. "Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L." International Journal of Molecular Sciences 24, no. 19: 14775. https://doi.org/10.3390/ijms241914775

APA Style

Zu, Q., Deng, X., Qu, Y., Chen, X., Cai, Y., Wang, C., Li, Y., Chen, Q., Zheng, K., Liu, X., & Chen, Q. (2023). Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L. International Journal of Molecular Sciences, 24(19), 14775. https://doi.org/10.3390/ijms241914775

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop