N6-methyladenosine mRNA methylation positively regulated the response of poplar to salt stress.

As the most abundant form of methylation modification in messenger RNA (mRNA), the distribution of N6-methyladenosine (m6A) has been preliminarily revealed in herbaceous plants under salt stress, but its function and mechanism in woody plants were still unknown. Here, we showed that global m6A levels increased during poplar response to salt stress. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) revealed that m6A significantly enriched in the coding sequence region and 3'-untranslated regions in poplar, by recognising the conserved motifs, AGACU, GGACA and UGUAG. A large number of differential m6A transcripts have been identified, and some have been proved involving in salt response and plant growth and development. Further combined analysis of MeRIP-seq and RNA-seq revealed that the m6A hypermethylated and enrich in the CDS region preferred to positively regulate expression abundance. Writer inhibitor, 3-deazaneplanocin A treatment increased the sensitivity of poplar to salt stress by reducing mRNA stability to regulate the expression of salt-responsive transcripts PagMYB48, PagGT2, PagNAC2, PagGPX8 and PagARF2. Furthermore, we verified that the methyltransferase PagFIP37 plays a positively role in the response of poplar to salt stress, overexpressed lines have stronger salt tolerance, while RNAi lines were more sensitive to salt, which relied on regulating mRNA stability in an m6A manner of salt-responsive transcripts PagMYB48, PagGT2, PagNAC2, PagGPX8 and PagARF2. Collectively, these results revealed the regulatory role of m6A methylation in poplar response to salt stress, and revealed the importance and mechanism of m6A methylation in the response of woody plants to salt stress for the first time.

Transcriptomic and Phenotypic Analyses Reveal the Molecular Mechanism of Dwarfing in Tetraploid Robinia pseudoacacia L.

Polyploid breeding techniques aid in the cultivation of new forestry cultivars, thus expanding the suite of strategies for the improvement of arboreal traits and innovation within the field of forestry. Compared to diploid Robinia pseudoacacia L. (black locust) 'D26-5①' (2×), its dwarfed homologous tetraploid 'D26-5②' (4×) variety has better application prospects in garden vegetation guardrails and urban landscape. However, the molecular mechanism of the generation and growth of this dwarf variety is still unclear. Here, plant growth and development as well as histological differences between the diploid and its autotetraploid were investigated. Levels of endogenous hormones at three different developmental stages (20, 40, and 70 days) of 2× and homologous 4× tissue culture plantlets were assessed, and it was found that the brassinosteroid (BR) contents of the former were significantly higher than the latter. Transcriptome sequencing data analysis of 2× and homologous 4× showed that differentially expressed genes (DEGs) were significantly enriched in plant hormone synthesis and signal transduction, sugar and starch metabolism, and the plant circadian rhythm pathway, which are closely related to plant growth and development. Therefore, these biological pathways may be important regulatory pathways leading to dwarfism and slow growth in tetraploids. Additionally, utilizing weighted gene coexpression network analysis (WGCNA), we identified three crucial differentially expressed genes (DEGs)-PRR5, CYP450, and SPA1-that potentially underlie the observed ploidy variation. This study provides a new reference for the molecular mechanism of dwarfism in dwarfed autotetraploid black locusts. Collectively, our results of metabolite analysis and comparative transcriptomics confirm that plant hormone signaling and the circadian rhythm pathway result in dwarfism in black locusts.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation.

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

SubPhaser: a robust allopolyploid subgenome phasing method based on subgenome-specific k-mers.

With advanced sequencing technology, dozens of complex polyploid plant genomes have been characterized. However, for many polyploid species, their diploid ancestors are unknown or extinct, making it impossible to unravel the subgenomes and genome evolution directly. We developed a novel subgenome-phasing algorithm, SubPhaser, specifically designed for a neoallopolyploid or a homoploid hybrid. SubPhaser first searches for the subgenome-specific sequence (k-mer), then assigns homoeologous chromosomes into subgenomes, and further provides tools to annotate and investigate specific sequences. SubPhaser works well on neoallopolyploids and homoploid hybrids containing subgenome-specific sequences like wheat, but fails on autopolyploids lacking subgenome-specific sequences like alfalfa, indicating that SubPhaser can phase neoallopolyploid/homoploid hybrids with high accuracy, sensitivity and performance. This highly accurate, highly sensitive, ancestral data free chromosome phasing algorithm, SubPhaser, offers significant application value for subgenome phasing in neoallopolyploids and homoploid hybrids, and for the subsequent exploration of genome evolution and related genetic/epigenetic mechanisms.

UV-B and UV-C radiation trigger both common and distinctive signal perceptions and transmissions in Pinus tabuliformis Carr.

In plants, ultraviolet (UV)-light is an important driver for growth and natural distribution, and is also a valuable tool for manipulating productivity as well as biotic interactions. Understanding of plant responses to different UV radiation is sparse, especially from a systems biology perspective and particularly for conifers. Here, we evaluated the physiological and transcriptomic responses to the short-term application of high-irradiance UV-B and UV-C waves on Pinus tabuliformis Carr., a major conifer in Northern China. By undertaking time-ordered gene coexpression network analyses and network comparisons incorporating physiological traits and gene expression variation, we uncovered communalities but also differences in P. tabuliformis responses to UV-B and UV-C. Both types of spectral bands caused a significant inhibition of photosynthesis, and conversely, the improvement of antioxidant capacity, flavonoid production and signaling pathways related to stress resistance, indicating a clear switch from predominantly primary metabolism to enhanced defensive metabolism in pine. We isolated distinct subnetworks for photoreceptor-mediated signal transduction, maximum quantum efficiency of photosystem II (Fv/Fm) regulation and flavonoid biosynthesis in response to UV-B and UV-C radiation. From these subnetworks, we further identified phototropins as potentially important elements in both UV-B and UV-C signaling and, for the first time, suggesting peptide hormones to be involved in promoting flavonoid biosynthesis against UV-B, while these hormones seem not to be implicated in the defense against UV-C exposure. The present study employed an effective strategy for disentangling the complex physiological and genetic regulatory mechanisms in a nonmodel plant species, and thus, provides a suitable reference for future functional evaluations and artificial UV-light mediated growing strategies in plant production.

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

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

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus.

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome.

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

Transcriptome-wide analysis of introgression-resistant regions reveals genetic divergence genes under positive selection in Populus trichocarpa.

Comparing gene expression patterns and genetic polymorphisms between populations is of central importance for understanding the origin and maintenance of biodiversity. Based on population-specific gene expression levels and allele frequency differences, we sought to identify population divergence (PD) genes across the introgression-resistant genomic regions of Populus trichocarpa. Genes containing highly diverged loci [i.e., genetic divergence (GD)] or showing expression divergence (ED) between populations were widely distributed in the genome and substantially enriched in functional categories related to stress responses, disease resistance, timing of flowering, cell cycle regulation, plant growth, and development. Nine genomic regions showing evidence of strong positive selection were overlapped with GD genes, which had significant differences between Oregon (a southernmost peripheral deme) and the other demes. However, we did not find evidence that genes under positive selection show an enrichment for ED. PD genes and genes under selection pertained to the same gene classes, such as SERINE/CYSTEINE PROTEASE, ABC TRANSPORTER, GLYCOSYLTRANSFERASE and other transferases. Our analysis also revealed that GD genes were polymorphic within the species (41.9 ± 3.66 biallelic variants per gene), as previously reported in herbaceous plants. By contrast, ED genes contained less genetic variants (10.73 ± 1.14) and were likely highly expressed. In addition, we found that trans- rather than cis-acting variants considerably contribute to the evolution of >90% PD genes. Overall, this study elucidates that cohorts of PD genes agree with the general attributes of known speciation genes and GD genes will provide substrates for positive selection to operate on.

Association genetics, geography and ecophysiology link stomatal patterning in Populus trichocarpa with carbon gain and disease resistance trade-offs.

Stomata are essential for diffusive entry of gases to support photosynthesis, but may also expose internal leaf tissues to pathogens. To uncover trade-offs in range-wide adaptation relating to stomata, we investigated the underlying genetics of stomatal traits and linked variability in these traits with geoclimate, ecophysiology, condensed foliar tannins and pathogen susceptibility in black cottonwood (Populus trichocarpa). Upper (adaxial) and lower (abaxial) leaf stomatal traits were measured from 454 accessions collected throughout much of the species range. We calculated broad-sense heritability (H(2) ) of stomatal traits and, using SNP data from a 34K Populus SNP array, performed a genome-wide association studies (GWAS) to uncover genes underlying stomatal trait variation. H(2) values for stomatal traits were moderate (average H(2) = 0.33). GWAS identified genes associated primarily with adaxial stomata, including polarity genes (PHABULOSA), stomatal development genes (BRASSINOSTEROID-INSENSITIVE 2) and disease/wound-response genes (GLUTAMATE-CYSTEINE LIGASE). Stomatal traits correlated with latitude, gas exchange, condensed tannins and leaf rust (Melampsora) infection. Latitudinal trends of greater adaxial stomata numbers and guard cell pore size corresponded with higher stomatal conductance (gs ) and photosynthesis (Amax ), faster shoot elongation, lower foliar tannins and greater Melampsora susceptibility. This suggests an evolutionary trade-off related to differing selection pressures across the species range. In northern environments, more adaxial stomata and larger pore sizes reflect selection for rapid carbon gain and growth. By contrast, southern genotypes have fewer adaxial stomata, smaller pore sizes and higher levels of condensed tannins, possibly linked to greater pressure from natural leaf pathogens, which are less significant in northern ecosystems.