TrancriptomeReconstructoR: data-driven annotation of complex transcriptomes

Global identification and regulatory network analysis reveal the significant roles of lncRNAs during anther and pollen development in Arabidopsis

KEY MESSAGE

A high-throughput sequencing identified 1283 lncRNAs in anthers at different stages in Arabidopsis and their relationship with protein-coding genes and miRNAs during anther and pollen development were analyzed. Long non-coding RNAs (lncRNAs) are important regulatory molecules involved in various biological processes. However, their roles in male reproductive development and interactions with miRNAs remained elusive. In this study, a high-throughput sequencing of anthers at different developmental stages in Arabidopsis identified 1283 lncRNAs including 524 differentially expressed lncRNAs (DELs). Most of these DELs exhibited positive correlations with the expression patterns of adjacent protein-coding genes. Weighted gene co-expression network analysis (WGCNA) revealed that protein-coding genes targeted by DELs in four modules related to the tetrad stage were associated with functions such as pollen wall formation, pollen germination, or pollen tube growth, respectively. Furthermore, five, 10, and 11 lncRNAs were predicted as miRNAs' endogenous target mimics (eTMs), precursors, and natural antisense transcripts of pri-miRNA, respectively. Remarkably, the lncRNA, host gene of ath-miR167a (ath-miR167aHG), predicted as the precursor of miR167a, was selected for function validation. Its overexpression resulted in the up-regulation of miR167a and the subsequent down-regulation of miR167a's target genes ARF6 and ARF8, demonstrating a functional interaction between ath-miR167aHG and miR167a. The transgenic plants showed delayed flowering, shorter filaments, abnormal anther dehiscence, and undeveloped siliques ultimately, suggesting a role of ath-miR167aHG in male reproductive development. Collectively, our research shed new light on the functions of lncRNAs in male reproductive development and uncovered the unique interactions between lncRNAs and miRNAs.

Efficient small fragment sequencing of human, cattle, and bison miRNA, small RNA, or csRNA-seq libraries using AVITI

BACKGROUND

Next-Generation Sequencing (NGS) catalyzed breakthroughs across various scientific domains. Illumina's sequencing by synthesis method has long been central to NGS, but new sequencing methods like Element Biosciences' AVITI technology are emerging. AVITI is reported to offer improved signal-to-noise ratios and cost reductions. However, its reliance on rolling circle amplification, which can be affected by polymer size, raises questions about its effectiveness in sequencing small RNAs (sRNAs) such as microRNAs (miRNAs), small nucleolar RNAs (snoRNAs), and many others. These sRNAs are crucial regulators of gene expression and involved in various biological processes. Additionally, capturing capped small RNAs (csRNA-seq) is a powerful method for mapping active or "nascent" RNA polymerase II transcription initiation in tissues and clinical samples.

RESULTS

Here, we report a new protocol for seamlessly sequencing short fragments on the AVITI and demonstrate that AVITI and Illumina sequencing technologies equivalently capture human, cattle (Bos taurus), and bison (Bison bison) sRNA or csRNA sequencing libraries, increasing confidence in both sequencing approaches. Additionally, analysis of generated nascent transcription start site (TSS) data for cattle and bison revealed inaccuracies in their current genome annotations, underscoring the potential and necessity to translate small and nascent RNA sequencing methodologies to livestock.

CONCLUSIONS

Our accelerated and optimized protocol bridges the advantages of AVITI sequencing with critical methods that rely on sequencing short fragments. This advance bolsters the utility of AVITI technology alongside traditional Illumina platforms, offering new opportunities for NGS applications.

HPC-T-Annotator: an HPC tool for de novo transcriptome assembly annotation

BACKGROUND

The availability of transcriptomic data for species without a reference genome enables the construction of de novo transcriptome assemblies as alternative reference resources from RNA-Seq data. A transcriptome provides direct information about a species' protein-coding genes under specific experimental conditions. The de novo assembly process produces a unigenes file in FASTA format, subsequently targeted for the annotation. Homology-based annotation, a method to infer the function of sequences by estimating similarity with other sequences in a reference database, is a computationally demanding procedure.

RESULTS

To mitigate the computational burden, we introduce HPC-T-Annotator, a tool for de novo transcriptome homology annotation on high performance computing (HPC) infrastructures, designed for straightforward configuration via a Web interface. Once the configuration data are given, the entire parallel computing software for annotation is automatically generated and can be launched on a supercomputer using a simple command line. The output data can then be easily viewed using post-processing utilities in the form of Python notebooks integrated in the proposed software.

CONCLUSIONS

HPC-T-Annotator expedites homology-based annotation in de novo transcriptome assemblies. Its efficient parallelization strategy on HPC infrastructures significantly reduces computational load and execution times, enabling large-scale transcriptome analysis and comparison projects, while its intuitive graphical interface extends accessibility to users without IT skills.

Trans2express - de novo transcriptome assembly pipeline optimized for gene expression analysis

BACKGROUND

As genomes of many eukaryotic species, especially plants, are large and complex, their de novo sequencing and assembly is still a difficult task despite progress in sequencing technologies. An alternative to genome assembly is the assembly of transcriptome, the set of RNA products of the expressed genes. While a bunch of de novo transcriptome assemblers exists, the challenges of transcriptomes (the existence of isoforms, the uneven expression levels across genes) complicates the generation of high-quality assemblies suitable for downstream analyses.

RESULTS

We developed Trans2express - a web-based tool and a pipeline of de novo hybrid transcriptome assembly and postprocessing based on rnaSPAdes with a set of subsequent filtrations. The pipeline was tested on Arabidopsis thaliana cDNA sequencing data obtained using Illumina and Oxford Nanopore Technologies platforms and three non-model plant species. The comparison of structural characteristics of the transcriptome assembly with reference Arabidopsis genome revealed the high quality of assembled transcriptome with 86.1% of Arabidopsis expressed genes assembled as a single contig. We tested the applicability of the transcriptome assembly for gene expression analysis. For both Arabidopsis and non-model species the results showed high congruence of gene expression levels and sets of differentially expressed genes between analyses based on genome and based on the transcriptome assembly.

CONCLUSIONS

We present Trans2express - a protocol for de novo hybrid transcriptome assembly aimed at recovering of a single transcript per gene. We expect this protocol to promote the characterization of transcriptomes and gene expression analysis in non-model plants and web-based tool to be of use to a wide range of plant biologists.

Efficient small fragment sequencing of human, cow, and bison miRNA, small RNA or csRNA-seq libraries using AVITI

Next-Generation Sequencing (NGS) catalyzed breakthroughs across various scientific domains. Illumina's sequencing by synthesis method has long been essential for NGS but emerging technologies like Element Biosciences' sequencing by avidity (AVITI) represent a novel approach. It has been reported that AVITI offers improved signal-to-noise ratios and cost reductions. However, the method relies on rolling circle amplification which can be impacted by polymer size, raising questions about its efficacy sequencing small RNAs (sRNA) molecules including microRNAs (miRNAs), piwi-interacting RNAs (piRNAs), and others that are crucial regulators of gene expression and involved in various biological processes. In addition, capturing capped small RNAs (csRNA-seq) has emerged as a powerful method to map active or "nascent" RNA polymerase II transcription initiation in tissues and clinical samples. Here, we report a new protocol for seamlessly sequencing short DNA fragments on the AVITI and demonstrate that AVITI and Illumina sequencing technologies equivalently capture human, cattle (Bos taurus) and the bison (Bison bison) sRNA or csRNA sequencing libraries, augmenting the confidence in both approaches. Additionally, analysis of generated nascent transcription start sites (TSSs) data for cattle and bison revealed inaccuracies in their current genome annotations and highlighted the possibility and need to translate small RNA sequencing methodologies to livestock. Our accelerated and optimized protocol therefore bridges the advantages of AVITI sequencing and critical methods that rely on sequencing short DNA fragments.

Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics

High-throughput detection of nascent RNA is critical for studies of transcription and much more challenging than that of mRNA. Recently, several massively parallel nascent RNA sequencing methods were established in eukaryotic cells. Here, we systematically compared 3 classes of methods on the same pure or crude nuclei preparations: GRO-seq for sequence nuclear run-on RNAs, pNET-seq for sequence RNA polymerase II-associated RNAs, and CB RNA-seq for sequence chromatin-bound (CB) RNAs in Arabidopsis (Arabidopsis thaliana). To improve the resolution of CB RNAs, 3'CB RNA-seq was established to sequence the 3' ends of CB RNAs. In addition, we modified pNET-seq to establish the Chromatin Native Elongation Transcript sequencing (ChrNET) method using chromatin as the starting material for RNA immunoprecipitation. Reproducibility, sensitivity and accuracy in detecting nascent transcripts, experimental procedures, and costs were analyzed, which revealed the strengths and weaknesses of each method. We found that pNET and GRO methods best detected active RNA polymerase II. CB RNA-seq is a simple and cost-effective alternative for nascent RNA studies, due to its high correlation with pNET-seq and GRO-seq. Compared with pNET, ChrNET has higher specificity for nascent RNA capture and lower sequencing cost. 3'CB is sensitive to transcription-coupled splicing. Using these methods, we identified 1,404 unknown transcripts, 4,482 unannotated splicing events, and 60 potential recursive splicing events. This comprehensive comparison of different nascent/chromatin RNA sequencing methods highlights the strengths of each method and serves as a guide for researchers aiming to select a method that best meets their study goals.

NICOTIANAMINE SYNTHASE activity affects nucleolar iron accumulation and impacts rDNA silencing and RNA methylation in Arabidopsis

In plant cells, a large pool of iron (Fe) is contained in the nucleolus, as well as in chloroplasts and mitochondria. A central determinant for intracellular distribution of Fe is nicotianamine (NA) generated by NICOTIANAMINE SYNTHASE (NAS). Here, we used Arabidopsis thaliana plants with disrupted NAS genes to study the accumulation of nucleolar iron and understand its role in nucleolar functions and more specifically in rRNA gene expression. We found that nas124 triple mutant plants, which contained lower quantities of the iron ligand NA, also contained less iron in the nucleolus. This was concurrent with the expression of normally silenced rRNA genes from nucleolar organizer regions 2 (NOR2). Notably, in nas234 triple mutant plants, which also contained lower quantities of NA, nucleolar iron and rDNA expression were not affected. In contrast, in both nas124 and nas234, specific RNA modifications were differentially regulated in a genotype dependent manner. Taken together, our results highlight the impact of specific NAS activities in RNA gene expression. We discuss the interplay between NA and nucleolar iron with rDNA functional organization and RNA methylation.

LncRNA FLAIL affects alternative splicing and represses flowering in Arabidopsis

How the noncoding genome affects cellular functions is a key biological question. A particular challenge is to distinguish the effects of noncoding DNA elements from long noncoding RNAs (lncRNAs) that coincide at the same loci. Here, we identified the flowering-associated intergenic lncRNA (FLAIL) in Arabidopsis through early flowering flail mutants. Expression of FLAIL RNA from a different chromosomal location in combination with strand-specific RNA knockdown characterized FLAIL as a trans-acting RNA molecule. FLAIL directly binds to differentially expressed target genes that control flowering via RNA-DNA interactions through conserved sequence motifs. FLAIL interacts with protein and RNA components of the spliceosome to affect target mRNA expression through co-transcriptional alternative splicing (AS) and linked chromatin regulation. In the absence of FLAIL, splicing defects at the direct FLAIL target flowering gene LACCASE 8 (LAC8) correlated with reduced mRNA expression. Double mutant analyses support a model where FLAIL-mediated splicing of LAC8 promotes its mRNA expression and represses flowering. Our study suggests lncRNAs as accessory components of the spliceosome that regulate AS and gene expression to impact organismal development.

Functional variation in the non-coding genome: molecular implications for food security

The growing world population, in combination with the anticipated effects of climate change, is pressuring food security. Plants display an impressive arsenal of cellular mechanisms conferring resilience to adverse environmental conditions, and humans rely on these mechanisms for stable food production. The elucidation of the molecular basis of the mechanisms used by plants to achieve resilience promises knowledge-based approaches to enhance food security. DNA sequence polymorphisms can reveal genomic regions that are linked to beneficial traits of plants. However, our ability to interpret how a given DNA sequence polymorphism confers a fitness advantage at the molecular level often remains poor. A key factor is that these polymorphisms largely localize to the enigmatic non-coding genome. Here, we review the functional impact of sequence variations in the non-coding genome on plant biology in the context of crop breeding and agricultural traits. We focus on examples of non-coding with particularly convincing functional support. Our survey combines findings that are consistent with the view that the non-coding genome contributes to cellular mechanisms assisting many plant traits. Understanding how DNA sequence polymorphisms in the non-coding genome shape plant traits at the molecular level offers a largely unexplored reservoir of solutions to address future challenges in plant growth and resilience.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

Unique features of transcription termination and initiation at closely spaced tandem human genes

The synthesis of RNA polymerase II (Pol2) products, which include messenger RNAs or long noncoding RNAs, culminates in transcription termination. How the transcriptional termination of a gene impacts the activity of promoters found immediately downstream of it, and which can be subject to potential transcriptional interference, remains largely unknown. We examined in an unbiased manner the features of the intergenic regions between pairs of 'tandem genes'-closely spaced (< 2 kb) human genes found on the same strand. Intergenic regions separating tandem genes are enriched with guanines and are characterized by binding of several proteins, including AGO1 and AGO2 of the RNA interference pathway. Additionally, we found that Pol2 is particularly enriched in this region, and it is lost upon perturbations affecting splicing or transcriptional elongation. Perturbations of genes involved in Pol2 pausing and R loop biology preferentially affect expression of downstream genes in tandem gene pairs. Overall, we find that features associated with Pol2 pausing and accumulation rather than those associated with avoidance of transcriptional interference are the predominant driving force shaping short tandem intergenic regions.

The Hda1 histone deacetylase limits divergent non-coding transcription and restricts transcription initiation frequency

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA and a divergent non-coding (DNC) transcript of unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. Using high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements, we identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC transcription. Nascent transcription profiling showed a genome-wide role of Hda1C in repression of DNC transcription. Live-cell imaging of transcription revealed that mutations in the Hda3 subunit increased the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC transcription regions, supporting DNC transcription repression by histone deacetylation. Our data support the interpretation that DNC transcription results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.