Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus

PagPXYs improve drought tolerance by regulating reactive oxygen species homeostasis in the cambium of Populus alba × P. glandulosa

PXY (Phloem intercalated with xylem) is a receptor kinase required for directional cell division during the development of plant vascular tissue. Drought stress usually affects plant stem cell division and differentiation thereby limiting plant growth. However, the role of PXY in cambial activities of woody plants under drought stress is unclear. In this study, we analyzed the biological functions of two PXY genes (PagPXYa and PagPXYb) in poplar growth and development and in response to drought stress in a hybrid poplar (Populus alba × P. glandulosa, '84K'). Expression analysis indicated that PagPXYs, similar to their orthologs PtrPXYs in Populus trichocarpa, are mainly expressed in the stem vascular system, and related to drought. Interestingly, overexpression of PagPXYa and PagPXYb in poplar did not have a significant impact on the growth status of transgenic plants under normal condition. However, when treated with 8 % PEG6000 or 100 mM H2O2, PagPXYa and PagPXYb overexpressing lines consistently exhibited more cambium cell layers, fewer xylem cell layers, and enhanced drought tolerance compared to the non-transgenic control '84K'. In addition, PagPXYs can alleviate the damage caused by H2O2 to the cambium under drought stress, thereby maintaining the cambial division activity of poplar under drought stress, indicating that PagPXYs play an important role in plant resistance to drought stress. This study provides a new insight for further research on the balance of growth and drought tolerance in forest trees.

Evaluation of novel promoters for vascular tissue-specific gene expression in Populus

Due to the extended generation cycle of trees, the breeding process for forest trees tends to be time-consuming. Genetic engineering has emerged as a viable approach to expedite the genetic breeding of forest trees. However, current genetic engineering techniques employed in forest trees often utilize continuous expression promoters such as CaMV 35S, which may result in unintended consequences by introducing genes into non-target tissues. Therefore, it is imperative to develop specific promoters for forest trees to facilitate targeted and precise design and breeding. In this study, we utilized single-cell RNA-Seq data and co-expression network analysis during wood formation to identify three vascular tissue-specific genes in poplar, PP2-A10, PXY, and VNS07, which are expressed in the phloem, cambium/expanding xylem, and mature xylem, respectively. Subsequently, we cloned the promoters of these three genes from '84K' poplar and constructed them into a vector containing the eyGFPuv visual selection marker, along with the 35S mini enhancer to drive GUS gene expression. Transgenic poplars expressing the ProPagPP2-A10::GUS, ProPagPXY::GUS, and ProPagVNS07::GUS constructs were obtained. To further elucidate the tissue specificity of these promoters, we employed qPCR, histochemical staining, and GUS enzyme activity. Our findings not only establish a solid foundation for the future utilization of these promoters to precisely express of specific functional genes in stems but also provide a novel perspective for the modular breeding of forest trees.

Integrated transcriptomic and CGAs analysis revealed IbGLK1 is a key transcription factor for chlorogenic acid accumulation in sweetpotato (Ipomoea batatas [L.] Lam.) blades

Sweetpotato blades are rich in the functional secondary metabolite chlorogenic acid (CGA), which deepen potential for effective utilization of the blade in industry. In this study, we evaluated the type and content of CGA in the blades of 16 sweetpotato genotypes and analyzed the correlation between CGA content and antioxidant capacity. Then we isolated and characterized IbGLK1, a GARP-type transcription factor, by comparative transcriptome analysis. A subcellular localization assay indicated that IbGLK1 is located in the nucleus. Overexpression and silencing of IbGLK1 in sweetpotato blade resulted in a 0.90-fold increase and 1.84-fold decrease, respectively, in CGA content compared to the control. Yeast one-hybrid and dual-luciferase assays showed that IbGLK1 binds and activates the promoters of IbHCT, IbHQT, IbC4H, and IbUGCT, resulting in the promotion of CGA biosynthesis. In conclusion, our study provides insights into a high-quality gene for the regulation of CGA metabolism and germplasm resources for breeding sweetpotato.

Genome-wide association study and network analysis of in vitro transformation in Populus trichocarpa support key roles of diverse phytohormone pathways and cross talk

Wide variation in amenability to transformation and regeneration (TR) among many plant species and genotypes presents a challenge to the use of genetic engineering in research and breeding. To help understand the causes of this variation, we performed association mapping and network analysis using a population of 1204 wild trees of Populus trichocarpa (black cottonwood). To enable precise and high-throughput phenotyping of callus and shoot TR, we developed a computer vision system that cross-referenced complementary red, green, and blue (RGB) and fluorescent-hyperspectral images. We performed association mapping using single-marker and combined variant methods, followed by statistical tests for epistasis and integration of published multi-omic datasets to identify likely regulatory hubs. We report 409 candidate genes implicated by associations within 5 kb of coding sequences, and epistasis tests implicated 81 of these candidate genes as regulators of one another. Gene ontology terms related to protein-protein interactions and transcriptional regulation are overrepresented, among others. In addition to auxin and cytokinin pathways long established as critical to TR, our results highlight the importance of stress and wounding pathways. Potential regulatory hubs of signaling within and across these pathways include GROWTH REGULATORY FACTOR 1 (GRF1), PHOSPHATIDYLINOSITOL 4-KINASE β1 (PI-4Kβ1), and OBF-BINDING PROTEIN 1 (OBP1).

An integrated metagenomic, metabolomic and transcriptomic survey of Populus across genotypes and environments

Bridging molecular information to ecosystem-level processes would provide the capacity to understand system vulnerability and, potentially, a means for assessing ecosystem health. Here, we present an integrated dataset containing environmental and metagenomic information from plant-associated microbial communities, plant transcriptomics, plant and soil metabolomics, and soil chemistry and activity characterization measurements derived from the model tree species Populus trichocarpa. Soil, rhizosphere, root endosphere, and leaf samples were collected from 27 different P. trichocarpa genotypes grown in two different environments leading to an integrated dataset of 318 metagenomes, 98 plant transcriptomes, and 314 metabolomic profiles that are supported by diverse soil measurements. This expansive dataset will provide insights into causal linkages that relate genomic features and molecular level events to system-level properties and their environmental influences.

GWAS supported by computer vision identifies large numbers of candidate regulators of in planta regeneration in Populus trichocarpa

Plant regeneration is an important dimension of plant propagation and a key step in the production of transgenic plants. However, regeneration capacity varies widely among genotypes and species, the molecular basis of which is largely unknown. Association mapping methods such as genome-wide association studies (GWAS) have long demonstrated abilities to help uncover the genetic basis of trait variation in plants; however, the performance of these methods depends on the accuracy and scale of phenotyping. To enable a large-scale GWAS of in planta callus and shoot regeneration in the model tree Populus, we developed a phenomics workflow involving semantic segmentation to quantify regenerating plant tissues over time. We found that the resulting statistics were of highly non-normal distributions, and thus employed transformations or permutations to avoid violating assumptions of linear models used in GWAS. We report over 200 statistically supported quantitative trait loci (QTLs), with genes encompassing or near to top QTLs including regulators of cell adhesion, stress signaling, and hormone signaling pathways, as well as other diverse functions. Our results encourage models of hormonal signaling during plant regeneration to consider keystone roles of stress-related signaling (e.g. involving jasmonates and salicylic acid), in addition to the auxin and cytokinin pathways commonly considered. The putative regulatory genes and biological processes we identified provide new insights into the biological complexity of plant regeneration, and may serve as new reagents for improving regeneration and transformation of recalcitrant genotypes and species.

Identification of environment-insensitive genes for oil content by combination of transcriptome and genome-wide association analysis in rapeseed

BACKGROUND

The primary objective of rapeseed breeding is to enhance oil content, which is predominantly influenced by environmental factors. However, the molecular mechanisms underlying the impact of these environmental factors on oil accumulation remain inadequately elucidated. In this study, we used transcriptome data from two higher (HOC) and two lower oil content (LOC) inbred lines at 35 days after pollination (DAP) to investigate genes exhibiting stable expression across three different environments. Meanwhile, a genome-wide association study (GWAS) was utilized to detect candidate genes exhibiting significant associations with seed oil content across three distinct environments.

RESULTS

The study found a total of 405 stable differentially expressed genes (DEGs), including 25 involved in lipid/fatty acid metabolism and 14 classified as transcription factors. Among these genes, BnBZIP10-A09, BnMYB61-A06, BnAPA1-A08, BnPAS2-A10, BnLCAT3-C05 and BnKASIII-C09 were also found to exhibit significant associations with oil content across multiple different environments based on GWAS of 50 re-sequenced semi-winter rapeseed inbred lines and previously reported intervals. Otherwise, we revealed the presence of additive effects among BnBZIP10-A09, BnKASIII-C09, BnPAS2-A10 and BnAPA1-A08, resulting in a significant increase in seed oil content. Meanwhile, the majority of these stable DEGs are interconnected either directly or indirectly through co-expression network analysis, thereby giving rise to an elaborate molecular network implicated in the potential regulation of seed oil accumulation and stability.

CONCLUSIONS

The combination of transcription and GWAS revealed that natural variation in six environment-insensitive gene regions exhibited significant correlations with seed oil content phenotypes. These results provide important molecular marker information for us to further improve oil content accumulation and stability in rapeseed.

Editing Metabolism, Sex, and Microbiome: How Can We Help Poplar Resist Pathogens?

Poplar (Populus) is a genus of woody plants of great economic value. Due to the growing economic importance of poplar, there is a need to ensure its stable growth by increasing its resistance to pathogens. Genetic engineering can create organisms with improved traits faster than traditional methods, and with the development of CRISPR/Cas-based genome editing systems, scientists have a new highly effective tool for creating valuable genotypes. In this review, we summarize the latest research data on poplar diseases, the biology of their pathogens and how these plants resist pathogens. In the final section, we propose to plant male or mixed poplar populations; consider the genes of the MLO group, transcription factors of the WRKY and MYB families and defensive proteins BbChit1, LJAMP2, MsrA2 and PtDef as the most promising targets for genetic engineering; and also pay attention to the possibility of microbiome engineering.

Overexpression of REDUCED WALL ACETYLATION C increases xylan acetylation and biomass recalcitrance in Populus

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Ectomycorrhizal symbiosis prepares its host locally and systemically for abiotic cue signaling

Tree growth and survival are dependent on their ability to perceive signals, integrate them, and trigger timely and fitted molecular and growth responses. While ectomycorrhizal symbiosis is a predominant tree-microbe interaction in forest ecosystems, little is known about how and to what extent it helps trees cope with environmental changes. We hypothesized that the presence of Laccaria bicolor influences abiotic cue perception by Populus trichocarpa and the ensuing signaling cascade. We submitted ectomycorrhizal or non-ectomycorrhizal P. trichocarpa cuttings to short-term cessation of watering or ozone fumigation to focus on signaling networks before the onset of any physiological damage. Poplar gene expression, metabolite levels, and hormone levels were measured in several organs (roots, leaves, mycorrhizas) and integrated into networks. We discriminated the signal responses modified or maintained by ectomycorrhization. Ectomycorrhizas buffered hormonal changes in response to short-term environmental variations systemically prepared the root system for further fungal colonization and alleviated part of the root abscisic acid (ABA) signaling. The presence of ectomycorrhizas in the roots also modified the leaf multi-omics landscape and ozone responses, most likely through rewiring of the molecular drivers of photosynthesis and the calcium signaling pathway. In conclusion, P. trichocarpa-L. bicolor symbiosis results in a systemic remodeling of the host's signaling networks in response to abiotic changes. In addition, ectomycorrhizal, hormonal, metabolic, and transcriptomic blueprints are maintained in response to abiotic cues, suggesting that ectomycorrhizas are less responsive than non-mycorrhizal roots to abiotic challenges.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.)

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

Plant-microbe interactions: Plant-exuded myo-inositol attracts specific bacterial taxa

Plants exude a plethora of metabolites that transform the microbiome composition. Initiated from genome-wide association studies of either a plant or a bacterium, two new studies dissect the impact of plant-secreted myo-inositol on recruitment of certain bacterial taxa by Arabidopsis.

Plant myo-inositol transport influences bacterial colonization phenotypes

Plant microbiomes are assembled and modified through a complex milieu of biotic and abiotic factors. Despite dynamic and fluctuating contributing variables, specific host metabolites are consistently identified as important mediators of microbial interactions. We combine information from a large-scale metatranscriptomic dataset from natural poplar trees and experimental genetic manipulation assays in seedlings of the model plant Arabidopsis thaliana to converge on a conserved role for transport of the plant metabolite myo-inositol in mediating host-microbe interactions. While microbial catabolism of this compound has been linked to increased host colonization, we identify bacterial phenotypes that occur in both catabolism-dependent and -independent manners, suggesting that myo-inositol may additionally serve as a eukaryotic-derived signaling molecule to modulate microbial activities. Our data suggest host control of this compound and resulting microbial behavior are important mechanisms at play surrounding the host metabolite myo-inositol.

Few-Shot Learning Enables Population-Scale Analysis of Leaf Traits in Populus trichocarpa

Plant phenotyping is typically a time-consuming and expensive endeavor, requiring large groups of researchers to meticulously measure biologically relevant plant traits, and is the main bottleneck in understanding plant adaptation and the genetic architecture underlying complex traits at population scale. In this work, we address these challenges by leveraging few-shot learning with convolutional neural networks to segment the leaf body and visible venation of 2,906 Populus trichocarpa leaf images obtained in the field. In contrast to previous methods, our approach (a) does not require experimental or image preprocessing, (b) uses the raw RGB images at full resolution, and (c) requires very few samples for training (e.g., just 8 images for vein segmentation). Traits relating to leaf morphology and vein topology are extracted from the resulting segmentations using traditional open-source image-processing tools, validated using real-world physical measurements, and used to conduct a genome-wide association study to identify genes controlling the traits. In this way, the current work is designed to provide the plant phenotyping community with (a) methods for fast and accurate image-based feature extraction that require minimal training data and (b) a new population-scale dataset, including 68 different leaf phenotypes, for domain scientists and machine learning researchers. All of the few-shot learning code, data, and results are made publicly available.

Validation of a metabolite-GWAS network for Populus trichocarpa family 1 UDP-glycosyltransferases

Metabolite genome-wide association studies (mGWASs) are increasingly used to discover the genetic basis of target phenotypes in plants such as Populus trichocarpa, a biofuel feedstock and model woody plant species. Despite their growing importance in plant genetics and metabolomics, few mGWASs are experimentally validated. Here, we present a functional genomics workflow for validating mGWAS-predicted enzyme-substrate relationships. We focus on uridine diphosphate-glycosyltransferases (UGTs), a large family of enzymes that catalyze sugar transfer to a variety of plant secondary metabolites involved in defense, signaling, and lignification. Glycosylation influences physiological roles, localization within cells and tissues, and metabolic fates of these metabolites. UGTs have substantially expanded in P. trichocarpa, presenting a challenge for large-scale characterization. Using a high-throughput assay, we produced substrate acceptance profiles for 40 previously uncharacterized candidate enzymes. Assays confirmed 10 of 13 leaf mGWAS associations, and a focused metabolite screen demonstrated varying levels of substrate specificity among UGTs. A substrate binding model case study of UGT-23 rationalized observed enzyme activities and mGWAS associations, including glycosylation of trichocarpinene to produce trichocarpin, a major higher-order salicylate in P. trichocarpa. We identified UGTs putatively involved in lignan, flavonoid, salicylate, and phytohormone metabolism, with potential implications for cell wall biosynthesis, nitrogen uptake, and biotic and abiotic stress response that determine sustainable biomass crop production. Our results provide new support for in silico analyses and evidence-based guidance for in vivo functional characterization.

Risk Assessment of Chlorogenic and Isochlorogenic Acids in Coffee By-Products

Chlorogenic and isochlorogenic acids are naturally occurring antioxidant dietary polyphenolic compounds found in high concentrations in plants, fruits, vegetables, coffee, and coffee by-products. The objective of this review was to assess the potential health risks associated with the oral consumption of coffee by-products containing chlorogenic and isochlorogenic acids, considering both acute and chronic exposure. An electronic literature search was conducted, revealing that 5-caffeoylquinic acid (5-CQA) and 3,5-dicaffeoylquinic acid (3,5-DCQA) are the major chlorogenic acids found in coffee by-products. Toxicological, pharmacokinetic, and clinical data from animal and human studies were available for the assessment, which indicated no significant evidence of toxic or adverse effects following acute oral exposure. The current state of knowledge suggests that long-term exposure to chlorogenic and isochlorogenic acids by daily consumption does not appear to pose a risk to human health when observed at doses within the normal range of dietary exposure. As a result, the intake of CQAs from coffee by-products can be considered reasonably safe.

Integrated Transcriptomic and Metabolomic Analyses Reveal the Effects of Grafting on Special Metabolites of Acanthopanax senticosus Leaves

Acanthopanax senticosus (A. senticosus) is a member of Acanthopanax Miq. and is used in traditional Chinese medicine, and it has been found that grafting technology can be used to alter plant metabolite composition and transcriptome characteristics. In this study, shoots of A. senticosus were grafted onto the rootstocks of the vigorous Acanthopanax sessiliflorus (A. sessiliflorus) to improve its varietal characteristics. In order to investigate the changes in metabolites and transcriptional patterns in grafted A. senticosus leaves (GSCL), fresh leaves were collected from 2-year-old grafted A. senticosus scions, while self-rooted seedling A. senticosus leaves (SCL) were used as controls to analyse the transcriptome and metabolome. Metabolic profiles and gene expression patterns were further identified and correlated in special metabolite target pathways. The content of chlorogenic acid and triterpenoids in the GSCL was higher than in the control, while the quercetin content was lower. All these metabolic changes were associated with changes in the expression pattern of transcripts. Our results revealed the transcriptome and metabolome characteristics of GSCL. This may help to improve leaf quality in A. senticosus cultivation, suggesting that it is feasible to improve the medicinal quality of GSCL through asexual propagation, but the long-term effects need further investigation. In conclusion, this dataset provides a useful resource for future studies on the effects of grafting on medicinal plants.

CRISPR/Cas9-based gene activation and base editing in Populus

The genus Populus has long been used for environmental, agroforestry and industrial applications worldwide. Today Populus is also recognized as a desirable crop for biofuel production and a model tree for physiological and ecological research. As such, various modern biotechnologies, including CRISPR/Cas9-based techniques, have been actively applied to Populus for genetic and genomic improvements for traits such as increased growth rate and tailored lignin composition. However, CRISPR/Cas9 has been primarily used as the active Cas9 form to create knockouts in the hybrid poplar clone "717-1B4" (P. tremula x P. alba clone INRA 717-1B4). Alternative CRISPR/Cas9-based technologies, e.g. those involving modified Cas9 for gene activation and base editing, have not been evaluated in most Populus species for their efficacy. Here we employed a deactivated Cas9 (dCas9)-based CRISPR activation (CRISPRa) technique to fine-tune the expression of two target genes, TPX2 and LecRLK-G which play important roles in plant growth and defense response, in hybrid poplar clone "717-1B4" and poplar clone "WV94" (P. deltoides "WV94"), respectively. We observed that CRISPRa resulted in 1.2-fold to 7.0-fold increase in target gene expression through transient expression in protoplasts and Agrobacterium-mediated stable transformation, demonstrating the effectiveness of dCas9-based CRISPRa system in Populus. In addition, we applied Cas9 nickase (nCas9)-based cytosine base editor (CBE) to precisely introduce premature stop codons via C-to-T conversion, with an efficiency of 13%-14%, in the target gene PLATZ which encodes a transcription factor involved in plant fungal pathogen response in hybrid poplar clone "717-1B4". Overall, we showcase the successful application of CRISPR/Cas-based technologies in gene expression regulation and precise gene engineering in two Populus species, facilitating the adoption of emerging genome editing tools in woody species.

Novel candidate genes for lignin structure identified through genome-wide association study of naturally varying Populus trichocarpa

Populus is a promising lignocellulosic feedstock for biofuels and bioproducts. However, the cell wall biopolymer lignin is a major barrier in conversion of biomass to biofuels. To investigate the variability and underlying genetic basis of the complex structure of lignin, a population of 409 three-year-old, naturally varying Populus trichocarpa genotypes were characterized by heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR). A subsequent genome-wide association study (GWAS) was conducted using approximately 8.3 million single nucleotide polymorphisms (SNPs), which identified 756 genes that were significantly associated (-log₁₀(p-value)>6) with at least one lignin phenotype. Several promising candidate genes were identified, many of which have not previously been reported to be associated with lignin or cell wall biosynthesis. These results provide a resource for gaining insights into the molecular mechanisms of lignin biosynthesis and new targets for future genetic improvement in poplar.

Enhanced genome-wide association reveals the role of YABBY11-NGATHA-LIKE1 in leaf serration development of Populus

Leaf margins are complex plant morphological features that contribute to leaf shape diversity, which affects plant structure, yield, and adaptation. Although several leaf margin regulators have been identified to date, the genetic basis of their natural variation has not been fully elucidated. In this study, we profiled two distinct leaf morphology types (serrated and smooth) using the persistent homology mathematical framework (PHMF) in two poplar species (Populus tomentosa and Populus simonii, respectively). A combined genome-wide association study (GWAS) and expression quantitative trait nucleotide (eQTN) mapping were applied to create a leaf morphology control module using data from P. tomentosa and P. simonii populations. Natural variation in leaf margins was associated with YABBY11 (YAB11) transcript abundance in poplar. In P. tomentosa, PtoYAB11 carries a premature stop codon (PtoYAB11PSC), resulting in the loss of its positive regulation of NGATHA-LIKE1 (PtoNGAL-1) and RIBULOSE BISPHOSPHATE CARBOXYLASE LARGE SUBUNIT (PtoRBCL). Overexpression of PtoYAB11PSC promoted serrated leaf margins, enlarged leaves, enhanced photosynthesis, and increased biomass. Overexpression of PsiYAB11 in P. tomentosa promoted smooth leaf margins, higher stomatal density, and greater light damage repair ability. In poplar, YAB11-NGAL1 is sensitive to environmental conditions, acts as a positive regulator of leaf margin serration, and may also link environmental signaling to leaf morphological plasticity.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

Exploring the role of plant lysin motif receptor-like kinases in regulating plant-microbe interactions in the bioenergy crop Populus

For plants, distinguishing between mutualistic and pathogenic microbes is a matter of survival. All microbes contain microbe-associated molecular patterns (MAMPs) that are perceived by plant pattern recognition receptors (PRRs). Lysin motif receptor-like kinases (LysM-RLKs) are PRRs attuned for binding and triggering a response to specific MAMPs, including chitin oligomers (COs) in fungi, lipo-chitooligosaccharides (LCOs), which are produced by mycorrhizal fungi and nitrogen-fixing rhizobial bacteria, and peptidoglycan in bacteria. The identification and characterization of LysM-RLKs in candidate bioenergy crops including Populus are limited compared to other model plant species, thus inhibiting our ability to both understand and engineer microbe-mediated gains in plant productivity. As such, we performed a sequence analysis of LysM-RLKs in the Populus genome and predicted their function based on phylogenetic analysis with known LysM-RLKs. Then, using predictive models, molecular dynamics simulations, and comparative structural analysis with previously characterized CO and LCO plant receptors, we identified probable ligand-binding sites in Populus LysM-RLKs. Using several machine learning models, we predicted remarkably consistent binding affinity rankings of Populus proteins to CO. In addition, we used a modified Random Walk with Restart network-topology based approach to identify a subset of Populus LysM-RLKs that are functionally related and propose a corresponding signal transduction cascade. Our findings provide the first look into the role of LysM-RLKs in Populus-microbe interactions and establish a crucial jumping-off point for future research efforts to understand specificity and redundancy in microbial perception mechanisms.

Lignin and Its Pathway-Associated Phytoalexins Modulate Plant Defense against Fungi

Fungi infections cause approximately 60-70% yield loss through diseases such as rice blast, powdery mildew, Fusarium rot, downy mildew, etc. Plants naturally respond to these infections by eliciting an array of protective metabolites to confer physical or chemical protection. Among plant metabolites, lignin, a phenolic compound, thickens the middle lamella and the secondary cell walls of plants to curtail fungi infection. The biosynthesis of monolignols (lignin monomers) is regulated by genes whose transcript abundance significantly improves plant defense against fungi. The catalytic activities of lignin biosynthetic enzymes also contribute to the accumulation of other defense compounds. Recent advances focus on modifying the lignin pathway to enhance plant growth and defense against pathogens. This review presents an overview of monolignol regulatory genes and their contributions to fungi immunity, as reported over the last five years. This review expands the frontiers in lignin pathway engineering to enhance plant defense.

Comparative transcriptomics and genomic analyses reveal differential gene expression related to Colletotrichum brevisporum resistance in papaya (Carica papaya L.)

Colletotrichum brevisporum is an important causal pathogen of anthracnose that seriously affects the fruit quality and yield of papaya (Carica papaya L.). Although many genes and biological processes involved in anthracnose resistance have been reported in other species, the molecular mechanisms involved in the response or resistance to anthracnose in post-harvest papaya fruits remain unclear. In this study, we compared transcriptome changes in the post-harvest fruits of the anthracnose-susceptible papaya cultivar Y61 and the anthracnose-resistant cultivar G20 following C. brevisporum inoculation. More differentially expressed genes (DEGs) and differentially expressed long non-coding RNAs (DElnRNAs) were identified in G20 than in Y61, especially at 24 h post-inoculation (hpi), suggesting a prompt activation of defense responses in G20 in the first 24 h after C. brevisporum inoculation. These DEGs were mainly enriched in plant-pathogen interaction, phenylpropanoid biosynthesis/metabolism, and peroxisome and flavonoid biosynthesis pathways in both cultivars. However, in the first 24 hpi, the number of DEGs related to anthracnose resistance was greater in G20 than in Y61, and changes in their expression levels were faster in G20 than in Y61. We also identified a candidate anthracnose-resistant gene cluster, which consisted of 12 genes, 11 in G20 and Y61, in response to C. brevisporum inoculation. Moreover, 529 resistance gene analogs were identified in papaya genome, most of which responded to C. brevisporum inoculation and were genetically different between papaya cultivars and wild-type populations. The total expression dose of the resistance gene analogs may help papaya resist C. brevisporum infection. This study revealed the mechanisms underlying different anthracnose resistance between the anthracnose-resistant and anthracnose-susceptible cultivars based on gene expression, and identified some potential anthracnose resistance-related candidate genes/major regulatory factors. Our findings provided potential targets for developing novel genetic strategies to overcome anthracnose in papaya.

Integrated transcriptomic and transgenic analyses reveal potential mechanisms of poplar resistance to Alternaria alternata infection

BACKGROUND

Populus davidiana × P. bollena is a species of poplar from northeastern China that is characterized by cold resistance and fast growth but now suffers from pathogen infections. Leaf blight caused by Alternaria alternata has become a common poplar disease that causes serious economic impacts, but the molecular mechanisms of resistance to A. alternata in P. davidiana × P. bollena are still unclear.

RESULTS

In this study, the transcriptomic response of P. davidiana × P. bollena to A. alternata infection was determined via RNA-Seq. Twelve cDNA libraries were generated from RNA isolated from three biological replicates at four time points (0, 2, 3, and 4 d post inoculation), and a total of 5,930 differentially expressed genes (DEGs) were detected (| log2 fold change |≥ 1 and FDR values < 0.05). Functional analysis revealed that the DEGs were mainly enriched for the "plant hormone signal transduction" pathway, followed by the "phenylpropanoid biosynthesis" pathway. In addition, DEGs that encode defense-related proteins and are related to ROS metabolism were also identified. Numerous transcription factors, such as the bHLH, WRKY and MYB families, were also induced by A. alternata infection. Among these DEGs, those related to JA biosynthesis and JA signal transduction were consistently activated. Therefore, the lipoxygenase gene PdbLOX2, which is involved in JA biosynthesis, was selected for functional characterization. Overexpression of PdbLOX2 enhanced the resistance of P. davidiana × P. bollena to A. alternata, whereas silencing this gene enhanced susceptibility to A. alternata infection.

CONCLUSIONS

These results provide new insight into the molecular mechanisms of poplar resistance to A. alternata infection and provide candidate genes for breeding resistant cultivars using genetic engineering.

Advanced Breeding for Biotic Stress Resistance in Poplar

Poplar is one of the most important forest trees because of its high economic value. Thanks to the fast-growing rate, easy vegetative propagation and transformation, and availability of genomic resources, poplar has been considered the model species for forest genetics, genomics, and breeding. Being a field-growing tree, poplar is exposed to environmental threats, including biotic stresses that are becoming more intense and diffused because of global warming. Current poplar farming is mainly based on monocultures of a few elite clones and the expensive and long-term conventional breeding programmes of perennial tree species cannot face current climate-change challenges. Consequently, new tools and methods are necessary to reduce the limits of traditional breeding related to the long generation time and to discover new sources of resistance. Recent advances in genomics, marker-assisted selection, genomic prediction, and genome editing offer powerful tools to efficiently exploit the Populus genetic diversity and allow enabling molecular breeding to support accurate early selection, increasing the efficiency, and reducing the time and costs of poplar breeding, that, in turn, will improve our capacity to face or prevent the emergence of new diseases or pests.

Transcriptome and Metabolome Analyses Reveal Differences in Terpenoid and Flavonoid Biosynthesis in Cryptomeria fortunei Needles Across Different Seasons

Cryptomeria fortunei (Chinese cedar) has outstanding medicinal value due to its abundant flavonoid and terpenoid contents. The metabolite contents of C. fortunei needles differ across different seasons. However, the biosynthetic mechanism of these differentially synthesized metabolites (DSMs) is poorly understood. To improve our understanding of this process, we performed integrated non-targeted metabolomic liquid chromatography and gas chromatography mass spectrometry (LC-MS and GC-MS), and transcriptomic analyses of summer and winter needles. In winter, the C. fortunei needle ultrastructure was damaged, and the chlorophyll content and F _v/F _m were significantly (p < 0.05) reduced. Based on GC-MS and LC-MS, we obtained 106 and 413 DSMs, respectively; based on transcriptome analysis, we obtained a total of 41.17 Gb of clean data and assembled 33,063 unigenes, including 14,057 differentially expressed unigenes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DSMs/DEGs were significantly (p < 0.05) enriched in many biosynthesis pathways, such as terpenoids, photosynthates, and flavonoids. Integrated transcriptomic and metabonomic analyses showed that seasonal changes have the greatest impact on photosynthesis pathways, followed by terpenoid and flavonoid biosynthesis pathways. In summer Chinese cedar (SCC) needles, DXS, DXR, and ispH in the 2-methyl-pentaerythritol 4-phosphate (MEP) pathway and GGPS were highly expressed and promoted the accumulation of terpenoids, especially diterpenoids. In winter Chinese cedar (WCC) needles, 9 genes (HCT, CHS, CHI, F3H, F3'H, F3'5'H, FLS, DFR, and LAR) involved in flavonoid biosynthesis were highly expressed and promoted flavonoid accumulation. This study broadens our understanding of the metabolic and transcriptomic changes in C. fortunei needles caused by seasonal changes and provides a reference regarding the adaptive mechanisms of C. fortunei and the extraction of its metabolites.

Integrative analysis of the metabolome and transcriptome reveals the molecular mechanism of chlorogenic acid synthesis in peach fruit

As the most abundant phenolic acid in peach fruit, chlorogenic acid (CGA) is an important entry point for the development of natural dietary supplements and functional foods. However, the metabolic and regulation mechanisms underlying its accumulation in peach fruits remain unclear. In this study, we evaluated the composition and content of CGAs in mature fruits of 205 peach cultivars. In peach fruits, three forms of CGA (52.57%), neochlorogenic acid (NCGA, 47.13%), and cryptochlorogenic acid (CCGA, 0.30%) were identified. During the growth and development of peach fruits, the content of CGAs generally showed a trend of rising first and then decreasing. Notably, the contents of quinic acid, shikimic acid, p-coumaroyl quinic acid, and caffeoyl shikimic acid all showed similar dynamic patterns to that of CGA, which might provide the precursor material basis for the accumulation of CGA in the later stage. Moreover, CGA, lignin, and anthocyanins might have a certain correlation and these compounds work together to maintain a dynamic balance. By the comparative transcriptome analysis, 8 structural genes (Pp4CL, PpCYP98A, and PpHCT) and 15 regulatory genes (PpMYB, PpWRKY, PpERF, PpbHLH, and PpWD40) were initially screened as candidate genes of CGA biosynthesis. Our findings preliminarily analyzed the metabolic and molecular regulation mechanisms of CGA biosynthesis in peach fruit, which provided a theoretical basis for developing high-CGA content peaches in future breeding programs.

Defining key metabolic roles in osmotic adjustment and ROS homeostasis in the recretohalophyte Karelinia caspia under salt stress

The recretohalophyte Karelinia caspia is of forage and medical value and can remediate saline soils. We here assess the contribution of primary/secondary metabolism to osmotic adjustment and ROS homeostasis in Karelinia caspia under salt stress using multi-omic approaches. Computerized phenomic assessments, tests for cellular osmotic changes and lipid peroxidation indicated that salt treatment had no detectable physical effect on K. caspia. Metabolomic analysis indicated that amino acids, saccharides, organic acids, polyamine, phenolic acids, and vitamins accumulated significantly with salt treatment. Transcriptomic assessment identified differentially expressed genes closely linked to the changes in above primary/secondary metabolites under salt stress. In particular, shifts in carbohydrate metabolism (TCA cycle, starch and sucrose metabolism, glycolysis) as well as arginine and proline metabolism were observed to maintain a low osmotic potential. Chlorogenic acid/vitamin E biosynthesis was also enhanced, which would aid in ROS scavenging in the response of K. caspia to salt. Overall, our findings define key changes in primary/secondary metabolism that are coordinated to modulate the osmotic balance and ROS homeostasis to contribute to the salt tolerance of K. caspia.

Comparative transcriptome and weighted correlation network analyses reveal candidate genes involved in chlorogenic acid biosynthesis in sweet potato

Chlorogenic acids (CGAs) are important secondary metabolites produced in sweet potato. However, the mechanisms of their biosynthesis and regulation remain unclear. To identify potential genes involved in CGA biosynthesis, analysis of the dynamic changes in CGA components and RNA sequencing were performed on young leaves (YL), mature leaves (ML), young stems (YS), mature stems (MS) and storage roots (SR). Accordingly, we found that the accumulation of six CGA components varied among the different tissues and developmental stages, with YS and YL recording the highest levels, while SR exhibited low levels. Moreover, the transcriptome analysis yielded 59,287 unigenes, 3,767 of which were related to secondary-metabolite pathways. The differentially expressed genes (DEGs) were identified based on CGA content levels by comparing the different samples, including ML vs. YL, MS vs. YS, SR vs. YL and SR vs. YS. A total of 501 common DEGs were identified, and these were mainly implicated in the secondary metabolites biosynthesis. Additionally, eight co-expressed gene modules were identified following weighted gene co-expression network analysis, while genes in darkgrey module were highly associated with CGA accumulation. Darkgrey module analysis revealed that 12 unigenes encoding crucial enzymes (PAL, 4CL, C4H, C3H and HCT/HQT) and 42 unigenes encoding transcription factors (MYB, bHLH, WD40, WRKY, ERF, MADS, GARS, bZIP and zinc finger protein) had similar expression patterns with change trends of CGAs, suggesting their potential roles in CGA metabolism. Our findings provide new insights into the biosynthesis and regulatory mechanisms of CGA pathway, and will inform future efforts to build a genetically improve sweet potato through the breeding of high CGA content varieties.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

A systems genetics approach reveals PbrNSC as a regulator of lignin and cellulose biosynthesis in stone cells of pear fruit

BACKGROUND

Stone cells in fruits of pear (Pyrus pyrifolia) negatively influence fruit quality because their lignified cell walls impart a coarse and granular texture to the fruit flesh.

RESULTS

We generate RNA-seq data from the developing fruits of 206 pear cultivars with a wide range of stone cell contents and use a systems genetics approach to integrate co-expression networks and expression quantitative trait loci (eQTLs) to characterize the regulatory mechanisms controlling lignocellulose formation in the stone cells of pear fruits. Our data with a total of 35,897 expressed genes and 974,404 SNPs support the identification of seven stone cell formation modules and the detection of 139,515 eQTLs for 3229 genes in these modules. Focusing on regulatory factors and using a co-expression network comprising 39 structural genes, we identify PbrNSC as a candidate regulator of stone cell formation. We then verify the function of PbrNSC in regulating lignocellulose formation using both pear fruit and Arabidopsis plants and further show that PbrNSC can transcriptionally activate multiple target genes involved in secondary cell wall formation.

CONCLUSIONS

This study generates a large resource for studying stone cell formation and provides insights into gene regulatory networks controlling the formation of stone cell and lignocellulose.

Genetic architecture of the metabolic pathway of salicylic acid biosynthesis in Populus

Salicylic acid (SA) is a vital hormone for adaptive responses to biotic and abiotic stresses, which facilitates growth-immunity trade-offs in plants. However, the genetic regulatory networks underlying the metabolic pathway of SA biosynthesis in perennial species remain unclear. Here, we integrated genome-wide association study (GWAS) with metabolite and expression profiling methodologies to dissect the genetic architecture of SA biosynthesis in Populus. First, we quantified nine intermediate metabolites of SA biosynthesis in 300 unrelated Populus tomentosa Carr. individuals. Then, we used a systematic genetic strategy to identify candidate genes for constructing the genetic regulatory network of SA biosynthesis. We focused on WRKY70, an efficient transcription factor, as the key causal gene in the regulatory network, and combined the novel genes coordinating the accumulation of SA. Finally, we identified eight GWAS signals and eight expression quantitative trait loci situated in a selective sweep, and showed the presence of large allele frequency differences among the three geographic populations, revealing that candidate genes subject to selection were involved in SA biosynthesis. This study provides an integrated strategy for dissecting the genetic architecture of the metabolic pathway of SA biosynthesis in Populus, thereby enhancing our understanding of genetic regulation of SA biosynthesis in trees, and accelerating marker-assisted breeding efforts toward high-resistance elite varieties of Populus.

Metabolomic Patterns of Septoria Canker Resistant and Susceptible Populus trichocarpa Genotypes 24 Hours Postinoculation

Sphaerulina musiva is an economically and ecologically important fungal pathogen that causes Septoria stem canker and leaf spot disease of Populus species. To bridge the gap between genetic markers and structural barriers previously found to be linked to Septoria canker disease resistance in poplar, we used hydrophilic interaction liquid chromatography and tandem mass spectrometry to identify and quantify metabolites involved with signaling and cell wall remodeling. Fluctuations in signaling molecules, organic acids, amino acids, sterols, phenolics, and saccharides in resistant and susceptible P. trichocarpa inoculated with S. musiva were observed. The patterns of 222 metabolites in the resistant host implicate systemic acquired resistance (SAR), cell wall apposition, and lignin deposition as modes of resistance to this hemibiotrophic pathogen. This pattern is consistent with the expected response to the biotrophic phase of S. musiva colonization during the first 24 h postinoculation. The fungal pathogen metabolized key regulatory signals of SAR, other phenolics, and precursors of lignin biosynthesis that were depleted in the susceptible host. This is the first study to characterize metabolites associated with the response to initial colonization by S. musiva between resistant and susceptible hosts.

The Ancient Salicoid Genome Duplication Event: A Platform for Reconstruction of De Novo Gene Evolution in Populus trichocarpa

Orphan genes are characteristic genomic features that have no detectable homology to genes in any other species and represent an important attribute of genome evolution as sources of novel genetic functions. Here, we identified 445 genes specific to Populus trichocarpa. Of these, we performed deeper reconstruction of 13 orphan genes to provide evidence of de novo gene evolution. Populus and its sister genera Salix are particularly well suited for the study of orphan gene evolution because of the Salicoid whole-genome duplication event which resulted in highly syntenic sister chromosomal segments across the Salicaceae. We leveraged this genomic feature to reconstruct de novo gene evolution from intergenera, interspecies, and intragenomic perspectives by comparing the syntenic regions within the P. trichocarpa reference, then P. deltoides, and finally Salix purpurea. Furthermore, we demonstrated that 86.5% of the putative orphan genes had evidence of transcription. Additionally, we also utilized the Populus genome-wide association mapping panel, a collection of 1,084 undomesticated P. trichocarpa genotypes to further determine putative regulatory networks of orphan genes using expression quantitative trait loci (eQTL) mapping. Functional enrichment of these eQTL subnetworks identified common biological themes associated with orphan genes such as response to stress and defense response. We also identify a putative cis-element for a de novo gene and leverage conserved synteny to describe evolution of a putative transcription factor binding site. Overall, 45% of orphan genes were captured in trans-eQTL networks.

PopulusPtERF85 Balances Xylem Cell Expansion and Secondary Cell Wall Formation in Hybrid Aspen

Secondary growth relies on precise and specialized transcriptional networks that determine cell division, differentiation, and maturation of xylem cells. We identified a novel role for the ethylene-induced Populus Ethylene Response Factor PtERF85 (Potri.015G023200) in balancing xylem cell expansion and secondary cell wall (SCW) formation in hybrid aspen (Populus tremula x tremuloides). Expression of PtERF85 is high in phloem and cambium cells and during the expansion of xylem cells, while it is low in maturing xylem tissue. Extending PtERF85 expression into SCW forming zones of woody tissues through ectopic expression reduced wood density and SCW thickness of xylem fibers but increased fiber diameter. Xylem transcriptomes from the transgenic trees revealed transcriptional induction of genes involved in cell expansion, translation, and growth. The expression of genes associated with plant vascular development and the biosynthesis of SCW chemical components such as xylan and lignin, was down-regulated in the transgenic trees. Our results suggest that PtERF85 activates genes related to xylem cell expansion, while preventing transcriptional activation of genes related to SCW formation. The importance of precise spatial expression of PtERF85 during wood development together with the observed phenotypes in response to ectopic PtERF85 expression suggests that PtERF85 contributes to the transition of fiber cells from elongation to secondary cell wall deposition.

Integration of genome wide association studies and co-expression networks reveal roles of PtoWRKY 42-PtoUGT76C1-1 in trans-zeatin metabolism and cytokinin sensitivity in poplar

Cytokinins are important for in vitro shoot regeneration in plants. Cytokinin N-glucosides are produced via an irreversible glycosylation pathway, which regulates the endogenous cytokinin content. Although cytokinin N-glucoside pathways have been uncovered in higher plants, no regulator has been identified to date. We performed a metabolome genome-wide association study (mGWAS), weighted gene co-expression network analysis (WGCNA), and expression quantitative trait nucleotide (eQTN) mappings to build a core triple genetic network (mGWAS-gene expression-phenotype) for the trans-zeatin N-glucoside (ZNG) metabolite using data from 435 unrelated Populus tomentosa individuals. Variation of the ZNG level in poplar is attributed to the differential transcription of PtoWRKY42, a member of WRKY multigene family group IIb. Functional analysis revealed that PtoWRKY42 negatively regulated ZNG accumulation by binding directly to the W-box of the UDP-glycosyltransferase 76C 1-1 (PtoUGT761-1) promoter. Also, PtoWRKY42 was strongly induced by leaf senescence, 6-BA, wounding, and salt stress, resulting in a reduced ZNG level. We identified PtoWRKY42, a negative regulator of cytokinin N-glucosides, which contributes to the natural variation in ZNG level and mediates ZNG accumulation by directly modulating the key glycosyltransferase gene PtoUGT76C1-1.

Transcriptional Regulation and Signaling of Developmental Programmed Cell Death in Plants

Developmental programmed cell death (dPCD) has multiple functions in plant growth and development, and is of great value for industrial production. Among them, wood formed by xylem dPCD is one of the most widely used natural materials. Therefore, it is crucial to explore the molecular mechanism of plant dPCD. The dPCD process is tightly regulated by genetic networks and is involved in the transduction of signaling molecules. Several key regulators have been identified in diverse organisms and individual PCD events. However, complex molecular networks controlling plant dPCD remain highly elusive, and the original triggers of this process are still unknown. This review summarizes the recent progress on the transcriptional regulation and signaling of dPCD during vegetative and reproductive development. It is hoped that this review will provide an overall view of the molecular regulation of dPCD in different developmental processes in plants and identify specific mechanisms for regulating these dPCD events. In addition, the application of plants in industrial production can be improved by manipulating dPCD in specific processes, such as xylogenesis.

CRISPR-Based Genome Editing Tools: Insights into Technological Breakthroughs and Future Challenges

Genome-editing (GE) is having a tremendous influence around the globe in the life science community. Among its versatile uses, the desired modifications of genes, and more importantly the transgene (DNA)-free approach to develop genetically modified organism (GMO), are of special interest. The recent and rapid developments in genome-editing technology have given rise to hopes to achieve global food security in a sustainable manner. We here discuss recent developments in CRISPR-based genome-editing tools for crop improvement concerning adaptation, opportunities, and challenges. Some of the notable advances highlighted here include the development of transgene (DNA)-free genome plants, the availability of compatible nucleases, and the development of safe and effective CRISPR delivery vehicles for plant genome editing, multi-gene targeting and complex genome editing, base editing and prime editing to achieve more complex genetic engineering. Additionally, new avenues that facilitate fine-tuning plant gene regulation have also been addressed. In spite of the tremendous potential of CRISPR and other gene editing tools, major challenges remain. Some of the challenges are related to the practical advances required for the efficient delivery of CRISPR reagents and for precision genome editing, while others come from government policies and public acceptance. This review will therefore be helpful to gain insights into technological advances, its applications, and future challenges for crop improvement.

Characterization and functional analysis of the Hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyl transferase (HCT) gene family in poplar

Hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyl transferase (HCT) divides the mass flux to H, G and S units in monolignol biosynthesis and affects lignin content. Ten HCT homologs were identified in the Populus trichocarpa (Torr. & Gray) genome. Both genome duplication and tandem duplication resulted in the expansion of HCT orthologs in Populus. Comprehensive analysis including motif analysis, phylogenetic analysis, expression profiles and co-expression analysis revealed the divergence and putative function of these candidate PoptrHCTs. PoptrHCT1 and 2 were identified as likely involved in lignin biosynthesis. PoptrHCT9 and 10- are likely to be involved in plant development and the response to cold stress. Similar functional divergence was also identified in Populus tomentosa Carr. Enzymatic assay of PtoHCT1 showed that PtoHCT1 was able to synthesize caffeoyl shikimate using caffeoyl-CoA and shikimic acid as substrates.

Genetic Architecture Underlying the Metabolites of Chlorogenic Acid Biosynthesis in Populus tomentosa

Chlorogenic acid (CGA) plays a crucial role in defense response, immune regulation, and the response to abiotic stress in plants. However, the genetic regulatory network of CGA biosynthesis pathways in perennial plants remains unclear. Here, we investigated the genetic architecture for CGA biosynthesis using a metabolite-based genome-wide association study (mGWAS) and expression quantitative trait nucleotide (eQTN) mapping in a population of 300 accessions of Populus tomentosa. In total, we investigated 204 SNPs which were significantly associated with 11 metabolic traits, corresponding to 206 genes, and were mainly involved in metabolism and cell growth processes of P. tomentosa. We identified 874 eQTNs representing 1066 genes, in which the expression and interaction of causal genes affected phenotypic variation. Of these, 102 genes showed significant signatures of selection in three geographical populations, which provided insights into the adaptation of CGA biosynthesis to the local environment. Finally, we constructed a genetic network of six causal genes that coordinately regulate CGA biosynthesis, revealing the multiple regulatory patterns affecting CGA accumulation in P. tomentosa. Our study provides a multiomics strategy for understanding the genetic basis underlying the natural variation in the CGA biosynthetic metabolites of Populus, which will enhance the genetic development of abiotic-resistance varieties in forest trees.

Transcriptome analysis and association mapping reveal the genetic regulatory network response to cadmium stress in Populus tomentosa

Long non-coding RNAs (lncRNAs) play essential roles in plant abiotic stress responses, but the response of lncRNA-mediated genetic networks to cadmium (Cd) treatment remain elusive in trees, the promising candidates for phytoremediation of Cd contamination. We identified 172 Cd-responsive lncRNAs and 295 differentially expressed target genes in the leaves of Cd-treated Populus tomentosa. Functional annotation revealed that these lncRNAs were involved in various processes, including photosynthesis, hormone regulation, and phenylalanine metabolism. Association studies identified 78 significant associations, representing 14 Cd-responsive lncRNAs and 28 target genes for photosynthetic and leaf physiological traits. Epistasis uncovered 83 pairwise interactions among these traits, revealing Cd-responsive lncRNA-mediated genetic networks for photosynthesis and leaf physiology in P. tomentosa. We focused on the roles of two Cd-responsive lncRNA-gene pairs, MSTRG.22608.1-PtoMYB73 and MSTRG.5634.1-PtoMYB27, in Cd tolerance of Populus, and detected insertions/deletions within lncRNAs as polymorphisms driving target gene expression. Genotype analysis of lncRNAs and heterologous overexpression of PtoMYB73 and PtoMYB27 in Arabidopsis indicated their effects on enhancing Cd tolerance, photosynthetic rate, and leaf growth, and the potential interaction mechanisms of PtoMYB73 with abiotic stresses. Our study identifies the genetic basis for the response of Populus to Cd treatment, facilitating genetic improvement of Cd tolerance in trees.

Integrated transcriptome, small RNA and degradome sequencing approaches proffer insights into chlorogenic acid biosynthesis in leafy sweet potato

Leafy sweet potato is rich in total phenolics (TP) which play key roles in health protection, the chlorogenic acid (CGA) constitutes the major components of phenolic compounds in leafy sweet potato. Unfortunately, the mechanism of CGA biosynthesis in leafy sweet potato is unclear. To dissect the mechanisms of CGA biosynthesis, we performed transcriptome, small RNA (sRNA) and degradome sequencing of one low-CGA content and one high-CGA content genotype at two stages. A total of 2,333 common differentially expressed genes (DEGs) were identified, and the enriched DEGs were related to photosynthesis, starch and sucrose metabolism and phenylpropanoid biosynthesis. The functional genes, such as CCR, CCoAOMT and HCT in the CGA biosynthetic pathway were down-regulated, indicating that the way to lignin was altered, and two possible CGA biosynthetic routes were hypothesized. A total of 38 DE miRNAs were identified, and 1,799 targets were predicated for 38 DE miRNAs by using in silico approaches. The target genes were enriched in lignin and phenylpropanoid catabolic processes. Transcription factors (TFs) such as apetala2/ethylene response factor (AP2/ERF) and Squamosa promoter binding protein-like (SPL) predicated in silico were validated by degradome sequencing. Association analysis of the DE miRNAs and transcriptome datasets identified that miR156 family negatively targeted AP2/ERF and SPL. Six mRNAs and six miRNAs were validated by qRT-PCR, and the results showed that the expression levels of the mRNAs and miRNAs were consistent with the sequencing data. This study established comprehensive functional genomic resources for the CGA biosynthesis, and provided insights into the molecular mechanisms involving in this process. The results also enabled the first perceptions of the regulatory roles of mRNAs and miRNAs, and offered candidate genes for leafy sweet potato improvements.

Transcriptional Regulation of Drought Response in Arabidopsis and Woody Plants

Within the context of global warming, long-living plants such as perennial woody species endure adverse conditions. Among all of the abiotic stresses, drought stress is one of the most detrimental stresses that inhibit plant growth and productivity. Plants have evolved multiple mechanisms to respond to drought stress, among which transcriptional regulation is one of the key mechanisms. In this review, we summarize recent progress on the regulation of drought response by transcription factor (TF) families, which include abscisic acid (ABA)-dependent ABA-responsive element/ABRE-binding factors (ABRE/ABF), WRKY, and Nuclear Factor Y families, as well as ABA-independent AP2/ERF and NAC families, in the model plant Arabidopsis. We also review what is known in woody species, particularly Populus, due to its importance and relevance in economic and ecological processes. We discuss opportunities for a deeper understanding of drought response in woody plants with the development of high-throughput omics analyses and advanced genome editing techniques.

Phylogenetic diversity of 200+ isolates of the ectomycorrhizal fungus Cenococcum geophilum associated with Populus trichocarpa soils in the Pacific Northwest, USA and comparison to globally distributed representatives

The ectomycorrhizal fungal symbiont Cenococcum geophilum is of high interest as it is globally distributed, associates with many plant species, and has resistance to multiple environmental stressors. C. geophilum is only known from asexual states but is often considered a cryptic species complex, since extreme phylogenetic divergence is often observed within nearly morphologically identical strains. Alternatively, C. geophilum may represent a highly diverse single species, which would suggest cryptic but frequent recombination. Here we describe a new isolate collection of 229 C. geophilum isolates from soils under Populus trichocarpa at 123 collection sites spanning a ~283 mile north-south transect in Western Washington and Oregon, USA (PNW). To further understanding of the phylogenetic relationships within C. geophilum, we performed maximum likelihood and Bayesian phylogenetic analyses to assess divergence within the PNW isolate collection, as well as a global phylogenetic analysis of 789 isolates with publicly available data from the United States, Japan, and European countries. Phylogenetic analyses of the PNW isolates revealed three distinct phylogenetic groups, with 15 clades that strongly resolved at >80% bootstrap support based on a GAPDH phylogeny and one clade segregating strongly in two principle component analyses. The abundance and representation of PNW isolate clades varied greatly across the North-South range, including a monophyletic group of isolates that spanned nearly the entire gradient at ~250 miles. A direct comparison between the GAPDH and ITS rRNA gene region phylogenies, combined with additional analyses revealed stark incongruence between the ITS and GAPDH gene regions, consistent with intra-species recombination between PNW isolates. In the global isolate collection phylogeny, 34 clades were strongly resolved using Maximum Likelihood and Bayesian approaches (at >80% MLBS and >0.90 BPP respectively), with some clades having intra- and intercontinental distributions. Together these data are highly suggestive of divergence within multiple cryptic species, however additional analyses such as higher resolution genotype-by-sequencing approaches are needed to distinguish potential species boundaries and the mode and tempo of recombination patterns.

Recent Advances in the Roles of HSFs and HSPs in Heat Stress Response in Woody Plants

A continuous increase in ambient temperature caused by global warming has been considered a worldwide threat. As sessile organisms, plants have evolved sophisticated heat shock response (HSR) to respond to elevated temperatures and other abiotic stresses, thereby minimizing damage and ensuring the protection of cellular homeostasis. In particular, for perennial trees, HSR is crucial for their long life cycle and development. HSR is a cell stress response that increases the number of chaperones including heat shock proteins (HSPs) to counter the negative effects on proteins caused by heat and other stresses. There are a large number of HSPs in plants, and their expression is directly regulated by a series of heat shock transcription factors (HSFs). Therefore, understanding the detailed molecular mechanisms of woody plants in response to extreme temperature is critical for exploring how woody species will be affected by climate changes. In this review article, we summarize the latest findings of the role of HSFs and HSPs in the HSR of woody species and discuss their regulatory networks and cross talk in HSR. In addition, strategies and programs for future research studies on the functions of HSFs and HSPs in the HSR of woody species are also proposed.