The Plant Genome Integrative Explorer Resource: PlantGenIE.org

The TELOMERE REPEAT BINDING proteins TRB4 and TRB5 function as transcriptional activators of PRC2-controlled genes to regulate plant development

Plant-specific transcriptional regulators called TELOMERE REPEAT BINDING proteins (TRBs) combine two DNA-binding domains, the GH1 domain, which binds to linker DNA and is shared with H1 histones, and the Myb/SANT domain, which specifically recognizes the telobox DNA-binding site motif. TRB1, TRB2, and TRB3 proteins recruit Polycomb group complex 2 (PRC2) to deposit H3K27me3 and JMJ14 to remove H3K4me3 at gene promoters containing telobox motifs to repress transcription. Here, we demonstrate that TRB4 and TRB5, two related paralogs belonging to a separate TRB clade conserved in spermatophytes, regulate the transcription of several hundred genes involved in developmental responses to environmental cues. TRB4 binds to several thousand sites in the genome, mainly at transcription start sites and promoter regions of transcriptionally active and H3K4me3-marked genes, but, unlike TRB1, it is not enriched at H3K27me3-marked gene bodies. However, TRB4 can physically interact with the catalytic components of PRC2, SWINGER, and CURLY LEAF (CLF). Unexpectedly, we show that TRB4 and TRB5 are required for distinctive phenotypic traits observed in clf mutant plants and thus function as transcriptional activators of several hundred CLF-controlled genes, including key flowering genes. We further demonstrate that TRB4 shares multiple target genes with TRB1 and physically and genetically interacts with members of both TRB clades. Collectively, these results reveal that TRB proteins engage in both positive and negative interactions with other members of the family to regulate plant development through both PRC2-dependent and -independent mechanisms.

Transcriptomics, proteomics, and metabolomics interventions prompt crop improvement against metal(loid) toxicity

The escalating challenges posed by metal(loid) toxicity in agricultural ecosystems, exacerbated by rapid climate change and anthropogenic pressures, demand urgent attention. Soil contamination is a critical issue because it significantly impacts crop productivity. The widespread threat of metal(loid) toxicity can jeopardize global food security due to contaminated food supplies and pose environmental risks, contributing to soil and water pollution and thus impacting the whole ecosystem. In this context, plants have evolved complex mechanisms to combat metal(loid) stress. Amid the array of innovative approaches, omics, notably transcriptomics, proteomics, and metabolomics, have emerged as transformative tools, shedding light on the genes, proteins, and key metabolites involved in metal(loid) stress responses and tolerance mechanisms. These identified candidates hold promise for developing high-yielding crops with desirable agronomic traits. Computational biology tools like bioinformatics, biological databases, and analytical pipelines support these omics approaches by harnessing diverse information and facilitating the mapping of genotype-to-phenotype relationships under stress conditions. This review explores: (1) the multifaceted strategies that plants use to adapt to metal(loid) toxicity in their environment; (2) the latest findings in metal(loid)-mediated transcriptomics, proteomics, and metabolomics studies across various plant species; (3) the integration of omics data with artificial intelligence and high-throughput phenotyping; (4) the latest bioinformatics databases, tools and pipelines for single and/or multi-omics data integration; (5) the latest insights into stress adaptations and tolerance mechanisms for future outlooks; and (6) the capacity of omics advances for creating sustainable and resilient crop plants that can thrive in metal(loid)-contaminated environments.

A cytokinin response factor PtCRF1 is involved in the regulation of wood formation in poplar

Wood formation is a complex developmental process under the control of multiple levels of regulatory transcriptional network and hormone signals in trees. It is well known that cytokinin (CK) signaling plays an important role in maintaining the activity of the vascular cambium. The CK response factors (CRFs) encoding a subgroup of AP2 transcription factors have been identified to mediate the CK-dependent regulation in different plant developmental processes. However, the functions of CRFs in wood development remain unclear. Here, we characterized the function of PtCRF1, a CRF transcription factor isolated from poplar, in the process of wood formation. The PtCRF1 is preferentially expressed in secondary vasculature, especially in vascular cambium and secondary phloem, and encodes a transcriptional activator. Overexpression of PtCRF1 in transgenic poplar plants led to a significant reduction in the cell layer number of vascular cambium. The development of wood tissue was largely promoted in the PtCRF1-overexpressing lines, while it was significantly compromised in the CRISPR/Cas9-generated double mutant plants of PtCRF1 and its closest homolog PtCRF2. The RNA sequencing (RNA-seq) and quantitative reverse transcription PCR (RT-qPCR) analyses showed that PtCRF1 repressed the expression of the typical CK-responsive genes. Furthermore, bimolecular fluorescence complementation assays revealed that PtCRF1 competitively inhibits the direct interactions between histidine phosphotransfer proteins and type-B response regulator by binding to PtHP protein. Collectively, these results indicate that PtCRF1 negatively regulates CK signaling and is required for woody cell differentiation in poplar.

Reference genes for Eucalyptus spp. under Beauveria bassiana inoculation and subsequently infestation by the galling wasp Leptocybe invasa

Relative gene expression analysis through RT-qPCR is an important molecular technique that helps understanding different molecular mechanisms, such as the plant defense response to insect pests. However, the use of RT-qPCR for gene expression analysis can be affected by factors that directly affect the reliability of the results. Among these factors, the appropriate choice of reference genes is crucial and can strongly impact RT-qPCR relative gene expression analyses, highlighting the importance in correctly choosing the most suitable genes for the success of the analysis. Thus, this study aimed to select and validate reference genes for relative gene expression studies through RT-qPCR in hybrids of Eucalyptus tereticornis × Eucalyptus camaldulensis (drought tolerant and susceptible to Leptocybe invasa) under conditions of inoculation by the Beauveria bassiana fungus and subsequent infestation by L. invasa. The expression level and stability of eleven candidate genes were evaluated. Stability was analyzed using the RefFinder tool, which integrates the geNorm, NormFinder, BestKeeper, and Delta-Ct algorithms. The selected reference genes were validated through the expression analysis of the transcriptional factor EcDREB2 (dehydration-responsive element-binding protein 2). For all treatments evaluated, EcPTB, EcPP2A-1, and EcEUC12 were the best reference genes. The triplets EcPTB/EcEUC12/EcUBP6, EcPP2A-1/EcEUC12/EcPTB, EcIDH/EcSAND/Ecα-TUB, EcPP2A-1/Ecα-TUB/EcPTB, and EcPP2A-1/EcUPL7/EcSAND were the best reference genes for the control plants, mother plants, plants inoculated with B. bassiana, plants infested with L. invasa, and plants inoculated with B. bassiana and subsequently infested with L. invasa, respectively. The best determined reference genes were used to normalize the RT-qPCR expression data for each experimental condition evaluated. The results emphasize the importance of this type of study to ensure the reliability of relative gene expression analyses. Furthermore, the findings of this study can be used as a basis for future research, comprising gene expression analysis of different eucalyptus metabolic pathways.

Genome-wide molecular evolution analysis of the GRF and GIF gene families in Plantae (Archaeplastida)

BACKGROUND

Plant growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) interact with each other and collectively have important regulatory roles in plant growth, development, and stress responses. Therefore, it is of great significance to explore the systematic evolution of GRF and GIF gene families. However, our knowledge and understanding of the role of GRF and GIF genes during plant evolution has been fragmentary.

RESULTS

In this study, a large number of genomic and transcriptomic datasets of algae, mosses, ferns, gymnosperms and angiosperms were used to systematically analyze the evolution of GRF and GIF genes during the evolution of plants. The results showed that GRF gene first appeared in the charophyte Klebsormidium nitens, whereas the GIF genes originated relatively early, and these two gene families were mainly expanded by segmental duplication events after plant terrestrialization. During the process of evolution, the protein sequences and functions of GRF and GIF family genes are relatively conservative. As cooperative partner, GRF and GIF genes contain the similar types of cis-acting elements in their promoter regions, which enables them to have similar transcriptional response patterns, and both show higher levels of expression in reproductive organs and tissues and organs with strong capacity for cell division. Based on protein-protein interaction analysis and verification, we found that the GRF-GIF protein partnership began to be established in pteridophytes and is highly conserved across different terrestrial plants.

CONCLUSIONS

These results provide a foundation for further exploration of the molecular evolution and biological functions of GRF and GIF genes.

PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar

Introduction

The biosynthesis of the secondary cell wall (SCW) is orchestrated by an intricate hierarchical transcriptional regulatory network. This network is initiated by first-layer master switches, SCW-NAC transcription factors, which in turn activate the second-layer master switches MYBs. These switches play a crucial role in regulating xylem specification and differentiation during SCW formation. However, the roles of most MYBs in woody plants are yet to be fully understood.

Methods

In this study, we identified and isolated the R2R3-MYB transcription factor, PtoMYB031, from Populus tomentosa. We explored its expression, mainly in xylem tissues, and its role as a transcriptional repressor in the nucleus. We used overexpression and RNA interference techniques in poplar, along with Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, to analyze the regulatory effects of PtoMYB031.

Results

Overexpression of PtoMYB031 in poplar significantly reduced lignin, cellulose, and hemicellulose content, and inhibited vascular development in stems, resulting in decreased SCW thickness in xylem tissues. Gene expression analysis showed that structural genes involved in SCW biosynthesis were downregulated in PtoMYB031-OE lines. Conversely, RNA interference of PtoMYB031 increased these compounds. Additionally, PtoMYB031 was found to recruit the repressor PtoZAT11, forming a transcriptional inhibition complex.

Discussion

Our findings provide new insights into how PtoMYB031, through its interaction with PtoZAT11, forms a complex that can suppress the expression of key regulatory genes, PtoWND1A and PtoWND2B, in SCW biosynthesis. This study enhances our understanding of the transcriptional regulation involved in SCW formation in poplar, highlighting the significant role of PtoMYB031.

PPGR: a comprehensive perennial plant genomes and regulation database

Perennial woody plants hold vital ecological significance, distinguished by their unique traits. While significant progress has been made in their genomic and functional studies, a major challenge persists: the absence of a comprehensive reference platform for collection, integration and in-depth analysis of the vast amount of data. Here, we present PPGR (Resource for Perennial Plant Genomes and Regulation; https://ngdc.cncb.ac.cn/ppgr/) to address this critical gap, by collecting, integrating, analyzing and visualizing genomic, gene regulation and functional data of perennial plants. PPGR currently includes 60 species, 847 million protein-protein/TF (transcription factor)-target interactions, 9016 transcriptome samples under various environmental conditions and genetic backgrounds. Noteworthy is the focus on genes that regulate wood production, seasonal dormancy, terpene biosynthesis and leaf senescence representing a wealth of information derived from experimental data, literature mining, public databases and genomic predictions. Furthermore, PPGR incorporates a range of multi-omics search and analysis tools to facilitate browsing and application of these extensive datasets. PPGR represents a comprehensive and high-quality resource for perennial plants, substantiated by an illustrative case study that demonstrates its capacity in unraveling gene functions and shedding light on potential regulatory processes.

Flexure wood formation via growth reprogramming in hybrid aspen involves jasmonates and polyamines and transcriptional changes resembling tension wood development

Stem bending in trees induces flexure wood but its properties and development are poorly understood. Here, we investigated the effects of low-intensity multidirectional stem flexing on growth and wood properties of hybrid aspen, and on its transcriptomic and hormonal responses. Glasshouse-grown trees were either kept stationary or subjected to several daily shakes for 5 wk, after which the transcriptomes and hormones were analyzed in the cambial region and developing wood tissues, and the wood properties were analyzed by physical, chemical and microscopy techniques. Shaking increased primary and secondary growth and altered wood differentiation by stimulating gelatinous-fiber formation, reducing secondary wall thickness, changing matrix polysaccharides and increasing cellulose, G- and H-lignin contents, cell wall porosity and saccharification yields. Wood-forming tissues exhibited elevated jasmonate, polyamine, ethylene and brassinosteroids and reduced abscisic acid and gibberellin signaling. Transcriptional responses resembled those during tension wood formation but not opposite wood formation and revealed several thigmomorphogenesis-related genes as well as novel gene networks including FLA and XTH genes encoding plasma membrane-bound proteins. Low-intensity stem flexing stimulates growth and induces wood having improved biorefinery properties through molecular and hormonal pathways similar to thigmomorphogenesis in herbaceous plants and largely overlapping with the tension wood program of hardwoods.

Integrating omics databases for enhanced crop breeding

Crop plant breeding involves selecting and developing new plant varieties with desirable traits such as increased yield, improved disease resistance, and enhanced nutritional value. With the development of high-throughput technologies, such as genomics, transcriptomics, and metabolomics, crop breeding has entered a new era. However, to effectively use these technologies, integration of multi-omics data from different databases is required. Integration of omics data provides a comprehensive understanding of the biological processes underlying plant traits and their interactions. This review highlights the importance of integrating omics databases in crop plant breeding, discusses available omics data and databases, describes integration challenges, and highlights recent developments and potential benefits. Taken together, the integration of omics databases is a critical step towards enhancing crop plant breeding and improving global food security.

Evolutionarily conserved protein motifs drive interactions between the plant nucleoskeleton and nuclear pores

The nucleoskeleton forms a filamentous meshwork under the nuclear envelope and contributes to the regulation of nuclear shape and gene expression. To understand how the Arabidopsis (Arabidopsis thaliana) nucleoskeleton physically connects to the nuclear periphery in plants, we investigated the Arabidopsis nucleoskeleton protein KAKU4 and sought for functional regions responsible for its localization at the nuclear periphery. We identified 3 conserved peptide motifs within the N-terminal region of KAKU4, which are required for intermolecular interactions of KAKU4 with itself, interaction with the nucleoskeleton protein CROWDED NUCLEI (CRWN), localization at the nuclear periphery, and nuclear elongation in differentiated tissues. Unexpectedly, we find these motifs to be present also in NUP82 and NUP136, 2 plant-specific nucleoporins from the nuclear pore basket. We further show that NUP82, NUP136, and KAKU4 have a common evolutionary history predating nonvascular land plants with KAKU4 mainly localizing outside the nuclear pore suggesting its divergence from an ancient nucleoporin into a new nucleoskeleton component. Finally, we demonstrate that both NUP82 and NUP136, through their shared N-terminal motifs, interact with CRWN and KAKU4 proteins revealing the existence of a physical continuum between the nuclear pore and the nucleoskeleton in plants.

Admixture mapping and selection scans identify genomic regions associated with stomatal patterning and disease resistance in hybrid poplars

Variation in fitness components can be linked in some cases to variation in key traits. Metric traits that lie at the intersection of development, defense, and ecological interactions may be expected to experience environmental selection, informing our understanding of evolutionary and ecological processes. Here, we use quantitative genetic and population genomic methods to investigate disease dynamics in hybrid and non-hybrid populations. We focus our investigation on morphological and ecophysiological traits which inform our understanding of physiology, growth, and defense against a pathogen. In particular, we investigate stomata, microscopic pores on the surface of a leaf that regulate gas exchange during photosynthesis and are sites of entry for various plant pathogens. Stomatal patterning traits were highly predictive of disease risk. Admixture mapping identified a polygenic basis of disease resistance. Candidate genes for stomatal and disease resistance map to the same genomic regions and experienced positive selection. Genes with functions to guard cell homeostasis, the plant immune system, components of constitutive defenses, and growth-related transcription factors were identified. Our results indicate positive selection acted on candidate genes for stomatal patterning and disease resistance, potentially acting in concert to structure their variation in naturally formed backcrossing hybrid populations.

Recent genome resequencing paraded COBRA-Like gene family roles in abiotic stress and wood formation in Poplar

A cell wall determines the mechanical properties of a cell, serves as a barrier against plant stresses, and allows cell division and growth processes. The COBRA-Like (COBL) gene family encodes a putative glycosylphosphatidylinositol (GPI)-anchored protein that controls cellulose deposition and cell progression in plants by contributing to the microfibril orientation of a cell wall. Despite being studied in different plant species, there is a dearth of the comprehensive global analysis of COBL genes in poplar. Poplar is employed as a model woody plant to study abiotic stresses and biomass production in tree research. Improved genome resequencing has enabled the comprehensive exploration of the evolution and functional capacities of PtrCOBLs (Poplar COBRA-Like genes) in poplar. Phylogeny analysis has discerned and classified PtrCOBLs into two groups resembling the Arabidopsis COBL family, and group I genes possess longer proteins but have fewer exons than group II. Analysis of gene structure and motifs revealed PtrCOBLs maintained a rather stable motif and exon-intron pattern across members of the same group. Synteny and collinearity analyses exhibited that the evolution of the COBL gene family was heavily influenced by gene duplication events. PtrCOBL genes have undergone both segmental duplication and tandem duplication, followed by purifying selection. Promotor analysis flaunted various phytohormone-, growth- and stress-related cis-elements (e.g., MYB, ABA, MeJA, SA, AuxR, and ATBP1). Likewise, 29 Ptr-miRNAs of 20 families were found targeting 11 PtrCOBL genes. PtrCOBLs were found localized at the plasma membrane and extracellular matrix, while gene ontology analysis showed their involvement in plant development, plant growth, stress response, cellulose biosynthesis, and cell wall biogenesis. RNA-seq datasets depicted the bulk of PtrCOBL genes expression being found in plant stem tissues and leaves, rendering mechanical strength and rejoinders to environmental cues. PtrCOBL2, 3, 10, and 11 manifested the highest expression in vasculature and abiotic stress, and resemblant expression trends were upheld by qRT-PCR. Co-expression network analysis identified PtrCOBL2 and PtrCOBL3 as hub genes across all abiotic stresses and wood developing tissues. The current study reports regulating roles of PtrCOBLs in xylem differentiating tissues, tension wood formation, and abiotic stress latency that lay the groundwork for future functional studies of the PtrCOBL genes in poplar breeding.

Carbohydrate dynamics in Populus trees under drought: An expression atlas of genes related to sensing, translocation, and metabolism across organs

In trees, nonstructural carbohydrates (NSCs) serve as long-term carbon storage and long-distance carbon transport from source to sink. NSC management in response to drought stress is key to our understanding of drought acclimation. However, the molecular mechanisms underlying these processes remain unclear. By combining a transcriptomic approach with NSC quantification in the leaves, stems, and roots of Populus alba under drought stress, we analyzed genes from 29 gene families related to NSC signaling, translocation, and metabolism. We found starch depletion across organs and accumulation of soluble sugars (SS) in the leaves. Activation of the trehalose-6-phosphate/SNF1-related protein kinase (SnRK1) signaling pathway across organs via the suppression of class I TREHALOSE-PHOSPHATE SYNTHASE (TPS) and the expression of class II TPS genes suggested an active response to drought. The expression of SnRK1α and β subunits, and SUCROSE SYNTHASE6 supported SS accumulation in leaves. The upregulation of active transporters and the downregulation of most passive transporters implied a shift toward active sugar transport and enhanced regulation over partitioning. SS accumulation in vacuoles supports osmoregulation in leaves. The increased expression of sucrose synthesis genes and reduced expression of sucrose degradation genes in the roots did not coincide with sucrose levels, implying local sucrose production for energy. Moreover, the downregulation of invertases in the roots suggests limited sucrose allocation from the aboveground organs. This study provides an expression atlas of NSC-related genes that respond to drought in poplar trees, and can be tested in tree improvement programs for adaptation to drought conditions.

Aspartyl tRNA-synthetase (AspRS) gene family enhances drought tolerance in poplar through BABA-PtrIBIs-PtrVOZ signaling module

BACKGROUND

Drought stress is a prevalent abiotic stress that significantly hinders the growth and development of plants. According to studies, β-aminobutyric acid (BABA) can influence the ABA pathway through the AtIBI1 receptor gene to enhance cold resistance in Arabidopsis. However, the Aspartate tRNA-synthetase (AspRS) gene family, which acts as the receptor for BABA, has not yet been investigated in poplar. Particularly, it is uncertain how the AspRS gene family (PtrIBIs)r can resist drought stress after administering various concentrations of BABA to poplar.

RESULTS

In this study, we have identified 12 AspRS family genes and noted that poplar acquired four PtrIBI pairs through whole genome duplication (WGD). We conducted cis-action element analysis and found a significant number of stress-related action elements on different PtrIBI genes promoters. The expression of most PtrIBI genes was up-regulated under beetle and mechanical damage stresses, indicating their potential role in responding to leaf damage stress. Our results suggest that a 50 mM BABA treatment can alleviate the damage caused by drought stress in plants. Additionally, via transcriptome sequencing, we observed that the partial up-regulation of BABA receptor genes, PtrIBI2/4/6/8/11, in poplars after drought treatment. We hypothesize that poplar responds to drought stress through the BABA-PtrIBIs-PtrVOZ coordinated ABA signaling pathway. Our research provides molecular evidence for understanding how plants respond to drought stress through external application of BABA.

CONCLUSIONS

In summary, our study conducted genome-wide analysis of the AspRS family of P. trichocarpa and identified 12 PtrIBI genes. We utilized genomics and bioinformatics to determine various characteristics of PtrIBIs such as chromosomal localization, evolutionary tree, gene structure, gene doubling, promoter cis-elements, and expression profiles. Our study found that certain PtrIBI genes are regulated by drought, beetle, and mechanical damage implying their crucial role in enhancing poplar stress tolerance. Additionally, we observed that external application of low concentrations of BABA increased plant drought resistance under drought stress. Through the BABA-PtrIBIs-PtrVOZ signaling module, poplar plants were able to transduce ABA signaling and regulate their response to drought stress. These results suggest that the PtrIBI genes in poplar have the potential to improve drought tolerance in plants through the topical application of low concentrations of BABA.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Salicylic acid metabolism and signalling coordinate senescence initiation in aspen in nature

Deciduous trees exhibit a spectacular phenomenon of autumn senescence driven by the seasonality of their growth environment, yet there is no consensus which external or internal cues trigger it. Senescence starts at different times in European aspen (Populus tremula L.) genotypes grown in same location. By integrating omics studies, we demonstrate that aspen genotypes utilize similar transcriptional cascades and metabolic cues to initiate senescence, but at different times during autumn. The timing of autumn senescence initiation appeared to be controlled by two consecutive "switches"; 1) first the environmental variation induced the rewiring of the transcriptional network, stress signalling pathways and metabolic perturbations and 2) the start of senescence process was defined by the ability of the genotype to activate and sustain stress tolerance mechanisms mediated by salicylic acid. We propose that salicylic acid represses the onset of leaf senescence in stressful natural conditions, rather than promoting it as often observed in annual plants.

Genome-Wide Analysis of the FBA Subfamily of the Poplar F-Box Gene Family and Its Role under Drought Stress

F-box proteins are important components of eukaryotic SCF E3 ubiquitin ligase complexes, which specifically determine protein substrate proteasomal degradation during plant growth and development, as well as biotic and abiotic stress. It has been found that the FBA (F-box associated) protein family is one of the largest subgroups of the widely prevalent F-box family and plays significant roles in plant development and stress response. However, the FBA gene family in poplar has not been systematically studied to date. In this study, a total of 337 F-box candidate genes were discovered based on the fourth-generation genome resequencing of P. trichocarpa. The domain analysis and classification of candidate genes revealed that 74 of these candidate genes belong to the FBA protein family. The poplar F-box genes have undergone multiple gene replication events, particularly in the FBA subfamily, and their evolution can be attributed to genome-wide duplication (WGD) and tandem duplication (TD). In addition, we investigated the P. trichocarpa FBA subfamily using the PlantGenIE database and quantitative real-time PCR (qRT-PCR); the results showed that they are expressed in the cambium, phloem and mature tissues, but rarely expressed in young leaves and flowers. Moreover, they are also widely involved in the drought stress response. At last, we selected and cloned PtrFBA60 for physiological function analysis and found that it played an important role in coping with drought stress. Taken together, the family analysis of FBA genes in P. trichocarpa provides a new opportunity for the identification of P. trichocarpa candidate FBA genes and elucidation of their functions in growth, development and stress response, thus demonstrating their utility in the improvement of P. trichocarpa.

From common gardens to candidate genes: exploring local adaptation to climate in red spruce

Local adaptation to climate is common in plant species and has been studied in a range of contexts, from improving crop yields to predicting population maladaptation to future conditions. The genomic era has brought new tools to study this process, which was historically explored through common garden experiments. In this study, we combine genomic methods and common gardens to investigate local adaptation in red spruce and identify environmental gradients and loci involved in climate adaptation. We first use climate transfer functions to estimate the impact of climate change on seedling performance in three common gardens. We then explore the use of multivariate gene-environment association methods to identify genes underlying climate adaptation, with particular attention to the implications of conducting genome scans with and without correction for neutral population structure. This integrative approach uncovered phenotypic evidence of local adaptation to climate and identified a set of putatively adaptive genes, some of which are involved in three main adaptive pathways found in other temperate and boreal coniferous species: drought tolerance, cold hardiness, and phenology. These putatively adaptive genes segregated into two 'modules' associated with different environmental gradients. This study nicely exemplifies the multivariate dimension of adaptation to climate in trees.

Orthologs of Human-Disease-Associated Genes in Plants Are Involved in Regulating Leaf Senescence

As eukaryotes, plants and animals have many commonalities on the genetic level, although they differ greatly in appearance and physiological habits. The primary goal of current plant research is to improve the crop yield and quality. However, plant research has a wider aim, exploiting the evolutionary conservatism similarities between plants and animals, and applying discoveries in the field of botany to promote zoological research that will ultimately serve human health, although very few studies have addressed this aspect. Here, we analyzed 35 human-disease-related gene orthologs in plants and characterized the genes in depth. Thirty-four homologous genes were found to be present in the herbaceous annual plant Arabidopsis thaliana and the woody perennial plant Populus trichocarpa, with most of the genes having more than two exons, including the ATM gene with 78 exons. More surprisingly, 27 (79.4%) of the 34 homologous genes in Arabidopsis were found to be senescence-associated genes (SAGs), further suggesting a close relationship between human diseases and cellular senescence. Protein-protein interaction network analysis revealed that the 34 genes formed two main subnetworks, and genes in the first subnetwork interacted with 15 SAGs. In conclusion, our results show that most of the 34 homologs of human-disease-associated genes in plants are involved in the leaf senescence process, suggesting that leaf senescence may offer a means to study the pathogenesis of human diseases and to screen drugs for the treat of diseases.

Populus endo-glucanase 16 localizes to the cell walls of developing tissues

The hemicelluloses comprise a group of matrix glycans that interact with cellulose microfibrils in plant cell walls and play important roles in establishing wall architecture. The structures of hemicelluloses are determined by carbohydrate-active enzymes (CAZymes) that synthesize, integrate, and break down these polymers. Specifically, endo-glucanase 16 (EG16) enzymes, which are related to the well-known xyloglucan endotransglycosylase/hydrolase (XTH) gene products in Glycoside Hydrolase Family 16 (GH16), have been implicated in the degradation of the β(1,4)-linked backbone of mixed-linkage β(1,3);β(1,4)-glucans (MLG) and xyloglucans. EG16 members are single-copy genes found in most plant clades but are absent from many eudicots, including the model plant Arabidopsis thaliana. Until recently, EG16 members had only been characterized in vitro, establishing their substrate specificity, protein structure, and phylogenetic history, but their biological function was unknown. Here we used a hybrid polar, Populus alba × Populus grandidentata (P39), as a model to examine EG16 expression, subcellular localization, and pheno- and chemotypes of EG16-downregulated P39 plants. Populus EG16 expression is strong in young tissues, but RNAi-mediated downregulation did not impact plant growth nor the fine structure of the hemicellulose xyloglucan, suggesting a restricted or currently unknown role in angiosperm physiology.

Identification and expression pattern analysis of PtCarA and PtCarB genes in Populus trichocarpa under different nitrogen treatments

Carbamoyl phosphate synthetase (CPS) catalyses the synthesis of ammonia carbamoyl phosphate (CP), which plays a key role in the biosynthesis of arginine and pyrimidine nucleotides. There are two subunits of the CPS enzyme in Populus trichocarpa, CarA (small subunit) and CarB (large subunit). Only when they coexist can CPS catalyse synthesis of CP. However, it is not clear how CPS responds to nitrogen (N) to affect arginine and pyrimidine nucleotide biosynthesis. In this study, bioinformatics methods were used to analyse the expression patterns of genes encoding CarA and CarB, and qRT-PCR and RNA-seq were used to investigate their molecular responses under different N concentrations. Phylogenetic analysis revealed that the phylogenetic trees of CarA and CarB had similar topologies. qRT-PCR showed that the PtCarA and PtCarB genes were regulated by N, while their N-regulated patterns differed in different tissues. The expression patterns of PtCarA and PtCarB show a significant positive correlation according to qRT-PCR and RNA-seq. The analysis of promoter cis-acting elements showed that the promoter regions of PtCarA1, PtCarA2 and PtCarB contained some identical cis-acting elements. According to analysis of the phylogenetic tree, expression patterns and promoter elements, we speculate that there might be coevolution among PtCarA1, PtCarA2 and PtCarB. This study provides valuable information for further understanding the function of CPS in poplar, especially for N response, and provides new ideas for studying the evolution of gene families related to heteromultimers.

Cone-setting in spruce is regulated by conserved elements of the age-dependent flowering pathway

Reproductive phase change is well characterized in angiosperm model species, but less studied in gymnosperms. We utilize the early cone-setting acrocona mutant to study reproductive phase change in the conifer Picea abies (Norway spruce), a gymnosperm. The acrocona mutant frequently initiates cone-like structures, called transition shoots, in positions where wild-type P. abies always produces vegetative shoots. We collect acrocona and wild-type samples, and RNA-sequence their messenger RNA (mRNA) and microRNA (miRNA) fractions. We establish gene expression patterns and then use allele-specific transcript assembly to identify mutations in acrocona. We genotype a segregating population of inbred acrocona trees. A member of the SQUAMOSA BINDING PROTEIN-LIKE (SPL) gene family, PaSPL1, is active in reproductive meristems, whereas two putative negative regulators of PaSPL1, miRNA156 and the conifer specific miRNA529, are upregulated in vegetative and transition shoot meristems. We identify a mutation in a putative miRNA156/529 binding site of the acrocona PaSPL1 allele and show that the mutation renders the acrocona allele tolerant to these miRNAs. We show co-segregation between the early cone-setting phenotype and trees homozygous for the acrocona mutation. In conclusion, we demonstrate evolutionary conservation of the age-dependent flowering pathway and involvement of this pathway in regulating reproductive phase change in the conifer P. abies.

Histone H2B.8 compacts flowering plant sperm through chromatin phase separation

Sperm chromatin is typically transformed by protamines into a compact and transcriptionally inactive state^1,2. Sperm cells of flowering plants lack protamines, yet they have small, transcriptionally active nuclei with chromatin condensed through an unknown mechanism^3,4. Here we show that a histone variant, H2B.8, mediates sperm chromatin and nuclear condensation in Arabidopsis thaliana. Loss of H2B.8 causes enlarged sperm nuclei with dispersed chromatin, whereas ectopic expression in somatic cells produces smaller nuclei with aggregated chromatin. This result demonstrates that H2B.8 is sufficient for chromatin condensation. H2B.8 aggregates transcriptionally inactive AT-rich chromatin into phase-separated condensates, which facilitates nuclear compaction without reducing transcription. Reciprocal crosses show that mutation of h2b.8 reduces male transmission, which suggests that H2B.8-mediated sperm compaction is important for fertility. Altogether, our results reveal a new mechanism of nuclear compaction through global aggregation of unexpressed chromatin. We propose that H2B.8 is an evolutionary innovation of flowering plants that achieves nuclear condensation compatible with active transcription.

Network-Based Analysis to Identify Hub Genes Involved in Spatial Root Response to Mechanical Constrains

Previous studies report that the asymmetric response, observed along the main poplar woody bent root axis, was strongly related to both the type of mechanical forces (compression or tension) and the intensity of force displacement. Despite a large number of targets that have been proposed to trigger this asymmetry, an understanding of the comprehensive and synergistic effect of the antistress spatially related pathways is still lacking. Recent progress in the bioinformatics area has the potential to fill these gaps through the use of in silico studies, able to investigate biological functions and pathway overlaps, and to identify promising targets in plant responses. Presently, for the first time, a comprehensive network-based analysis of proteomic signatures was used to identify functions and pivotal genes involved in the coordinated signalling pathways and molecular activities that asymmetrically modulate the response of different bent poplar root sectors and sides. To accomplish this aim, 66 candidate proteins, differentially represented across the poplar bent root sides and sectors, were grouped according to their abundance profile patterns and mapped, together with their first neighbours, on a high-confidence set of interactions from STRING to compose specific cluster-related subnetworks (I–VI). Successively, all subnetworks were explored by a functional gene set enrichment analysis to identify enriched gene ontology terms. Subnetworks were then analysed to identify the genes that are strongly interconnected with other genes (hub gene) and, thus, those that have a pivotal role in the bent root asymmetric response. The analysis revealed novel information regarding the response coordination, communication, and potential signalling pathways asymmetrically activated along the main root axis, delegated mainly to Ca2+ (for new lateral root formation) and ROS (for gravitropic response and lignin accumulation) signatures. Furthermore, some of the data indicate that the concave side of the bent sector, where the mechanical forces are most intense, communicates to the other (neighbour and distant) sectors, inducing spatially related strategies to ensure water uptake and accompanying cell modification. This information could be critical for understanding how plants maintain and improve their structural integrity—whenever and wherever it is necessary—in natural mechanical stress conditions.

Laccaria bicolor pectin methylesterases are involved in ectomycorrhiza development with Populus tremula × Populus tremuloides

The development of ectomycorrhizal (ECM) symbioses between soil fungi and tree roots requires modification of root cell walls. The pectin-mediated adhesion between adjacent root cells loosens to accommodate fungal hyphae in the Hartig net, facilitating nutrient exchange between partners. We investigated the role of fungal pectin modifying enzymes in Laccaria bicolor for ECM formation with Populus tremula × Populus tremuloides. We combine transcriptomics of cell-wall-related enzymes in both partners during ECM formation, immunolocalisation of pectin (Homogalacturonan, HG) epitopes in different methylesterification states, pectin methylesterase (PME) activity assays and functional analyses of transgenic L. bicolor to uncover pectin modification mechanisms and the requirement of fungal pectin methylesterases (LbPMEs) for ECM formation. Immunolocalisation identified remodelling of pectin towards de-esterified HG during ECM formation, which was accompanied by increased LbPME1 expression and PME activity. Overexpression or RNAi of the ECM-induced LbPME1 in transgenic L. bicolor lines led to reduced ECM formation. Hartig Nets formed with LbPME1 RNAi lines were shallower, whereas those formed with LbPME1 overexpressors were deeper. This suggests that LbPME1 plays a role in ECM formation potentially through HG de-esterification, which initiates loosening of adjacent root cells to facilitate Hartig net formation.

Enhanced lignin synthesis and ecotypic variation in defense-related gene expression in response to shade in Norway spruce

During the growth season, northern forests in Sweden daily receive more hours of far-red (FR)-enriched light or twilight (shade) as compared to southern forests. Norway spruce (shade-tolerant) are adapted to latitudinal variation in twilight characterized by a northward increase in FR requirement to maintain growth. Shade is a stressful condition that affects plant growth and increases plant's susceptibility to pathogen attack. Lignin plays a central role in plant defense and its metabolism is regulated by light wavelength composition (light quality). In the current work, we studied regulation of lignin synthesis and defense-related genes (growth-defense trade-offs) in response to shade in Norway spruce. In most angiosperms, light promotes lignin synthesis, whereas shade decreases lignin production leading to weaker stem, which may make plants more disease susceptible. In contrast, enhanced lignin synthesis was detected in response to shade in Norway spruce. We detected a higher number of immunity/defense-related genes up-regulated in northern populations as compared to south ones in response to shade. Enhanced lignin synthesis coupled with higher defense-related gene expression can be interpreted as an adaptive strategy for better survival in northern populations. Findings will contribute to ensuring deployment of well-adapted genetic material and identifying tree families with enhanced disease resistance.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

Genome-Wide Classification and Evolutionary Analysis Reveal Diverged Patterns of Chalcone Isomerase in Plants

Flavonoids as a class of important secondary metabolites are widely present in land plants, and chalcone isomerase (CHI) is the key rate-limiting enzyme that participates in catalyzing the stereospecific isomerization of chalcones to yield their corresponding flavanones. However, the phylogenetic dynamics and functional divergence of CHI family genes during the evolutionary path of green plants remains poorly understood. Here, a total of 122 CHI genes were identified by performing a genome-wide survey of 15 representative green plants from the most ancestral basal plant chlorophyte algae to higher angiosperm plants. Phylogenetic, orthologous groups (OG) classification, and genome structure analysis showed that the CHI family genes have evolved into four distinct types (types I–IV) containing eight OGs after gene duplication, and further studies indicated type III CHIs consist of three subfamilies (FAP1, FAP2, and FAP3). The phylogeny showed FAP3 CHIs as an ancestral out-group positioned on the outer layers of the main branch, followed by type IV CHIs, which are placed in an evolutionary intermediate between FAP3 CHIs and bona fide CHIs (including type I and type II). The results imply a potential intrinsic evolutionary connection between CHIs existing in the green plants. The amino acid substitutions occurring in several residues have potentially affected the functional divergence between CHI proteins. This is supported by the analysis of transcriptional divergence and cis-acting element analysis. Evolutionary dynamics analyses revealed that the differences in the total number of CHI family genes in each plant are primarily attributed to the lineage-specific expansion by natural selective forces. The current studies provide a deeper understanding of the phylogenetic relationships and functional diversification of CHI family genes in green plants, which will guide further investigation on molecular characteristics and biological functions of CHIs.

Advances in understanding Norway spruce natural resistance to needle bladder rust infection: transcriptional and secondary metabolites profiling

Background

Needle rust caused by the fungus Chrysomyxa rhododendri causes significant growth decline and increased mortality of young Norway spruce trees in subalpine forests. Extremely rare trees with enhanced resistance represent promising candidates for practice-oriented reproduction approaches. They also enable the investigation of tree molecular defence and resistance mechanisms against this fungal disease. Here, we combined RNA-Seq, RT-qPCR and secondary metabolite analyses during a period of 38 days following natural infection to investigate differences in constitutive and infection-induced defence between the resistant genotype PRA-R and three susceptible genotypes.

Results

Gene expression and secondary metabolites significantly differed among genotypes from day 7 on and revealed already known, but also novel candidate genes involved in spruce molecular defence against this pathogen. Several key genes related to (here and previously identified) spruce defence pathways to needle rust were differentially expressed in PRA-R compared to susceptible genotypes, both constitutively (in non-symptomatic needles) and infection-induced (in symptomatic needles). These genes encoded both new and well-known antifungal proteins such as endochitinases and chitinases. Specific genetic characteristics concurred with varying phenolic, terpene, and hormone needle contents in the resistant genotype, among them higher accumulation of several flavonoids (mainly kaempferol and taxifolin), stilbenes, geranyl acetone, α-ionone, abscisic acid and salicylic acid.

Conclusions

Combined transcriptional and metabolic profiling of the Norway spruce defence response to infection by C. rhododendri in adult trees under subalpine conditions confirmed the results previously gained on artificially infected young clones in the greenhouse, both regarding timing and development of infection, and providing new insights into genes and metabolic pathways involved. The comparison of genotypes with different degrees of susceptibility proved that several of the identified key genes are differently regulated in PRA-R, and that the resistant genotype combines a strong constitutive defence with an induced response in infected symptomatic needles following fungal invasion. Genetic and metabolic differences between the resistant and susceptible genotypes indicated a more effective hypersensitive response (HR) in needles of PRA-R that prevents penetration and spread of the rust fungus and leads to a lower proportion of symptomatic needles as well as reduced symptom development on the few affected needles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08661-y.

Transcriptomic changes during the establishment of long-term methyl jasmonate-induced resistance in Norway spruce

Norway spruce (Picea abies) is an economically and ecologically important tree species that grows across northern and central Europe. Treating Norway spruce with jasmonate has long-lasting beneficial effects on tree resistance to damaging pests, such as the European spruce bark beetle Ips typographus and its fungal associates. The (epi)genetic mechanisms involved in such long-lasting jasmonate induced resistance (IR) have gained much recent interest but remain largely unknown. In this study, we treated 2-year-old spruce seedlings with methyl jasmonate (MeJA) and challenged them with the I. typographus vectored necrotrophic fungus Grosmannia penicillata. MeJA treatment reduced the extent of necrotic lesions in the bark 8 weeks after infection and thus elicited long-term IR against the fungus. The transcriptional response of spruce bark to MeJA treatment was analysed over a 4-week time course using mRNA-seq. This analysis provided evidence that MeJA treatment induced a transient upregulation of jasmonic acid, salicylic acid and ethylene biosynthesis genes and downstream signalling genes. Our data also suggests that defence-related genes are induced while genes related to growth are repressed by methyl jasmonate treatment. These results provide new clues about the potential underpinning mechanisms and costs associated with long-term MeJA-IR in Norway spruce.

Diversification of SEC15a and SEC15b isoforms of an exocyst subunit in seed plants is manifested in their specific roles in Arabidopsis sporophyte and male gametophyte

The exocyst complex is an octameric evolutionarily conserved tethering complex engaged in the regulation of polarized secretion in eukaryotic cells. Here, we focus on the systematic comparison of two isoforms of the SEC15 exocyst subunit, SEC15a and SEC15b. We infer that SEC15 gene duplication and diversification occurred in the common ancestor of seed plants (Spermatophytes). In Arabidopsis, SEC15a represents the main SEC15 isoform in the male gametophyte, and localizes to the pollen tube tip at the plasma membrane. Although pollen tubes of sec15a mutants are impaired, sporophytes show no phenotypic deviations. Conversely, SEC15b is the dominant isoform in the sporophyte and localizes to the plasma membrane in root and leaf cells. Loss-of-function sec15b mutants exhibit retarded elongation of hypocotyls and root hairs, a loss of apical dominance, dwarfed plant stature and reduced seed coat mucilage formation. Surprisingly, the sec15b mutants also exhibit compromised pollen tube elongation in vitro, despite its very low expression in pollen, suggesting a non-redundant role for the SEC15b isoform there. In pollen tubes, SEC15b localizes to distinct cytoplasmic structures. Reciprocally to this, SEC15a also functions in the sporophyte, where it accumulates at plasmodesmata. Importantly, although overexpressed SEC15a could fully complement the sec15b phenotypic deviations in the sporophyte, the pollen-specific overexpression of SEC15b was unable to fully compensate for the loss of SEC15a function in pollen. We conclude that the SEC15a and SEC15b isoforms evolved in seed plants, with SEC15a functioning mostly in pollen and SEC15b functioning mostly in the sporophyte.

FLOWERING LOCUS T paralogs control the annual growth cycle in Populus trees

In temperate and boreal regions, perennials adapt their annual growth cycle to the change of seasons. These adaptations ensure survival in harsh environmental conditions, allowing growth at different latitudes and altitudes, and are therefore tightly regulated. Populus tree species cease growth and form terminal buds in autumn when photoperiod falls below a certain threshold.¹ This is followed by establishment of dormancy and cold hardiness over the winter. At the center of the photoperiodic pathway in Populus is the gene FLOWERING LOCUS T2 (FT2), which is expressed during summer and harbors significant SNPs in its locus associated with timing of bud set.^1-4 The paralogous gene FT1, on the other hand, is hyper-induced in chilling buds during winter.^3,5 Even though its function is so far unknown, it has been suggested to be involved in the regulation of flowering and the release of winter dormancy.^3,5 In this study, we employ CRISPR-Cas9-mediated gene editing to individually study the function of the FT-like genes in Populus trees. We show that while FT2 is required for vegetative growth during spring and summer and regulates the entry into dormancy, expression of FT1 is absolutely required for bud flush in spring. Gene expression profiling suggests that this function of FT1 is linked to the release of winter dormancy rather than to the regulation of bud flush per se. These data show how FT duplication and sub-functionalization have allowed Populus trees to regulate two completely different and major developmental control points during the yearly growth cycle.

Leaf necrosis resulting from downregulation of poplar glycosyltransferase UGT72A2

Reactive species (RS) causing oxidative stress are unavoidable by-products of various plant metabolic processes, such as photosynthesis, respiration or photorespiration. In leaves, flavonoids scavenge RS produced during photosynthesis and protect plant cells against deleterious oxidative damages. Their biosynthesis and accumulation are therefore under tight regulation at the cellular level. Glycosylation has emerged as an essential biochemical reaction in the homeostasis of various specialized metabolites such as flavonoids. This article provides a functional characterization of the Populus tremula x P. alba (poplar) UGT72A2 coding for a UDP-glycosyltransferase that is localized in the chloroplasts. Compared with the wild type, transgenic poplar lines with decreased expression of UGT72A2 are characterized by reduced growth and oxidative damages in leaves, as evidenced by necrosis, higher content of glutathione and lipid peroxidation products as well as diminished soluble peroxidase activity and NADPH to NADP+ ratio under standard growing conditions. They furthermore display lower pools of phenolics, anthocyanins and total flavonoids but higher proanthocyanidins content. Promoter analysis revealed the presence of cis-elements involved in photomorphogenesis, chloroplast biogenesis and flavonoid biosynthesis. The UGT72A2 is regulated by the poplar MYB119, a transcription factor known to regulate the flavonoid biosynthesis pathway. Phylogenetic analysis and molecular docking suggest that UGT72A2 could glycosylate flavonoids; however, the actual substrate(s) was not consistently evidenced with either in vitro assays nor analyses of glycosylated products in leaves of transgenic poplar overexpressing or downregulated for UGT72A2. This article provides elements highlighting the importance of flavonoid glycosylation regarding protection against oxidative stress in poplar leaves and raises new questions about the link between this biochemical reaction and regulation of the redox homeostasis system.

Transcriptional cascades in the regulation of 2-AP biosynthesis under Zn supply in fragrant rice

Transcription factors (TFs) regulate gene expression to control certain genetic programs, such as growth and development, phytohormone regulation, and environmental stresses. 2-acetyl-1-pyrroline (2-AP) is the key element involved in aroma biosynthesis pathway, and the application of micronutrients can increase the 2-AP levels. However, little is known about the micronutrient-induced TFs involved in 2-AP biosynthesis. Here, we identify a number of TF families in two fragrant rice varieties, "Meixiangzhan-2" (M) and "Xiangyaxiangzhan" (X), in response to Zinc (Zn) application through transcriptomic analysis. A total of ~678 TFs were identified and grouped into 26 TF families, each of which was found to be involved in numerous signaling pathways. The WRKY TF family was found to be the most abundant, followed by bHLH and MYB. Furthermore, members of the WRKY, bHLH, MYB, ERF, HSF, MADS-box, NFY, and AP2 TF families were significantly upregulated and may be involved in the transcriptional regulation of aroma biosynthesis. In brief, this study enhances our understanding of the molecular mechanism of 2-AP biosynthesis and highlights the key TFs potentially involved in the production of aroma in fragrant rice.

DUF1005 Family Identification, Evolution Analysis in Plants, and Primary Root Elongation Regulation of CiDUF1005 From Caragana intermedia

Proteins with a domain of unknown function (DUF) represent a number of gene families that encode functionally uncharacterized proteins in eukaryotes. In particular, members of the DUF1005 family in plants have a 411-amino-acid conserved domain, and this family has not been described previously. In this study, a total of 302 high-confidence DUF1005 family members were identified from 58 plant species, and none were found in the four algae that were selected. Thus, this result showed that DUF1005s might belong to a kind of plant-specific gene family, and this family has not been evolutionarily expanded. Phylogenetic analysis showed that the DUF1005 family genes could be classified into four subgroups in 58 plant species. The earliest group to emerge was Group I, including a total of 100 gene sequences, and this group was present in almost all selected species spanning from mosses to seed plants. Group II and Group III, with 69 and 74 members, respectively, belong to angiosperms. Finally, with 59 members, Group IV was the last batch of genes to emerge, and this group is unique to dicotyledons. Expression pattern analysis of the CiDUF1005, a member of the DUF1005 family from Caragana intermedia, showed that CiDUF1005 genes were differentially regulated under various treatments. Compared to the wild type, transgenic lines with heterologous CiDUF1005 expression in Arabidopsis thaliana had longer primary roots and more lateral roots. These results expanded our knowledge of the evolution of the DUF1005 family in plants and will contribute to elucidating biological functions of the DUF1005 family in the future.

miR160 Interacts in vivo With Pinus pinaster AUXIN RESPONSE FACTOR 18 Target Site and Negatively Regulates Its Expression During Conifer Somatic Embryo Development

MicroRNAs (miRNAs) are key regulators of several plant developmental processes including embryogenesis. Most miRNA families are conserved across major groups of plant species, but their regulatory roles have been studied mainly in model species like Arabidopsis and other angiosperms. In gymnosperms, miRNA-dependent regulation has been less studied since functional approaches in these species are often difficult to establish. Given the fundamental roles of auxin signaling in somatic embryogenesis (SE) induction and embryo development, we investigated a previously predicted interaction between miR160 and a putative target encoding AUXIN RESPONSE FACTOR 18 in Pinus pinaster (PpARF18) embryonic tissues. Phylogenetic analysis of AUXIN RESPONSE FACTOR 18 (ARF18) from Pinus pinaster and Picea abies, used here as a model system of conifer embryogenesis, showed their close relatedness to AUXIN RESPONSE FACTOR (ARF) genes known to be targeted by miR160 in other species, including Arabidopsis ARF10 and ARF16. By using a luciferase (LUC) reporter system for miRNA activity in Arabidopsis protoplasts, we have confirmed that P. pinaster miR160 (ppi-miR160) interacts in vivo with PpARF18 target site. When the primary miR160 from P. pinaster was overexpressed in protoplasts under non-limiting levels of ARGONAUTE1, a significant increase of miR160 target cleavage activity was observed. In contrast, co-expression of the primary miRNA and the target mimic MIM160 led to a decrease of miR160 activity. Our results further support that this interaction is functional during consecutive stages of SE in the conifer model P. abies. Expression analyses conducted in five stages of development, from proembryogenic masses (PEMs) to the mature embryo, show that conifer ARF18 is negatively regulated by miR160 toward the fully developed mature embryo when miR160 reached its highest expression level. This study reports the first in vivo validation of a predicted target site of a conifer miRNA supporting the conservation of miR160 interaction with ARF targets in gymnosperms. The approach used here should be useful for future characterization of miRNA functions in conifer embryogenesis.

Evolution of chlorophyll degradation is associated with plant transition to land

Chlorophyll, the central pigment of photosynthesis, is highly photo‐active and degraded enzymatically during leaf senescence. Merging comparative genomics and metabolomics, we evaluate the extent to which the chlorophyll detoxification pathway has evolved in Viridiplantae. We argue that cytosolic detoxification of phyllobilins in particular was a critical process to the green lineage’s transition to land.

Populus SVL Acts in Leaves to Modulate the Timing of Growth Cessation and Bud Set

SHORT VEGETATIVE PHASE (SVP) is an important regulator of FLOWERING LOCUS T (FT) in the thermosensory pathway of Arabidopsis. It is a negative regulator of flowering and represses FT transcription. In poplar trees, FT2 is central for the photoperiodic control of growth cessation, which also requires the decrease of bioactive gibberellins (GAs). In angiosperm trees, genes similar to SVP, sometimes named DORMANCY-ASSOCIATED MADS-BOX genes, control temperature-mediated bud dormancy. Here we show that SVL, an SVP ortholog in aspen trees, besides its role in controlling dormancy through its expression in buds, is also contributing to the regulation of short day induced growth cessation and bud set through its expression in leaves. SVL is upregulated during short days in leaves and binds to the FT2 promoter to repress its transcription. It furthermore decreases the amount of active GAs, whose downregulation is essential for growth cessation, by repressing the transcription of GA20 oxidase. Finally, the SVL protein is more stable in colder temperatures, thus integrating the temperature signal into the response. We conclude that the molecular function of SVL in the photoperiodic pathway has been conserved between Arabidopsis and poplar trees, albeit the physiological process it controls has changed. SVL is thus both involved in regulating the photoperiod response in leaves, modulating the timing of growth cessation and bud set, and in the subsequent temperature regulation of dormancy in the buds.

pHBMT1, a BAHD-family monolignol acyltransferase, mediates lignin acylation in poplar

Poplar (Populus) lignin is naturally acylated with p-hydroxybenzoate ester moieties. However, the enzyme(s) involved in the biosynthesis of the monolignol-p-hydroxybenzoates have remained largely unknown. Here, we performed an in vitro screen of the Populus trichocarpa BAHD acyltransferase superfamily (116 genes) using a wheatgerm cell-free translation system and found five enzymes capable of producing monolignol-p-hydroxybenzoates. We then compared the transcript abundance of the five corresponding genes with p-hydroxybenzoate concentrations using naturally occurring unrelated genotypes of P. trichocarpa and revealed a positive correlation between the expression of p-hydroxybenzoyl-CoA monolig-nol transferase (pHBMT1, Potri.001G448000) and p-hydroxybenzoate levels. To test whether pHBMT1 is responsible for the biosynthesis of monolignol-p-hydroxybenzoates, we overexpressed pHBMT1 in hybrid poplar (Populus alba × P. grandidentata) (35S::pHBMT1 and C4H::pHBMT1). Using three complementary analytical methods, we showed that there was an increase in soluble monolignol-p-hydroxybenzoates and cell-wall-bound monolignol-p-hydroxybenzoates in the poplar transgenics. As these pendent groups are ester-linked, saponification releases p-hydroxybenzoate, a precursor to parabens that are used in pharmaceuticals and cosmetics. This identified gene could therefore be used to engineer lignocellulosic biomass with increased value for emerging biorefinery strategies.

Norway spruce deploys tissue-specific responses during acclimation to cold

Climate change in the conifer-dominated boreal forest is expected to lead to warmer but more dynamic winter air temperatures, reducing the depth and duration of snow cover and lowering winter soil temperatures. To gain insight into the mechanisms that have enabled conifers to dominate extreme cold environments, we performed genome-wide RNA-Seq analysis from needles and roots of non-dormant two-year Norway spruce (Picea abies (L.) H. Karst), and contrasted these response to herbaceous model Arabidopsis We show that the main transcriptional response of Norway spruce needles exposed to cold was delayed relative to Arabidopsis, and this delay was associated with slower development of freezing tolerance. Despite this difference in timing, Norway spruce principally utilizes early response transcription factors (TFs) belonging to the same gene families as Arabidopsis, indicating broad evolutionary conservation of cold response networks. In keeping with their different metabolic and developmental states, needles and root of Norway spruce showed contrasting results. Regulatory network analysis identified both conserved TFs with known roles in cold acclimation (e.g. homologs of ICE1, AKS3, and of the NAC and AP2/ERF superfamilies), but also a root-specific bHLH101 homolog, providing functional insights into cold stress response strategies in Norway spruce.

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Functional Analysis of Poplar Sombrero-Type NAC Transcription Factors Yields a Strategy to Modify Woody Cell Wall Properties

Woody cells generate lignocellulosic biomass, which is a promising sustainable bioresource for wide industrial applications. Woody cell differentiation in vascular plants, including the model plant poplar (Populus trichocarpa), is regulated by a set of NAC family transcription factors, the VASCULAR-RELATED NAC-DOMAIN (VND), NAC SECONDARY CELL WALL THICKENING PROMOTING FACTOR (NST)/SND, and SOMBRERO (SMB) (VNS)-related proteins, but the precise contributions of each VNS protein to wood quality are unknown. Here, we performed a detailed functional analysis of the poplar SMB-type VNS proteins PtVNS13-PtVNS16. PtVNS13-PtVNS16 were preferentially expressed in the roots of young poplar plantlets, similar to the Arabidopsis thalianaSMB gene. PtVNS13 and PtVNS14, as well as the NST-type PtVNS11, suppressed the abnormal root cap phenotype of the Arabidopsis sombrero-3 mutant, whereas the VND-type PtVNS07 gene did not, suggesting a functional gap between SMB- or NST-type VNS proteins and VND-type VNS proteins. Overexpressing PtVNS13-PtVNS16 in Arabidopsis seedlings and poplar leaves induced ectopic xylem-vessel-like cells with secondary wall deposition, and a transient expression assay showed that PtVNS13-16 transactivated woody-cell-related genes. Interestingly, although any VNS protein rescued the pendant stem phenotype of the Arabidopsis nst1-1 nst3-1 mutant, the resulting inflorescence stems exhibited distinct cell wall properties: poplar VNS genes generated woody cell walls with higher enzymatic saccharification efficiencies compared with Arabidopsis VNS genes. Together, our data reveal clear functional diversity among VNS proteins in woody cell differentiation and demonstrate a novel VNS-based strategy for modifying woody cell wall properties toward enhanced utilization of woody biomass.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BACKGROUND

Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.

RESULTS

We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.

CONCLUSIONS

qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Cytotype and genotype predict mortality and recruitment in Colorado quaking aspen (Populus tremuloides)

Species responses to climate change depend on environment, genetics, and interactions among these factors. Intraspecific cytotype (ploidy level) variation is a common type of genetic variation in many species. However, the importance of intraspecific cytotype variation in determining demography across environments is poorly known. We studied quaking aspen (Populus tremuloides), which occurs in diploid and triploid cytotypes. This widespread tree species is experiencing contractions in its western range, which could potentially be linked to cytotype-dependent drought tolerance. We found that interactions between cytotype and environment drive mortality and recruitment across 503 plots in Colorado. Triploids were more vulnerable to mortality relative to diploids and had reduced recruitment on more drought-prone and disturbed plots relative to diploids. Furthermore, there was substantial genotype-dependent variation in demography. Thus, cytotype and genotype variation are associated with decline in this foundation species. Future assessment of demographic responses to climate change will benefit from knowledge of how genetic and environmental mosaics interact to determine species' ecophysiology and demography.

An atlas of the Norway spruce needle seasonal transcriptome

Boreal conifers possess a tremendous ability to survive and remain evergreen during harsh winter conditions and resume growth during summer. This is enabled by coordinated regulation of major cellular functions at the level of gene expression, metabolism, and physiology. Here we present a comprehensive characterization of the annual changes in the global transcriptome of Norway spruce (Picea abies) needles as a resource to understand needle development and acclimation processes throughout the year. In young, growing needles (May 15 until June 30), cell walls, organelles, etc., were formed, and this developmental program heavily influenced the transcriptome, explained by over-represented Gene Ontology (GO) categories. Later changes in gene expression were smaller but four phases were recognized: summer (July-August), autumn (September-October), winter (November-February), and spring (March-April), where over-represented GO categories demonstrated how the needles acclimated to the various seasons. Changes in the seasonal global transcriptome profile were accompanied by differential expression of members of the major transcription factor families. We present a tentative model of how cellular activities are regulated over the year in needles of Norway spruce, which demonstrates the value of mining this dataset, accessible in ConGenIE together with advanced visualization tools.

Genetic diversity of Norway spruce ecotypes assessed by GBS-derived SNPs

We investigated the genetic structure of three phenotypically distinct ecotypic groups of Norway spruce (Picea abies) belonging to three elevational classes; namely, low- (acuminata), medium- (europaea), and high-elevation (obovata) form, each represented by 150 trees. After rigorous filtering, we used 1916 Genotyping-by-Sequencing generated SNPs for analysis. Outputs from three multivariate analysis methods (Bayesian clustering algorithm implemented in STRUCTURE, Principal Component Analysis, and the Discriminant Analysis of Principal Components) indicated the presence of a distinct genetic cluster representing the high-elevation ecotypic group. Our findings bring a vital message to forestry practice affirming that artificial transfer of forest reproductive material, especially for stands under harsh climate conditions, should be considered with caution.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

Genome-Wide Identification, Characterization and Expression Analysis of Soybean CHYR Gene Family

The CHYR (CHY ZINC-FINGER AND RING FINGER PROTEIN) proteins have been functionally characterized in iron regulation and stress response in Arabidopsis, rice and Populus. However, their roles in soybean have not yet been systematically investigated. Here, in this study, 16 GmCHYR genes with conserved Zinc_ribbon, CHY zinc finger and Ring finger domains were obtained and divided into three groups. Moreover, additional 2-3 hemerythrin domains could be found in the N terminus of Group III. Phylogenetic and homology analysis of CHYRs in green plants indicated that three groups might originate from different ancestors. Expectedly, GmCHYR genes shared similar conserved domains/motifs distribution within the same group. Gene expression analysis uncovered their special expression patterns in different soybean tissues/organs and under various abiotic stresses. Group I and II members were mainly involved in salt and alkaline stresses. The expression of Group III members was induced/repressed by dehydration, salt and alkaline stresses, indicating their diverse roles in response to abiotic stress. In conclusion, our work will benefit for further revealing the biological roles of GmCHYRs.

Prior secondary cell wall formation is required for gelatinous layer deposition and posture control in gravi-stimulated aspen

Cell walls, especially secondary cell walls (SCWs), maintain cell shape and reinforce wood, but their structure and shape can be altered in response to gravity. In hardwood trees, tension wood is formed along the upper side of a bending stem and contains wood fiber cells that have a gelatinous layer (G-layer) inside the SCW. In a previous study, we generated nst/snd quadruple-knockout aspens (Populus tremula × Populus tremuloides), in which SCW formation was impaired in 99% of the wood fiber cells. In the present study, we produced nst/snd triple-knockout aspens, in which a large number of wood fibers had thinner SCWs than the wild type (WT) and some had no SCW. Because SCW layers are always formed prior to G-layer deposition, the nst/snd mutants raise interesting questions of whether the mutants can form G-layers without SCW and whether they can control their postures in response to changes in gravitational direction. The nst/snd mutants and the WT plants showed growth eccentricity and vessel frequency reduction when grown on an incline, but the triple mutants recovered their upright growth only slightly, and the quadruple mutants were unable to maintain their postures. The mutants clearly showed that the G-layers were formed in SCW-containing wood fibers but not in those lacking the SCW. Our results indicate that SCWs are essential for G-layer formation and posture control. Furthermore, each wood fiber cell may be able to recognize its cell wall developmental stage to initiate the formation of the G-layer as a response to gravistimulation.

X-ray computed tomography for 3D plant imaging

X-ray computed tomography (CT) is a valuable tool for 3D imaging of plant tissues and organs. Applications include the study of plant development and organ morphogenesis, as well as modeling of transport processes in plants. Some challenges remain, however, including attaining higher contrast for easier quantification, increasing the resolution for imaging subcellular features, and decreasing image acquisition and processing time for high-throughput phenotyping. In addition, phase contrast, multispectral, dark-field, soft X-ray, and time-resolved imaging are emerging. At the same time, a large amount of 3D image data are becoming available, posing challenges for data management. We review recent advances in the area of X-ray CT for plant imaging, and describe opportunities for using such images for studying transport processes in plants.