Chromosome-scale de novo genome assembly and annotation of three representative Casuarina species: C. equisetifolia, C. glauca, and C. cunninghamiana.

Australian pine (Casuarina spp.) is extensively planted in tropical and subtropical regions for wood production, shelterbelts, environmental protection, and ecological restoration due to their superior biological characteristics, such as rapid growth, wind and salt tolerance, and nitrogen fixation. To analyze the genomic diversity of Casuarina, we sequenced the genomes and constructed de novo genome assemblies of the three most widely planted Casuarina species: C. equisetifolia, C. glauca, and C. cunninghamiana. We generated chromosome-scale genome sequences using both Pacific Biosciences (PacBio) Sequel sequencing and chromosome conformation capture technology (Hi-C). The total genome sizes for C. equisetifolia, C. glauca, and C. cunninghamiana are 268 942 579 bp, 296 631 783 bp, and 293 483 606 bp, respectively, of which 25.91, 27.15, and 27.74% were annotated as repetitive sequences. We annotated 23 162, 24 673, and 24 674 protein-coding genes in C. equisetifolia, C. glauca, and C. cunninghamiana, respectively. We then collected branchlets from male and female individuals for whole-genome bisulfite sequencing (BS-seq) to explore the epigenetic regulation of sex determination in these three species. Transcriptome sequencing (RNA-seq) revealed differential expression of phytohormone-related genes between male and female plants. In summary, we generated three chromosome-level genome assemblies and comprehensive DNA methylation and transcriptome datasets from both male and female material for three Casuarina species, providing a basis for the comprehensive investigation of genomic diversity and functional gene discovery of Casuarina in the future.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa.

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

Whole-genome characterization of chronological age-associated changes in methylome and circular RNAs in moso bamboo (Phyllostachys edulis) from vegetative to floral growth.

In mammals, DNA methylation is associated with aging. However, age-related DNA methylation changes during phase transitions largely remain unstudied in plants. Moso bamboo (Phyllostachys edulis) requires a very long time to transition from the vegetative to the floral phase. To comprehensively investigate the association of DNA methylation with aging, we present here single-base-resolution DNA methylation profiles using both high-throughput bisulfite sequencing and single-molecule nanopore-based DNA sequencing, covering the long period of vegetative growth and transition to flowering in moso bamboo. We discovered that CHH methylation gradually accumulates from vegetative to reproductive growth in a time-dependent fashion. Differentially methylated regions, correlating with chronological aging, occurred preferentially at both transcription start sites and transcription termination sites. Genes with CG methylation changes showed an enrichment of Gene Ontology (GO) categories in 'vegetative to reproductive phase transition of meristem'. Combining methylation data with mRNA sequencing revealed that DNA methylation in promoters, introns and exons may have different roles in regulating gene expression. Finally, circular RNA (circRNA) sequencing revealed that the flanking introns of circRNAs are hypermethylated and enriched in long terminal repeat (LTR) retrotransposons. Together, the observations in this study provide insights into the dynamic DNA methylation and circRNA landscapes, correlating with chronological age, which paves the way to study further the impact of epigenetic factors on flowering in moso bamboo.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro).

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.

Genome-wide profiling of circular RNAs, alternative splicing, and R-loops in stem-differentiating xylem of Populus trichocarpa.

Circular RNAs (circRNAs) are a recently discovered type of non-coding RNA derived from pre-mRNAs. R-loops consist of a DNA:RNA hybrid and the associated single-stranded DNA. In Arabidopsis thaliana, circRNA:DNA R-loops regulate alternative splicing (AS) of SEPALLATA3 (SEP3). However, the occurrence and functions of circRNAs and R-loops in Populus trichocarpa are largely unexplored. Here, we performed circRNA-enriched sequencing in the stem-differentiating xylem (SDX) of P. trichocarpa and identified 2,742 distinct circRNAs, including circ-CESA4, circ-IRX7, and circ-GUX1, which are generated from genes involved in cellulose, and hemicellulose biosynthesis, respectively. To investigate the roles of circRNAs in modulating alternative splicing (AS), we detected 7,836 AS events using PacBio Iso-Seq and identified 634 circRNAs that overlapped with 699 AS events. Furthermore, using DNA:RNA hybrid immunoprecipitation followed by sequencing (DRIP-seq), we identified 8,932 R-loop peaks that overlapped with 181 circRNAs and 672 AS events. Notably, several SDX-related circRNAs overlapped with R-loop peaks, pointing to their possible roles in modulating AS in SDX. Indeed, overexpressing circ-IRX7 increased the levels of R-loop structures and decreased the frequency of intron retention in linear IRX7 transcripts. This study provides a valuable R-loop atlas resource and uncovers the interplay between circRNAs and AS in SDX of P. trichocarpa.

De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth.

Casuarina equisetifolia (C. equisetifolia), a conifer-like angiosperm with resistance to typhoon and stress tolerance, is mainly cultivated in the coastal areas of Australasia. C. equisetifolia, making it a valuable model to study secondary growth associated genes and stress-tolerance traits. However, the genome sequence is unavailable and therefore wood-associated growth rate and stress resistance at the molecular level is largely unexplored. We therefore constructed a high-quality draft genome sequence of C. equisetifolia by a combination of Illumina second-generation sequencing reads and Pacific Biosciences single-molecule real-time (SMRT) long reads to advance the investigation of this species. Here, we report the genome assembly, which contains approximately 300 megabases (Mb) and scaffold size of N50 is 1.06 Mb. Additionally, gene annotation, assisted by a combination of prediction and RNA-seq data, generated 29 827 annotated protein-coding genes and 1983 non-coding genes, respectively. Furthermore, we found that the total number of repetitive sequences account for one-third of the genome assembly. Here we also construct the genome-wide map of DNA modification, such as two novel forms N6 -adenine (6mA) and N4-methylcytosine (4mC) at the level of single-nucleotide resolution using single-molecule real-time (SMRT) sequencing. Interestingly, we found that 17% of 6mA modification genes and 15% of 4mC modification genes also included alternative splicing events. Finally, we investigated cellulose, hemicellulose, and lignin-related genes, which were associated with secondary growth and contained different DNA modifications. The high-quality genome sequence and annotation of C. equisetifolia in this study provide a valuable resource to strengthen our understanding of the diverse traits of trees.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species.

MOTIVATION: MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species. RESULTS: In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms. AVAILABILITY AND IMPLEMENTATION: The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-Wide Profiling of Circular RNAs in the Rapidly Growing Shoots of Moso Bamboo (Phyllostachys edulis).

Circular RNAs, including circular exonic RNAs (circRNA), circular intronic RNAs (ciRNA) and exon-intron circRNAs (EIciRNAs), are a new type of noncoding RNAs. Growing shoots of moso bamboo (Phyllostachys edulis) represent an excellent model of fast growth and their circular RNAs have not been studied yet. To understand the potential regulation of circular RNAs, we systematically characterized circular RNAs from eight different developmental stages of rapidly growing shoots. Here, we identified 895 circular RNAs including a subset of mutually inclusive circRNA. These circular RNAs were generated from 759 corresponding parental coding genes involved in cellulose, hemicellulose and lignin biosynthetic process. Gene co-expression analysis revealed that hub genes, such as DEFECTIVE IN RNA-DIRECTED DNA METHYLATION 1 (DRD1), MAINTENANCE OF METHYLATION (MOM), dicer-like 3 (DCL3) and ARGONAUTE 1 (AGO1), were significantly enriched giving rise to circular RNAs. The expression level of these circular RNAs presented correlation with its linear counterpart according to transcriptome sequencing. Further protoplast transformation experiments indicated that overexpressing circ-bHLH93 generating from transcription factor decreased its linear transcript. Finally, the expression profiles suggested that circular RNAs may have interplay with miRNAs to regulate their cognate linear mRNAs, which was further supported by overexpressing miRNA156 decreasing the transcript of circ-TRF-1 and linear transcripts of TRF-1. Taken together, the overall profile of circular RNAs provided new insight into an unexplored category of long noncoding RNA regulation in moso bamboo.

PRAPI: post-transcriptional regulation analysis pipeline for Iso-Seq.

Summary: The single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) based on Pacific Bioscience (PacBio) platform has received increasing attention for its ability to explore full-length isoforms. Thus, comprehensive tools for Iso-Seq bioinformatics analysis are extremely useful. Here, we present a one-stop solution for Iso-Seq analysis, called PRAPI to analyze alternative transcription initiation (ATI), alternative splicing (AS), alternative cleavage and polyadenylation (APA), natural antisense transcripts (NAT), and circular RNAs (circRNAs) comprehensively. PRAPI is capable of combining Iso-Seq full-length isoforms with short read data, such as RNA-Seq or polyadenylation site sequencing (PAS-seq) for differential expression analysis of NAT, AS, APA and circRNAs. Furthermore, PRAPI can annotate new genes and correct mis-annotated genes when gene annotation is available. Finally, PRAPI generates high-quality vector graphics to visualize and highlight the Iso-Seq results. Availability and implementation: The Dockerfile of PRAPI is available at http://www.bioinfor.org/tool/PRAPI. Contact: lfgu@fafu.edu.cn.

Comprehensive profiling of rhizome-associated alternative splicing and alternative polyadenylation in moso bamboo (Phyllostachys edulis).

Moso bamboo (Phyllostachys edulis) represents one of the fastest-spreading plants in the world, due in part to its well-developed rhizome system. However, the post-transcriptional mechanism for the development of the rhizome system in bamboo has not been comprehensively studied. We therefore used a combination of single-molecule long-read sequencing technology and polyadenylation site sequencing (PAS-seq) to re-annotate the bamboo genome, and identify genome-wide alternative splicing (AS) and alternative polyadenylation (APA) in the rhizome system. In total, 145 522 mapped full-length non-chimeric (FLNC) reads were analyzed, resulting in the correction of 2241 mis-annotated genes and the identification of 8091 previously unannotated loci. Notably, more than 42 280 distinct splicing isoforms were derived from 128 667 intron-containing full-length FLNC reads, including a large number of AS events associated with rhizome systems. In addition, we characterized 25 069 polyadenylation sites from 11 450 genes, 6311 of which have APA sites. Further analysis of intronic polyadenylation revealed that LTR/Gypsy and LTR/Copia were two major transposable elements within the intronic polyadenylation region. Furthermore, this study provided a quantitative atlas of poly(A) usage. Several hundred differential poly(A) sites in the rhizome-root system were identified. Taken together, these results suggest that post-transcriptional regulation may potentially have a vital role in the underground rhizome-root system.

Transcriptome characterization of moso bamboo (Phyllostachys edulis) seedlings in response to exogenous gibberellin applications.

BACKGROUND: Moso bamboo (Phyllostachys edulis) is a well-known bamboo species of high economic value in the textile industry due to its rapid growth. Phytohormones, which are master regulators of growth and development, serve as important endogenous signals. However, the mechanisms through which phytohormones regulate growth in moso bamboo remain unknown to date. RESULTS: Here, we reported that exogenous gibberellins (GA) applications resulted in a significantly increased internode length and lignin condensation. Transcriptome sequencing revealed that photosynthesis-related genes were enriched in the GA-repressed gene class, which was consistent with the decrease in leaf chlorophyll concentrations and the lower rate of photosynthesis following GA treatment. Exogenous GA applications on seedlings are relatively easy to perform, thus we used 4-week-old whole seedlings of bamboo for GA- treatment followed by high throughput sequencing. In this study, we identified 932 cis-nature antisense transcripts (cis-NATs), and 22,196 alternative splicing (AS) events in total. Among them, 42 cis-nature antisense transcripts (cis-NATs) and 442 AS events were differentially expressed upon exposure to exogenous GA3, suggesting that post-transcriptional regulation might be also involved in the GA3 response. Targets of differential expression of cis-NATs included genes involved in hormone receptor, photosynthesis and cell wall biogenesis. For example, LAC4 and its corresponding cis-NATs were GA3-induced, and may be involved in the accumulation of lignin, thus affecting cell wall composition. CONCLUSIONS: This study provides novel insights illustrating how GA alters post-transcriptional regulation and will shed light on the underlying mechanism of growth modulated by GA in moso bamboo.

Genome-wide analysis and transcriptomic profiling of the auxin biosynthesis, transport and signaling family genes in moso bamboo (Phyllostachys heterocycla).

BACKGROUND: Auxin is essential for plant growth and development. Although substantial progress has been made in understanding auxin pathways in model plants such as Arabidopsis and rice, little is known in moso bamboo which is famous for its fast growth resulting from the rapid cell elongation and division. RESULTS: Here we showed that exogenous auxin has strong effects on crown and primary roots. Genes involved in auxin action, including 13 YUCCA (YUC) genes involved in auxin synthesis, 14 PIN-FORMED/PIN-like (PIN/PILS) and 7 AUXIN1/LIKE-AUX1 (AUX1/LAX) members involved in auxin transport, 10 auxin receptors (AFB) involved in auxin perception, 43 auxin/indole-3-aceticacid (AUX/IAA) genes, and 41 auxin response factors (ARF) involved in auxin signaling were identified through genome-wide analysis. Phylogenetic analysis of these genes from Arabidopsis, Oryza sativa and bamboo revealed that auxin biosynthesis, transport, and signaling pathways are conserved in these species. A comprehensive study of auxin-responsive genes using RNA sequencing technology was performed, and the results also supported that moso bamboo shared a conserved regulatory mechanism for the expression of auxin pathway genes; meanwhile it harbors its own specific properties. CONCLUSIONS: In summary, we generated an overview of the auxin pathway in bamboo, which provides information for uncovering the precise roles of auxin pathway in this important species in the future.

Telomere-to-telomere genome assembly of Oldenlandia diffusa.

We report the complete telomere-to-telomere genome assembly of Oldenlandia diffusa which renowned in traditional Chinese medicine, comprising 16 chromosomes and spanning 499.7 Mb. The assembly showcases 28 telomeres and minimal gaps, with a total of only five. Repeat sequences constitute 46.41% of the genome, and 49,701 potential protein-coding genes have been predicted. Compared with O. corymbosa, O. diffusa exhibits chromosome duplication and fusion events, diverging 20.34 million years ago. Additionally, a total of 11 clusters of terpene synthase have been identified. The comprehensive genome sequence, gene catalog, and terpene synthase clusters of O. diffusa detailed in this study will significantly contribute to advancing research in this species' genetic, genomic, and pharmacological aspects.

Multi-omics of Circular RNAs and Their Responses to Hormones in Moso Bamboo (Phyllostachys edulis).

Circular RNAs (circRNAs) are endogenous non-coding RNAs with covalently closed structures, which have important functions in plants. However, their biogenesis, degradation, and function upon treatment with gibberellins (GAs) and auxins (1-naphthaleneacetic acid, NAA) remain unknown. Here, we systematically identified and characterized the expression patterns, evolutionary conservation, genomic features, and internal structures of circRNAs using RNase R-treated libraries from moso bamboo (Phyllostachys edulis) seedlings. Moreover, we investigated the biogenesis of circRNAs dependent on both cis- and trans-regulation. We explored the function of circRNAs, including their roles in regulating microRNA (miRNA)-related genes and modulating the alternative splicing of their linear counterparts. Importantly, we developed a customized degradome sequencing approach to detect miRNA-mediated cleavage of circRNAs. Finally, we presented a comprehensive view of the participation of circRNAs in the regulation of hormone metabolism upon treatment of bamboo seedlings with GA and NAA. Collectively, our study provides insights into the biogenesis, function, and miRNA-mediated degradation of circRNAs in moso bamboo.

The pattern of DNA methylation alteration, and its association with the expression changes of non-coding RNAs and mRNAs in Moso bamboo under abiotic stress.

Epigenetic changes play an important role in plant growth and development and in stress response. However, DNA methylation pattern and its relationship with the expression changes of non-coding RNAs and mRNAs of Moso bamboo in response to abiotic stress is still largely unknown. In this work, we used whole-genome bisulfite sequencing in combination with whole-transcriptome sequencing to analyze the DNA methylation and transcription patterns of mRNAs and non-coding RNAs in Moso bamboo under abiotic stresses such as cold, heat, ultraviolet (UV) and salinity. We found that CHH methylation in the promoter region was positively correlated with gene expression, while CHG and CHH methylations in the gene body regions were negatively associated with gene expression. Moreover, CG and CHG methylations in the promoter regions were negatively correlated with the transcript abundance of long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs). Similarly, the methylation levels of three contexts in the genic regions were negatively correlated with the transcript abundance of lncRNAs and miRNAs but positively correlated with that of circRNAs. In addition, we suggested that the reduction of 21-nt and 24-nt small interfering RNA (siRNA) expression tended to increase methylation levels in the genic regions. We found that stress-responsive genes such as CRPK1, HSFB2A and CIPK were differentially methylated and expressed. Our results also proposed that DNA methylation may regulate the expression of the transcription factors (TFs) and plant hormone signalling genes such as IAA9, MYC2 and ERF110 in response to abiotic stress. This study firstly reports the abiotic stress-responsive DNA methylation pattern and its involvement of expression of coding RNAs and non-coding RNAs in Moso bamboo. The results expand the knowledge of epigenetic mechanisms in Moso bamboo under abiotic stress and support in-depth deciphering of the function of specific non-coding RNAs in future studies.

Comprehensive Transcriptome Analysis of Stem-Differentiating Xylem Upon Compression Stress in Cunninghamia Lanceolata.

Compression wood (CW) in gymnosperm brings great difficulties to wood industry using wood as raw materials since CW presents special wood structure and have different physical and chemical properties from those of normal wood (NW). Chinese fir (Cunninghamia lanceolata) is widely distributed in China. However, global transcriptome profiling of coding and long non-coding RNA in response to compression stress has not been reported in the gymnosperm species. In this study, we revealed that CW in Chinese fir exhibited distinct morphology and cytology properties compared with those of NW, including high lignin content, thick and round tracheid cells. Furthermore, we combined both PacBio long-read SMRT sequencing (Iso-Seq) and Illumina short-read RNA-Seq to reveal the transcriptome in stem-differentiating xylem (SDX) under different time points (2, 26, and 74 h) upon compression stress in NW, CW, and OW (opposite wood), respectively. Iso-Seq was successfully assembled into 41,253 de-novo full-length transcriptome reference (average length 2,245 bp). Moreover, there were striking differences in expression upon compression stress, which were involved 13 and 7 key enzyme genes in the lignin and cellulose synthesis, respectively. Especially, we revealed 11 secondary growth-related transcription factors show differential expression under compression stress, which was further validated by qRT-PCR. Finally, the correlation between 6,533 differentially expressed coding genes and 372 differentially expressed long non-coding RNAs (lncRNAs) indicates that these lncRNAs may affect cell wall biogenesis and xyloglucan metabolism. In conclusion, our results provided comprehensive cytology properties and full-length transcriptome profiling of wood species upon compression stress. Especially we explored candidate genes, including both coding and long non-coding genes, and provided a theoretical basis for further research on the formation mechanism of CW in gymnosperm Chinese fir.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

Quantitative profiling of N6-methyladenosine at single-base resolution in stem-differentiating xylem of Populus trichocarpa using Nanopore direct RNA sequencing.

There are no comprehensive methods to identify N6-methyladenosine (m6A) at single-base resolution for every single transcript, which is necessary for the estimation of m6A abundance. We develop a new pipeline called Nanom6A for the identification and quantification of m6A modification at single-base resolution using Nanopore direct RNA sequencing based on an XGBoost model. We validate our method using methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and m6A-sensitive RNA-endoribonuclease-facilitated sequencing (m6A-REF-seq), confirming high accuracy. Using this method, we provide a transcriptome-wide quantification of m6A modification in stem-differentiating xylem and reveal that different alternative polyadenylation (APA) usage shows a different ratio of m6A.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses.

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.