The genome of black cottonwood, Populus trichocarpa (Torr. & Gray)

Differential prolyl hydroxylation by six Physcomitrella prolyl-4 hydroxylases

Hydroxylation of prolines to 4-trans-hydroxyproline (Hyp) is mediated by prolyl-4 hydroxylases (P4Hs). In plants, Hyps occur in Hydroxyproline-rich glycoproteins (HRGPs), and are frequently O-glycosylated. While both modifications are important, e.g. for cell wall stability, they are undesired in plant-made pharmaceuticals. Sequence motifs for prolyl-hydroxylation were proposed but did not include data from mosses, such as Physcomitrella. We identified six moss P4Hs by phylogenetic reconstruction. Our analysis of 73 Hyps in 24 secretory proteins from multiple mass spectrometry datasets revealed that prolines near other prolines, alanine, serine, threonine and valine were preferentially hydroxylated. About 95 % of Hyps were predictable with combined established methods. In our data, AOV was the most frequent pattern. A combination of 443 AlphaFold models and MS data with 3000 prolines found Hyps mainly on protein surfaces in disordered regions. Moss-produced human erythropoietin (EPO) exhibited O-glycosylation with arabinose chains on two Hyps. This modification was significantly reduced in a p4h1 knock-out (KO) Physcomitrella mutant. Quantitative proteomics with different p4h mutants revealed specific changes in protein amounts, and a modified prolyl-hydroxylation pattern, suggesting a differential function of the Physcomitrella P4Hs. Quantitative RT-PCR revealed a differential effect of single p4h KOs on the expression of the other five p4h genes, suggesting a partial compensation of the mutation. AlphaFold-Multimer models for Physcomitrella P4H1 and its target EPO peptide superposed with the crystal structure of Chlamydomonas P4H1 suggested significant amino acids in the active centre of the enzyme and revealed differences between P4H1 and the other Physcomitrella P4Hs.

C2H2 zinc finger protein PagIDD15A regulates secondary wall thickening and lignin biosynthesis in poplar

Wood production is largely determined by the activity of cambial cell proliferation, and the secondary cell wall (SCW) thickening of xylem cells determines the wood property. In this study, we identified an INDETERMINATE DOMAIN (IDD) type C2H2 zinc finger transcription factor PagIDD15A as a regulator of wood formation in Populus alba × Populus glandulosa. Downregulation of PagIDD15A expression by RNA interference (RNAi) inhibited xylem development and xylem cell secondary wall thickening. RNA-seq analysis showed that PagPAL1, PagCCR2 and PagCCoAOMT1 were downregulated in the differentiating xylem of the PagIDD15A-RNAi transgenic plants, showing that PagIDD15A may regulate SCW biosynthesis through inhibiting lignin biosynthesis. The downregulation of PagVND6-B2, PagMYB10 and PagMYC4 and upregulation of PagWRKY12 in the differentiating xylem of RNAi transgenic plants suggest that PagIDD15A may also regulate these transcription factor (TF) genes to affect SCW thickening. RT-qPCR analysis in the phloem-cambium of RNAi transgenic demonstrates that PagIDD15A may regulate the expression of the genes associated with cell proliferation, including, PagSHR (SHORTROOT), PagSCR (SCARECROW), PagCYCD3;1 (CYCLIN D3;1) and PagSMR4 (SIAMESE-RELATED4), to affect the cambial activity. This study provides the knowledge of the IDD-type C2H2 zinc finger protein in regulating wood formation.

Genome-wide investigation and analysis of NAC transcription factor family in Populus tomentosa and expression analysis under salt stress

The NAC transcription factor family is one of the largest families of TFs in plants, and members of NAC gene family play important roles in plant growth and stress response. Recent release of the haplotype-resolved genome assembly of P. tomentosa provide a platform for NAC protein genome-wide analysis. A total of 270 NAC genes were identified and a comprehensive overview of the PtoNAC gene family is presented, including gene promoter, structure and conserved motif analyses, chromosome localization and collinearity analysis, protein phylogeny, expression pattern, and interaction analysis. The results indicate that protein length, molecular weight, and theoretical isoelectric points of the NAC TF family vary, while gene structure and motif are relatively conserved. Chromosome mapping analysis showed that the P. tomentosa NAC genes are unevenly distributed on 19 chromosomes. The interchromosomal evolutionary results indicate 12 pairs of tandem and 280 segmental duplications. Segmental duplication is possibly related to amplification of P. tomentosa NAC gene family. Expression patterns of 35 PtoNAC genes from P. tomentosa subgroup were analysed under high salinity, and seven NAC genes were induced by this treatment. Promoter and protein interaction network analyses showed that PtoNAC genes are closely associated with growth, development, and abiotic and biotic stress, especially salt stress. These results provide a meaningful reference for follow-up studies of the functional characteristics of NAC genes in the mechanism of stress response and their potential roles in development of P. tomentosa.

Forest genomics in the Caucasus through the lens of its dominant tree species - Fagus orientalis

The last glacial period is known to have greatly influenced the demographic history of temperate forest trees, with important range contractions and post-glacial expansions that led to the formation of multiple genetic lineages and secondary contact zones in the Northern Hemisphere. These dynamics have been extensively studied for European and North American species but are still poorly understood in other temperate regions of rich biodiversity such as the Caucasus. Our study helps filling that gap by deciphering the genomic landscapes of F. orientalis across the South Caucasus. The use of genome-wide data confirmed a past demographic history strongly influenced by the Last Glacial Maximum, revealing two disjunct glacial refugia in the Colchis and Hyrcanian regions. The resulting patterns of genetic diversity, load and differentiation are not always concordant across the region, with genetic load pinpointing the location of the glacial refugia more efficiently than genetic diversity alone. The Hyrcanian forests show depleted genetic diversity and substantial isolation, even if long-distance gene flow is still present with the main centre of diversity in the Greater Caucasus. Finally, we characterize a strong heterogeneity of genetic diversity and differentiation along the species chromosomes, with noticeably a first chromosome showing low diversity and weak differentiation.

Genome-wide characterization of RR gene family members in Zanthoxylum armatum and the subsequent functional characterization of the C-type RR

Response Regulators (RRs) are crucial regulators in plant development and stress responses, comprising A-type, B-type, C-type, and pseudo-RR subfamilies. However, previous studies have often focused on specific subfamilies, which restricts our understanding of the complete RR gene family. In this study, we conducted a comprehensive analysis of 63 RR members from Zanthoxylum armatum, using phylogenetic relationships, motif composition, cis-acting elements, gene duplication and collinearity analyses. Segmental repeats among ZaRR genes enhanced the various environmental adaptabilities of Z. armatum, and the B-type ZaRR exhibited significant collinearity with the RRs in P. trichocarpa and C. sinensis. Cis-element analysis indicated ZaRRs play a significant role in abiotic stress and phytohormone pathways, particularly in light, drought, cold, abscisic acid (ABA) and salicylic acid (SA) responses. Abundant Ethylene Response Factor (ERF) and reproduction-associated binding sites in ZaRR promoters suggested their roles in stress and reproductive processes. A-type ZaRRs were implicated in plant vegetative and reproductive growth, whereas B-type ZaRRs contributed to both growth and stress responses. C-type ZaRRs were associated with plant reproductive growth, whereas pseudo-RRs may function in plant stress responses, such as water logging, cold, and response to ethylene (ETH), SA, and jasmonic acid (JA). Ectopic expression of ZaRR24, a C-type RR, inhibits growth, induces early flowering, and shortens fruit length in Arabidopsis. ZaRR24 overexpression also affected the expression of A- and B-type RRs, as well as floral meristem and organ identity genes. These findings establish a solid and comprehensive foundation for RR gene research in Z. armatum, and provide a platform for investigating signal transduction in other woody plants.

Rare variations within the serine/arginine-rich splicing factor PtoRSZ21 modulate stomatal size to determine drought tolerance in Populus

Rare variants contribute significantly to the 'missing heritability' of quantitative traits. The genome-wide characteristics of rare variants and their roles in environmental adaptation of woody plants remain unexplored. Utilizing genome-wide rare variant association study (RVAS), expression quantitative trait loci (eQTL) mapping, genetic transformation, and molecular experiments, we explored the impact of rare variants on stomatal morphology and drought adaptation in Populus. Through comparative analysis of five world-wide Populus species, we observed the influence of mutational bias and adaptive selection on the distribution of rare variants. RVAS identified 75 candidate genes correlated with stomatal size (SS)/stomatal density (SD), and a rare haplotype in the promoter of serine/arginine-rich splicing factor PtoRSZ21 emerged as the foremost association signal governing SS. As a positive regulator of drought tolerance, PtoRSZ21 can recruit the core splicing factor PtoU1-70K to regulate alternative splicing (AS) of PtoATG2b (autophagy-related 2). The rare haplotype PtoRSZ21hap2 weakens binding affinity to PtoMYB61, consequently affecting PtoRSZ21 expression and SS, ultimately resulting in differential distribution of Populus accessions in arid and humid climates. This study enhances the understanding of regulatory mechanisms that underlie AS induced by rare variants and might provide targets for drought-tolerant varieties breeding in Populus.

Genome-wide and functional analysis of late embryogenesis abundant (LEA) genes during dormancy and sprouting periods of kernel consumption apricots (P. armeniaca L. × P. sibirica L.)

Late embryogenesis abundant (LEA) proteins play a crucial role in protecting cells from stress, making them potential contributors to abiotic stress tolerance. This study focuses on apricot (P. armeniaca L. × P. sibirica L.), where a comprehensive genome-wide analysis identified 54 LEA genes, categorized into eight subgroups based on phylogenetic relationships. Synteny analysis revealed 14 collinear blocks containing LEA genes between P. armeniaca × P. sibirica and Arabidopsis thaliana, with an additional 9 collinear blocks identified between P. armeniaca × P. sibirica and poplar. Examination of gene structure and conserved motifs indicated that these subgroups exhibit consistent exon-intron patterns and shared motifs. The expansion and duplication of LEA genes in P. armeniaca × P. sibirica were driven by whole-genome duplication (WGD), segmental duplication, and tandem duplication events. Expression analysis, utilizing RNA-seq data and quantitative real-time RT-PCR (qRT-PCR), indicated induction of PasLEA2-20, PasLEA3-2, PasLEA6-1, Pasdehydrin-3, and Pasdehydrin-5 in flower buds during dormancy and sprouting phases. Coexpression network analysis linked LEA genes with 15 cold-resistance genes. Remarkably, during the four developmental stages of flower buds in P. armeniaca × P. sibirica - physiological dormancy, ecological dormancy, sprouting period, and germination stage - the expression patterns of all PasLEAs coexpressed with cold stress-related genes remained consistent. Protein-protein interaction networks, established using Arabidopsis orthologs, emphasized connections between PasLEA proteins and cold resistance pathways. Overexpression of certain LEA genes in yeast and Arabidopsis conferred advantages under cold stress, including increased pod length, reduced bolting time and flowering time, improved survival and seed setting rates, elevated proline accumulation, and enhanced antioxidative enzymatic activities. Furthermore, these overexpressed plants exhibited upregulation of genes related to flower development and cold resistance. The Y1H assay confirmed that PasGBF4 and PasDOF3.5 act as upstream regulatory factors by binding to the promoter region of PasLEA3-2. PasDOF2.4, PasDnaJ2, and PasAP2 were also found to bind to the promoter of Pasdehydrin-3, regulating the expression levels of downstream genes. This comprehensive study explores the evolutionary relationships among PasLEA genes, protein interactions, and functional analyses during various stages of dormancy and sprouting in P. armeniaca × P. sibirica. It offers potential targets for enhancing cold resistance and manipulating flower bud dormancy in this apricot hybrid.

Evolutionary divergence of CXE gene family in green plants unveils that PtoCXEs overexpression reduces fungal colonization in transgenic Populus

Plant enzymes significantly contribute to the rapidly diversified metabolic repertoire since the colonization of land by plants. Carboxylesterase is just one of the ubiquitous, multifunctional and ancient enzymes that has particularly diversified during plant evolution. This study provided a status on the carboxylesterase landscape within Viridiplantae. A total of 784 carboxylesterases were identified from the genome of 31 plant species representing nine major lineages of sequenced Viridiplantae and divided into five clades based on phylogenetic analysis. Clade I carboxylesterase genes may be of bacterial origin and then expanded and diversified during plant evolution. Clade II was first gained in the ancestor of bryophytes after colonization of land by plants, Clade III and Clade IV in ferns which were considered the most advanced seedless vascular plants, while Clade V was gained in seed plants. To date, the functions of carboxylesterase genes in woody plants remain unclear. In this study, 51 carboxylesterase genes were identified from the genome of Populus trichocarpa and further divided into eight classes. Tandem and segmental duplication events both contributed to the expansion of carboxylesterase genes in Populus. Although carboxylesterase genes were proven to enhance resistance to pathogens in many herbaceous species, relevant researches on forest trees are still needed. In this study, pathogen incubation assays showed that overexpressing of six Class VI carboxylesterases in Populus tomentosa, to a greater or lesser degree, reduced colonization of detached leaves by fungus Cytospora chrysosperma. A significant difference was also found in functional divergence patterns for genes derived from different gene duplication events. Functional differentiation of duplicated carboxylesterase genes in Populus was proved for the first time by in vivo physiological analysis. The identification of the potentially anti-fungal PtoCXE06 gene also laid a theoretical foundation for promoting the genetic improvement of disease-resistance traits in forest trees.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

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.

Enabling Lignin Valorization Through Integrated Advances in Plant Biology and Biorefining

Despite lignin having long been viewed as an impediment to the processing of biomass for the production of paper, biofuels, and high-value chemicals, the valorization of lignin to fuels, chemicals, and materials is now clearly recognized as a critical element for the lignocellulosic bioeconomy. However, the intended application for lignin will likely require a preferred lignin composition and form. To that end, effective lignin valorization will require the integration of plant biology, providing optimal feedstocks, with chemical process engineering, providing efficient lignin transformations. Recent advances in our understanding of lignin biosynthesis have shown that lignin structure is extremely diverse and potentially tunable, while simultaneous developments in lignin refining have resulted in the development of several processes that are more agnostic to lignin composition. Here, we review the interface between in planta lignin design and lignin processing and discuss the advances necessary for lignin valorization to become a feature of advanced biorefining.

The fnr-like mutants confer isoxaben tolerance by initiating mitochondrial retrograde signalling

Isoxaben is a pre-emergent herbicide used to control broadleaf weeds. While the phytotoxic mechanism is not completely understood, isoxaben interferes with cellulose synthesis. Certain mutations in cellulose synthase complex proteins can confer isoxaben tolerance; however, these mutations can cause compromised cellulose synthesis and perturbed plant growth, rendering them unsuitable as herbicide tolerance traits. We conducted a genetic screen to identify new genes associated with isoxaben tolerance by screening a selection of Arabidopsis thaliana T-DNA mutants. We found that mutations in a FERREDOXIN-NADP(+) OXIDOREDUCTASE-LIKE (FNRL) gene enhanced tolerance to isoxaben, exhibited as a reduction in primary root stunting, reactive oxygen species accumulation and ectopic lignification. The fnrl mutant did not exhibit a reduction in cellulose levels following exposure to isoxaben, indicating that FNRL operates upstream of isoxaben-induced cellulose inhibition. In line with these results, transcriptomic analysis revealed a highly reduced response to isoxaben treatment in fnrl mutant roots. The fnrl mutants displayed constitutively induced mitochondrial retrograde signalling, and the observed isoxaben tolerance is partially dependent on the transcription factor ANAC017, a key regulator of mitochondrial retrograde signalling. Moreover, FNRL is highly conserved across all plant lineages, implying conservation of its function. Notably, fnrl mutants did not show a growth penalty in shoots, making FNRL a promising target for biotechnological applications in breeding isoxaben tolerance in crops.

RNA m6A detection using raw current signals and basecalling errors from Nanopore direct RNA sequencing reads

MOTIVATION

Nanopore direct RNA sequencing (DRS) enables the detection of RNA N6-methyladenosine (m6A) without extra laboratory techniques. A number of supervised or comparative approaches have been developed to identify m6A from Nanopore DRS reads. However, existing methods typically utilize either statistical features of the current signals or basecalling-error features, ignoring the richer information of the raw signals of DRS reads.

RESULTS

Here, we propose RedNano, a deep-learning method designed to detect m6A from Nanopore DRS reads by utilizing both raw signals and basecalling errors. RedNano processes the raw-signal feature and basecalling-error feature through residual networks. We validated the effectiveness of RedNano using synthesized, Arabidopsis, and human DRS data. The results demonstrate that RedNano surpasses existing methods by achieving higher area under the ROC curve (AUC) and area under the precision-recall curve (AUPRs) in all three datasets. Furthermore, RedNano performs better in cross-species validation, demonstrating its robustness. Additionally, when detecting m6A from an independent dataset of Populus trichocarpa, RedNano achieves the highest AUC and AUPR, which are 3.8%-9.9% and 5.5%-13.8% higher than other methods, respectively.

AVAILABILITY AND IMPLEMENTATION

The source code of RedNano is freely available at https://github.com/Derryxu/RedNano.

Genome-Wide Identification of TLP Gene Family in Populus trichocarpa and Functional Characterization of PtTLP6, Preferentially Expressed in Phloem

Thaumatin-like proteins (TLPs) in plants are involved in diverse biotic and abiotic stresses, including antifungal activity, low temperature, drought, and high salinity. However, the roles of the TLP genes are rarely reported in early flowering. Here, the TLP gene family was identified in P. trichocarpa. The 49 PtTLP genes were classified into 10 clusters, and gene structures, conserved motifs, and expression patterns were analyzed in these PtTLP genes. Among 49 PtTLP genes, the PtTLP6 transcription level is preferentially high in stems, and GUS staining signals were mainly detected in the phloem tissues of the PtTLP6pro::GUS transgenic poplars. We generated transgenic Arabidopsis plants overexpressing the PtTLP6 gene, and its overexpression lines showed early flowering phenotypes. However, the expression levels of main flowering regulating genes were not significantly altered in these PtTLP6-overexpressing plants. Our data further showed that overexpression of the PtTLP6 gene led to a reactive oxygen species (ROS) burst in Arabidopsis, which might advance the development process of transgenic plants. In addition, subcellular localization of PtTLP6-fused green fluorescent protein (GFP) was in peroxisome, as suggested by tobacco leaf transient transformation. Overall, this work provides a comprehensive analysis of the TLP gene family in Populus and an insight into the role of TLPs in woody plants.

Nanoplastic toxicity induces metabolic shifts in Populus × euramericana cv. '74/76' revealed by multi-omics analysis

There is increasing global concern regarding the pervasive issue of plastic pollution. We investigated the response of Populus × euramericana cv. '74/76' to nanoplastic toxicity via phenotypic, microanatomical, physiological, transcriptomic, and metabolomic approaches. Polystyrene nanoplastics (PS-NPs) were distributed throughout the test plants after the application of PS-NPs. Nanoplastics principally accumulated in the roots; minimal fractions were translocated to the leaves. In leaves, however, PS-NPs easily penetrated membranes and became concentrated in chloroplasts, causing thylakoid disintegration and chlorophyll degradation. Finally, oxidant damage from the influx of PS-NPs led to diminished photosynthesis, stunted growth, and etiolation and/or wilting. By integrating dual-omics data, we found that plants could counteract mild PS-NP-induced oxidative stress through the antioxidant enzyme system without initiating secondary metabolic defense mechanisms. In contrast, severe PS-NP treatments promoted a shift in metabolic pattern from primary metabolism to secondary metabolic defense mechanisms, an effect that was particularly pronounced during the upregulation of flavonoid biosynthesis. Our findings provide a useful framework from which to further clarify the roles of key biochemical pathways in plant responses to nanoplastic toxicity. Our work also supports the development of effective strategies to mitigate the environmental risks of nanoplastics by biologically immobilizing them in contaminated lands.

Genome-Wide Analysis of C/S1-bZIP Subfamilies in Populus tomentosa and Unraveling the Role of PtobZIP55/21 in Response to Low Energy

C/S1 basic leucine zipper (bZIP) transcription factors are essential for plant survival under energy deficiency. However, studies on the responses of C/S1-bZIPs to low energy in woody plants have not yet been reported. In this study, members of C/S1-bZIP subfamilies in Populus tomentosa were systematically analyzed using bioinformatic approaches. Four C-bZIPs and 10 S1-bZIPs were identified, and their protein properties, phylogenetic relationships, gene structures, conserved motifs, and uORFs were systematically investigated. In yeast two-hybrid assays, direct physical interactions between C-bZIP and S1-bZIP members were observed, highlighting their potential functional synergy. Moreover, expression profile analyses revealed that low energy induced transcription levels of most C/S1-bZIP members, with bZIP55 and bZIP21 (a homolog of bZIP55) exhibiting particularly significant upregulation. When the expression of bZIP55 and bZIP21 was co-suppressed using artificial microRNA mediated gene silencing in transgenic poplars, root growth was promoted. Further analyses revealed that bZIP55/21 negatively regulated the root development of P. tomentosa in response to low energy. These findings provide insights into the molecular mechanisms by which C/S1-bZIPs regulate poplar growth and development in response to energy deprivation.

Assembly and analysis of the genome of Notholithocarpus densiflorus

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to "sudden oak death" (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Genetic underpinnings of arthropod community distributions in Populus trichocarpa

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Transcription factor PagMYB31 positively regulates cambium activity and negatively regulates xylem development in poplar

Wood formation involves consecutive developmental steps, including cell division of vascular cambium, xylem cell expansion, secondary cell wall (SCW) deposition, and programmed cell death. In this study, we identified PagMYB31 as a coordinator regulating these processes in Populus alba × Populus glandulosa and built a PagMYB31-mediated transcriptional regulatory network. PagMYB31 mutation caused fewer layers of cambial cells, larger fusiform initials, ray initials, vessels, fiber and ray cells, and enhanced xylem cell SCW thickening, showing that PagMYB31 positively regulates cambial cell proliferation and negatively regulates xylem cell expansion and SCW biosynthesis. PagMYB31 repressed xylem cell expansion and SCW thickening through directly inhibiting wall-modifying enzyme genes and the transcription factor genes that activate the whole SCW biosynthetic program, respectively. In cambium, PagMYB31 could promote cambial activity through TRACHEARY ELEMENT DIFFERENTIATION INHIBITORY FACTOR (TDIF)/PHLOEM INTERCALATED WITH XYLEM (PXY) signaling by directly regulating CLAVATA3/ESR-RELATED (CLE) genes, and it could also directly activate WUSCHEL HOMEOBOX RELATED4 (PagWOX4), forming a feedforward regulation. We also observed that PagMYB31 could either promote cell proliferation through the MYB31-MYB72-WOX4 module or inhibit cambial activity through the MYB31-MYB72-VASCULAR CAMBIUM-RELATED MADS2 (VCM2)/PIN-FORMED5 (PIN5) modules, suggesting its role in maintaining the homeostasis of vascular cambium. PagMYB31 could be a potential target to manipulate different developmental stages of wood formation.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

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).

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

Genome-Wide Identification and Analysis of the EIN3/EIL Transcription Factor Gene Family in Doubled Haploid (DH) Poplar

Ethylene (ET) is an important phytohormone that regulates plant growth, development and stress responses. The ethylene-insensitive3/ethylene-insensitive3-like (EIN3/EIL) transcription factor family, as a key regulator of the ET signal transduction pathway, plays an important role in regulating the expression of ET-responsive genes. Although studies of EIN3/EIL family members have been completed in many species, their role in doubled haploid (DH) poplar derived from another culture of diploid Populus simonii × P. nigra (donor tree, DT) remains ambiguous. In this study, a total of seven EIN3/EIL gene family members in the DH poplar genome were identified. Basic physical and chemical property analyses of these genes were performed, and these proteins were predicted to be localized to the nucleus. According to the phylogenetic relationship, EIN3/EIL genes were divided into two groups, and the genes in the same group had a similar gene structure and conserved motifs. The expression patterns of EIN3/EIL genes in the apical buds of different DH poplar plants were analyzed based on transcriptome data. At the same time, the expression patterns of PsnEIL1, PsnEIN3, PsnEIL4 and PsnEIL5 genes in different tissues of different DH plants were detected via RT-qPCR, including the apical buds, young leaves, functional leaves, xylem, cambium and roots. The findings presented above indicate notable variations in the expression levels of PsnEIL genes across various tissues of distinct DH plants. Finally, the PsnEIL1 gene was overexpressed in DT, and the transgenic plants showed a dwarf phenotype, indicating that the PsnEIL1 gene was involved in regulating the growth and development of poplar. In this study, the EIN3/EIL gene family of DH poplar was analyzed and functionally characterized, which provides a theoretical basis for the future exploration of the EIN3/EIL gene function.

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.

A multi-omic survey of black cottonwood tissues highlights coordinated transcriptomic and metabolomic mechanisms for plant adaptation to phosphorus deficiency

Introduction

Phosphorus (P) deficiency in plants creates a variety of metabolic perturbations that decrease photosynthesis and growth. Phosphorus deficiency is especially challenging for the production of bioenergy feedstock plantation species, such as poplars (Populus spp.), where fertilization may not be practically or economically feasible. While the phenotypic effects of P deficiency are well known, the molecular mechanisms underlying whole-plant and tissue-specific responses to P deficiency, and in particular the responses of commercially valuable hardwoods, are less studied.

Methods

We used a multi-tissue and multi-omics approach using transcriptomic, proteomic, and metabolomic analyses of the leaves and roots of black cottonwood (Populus trichocarpa) seedlings grown under P-deficient (5 µM P) and replete (100 µM P) conditions to assess this knowledge gap and to identify potential gene targets for selection for P efficiency.

Results

In comparison to seedlings grown at 100 µM P, P-deficient seedlings exhibited reduced dry biomass, altered chlorophyll fluorescence, and reduced tissue P concentrations. In line with these observations, growth, C metabolism, and photosynthesis pathways were downregulated in the transcriptome of the P-deficient plants. Additionally, we found evidence of strong lipid remodeling in the leaves. Metabolomic data showed that the roots of P-deficient plants had a greater relative abundance of phosphate ion, which may reflect extensive degradation of P-rich metabolites in plants exposed to long-term P-deficiency. With the notable exception of the KEGG pathway for Starch and Sucrose Metabolism (map00500), the responses of the transcriptome and the metabolome to P deficiency were consistent with one another. No significant changes in the proteome were detected in response to P deficiency.

Discussion and conclusion

Collectively, our multi-omic and multi-tissue approach enabled the identification of important metabolic and regulatory pathways regulated across tissues at the molecular level that will be important avenues to further evaluate for P efficiency. These included stress-mediating systems associated with reactive oxygen species maintenance, lipid remodeling within tissues, and systems involved in P scavenging from the rhizosphere.

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.

Genome-wide analysis of the mulberry (Morus abla L.) GH9 gene family and the functional characterization of MaGH9B6 during the development of the abscission zone

Plant glycoside hydrolase family 9 genes (GH9s) are widely distributed in plants and involved in a variety of cellular and physiological processes. In the current study, nine GH9 genes were identified in the mulberry and were divided into two subfamilies based on the phylogenetic analysis. Conserved motifs and gene structure analysis suggested that the evolution of the two subfamilies is relatively conserved and the glycoside hydrolase domain almost occupy the entire coding region of the GH9s gene. Only segmental duplication has played a role in the expansion of gene family. Collinearity analysis showed that mulberry GH9s had the closest relationship with poplar GH9s. MaGH9B1, MaGH9B6, MaGH9B5, and MaGH9B3 were detected to have transcript accumulation in the stalk of easy-to drop mature fruit drop, suggesting that these could play a role in mulberry fruit drop. Multiple cis-acting elements related to plant hormones and abiotic stress responses were found in the mulberry GH9 promoter regions and showed different activities under exogenous abscisic acid (ABA) and 2,4- dichlorophenoxyacetic acid (2,4-D) stresses. We found that the lignin content in the fruit stalk decreased with the formation of the abscission zone (AZ), which could indirectly reflect the formation process of the AZ. These results provide a theoretical basis for further research on the role of GH9s in mulberry abscission.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Comparative analysis of PLATZ transcription factors in six poplar species and analysis of the role of PtrPLATZ14 in leaf development

Plant AT-rich sequence and zinc-binding (PLATZ) proteins are a class of plant-specific transcription factor that play a crucial role in plant growth, development, and stress response. However, the evolutionary relationship of the PLATZ gene family across the Populus genus and the biological functions of the PLATZ protein require further investigation. In this study, we identified 133 PLATZ genes from six Populus species belonging to four Populus sections. Synteny analysis of the PLATZ gene family indicated that whole genome duplication events contributed to the expansion of the PLATZ family. Among the nine paralogous pairs, the protein structure of PtrPLATZ14/18 pair exhibited significant differences with others. Through gene expression patterns and co-expression networks, we discovered divergent expression patterns and sub-networks, and found that the members of pair PtrPLATZ14/18 might play different roles in the regulation of macromolecule biosynthesis and modification. Furthermore, we found that PtrPLATZ14 regulates poplar leaf development by affecting cell size control genes PtrGRF/GIF and PtrTCP. In conclusion, our study provides a theoretical foundation for exploring the evolution relationships and functions of the PLATZ gene family within Populus species and provides insights into the function and potential mechanism of PtrPLATZ14 in leaf morphology that were diverse across the Populus genus.

Suberin deficiency and its effect on the transport physiology of young poplar roots

The precise functions of suberized apoplastic barriers in root water and nutrient transport physiology have not fully been elucidated. While lots of research has been performed with mutants of Arabidopsis, little to no data are available for mutants of agricultural crop or tree species. By employing a combined set of physiological, histochemical, analytical, and transport physiological methods as well as RNA-sequencing, this study investigated the implications of remarkable CRISPR/Cas9-induced suberization defects in young roots of the economically important gray poplar. While barely affecting overall plant development, contrary to literature-based expectations significant root suberin reductions of up to 80-95% in four independent mutants were shown to not evidently affect the root hydraulic conductivity during non-stress conditions. In addition, subliminal iron deficiency symptoms and increased translocation of a photosynthesis inhibitor as well as NaCl highlight the involvement of suberin in nutrient transport physiology. The multifaceted nature of the root hydraulic conductivity does not allow drawing simplified conclusions such as that the suberin amount must always be correlated with the water transport properties of roots. However, the decreased masking of plasma membrane surface area could facilitate the uptake but also leakage of beneficial and harmful solutes.

Genome-Wide Identification and Expression Profiling of Velvet Complex Transcription Factors in Populus alba × Populus glandulosa

The discovery of new genes with novel functions is a major driver of adaptive evolutionary innovation in plants. Especially in woody plants, due to genome expansion, new genes evolve to regulate the processes of growth and development. In this study, we characterized the unique VeA transcription factor family in Populus alba × Populus glandulosa, which is associated with secondary metabolism. Twenty VeA genes were characterized systematically on their phylogeny, genomic distribution, gene structure and conserved motif, promoter binding site, and expression profiling. Furthermore, through ChIP-qPCR, Y1H, and effector-reporter assays, it was demonstrated that PagMYB128 directly regulated PagVeA3 to influence the biosynthesis of secondary metabolites. These results provide a basis for further elucidating the function of VeAs gene in poplar and its genetic regulation mechanism.

Comparative genomic analysis of the Growth-Regulating Factors-Interacting Factors (GIFs) in six Salicaceae species and functional analysis of PeGIF3 reveals their regulatory role in Populus heteromorphic leaves

BACKGROUND

The growth-regulating factor-interacting factor (GIF) gene family plays a vital role in regulating plant growth and development, particularly in controlling leaf, seed, and root meristem homeostasis. However, the regulatory mechanism of heteromorphic leaves by GIF genes in Populus euphratica as an important adaptative trait of heteromorphic leaves in response to desert environment remains unknown.

RESULTS

This study aimed to identify and characterize the GIF genes in P. euphratica and other five Salicaceae species to investigate their role in regulating heteromorphic leaf development. A total of 27 GIF genes were identified and characterized across six Salicaceae species (P. euphratica, Populus pruinose, Populus deltoides, Populus trichocarpa, Salix sinopurpurea, and Salix suchowensis) at the genome-wide level. Comparative genomic analysis among these species suggested that the expansion of GIFs may be derived from the specific Salicaceae whole-genome duplication event after their divergence from Arabidopsis thaliana. Furthermore, the expression data of PeGIFs in heteromorphic leaves, combined with functional information on GIF genes in Arabidopsis, indicated the role of PeGIFs in regulating the leaf development of P. euphratica, especially PeGIFs containing several cis-acting elements associated with plant growth and development. By heterologous expression of the PeGIF3 gene in wild-type plants (Col-0) and atgif1 mutant of A. thaliana, a significant difference in leaf expansion along the medial-lateral axis, and an increased number of leaf cells, were observed between the overexpressed plants and the wild type.

CONCLUSION

PeGIF3 enhances leaf cell proliferation, thereby resulting in the expansion of the central-lateral region of the leaf. The findings not only provide global insights into the evolutionary features of Salicaceae GIFs but also reveal the regulatory mechanism of PeGIF3 in heteromorphic leaves of P. euphratica.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

Transcriptomic Analysis for Diurnal Temperature Differences Reveals Gene-Regulation-Network Response to Accumulation of Bioactive Ingredients of Protocorm-like Bodies in Dendrobium officinale

Dendrobium officinale Kimura et Migo (D. officinale) is one of the most important traditional Chinese medicinal herbs, celebrated for its abundant bioactive ingredients. This study demonstrated that the diurnal temperature difference (DIF) (T1: 13/13 °C, T2: 25/13 °C, and T3: 25/25 °C) was more favorable for high chlorophyll, increased polysaccharide, and total flavonoid contents compared to constant temperature treatments in D. officinale PLBs. The transcriptome analysis revealed 4251, 4404, and 4536 differentially expressed genes (DEGs) in three different comparisons (A: 25/13 °C vs. 13/13 °C, B: 13/13 °C vs. 25/25 °C, and C: 25/13 °C vs. 25/25 °C, respectively). The corresponding up-/down-regulated DEGs were 1562/2689, 2825/1579, and 2310/2226, respectively. GO and KEGG enrichment analyses of DEGs showed that the pathways of biosynthesis of secondary metabolites, carotenoid biosynthesis, and flavonoid biosynthesis were enriched in the top 20; further analysis of the sugar- and flavonol-metabolism pathways in D. officinale PLBs revealed that the DIF led to a differential gene expression in the enzymes linked to sugar metabolism, as well as to flavonol metabolism. Certain key metabolic genes related to ingredient accumulation were identified, including those involved in polysaccharide metabolism (SUS, SUT, HKL1, HGL, AMY1, and SS3) and flavonol (UGT73C and UGT73D) metabolism. Therefore, these findings indicated that these genes may play an important role in the regulatory network of the DIF in the functional metabolites of D. officinale PLBs. In a MapMan annotation of abiotic stress pathways, the DEGs with significant changes in their expression levels were mainly concentrated in the heat-stress pathways, including heat-shock proteins (HSPs) and heat-shock transcription factors (HSFs). In particular, the expression levels of HSP18.2, HSP70, and HSF1 were significantly increased under DIF treatment, which suggested that HSF1, HSP70 and HSP18.2 may respond to the DIF. In addition, they can be used as candidate genes to study the effect of the DIF on the PLBs of D. officinale. The results of our qPCR analysis are consistent with those of the transcriptome-expression analysis, indicating the reliability of the sequencing. The results of this study revealed the transcriptome mechanism of the DIF on the accumulation of the functional metabolic components of D. officinale. Furthermore, they also provide an important theoretical basis for improving the quality of D. officinale via the DIF in production.

Analysis of the UDP-Glucosyltransferase (UGT) Gene Family and Its Functional Involvement in Drought and Salt Stress Tolerance in Phoebe bournei

Uridine diphosphate glycosyltransferases (UDP-GTs, UGTs), which are regulated by UGT genes, play a crucial role in glycosylation. In vivo, the activity of UGT genes can affect the availability of metabolites and the rate at which they can be eliminated from the body. UGT genes can exert their regulatory effects through mechanisms such as post-transcriptional modification, substrate subtype specificity, and drug interactions. Phoebe bournei is an economically significant tree species that is endemic to southern China. Despite extensive studies on the UGT gene family in various species, a comprehensive investigation of the UGT family in P. bournei has not been reported. Therefore, we conducted a systematic analysis to identify 156 UGT genes within the entire P. bournei genome, all of which contained the PSPG box. The PbUGT family consists of 14 subfamilies, consistent with Arabidopsis thaliana. We observed varying expression levels of PbUGT genes across different tissues in P. bournei, with the following average expression hierarchy: leaf > stem xylem > stem bark > root xylem > root bark. Covariance analysis revealed stronger covariance between P. bournei and closely related species. In addition, we stressed the seedlings with 10% NaCl and 10% PEG-6000. The PbUGT genes exhibited differential expression under drought and salt stresses, with specific expression patterns observed under each stress condition. Our findings shed light on the transcriptional response of PbUGT factors to drought and salt stresses, thereby establishing a foundation for future investigations into the role of PbUGT transcription factors.

Comprehensive Evolutionary Analysis of the SMXL Gene Family in Rosaceae: Further Insights into Its Origin, Expansion, Diversification, and Role in Regulating Pear Branching

SMXL genes constitute a conserved gene family that is ubiquitous in angiosperms and involved in regulating various plant processes, including branching, leaf elongation, and anthocyanin biosynthesis, but little is known about their molecular functions in pear branching. Here, we performed genome-wide identification and investigation of the SMXL genes in 16 angiosperms and analyzed their phylogenetics, structural features, conserved motifs, and expression patterns. In total, 121 SMXLs genes were identified and were classified into four groups. The number of non-redundant SMXL genes in each species varied from 3 (Amborella trichopoda Baill.) to 18 (Glycine max Merr.) and revealed clear gene expansion events over evolutionary history. All the SMXL genes showed conserved structures, containing no more than two introns. Three-dimensional protein structure prediction revealed distinct structures between but similar structures within groups. A quantitative real-time PCR analysis revealed different expressions of 10 SMXL genes from pear branching induced by fruit-thinning treatment. Overall, our study provides a comprehensive investigation of SMXL genes in the Rosaceae family, especially pear. The results offer a reference for understanding the evolutionary history of SMXL genes and provide excellent candidates for studying fruit tree branching regulation, and in facilitating pear pruning and planting strategies.

Host-Specific and Homologous Pairs of Melampsora larici-populina Effectors Unveil Novel Nicotiana benthamiana Stromule Induction Factors

The poplar rust fungus Melampsora larici-populina is part of one of the most devastating group of fungi (Pucciniales) and causes important economic losses to the poplar industry. Because M. larici-populina is a heteroecious obligate biotroph, its spread depends on its ability to carry out its reproductive cycle through larch and then poplar parasitism. Genomic approaches have identified more than 1,000 candidate secreted effector proteins (CSEPs) from the predicted secretome of M. larici-populina that are potentially implicated in the infection process. In this study, we selected CSEP pairs (and one triplet) among CSEP gene families that share high sequence homology but display specific gene expression profiles among the two distinct hosts. We determined their subcellular localization by confocal microscopy through expression in the heterologous plant system Nicotiana benthamiana. Five out of nine showed partial or complete chloroplastic localization. We also screened for potential protein interactors from larch and poplar by yeast two-hybrid assays. One pair of CSEPs and the triplet shared common interactors, whereas the members of the two other pairs did not have common targets from either host. Finally, stromule induction quantification revealed that two pairs and the triplet of CSEPs induced stromules when transiently expressed in N. benthamiana. The use of N. benthamiana eds1 and nrg1 knockout lines showed that CSEPs can induce stromules through an eds1-independent mechanism. However, CSEP homologs shared the same impact on stromule induction and contributed to discovering a new stromule induction cascade that can be partially and/or fully independent of eds1. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome-Wide Characterization and Expression Analyses of Major Latex Protein Gene Family in Populus simonii × P. nigra

Major latex proteins, or MLPs, are crucial to plants' capacity to grow, develop, and endure biotic and abiotic stresses. The MLP gene family has been found in numerous plants, but little is known about its role in Populus simonii × P. nigra. This study discovered and assessed 43 PtMLP genes that were unevenly dispersed throughout 12 chromosomes in terms of their physicochemical characteristics, gene structure, conserved motifs, and protein localization. Based on their phylogeny and protein structural characteristics, three separate subclasses of PtMLP family were identified. Segmental and tandem duplication were found to be essential variables in the expansion of the PtMLP genes. The involvement of the PtMLP genes in growth and development, as well as in the responses to different hormones and stresses, was demonstrated by cis-regulatory element prediction. The PtMLP genes showed varying expression patterns in various tissues and under different conditions (cold, salt, and drought stress), as demonstrated in RNA-Seq databases, suggesting that PsnMLP may have different functions. Following the further investigation of the genes demonstrating notable variations in expression before and after the application of three stresses, PsnMLP5 was identified as a candidate gene. Subsequent studies revealed that PsnMLP5 could be induced by ABA treatment. This study paves the way for further investigations into the MLP genes' functional mechanisms in response to abiotic stressors, as well as the ways in which they can be utilized in poplar breeding for improved stress tolerance.

Detection of Hybrids in Willows (Salix, Salicaceae) Using Genome-Wide DArTseq Markers

The genus Salix, comprising some 400-500 species, is important in various alluvial or wet habitats of the northern hemisphere. It is a promising crop for applications such as biomass production, biofuels, or environmental projects. Clear species delimitation is crucial in ecology, biotechnology, and horticulture. DArTseq markers, a genome-wide technique, were tested for species and hybrid identification. A total of 179 willow samples were analysed, including six species of Salix subgen. Salix and four species of Salix subgen. Vetrix, including those used in biomass crop production, representing important European taxa. Identification of species-specific markers, clustering analyses (principal coordinate analysis, neighbor-joining) and Bayesian methods (Structure) unambiguously identified putative hybrids. In addition to demonstrating the high efficiency of DArT-seq markers in identifying willow hybrids, we also opened-up new questions about hybridisation processes and systematics. We detected unidirectional hybridisation between S. alba and S. fragilis, forming backcross hybrids, and we rejected the hypothesis that S. fragilis does not occur naturally in Europe. Further, the isolated position of Salix triandra within the genus was confirmed.

Characterization of the complete chloroplast genome of 'Quanhong poplar' (Populus deltoides W. Bartram ex Humphry Marshall, 2011)

The color of the leaves is one of the most important factors for horticultural crops that are considered by breeders, and is also attracting more and more attention from economists and academics. 'Quanhong poplar' (QHP), a rare, bright reddish-purple color-leaf cultivar that has been widely cultivated in China as a landscape tree, is a very precious color-leaf cultivar. In the present study, a reference-based assembly was performed using whole-genome sequencing data to characterize the chloroplast genome of 'QHP'. The total chloroplast genome size of 'QHP' is 156,950 bp, which is divided into two inverted repeat structures of 27,649 bp each, a small single-copy region of 16,563 bp, and a large single-copy region (LSC) of 85,089 bp. From the chloroplast genome, 130 genes have been predicted, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. A chloroplast genome containing 36.68% GC content was detected in 'QHP'. Three SNP sites have been developed between 'QHP' and Populus deltoides Zhonglin 2025. Based on the phylogenetic analysis of chloroplast genomes reported for Populus, the chloroplast of 'QHP' is closest to several strains of Populus deltoides.

The Genome of Vitis zhejiang-adstricta Strengthens the Protection and Utilization of the Endangered Ancient Grape Endemic to China

Vitis zhejiang-adstricta (V. zhejiang-adstricta) is one of the most important and endangered wild grapes. It is a national key protected wild, rare and endangered ancient grape endemic to China and used as a candidate material for resistance breeding owing to its excellent significant disease resistance. Here, we present a high-quality chromosome-level assembly of V. zhejiang-adstricta (IB-VB-01), comprising 506.66 Mb assembled into 19 pseudo-chromosomes. The contig N50 length is 3.91 Mb with 31,196 annotated protein-coding genes. Comparative genome and evolutionary analyses illustrated that V. zhejiang-adstricta has a specific position in the evolution of East Asian Vitis and shared a common ancestor with Vitis vinifera during the divergence of the two species about 10.42 (between 9.34 and 11.12) Mya. The expanded gene families compared with those in plants were related to disease resistance, and constructed gene families were related to plant growth and primary metabolism. With the analysis of gene family expansion and contraction, the evolution of environmental adaptability and especially the NBS-LRR gene family of V. zhejiang-adstricta was elucidated based on the pathways of resistance genes (R genes), unique genes and structural variations. The near-complete and accurate diploid V. zhejiang-adstricta reference genome obtained herein serves as an important complement to wild grape genomes and will provide valuable genomic resources for investigating the genomic architecture of V. zhejiang-adstricta as well as for improving disease resistance breeding strategies in grape.

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

Phylogenomics reveal Populusgonggaensis as a hybrid between P.lasiocarpa and P.cathayana (Salicaceae)

High levels of intra-specific polymorphism and frequent hybridisation make it difficult to define species and correctly apply their scientific names. Populus L. is a challenging genus with plentiful natural and artificial hybrids. This study is a part of the project 'Flora of Pan-Himalaya' and aims to determine the taxonomic identity of P.gonggaensis N. Chao & J.R. He and to find out whether it is of hybrid origin. Whole-genome sequencing data were obtained from 57 samples. The SNP matrix was developed for phylogenetic reconstruction, ABBA-BABA statistics, PCA and ADMIXTURE analysis. The results indicate that P.gonggaensis is a spontaneous hybrid between P.lasiocarpa and P.cathayana. This study points out the importance of SNP data and comprehensive analyses for discovering the potential interspecific hybridisation and clarifies the usage of the name. In addition, the lectotype of P.gonggaensis was designated.

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.

Integrative Analysis of Oleosin Genes Provides Insights into Lineage-Specific Family Evolution in Brassicales

Oleosins (OLEs) are a class of small but abundant structural proteins that play essential roles in the formation and stabilization of lipid droplets (LDs) in seeds of oil crops. Despite the proposal of five oleosin clades (i.e., U, SL, SH, T, and M) in angiosperms, their evolution in eudicots has not been well-established. In this study, we employed Brassicales, an economically important order of flowering plants possessing the lineage-specific T clade, as an example to address this issue. Three to 10 members were identified from 10 species representing eight plant families, which include Caricaceae, Moringaceae, Akaniaceae, Capparaceae, and Cleomaceae. Evolutionary and reciprocal best hit-based homologous analyses assigned 98 oleosin genes into six clades (i.e., U, SL, SH, M, N, and T) and nine orthogroups (i.e., U1, U2, SL, SH1, SH2, SH3, M, N, and T). The newly identified N clade represents an ancient group that has already appeared in the basal angiosperm Amborella trichopoda, which are constitutively expressed in the tree fruit crop Carica papaya, including pulp and seeds of the fruit. Moreover, similar to Clade N, the previously defined M clade is actually not Lauraceae-specific but an ancient and widely distributed group that diverged before the radiation of angiosperm. Compared with A. trichopoda, lineage-specific expansion of the family in Brassicales was largely contributed by recent whole-genome duplications (WGDs) as well as the ancient γ event shared by all core eudicots. In contrast to the flower-preferential expression of Clade T, transcript profiling revealed an apparent seed/embryo/endosperm-predominant expression pattern of most oleosin genes in Arabidopsis thaliana and C. papaya. Moreover, the structure and expression divergence of paralogous pairs was frequently observed, and a good example is the lineage-specific gain of an intron. These findings provide insights into lineage-specific family evolution in Brassicales, which facilitates further functional studies in nonmodel plants such as C. papaya.

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.

Heterologous expression of Arabidopsis laccase2, laccase4 and peroxidase52 driven under developing xylem specific promoter DX15 improves saccharification in populus

BACKGROUND

Secondary cell wall holds considerable potential as it has gained immense momentum to replace the lignocellulosic feedstock into fuels. Lignin one of the components of secondary cell wall tightly holds the polysaccharides thereby enhancing the recalcitrance and complexity in the biomass. Laccases (LAC) and peroxidases (PRX) are the major phenyl-oxidases playing key functions during the polymerization of monolignols into lignin. Yet, the functions of laccase and peroxidases gene families remained largely unknown. Hence, the objective of this conducted study is to understand the role of specific LAC and PRX in Populus wood formation and to further investigate how the altered Lac and Prx expression affects biomass recalcitrance and plant growth. This study of heterologous expression of Arabidopsis Lac and Prx genes was conducted in poplar to avoid any otherwise occurring co-suppression mechanism during the homologous overexpression of highly expressed native genes. In the pursuit of optimizing lignocellulosic biomass for biofuel production, the present study focuses on harnessing the enzymatic potential of Arabidopsis thaliana Laccase2, Laccase4, and Peroxidase52 through heterologous expression.

RESULTS

We overexpressed selected Arabidopsis laccase2 (AtLac2), laccase4 (AtLac4), and peroxidase52 (AtPrx52) genes, based on their high transcript expression respective to the differentiating xylem tissues in the stem, in hybrid poplar (cv. 717) expressed under the developing xylem tissue-specific promoter, DX15 characterized the transgenic populus for the investigation of growth phenotypes and recalcitrance efficiency. Bioinformatics analyses conducted on AtLac2 and AtLac4 and AtPrx52, revealed the evolutionary relationship between the laccase gene and peroxidase gene homologs, respectively. Transgenic poplar plant lines overexpressing the AtLac2 gene (AtLac2-OE) showed an increase in plant height without a change in biomass yield as compared to the controls; whereas, AtLac4-OE and AtPrx52-OE transgenic lines did not show any such observable growth phenotypes compared to their respective controls. The changes in the levels of lignin content and S/G ratios in the transgenic poplar resulted in a significant increase in the saccharification efficiency as compared to the control plants.

CONCLUSIONS

Overall, saccharification efficiency was increased by 35-50%, 21-42%, and 8-39% in AtLac2-OE, AtLac4-OE, and AtPrx52-OE transgenic poplar lines, respectively, as compared to their controls. Moreover, the bioengineered plants maintained normal growth and development, underscoring the feasibility of this approach for biomass improvement without compromising overall plant fitness. This study also sheds light on the potential of exploiting regulatory elements of DX15 to drive targeted expression of lignin-modifying enzymes, thereby providing a promising avenue for tailoring biomass for improved biofuel production. These findings contribute to the growing body of knowledge in synthetic biology and plant biotechnology, offering a sustainable solution to address the challenges associated with lignocellulosic biomass recalcitrance.

Gender specific SNP markers in Coscinium fenestratum (Gaertn.) Colebr. for resource augmentation

BACKGROUND

Unregulated extraction of highly traded medicinal plant species results in drastic decline of the natural resources and alters viable sex ratio of populations. Conservation and long-term survival of such species, require gender specific restoration programs to ensure reproductive success. However, it is often difficult to differentiate sex of individuals before reaching reproductive maturity. C. fenestratum is one of the medicinally important and overexploited dioecious woody liana, with a reproductive maturity of 15 years. Currently, no information is available to identify sex of C. fenestratum in seedling stage while augmenting the resources. Thus, the current study envisages to utilize transcriptomics approach for gender differentiation which is imperative for undertaking viable resource augmentation programmes.

METHODS AND RESULTS

Gender specific SNPs with probable role in sexual reproduction/sex determination was located using comparative transcriptomics approach (sampling male and female individuals), alongside gene ontology and annotation. Nine sets of primers were synthesized from 7 transcripts (involved in sexual reproduction/other biological process) containing multiple SNP variants. Out of the nine primer pairs, only one SNP locus with no available information of its role in reproduction, showed consistent and accurate results (males-heterozygous and females-homozygous), in the analyzed 40 matured individuals of known sexes. Thus validated the efficiency of this SNP marker in differentiating male and female individuals.

CONCLUSIONS

The study could identify SNPs linked to the loci with apparent role in gender differentiation. This SNP marker can be used for early sexing of seedlings for in-situ conservation and resource augmentation of C. fenestratum in Kerala, India.