Secondary metabolites of the flavonoid pathway participate in plant defense, and bHLH and MYB transcription factors regulate the synthesis of these metabolites. Here, we define the regulatory mechanisms in response to pathogens. Two transcription factors from Populus alba var. pyramidalis, PalbHLH1 and PalMYB90, were overexpressed together in poplar, and transcriptome analysis revealed differences in response to pathogen infection. The transgenic plants showed elevated levels of several key flavonoid pathway components: total phenols, proanthocyanidins (PAs), and anthocyanins and intermediates quercetin and kaempferol. Furthermore, PalbHLH1 and PalMYB90 overexpression in poplar enhanced antioxidase activities and H₂O₂ release and also increased resistance to Botrytis cinerea and Dothiorella gregaria infection. Gene expression profile analysis showed most genes involved in the flavonoid biosynthesis pathway or antioxidant response to be upregulated in MYB90/bHLH1-OE poplar, but significant differential expression occurred in response to pathogen infection. Specifically, expression of PalF3H (flavanone 3-hydroxylase), PalDFR (dihydroflavonol 4-seductase), PalANS (anthocyanin synthase), and PalANR (anthocyanin reductase), which function in initial, middle, and final steps of anthocyanin and PA biosynthesis, respectively, was significantly upregulated in D. gregaria-infected MYB90/bHLH1-OE poplar. Our results highlight that PalbHLH1 and PalMYB90 function as transcriptional activators of flavonoid pathway secondary-metabolite synthesis genes, with differential mechanisms in response to bacterial or fungal infection.

Copy number variation (CNV) represents an important source of genetic divergence that can produce drastic phenotypic differences and may therefore be subject to selection during domestication and environmental adaptation. To investigate the evolutionary dynamics of CNV in the yak genome, we used a read depth approach to detect CNV based on genome resequencing data from 14 wild and 65 domestic yaks and determined CNV regions related to domestication and adaptations to high-altitude. We identified 2,634 CNV regions (CNVRs) comprising a total of 153 megabases (5.7 % of the yak genome) and 3,879 overlapping annotated genes. Comparison between domestic and wild yak populations identified 121 potentially selected CNVRs, harboring genes related to neuronal development, reproduction, nutrition and energy metabolism. In addition, we found 85 CNVRs that are significantly different between domestic yak living in high- and low-altitude areas, including three genes related to hypoxia response and six related to immune defense. This analysis shows that genic CNVs may play an important role in phenotypic changes during yak domestication and adaptation to life at high-altitude. We present the first refined CNV map for yak along with comprehensive genomic analysis of yak CNV. Our results provide new insights into the genetic basis of yak domestication and adaptation to living in a high-altitude environment, as well as a valuable genetic resource that will facilitate future CNV association studies of important traits in yak and other bovid species.

Functional divergence in paralogs is an important genetic source of evolutionary innovation. Actin-depolymerizing factors (ADFs) are among the most important actin binding proteins and are involved in generating and remodeling actin cytoskeletal architecture via their conserved F-actin severing or depolymerizing activity. In plants, ADFs coevolved with actin, but their biochemical properties are diverse. Unfortunately, the biochemical function of most plant ADFs and the potential mechanisms of their functional divergence remain unclear. Here, in vitro biochemical analyses demonstrated that all 11 ADF genes in Arabidopsis thaliana exhibit opposing biochemical properties. Subclass III ADFs evolved F-actin bundling (B-type) function from conserved F-actin depolymerizing (D-type) function, and subclass I ADFs have enhanced D-type function. By tracking historical mutation sites on ancestral proteins, several fundamental amino acid residues affecting the biochemical functions of these proteins were identified in Arabidopsis and various plants, suggesting that the biochemical divergence of ADFs has been conserved during the evolution of angiosperm plants. Importantly, N-terminal extensions on subclass III ADFs that arose from intron-sliding events are indispensable for the alteration of D-type to B-type function. We conclude that the evolution of these N-terminal extensions and several conserved mutations produced the diverse biochemical functions of plant ADFs from a putative ancestor.

Drought severely limits plant development and growth; accordingly, plants have evolved strategies to prevent water loss and adapt to water deficit conditions. However, experimental cases that corroborate these evolutionary processes are limited. The LACCASEs (LACs) family is involved in various plant development and growth processes. Here, we performed an evolutionary analysis ofLACsfromPopulus euphraticaand characterized the functions ofLACsinArabidopsisand poplar. The results showed that inPeuLACs, multiple gene duplications led to apparent functional redundancy as the result of various selective pressures. Among them,PeuLAC2underwent strong positive selection. Heterologous expression analyses showed that the overexpression ofPeuLAC2alters the xylem structure of plants, including thickening the secondary cell wall (SCW) and increasing the fiber cell length and stem tensile strength. Altogether, these changes improve the water transport capacity of plants. The analysis of the physiological experimental results showed thatPeuLAC2-OE lines exhibited a stronger antioxidant response and greater drought tolerance than WT. Three genes screened by transcriptome analysis,NAC025,BG1, andUGT, that are associated with SCW synthesis and drought stress were all upregulated in thePeuLAC2-OE lines, implying that the overexpression ofPeuLAC2thickened the SCW, improved the water transport capacity of the plant, and further enhanced its drought tolerance. Our study highlights that genes that have undergone adaptive evolution may participate in the development of adaptive traits inP. euphraticaand thatPeuLAC2could be a candidate gene for molecular genetic breeding in trees.

4-Coumarate:CoA ligase (4CL) genes are critical for the biosynthesis of plant phenylpropanoids. Here we identified 20 4CL genes in the genomes of two desert poplars (Populus euphratica and P. pruinosa) and salt-sensitive congener (P. trichocarpa), but 12 in Salix suchowensis (Salix willow). Phylogenetic analyses clustered all Salicaceae 4CL genes into two clades, and one of them (corresponding to the 4CL-like clade from Arabidopsis) showed signals of adaptive evolution, with more genes retained in Populus than Salix and Arabidopsis. We also found that 4CL12 (in 4CL-like clade) showed positive selection along the two desert poplar lineages. Transcriptional profiling analyses indicated that the expression of 4CL2, 4CL11, and 4CL12 changed significantly in one or both desert poplars in response to salt stress compared to that of in P. trichocarpa. Our results suggest that the evolution of the 4CL genes may have contributed to the development of salt tolerance in the two desert poplars.