Benchmarking sequencing methods and tools that facilitate the study of alternative polyadenylation

BACKGROUND

Alternative cleavage and polyadenylation (APA), an RNA processing event, occurs in over 70% of human protein-coding genes. APA results in mRNA transcripts with distinct 3' ends. Most APA occurs within 3' UTRs, which harbor regulatory elements that can impact mRNA stability, translation, and localization.

RESULTS

APA can be profiled using a number of established computational tools that infer polyadenylation sites from standard, short-read RNA-seq datasets. Here, we benchmarked a number of such tools-TAPAS, QAPA, DaPars2, GETUTR, and APATrap- against 3'-Seq, a specialized RNA-seq protocol that enriches for reads at the 3' ends of genes, and Iso-Seq, a Pacific Biosciences (PacBio) single-molecule full-length RNA-seq method in their ability to identify polyadenylation sites and quantify polyadenylation site usage. We demonstrate that 3'-Seq and Iso-Seq are able to identify and quantify the usage of polyadenylation sites more reliably than computational tools that take short-read RNA-seq as input. However, we find that running one such tool, QAPA, with a set of polyadenylation site annotations derived from small quantities of 3'-Seq or Iso-Seq can reliably quantify variation in APA across conditions, such asacross genotypes, as demonstrated by the successful mapping of alternative polyadenylation quantitative trait loci (apaQTL).

CONCLUSIONS

We envisage that our analyses will shed light on the advantages of studying APA with more specialized sequencing protocols, such as 3'-Seq or Iso-Seq, and the limitations of studying APA with short-read RNA-seq. We provide a computational pipeline to aid in the identification of polyadenylation sites and quantification of polyadenylation site usages using Iso-Seq data as input.

Neuregulin1 and ErbB expression in the uninjured and regenerating olfactory mucosa

Neuregulin1, a protein involved in signaling through the ErbB receptors, is required for the proper development of multiple organ systems. A complete understanding of the expression profile of Neuregulin1 is complicated by the presence of multiple isoform variants that result from extensive alternative splicing. Remarkably, these numerous protein products display a wide range of divergent functional roles, making the characterization of tissue-specific isoforms critical to understanding signaling. Recent evidence suggests an important role for Neuregulin1 signaling during olfactory epithelium development and regeneration. In order to understand the physiological consequences of this signaling, we sought to identify the isoform-specific and cell type-specific expression pattern of Neuregulin1 in the adult olfactory mucosa using a combination of RT-qPCR, FACS, and immunohistochemistry. To complement this information, we also analyzed the cell-type specific expression patterns of the ErbB receptors using immunohistochemistry. We found that multiple Neuregulin1 isoforms, containing predominantly the Type I and Type III N-termini, are expressed in the uninjured olfactory mucosa. Specifically, we found that Type III Neuregulin1 is highly expressed in mature olfactory sensory neurons and Type I Neuregulin1 is highly expressed in duct gland cells. Surprisingly, the divergent localization of these Neuregulin isoforms and their corresponding ErbB receptors does not support a role for active signaling during normal turnover and maintenance of the olfactory mucosa. Conversely, we found that injury to the olfactory epithelium specifically upregulates the Neuregulin1 Type I isoform bringing the expression pattern adjacent to cells expressing both ErbB2 and ErbB3 which is compatible with active signaling, supporting a functional role for Neuregulin1 specifically during regeneration.