Divergence of duplicated genes by repeated partitioning of splice forms and subcellular localization

Physiological function and regulation of ascorbate peroxidase isoforms

Ascorbate peroxidase (APX) reduces H2O2 to H2O by utilizing ascorbate as a specific electron donor and constitutes the ascorbate-glutathione cycle in organelles of plants including chloroplasts, cytosol, mitochondria, and peroxisomes. It has been almost 40 years since APX was discovered as an important plant-specific H2O2-scavenging enzyme, during which time many research groups have conducted molecular physiological analyses. It is now clear that APX isoforms function not only just as antioxidant enzymes but also as important factors in intracellular redox regulation through the metabolism of reactive oxygen species. The function of APX isoforms is regulated at multiple steps, from the transcriptional level to post-translational modifications of enzymes, thereby allowing them to respond flexibly to ever-changing environmental factors and physiological phenomena such as cell growth and signal transduction. In this review, we summarize the physiological functions and regulation mechanisms of expression of each APX isoform.

Ascorbate peroxidase in fruits and modulation of its activity by reactive species

Ascorbate peroxidase (APX) is one of the enzymes of the ascorbate-glutathione cycle and is the key enzyme that breaks down H2O2 with the aid of ascorbate as an electron source. APX is present in all photosynthetic eukaryotes from algae to higher plants and, at the cellular level, it is localized in all subcellular compartments where H2O2 is generated, including the apoplast, cytosol, plastids, mitochondria, and peroxisomes, either in soluble form or attached to the organelle membranes. APX activity can be modulated by various post-translational modifications including tyrosine nitration, S-nitrosation, persulfidation, and S-sulfenylation. This allows the connection of H2O2 metabolism with other relevant signaling molecules such as NO and H2S, thus building a complex coordination system. In both climacteric and non-climacteric fruits, APX plays a key role during the ripening process and during post-harvest, since it participates in the regulation of both H2O2 and ascorbate levels affecting fruit quality. Currently, the exogenous application of molecules such as NO, H2S, H2O2, and, more recently, melatonin is seen as a new alternative to maintain and extend the shelf life and quality of fruits because they can modulate APX activity as well as other antioxidant systems. Therefore, these molecules are being considered as new biotechnological tools to improve crop quality in the horticultural industry.

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.

Kinase regulators evolved into two families by gain and loss of ability to bind plant steroid receptors

All biological functions evolve by fixing beneficial mutations and removing deleterious ones. Therefore, continuously fixing and removing the same essential function to separately diverge monophyletic gene families sounds improbable. Yet, here we report that brassinosteroid insensitive1 kinase inhibitor1 (BKI1)/membrane-associated kinase regulators (MAKRs) regulating a diverse function evolved into BKI1 and MAKR families from a common ancestor by respectively enhancing and losing ability to bind brassinosteroid receptor brassinosteroid insensitive1 (BRI1). The BKI1 family includes BKI1, MAKR1/BKI1-like (BKL) 1, and BKL2, while the MAKR family contains MAKR2-6. Seedless plants contain only BKL2. In seed plants, MAKR1/BKL1 and MAKR3, duplicates of BKL2, gained and lost the ability to bind BRI1, respectively. In angiosperms, BKL2 lost the ability to bind BRI1 to generate MAKR2, while BKI1 and MAKR6 were duplicates of MAKR1/BKL1 and MAKR3, respectively. In dicots, MAKR4 and MAKR5 were duplicates of MAKR3 and MAKR2, respectively. Importantly, BKI1 localized in the plasma membrane, but BKL2 localized to the nuclei while MAKR1/BKL1 localized throughout the whole cell. Importantly, BKI1 strongly and MAKR1/BKL1 weakly inhibited plant growth, but BKL2 and the MAKR family did not inhibit plant growth. Functional study of the chimeras of their N- and C-termini showed that only the BKI1 family was partially reconstructable, supporting stepwise evolution by a seesaw mechanism between their C- and N-termini to alternately gain an ability to bind and inhibit BRI1, respectively. Nevertheless, the C-terminal BRI1-interacting motif best defines the divergence of BKI1/MAKRs. Therefore, BKI1 and MAKR families evolved by gradually gaining and losing the same function, respectively, extremizing divergent evolution and adding insights into gene (BKI1/MAKR) duplication and divergence.

Ascorbate-Glutathione Cycle Genes Families in Euphorbiaceae: Characterization and Evolutionary Analysis

Ascorbate peroxidase (APX), Monodehydroascorbate Reductase (MDAR), Dehydroascorbate Reductase (DHAR) and Glutathione Reductase (GR) enzymes participate in the ascorbate-glutathione cycle, which exerts a central role in the antioxidant metabolism in plants. Despite the importance of this antioxidant system in different signal transduction networks related to development and response to environmental stresses, the pathway has not yet been comprehensively characterized in many crop plants. Among different eudicotyledons, the Euphorbiaceae family is particularly diverse with some species highly tolerant to drought. Here the APX, MDAR, DHAR, and GR genes in Ricinus communis, Jatropha curcas, Manihot esculenta, and Hevea brasiliensis were identified and characterized. The comprehensive phylogenetic and genomic analyses allowed the classification of the genes into different classes, equivalent to cytosolic, peroxisomal, chloroplastic, and mitochondrial enzymes, and revealed the duplication events that contribute to the expansion of these families within plant genomes. Due to the high drought stress tolerance of Ricinus communis, the expression patterns of ascorbate-glutathione cycle genes in response to drought were also analyzed in leaves and roots, indicating a differential expression during the stress. Altogether, these data contributed to the characterization of the expression pattern and evolutionary analysis of these genes, filling the gap in the proposed functions of core components of the antioxidant mechanism during stress response in an economically relevant group of plants.

Genome-wide, evolutionary, and functional analyses of ascorbate peroxidase (APX) family in Poaceae species

Ascorbate peroxidases (APXs) are heme peroxidases involved in the control of hydrogen peroxide levels and signal transduction pathways related to development and stress responses. Here, a total of 238 APX, 30 APX-related (APX-R), and 34 APX-like (APX-L) genes were identified from 24 species from the Poaceae family. Phylogenetic analysis of APX indicated five distinct clades, equivalent to cytosolic (cAPX), peroxisomal (pAPX), mitochondrial (mitAPX), stromal (sAPX), and thylakoidal (tAPX) isoforms. Duplication events contributed to the expansion of this family and the divergence times. Different from other APX isoforms, the emergence of Poaceae mitAPXs occurred independently after eudicot and monocot divergence. Our results showed that the constitutive silencing of mitAPX genes is not viable in rice plants, suggesting that these isoforms are essential for rice regeneration or development. We also obtained rice plants silenced individually to sAPX isoforms, demonstrating that, different to plants double silenced to both sAPX and tAPX or single silenced to tAPX previously obtained, these plants do not show changes in the total APX activity and hydrogen peroxide content in the shoot. Among rice plants silenced to different isoforms, plants silenced to cAPX showed a higher decrease in total APX activity and an increase in hydrogen peroxide levels. These results suggest that the cAPXs are the main isoforms responsible for regulating hydrogen peroxide levels in the cell, whereas in the chloroplast, this role is provided mainly by the tAPX isoform. In addition to broadening our understanding of the core components of the antioxidant defense in Poaceae species, the present study also provides a platform for their functional characterization.

Chloroplastic ascorbate peroxidases targeted to stroma or thylakoid membrane: The chicken or egg dilemma

Ascorbate peroxidases (APXs) are heme peroxidases that remove hydrogen peroxide in different subcellular compartments with concomitant ascorbate cycling. Here, we analysed and discussed phylogenetic and molecular features of the APX family. Ancient APX originated as a soluble stromal enzyme, and early during plant evolution, acquired both chloroplast-targeting and mitochondrion-targeting sequences and an alternative splicing mechanism whereby it could be expressed as a soluble or thylakoid membrane-bound enzyme. Later, independent duplication and neofunctionalization events in some angiosperm groups resulted in individual genes encoding stromal, thylakoidal and mitochondrial isoforms. These data reaffirm the complexity of plant antioxidant defenses that allow diverse plant species to acquire new means to adapt to changing environmental conditions.

Global Patterns of subgenome evolution in organelle-targeted genes of six allotetraploid angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary “mismatches” between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Expression of Two Rye CENH3 Variants and Their Loading into Centromeres

Gene duplication and the preservation of both copies during evolution is an intriguing evolutionary phenomenon. Their preservation is related to the function they perform. The central component of centromere specification and function is the centromere-specific histone H3 (CENH3). Some cereal species (maize, rice) have one copy of the gene encoding this protein, while some (wheat, barley, rye) have two. Therefore, they represent a good model for a comparative study of the functional activity of the duplicated CENH3 genes and their protein products. We determined the organization of the CENH3 locus in rye (Secale cereale L.) and identified the functional motifs in the vicinity of the CENH3 genes. We compared the expression of these genes at different stages of plant development and the loading of their products, the CENH3 proteins, into nucleosomes during mitosis and meiosis. Using extended chromatin fibers, we revealed patterns of loading CENH3 proteinsinto polynucleosomal domains in centromeric chromatin. Our results indicate no sign of neofunctionalization, subfunctionalization or specialization in the gene copies. The influence of negative selection on the coding part of the genes led them to preserve their conserved function. The advantage of having two functional genes appears as the gene-dosage effect.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

CropPAL for discovering divergence in protein subcellular location in crops to support strategies for molecular crop breeding

Agriculture faces increasing demand for yield, higher plant-derived protein content and diversity while facing pressure to achieve sustainability. Although the genomes of many of the important crops have been sequenced, the subcellular locations of most of the encoded proteins remain unknown or are only predicted. Protein subcellular location is crucial in determining protein function and accumulation patterns in plants, and is critical for targeted improvements in yield and resilience. Integrating location data from over 800 studies for 12 major crop species into the cropPAL2020 data collection showed that while >80% of proteins in most species are not localised by experimental data, combining species data or integrating predictions can help bridge gaps at similar accuracy. The collation and integration of over 61 505 experimental localisations and more than 6 million predictions showed that the relative sizes of the protein catalogues located in different subcellular compartments are comparable between crops and Arabidopsis. A comprehensive cross-species comparison showed that between 50% and 80% of the subcellulomes are conserved across species and that conservation only depends to some degree on the phylogenetic relationship of the species. Protein subcellular locations in major biosynthesis pathways are more often conserved than in metabolic pathways. Underlying this conservation is a clear potential for subcellular diversity in protein location between species by means of gene duplication and alternative splicing. Our cropPAL data set and search platform (https://crop-pal.org) provide a comprehensive subcellular proteomics resource to drive compartmentation-based approaches for improving yield, protein composition and resilience in future crop varieties.

Generating and characterizing single- and multigene mutants of the Rubisco small subunit family in Arabidopsis

The primary CO2-fixing enzyme Rubisco limits the productivity of plants. The small subunit of Rubisco (SSU) can influence overall Rubisco levels and catalytic efficiency, and is now receiving increasing attention as a potential engineering target to improve the performance of Rubisco. However, SSUs are encoded by a family of nuclear rbcS genes in plants, which makes them challenging to engineer and study. Here we have used CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9] and T-DNA insertion lines to generate a suite of single and multiple gene knockout mutants for the four members of the rbcS family in Arabidopsis, including two novel mutants 2b3b and 1a2b3b. 1a2b3b contained very low levels of Rubisco (~3% relative to the wild-type) and is the first example of a mutant with a homogenous Rubisco pool consisting of a single SSU isoform (1B). Growth under near-outdoor levels of light demonstrated Rubisco-limited growth phenotypes for several SSU mutants and the importance of the 1A and 3B isoforms. We also identified 1a1b as a likely lethal mutation, suggesting a key contributory role for the least expressed 1B isoform during early development. The successful use of CRISPR/Cas here suggests that this is a viable approach for exploring the functional roles of SSU isoforms in plants.

Enzyme evolution in natural products biosynthesis: target- or diversity-oriented?

Natural product biosynthesis (NPB) is the Panda's Thumb of evolutionary biochemistry. Arm races between organisms, and ever-changing environments, result in relentless innovation. This review focusses on enzyme evolution in NPB. First, we review cases of de novo emergence, whereby a completely new enzymatic activity arose in a ligand-binding protein, or a new enzyme emerged including a completely new scaffold. Second, we briefly review the current models for enzyme evolution, and how they explain the inherent promiscuity of NPB enzymes and their tendency to produce multiple related products. We thus suggest that NPB enzymes a priori evolved to generate a specific product; they are, however, trapped in a multifunctional, generalist evolutionary state and thereby produce a diversity of products.

Regulation plays a multifaceted role in the retention of gene duplicates

A gene duplication can lead to all sorts of problems in a cell. However, it can also lead to all sorts of benefits. Beneficial or not, the gene duplicates might be kept in the genome because of several different reasons. For instance, if natural selection works towards optimizing one function of a gene at the expense of another, then gene duplication could resolve this conflict by separating the functions in two genes. Here, we outline evolutionary incentives to keep a duplicated gene in the genome, focusing on divergence in expression and trade-off resolution as featured in a new and exciting paper published in this edition of PLOS Biology.