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 co-evolved 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 co-opted the ancient CB5D, forming a modern cytochrome P450 monooxygenase system for aromatic ring meta-hydroxylation, enabling the re-emergence of S-lignin biosynthesis in angiosperms.

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.

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.

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.

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.

Ectopic expression of Populus MYB10 promotes secondary cell wall thickening and inhibits anthocyanin accumulation.

Secondary cell wall (SCW) formation is regulated by a multilevel transcriptional regulatory network, in which MYB transcription factors (TFs) play key roles. In woody plants, hundreds of MYB TFs have been identified, most of which have unknown functions in wood SCW biosynthesis. Here, we characterized the function of a Populus MYB gene, PtoMYB10. PtoMYB10 was found to encode an R2R3-MYB TF and exhibit dominant expression in xylem tissues. PtoMYB10 was determined to be located in the nucleus with the ability to activate transcription. Overexpression of PtoMYB10 in Populus resulted in a drastic increase in SCW thickening in xylem fiber cells as well as ectopic deposition of lignin in cortex cells. The expression of genes associated with lignin biosynthesis was induced in PtoMYB10 overexpressing plants, whereas repressed gene expression was found with the anthocyanin biosynthesis pathway. Lignin and anthocyanin are both produced from metabolites of the phenylpropanoid pathway. Accordingly, the anthocyanin content of Populus overexpressing PtoMYB10 decreased by more than 68%. These results indicate that PtoMYB10 can positively regulate xylary fiber SCW thickening, accompanied by the reprogramming of phenylpropanoid metabolism, which redirects metabolic flux from anthocyanin biosynthesis to monolignol biosynthesis.

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.

NRAMP1 promotes iron uptake at the late stage of iron deficiency in poplars.

Iron (Fe) is an essential micronutrient for plant survival and proliferation. Plants have evolved complex mechanisms to maintain Fe homeostasis in response to Fe deficiency. In this study, we evaluated the physiological, biochemical and transcriptomic differences between poplars grown under Fe-sufficient and Fe-deficient conditions to elucidate the mechanistic responses of poplars to Fe deficiency. Our results revealed that chlorophyll synthesis and photosynthesis were inhibited under Fe-deficient conditions. The inhibition of these pathways caused chlorosis and reduced shoot growth. Although both photosynthetic systems (PSI and PSII) were inhibited under Fe limitation, PSI was affected more severely and earlier than PSII. Fe deficiency also promoted root growth and increased the accumulation of divalent metal ions in roots. IRT1 and NRAMP1 are both Fe2+ transporters for iron uptake in Arabidopsis. In this study, however, only NRAMP1 was induced to promote Fe2+ uptake in roots at the late stage of Fe deficiency response. It indicated that NRAMP1, rather than the more well-known IRT1, might be a major Fe2+ transporter at the late stage of Fe-deficiency in poplars.

De novo assembly of white poplar genome and genetic diversity of white poplar population in Irtysh River basin in China.

The white poplar (Populus alba) is widely distributed in Central Asia and Europe. There are natural populations of white poplar in Irtysh River basin in China. It also can be cultivated and grown well in northern China. In this study, we sequenced the genome of P. alba by single-molecule real-time technology. De novo assembly of P. alba had a genome size of 415.99 Mb with a contig N50 of 1.18 Mb. A total of 32,963 protein-coding genes were identified. 45.16% of the genome was annotated as repetitive elements. Genome evolution analysis revealed that divergence between P. alba and Populus trichocarpa (black cottonwood) occurred ~5.0 Mya (3.0, 7.1). Fourfold synonymous third-codon transversion (4DTV) and synonymous substitution rate (ks) distributions supported the occurrence of the salicoid WGD event (~ 65 Mya). Twelve natural populations of P. alba in the Irtysh River basin in China were sequenced to explore the genetic diversity. Average pooled heterozygosity value of P. alba populations was 0.170±0.014, which was lower than that in Italy (0.271±0.051) and Hungary (0.264±0.054). Tajima's D values showed a negative distribution, which might signify an excess of low frequency polymorphisms and a bottleneck with later expansion of P. alba populations examined.

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.

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.

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.

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.

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.

Genome-Wide Analysis of the Glutathione S-Transferase Gene Family in Capsella rubella: Identification, Expression, and Biochemical Functions.

Extensive subfunctionalization might explain why so many genes have been maintained after gene duplication, which provides the engine for gene family expansion. However, it is still a particular challenge to trace the evolutionary dynamics and features of functional divergences in a supergene family over the course of evolution. In this study, we identified 49 Glutathione S-transferase (GST) genes from the Capsella rubella, a close relative of Arabidopsis thaliana and a member of the mustard family. Capsella GSTs can be categorized into eight classes, with tau and phi GSTs being the most numerous. The expansion of the two classes mainly occurs through tandem gene duplication, which results in tandem-arrayed gene clusters on chromosomes. By integrating phylogenetic analysis, expression patterns, and biochemical functions of Capsella and Arabidopsis GSTs, functional divergence, both in gene expression and enzymatic properties, were clearly observed in paralogous gene pairs in Capsella (even the most recent duplicates), and orthologous GSTs in Arabidopsis/Capsella. This study provides functional evidence for the expansion and organization of a large gene family in closely related species.

Role of UDP-Glucuronic Acid Decarboxylase in Xylan Biosynthesis in Arabidopsis.

UDP-xylose (UDP-Xyl) is the Xyl donor used in the synthesis of major plant cell-wall polysaccharides such as xylan (as a backbone-chain monosaccharide) and xyloglucan (as a branching monosaccharide). The biosynthesis of UDP-Xyl from UDP-glucuronic acid (UDP-GlcA) is irreversibly catalyzed by UDP-glucuronic acid decarboxylase (UXS). Until now, little has been known about the physiological roles of UXS in plants. Here, we report that AtUXS1, AtUXS2, and AtUXS4 are located in the Golgi apparatus whereas AtUXS3, AtUXS5, and AtUXS6 are located in the cytosol. Although all six single AtUXS T-DNA mutants and the uxs1 usx2 uxs4 triple mutant show no obvious phenotype, the uxs3 uxs5 uxs6 triple mutant has an irregular xylem phenotype. Monosaccharide analysis showed that Xyl levels decreased in uxs3 uxs5 uxs6 and linkage analysis confirmed that the xylan content in uxs3 xus5 uxs6 declined, indicating that UDP-Xyl from cytosol AtUXS participates in xylan synthesis. Gel-permeation chromatography showed that the molecular weight of non-cellulosic polysaccharides in the triple mutants, mainly composed of xylans, is lower than that in the wild type, suggesting an effect on the elongation of the xylan backbone. Upon saccharification treatment stems of the uxs3 uxs5 uxs6 triple mutants released monosaccharides with a higher efficiency than those of the wild type. Taken together, our results indicate that the cytosol UXS plays a more important role than the Golgi-localized UXS in xylan biosynthesis.

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.

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.

Weak crossability barrier but strong juvenile selection supports ecological speciation of the hybrid pine Pinus densata on the Tibetan plateau.

Determining how a new hybrid lineage can achieve reproductive isolation is a key to understanding the process and mechanisms of homoploid hybrid speciation. Here, we evaluated the degree and nature of reproductive isolation between the ecologically successful hybrid species Pinus densata and its parental species P. tabuliformis and P. yunnanensis. We performed interspecific crosses among the three species to assess their crossability. We then conducted reciprocal transplantation experiments to evaluate their fitness differentiation, and to examine how natural populations representing different directions of introgression differ in adaptation. The crossing experiments revealed weak genetic barriers among the species. The transplantation trials showed manifest evidence of local adaptation as the three species all performed best in their native habitats. Pinus densata populations from the western edge of its distribution have evolved a strong local adaptation to the specific habitat in that range; populations representing different directions of introgressants with the two parental species all showed fitness disadvantages in this P. densata habitat. These observations illustrate that premating isolation through selection against immigrants from other habitat types or postzygotic isolation through selection against backcrosses between the three species is strong. Thus, ecological selection in combination with endogenous components and geographic isolation has likely played a significant role in the speciation of P. densata.