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

PtrVINV2 is dispensable for cellulose synthesis but essential for salt tolerance in Populus trichocarpa Torr. and Gray

Invertase (EC.3.2.1.26), a key enzyme in sucrose breakdown, is crucial for cellulose synthesis. However, the function of the vacuolar invertase (VINV) in woody plants remains unclear. In this study, transgenic lines of Populus trichocarpa Torr. and Gray were generated to investigate the role of PtrVINV2 in wood formation and under high salinity stress. Compared to wild-type (WT), VINV activity in the developing xylem of knockout lines was reduced, resulting in a decrease in lignin content and an increase in hemicellulose content, while cellulose content remained unaffected. These changes in structural carbohydrate content were accompanied by reductions in xylem width and fibre cell wall thickness. The overexpression lines of the developing xylem exhibited opposite trends. Transcriptome analyses of developing xylem indicated that the expression level of PtrVINV2 affects the expression of genes involved in hemicellulose and lignin biosynthesis pathways, such as AXS, UAMs, HCT, COMT, CAD and peroxidases, while CesA expression remained unaffected. WGCNA analysis revealed that Potri.001G219100, Potri.009G106600 and Potri.002G081000 serve as 'hub' transcription factor genes within the structural/non-structural carbohydrate modules of PtrVINV2 transgenic lines, potentially involved in plant salt tolerance. Additionally, under 200 mmol/L NaCl treatment, the knockout lines exhibited increased salt sensitivity compared to WT. This increased sensitivity was accompanied by elevated activities of SOD, CAT and MDA, as well as higher sucrose content and reduced contents of glucose and fructose. The findings indicate that although PtrVINV2 is not essential for cellulose synthesis, it enhances salt tolerance in poplar and presents a promising candidate gene for breeding salt-tolerant poplar.

Genome-wide identification, expression analysis, and stress response analysis of the RdbZIP gene family in Rhododendron delavayi

BACKGROUND

Basic leucine zipper (bZIP) gene family members represent one of the most diverse and largest groups of transcription factors in eukaryotes. Research has demonstrated that bZIP transcription factors play crucial roles not only in plant growth and development but also in response to various abiotic stresses. However, studies focusing on bZIP factors in Rhododendron delavayi (RdbZIPs) remain limited.

RESULT

In this study, a total of 59 RdbZIPs were identified using bioinformatics approaches, and these could be classified into 13 subfamilies based on the genomic data of R. delavayi. Members of the same RdbZIP subfamily exhibited similar gene structures and conserved motifs, and were unevenly distributed across the 13 chromosomes of R. delavayi. Collinearity analysis revealed a total of 20 duplication events, comprising 3 pairs of tandem duplications and 17 pairs of segmental duplications. Additionally, cis-acting element analysis indicated that RdbZIP family members may be involved in various biological processes, including transcription, development, hormone regulation, and responses to biotic and abiotic stresses. Transcriptomic analysis revealed that RdbZIP family genes were highly expressed in reproductive tissues. RT-qPCR expression analysis revealed that many selected RdbZIP genes were significantly upregulated under high salinity and drought conditions, suggesting their potential involvement in stress-responsive regulatory networks.

CONCLUSION

This study provides the first comprehensive characterization of the bZIP transcription factor family in Rhododendron delavayi, laying a foundational framework for functional studies of individual RdbZIP genes. The results highlight the pivotal role of RdbZIP genes in abiotic stress tolerance, which is crucial for understanding the adaptive mechanisms of R. delavayi. Future research should focus on the functional validation of key RdbZIP genes and elucidation of their regulatory pathways, which may contribute to the genetic improvement of Rhododendron species under adverse environmental conditions.

CLINICAL TRIAL

Not applicable.

PagWRKY11 regulates leaf morphology and salt sensitivity in Populus alba×P.glandulosa

WRKY transcription factors (TFs) are key regulators of plant tissue morphogenesis, defense responses, and metabolic regulation. However, the functions for most of WRKY TFs in 84K poplar (Populus alba × P. glandulosa) in regulating leaf morphology and responding to salt stress are remain unclear. In this study, overexpressing PagWRKY11 poplars were generated. Phenotypic analysis revealed that transgenic poplar leaves were narrower and smoother compared to the traditionally elliptical and relatively rough leaves of wild-type (WT) plants. Then, the apical buds of transgenic poplars were sharp and elongated, with young leaves pointing upwards and inwards, whereas WT buds were rounder with smoother young leaves pointing downwards. Functional analysis indicated that under salt stress, the activities of SOD and POD enzymes and the expression of their encoding genes were significantly lower in transgenic poplars compared to WT. Conversely, the accumulation of H2O2 and MDA was significantly higher. These results suggest that overexpression lines of PagWRKY11 increase salt sensitivity by downregulating the expression of antioxidant enzyme genes. Meanwhile, overexpression of PagWRKY11 increased the natural water loss rate of poplar leaves, and negatively regulated salt stress by affecting water retention. In addition, yeast one-hybrid assays showed that PagWRKY11 binds specifically to W-box elements. These results provide a theoretical basis for further exploration of the molecular mechanisms by which PagWRKY11 regulates leaf morphogenesis and stress responses, and offer new potential strategies for resistance breeding.

Evolution of NAC transcription factors from early land plants to domesticated crops

NAC [NO APICAL MERISTEM (NAM), ARABIDOPSIS TRANSCRIPTION ACTIVATOR FACTOR 1/2 (ATAF1/2), and CUP-SHAPED COTYLEDON (CUC2)] transcription factors are key regulators of plant growth, development, and stress responses but were also crucial players during land plant adaptation and crop domestication. Using representative members of green algae, bryophytes, lycophytes, gymnosperms, and angiosperms, we expanded the evolutionary history of NAC transcription factors to unveil the relationships among members of this gene family. We found a massive increase in the number of NAC transcription factors from green algae to lycophytes and an even larger increase in flowering plants. Many of the NAC clades arose later during evolution since we found eudicot- and monocot-specific clades. Cis-elements analysis in NAC promoters showed the presence of abiotic and biotic stress as well as hormonal response elements, which indicate the ancestral function of NAC transcription factor genes in response to environmental stimuli and in plant development. At the transcriptional level, the expression of NAC transcription factors was low or absent in male reproduction, particularly mature pollen, across the plant kingdom. We also identified NAC genes with conserved expression patterns in response to heat stress in Marchantia polymorpha and Oryza sativa. Our study provides further evidence that transcriptional mechanisms associated with stress responses and development emerged early during plant land adaptation and are still conserved in flowering plants and domesticated crops.

Generation and Functional Characteristics of CRISPR/Cas9-Edited PtrPHOTs Triple-Gene Mutants in Poplar

Phototropins (PHOTs), as blue light receptors, play a pivotal role in plant light signal perception and adaptive regulation, yet their functional characteristics in trees remain poorly understood. In this study, the PHOT gene family was identified in Populus trichocarpa, and it included three members, PtrPHOT1, PtrPHOT2.1, and PtrPHOT2.2, all of which were highly expressed in mature leaves. Using CRISPR/Cas9 gene editing technology, triple-gene mutations in the PtrPHOT1/2.1/2.2 (PtrPHOTs) were generated, providing initial insights into the functions of PHOTs in trees. Compared to the wild type (WT), triple-gene ptrphots mutants displayed curved and wrinkled leaves, reduced leaf area, and delayed phototropic responses, indicating the central role of PHOTs in blue light signal perception. The stomatal aperture recovery rate in mutants was only 40% of that observed in WT, accompanied by significant downregulation of the BLUS1 gene transcription levels, confirming the conservation of the PHOT-BLUS1-H⁺-ATPase signaling axis in stomatal regulation. Transcriptome of triple-gene ptrphots mutants revealed 1413 differentially expressed genes, of which were enriched in auxin response (upregulation of SAUR family genes), jasmonic acid (downregulation of JAZ genes), and light signaling pathways, suggesting that PHOTs could regulate plant adaptability by integrating light signals and hormone homeostasis. Overall, this study achieved the knockouts of three PtrPHOTs family genes, and characteristics of triple-gene ptrphots mutants elucidated the multifunctional roles of PHOTs in leaf development, phototropism, and stomatal movement in poplar. Our work provides a foundation for deciphering light signaling networks and molecular breeding in woody plants.

Transcriptomic analysis of peaches and nectarines reveals alternative mechanism for trichome formation

Trichomes in Prunus persica (L.) Batsch are crucial specialized structures that play a protective role against both biotic and abiotic stresses. The fruits with and without trichomes are respectively named as peach and nectarine. At the genetic level, the formation of trichome in peach is controlled by a single gene, PpMYB25, at the G locus. Peach (GG or Gg) is dominant to nectarine (gg), but such regulatory role was reported in a small-scale accession. In this study, we performed large-scale genotype and phenotype screening on 295 accessions. Almost all accessions supported the casual relationship between trichome formation and PpMYB25. However, a peach to nectarine mutant, named Maravilha Nectarine Mutant (MN), was discovered to possess a putative functional PpMYB25 gene sequence (Gg) but revealed nectarine phenotype. Comparative transcriptomic analyses revealed that PpMYB25 transcript was absent in MN. Correlation analyses also demonstrated that the PpMYB25-mediated regulatory network was abolished in MN. In summary, our results demonstrated an alternative mechanism beyond genetic regulation on trichome formation.

PtrCWINV3 encoding a cell wall invertase regulates carbon flow to wood in Populus trichocarpa

Cell wall invertase (CWINV) catalyzes the hydrolysis of sucrose into glucose and fructose in the apoplastic unloading pathway, with carbon sources provided for sink tissues. However, its role in wood formation remains undetermined. Therefore, transgenic lines overexpressing PtrCWINV3 or with knocked-out PtrCWINV3 expression were generated in Populus trichocarpa. Compared with wild type, the PtrCWINV3-knockout lines showed decreased CWINV activity (by 7.4 %-10.8 %), which resulted in a 1.5 %-1.8 % decrease in cellulose content, a 0.82 %-0.98 % decrease in hemicellulose content, and an increase in lignin content (by 2.9 %-4.7 %). These changes in structural carbohydrate contents were accompanied with anomalies in the late stages of secondary xylem development, characterized by reduced width of the secondary xylem, fewer cell layers in secondary xylem, and thinner fiber cell walls. The lines overexpressing PtrCWINV3 under the control of the DX15 promoter in the developing xylem showed the opposite phenotype. Transcriptome data from the developing xylem indicated that PtrCWINV3 regulated the expression of genes involved in the biosynthesis of cellulose (CesA, EG, and CB), hemicellulose/pectin (UGD, AXS, GATL, UAM, PAE, and GAUT), and starch (GBSS), which suggested its involvement in multiple polysaccharide metabolic pathways. Ultimately, this facilitated the synthesis of structural carbohydrate components such as cellulose and hemicellulose, which promoted the later stages of secondary xylem development. These findings not only demonstrate the significant role of CWINV activity in wood formation, but also highlight an excellent candidate gene for breeding new poplar varieties with high cellulose and low lignin contents.

Construction of multi-targeted CRISPR libraries in tomato to overcome functional redundancy at genome-scale level

Genetic variance is vital for breeding programs and mutant screening, yet traditional mutagenesis methods wrestle with genetic redundancy and a lack of specificity in gene targeting. CRISPR-Cas9 offers precise, site-specific gene editing, but its application in crop improvement has been limited by scalability challenges. In this study, we develop genome-wide multi-targeted CRISPR libraries in tomato, enhancing the scalability of CRISPR gene editing in crops and addressing the challenges of redundancy while maintaining its precision. We design 15,804 unique single guide RNAs (sgRNAs), each targeting multiple genes within the same gene families. These sgRNAs are classified into 10 sub-libraries based on gene function. We generate approximately 1300 independent CRISPR lines and successfully identify mutants with distinct phenotypes related to fruit development, fruit flavor, nutrient uptake, and pathogen response. Additionally, we develop CRISPR-GuideMap, a double-barcode tagging system to enable large-scale sgRNA tracking in generated plants. Our results demonstrate that multi-targeted CRISPR libraries are scalable and effective for large-scale gene editing and offer an approach to overcome gene functional redundancy in basic plant research and crop breeding.

Four near-complete genome assemblies reveal the landscape and evolution of centromeres in Salicaceae

BACKGROUND

Centromeres play a crucial role in maintaining genomic stability during cell division. They are typically composed of large arrays of tandem satellite repeats, which hinder high-quality assembly and complicate our efforts to understand their evolution across species. Here, we use long-read sequencing to generate near-complete genome assemblies for two Populus and two Salix species belonging to the Salicaceae family and characterize the genetic and epigenetic landscapes of their centromeres.

RESULTS

The results show that only limited satellite repeats are present as centromeric components in these species, while most of them are located outside the centromere but exhibit a homogenized structure similar to that of the Arabidopsis centromeres. Instead, the Salicaceae centromeres are mainly composed of abundant transposable elements, including CRM and ATHILA, while LINE elements are exclusively discovered in the poplar centromeres. Comparative analysis reveals that these centromeric repeats are extensively expanded and interspersed with satellite arrays in a species-specific and chromosome-specific manner, driving rapid turnover of centromeres both in sequence compositions and genomic locations in the Salicaceae.

CONCLUSIONS

Our results highlight the dynamic evolution of diverse centromeric landscapes among closely related species mediated by satellite homogenization and widespread invasions of transposable elements and shed further light on the role of centromere in genome evolution and species diversification.

Effect of temperature on circadian clock functioning of trees in the context of global warming

Plant survival in a warmer world requires the timely adjustment of biological processes to cyclical changes in the new environment. Circadian oscillators have been proposed to contribute to thermal adaptation and plasticity. However, the influence of temperature on circadian clock performance and its impact on plant behaviour in natural ecosystems are not well-understood. We combined bioinformatics, molecular biology and ecophysiology to investigate the effects of increasing temperatures on the functioning of the circadian clock in two closely related tree species from Patagonian forests that constitute examples of adaptation to different thermal environments based on their altitudinal profiles. Nothofagus pumilio, the species from colder environments, showed a major rearrangement of its transcriptome and reduced ability to maintain rhythmicity at high temperatures compared with Nothofagus obliqua, which inhabits warmer zones. In altitude-swap experiments, N. pumilio, but not N. obliqua, showed limited oscillator function in warmer zones of the forest, and reduced survival and growth. Our findings show that interspecific differences in the influence of temperature on circadian clock performance are associated with preferred thermal niches, and to thermal plasticity of seedlings in natural environments, highlighting the potential role of a resonating oscillator in ecological adaptation to a warming environment.

TCP23-WRKY15 module negatively regulates lignin deposition and xylem development of wood formation in Populus

Secondary wall, a critical component of wood, is influenced by multiple factors during its formation. The TCP family encodes plant-specific transcription factors (TFs) that play key roles in multiple aspects of plant development. In this study, we identified all TCP TFs in five poplar species and analyzed their evolutionary relationships, gene structures, tissue-specific expression patterns, and potential interactions with microRNAs. Additionally, we screened for TCP proteins associated with secondary wall development that are independent of miRNA regulation. Three candidate TFs were identified, with TCP23 showing high conservation across poplar species and the highest expression levels in the xylem of Populus trichocarpa and Populus wilsonii. The overexpression of TCP23 in poplar inhibited the expression of MYB TFs and structural genes involved in xylem biosynthesis, thereby reducing the lignin content within the stems. By contrast, CRISPR/Cas9-mediated knockout of TCP23 resulted in the opposite effect. Furthermore, we successfully identified WRKY15 as an interaction partner of TCP23 via a yeast two-hybrid library and demonstrated that TCP23 negatively regulates lignin synthesis and xylem development by enhancing the inhibitory function of WRKY15. Our study provides new insights into the transcriptional regulatory mechanisms underlying secondary wall formation.

An integrated transcriptome, metabolome, and microbiome dataset of Populus under nutrient-poor conditions

The rhizosphere microbiota recruited by plants contributes significantly to maintaining host productivity and resisting stress. However, the genetic mechanisms by which plants regulate this recruitment process remain largely unclear. Here, we generated a comprehensive dataset, including 27 root transcriptomes, 27 root metabolomes, and 54 bulk or rhizosphere soil 16S rRNA amplicons across nine poplar species from four sections grown in nutrient-poor natural soil, along with eleven growth phenotype data. We provided a thorough description of this dataset, followed by a comprehensive co-expression network analysis example that broke down the wall of the four-way relationship between plant gene-metabolite-microbe-phenotype, thus identifying the links between plant gene expression, metabolite accumulation, growth behavior, and rhizosphere microbiome variation under nutrient-poor conditions. Overall, this dataset enhances our understanding of plant and microbe interactions, offering valuable strategies and novel insights for resolving how plants regulate rhizosphere microbial compositions and functions, thereby improving host fitness, which will benefit future research.

Bioengineered poplar fibres via PagGLR2.8 editing: A synergistic design for high-performance biocomposites

The urgent need to replace petroleum-derived materials with sustainable alternatives drives innovation at the nexus of plant biotechnology and materials science. Here, we engineered Populus alba × P. glandulosa '84 K' through CRISPR-Cas9-mediated knockout of PagGLR2.8, a glutamate receptor gene regulating vascular development, to investigate its role in fibre biosynthesis and composite performance. Knockout of PagGLR2.8 improved the quality of poplar fibre by altering the structure and development mode of poplar vascular tissue. Our study established the relationship between fibre quantity and structure and the performance of polylactic acid (PLA) composites. The mechanical and fire-resistance properties of these transgenic plant fibres/PLA composites significantly outperformed those of pure PLA, demonstrating the potential of phloem fibres to reinforce toughened composites. Notably, we also evaluated flammability and dripping behaviours, with findings indicating that our optimised fibre/PLA composites exhibit superior strengths, modulus, fire resistance, and anti-dripping, surpassing those of PLA. This research unveils a groundbreaking approach to regulating composite properties through genetic manipulation and highlights the promising potential of plant-derived materials in enriching forest resources and advancing the sustainable utilisation of poplar fibres and polymers.

idopNetwork Analysis of Salt-Responsive Transcriptomes Reveals Hub Regulatory Modules and Genes in Populus euphratica

Euphrates poplar (Populus euphratica) is known as a system model to study the genomic mechanisms underlying the salt resistance of woody species. To characterize how dynamic gene regulatory networks (GRNs) drive the defense response of this species to salt stress, we performed mRNA sequencing of P. euphratica roots under short-term (ST) and long-term (LT) salt stress treatments across multiple time points. Comparisons of these transcriptomes revealed the diverged gene expression patterns between the ST and LT treated samples. Based on the informative, dynamic, omnidirectional, and personalized networks model (idopNetwork), inter- and intra-module networks were constructed across different time points for both the ST and LT groups. Through the analysis of the inter-module network, we identified module 4 as the hub, containing the largest number of genes. Further analysis of the gene network within module 4 revealed that gene XM_011048240.1 had the most prominent interactions with other genes. Under short-term salt stress, gene interactions within the network were predominantly promoted, whereas under long-term stress, these interactions shifted towards inhibition. As for the gene ontology (GO) annotation of differentially expressed genes, the results suggest that P. euphratica may employ distinct response mechanisms during the early and late stages of salt stress. Taking together, these results offer valuable insights into the regulatory mechanism involved in P. euphratica's stress response, advancing our understanding of complex biological processes.

Transposon-Associated Small RNAs Involved in Plant Defense in Poplar

Utilizing high-throughput Illumina sequencing, we examined how small RNA (sRNA) profiles vary in Chinese white poplar (Populus tomentosa) across two pivotal infection stages by the rust fungus Melampsora larici-populina: the biotrophic growth phase (T02; 48 h post infection) and the urediniospore development and dispersal phase (T03; 168 h), both essential for plant colonization and prolonged biotrophic engagement. Far exceeding random expectations, siRNA clusters predominantly arose from transposon regions, with pseudogenes also contributing significantly, and infection-stage-specific variations were notably tied to these transposon-derived siRNAs. As the infection advanced, clusters of 24 nt siRNAs in transposon and intergenic regions exhibited pronounced abundance shifts. An analysis of targets indicated that Populus sRNAs potentially regulate 95% of Melampsora larici-populina genes, with pathogen effector genes showing heightened targeting by sRNAs during the biotrophic and urediniospore phases compared to controls, pointing to selective sRNA-target interactions. In contrast to conserved miRNAs across plant species, Populus-specific miRNAs displayed a markedly greater tendency to target NB-LRR genes. These observations collectively highlight the innovative roles of sRNAs in plant defense, their evolutionary roots, and their dynamic interplay with pathogen coevolution.

Genome-Wide Identification of the Dof Gene Family and Functional Analysis of PeSCAP1 in Regulating Guard Cell Maturation in Populus euphratica

DNA-binding with one finger (Dof) transcription factors plays critical roles in regulating plant growth and development, as well as modulating responses to biotic and abiotic stresses. While the biological characteristics of the Dof family have been explored across various species, their functions in Populus euphratica remain largely uncharacterized. In this study, we identified 43 PeDof family genes through a genome-wide approach, revealing a total of 10 conserved motifs across all family members. Predictions of cis-acting elements indicated that Dof genes are involved in light signaling, hormone signaling, and stress responses. Phylogenetic analysis classified the 43 Dof genes of P. euphratica into six distinct groups, with genes within the same group exhibiting relatively conserved structures. Expression pattern analyses demonstrated significant regulation of PeDof genes by drought stress, with their expression also being influenced by environmental conditions during seed germination. Furthermore, we identified the Dof gene PeSCAP1, which plays a conserved role in regulating guard cell maturation, underscoring the importance of stomatal morphology and function in leaf water retention. This study enhances our understanding of the role of Dofs in abiotic stress responses and provides valuable insights into their function in Populus euphratica.

Bioinformatics Analysis Reveals the Evolutionary Characteristics of the Phoebe bournei ARF Gene Family and Its Expression Patterns in Stress Adaptation

Auxin response factors (ARFs) are pivotal transcription factors that regulate plant growth, development, and stress responses. Yet, the genomic characteristics and functions of ARFs in Phoebe bournei remain undefined. In this study, 25 PbARF genes were identified for the first time across the entire genome of P. bournei. Phylogenetic analysis categorized these genes into five subfamilies, with members of each subfamily displaying similar conserved motifs and gene structures. Notably, Classes III and V contained the largest number of members. Collinearity analysis suggested that segmental duplication events were the primary drivers of PbARF gene family expansion. Structural analysis revealed that all PbARF genes possess a conserved B3 binding domain and an auxin response element, while additional motifs varied among different classes. Promoter cis-acting element analysis revealed that PbARF genes are extensively involved in hormonal responses-particularly to abscisic acid and jasmonic acid and abiotic stresses-as well as abiotic stresses, including heat, drought, light, and dark. Tissue-specific expression analysis showed that PbARF25, PbARF23, PbARF19, PbARF22, and PbARF20 genes (class III), and PbARF18 and PbARF11 genes (class V) consistently exhibited high expression levels in the five tissues. In addition, five representative PbARF genes were analyzed using qRT-PCR. The results demonstrated significant differences in the expression of PbARF genes under various abiotic stress conditions (drought, salt stress, light, and dark), indicating their important roles in stress response. This study laid a foundation for elucidating the molecular evolution mechanism of ARF genes in P. bournei and for determining the candidate genes for stress-resistance breeding.

Integrated genomic and transcriptomic analyses reveal the genetic and molecular mechanisms underlying hawthorn peel color and seed hardness diversity

Hawthorn (Crataegus pinnatifida) fruit peel color and seed hardness are key traits that significantly impact economic value. We present here the high-quality chromosome-scale genomes of two cultivars, including the hard-seed, yellow-peel C. pinnatifida "Jinruyi" (JRY) and the soft-seed, red-peel C. pinnatifida "Ruanzi" (RZ). The assembled genomes comprising 17 chromosomes are 809.1 Mb and 760.5 Mb in size, achieving scaffold N50 values of 48.5 Mb and 46.8 Mb for JRY and RZ, respectively. Comparative genomic analysis identifies 3.6-3.8 million single nucleotide polymorphisms, 8.5-9.3 million insertions/deletions, and approximately 30 Mb of presence/absence variations across different hawthorn genomes. Through integrating differentially expressed genes and accumulated metabolites, we filter candidate genes CpMYB114 and CpMYB44 associated with differences in hawthorn fruit peel color and seed hardness, respectively. Functional validation confirms that the CpMYB114-CpANS regulates anthocyanin biosynthesis in hawthorn peels, contributing to the observed variation in peel color. CpMYB44-CpCOMT is significantly upregulated in JRY and is verified to promote lignin biosynthesis, resulting in the distinction in seed hardness. Overall, this study reveals the new insights into understanding of distinct peel pigmentation and seed hardness in hawthorn and provides an abundant resource for molecular breeding.

SOI: robust identification of orthologous synteny with the Orthology Index and broad applications in evolutionary genomics

With the explosive growth of whole-genome datasets, accurate detection of orthologous synteny has become crucial for reconstructing evolutionary history. However, current methods for identifying orthologous synteny face great limitations, particularly in scaling with varied polyploidy histories and accurately removing out-paralogous synteny. In this study, we developed a scalable and robust approach, based on the Orthology Index (OI), to effectively identify orthologous synteny. Our evaluation across a large-scale empirical dataset with diverse polyploidization events demonstrated the high reliability and robustness of the OI method. Simulation-based benchmarks further validated the accuracy of our method, showing its superior performance against existing methods across a wide range of scenarios. Additionally, we explored its broad applications in reconstructing the evolutionary histories of plant genomes, including the inference of polyploidy, identification of reticulation, and phylogenomics. In conclusion, OI offers a robust, interpretable, and scalable approach for identifying orthologous synteny, facilitating more accurate and efficient analyses in plant evolutionary genomics.

Identification of phloem-specific proteinaSEOus structure heterogeneity in sieve element of Populus trichocarpa

Phloem, an exceptional plant vascular tissue, facilitates the transport of photoassimilates, RNAs, and other signaling substances from the leaves to the roots throughout the plant. Among the specialized phloem cells are the conductive sieve elements (SEs), which are unique in that they remain alive despite lacking several cell organelles, including the nucleus, plastids, and most mitochondria. These SEs contain a specific proteinaceous structure composed of phloem-specific proteins (P-proteins), whose function is not yet fully understood. Various P-proteins have been characterized in broad range of model species, including Arabidopsis thaliana, and reported in Fabaceae and Cucurbitaceae plants. To date, only one P-protein has been identified in the model tree species Populus trichocarpa. Given the presence of multiple P-protein encoding genes across numerous plant species, we hypothesized the existence of multiple such genes in the Populus genome. Our genomic analysis uncovered 12 genes being potential orthologues to one of A. thaliana P-protein - SEOR (sieve element occlusion-related) genes, which may contribute to the proteinaceous structures observed in differentiating sieve elements. Our transcriptomic and proteomic analyses confirmed the expression of at least seven of these genes, indicating that the protein structure visible in mature sieve elements in P. trichocarpa may be heterogeneous.

Chromosome-scale and haplotype-resolved genome assembly of Populus trichocarpa

Populus trichocarpa, a pivotal model organism for woody transgenic research, not only garners substantial scientific interest but plays an integral role in forestry economics. Previous genomic assemblies of P. trichocarpa predominantly treated its heterozygous genome as homozygous, thereby neglecting crucial haplotypic diversity. Leveraging the high-fidelity (HiFi) sequencing capabilities of PacBio sequencing and the chromosome conformation capture insights provided by Illumina's Hi-C technique, this study is the first to achieve a near telomere-to-telomere assembly of both paternal and maternal haplotypes in P. trichocarpa. Comparative genomic analysis between these haplotypes has uncovered several allelic variants and pathways critical for trait determination through allele-specific expression. Furthermore, utilizing RNA-seq data from multiple tissues, this investigation has detailed the tissue-specific expression patterns of the leucine-rich repeat gene family, which are essential in mediating plant signal transduction and developmental regulation. Our results not only illuminate the functional genomics landscape of P. trichocarpa but also provide invaluable theoretical underpinnings for the genetic improvement of woody plants and a robust framework for exploring genetic variability and allelic expression disparities in arboreal species.

Chromosomal-Level Genome Suggests Adaptive Constraints Leading to the Historical Population Decline in an Extremely Endangered Plant

Increased human activity and climate change have significantly impacted wild habitats and increased the number of endangered species. Exploring evolutionary history and predicting adaptive potential using genomic data will facilitate species conservation and biodiversity recovery. Here, we examined the genome evolution of a critically endangered tree Pellacalyx yunnanensis, a plant species with extremely small populations (PSESP) that is narrowly distributed in Xishuangbanna, China. The species has neared extinction due to economic exploitation in recent decades. We assembled a chromosome-level genome of 334 Mb, with the N50 length of 20.5 Mb. Using the genome, we discovered that P. yunnanensis has undergone several population size reductions, leading to excess deleterious mutations. The species may possess low adaptive potential due to reduced genetic diversity and the loss of stress-responsive genes. We estimate that P. yunnanensis is the basal species of its genus and diverged from its relatives during global cooling, suggesting it was stranded in unsuitable environments during periods of dramatic climate change. In particular, the loss of seed dormancy leads to germination under unfavourable conditions and reproduction challenges. This dormancy loss may have occurred through genetic changes that suppress ABA signalling and the loss of genes involved in seed maturation. The high-quality genome has also enabled us to reveal phenotypic trait evolution in Rhizophoraceae and identify divergent adaptation to intertidal and inland habitats. In summary, our study elucidates mechanisms underlying the decline and evaluates the adaptive potential of P. yunnanensis to future climate change, informing future conservation efforts.

Optimising Exome Captures in Species With Large Genomes Using Species-Specific Repetitive DNA Blocker

Large and highly repetitive genomes are common. However, research interests usually lie within the non-repetitive parts of the genome, as they are more likely functional, and can be used to answer questions related to adaptation, selection and evolutionary history. Exome capture is a cost-effective method for providing sequencing data from protein-coding parts of the genes. C0t-1 DNA blockers consist of repetitive DNA and are used in exome captures to prevent the hybridisation of repetitive DNA sequences to capture baits or bait-bound genomic DNA. Universal blockers target repetitive regions shared by many species, while species-specific c0t-1 DNA is prepared from the DNA of the studied species, thus perfectly matching the repetitive DNA contents of the species. So far, the use of species-specific c0t-1 DNA has been limited to a few model species. Here, we evaluated the performance of blocker treatments in exome captures of Pinus sylvestris, a widely distributed conifer species with a large (> 20 Gbp) and highly repetitive genome. We compared treatment with a commercial universal blocker to treatments with species-specific c0t-1 (30,000 and 60,000 ng). Species-specific c0t-1 captured more unique exons than the initial set of targets leading to increased SNP discovery and reduced sequencing of tandem repeats compared to the universal blocker. Based on our results, we recommend optimising exome captures using at least 60,000 ng of species-specific c0t-1 DNA. It is relatively easy and fast to prepare and can also be used with existing bait set designs.

Genome-Wide Identification and Expression of the ERF Gene Family in Populus trichocarpa and Their Responses to Nitrogen and Abiotic Stresses

The ethylene response factor (ERF) family is a prominent plant-specific transcription factor family, which plays a crucial role in modulating plant growth and stress tolerance. In this study, a total of 210 ERFs were identified in Populus trichocarpa, comprising 29 AP2 (APETALA2) subfamily members, 176 ERF subfamily members, and 5 RAV (related to ABI3/VP1) subfamily members. The duplication events of the PtERF family members exclusively occurred within the subfamilies. A total of 168 duplication pairs were found among 161 PtERF genes, and all of them were fragment duplications. Gene structure analysis revealed that most ERF subfamily members only had one exon without introns, the AP2 subfamily members had six or more introns and exons, and RAV subfamily members lacked introns except for PtERF102. Considerable cis-acting elements associated with plant growth and development, stress response, hormone response, and light response were detected in the promoters of PtERF genes. The expression levels of PtERFs were highest in roots across tissues and in winter among seasons. Furthermore, the nitrate and urea stimulated the expression of PtERF genes. The co-expression network analysis based on PtERFs indicated their potential roles in hormone signaling, acyltransferase activity, and response to chemicals. This study provides novel insights into investigating the role of PtERFs in environmental stress in poplar species.

Identification of a root-specific expression promoter in poplar and its application in genetic engineering for cadmium phytoremediation

KEY MESSAGE

A promoter, PRSEP7, was identified and confirmed to be specifically expressed in poplar roots. Poplar PRSEP7::CadWp transgenic lines showed high phytoremediation of Cd(II)-contaminated WPM and soil. Cadmium ions (Cd(II)) are heavy metals that are difficult for organisms to decompose in our natural environment. The generation of plants by genetic engineering with a high ability to phytoremediate Cd(II) from the soil is an ideal biological remediation strategy. Here, we identified and confirmed a promoter, PRSEP7, that is specifically expressed in poplar (Populus L.) roots. The promoter of PRSEP7 was then used to construct the poplar root expression vector 2301S-root. The CadW gene encoding a carbonic anhydrase (CA) was reported to play important roles in the phytoremediation of Cd(II) in microorganisms in a previous study. The sequence of CadW was optimized for plants, and the resulting gene CadWp also showed high activity for sequestration of Cd(II). CadWp was then introduced to 2301S-root to generate the PRSEP7::CadWp construct. This construct was used to transform poplar via Agrobacterium-mediated transformation. A number of stable transgenic poplar lines were generated, and two lines were randomly selected to test their ability to phytoremediate Cd(II). With several parameter measurements, the two transgenic lines showed high phytoremediation of Cd(II) under multiple growth conditions. Overall, we generated elite plant materials for the phytoremediation of Cd(II) in this study.

Genome-wide identification of MYB gene family and exploration of selenium metabolism-related candidates in paper mulberry (Broussonetia papyrifera)

KEY MESSAGE

Genome-wide identified 144 MYB family members in B. papyrifera. Integrated correlation analysis and target gene-binding motif prediction indicate that BpMYB135 is vital in regulating selenium metabolism. Selenium is an essential micronutrient for maintaining the health of humans and animals. Broussonetia papyrifera, a forage tree with high nutritional value, exhibits a remarkable ability to accumulate selenium. Although previous studies have preliminarily unfolded the molecular mechanisms underlying selenium accumulation, the roles of transcription factors in regulating selenium uptake and transformation remain poorly understood. This study used various strategies including bioinformatic, physiological, and molecular experiments to explore candidates regarding Se metabolism. Briefly, 144 MYB transcription factor family members were identified and classified into four types (R1, R2R3, R1R2R3, and R4), with phylogenetic analysis further dividing them into 58 subfamilies. The promoters of those BpMYBs contain numerous cis-acting elements associated with plant growth, development, and stress response. qRT-PCR assay confirmed 8 of 15 BpMYBs exhibit a remarkable correlation with selenium content at the threshold absolute value of 0.5. Additionally, foliar application of exogenous abscisic acid (ABA), methyl jasmonate (MeJA), and salicylic acid (SA) reveals different response patterns of BpMYBs. The subcellular localization assay simultaneously verifies that the candidate BpMYB135 functions within the nucleus. Overall, this funding highlights the potential regulatory mechanisms of selenium metabolism in B. papyrifera, providing a foundation for improving its forage value through genetic modification.

A special short-wing petal faba genome and genetic dissection of floral and yield-related traits accelerate breeding and improvement of faba bean

BACKGROUND

A comprehensive study of the genome and genetics of superior germplasms is fundamental for crop improvement. As a widely adapted protein crop with high yield potential, the improvement in breeding and development of the seeds industry of faba bean have been greatly hindered by its giant genome size and high outcrossing rate.

RESULTS

To fully explore the genomic diversity and genetic basis of important agronomic traits, we first generate a de novo genome assembly and perform annotation of a special short-wing petal faba bean germplasm (VF8137) exhibiting a low outcrossing rate. Comparative genome and pan-genome analyses reveal the genome evolution characteristics and unique pan-genes among the three different faba bean genomes. In addition, the genome diversity of 558 accessions of faba bean germplasm reveals three distinct genetic groups and remarkable genetic differences between the southern and northern germplasms. Genome-wide association analysis identifies several candidate genes associated with adaptation- and yield-related traits. We also identify one candidate gene related to short-wing petals by combining quantitative trait locus mapping and bulked segregant analysis. We further elucidate its function through multiple lines of evidence from functional annotation, sequence variation, expression differences, and protein structure variation.

CONCLUSIONS

Our study provides new insights into the genome evolution of Leguminosae and the genomic diversity of faba bean. It offers valuable genomic and genetic resources for breeding and improvement of faba bean.

Comparative genomics provides insights into the biogeographic and biochemical diversity of meliaceous species

Meliaceous plants such as Azadirachta indica (neem) and Melia azedarach (chinaberry) contain large amounts of limonoids with unique anti-insect activities. However, genes responsible for downstream modifications of limonoids are not well known. Here, we improve the genome assemblies of neem and chinaberry to the telomere-to-telomere (T2T) level. Allopatric speciation of the two plants is confirmed by the lineage-specific inversion of chromosome 12 in the neem lineage. We further identify two BAHD-acetyltransferases (ATs) in chinaberry (MaAT8824 and MaAT1704) that catalyse acetylation at both the C-12 and C-3 hydroxyl groups of limonoids, whereas the syntenic neem copy (AiAT0635) does not possess this activity. A critical N-terminal region (SAGAVP) is crucial for the acetylation of AiAT0635, and swapping it into the MaAT8824 version (CHRSSG) can endow it with acetylation activity. Our improved genome assemblies provide insights into allopatric speciation of neem, as well as limonoid biosynthesis and chemical diversity in meliaceous plants.

Genomic, transcriptomic, and metabolomic analyses reveal convergent evolution of oxime biosynthesis in Darwin's orchid

Angraecum sesquipedale, also known as Darwin's orchid, possesses an exceptionally long nectar spur. Charles Darwin predicted the orchid to be pollinated by a hawkmoth with a correspondingly long proboscis, later identified as Xanthopan praedicta. In this plant-pollinator interaction, the A. sesquipedale flower emits a complex blend of scent compounds dominated by diurnally regulated oximes (R1R2C = N-OH) to attract crepuscular and nocturnal pollinators. The molecular mechanism of oxime biosynthesis remains unclear in orchids. Here, we present the chromosome-level genome of A. sesquipedale. The haploid genome size is 2.10 Gb and represents 19 pseudochromosomes. Cytochrome P450 encoding genes of the CYP79 family known to be involved in oxime biosynthesis in seed plants are not present in the A. sesquipedale genome nor the genomes of other members of the orchid family. Metabolomic analysis of the A. sesquipedale flower revealed a substantial release of oximes at dusk during the blooming stage. By integrating metabolomic and transcriptomic correlation approaches, flavin-containing monooxygenases (FMOs) encoded by six tandem-repeat genes in the A. sesquipedale genome are identified as catalyzing the formation of oximes present. Further in vitro and in vivo assays confirm the function of FMOs in the oxime biosynthesis. We designate these FMOs as orchid oxime synthases 1-6. The evolutionary aspects related to the CYP79 gene losses and neofunctionalization of FMO-catalyzed biosynthesis of oximes in Darwin's orchid provide new insights into the convergent evolution of biosynthetic pathways.

Exon disruptive variants in Populus trichocarpa associated with wood properties exhibit distinct gene expression patterns

Forest trees may harbor naturally occurring exon disruptive variants (DVs) in their gene sequences, which potentially impact important ecological and economic phenotypic traits. However, the abundance and molecular regulation of these variants remain largely unexplored. Here, 24,420 DVs were identified by screening 1014 Populus trichocarpa full genomes. The identified DVs were predominantly heterozygous with allelic frequencies below 5% (only 26% of DVs had frequencies greater than 5%). Using common garden-grown trees, DVs were assessed for gene expression variation in the developing xylem, revealing that their gene expression can be significantly altered, particularly for homozygous DVs (in the range of 27%-38% of cases depending on the studied common garden). DVs were further investigated for their correlations with 13 wood quality traits, revealing that, among the 148 discovered DV associations, 15 correlated with more than one wood property and six genes had more than one DV in their coding sequences associated with wood traits. Approximately one-third of DVs correlated with wood property variation also showed significant gene expression variation, confirming their non-spurious impact. These findings offer potential avenues for targeted introduction of homozygous mutations using tree biotechnology, and while the exact mechanisms by which DVs may directly influence wood formation remain to be unraveled, this study lays the groundwork for further investigation.

Chromosome-level de novo genome unveils the evolution of Gleditsia sinensis and thorns development

Gleditsia sinensis Lam. (G. sinensis) as an important species within the Leguminosae family, has been utilized in Chinese medicine for centuries, and its thorns serve as a chief medicinal ingredient. The absence of a comprehensive genome database has hindered its in-depth research. In this investigation, a chromosome-level de novo genome assembly of G. sinensis 'Yulin No.1' was achieved, which harbors a 786.13 Mb sized genome with 36,408 protein-coding genes and experiences two WGD events. The comparative and evolutionary analysis unveiled the close phylogenetic relationship between G. sinensis and eight other Leguminosae species. The WGCNA and gene family analysis further indicated that GsinMYB was involved in the development of thorns. This investigation offered a high-level genome of G. sinensis, facilitating comparisons in Leguminosae species evolution and functional elucidation. It also provided key insights for further research on the molecular regulation mechanisms of thorn development in plants and the molecular breeding of G. sinensis.

Functional Genomics of Legumes in Bulgaria-Advances and Future Perspectives

Members of the Leguminosae family are important crops that provide food, animal feed and vegetable oils. Legumes make a substantial contribution to sustainable agriculture and the nitrogen cycle through their unique ability to fix atmospheric nitrogen in agricultural ecosystems. Over the past three decades, Medicago truncatula and Lotus japonicus have emerged as model plants for genomic and physiological research in legumes. The advancement of innovative molecular and genetic tools, particularly insertional mutagenesis using the retrotransposon Tnt1, has facilitated the development of extensive mutant collections and enabled precise gene tagging in plants for the identification of key symbiotic and developmental genes. Building on these resources, twelve years ago, our research team initiated the establishment of a platform for functional genomic studies of legumes in Bulgaria. In the framework of this initiative, we conducted systematic sequencing of selected mutant lines and identified genes involved in plant growth and development for detailed functional characterization. This review summarizes our findings on the functions of selected genes involved in the growth and development of the model species, discusses the molecular mechanisms underlying important developmental processes and examines the potential for the translation of this fundamental knowledge to improve commercially important legume crops in Bulgaria and globally.

Genome-wide identification, characterization and expression analysis of tubulin gene family in Populus deltoides

BACKGROUND

Tubulin proteins, the main components of microtubules in all eukaryotes, are involved in numerous aspects of plant morphogenesis and adaptation to the environment. In woody plants, microtubules are closely associated with the orientation of cellulose microfibril deposition in the secondary xylem cells and thereby exert an influence on the strength and flexibility of wood. Three major types of tubulin proteins-α-, β- and γ-tubulin-are ubiquitously present in all flowering plants, with α- and β- tubulin serving as basic subunits of microtubules and γ-tubulin directing microtubule nucleation. Compared with herbaceous plants, information on tubulin gene family has been limited in forest trees. This study aimed to characterize the tubulin gene family in the model forest tree Populus deltoides.

RESULT

Based on the whole genome sequence of P. deltoides, 25 PdTubulin genes were identified, including 6 PdTUAs, 17 PdTUBs, and 2 PdTUBGs were identified, with an uneven distribution across 14 chromosomes. Unlike Arabidopsis, which has only three pairs of tubulin paralogs, nearly all PdTubulin were paralogous duplicates, primarily generated by p-whole genome duplication (WGD), γ-WGD, or segmental duplication, indicating multiple rounds of gene family expansion. After the duplication events, the number of TUA genes in Populus was more strictly constrained compared to TUB genes. All paralogous and orthologous tubulin pairs have been under strong purifying selection. Expression analysis revealed that each PdTubulin gene was preferentially expressed in one of three organs: root, leaf, or stem. 5 PdTUB paralogs exhibited similar expression patterns, suggesting potential redundancy. Additionally, expression analysis in male and female floral buds across developmental stages indicated that different members might be involved in sex-specific processes.

Divergent evolutionary paces among eudicot plants revealed by simultaneously duplicated genes produced billions of years ago

Polyploidization often occurs more than once along an evolutionary lineage to form extant plants. Major core eudicot plants share a whole-genome triplication (ceWGT), through which thousands of simultaneously duplicated genes are retained in extant genomes, providing a valuable starting line to check the difference in their evolutionary paces. Here, by characterizing the synonymous nucleotide substitutions (Ks) between these duplicates from 28 representative plants from 21 families, we checked the various evolutionary rates among plants among plants subjected to different rounds of extra polyploidization events. We found up to 68.04% difference in evolutionary rates among the selected plants. A statistical correlation analysis (correlation coefficient =0.57, at significant level = 0.01) indicated that plants affected by extra polyploidies have evolved faster than plants without such extra polyploidies showing that (additional) polyploidization has resulted in elevated genetic diversity. Comparing the plants affected by additional polyploidization and plants without it, the duplicated genes produced by the ceWGT and retained in extant genomes have gathered 4.75% more nucleotide substitutions in the former plants. By identifying the fast- and slowly evolving genes, we showed that genes evolving at divergent rates were often related to different evolutionary paths. By performing correction to evolutionary rates using a genome-scale approach, we revised the estimated timing of key evolutionary events. The present effort exploited the simultaneously duplicated genes produced by the shared polyploidization and help deepen the understanding of the role of polyploidization, especially its long-term effect in plant evolution and biological innovation.

Genome-wide identification and characterization of the bZIP family in the Mangrove Plant Kandelia obovata and its role in response to stress

BACKGROUND

The basic leucine zipper (bZIP) transcription factors play crucial roles in plant growth, development, and responses to environmental changes. The mangrove plant Kandelia obovata, native to subtropical and tropical coastal intertidal zones, has evolved various adaptive mechanisms to cope with unstable muddy substrates, tidal fluctuations, saltwater intrusion, and intense ultraviolet radiation. This study aims to provide a comprehensive characterization of the bZIP gene family in K. obovata and investigate its functional roles in response to environmental stresses.

RESULTS

In the K. obovata genome, 66 bZIP genes were identified and named KobZIP1 to KobZIP66, categorized based on their chromosomal locations. These KobZIP genes exhibited diversity in physicochemical properties, such as protein length, molecular weight, and isoelectric point, and were all predicted to localize to the nucleus. Phylogenetic and structural analyses classified the KobZIP genes into 12 subfamilies, with subfamily A containing the majority of members. Gene structure analysis revealed variations in the number and position of exons and introns among subfamilies, reflecting their evolutionary history and potential functional diversity. Conserved motif analysis showed that all bZIP family members contained motifs in the basic and hinge regions, with subfamily D displaying the greatest motif diversity. Promoter region analysis identified various cis-acting elements associated with responses to phytohormones (ABA, GA, ET, IAA, MeJA, SA) and environmental stress. The expression patterns of KobZIP genes across different tissues and under various abiotic stress conditions were analyzed using transcriptomic data and experimental validation.

CONCLUSION

This study provides a comprehensive characterization and functional analysis of the bZIP gene family in K. obovata, offering new insights into their roles in plant development and environmental adaptation. The expression profiles of KobZIP genes during root development and post-embryonic stages, along with their responses to ABA, low temperature, and salt stress, underscore their potential significance in the adaptation of mangrove plants to the intertidal environment.

Evolutionary signatures of the erosion of sexual reproduction genes in domesticated cassava (Manihot esculenta)

Centuries of clonal propagation in cassava (Manihot esculenta) have reduced sexual recombination, leading to the accumulation of deleterious mutations. This has resulted in both inbreeding depression affecting yield and a significant decrease in reproductive performance, creating hurdles for contemporary breeding programs. Cassava is a member of the Euphorbiaceae family, including notable species such as rubber tree (Hevea brasiliensis) and poinsettia (Euphorbia pulcherrima). Expanding upon preliminary draft genomes, we annotated 7 long-read genome assemblies and aligned a total of 52 genomes, to analyze selection across the genome and the phylogeny. Through this comparative genomic approach, we identified 48 genes under relaxed selection in cassava. Notably, we discovered an overrepresentation of floral expressed genes, especially focused at 6 pollen-related genes. Our results indicate that domestication and a transition to clonal propagation have reduced selection pressures on sexually reproductive functions in cassava leading to an accumulation of mutations in pollen-related genes. This relaxed selection and the genome-wide deleterious mutations responsible for inbreeding depression are potential targets for improving cassava breeding, where the generation of new varieties relies on recombining favorable alleles through sexual reproduction.

Comparative transcriptomic and phenotypic analysis of monoclonal and polyclonal Populus deltoides genotypes

Populus species are highly valued for bioenergy and bioproducts due to their rapid growth and productivity. Polyclonal plantings, or mixtures of Populus clones, have shown the potential to enhance resource utilization and productivity, likely due to phenotypic differences arising from niche differentiation. In this study, we investigated gene expression and productivity in monoclonal and polyclonal stands of P. deltoides. Phenotypic results showed that polyclonal plots exhibited higher leaf area index (LAI; p < 0.01, 2.96 ± 0.057 m2) and total biomass (p < 0.01, 2.74 ± 0.06) compared to monoclonal plots, indicating superior productivity. RNA sequencing revealed upregulation of key genes such as exocyst subunit exo70 family protein H7 (EXO70H7), NDH-dependent cyclic electron flow 5 (NDF5), and expansin-like A3 (EXLA3). We also observed enrichment in phenylalanine metabolism and other secondary metabolic pathways in clone S7C8. Phenotypic results, upregulated genes and enriched biological pathways identified in this study may explain the enhanced productivity, increased nitrate content, and expanded canopy in polyclonal plantings. Overall, this study provides a foundation for future research to enhance forest productivity by linking molecular mechanisms to practical applications in field plantings.

Comprehensive analysis of the multi-rings mitochondrial genome of Populus tomentosa

BACKGROUND

Populus tomentosa, known as Chinese white poplar, is indigenous and distributed across large areas of China, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. However, limited accessibility to the mitochondrial (mt) genome of P. tomentosa impedes phylogenetic and population genetic analyses and restricts functional gene research in Salicaceae family.

RESULTS

Single-molecule real-time (SMRT) sequencing technology was used to sequence, assemble, and annotate the mt genome of P. tomentosa. This genome has a complex structure composed of four circular molecules ranging from 153,004 to 330,873 base pairs (bp). Each of these four circular molecules contains unique gene sequences that constitute the mt genome of P. tomentosa. The mt genome comprises 69 functional genes, including 38 protein-coding genes (PCGs), 26 tRNA genes, and 5 rRNA genes. After removing duplications, 19 different tRNA coding genes remain, though only 10 amino acids can be recognized. The noncoding region constitutes 93.38% of the mt genome, comprising a large number of repetitive sequences, gene spacer regions, and insertion from chloroplast sequences. Specifically, 40 chloroplast-derived sequences, with a total length of 24,381 bp, were identified in P. tomentosa.

CONCLUSIONS

In the current study, the results provide mitochondrial genomic evidence for the maternal origin of P. tomentosa and enhance understanding of the gene dialog between organelle genomes, contributing to the conservation and utilization of the genetic resources of P. tomentosa.

Comparative analysis of HKTs in six poplar species and functional characterization of PyHKTs in stress-affected tissues

Plant HKTs (High-affinity K+ transporters) are essential transporters for ion transport and homeostasis and play crucial roles in plant growth and stress responses. However, the evolution of HKTs in Populus species and their functions require further investigation. In this study, we identified 16 HKTs from six Populus species. All poplar HKTs were classified as Class I HKTs because of their physiological relationships and the conservation of amino acids in key structures, which aligns with their conserved evolutionary coding sequences. The analysis of the protein domains, motifs and gene structures of 16 poplar HKTs revealed consistent conservation, with the exception of two members. The number of homologs and their chromosome locations indicated the differentiation of HKTs during poplar evolution and adaptation. Poplar HKTs can be classified into two subgroups on the basis of their physiological relationships and distinct protein structures. Gene expression pattern analysis revealed that poplar HKTs presented relatively high expression levels in roots and stems under salt stress. Furthermore, cis-element analysis and protein interaction predictions provide insights into the functions of HKTs under salt stress through the activation of ion transporters, proline content, and ATPases regulated by hormonal signals and MYB transcription factors. In conclusion, our research established a theoretical framework for investigating the evolutionary relationships and functional roles of HKTs in Populus species and offered valuable insights into the functions and underlying mechanisms of poplar HKTs in specific tissues under various stress conditions.

Identification of the xyloglucan endotransglycosylase/hydrolase genes and the role of PagXTH12 in drought resistance in poplar

The xyloglucan endotransglycosylase/hydrolase (XTH) gene family plays a crucial role in plant cell wall remodeling, facilitating growth and structural changes. However, the divergence of paralogous genes among different species of Populus remains inadequately understood. This study investigates the phylogenetic relationships and expression characteristics of XTH genes in two Populus species: Populus trichocarpa and Populus alba × P. glandulosa '84K'. Forty-one XTHs were identified in P. trichocarpa and 38 and 33 members in the subgenome A and G of '84K' poplar, respectively. Gene expression analysis demonstrated differences among paralogous genes within the same subgenome and between orthologous genes across species, likely influenced by variations in promoter regions. Notably, XTH12 showed a specific response to drought stress among various abiotic stresses. In a population of 549 Populus individuals, functional SNPs in XTH12's coding region did not affect its conserved ExDxE catalytic site, highlighting its irreplaceable function. Furthermore, validation through qRT-PCR and ProPagXTH12::GUS activity, alongside PagXTH12-overexpression poplar lines, substantiated the role of PagXTH12 in modulating the balance between plant biomass and drought resistance. Overall, this research provides valuable insights into the biological functions of XTHs in plant environmental adaptability and offers strategies for targeted regulation of tree growth and stress resistance.

The haplotype-phased genome assembly facilitated the deciphering of the bud dormancy-related QTLs in Prunus mume

Bud dormancy is a vital physiological process in woody perennials, facilitating their adaptation to seasonal environmental changes. Satisfying genotype-specific chilling requirements (CR) and heat requirements (HR) through exposure to specific chilling and warm temperatures is essential for dormancy release and the subsequent resumption of growth. The genetic mechanisms regulating bud dormancy traits in Prunus mume remain unclear. In this study, we first assembled the genome of 'Nanko', the leading P. mume cultivar in Japan, in a haplotype-resolved manner. Using an F1 segregating population from a cross between 'Nanko' (high-chill) and 'SC' (low-chill), a cultivar adapted to subtropical conditions, we identified quantitative trait loci (QTLs) for vegetative bud dormancy traits on chromosome 4 (LG4 QTLs) in the 'Nanko' genome and for CR and HR on chromosome 7 (LG7 QTL) in the 'SC' genome. A notable 5.6 Mb chromosome inversion was overlapped with LG4 QTL interval in one of the 'Nanko' haplotypes. We also identified candidate genes based on haplotyping, differential expression between the parents or the presence of trait-correlated variants in coding regions. Notably, genes such as PmuMAIN, PmuNAC2, PmuDOG1, PmuSUI1, PmuATG8CL, PmubZIP44, and PmuSAUR50 were identified. This study provides valuable insights into the genetic regulation of vegetative bud dormancy in Prunus species.

High-quality genome of Firmiana hainanensis provides insights into the evolution of Malvaceae subfamilies and the mechanism of their wood density formation

The Malvaceae family, the most diverse family in the order Malvales, consists of nine subfamilies. Within the Firmiana genus of the Sterculioideae subfamily, most species are considered globally vulnerable, yet their genomes remain unexplored. Here, we present a chromosome-level genome assembly for a representative Firmiana species, F. hainanensis, 2n = 40, totaling 1536 Mb. Phylogenomic analysis shows that F. hainanensis and Durio zibethinus have the closest evolutionary relationship, with an estimated divergence time of approximately 21 MYA and distinct polyploidization events in their histories. Evolutionary trajectory analyses indicate that fissions and fusions may play a crucial role in chromosome number variation (2n = 14 to 2n = 96). Analysis of repetitive elements among Malvaceae reveals that the Tekay subfamily (belonging to the Gypsy group) contributes to variation in genome size (ranging from 324 Mb to 1620 Mb). Additionally, genes associated with P450, peroxidase, and microtubules, and thereby related to cell wall biosynthesis, are significantly contracted in F. hainanensis, potentially leading to its lower wood density relative to Hopea hainanensis. Overall, our study provides insights into the evolution of chromosome number, genome size, and the genetic basis of cell wall biosynthesis in Malvaceae species.

Genome-wide identification and expression analysis of the BAHD gene family in Leonurus japonicus

Acylation represents a pivotal biochemical process that is instrumental in the modification of secondary metabolites throughout the growth and developmental stages of plants. The BAHD acyltransferase family within the plant kingdom predominantly utilizes coenzyme A thioester as the acyl donor, while employing alcohol or amine compounds as the acceptor substrates to facilitate acylation reactions. Using bioinformatics approaches, the LjBAHD gene family members in the genome of Leonurus japonicus (L. japonicus) were identified and characterized including gene structure, conserved motifs, cis-acting elements, and potential gene functions. To elucidate the roles of BAHD genes in various tissues of L. japonicus, the expression profiles of LjBAHD family members across different organs were scrutinized. Under drought stress treatment, some LjBAHDs were upregulation, suggesting their potential involvement in drought response. Notably, a detailed study was conducted on a specific HCT gene (i.e., LjBAHD25) within the BAHD gene family. Analysis of its expression patterns suggested a role for LjBAHD25 in the phenylpropanoid metabolism pathway in L. japonicus, contributing to the biosynthesis of secondary metabolites with unique bioactivity. The findings of this study have established a scientific foundation for the subsequent development and functional validation of the BAHD gene family in L. japonicus.

Evolution and functional characterization of Populus salt stress-responsive calcineurin B-like protein-interacting protein kinases

KEY MESSAGE

Identification of salt-responsive calcineurin B-like protein-interacting protein kinases (CIPKs) in Populus. Calcineurin B-like protein-interacting protein kinases (CIPKs) play vital roles in plant growth and abiotic stress responses. Currently, the regulatory mechanisms underlying these processes mediated by CIPK proteins are not completely understood in woody species. This study provided the first systematic analysis of 31 Populus CIPK genes and investigated their evolutionary relationships, gene structures, motif compositions, and salt stress responses. A total of 11 pairs of paralogous PtCIPK genes were identified, of which three pairs may be resulted from whole genome duplication, and two pairs that may be created by tandem duplications. RT-qPCR analysis revealed that 93.5% (29/31) genes showed altered expression levels in roots after salt treatment. Ectopic expression of PdCIPK21 or PdCIPK31 in Arabidopsis resulted in significant increases of seed germination, root elongation and fresh weight under salt stress conditions. Cytological observation revealed that PdCIPK21/31 overexpression lines showed increased number, lumen area and cell wall thickness of xylem vessels, and higher lignin content in stems compared with the wild type, with decreased sensitivity to long-term salt stress treatment. Our results suggest that PdCIPK21/31 serve as candidate genes for improving wood production and enhancing salt tolerance of tree species.

Unraveling site-specific seed formation abnormalities in Picea neoveitchii Mast. trees via widely metabolomic and transcriptomic analysis

Picea neoveitchii Mast. is a rare and threatened species of evergreen coniferous tree in China, commonly facing issues such as damaged seeds, abnormal seed growth, and empty seed shells. These abnormalities vary by location; unfortunately, the reasons behind these inconsistencies are completely unknown. This study compared seeds from two 150-year-old trees located in Taibai (Shaanxi province, TB150) and Zhouqu (Gansu province, ZQ150). The results showed significant differences in 43 metabolites and hormone levels, with higher levels of indole-3-acetic acid (IAA), methyl jasmonate (MeJA), and brassinosteroid (BR) in ZQ150, which were associated with more viable seeds. In contrast, TB150 exhibited more damaged seeds and empty seed shells due to higher abscisic acid (ABA) levels. Moreover, to further investigate these inconsistencies, we performed de-novo transcriptomic assembly and functional annotation of unigenes using high-throughput sequencing. A total of 2,355 differentially expressed unigenes were identified between TB150 and ZQ150, with 1,280 upregulated and 1,075 downregulated. Hormone signaling and sugar metabolism-related unigenes were further examined for their role in seed development. ZQ150 increased the number of normal seeds by enhancing endogenous IAA levels and upregulating auxin signaling and sugar metabolism-related genes. Conversely, TB150 showed more empty seed shells, correlated with elevated ABA levels and the activation of ABA signaling genes. We hypothesize that enhanced IAA levels and the upregulation of sugar metabolism and auxin signaling genes promote normal seed development.

COBRA-LIKE4 modulates cellulose synthase velocity and facilitates cellulose deposition in the secondary cell wall

Cellulose is a critical component of secondary cell walls (CWs) and woody tissues of plants. Cellulose synthase (CESA) complexes (CSCs) produce cellulose as they move within the plasma membrane, extruding glucan chains into the CW that coalesce and often crystallize into cellulose fibrils. Here we examine COBRA-LIKE4 (COBL4), a GPI-anchored protein on the outer leaflet of the plasma membrane that is required for normal cellulose deposition in secondary CWs. Characterization of the Arabidopsis (Arabidopsis thaliana) cobl4 mutant alleles called irregular xylem6, irx6-2 and irx6-3, showed reduced α-cellulose content and lower crystallinity, supporting a role for COBL4 in maintaining cellulose quantity and quality. In live-cell imaging, mNeon Green-tagged CESA7 moved in the plasma membrane at higher speeds in the irx6-2 background compared to wild-type. To test conservation of COBL4 function between herbaceous and woody plants, poplar (Populus trichocarpa) COBL4 homologs PtCOBL4a and PtCOBL4b were transformed into, and rescued, the Arabidopsis irx6 mutants. Using the Arabidopsis secondary CW-inducible VND7-GR system to study poplar COBL4 dynamics, YFP-tagged PtCOBL4a localized to the plasma membrane in regions of high cellulose deposition in secondary CW bands. As predicted for a lipid-linked protein, COBL4 was more mobile in the plane of the plasma membrane than CESA7 or a control plasma membrane marker. Following programmed cell death, COBL4 anchored to the secondary CW bands. These data support a role for COBL4 as a modulator of cellulose organization in the secondary CW, influencing cellulose production, and CSC velocity at the plasma membrane.

Genome-Wide Identification of NLP Gene Families and Haplotype Analysis of SiNLP2 in Foxtail Millet (Setaria italica)

Nitrogen is a critical factor in plant growth, development, and crop yield. NODULE-INCEPTION-like proteins (NLPs), which are plant-specific transcription factors, function as nitrate sensors and play a vital role in the nitrogen response of plants. However, the genome-wide identification of the NLP gene family, the elucidation of the underlying molecular mechanism governing nitrogen response, and haplotype mining remain elusive in millet. In this study, we identified seven members of the NLP gene family in the millet genome and systematically analyzed their physicochemical properties. Evolutionary tree analysis indicated that SiNLP members can be classified into three subgroups, with NLP members from the same species preferentially grouped together within each subgroup. Analysis of gene structure characteristics revealed that all SiNLP members contained 10 conserved motifs, as well as the RWP-RK and PB1 domains, indicating that these motifs and domains have been relatively conserved throughout evolution. Additionally, we identified a significant abundance of response elements related to hormones, stress, growth, and development within the promoter regions of SiNLP members, suggesting that these members are involved in regulating diverse physiological processes in millet. Transcriptome data under low-nitrogen conditions showed significant differences in the expression profiles of SiNLP2 and SiNLP4 compared to the other members. RNA-seq and qRT-PCR results demonstrated that SiNLP2 significantly responds to low-nitrogen stress. Notably, we found that SiNLP2 is involved in nitrogen pathways by regulating the expression of the SiNAR2.1A, SiNAR2.1B, SiNRT1.1, and SiNR2 genes. More importantly, we identified an elite haplotype, Hap2, of SiNLP2, which is gradually being utilized in the breeding process. Our study established a foundation for a comprehensive understanding of the SiNLP gene family and provided gene resources for variety improvement and marker-assisted selection breeding.