Functional and structural profiles of GST gene family from three Populus species reveal the sequence-function decoupling of orthologous genes.

A common assumption in comparative genomics is that orthologous genes are functionally more similar than paralogous genes. However, the validity of this assumption needs to be assessed using robust experimental data. We conducted tissue-specific gene expression and protein function analyses of orthologous groups within the glutathione S-transferase (GST) gene family in three closely related Populus species: Populus trichocarpa, Populus euphratica and Populus yatungensis. This study identified 21 GST orthologous groups in the three Populus species. Although the sequences of the GST orthologous groups were highly conserved, the divergence in enzymatic functions was prevalent. Through site-directed mutagenesis of orthologous proteins, this study revealed that nonsynonymous substitutions at key amino acid sites played an important role in the divergence of enzymatic functions. In particular, a single amino acid mutation (Arg39→Trp39) contributed to P. euphratica PeGSTU30 possessing high enzymatic activity via increasing the hydrophobicity of the active cavity. This study provided experimental evidence showing that orthologues belonging to the gene family have functional divergences. The nonsynonymous substitutions at a few amino acid sites resulted in functional divergence of the orthologous genes. Our findings provide new insights into the evolution of orthologous genes in closely related species.

Evolution and Function of the Populus SABATH Family Reveal That a Single Amino Acid Change Results in a Substrate Switch.

Evolutionary mechanisms of substrate specificities of enzyme families remain poorly understood. Plant SABATH methyltransferases catalyze methylation of the carboxyl group of various low molecular weight metabolites. Investigation of the functional diversification of the SABATH family in plants could shed light on the evolution of substrate specificities in this enzyme family. Previous studies identified 28 SABATH genes from the Populus trichocarpa genome. In this study, we re-annotated the Populus SABATH gene family, and performed molecular evolution, gene expression and biochemical analyses of this large gene family. Twenty-eight Populus SABATH genes were divided into three classes with distinct divergences in their gene structure, expression responses to abiotic stressors and enzymatic properties of encoded proteins. Populus class I SABATH proteins converted IAA to methyl-IAA, class II SABATH proteins converted benzoic acid (BA) and salicylic acid (SA) to methyl-BA and methyl-SA, while class III SABATH proteins converted farnesoic acid (FA) to methyl-FA. For Populus class II SABATH proteins, both forward and reverse mutagenesis studies showed that a single amino acid switch between PtSABATH4 and PtSABATH24 resulted in substrate switch. Our findings provide new insights into the evolution of substrate specificities of enzyme families.

Divergence in Enzymatic Activities in the Soybean GST Supergene Family Provides New Insight into the Evolutionary Dynamics of Whole-Genome Duplicates.

Whole-genome duplication (WGD), or polyploidy, is a major force in plant genome evolution. A duplicate of all genes is present in the genome immediately following a WGD event. However, the evolutionary mechanisms responsible for the loss of, or retention and subsequent functional divergence of polyploidy-derived duplicates remain largely unknown. In this study we reconstructed the evolutionary history of the glutathione S-transferase (GST) gene family from the soybean genome, and identified 72 GST duplicated gene pairs formed by a recent Glycine-specific WGD event occurring approximately 13 Ma. We found that 72% of duplicated GST gene pairs experienced gene losses or pseudogenization, whereas 28% of GST gene pairs have been retained in the soybean genome. The GST pseudogenes were under relaxed selective constraints, whereas functional GSTs were subject to strong purifying selection. Plant GST genes play important roles in stress tolerance and detoxification metabolism. By examining the gene expression responses to abiotic stresses and enzymatic properties of the ancestral and current proteins, we found that polyploidy-derived GST duplicates show the divergence in enzymatic activities. Through site-directed mutagenesis of ancestral proteins, this study revealed that nonsynonymous substitutions of key amino acid sites play an important role in the divergence of enzymatic functions of polyploidy-derived GST duplicates. These findings provide new insights into the evolutionary and functional dynamics of polyploidy-derived duplicate genes.

Functional divergence of the glutathione S-transferase supergene family in Physcomitrella patens reveals complex patterns of large gene family evolution in land plants.

Plant glutathione S-transferases (GSTs) are multifunctional proteins encoded by a large gene family that play major roles in the detoxification of xenobiotics and oxidative stress metabolism. To date, studies on the GST gene family have focused mainly on vascular plants (particularly agricultural plants). In contrast, little information is available on the molecular characteristics of this large gene family in nonvascular plants. In addition, the evolutionary patterns of this family in land plants remain unclear. In this study, we identified 37 GST genes from the whole genome of the moss Physcomitrella patens, a nonvascular representative of early land plants. The 37 P. patens GSTs were divided into 10 classes, including two new classes (hemerythrin and iota). However, no tau GSTs were identified, which represent the largest class among vascular plants. P. patens GST gene family members showed extensive functional divergence in their gene structures, gene expression responses to abiotic stressors, enzymatic characteristics, and the subcellular locations of the encoded proteins. A joint phylogenetic analysis of GSTs from P. patens and other higher vascular plants showed that different class GSTs had distinct duplication patterns during the evolution of land plants. By examining multiple characteristics, this study revealed complex patterns of evolutionary divergence among the GST gene family in land plants.

Discovery of specific catalytic activity toward IAA/FA by LaSABATHs based on genome-wide phylogenetic and enzymatic analysis of SABATH gene family from Larix kaempferi.

The SABATH methyltransferases catalyze methylation of small-molecule metabolites, which participate in plant growth, development and defense response. Given lack of genome-wide studies on gymnosperms SABATH family, the formation and functional differentiation mechanism of the Larix kaempferi SABATH gene family was systematically and exhaustively explored by analyzing gene sequence characteristics, phylogenetic relationship, expression pattern, and enzyme activities. Phylogenetic analysis showed that 247 SABATH genes from 14 land plants were divided into 4 clades, and lineage-specific gene duplication events were important factors that contributed to the evolution of the SABATH gene family in gymnosperms and angiosperms. Substrate specificity analysis of 18 Larix SABATH proteins showed that LaSABATHs could catalyze O-methylation of indole-3-acetic acid (IAA) and farnesic acid (FA), N-methylation of theobromine, and S-methylation of thiobenzoic acid. Furthermore, only LaSABATH2 and LaSABATH29 could catalyze O-methylation of FA, and only LaSABATH30 could catalyze O-methylation of IAA. Homology modeling and molecular docking studies showed the hydrogen bond formed between the His188 of LaSABATH30 and IAA and the noticeable hydrophobic IAA-binding pocket may be helpful for IAA methylation. In this study, identification of proteins with significant specific catalytic activity toward FA and IAA provided high-quality candidate genes for forest genetics and breeding.

The complete chloroplast genome of Chinese endemic species Abies ferreana (Pinaceae) and its phylogenetic analysis.

Abies ferreana Bordères & Gaussen 1947 is endemic to China, where it is distributed at 3300-4000 meters in the mountains of Southwest Sichuan and Northwest Yunnan. In this study, the complete chloroplast genome of A. ferreana was reconstructed by de novo assembly using whole-genome sequencing data. The complete chloroplast genome of A. ferreana was 120,049 bp in length with a GC content of 37.9%. A total of 113 genes were identified, including 4 rRNA genes, 35 tRNA genes, and 74 protein-coding genes. Among these, 14 genes contain introns. In the phylogenetic tree with 12 other species of Abies, A. ferreana and Abies fanjingshanensis W. L. Huang et al. 1984 were grouped into the same branch, with a bootstrap value of 100%. The complete chloroplast genome of A. ferreana provides potential genetic resources for further Abies evolutionary and genomic studies.

Molecular evolution and functional characterization of chitinase gene family in Populus trichocarpa.

Chitinases, the chitin-degrading enzymes, have been shown to play important role in defense against the chitin-containing fungal pathogens. In this study, we identified 48 chitinase-coding genes from the woody model plant Populus trichocarpa. Based on phylogenetic analysis, the Populus chitinases were classified into seven groups. Different gene structures and protein domain architectures were found among the seven Populus chitinase groups. Selection pressure analysis indicated that all the seven groups are under purifying selection. Phylogenetic analysis combined with chromosome location analysis showed that Populus chitinase gene family mainly expanded through tandem duplication. The Populus chitinase gene family underwent marked expression divergence and is inducibly expressed in response to treatments, such as chitosan, chitin, salicylic acid and methyl jasmonate. Protein enzymatic activity analysis showed that Populus chitinases had activity towards both chitin and chitosan. By integrating sequence characteristic, phylogenetic, selection pressure, gene expression and protein activity analysis, this study shed light on the evolution and function of chitinase family in poplar.

Genome-wide analysis of superoxide dismutase genes in Larix kaempferi.

Superoxide dismutase is a key enzyme that scavenges superoxide anion and plays vital roles in plant antioxidant system. This study identified six SOD genes from the deciduous conifer Larix kaempferi, which is widely distributed across the cooler regions of the northern hemisphere. These six SOD genes were classified into three types: Cu/Zn-SOD (LkSOD1, 2, 3 and 4), Fe-SOD (LkSOD5) and Mn-SODs (LkSOD6). Three Cu/Zn-SOD proteins (LkSOD1, 3 and 4) were cytosolic-localized, while the other three proteins (LkSOD2, 5 and 6) were chloroplast-localized. Larix SOD proteins displayed catalytic activities toward superoxide anion, and retained >55% of its maximum enzymatic activity between 10 °C and 40 °C. Over expressions of three Larix SOD genes (LkSOD2, 4 and 6) in Arabidopsis thaliana, respectively, showed increased germination rates and root lengths during salt stress. LkSOD5 gene could rescue pale green and dwarf phenotype of Arabidopsis atfsd2-2 mutant. Taken together, this study provided comprehensive insight into the gymnosperm SOD gene family.

Genome-wide profiling of expression and biochemical functions of the Medicago glutathione S-transferase gene family.

Glutathione S-transferases are ubiquitous enzyme in plants, playing vital roles in several physiological and developmental processes. In this study we identified 73 GST genes from the genome of Medicago truncatula. The Medicago GSTs were divided to eight classes with tau and phi being the most numerous. Six clusters were found on four Medicago chromosomes. The local gene duplication mainly contributed to the expansion of this large gene family. Functional divergence was found in their gene structures, gene expression patterns, and enzyme properties. A genomic comparative analysis revealed lineage-specific loss/gain events between Medicago and Glycine. This study offered new insights into the evolution of gene family between closely related species.

The complete chloroplast genome sequence of Populus wilsonii and its phylogenetic analysis.

The complete chloroplast genome of Populus wilsonii was reconstructed by reference-based assembly using whole-genome sequencing data. The total chloroplast genome size of P. wilsonii was 158,080 bp in length, including a pair of inverted repeat regions (IRs) of 27,749 bp each, a large single-copy region (LSC) of 85,949 bp and a small single-copy region (SSC) of 16,633 bp. A total of 133 genes were predicted from the chloroplast genome, including 86 protein-coding genes, 39 tRNA genes and eight rRNA genes. Among these genes, 20 genes occurred in IRs, containing nine protein-coding genes, seven tRNA genes and four rRNA genes. The GC content of P. wilsonii chloroplast genome was 36.6%. The phylogenetic analysis with 15 other species showed that P. wilsonii was closely clustered with Populus cathayana. The complete chloroplast genome of P. wilsonii provides new insights into Populus evolutionary and genomic studies.