Small RNA library construction from minute biological samples

PBOX-sRNA-seq uncovers novel features of miRNA modification and identifies selected 5'-tRNA fragments bearing 2'-O-modification

The concomitant cloning of RNA degradation products is a major concern in standard small RNA-sequencing practices. This not only complicates the characterization of bona fide sRNAs but also hampers cross-batch experimental replicability and sometimes even results in library construction failure. Given that all types of plant canonical small RNAs possess the 3' end 2'-O-methylation modification, a new small RNA sequencing (sRNA-seq) method, designated as PBOX-sRNA-seq, has been developed specifically to capture this modification. PBOX-sRNA-seq, as its name implies, relies on the sequential treatment of RNA samples with phenylboronic acid-polyacrylamide gel electrophoresis (PBA-PAGE) and sodium periodate (NaIO4) oxidation, before sRNA library construction and sequencing. PBOX-sRNA-seq outperformed separate treatments (i.e. PBA-PAGE only or NaIO4 only) in terms of the depletion of unmethylated RNA species and capture 2'-O-modified sRNAs with extra-high purity. Using PBOX-sRNA-seq, we discovered that nascent miRNA-5p/-3p duplexes may undergo mono-cytidylation/uridylation before 2'-O-methylation. We also identified two highly conserved types of 5'-tRNA fragments (tRF) bearing HEN1-independent 2'-O modification (mainly the 13-nt tRF-5aAla and the 26-nt tRF-5bGly). We believe that PBOX-sRNA-seq is powerful for both qualitative and quantitative analyses of sRNAs in plants and piRNAs in animals.

RNase H-based analysis of synthetic mRNA 5' cap incorporation

Advances in mRNA synthesis and lipid nanoparticles technologies have helped make mRNA therapeutics and vaccines a reality. The 5' cap structure is a crucial modification required to functionalize synthetic mRNA for efficient protein translation in vivo and evasion of cellular innate immune responses. The extent of 5' cap incorporation is one of the critical quality attributes in mRNA manufacturing. RNA cap analysis involves multiple steps: generation of predefined short fragments from the 5' end of the kilobase-long synthetic mRNA molecules using RNase H, a ribozyme or a DNAzyme, enrichment of the 5' cleavage products, and LC-MS intact mass analysis. In this paper, we describe (1) a framework to design site-specific RNA cleavage using RNase H; (2) a method to fluorescently label the RNase H cleavage fragments for more accessible readout methods such as gel electrophoresis or high-throughput capillary electrophoresis; (3) a simplified method for post-RNase H purification using desthiobiotinylated oligonucleotides and streptavidin magnetic beads followed by elution using water. By providing a design framework for RNase H-based RNA 5' cap analysis using less resource-intensive analytical methods, we hope to make RNA cap analysis more accessible to the scientific community.

Improved Boronate Affinity Electrophoresis by Optimization of the Running Buffer for a Single-step Separation of piRNA from Mouse Testis Total RNA

Here we examined optimization of the running buffer in boronate affinity electrophoresis for improved separation of PIWI-interacting RNA (piRNA) with 2'-O-methylated ribose in 3'-terminal nucleotide. The use of Good's buffer, such as HEPES, significantly increased the separation efficiency for piRNA over normal RNA with free 3'-terminal ribose, and retained an ability to resolve the difference by at least 4-nucleotide lengths in the target piRNAs. We also demonstrated a single-step separation of piRNA from mouse testis total RNA.

Boronate affinity electrophoresis for the purification and analysis of cofactor-modified RNAs

RNA modifications are widely distributed in Nature, and their thorough analysis helps answering fundamental biological questions. Nowadays, mass spectrometry or deep-sequencing methods are often used for the analysis. With the raising number of newly discovered RNA modifications, such as the 5'-NAD cap in Escherichia coli, there is an important need for new, less complex and fast analytical tools to analyze the occurrence, amount, and distribution of modified RNAs in cells. To accomplish this task, we have revisited the previously developed affinity gel electrophoresis principles and copolymerized acryloylaminophenyl boronic acid (APB) in standard denaturing polyacrylamide gels to retard the NAD- or FAD-modified RNAs compared to the unmodified RNAs in the gels. The boronyl groups inside the gel form relatively stable complexes with 1,2-cis diols, occurring naturally at the 3'-end of RNA, and also in the nicotinamide riboside of NAD-modified RNA at the 5'-end. The transient formation of diesters between the immobilized boronic acid and the diols causes lower mobility of the modified RNAs, compared to unmodified RNAs, resulting in two distinct bands for one RNA sequence. We used APB affinity gel electrophoresis to preparatively purify in vitro transcribed NAD-RNA from triphosphorylated RNA, to study the enzyme kinetics of the NAD-RNA decapping enzyme NudC, and to determine the NAD modification ratios of various cellular sRNAs. In summary, APB affinity gels can be used to study cofactor-modified RNAs with low amounts of material, and to rapidly screen for their occurrence in total RNA while avoiding complex sample treatments.

Identification of extracellular miRNA in archived serum samples by next-generation sequencing from RNA extracted using multiple methods

miRNAs act as important regulators of gene expression by promoting mRNA degradation or by attenuating protein translation. Since miRNAs are stably expressed in bodily fluids, there is growing interest in profiling these miRNAs, as it is minimally invasive and cost-effective as a diagnostic matrix. A technical hurdle in studying miRNA dynamics is the ability to reliably extract miRNA as small sample volumes and low RNA abundance create challenges for extraction and downstream applications. The purpose of this study was to develop a pipeline for the recovery of miRNA using small volumes of archived serum samples. The RNA was extracted employing several widely utilized RNA isolation kits/methods with and without addition of a carrier. The small RNA library preparation was carried out using Illumina TruSeq small RNA kit and sequencing was carried out using Illumina platform. A fraction of five microliters of total RNA was used for library preparation as quantification is below the detection limit. We were able to profile miRNA levels in serum from all the methods tested. We found out that addition of nucleic acid based carrier molecules had higher numbers of processed reads but it did not enhance the mapping of any miRBase annotated sequences. However, some of the extraction procedures offer certain advantages: RNA extracted by TRIzol seemed to align to the miRBase best; extractions using TRIzol with carrier yielded higher miRNA-to-small RNA ratios. Nuclease free glycogen can be carrier of choice for miRNA sequencing. Our findings illustrate that miRNA extraction and quantification is influenced by the choice of methodologies. Addition of nucleic acid- based carrier molecules during extraction procedure is not a good choice when assaying miRNA using sequencing. The careful selection of an extraction method permits the archived serum samples to become valuable resources for high-throughput applications.

Conserved piRNA Expression from a Distinct Set of piRNA Cluster Loci in Eutherian Mammals

The Piwi pathway is deeply conserved amongst animals because one of its essential functions is to repress transposons. However, many Piwi-interacting RNAs (piRNAs) do not base-pair to transposons and remain mysterious in their targeting function. The sheer number of piRNA cluster (piC) loci in animal genomes and infrequent piRNA sequence conservation also present challenges in determining which piC loci are most important for development. To address this question, we determined the piRNA expression patterns of piC loci across a wide phylogenetic spectrum of animals, and reveal that most genic and intergenic piC loci evolve rapidly in their capacity to generate piRNAs, regardless of known transposon silencing function. Surprisingly, we also uncovered a distinct set of piC loci with piRNA expression conserved deeply in Eutherian mammals. We name these loci Eutherian-Conserved piRNA cluster (ECpiC) loci. Supporting the hypothesis that conservation of piRNA expression across ~100 million years of Eutherian evolution implies function, we determined that one ECpiC locus generates abundant piRNAs antisense to the STOX1 transcript, a gene clinically associated with preeclampsia. Furthermore, we confirmed reduced piRNAs in existing mouse mutations at ECpiC-Asb1 and -Cbl, which also display spermatogenic defects. The Asb1 mutant testes with strongly reduced Asb1 piRNAs also exhibit up-regulated gene expression profiles. These data indicate ECpiC loci may be specially adapted to support Eutherian reproduction.

Animal genomes from flies to humans contain many hundreds of non-coding elements called Piwi-interacting RNAs (piRNAs) cluster loci (piC loci). Some of these elements generate piRNAs that direct the silencing of transposable elements, which are pervasive genetic parasites. However, we lack an understanding of the targeting function for the remaining bulk of piRNAs because their loci are not complementarity to transposable elements. In addition, the field does not know if all piC loci are quickly evolving, or if some piC loci might be deeply conserved in piRNA expression, an indication of its potentially functional importance. Our study confirms the highly rapid evolution in piRNA expression capacity for the majority of piC loci in flies and mammals, with many clade- and species-specific piC loci expression patterns. In spite of this, we also discover a cohort of piC loci that are deeply conserved in piRNA expression from the human to the dog, a significantly broad phylogenetic spectrum of eutherian mammals. However, this conservation in piRNA expression ends at non-eutherian mammals like marsupials and monotremes. Existing mutations in two of these Eutherian-Conserved piC (ECpiC) loci impair mouse reproduction and abrogate piRNA production. Therefore, we suggest these ECpiC loci are conserved for piRNA expression due to their important function in eutherian reproduction and stand out as prime candidates for future genetic studies.

Transposable element dynamics and PIWI regulation impacts lncRNA and gene expression diversity in Drosophila ovarian cell cultures

Piwi proteins and Piwi-interacting RNAs (piRNAs) repress transposable elements (TEs) from mobilizing in gonadal cells. To determine the spectrum of piRNA-regulated targets that may extend beyond TEs, we conducted a genome-wide survey for transcripts associated with PIWI and for transcripts affected by PIWI knockdown in Drosophila ovarian somatic sheet (OSS) cells, a follicle cell line expressing the Piwi pathway. Despite the immense sequence diversity among OSS cell piRNAs, our analysis indicates that TE transcripts are the major transcripts associated with and directly regulated by PIWI. However, several coding genes were indirectly regulated by PIWI via an adjacent de novo TE insertion that generated a nascent TE transcript. Interestingly, we noticed that PIWI-regulated genes in OSS cells greatly differed from genes affected in a related follicle cell culture, ovarian somatic cells (OSCs). Therefore, we characterized the distinct genomic TE insertions across four OSS and OSC lines and discovered dynamic TE landscapes in gonadal cultures that were defined by a subset of active TEs. Particular de novo TEs appeared to stimulate the expression of novel candidate long noncoding RNAs (lncRNAs) in a cell lineage-specific manner, and some of these TE-associated lncRNAs were associated with PIWI and overlapped PIWI-regulated genes. Our analyses of OSCs and OSS cells demonstrate that despite having a Piwi pathway to suppress endogenous mobile elements, gonadal cell TE landscapes can still dramatically change and create transcriptome diversity.

Quantitative genome-wide enhancer activity maps for five Drosophila species show functional enhancer conservation and turnover during cis-regulatory evolution

Phenotypic differences between closely related species are thought to arise primarily from changes in gene expression due to mutations in cis-regulatory sequences (enhancers). However, it has remained unclear how frequently mutations alter enhancer activity or create functional enhancers de novo. Here we use STARR-seq, a recently developed quantitative enhancer assay, to determine genome-wide enhancer activity profiles for five Drosophila species in the constant trans-regulatory environment of Drosophila melanogaster S2 cells. We find that the functions of a large fraction of D. melanogaster enhancers are conserved for their orthologous sequences owing to selection and stabilizing turnover of transcription factor motifs. Moreover, hundreds of enhancers have been gained since the D. melanogaster-Drosophila yakuba split about 11 million years ago without apparent adaptive selection and can contribute to changes in gene expression in vivo. Our finding that enhancer activity is often deeply conserved and frequently gained provides functional insights into regulatory evolution.