Cytochrome b5 diversity in green lineages preceded the evolution of syringyl lignin biosynthesis.

Lignin production marked a milestone in vascular plant evolution, and the emergence of syringyl (S) lignin is lineage specific. S-lignin biosynthesis in angiosperms, mediated by ferulate 5-hydroxylase (F5H, CYP84A1), has been considered a recent evolutionary event. F5H uniquely requires the cytochrome b5 protein CB5D as an obligatory redox partner for catalysis. However, it remains unclear how CB5D functionality originated and whether it coevolved with F5H. We reveal here the ancient evolution of CB5D-type function supporting F5H-catalyzed S-lignin biosynthesis. CB5D emerged in charophyte algae, the closest relatives of land plants, and is conserved and proliferated in embryophytes, especially in angiosperms, suggesting functional diversification of the CB5 family before terrestrialization. A sequence motif containing acidic amino residues in Helix 5 of the CB5 heme-binding domain contributes to the retention of CB5D function in land plants but not in algae. Notably, CB5s in the S-lignin-producing lycophyte Selaginella lack these residues, resulting in no CB5D-type function. An independently evolved S-lignin biosynthetic F5H (CYP788A1) in Selaginella relies on NADPH-dependent cytochrome P450 reductase as sole redox partner, distinct from angiosperms. These results suggest that angiosperm F5Hs coopted the ancient CB5D, forming a modern cytochrome P450 monooxygenase system for aromatic ring meta-hydroxylation, enabling the reemergence of S-lignin biosynthesis in angiosperms.

Divergence of active site motifs among different classes of Populus glutaredoxins results in substrate switches.

Enzymes are essential components of all biological systems. The key characteristics of proteins functioning as enzymes are their substrate specificities and catalytic efficiencies. In plants, most genes encoding enzymes are members of large gene families. Within such families, the contributions of active site motifs to the functional divergence of duplicate genes have not been well elucidated. In this study, we identified 41 glutaredoxin (GRX) genes in the Populus trichocarpa genome. GRXs are ubiquitous enzymes in plants that play important roles in developmental and stress tolerance processes. In poplar, GRX genes were divided into four classes based on clear differences in gene structure and expression pattern, subcellular localization, enzymatic activity, and substrate specificity of the encoded proteins. Using site-directed mutagenesis, this study revealed that the divergence of the active site motif among different classes of GRX proteins resulted in substrate switches and thus provided new insights into the molecular evolution of these important plant enzymes.

The Larix kaempferi genome reveals new insights into wood properties.

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L.

BACKGROUNDS: Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae. RESULTS: In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists. CONCLUSIONS: The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.

The class II KNOX transcription factors KNAT3 and KNAT7 synergistically regulate monolignol biosynthesis in Arabidopsis.

The function of the transcription factor KNOTTED ARABIDOPSIS THALIANA7 (KNAT7) is still unclear since it appears to be either a negative or a positive regulator for secondary cell wall deposition with its loss-of-function mutant displaying thicker interfascicular and xylary fiber cell walls but thinner vessel cell walls in inflorescence stems. To explore the exact function of KNAT7, class II KNOTTED1-LIKE HOMEOBOX (KNOX II) genes in Arabidopsis including KNAT3, KNAT4, and KNAT5 were studied together. By chimeric repressor technology, we found that both KNAT3 and KNAT7 repressors exhibited a similar dwarf phenotype. Both KNAT3 and KNAT7 genes were expressed in the inflorescence stems and the knat3 knat7 double mutant exhibited a dwarf phenotype similar to the repressor lines. A stem cross-section of knat3 knat7 displayed an enhanced irregular xylem phenotype as compared with the single mutants, and its cell wall thickness in xylem vessels and interfascicular fibers was significantly reduced. Analysis of cell wall chemical composition revealed that syringyl lignin was significantly decreased while guaiacyl lignin was increased in the knat3 knat7 double mutant. Coincidently, the knat3 knat7 transcriptome showed that most lignin pathway genes were activated, whereas the syringyl lignin-related gene Ferulate 5-Hydroxylase (F5H) was down-regulated. Protein interaction analysis revealed that KNAT3 and KNAT7 can form a heterodimer, and KNAT3, but not KNAT7, can interact with the key secondary cell wall formation transcription factors NST1/2, which suggests that the KNAT3-NST1/2 heterodimer complex regulates F5H to promote syringyl lignin synthesis. These results indicate that KNAT3 and KNAT7 synergistically work together to promote secondary cell wall biosynthesis.

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.

Three UDP-xylose transporters participate in xylan biosynthesis by conveying cytosolic UDP-xylose into the Golgi lumen in Arabidopsis.

UDP-xylose (UDP-Xyl) is synthesized by UDP-glucuronic acid decarboxylases, also termed UDP-Xyl synthases (UXSs). The Arabidopsis genome encodes six UXSs, which fall into two groups based upon their subcellular location: the Golgi lumen and the cytosol. The latter group appears to play an important role in xylan biosynthesis. Cytosolic UDP-Xyl is transported into the Golgi lumen by three UDP-Xyl transporters (UXT1, 2, and 3). However, while single mutants affected in the UDP-Xyl transporter 1 (UXT1) showed a substantial reduction in cell wall xylose content, a double mutant affected in UXT2 and UXT3 had no obvious effect on cell wall xylose deposition. This prompted us to further investigate redundancy among the members of the UXT family. Multiple uxt mutants were generated, including a triple mutant, which exhibited collapsed vessels and reduced cell wall thickness in interfascicular fiber cells. Monosaccharide composition, molecular weight, nuclear magnetic resonance, and immunolabeling studies demonstrated that both xylan biosynthesis (content) and fine structure were significantly affected in the uxt triple mutant, leading to phenotypes resembling those of the irx mutants. Pollination was also impaired in the uxt triple mutant, likely due to reduced filament growth and anther dehiscence caused by alterations in the composition of the cell walls. Moreover, analysis of the nucleotide sugar composition of the uxt mutants indicated that nucleotide sugar interconversion is influenced by the cytosolic UDP-Xyl pool within the cell. Taken together, our results underpin the physiological roles of the UXT family in xylan biosynthesis and provide novel insights into the nucleotide sugar metabolism and trafficking in plants.

Structural insight into the type-II mitochondrial NADH dehydrogenases.

The single-component type-II NADH dehydrogenases (NDH-2s) serve as alternatives to the multisubunit respiratory complex I (type-I NADH dehydrogenase (NDH-1), also called NADH:ubiquinone oxidoreductase; EC 1.6.5.3) in catalysing electron transfer from NADH to ubiquinone in the mitochondrial respiratory chain. The yeast NDH-2 (Ndi1) oxidizes NADH on the matrix side and reduces ubiquinone to maintain mitochondrial NADH/NAD(+) homeostasis. Ndi1 is a potential therapeutic agent for human diseases caused by complex I defects, particularly Parkinson's disease, because its expression restores the mitochondrial activity in animals with complex I deficiency. NDH-2s in pathogenic microorganisms are viable targets for new antibiotics. Here we solve the crystal structures of Ndi1 in its substrate-free, NADH-, ubiquinone- and NADH-ubiquinone-bound states, to help understand the catalytic mechanism of NDH-2s. We find that Ndi1 homodimerization through its carboxy-terminal domain is critical for its catalytic activity and membrane targeting. The structures reveal two ubiquinone-binding sites (UQ(I) and UQ(II)) in Ndi1. NADH and UQ(I) can bind to Ndi1 simultaneously to form a substrate-protein complex. We propose that UQ(I) interacts with FAD to act as an intermediate for electron transfer, and that NADH transfers electrons through this FAD-UQ(I) complex to UQ(II). Together our data reveal the regulatory and catalytic mechanisms of Ndi1 and may facilitate the development or targeting of NDH-2s for potential therapeutic applications.

Subcellular Relocalization and Positive Selection Play Key Roles in the Retention of Duplicate Genes of Populus Class III Peroxidase Family.

Gene duplication is the primary source of new genes and novel functions. Over the course of evolution, many duplicate genes lose their function and are eventually removed by deletion. However, some duplicates have persisted and evolved diverse functions. A particular challenge is to understand how this diversity arises and whether positive selection plays a role. In this study, we reconstructed the evolutionary history of the class III peroxidase (PRX) genes from the Populus trichocarpa genome. PRXs are plant-specific enzymes that play important roles in cell wall metabolism and in response to biotic and abiotic stresses. We found that two large tandem-arrayed clusters of PRXs evolved from an ancestral cell wall type PRX to vacuole type, followed by tandem duplications and subsequent functional specification. Substitution models identified seven positively selected sites in the vacuole PRXs. These positively selected sites showed significant effects on the biochemical functions of the enzymes. We also found that positive selection acts more frequently on residues adjacent to, rather than directly at, a critical active site of the enzyme, and on flexible regions rather than on rigid structural elements of the protein. Our study provides new insights into the adaptive molecular evolution of plant 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.

Methylation mediated by an anthocyanin, O-methyltransferase, is involved in purple flower coloration in Paeonia.

Anthocyanins are major pigments in plants. Methylation plays a role in the diversity and stability of anthocyanins. However, the contribution of anthocyanin methylation to flower coloration is still unclear. We identified two homologous anthocyanin O-methyltransferase (AOMT) genes from purple-flowered (PsAOMT) and red-flowered (PtAOMT) Paeonia plants, and we performed functional analyses of the two genes in vitro and in vivo. The critical amino acids for AOMT catalytic activity were studied by site-directed mutagenesis. We showed that the recombinant proteins, PsAOMT and PtAOMT, had identical substrate preferences towards anthocyanins. The methylation activity of PsAOMT was 60 times higher than that of PtAOMT in vitro. Interestingly, this vast difference in catalytic activity appeared to result from a single amino acid residue substitution at position 87 (arginine to leucine). There were significant differences between the 35S::PsAOMT transgenic tobacco and control flowers in relation to their chromatic parameters, which further confirmed the function of PsAOMT in vivo. The expression levels of the two homologous AOMT genes were consistent with anthocyanin accumulation in petals. We conclude that AOMTs are responsible for the methylation of cyanidin glycosides in Paeonia plants and play an important role in purple coloration in Paeonia spp.

Structural and functional evolution of positively selected sites in pine glutathione S-transferase enzyme family.

Phylogenetic analyses have identified positive selection as an important driver of protein evolution, both structural and functional. However, the lack of appropriate combined functional and structural assays has generally hindered attempts to elucidate patterns of positively selected sites and their effects on enzyme activity and substrate specificity. In this study we investigated the evolutionary divergence of the glutathione S-transferase (GST) family in Pinus tabuliformis, a pine that is widely distributed from northern to central China, including cold temperate and drought-stressed regions. GSTs play important roles in plant stress tolerance and detoxification. We cloned 44 GST genes from P. tabuliformis and found that 26 of the 44 belong to the largest (Tau) class of GSTs and are differentially expressed across tissues and developmental stages. Substitution models identified five positively selected sites in the Tau GSTs. To examine the functional significance of these positively selected sites, we applied protein structural modeling and site-directed mutagenesis. We found that four of the five positively selected sites significantly affect the enzyme activity and specificity; thus their variation broadens the GST family substrate spectrum. In addition, positive selection has mainly acted on secondary substrate binding sites or sites close to (but not directly at) the primary substrate binding site; thus their variation enables the acquisition of new catalytic functions without compromising the protein primary biochemical properties. Our study sheds light on selective aspects of the functional and structural divergence of the GST family in pine and other organisms.

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.

Molecular evolution and expression divergence of the Populus polygalacturonase supergene family shed light on the evolution of increasingly complex organs in plants.

Plant polygalacturonases (PGs) are involved in cell separation processes during many stages of plant development. Investigation into the diversification of this large gene family in land plants could shed light on the evolution of structural development. We conducted whole-genome annotation, molecular evolution and gene expression analyses of PG genes in five species of land plant: Populus, Arabidopsis, rice, Selaginella and Physcomitrella. We identified 75, 44, 16 and 11 PG genes from Populus, rice, Selaginella and Physcomitrella genomes, respectively, which were divided into three classes. We inferred rapid expansion of class I PG genes in Populus, Arabidopsis and rice, while copy numbers of classes II and III PG genes were relatively conserved in all five species. Populus, Arabidopsis and rice class I PG genes were under more relaxed selection constraints than class II PG genes, while this selective pressure divergence was not observed in Selaginella and Physcomitrella PG families. In addition, class I PG genes underwent marked expression divergence in Populus, rice and Selaginella. Our results suggest that PG gene expansion occurred after the divergence of the lycophytes and euphyllophytes, and this expansion was likely paralleled by the evolution of increasingly complex organs in land plants.

CRISPR/Cas9 mutated p-coumaroyl shikimate 3'-hydroxylase 3 gene in Populus tomentosa reveals lignin functioning on supporting tree upright.

The lignin plays one of the most important roles in plant secondary metabolism. However, it is still unclear how lignin can contribute to the impressive height of wood growth. In this study, C3'H, a rate-limiting enzyme of the lignin pathway, was used as the target gene. C3'H3 was knocked out by CRISPR/Cas9 in Populus tomentosa. Compared with wild-type popular trees, c3'h3 mutants exhibited dwarf phenotypes, collapsed xylem vessels, weakened phloem thickening, decreased hydraulic conductivity and photosynthetic efficiency, and reduced auxin content, except for reduced total lignin content and significantly increased H-subunit lignin. In the c3'h3 mutant, the flavonoid biosynthesis genes CHS, CHI, F3H, DFR, ANR, and LAR were upregulated, and flavonoid metabolite accumulations were detected, indicating that decreasing the lignin biosynthesis pathway enhanced flavonoid metabolic flux. Furthermore, flavonoid metabolites, such as naringenin and hesperetin, were largely increased, while higher hesperetin content suppressed plant cell division. Thus, studying the c3'h3 mutant allows us to deduce that lignin deficiency suppresses tree growth and leads to the dwarf phenotype due to collapsed xylem and thickened phloem, limiting material exchanges and transport.

Transcriptome-wide identification and characterization of miRNAs from Pinus densata.

BACKGROUND: MicroRNAs (miRNAs) play key roles in diverse developmental processes, nutrient homeostasis and responses to biotic and abiotic stresses. The biogenesis and regulatory functions of miRNAs have been intensively studied in model angiosperms, such as Arabidopsis thaliana, Oryza sativa and Populus trichocarpa. However, global identification of Pinus densata miRNAs has not been reported in previous research. RESULTS: Here, we report the identification of 34 conserved miRNAs belonging to 25 miRNA families from a P. densata mRNA transcriptome database using local BLAST and MIREAP programs. The primary and/or precursor sequences of 29 miRNAs were further confirmed by RT-PCR amplification and subsequent sequencing. The average value of the minimal folding free energy indexes of the 34 miRNA precursors was 0.92. Nineteen (58%) mature miRNAs began with a 5' terminal uridine residue. Analysis of miRNA precursors showed that 19 mature miRNAs were novel members of 14 conserved miRNA families, of which 17 miRNAs were further validated by subcloning and sequencing. Using real-time quantitative RT-PCR, we found that the expression levels of 7 miRNAs were more than 2-fold higher in needles than in stems. In addition, 72 P. densata mRNAs were predicted to be targets of 25 miRNA families. Four target genes, including a nodal modulator 1-like protein gene, two GRAS family transcription factor protein genes and one histone deacetylase gene, were experimentally verified to be the targets of 3 P. densata miRNAs, pde-miR162a, pde-miR171a and pde-miR482a, respectively. CONCLUSIONS: This study led to the discovery of 34 conserved miRNAs comprising 25 miRNA families from Pinus densata. These results lay a solid foundation for further studying the regulative roles of miRNAs in the development, growth and responses to environmental stresses in P. densata.

Demography and speciation history of the homoploid hybrid pine Pinus densata on the Tibetan Plateau.

Pinus densata is an ecologically successful homoploid hybrid that inhabits vast areas of heterogeneous terrain on the south-eastern Tibetan Plateau as a result of multiple waves of colonization. Its region of origin, route of colonization onto the plateau and the directions of introgression with its parental species have previously been defined, but little is known about the isolation and divergence history of its populations. In this study, we surveyed nucleotide polymorphism over eight nuclear loci in 19 representative populations of P. densata and its parental species. Using this information and coalescence simulations, we assessed the historical changes in its population size, gene flow and divergence in time and space. The results indicate a late Miocene origin for P. densata associated with the recent uplift of south-eastern Tibet. The subsequent differentiation between geographical regions of this species began in the late Pliocene and was induced by regional topographical changes and Pleistocene glaciations. The ancestral P. densata population had a large effective population size but the central and western populations were established by limited founders, suggesting that there were severe bottlenecks during the westward migration out of the ancestral hybrid zone. After separating from their ancestral populations, population expansion occurred in all geographical regions especially in the western range. Gene flow in P. densata was restricted to geographically neighbouring populations, resulting in significant differentiation between regional groups. The new information on the divergence and demographic history of P. densata reported herein enhances our understanding of its speciation process on the Tibetan Plateau.

Extensive functional diversification of the Populus glutathione S-transferase supergene family.

Identifying how genes and their functions evolve after duplication is central to understanding gene family radiation. In this study, we systematically examined the functional diversification of the glutathione S-transferase (GST) gene family in Populus trichocarpa by integrating phylogeny, expression, substrate specificity, and enzyme kinetic data. GSTs are ubiquitous proteins in plants that play important roles in stress tolerance and detoxification metabolism. Genome annotation identified 81 GST genes in Populus that were divided into eight classes with distinct divergence in their evolutionary rate, gene structure, expression responses to abiotic stressors, and enzymatic properties of encoded proteins. In addition, when all the functional parameters were examined, clear divergence was observed within tandem clusters and between paralogous gene pairs, suggesting that subfunctionalization has taken place among duplicate genes. The two domains of GST proteins appear to have evolved under differential selective pressures. The C-terminal domain seems to have been subject to more relaxed functional constraints or divergent directional selection, which may have allowed rapid changes in substrate specificity, affinity, and activity, while maintaining the primary function of the enzyme. Our findings shed light on mechanisms that facilitate the retention of duplicate genes, which can result in a large gene family with a broad substrate spectrum and a wide range of reactivity toward different substrates.

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.