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

Chromosome-scale genome assembly of sweet tea (Lithocarpus polystachyus Rehder)

Lithocarpus, with >320 species, is the second largest genus of Fagaceae. However, the lack of a reference genome limits the molecular biology and functional study of Lithocarpus species. Here, we report the chromosome-scale genome assembly of sweet tea (Lithocarpus polystachyus Rehder), the first Lithocarpus species to be sequenced to date. Sweet tea has a 952-Mb genome, with a 21.4-Mb contig N50 value and 98.6% complete BUSCO score. In addition, the per-base consensus accuracy and completeness of the genome were estimated at 60.6 and 81.4, respectively. Genome annotation predicted 37,396 protein-coding genes, with repetitive sequences accounting for 64.2% of the genome. The genome did not undergo whole-genome duplication after the gamma (γ) hexaploidy event. Phylogenetic analysis showed that sweet tea diverged from the genus Quercus approximately at 59 million years ago. The high-quality genome assembly and gene annotation resources enrich the genomics of sweet tea, and will facilitate functional genomic studies in sweet tea and other Fagaceae species.

Three-dimensional chromatin architecture in plants - General features and novelties

Research on the three-dimensional (3D) structure of the genome and its distribution within the nuclear space has made a big leap in the last two decades. Work in the animal field has led to significant advances in our general understanding on eukaryotic genome organization. This did not only bring along insights into how the 3D genome interacts with the epigenetic landscape and the transcriptional machinery but also how 3D genome architecture is relevant for fundamental developmental processes, such as cell differentiation. In parallel, the 3D organization of plant genomes have been extensively studied, which resulted in both congruent and novel findings, contributing to a more complete view on how eukaryotic genomes are organized in multiple dimensions. Plant genomes are remarkably diverse in size, composition, and ploidy. Furthermore, as intrinsically sessile organisms without the possibility to relocate to more favorable environments, plants have evolved an elaborate epigenetic repertoire to rapidly respond to environmental challenges. The diversity in genome organization and the complex epigenetic programs make plants ideal study subjects to acquire a better understanding on universal features and inherent constraints of genome organization. Furthermore, considering a wide range of species allows us to study the evolutionary crosstalk between the various levels of genome architecture. In this article, we aim at summarizing important findings on 3D genome architecture obtained in various plant species. These findings cover many aspects of 3D genome organization on a wide range of levels, from gene loops to topologically associated domains and to global 3D chromosome configurations. We present an overview on plant 3D genome organizational features that resemble those in animals and highlight facets that have only been observed in plants to date.

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

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

Above- and belowground fungal biodiversity of Populus trees on a continental scale

Understanding drivers of terrestrial fungal communities over large scales is an important challenge for predicting the fate of ecosystems under climate change and providing critical ecological context for bioengineering plant-microbe interactions in model systems. We conducted an extensive molecular and microscopy field study across the contiguous United States measuring natural variation in the Populus fungal microbiome among tree species, plant niche compartments and key symbionts. Our results show clear biodiversity hotspots and regional endemism of Populus-associated fungal communities explained by a combination of climate, soil and geographic factors. Modelling climate change impacts showed a deterioration of Populus mycorrhizal associations and an increase in potentially pathogenic foliar endophyte diversity and prevalence. Geographic differences among these symbiont groups in their sensitivity to environmental change are likely to influence broader forest health and ecosystem function. This dataset provides an above- and belowground atlas of Populus fungal biodiversity at a continental scale.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Integrating genomic and multiomic data for Angelica sinensis provides insights into the evolution and biosynthesis of pharmaceutically bioactive compounds

Angelica sinensis roots (Angelica roots) are rich in many bioactive compounds, including phthalides, coumarins, lignans, and terpenoids. However, the molecular bases for their biosynthesis are still poorly understood. Here, an improved chromosome-scale genome for A. sinensis var. Qinggui1 is reported, with a size of 2.16 Gb, contig N50 of 4.96 Mb and scaffold N50 of 198.27 Mb, covering 99.8% of the estimated genome. Additionally, by integrating genome sequencing, metabolomic profiling, and transcriptome analysis of normally growing and early-flowering Angelica roots that exhibit dramatically different metabolite profiles, the pathways and critical metabolic genes for the biosynthesis of these major bioactive components in Angelica roots have been deciphered. Multiomic analyses have also revealed the evolution and regulation of key metabolic genes for the biosynthesis of pharmaceutically bioactive components; in particular, TPSs for terpenoid volatiles, ACCs for malonyl CoA, PKSs for phthalide, and PTs for coumarin biosynthesis were expanded in the A. sinensis genome. These findings provide new insights into the biosynthesis of pharmaceutically important compounds in Angelica roots for exploration of synthetic biology and genetic improvement of herbal quality.

A multitiered haplotype strategy to enhance phased assembly and fine mapping of a disease resistance locus

Fine mapping of quantitative trait loci (QTL) to dissect the genetic basis of traits of interest is essential to modern breeding practice. Here, we employed a multitiered haplotypic marker system to increase fine mapping accuracy by constructing a chromosome-level, haplotype-resolved parental genome, accurate detection of recombination sites, and allele-specific characterization of the transcriptome. In the first tier of this system, we applied the preexisting panel of 2,000 rhAmpSeq core genome markers that is transferable across the entire Vitis genus and provides a genomic resolution of 200 kb to 1 Mb. The second tier consisted of high-density haplotypic markers generated from Illumina skim sequencing data for samples enriched for relevant recombinations, increasing the potential resolution to hundreds of base pairs. We used this approach to dissect a novel Resistance to Plasmopara viticola-33 (RPV33) locus conferring resistance to grapevine downy mildew, narrowing the candidate region to only 0.46 Mb. In the third tier, we used allele-specific RNA-seq analysis to identify a cluster of 3 putative disease resistance RPP13-like protein 2 genes located tandemly in a nonsyntenic insertion as candidates for the disease resistance trait. In addition, combining the rhAmpSeq core genome haplotype markers and skim sequencing-derived high-density haplotype markers enabled chromosomal-level scaffolding and phasing of the grape Vitis × doaniana 'PI 588149' assembly, initially built solely from Pacific Biosciences (PacBio) high-fidelity (HiFi) reads, leading to the correction of 16 large-scale phasing errors. Our mapping strategy integrates high-density, phased genetic information with individual reference genomes to pinpoint the genetic basis of QTLs and will likely be widely adopted in highly heterozygous species.

Application and prospects of single-cell and spatial omics technologies in woody plants

Over the past decade, high-throughput sequencing and high-resolution single-cell transcriptome sequencing technologies have undergone rapid development, leading to significant breakthroughs. Traditional molecular biology methods are limited in their ability to unravel cellular-level heterogeneity within woody plant tissues. Consequently, techniques such as single-cell transcriptomics, single-cell epigenetics, and spatial transcriptomics are rapidly gaining popularity in the study of woody plants. In this review, we provide a comprehensive overview of the development of these technologies, with a focus on their applications and the challenges they present in single-cell transcriptome research in woody plants. In particular, we delve into the similarities and differences among the results of current studies and analyze the reasons behind these differences. Furthermore, we put forth potential solutions to overcome the challenges encountered in single-cell transcriptome applications in woody plants. Finally, we discuss the application directions of these techniques to address key challenges in woody plant research in the future.

Identification of Functional Brassinosteroid Receptor Genes in Oaks and Functional Analysis of QmBRI1

Brassinosteroids (BRs) play important regulatory roles in plant growth and development, with functional BR receptors being crucial for BR recognition or signaling. Although functional BR receptors have been extensively studied in herbaceous plants, they remain largely under-studied in forest tree species. In this study, nine BR receptors were identified in three representative oak species, of which BRI1s and BRL1s were functional BR receptors. Dispersed duplications were a driving force for oak BR receptor expansion, among which the Brassinosteroid-Insensitive-1 (BRI1)-type genes diverged evolutionarily from most rosids. In oak BRI1s, we identified that methionine in the conserved Asn-Gly-Ser-Met (NGSM) motif was replaced by isoleucine and that the amino acid mutation occurred after the divergence of Quercus and Fagus. Compared with QmBRL1, QmBRI1 was relatively highly expressed during BR-induced xylem differentiation and in young leaves, shoots, and the phloem and xylem of young stems of Quercus mongolica. Based on Arabidopsis complementation experiments, we proved the important role of QmBRI1 in oak growth and development, especially in vascular patterning and xylem differentiation. These findings serve as an important supplement to the findings of the structural, functional and evolutionary studies on functional BR receptors in woody plants and provide the first example of natural mutation occurring in the conserved BR-binding region (NGSM motif) of angiosperm BRI1s.

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

In vitro induction of tetraploidy and its effects on phenotypic variations in Populus hopeiensis

BACKGROUND

Artificial induction of polyploidy is the most common and effective way to improve the biological properties of Populus and develop new varieties of this tree. In this study, in order to confirm and expand earlier findings, we established a protocol using colchicine and based on an efficient shoot regeneration system of leaf blades to induce tetraploidy in vitro in three genotypes from diploid Populus hopeiensis. The stomatal characteristics, leaf blade size, and growth were evaluated for diploids and tetraploids of three genotypes.

RESULTS

We found that genotype, preculture duration, colchicine concentration, and colchicine exposure time had highly significant effects on the tetraploid induction rate. The optimal protocol for inducing tetraploidy in P. hopeiensis was to preculture leaf blades for 7 days and then treat them for 4 days with 40 mg/L colchicine. The tetraploid induction rates of genotypes BT1, BT3, and BT8 were 21.2, 11.4 and 16.7%, respectively. A total of 136 tetraploids were identified by flow cytometry analysis and somatic chromosome counting. The stomatal length, width, and density of leaf blades significantly differed between diploid and tetraploid plants. Compared with their diploid counterparts, the tetraploids produced larger leaf blades and had a slower growth rate. Our findings further document the modified morphological characteristics of P. hopeiensis following whole-genome duplication (e.g., induced tetraploidy).

CONCLUSIONS

We established a protocol for in vitro induction of tetraploidy from diploid leaf blades treated with colchicine, which can be applied to different genotypes of P. hopeiensis.

Evolution and diversification of the O-methyltransferase (OMT) gene family in Solanaceae

O-methyltransferases (OMTs) are a group of enzymes involved in several fundamental biological processes in plants, including lignin biosynthesis, pigmentation, and aroma production. Despite the intensive investigation of the role of OMTs in plant secondary metabolism, the evolution and diversification of this gene family in Solanaceae remain poorly understood. Here, we conducted a genome-wide survey of OMT genes in six Solanaceae species, reconstructing gene phylogenetic trees, predicting the potential involvement in biological processes, and investigating the exon/intron structure and chromosomal location. We identified 57 caffeoyl-CoA OMTs (CCoAOMTs) and 196 caffeic acid OMTs (COMTs) in the studied species. We observed a conserved gene block on chromosome 2 that consisted of tandem duplicated copies of OMT genes. Our results suggest that the expansion of the OMT gene family in Solanaceae was driven by whole genome duplication, segmental duplication, and tandem duplication, with multiple genes being retained by neofunctionalization and subfunctionalization. This study represents an essential first step in unraveling the evolutionary history of OMTs in Solanaceae. Our findings deepen our understanding of the crucial role of OMTs in several biological processes and highlight their significance as potential biotechnological targets.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

Heterologous expression of PtAAS1 reveals the metabolic potential of the common plant metabolite phenylacetaldehyde for auxin synthesis in planta

Aromatic aldehydes and amines are common plant metabolites involved in several specialized metabolite biosynthesis pathways. Recently, we showed that the aromatic aldehyde synthase PtAAS1 and the aromatic amino acid decarboxylase PtAADC1 contribute to the herbivory-induced formation of volatile 2-phenylethanol and its glucoside 2-phenylethyl-β-D-glucopyranoside in Populus trichocarpa. To unravel alternative metabolic fates of phenylacetaldehyde and 2-phenylethylamine beyond alcohol and alcohol glucoside formation, we heterologously expressed PtAAS1 and PtAADC1 in Nicotiana benthamiana and analyzed plant extracts using untargeted LC-qTOF-MS and targeted LC-MS/MS analysis. While the metabolomes of PtAADC1-expressing plants did not significantly differ from those of control plants, expression of PtAAS1 resulted in the accumulation of phenylacetic acid (PAA) and PAA-amino acid conjugates, identified as PAA-aspartate and PAA-glutamate. Herbivory-damaged poplar leaves revealed significantly induced accumulation of PAA-Asp, while levels of PAA remained unaltered upon herbivory. Transcriptome analysis showed that members of auxin-amido synthetase GH3 genes involved in the conjugation of auxins with amino acids were significantly upregulated upon herbivory in P. trichocarpa leaves. Overall, our data indicates that phenylacetaldehyde generated by poplar PtAAS1 serves as a hub metabolite linking the biosynthesis of volatile, non-volatile herbivory-induced specialized metabolites, and phytohormones, suggesting that plant growth and defense can be balanced on a metabolic level.

Functional analysis of Salix purpurea genes support roles for ARR17 and GATA15 as master regulators of sex determination

The Salicaceae family is of growing interest in the study of dioecy in plants because the sex determination region (SDR) has been shown to be highly dynamic, with differing locations and heterogametic systems between species. Without the ability to transform and regenerate Salix in tissue culture, previous studies investigating the mechanisms regulating sex in the genus Salix have been limited to genome resequencing and differential gene expression, which are mostly descriptive in nature, and functional validation of candidate sex determination genes has not yet been conducted. Here, we used Arabidopsis to functionally characterize a suite of previously identified candidate genes involved in sex determination and sex dimorphism in the bioenergy shrub willow Salix purpurea. Six candidate master regulator genes for sex determination were heterologously expressed in Arabidopsis, followed by floral proteome analysis. In addition, 11 transcription factors with predicted roles in mediating sex dimorphism downstream of the SDR were tested using DAP-Seq in both male and female S. purpurea DNA. The results of this study provide further evidence to support models for the roles of ARR17 and GATA15 as master regulator genes of sex determination in S. purpurea, contributing to a regulatory system that is notably different from that of its sister genus Populus. Evidence was also obtained for the roles of two transcription factors, an AP2/ERF family gene and a homeodomain-like transcription factor, in downstream regulation of sex dimorphism.

Genome-wide identification, characterization, and expression analysis of CCT transcription factors in poplar

The CCT [CONSTANS (CO), CO-like, and TIMING OF CAB EXPRESSION1 (TOC1)] gene family is involved in photoperiodic flowering and adaptation to different environments. In this study, 39 CCT family genes from the poplar genome were identified and characterized, including 18 COL, 7 PRR, and 14 CMF TFs. Phylogenetics analysis showed that the PtrCCT gene family could be classified into five classes (Classes I-V) that have close relationships with Arabidopsis thaliana. Eight pairs of PtrCCTs had collinear relationships through interchromosomal synteny analysis in poplar, suggesting segmental duplication played a vital role in the expansion of the poplar CCT gene family. Besides, synteny analyses of the CCT members among poplar and different species provided more clues for PtrCCT gene family evolution. Cis-acting elements in the promoters of PtrCCTs predicted their involvement in light responses, hormone responses, biotic/abiotic stress responses, and plant growth and development. Eight members of the PpnCCT gene family were differentially expressed in the apical buds and leaves of triploid poplar compared to diploids. We then focused on PpnCCT39 upregulated in triploid poplars and showed that PpnCCT39 was localized in the nucleus, chloroplast, and cytoplasm and could interact with CLPP1 in the chloroplast. Overexpression of PpnCCT39 in poplar increased chlorophyll contents and enhanced photosynthetic rate. This study provided comprehensive information for the CCT gene family and set up a basis for its function identification in poplar.

Engineering custom morpho- and chemotypes of Populus for sustainable production of biofuels, bioproducts, and biomaterials

Humans have been modifying plant traits for thousands of years, first through selection (i.e., domestication) then modern breeding, and in the last 30 years, through biotechnology. These modifications have resulted in increased yield, more efficient agronomic practices, and enhanced quality traits. Precision knowledge of gene regulation and function through high-resolution single-cell omics technologies, coupled with the ability to engineer plant genomes at the DNA sequence, chromatin accessibility, and gene expression levels, can enable engineering of complex and complementary traits at the biosystem level. Populus spp., the primary genetic model system for woody perennials, are among the fastest growing trees in temperate zones and are important for both carbon sequestration and global carbon cycling. Ample genomic and transcriptomic resources for poplar are available including emerging single-cell omics datasets. To expand use of poplar outside of valorization of woody biomass, chassis with novel morphotypes in which stem branching and tree height are modified can be fabricated thereby leading to trees with altered leaf to wood ratios. These morphotypes can then be engineered into customized chemotypes that produce high value biofuels, bioproducts, and biomaterials not only in specific organs but also in a cell-type-specific manner. For example, the recent discovery of triterpene production in poplar leaf trichomes can be exploited using cell-type specific regulatory sequences to synthesize high value terpenes such as the jet fuel precursor bisabolene specifically in the trichomes. By spatially and temporally controlling expression, not only can pools of abundant precursors be exploited but engineered molecules can be sequestered in discrete cell structures in the leaf. The structural diversity of the hemicellulose xylan is a barrier to fully utilizing lignocellulose in biomaterial production and by leveraging cell-type-specific omics data, cell wall composition can be modified in a tailored and targeted specific manner to generate poplar wood with novel chemical features that are amenable for processing or advanced manufacturing. Precision engineering poplar as a multi-purpose sustainable feedstock highlights how genome engineering can be used to re-imagine a crop species.

Hidden prevalence of deletion-inversion bi-alleles in CRISPR-mediated deletions of tandemly arrayed genes in plants

Tandemly arrayed genes (TAGs) with functional redundancy and chromosomal linkage constitute 14 ~ 35% in sequenced plant genomes. The multiplex CRISPR system is the tool of choice for creating targeted TAG deletions. Here, we show that up to ~80% of CRISPR-mediated TAG knockout alleles in Arabidopsis and rice are deletion-inversion (delinver) bi-alleles, which are easily misidentified as homozygous deletion alleles by routine PCR-based genotyping. This can lead to misinterpretation of experimental data and production of progenies with genetic heterogeneity in an unnoticed manner. In ~2,650 transgenic events, delinver mutation frequencies are predominantly correlated with deletion frequencies but unrelated to chromosomal locations or deletion sizes. Delinver mutations also occur frequently at genomic non-TAG loci during multiplexed CRISPR editing. Our work raises the alarm about delinver mutations as common unwanted products of targeted TAG deletions in plants and helps prevent false interpretation of plant TAG functions due to this hidden genotype issue.

Genome-Wide Identification and Expression Analysis of bZIP Family Genes in Stevia rebaudiana

The basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors widely distributed in eukaryotic organisms. In plants, they are not only involved in growth and development, defense and stress responses and regulation of physiological processes but also play a pivotal role in regulating secondary metabolism. To explore the function related to the bZIP gene family in Stevia rebaudiana Bertoni, we identified 105 SrbZIP genes at the genome-wide level and classified them into 12 subfamilies using bioinformation methods. Three main classes of cis-acting elements were found in the SrbZIP promoter regions, including development-related elements, defense and stress-responsive elements and phytohormone-responsive elements. Through protein-protein interaction network of 105 SrbZIP proteins, SrbZIP proteins were mainly classified into four major categories: ABF2/ABF4/ABI5 (SrbZIP51/SrbZIP38/SrbZIP7), involved in phytohormone signaling, GBF1/GBF3/GBF4 (SrbZIP29/SrbZIP63/SrbZIP60) involved in environmental signaling, AREB3 (SrbZIP88), PAN (SrbZIP12), TGA1 (SrbZIP69), TGA4 (SrbZIP82), TGA7 (SrbZIP31), TGA9 (SrbZIP95), TGA10 (SrbZIP79) and HY5 (SrbZIP96) involved in cryptochrome signaling, and FD (SrbZIP72) promoted flowering. The transcriptomic data showed that SrbZIP genes were differentially expressed in six S. rebaudiana cultivars ('023', '110', 'B1188', '11-14', 'GP' and 'GX'). Moreover, the expression levels of selected 15 SrbZIP genes in response to light, abiotic stress (low temperature, salt and drought), phytohormones (methyl jasmonate, gibberellic acid and salicylic acid) treatment and in different tissues were analyzed utilizing qRT-PCR. Some SrbZIP genes were further identified to be highly induced by factors affecting glycoside synthesis. Among them, three SrbZIP genes (SrbZIP54, SrbZIP63 and SrbZIP32) were predicted to be related to stress-responsive terpenoid synthesis in S. rebaudiana. The protein-protein interaction network expanded the potential functions of SrbZIP genes. This study firstly provided the comprehensive genome-wide report of the SrbZIP gene family, laying a foundation for further research on the evolution, function and regulatory role of the bZIP gene family in terpenoid synthesis in S. rebaudiana.

RNAseq based variant dataset in a black poplar association panel

OBJECTIVE

Black poplar (Populus nigra L.) is a species native to Eurasia with a wide distribution area. It is an ecologically important species from riparian ecosystems, that is used as a parent of interspecific (P. deltoides x P. nigra) cultivated poplar hybrids. Variant detection from transcriptomics sequences of 241 P. nigra individuals, sampled in natural populations from 11 river catchments (in four European countries) is described here. These data provide new valuable resources for population structure analysis, population genomics and genome-wide association studies.

DATA DESCRIPTION

We generated transcriptomics data from a mixture of young differentiating xylem and cambium tissues of 480 Populus nigra trees sampled in a common garden experiment located at Orléans (France), corresponding to 241 genotypes (2 clonal replicates per genotype, at maximum) by using RNAseq technology. We launched on the resulting sequences an in-silico pipeline that allowed us to obtain 878,957 biallelic polymorphisms without missing data. More than 99% of these positions are annotated and 98.8% are located on the 19 chromosomes of the P. trichocarpa reference genome. The raw RNAseq sequences are available at the NCBI Sequence Read Archive SPR188754 and the variant dataset at the Recherche Data Gouv repository under https://doi.org/10.15454/8DQXK5 .

Genome-wide identification of GH9 gene family and the assessment of its role during fruit abscission zone formation in Vaccinium ashei

KEY MESSAGE

The genomic location and stage-specific expression pattern of GH9 genes reveal their critical roles during fruit abscission zone formation in Vaccinium ashei. Glycosyl hydrolase family 9 (GH9) cellulases play a crucial role in both cellulose synthesis and hydrolysis during plant growth and development. Despite this importance, there is currently no study on the involvement of GH9-encoding genes, specifically VaGH9s, in abscission zone formation of rabbiteye blueberries (Vaccinium ashei). In this study, we identified a total of 61 VaGH9s in the genome, which can be classified into 3 subclasses based on conserved motifs and domains, gene structures, and phylogenetic analyses. Our synteny analysis revealed that VaGH9s are more closely related to the GH9s of Populus L. than to those of Arabidopsis, Vitis vinifera, and Citrus sinensis. In silico structural analysis predicted that most of VaGH9s are hydrophilic, and localized in cell membrane and/or cell wall, and the variable sets of cis-acting regulatory elements and functional diversity with four categories of stress response, hormone regulation, growth and development, and transcription factor-related elements are present in the promoter sequence of VaGH9s genes. Transcriptomic analysis showed that there were 22 differentially expressed VaGH9s in fruit abscission zone tissue at the veraison stage, and the expression of VaGH9B2 and VaGH9C10 was continuously increased during fruit maturation, which were in parallel with the increasing levels of cellulase activity and oxidative stress indicators, suggesting that they are involved in the separation stage of fruit abscission in Vaccinium ashei. Our work identified 22 VaGH9s potentially involved in different stages of fruit abscission and would aid further investigation into the molecular regulation of abscission in rabbiteye blueberries fruit.

Genome-wide identification and expression analysis of Hsp70 family genes in Cassava (Manihot esculenta Crantz)

Hsp70 proteins function as molecular chaperones, regulating various cellular processes in plants. In this study, a genome-wide analysis led to the identification of 22 Hsp70 (MeHsp70) genes in cassava. Phylogenetic relationship studies with other Malpighiales genomes (Populus trichocarpa, Ricinus communis and Salix purpurea) classified MeHsp70 proteins into eight groups (Ia, Ib, Ic, Id, Ie, If, IIa and IIb). Promoter analysis of MeHsp70 genes revealed the presence of tissue-specific, light, biotic and abiotic stress-responsive cis-regulatory elements showing their functional importance in cassava. Meta-analysis of publically available RNA-seq transcriptome datasets showed constitutive, tissue-specific, biotic and abiotic stress-specific expression patterns among MeHsp70s in cassava. Among 22 Hsp70, six MeHsp70s viz., MecHsp70-3, MecHsp70-6, MeBiP-1, MeBiP-2, MeBiP-3 and MecpHsp70-2 displayed constitutive expression, while three MecHsp70s were induced under both drought and cold stress conditions. Five MeHsp70s, MecHsp70-7, MecHsp70-11, MecHsp70-12, MecHsp70-13, and MecHsp70-14 were induced under drought stress conditions. We predicted that 19 MeHsp70 genes are under the regulation of 24 miRNAs. This comprehensive genome-wide analysis of the Hsp70 gene family in cassava provided valuable insights into their functional roles and identified various potential Hsp70 genes associated with stress tolerance and adaptation to environmental stimuli.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03760-3.

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

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

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

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

Comprehensive Analysis and Characterization of the GATA Gene Family, with Emphasis on the GATA6 Transcription Factor in Poplar

GATA transcription factors are ubiquitously present in eukaryotic organisms and play a crucial role in multiple biological processes, such as plant growth, stress response, and hormone signaling. However, the study of GATA factors in poplar is currently limited to a small number of proteins, despite their evident functional importance. In this investigation, we utilized the most recent genome annotation and stringent criteria to identify 38 GATA transcription factor genes in poplar. Subsequently, we conducted a comprehensive analysis of this gene family, encompassing phylogenetic classification, protein characterization, analysis of promoter cis-acting elements, and determination of chromosomal location. Our examination of gene duplication events indicated that both tandem and segmental duplications have contributed to the expansion of the GATA gene family in poplar, with segmental duplication potentially being a major driving force. By performing collinearity analysis of genes across six different species, we identified 74 pairs of co-linear genes, which provide valuable insights for predicting gene functions from a comparative genomics perspective. Furthermore, through the analysis of gene expression patterns, we identified five GATA genes that exhibited differential expression in leaf-stem-root tissues and eight genes that were responsive to salt stress. Of particular interest was GATA6, which displayed strong induction by salt stress and overlapped between the two gene sets. We discovered that GATA6 encodes a nuclear-localized protein with transcription activation activity, which is continuously induced by salt stress in leaf and root tissues. Moreover, we constructed a co-expression network centered around GATA6, suggesting the potential involvement of these genes in the growth, development, and response to abiotic stress processes in poplar through cell transport systems and protein modification mechanisms, such as vesicle-mediated transport, intracellular transport, ubiquitination, and deubiquitination. This research provides a foundation for further exploration of the functions and mechanisms of GATA transcription factors in poplar.

Integrative analysis of metabolome, proteome, and transcriptome for identifying genes influencing total lignin content in Populus trichocarpa

Lignin, a component of plant cell walls, possesses significant research potential as a renewable energy source to replace carbon-based products and as a notable pollutant in papermaking processes. The monolignol biosynthetic pathway has been elucidated and it is known that not all monolignol genes influence the total lignin content. However, it remains unclear which monolignol genes are more closely related to the total lignin content and which potential genes influence the total lignin content. In this study, we present a combination of t-test, differential gene expression analysis, correlation analysis, and weighted gene co-expression network analysis to identify genes that regulate the total lignin content by utilizing multi-omics data from transgenic knockdowns of the monolignol genes that includes data related to the transcriptome, proteome, and total lignin content. Firstly, it was discovered that enzymes from the PtrPAL, Ptr4CL, PtrC3H, and PtrC4H gene families are more strongly correlated with the total lignin content. Additionally, the co-downregulation of three genes, PtrC3H3, PtrC4H1, and PtrC4H2, had the greatest impact on the total lignin content. Secondly, GO and KEGG analysis of lignin-related modules revealed that the total lignin content is not only influenced by monolignol genes, but also closely related to genes involved in the "glutathione metabolic process", "cellular modified amino acid metabolic process" and "carbohydrate catabolic process" pathways. Finally, the cinnamyl alcohol dehydrogenase genes CAD1, CADL3, and CADL8 emerged as potential contributors to total lignin content. The genes HYR1 (UDP-glycosyltransferase superfamily protein) and UGT71B1 (UDP-glucosyltransferase), exhibiting a close relationship with coumarin, have the potential to influence total lignin content by regulating coumarin metabolism. Additionally, the monolignol genes PtrC3H3, PtrC4H1, and PtrC4H2, which belong to the cytochrome P450 genes, may have a significant impact on the total lignin content. Overall, this study establishes connections between gene expression levels and total lignin content, effectively identifying genes that have a significant impact on total lignin content and offering novel perspectives for future lignin research endeavours.

Optimization of High-Efficiency Tissue Culture Regeneration Systems in Gray Poplar

A series of tissue culture regeneration protocols were conducted on gray poplar (P. tremula × P. alba) to select the most efficient callus induction medium, adventitious shoot differentiation medium, shoot elongation medium and rooting medium, which laid the foundation for the optimization of genetic transformation technology for gray poplar. The results showed that the Woody Plant Medium (WPM) supplemented with 0.10 mg L-1 kinetin (KT) and 1.00 mg L-1 2,4-dichlorophenoxyacetic acid (2,4-D) was the most suitable medium for callus induction. The callus induction rates of different tissues were greater than 85.7%. The optimal adventitious shoot differentiation medium was the WPM supplemented with 0.02 mg L-1 thidiazuron (TDZ), and the adventitious shoot differentiation rates of young tissues were 22.2-41.9%. The optimal direct differentiation medium was the Murashige and Skoog (MS) medium supplemented with 0.20 mg L-1 6-benzylaminopurine (6-BA), 0.10 mg L-1 indole butyric acid (IBA) and 0.001 mg L-1 TDZ, and the differentiation rate of adventitious shoots was greater than 94%. The best shoot elongation medium for adventitious shoots was the MS medium with 0.10 mg L-1 naphthylacetic acid (NAA). After 45 days of cultivation in the MS medium with 0.10 mg L-1 NAA, the average plant height was 1.8 cm, and the average number of elongated adventitious shoots was 11 per explant. The 1/2 MS medium with 0.10 mg L-1 NAA showed the best performance for rooting, and later, shoot growth. The direct shoot induction pathway can induce adventitious shoots much faster than the indirect adventitious shoot induction pathway can, and the time cost via the direct adventitious shoot induction pathway can be shortened by 2-6 weeks compared to that of the indirect shoot induction pathway.

PtrABR1 Increases Tolerance to Drought Stress by Enhancing Lateral Root Formation in Populus trichocarpa

Roots are the main organ for water uptake and the earliest part of a plant's response to drought, making them of great importance to our understanding of the root system's response to drought. However, little is known about the underlying molecular mechanisms that control root responses to drought stress. Here, we identified and functionally characterized the AP2/ERF family transcription factor (TF) PtrABR1 and the upstream target gene zinc-finger protein TF PtrYY1, which respond to drought stress by promoting the growth and development of lateral roots in Populus trichocarpa. A root-specific induction of PtrABR1 under drought stress was explored. The overexpression of PtrABR1 (PtrABR1-OE) promoted root growth and development, thereby increasing tolerance to drought stress. In addition, PtrYY1 is directly bound to the promoter of PtrABR1 under drought stress, and the overexpression of PtrYY1 (PtrYY1-OE) promoted lateral root growth and development and increased tolerance to drought stress. An RNA-seq analysis of PtrABR1-OE with wild-type (WT) poplar identified PtrGH3.6 and PtrPP2C44, which share the same pattern of expression changes as PtrABR1. A qRT-PCR and cis-element analysis further suggested that PtrGH3.6 and PtrPP2C44 may act as potential downstream targets of PtrABR1 genes in the root response pathway to drought stress. In conclusion, these results reveal a novel drought regulatory pathway in which PtrABR1 regulates the network through the upstream target gene PtrYY1 and the potential downstream target genes PtrGH3.6 and PtrPP2C44, thereby promoting root growth and development and improving tolerance to drought stress.

Identification of Differentially Expressed lncRNAs in Response to Blue Light and Expression Pattern Analysis of Populus tomentosa Hybrid Poplar 741

Poplar is an important shelterbelt, timber stand, and city tree species that has been the focus of forestry research. The regulatory role of the long non-coding RNA molecule (lncRNA; length > 200 nt) has been a research hotspot in plants. In this study, seedlings of 741 poplar were irradiated with LED blue and white light, and the Illumina HiSeq 2000 sequencing platform was used to identify lncRNAs. |logFC| > 1 and p < 0.05 were considered to indicate differentially expressed lncRNAs, and nine differentially expressed lncRNAs were screened, the target genes of which were predicted, and three functionally annotated target genes were obtained. The differentially expressed lncRNAs were identified as miRNA targets. Six lncRNAs were determined to be target sites for twelve mRNAs in six miRNA families. LncRNAs and their target genes, including lncRNA MSTRG.20413.1-ptc-miR396e-5p-GRF9, were verified using quantitative real-time polymerase chain reaction analysis, and the expression patterns were analyzed. The analysis showed that the ptc-miR396e-5p expression was downregulated, while lncRNA MSTRG.20413.1 and GRF9 expression was upregulated, after blue light exposure. These results indicate that lncRNAs interact with miRNAs to regulate gene expression and affect plant growth and development.

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

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

Seed Treatment with Electromagnetic Field Induces Different Effects on Emergence, Growth and Profiles of Biochemical Compounds in Seven Half-Sib Families of Silver Birch

In the context of climate change, strategies aimed at enhancing trees' resistance to biotic and abiotic stress are particularly relevant. We applied an electromagnetic field (EMF) seed treatment to observe changes in the establishment and content of biochemical compounds in silver birch seedlings induced by a short (1 min) seed exposure to a physical stressor. The impact of EMF treatment was evaluated on seedling emergence and growth of one-year-old and two-year-old seedlings from seven half-sib families of silver birch. The effects on numerous biochemical parameters in seedling leaves, such as total phenolic content (TPC), total flavonoid content (TFC), amounts of photosynthetic pigments, total soluble sugars (TSS), level of lipid peroxidation level, antioxidant activity and activity of antioxidant enzymes, were compared using spectrophotometric methods. The results indicated that, in one-year-old seedlings, two of seven (60th and 73rd) half-sib families exhibited a positive response to seed treatment with EMFs in nearly all analyzed parameters. For example, in the 60th family, seed treatment with EMFs increased the percentage of emergence by 3 times, one-year-old seedling height by 71%, leaf TPC by 47%, antioxidant activity by 2 times and amount of chlorophyll a by 4.6 times. Meanwhile, the other two (86th and 179th) families exhibited a more obvious positive response to EMF in two-year-old seedlings as compared to one-year-old seedling controls. The results revealed that short-term EMF treatment of silver birch seeds can potentially be used to improve seedling emergence and growth and increase the content of secondary metabolites, antioxidant capacity and photosynthetic pigments. Understanding of the impact of EMFs as well as the influence of genetic differences on tree responses can be significant for practical applications in forestry. Genetic selection of plant genotypes that exhibit positive response trends can open the way to improve the quality of forest stands.

DIENELACTONE HYDROLASE LIKE PROTEIN1 negatively regulates the KAI2-ligand pathway in Marchantia polymorpha

Karrikins are smoke-derived butenolides that induce seed germination and photomorphogenesis in a wide range of plants.1,2,3 KARRIKIN INSENSITIVE2 (KAI2), a paralog of a strigolactone receptor, perceives karrikins or their metabolized products in Arabidopsis thaliana.4,5,6,7 Furthermore, KAI2 is thought to perceive an unidentified plant hormone, called KAI2 ligand (KL).8,9 KL signal is transduced via the interaction between KAI2, MORE AXILLARY GROWTH2 (MAX2), and SUPPRESSOR of MORE AXILLARY GROWTH2 1 LIKE family proteins (SMXLs), followed by the degradation of SMXLs.4,7,10,11,12,13,14 This signaling pathway is conserved both in A. thaliana and the bryophyte Marchantia polymorpha.14 Although the KL signaling pathway is well characterized, the KL metabolism pathways remain poorly understood. Here, we show that DIENELACTONE HYDROLASE LIKE PROTEIN1 (DLP1) is a negative regulator of the KL pathway in M. polymorpha. The KL signal induces DLP1 expression. DLP1 overexpression lines phenocopied the Mpkai2a and Mpmax2 mutants, while dlp1 mutants phenocopied the Mpsmxl mutants. Mutations in the KL signaling genes largely suppressed these phenotypes, indicating that DLP1 acts upstream of the KL signaling pathway, although DLP1 also has KL pathway-independent functions. DLP1 exhibited enzymatic activity toward a potential substrate, suggesting the possibility that DLP1 works through KL inactivation. Investigation of DLP1 homologs in A. thaliana revealed that they do not play a major role in the KL pathway, suggesting different mechanisms for the KL signal regulation. Our findings provide new insights into the regulation of the KL signal in M. polymorpha and the evolution of the KL pathway in land plants.

The complete chloroplast genome of Salix lindleyana (salicaceae), a plateau plant species

Salix lindleyana Wallich ex Andersson 1851 is a species of genus Salix which mainly grows on mountains above 3000 m at sea level in Qinghai-Tibetan Plateau (including the Himalayas and Hengduan Mountains). To determine its phylogenetic position within Salix, we reconstructed S. lindleyana complete chloroplast (cp) genome sequence by de novo assembly using whole-genome sequencing data. The completed chloroplast genome was 155,304 bp, with a total GC content of 36.7%. It had a very typical tetrad structure, including a large single-copy (LSC) region of 84,539 bp, a small single-copy (SSC) region of 16,161 bp, and two inverted repeats (IR) regions of 27,302 bp. A total of 132 functional genes were distributed in the chloroplast genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis showed that S. lindleyana was clustered with Salix dasyclados Wimmer 1849 and Salix variegata Franchet 1887. The complete chloroplast genome of S. lindleyana provides potential genetic resources for further phylogenetic studies.

Biotechnological Potential of the Stress Response and Plant Cell Death Regulators Proteins in the Biofuel Industry

Production of biofuel from lignocellulosic biomass is relatively low due to the limited knowledge about natural cell wall loosening and cellulolytic processes in plants. Industrial separation of cellulose fiber mass from lignin, its saccharification and alcoholic fermentation is still cost-ineffective and environmentally unfriendly. Assuming that the green transformation is inevitable and that new sources of raw materials for biofuels are needed, we decided to study cell death-a natural process occurring in plants in the context of reducing the recalcitrance of lignocellulose for the production of second-generation bioethanol. "Members of the enzyme families responsible for lysigenous aerenchyma formation were identified during the root hypoxia stress in Arabidopsis thaliana cell death mutants. The cell death regulatory genes, LESION SIMULATING DISEASE 1 (LSD1), PHYTOALEXIN DEFICIENT 4 (PAD4) and ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) conditionally regulate the cell wall when suppressed in transgenic aspen. During four years of growth in the field, the following effects were observed: lignin content was reduced, the cellulose fiber polymerization degree increased and the growth itself was unaffected. The wood of transgenic trees was more efficient as a substrate for saccharification, alcoholic fermentation and bioethanol production. The presented results may trigger the development of novel biotechnologies in the biofuel industry.

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis

BACKGROUND

The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis.

RESULTS

In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR.

CONCLUSIONS

This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

Plant myo-inositol transport influences bacterial colonization phenotypes

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

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient

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

GWAS identifies candidate genes controlling adventitious rooting in Populus trichocarpa

Adventitious rooting (AR) is critical to the propagation, breeding, and genetic engineering of trees. The capacity for plants to undergo this process is highly heritable and of a polygenic nature; however, the basis of its genetic variation is largely uncharacterized. To identify genetic regulators of AR, we performed a genome-wide association study (GWAS) using 1148 genotypes of Populus trichocarpa. GWASs are often limited by the abilities of researchers to collect precise phenotype data on a high-throughput scale; to help overcome this limitation, we developed a computer vision system to measure an array of traits related to adventitious root development in poplar, including temporal measures of lateral and basal root length and area. GWAS was performed using multiple methods and significance thresholds to handle non-normal phenotype statistics and to gain statistical power. These analyses yielded a total of 277 unique associations, suggesting that genes that control rooting include regulators of hormone signaling, cell division and structure, reactive oxygen species signaling, and other processes with known roles in root development. Numerous genes with uncharacterized functions and/or cryptic roles were also identified. These candidates provide targets for functional analysis, including physiological and epistatic analyses, to better characterize the complex polygenic regulation of AR.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

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

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

Apoplastic barriers of Populus × canescens roots in reaction to different cultivation conditions and abiotic stress treatments

Populus is an important tree genus frequently cultivated for economical purposes. However, the high sensitivity of poplars towards water deficit, drought, and salt accumulation significantly affects plant productivity and limits biomass yield. Various cultivation and abiotic stress conditions have been described to significantly induce the formation of apoplastic barriers (Casparian bands and suberin lamellae) in roots of different monocotyledonous crop species. Thus, this study aimed to investigate to which degree the roots of the dicotyledonous gray poplar (Populus × canescens) react to a set of selected cultivation conditions (hydroponics, aeroponics, or soil) and abiotic stress treatments (abscisic acid, oxygen deficiency) because a differing stress response could potentially help in explaining the observed higher stress susceptibility. The apoplastic barriers of poplar roots cultivated in different environments were analyzed by means of histochemistry and gas chromatography and compared to the available literature on monocotyledonous crop species. Overall, dicotyledonous poplar roots showed only a remarkably low induction or enhancement of apoplastic barriers in response to the different cultivation conditions and abiotic stress treatments. The genetic optimization (e.g., overexpression of biosynthesis key genes) of the apoplastic barrier development in poplar roots might result in more stress-tolerant cultivars in the future.

Transcript-wide identification and expression pattern analysis to comprehend the roles of AP2/ERF genes under development and abiotic stress in Trichosanthes kirilowii

BACKGROUND

The APETALA 2/ ethylene-responsive element binding factors (AP2/ERF), are thought to be associated with plant abiotic stress response, and involved in some plant hormone signaling pathways. Trichosanthes kirilowii is an important edible and medicinal crop, so far no research has been conducted on the TkAP2/ERF genes.

RESULT

In this study, a total of 135 TkERFs were identified, these genes were divided into 4 subfamilies and clustered into 13 groups. Moreover, 37 paralogous pairs were identified, with only two having Ka/Ks values greater than 1, proving that most TkERF genes underwent purifying selection during evolution. Co-expression networks constructed using transcriptome data at various flowering stages revealed that 50, 64, and 67 AP2/ERF genes correlated with members of the ethylene, gibberellin, and abscisic acid signaling pathways, respectively. When tissue cultured seedlings were treated with ETH, GA3 and ABA, 11, 12 and 17 genes were found to be up-regulated, respectively, suggesting that some members of the TkERF gene family may be involved in plant hormone signaling pathways. And under 4 ℃, PEG and NaCl treatment, 15, 20 and 19 genes were up-regulated, respectively, this suggested that these selected genes might be involved in plant abiotic stresses.

CONCLUSIONS

Overall, we identified 135 AP2/ERF family members, a comprehensive analysis of AP2/ERF gene expression patterns by RNA-seq and qRT-PCR showed that they played important roles in flower development and abiotic stress. This study provided a theoretical basis for the functional study of TkAP2/ERF genes and the genetic improvement of T. kirilowii.

Chromosome-level genome assembly of the Asian aspen Populus davidiana Dode

The genome of Populus davidiana, a keystone aspen species, has been sequenced to improve our understanding of the evolutionary and functional genomics of the Populus genus. The Hi-C scaffolding genome assembly resulted in a 408.1 Mb genome with 19 pseudochromosomes. The BUSCO assessment revealed that 98.3% of the genome matched the embryophytes dataset. A total of 31,862 protein-coding sequences were predicted, of which 31,619 were functionally annotated. The assembled genome was composed of 44.9% transposable elements. These findings provide new knowledge about the characteristics of the P. davidiana genome and will facilitate comparative genomics and evolutionary research on the genus Populus.

Welcome to the big leaves: Best practices for improving genome annotation in non-model plant genomes

Premise

Robust standards to evaluate quality and completeness are lacking in eukaryotic structural genome annotation, as genome annotation software is developed using model organisms and typically lacks benchmarking to comprehensively evaluate the quality and accuracy of the final predictions. The annotation of plant genomes is particularly challenging due to their large sizes, abundant transposable elements, and variable ploidies. This study investigates the impact of genome quality, complexity, sequence read input, and method on protein-coding gene predictions.

Methods

The impact of repeat masking, long-read and short-read inputs, and de novo and genome-guided protein evidence was examined in the context of the popular BRAKER and MAKER workflows for five plant genomes. The annotations were benchmarked for structural traits and sequence similarity.

Results

Benchmarks that reflect gene structures, reciprocal similarity search alignments, and mono-exonic/multi-exonic gene counts provide a more complete view of annotation accuracy. Transcripts derived from RNA-read alignments alone are not sufficient for genome annotation. Gene prediction workflows that combine evidence-based and ab initio approaches are recommended, and a combination of short and long reads can improve genome annotation. Adding protein evidence from de novo assemblies, genome-guided transcriptome assemblies, or full-length proteins from OrthoDB generates more putative false positives as implemented in the current workflows. Post-processing with functional and structural filters is highly recommended.

Discussion

While the annotation of non-model plant genomes remains complex, this study provides recommendations for inputs and methodological approaches. We discuss a set of best practices to generate an optimal plant genome annotation and present a more robust set of metrics to evaluate the resulting predictions.

Genome-Wide Analysis of Strictosidine Synthase-like Gene Family Revealed Their Response to Biotic/Abiotic Stress in Poplar

The strictosidine synthase-like (SSL) gene family is a small plant immune-regulated gene family that plays a critical role in plant resistance to biotic/abiotic stresses. To date, very little has been reported on the SSL gene in plants. In this study, a total of thirteen SSLs genes were identified from poplar, and these were classified into four subgroups based on multiple sequence alignment and phylogenetic tree analysis, and members of the same subgroup were found to have similar gene structures and motifs. The results of the collinearity analysis showed that poplar SSLs had more collinear genes in the woody plants Salix purpurea and Eucalyptus grandis. The promoter analysis revealed that the promoter region of PtrSSLs contains a large number of biotic/abiotic stress response elements. Subsequently, we examined the expression patterns of PtrSSLs following drought, salt, and leaf blight stress, using RT-qPCR to validate the response of PtrSSLs to biotic/abiotic stresses. In addition, the prediction of transcription factor (TF) regulatory networks identified several TFs, such as ATMYB46, ATMYB15, AGL20, STOP1, ATWRKY65, and so on, that may be induced in the expression of PtrSSLs in response to adversity stress. In conclusion, this study provides a solid basis for a functional analysis of the SSL gene family in response to biotic/abiotic stresses in poplar.

Oligo-FISH of Populus simonii Pachytene Chromosomes Improves Karyotyping and Genome Assembly

Poplar was one of the first woody species whose individual chromosomes could be identified using chromosome specific painting probes. Nevertheless, high-resolution karyotype construction remains a challenge. Here, we developed a karyotype based on the meiotic pachytene chromosome of Populus simonii which is a Chinese native species with many excellent traits. This karyotype was anchored by oligonucleotide (oligo)-based chromosome specific painting probes, a centromere-specific repeat (Ps34), ribosomal DNA, and telomeric DNA. We updated the known karyotype formula for P. simonii to 2n = 2x = 38 = 26m + 8st + 4t and the karyotype was 2C. The fluorescence in situ hybridization (FISH) results revealed some errors in the current P. simonii genome assembly. The 45S rDNA loci were located at the end of the short arm of chromosomes 8 and 14 by FISH. However, they were assembled on pseudochromosomes 8 and 15. In addition, the Ps34 loci were distributed in every centromere of the P. simonii chromosome in the FISH results, but they were only found to be present in pseudochromosomes 1, 3, 6, 10, 16, 17, 18, and 19. Our results reveal that pachytene chromosomes oligo-FISH is a powerful tool for constructing high-resolution karyotypes and improving the quality of genome assembly.