Ribosome profiling reveals translational regulation of mammalian cells in response to hypoxic stress

Nuclear release of eIF1 restricts start-codon selection during mitosis

Regulated start-codon selection has the potential to reshape the proteome through the differential production of upstream open reading frames, canonical proteins, and alternative translational isoforms1-3. However, conditions under which start codon selection is altered remain poorly defined. Here, using transcriptome-wide translation-initiation-site profiling4, we reveal a global increase in the stringency of start-codon selection during mammalian mitosis. Low-efficiency initiation sites are preferentially repressed in mitosis, resulting in pervasive changes in the translation of thousands of start sites and their corresponding protein products. This enhanced stringency of start-codon selection during mitosis results from increased association between the 40S ribosome and the key regulator of start-codon selection, eIF1. We find that increased eIF1-40S ribosome interaction during mitosis is mediated by the release of a nuclear pool of eIF1 upon nuclear envelope breakdown. Selectively depleting the nuclear pool of eIF1 eliminates the change to translational stringency during mitosis, resulting in altered synthesis of thousands of protein isoforms. In addition, preventing mitotic translational rewiring results in substantially increased cell death and decreased mitotic slippage in cells that experience a mitotic delay induced by anti-mitotic chemotherapies. Thus, cells globally control stringency of translation initiation, which has critical roles during the mammalian cell cycle in preserving mitotic cell physiology.

RNA G-quadruplexes and stress: emerging mechanisms and functions

RNA G-quadruplexes (rG4s) are noncanonical secondary structures formed by guanine-rich sequences that are found in different regions of RNA molecules. These structures have been implicated in diverse biological processes, including translation, splicing, and RNA stability. Recent studies have suggested that rG4s play a role in the cellular response to stress. This review summarizes the current knowledge on rG4s under stress, focusing on their formation, regulation, and potential functions in stress response pathways. We discuss the molecular mechanisms that regulate the formation of rG4 under different stress conditions and the impact of these structures on RNA metabolism, gene expression, and cell survival. Finally, we highlight the potential therapeutic implications of targeting rG4s for the treatment of stress-related diseases through modulating cell survival.

Ribosome Profiling and RNA Sequencing Reveal Translation and Transcription Regulation under Acute Heat Stress in Rainbow Trout (Oncorhynchus mykiss, Walbaum, 1792) Liver

Rainbow trout (Oncorhynchus mykiss, Walbaum, 1792) is an important economic cold-water fish that is susceptible to heat stress. To date, the heat stress response in rainbow trout is more widely understood at the transcriptional level, while little research has been conducted at the translational level. To reveal the translational regulation of heat stress in rainbow trout, in this study, we performed a ribosome profiling assay of rainbow trout liver under normal and heat stress conditions. Comparative analysis of the RNA-seq data with the ribosome profiling data showed that the folding changes in gene expression at the transcriptional level are moderately correlated with those at the translational level. In total, 1213 genes were significantly altered at the translational level. However, only 32.8% of the genes were common between both levels, demonstrating that heat stress is coordinated across both transcriptional and translational levels. Moreover, 809 genes exhibited significant differences in translational efficiency (TE), with the TE of these genes being considerably affected by factors such as the GC content, coding sequence length, and upstream open reading frame (uORF) presence. In addition, 3468 potential uORFs in 2676 genes were identified, which can potentially affect the TE of the main open reading frames. In this study, Ribo-seq and RNA-seq were used for the first time to elucidate the coordinated regulation of transcription and translation in rainbow trout under heat stress. These findings are expected to contribute novel data and theoretical insights to the international literature on the thermal stress response in fish.

Nuclear release of eIF1 globally increases stringency of start-codon selection to preserve mitotic arrest physiology

Regulated start-codon selection has the potential to reshape the proteome through the differential production of uORFs, canonical proteins, and alternative translational isoforms. However, conditions under which start-codon selection is altered remain poorly defined. Here, using transcriptome-wide translation initiation site profiling, we reveal a global increase in the stringency of start-codon selection during mammalian mitosis. Low-efficiency initiation sites are preferentially repressed in mitosis, resulting in pervasive changes in the translation of thousands of start sites and their corresponding protein products. This increased stringency of start-codon selection during mitosis results from increased interactions between the key regulator of start-codon selection, eIF1, and the 40S ribosome. We find that increased eIF1-40S ribosome interactions during mitosis are mediated by the release of a nuclear pool of eIF1 upon nuclear envelope breakdown. Selectively depleting the nuclear pool of eIF1 eliminates the changes to translational stringency during mitosis, resulting in altered mitotic proteome composition. In addition, preventing mitotic translational rewiring results in substantially increased cell death and decreased mitotic slippage following treatment with anti-mitotic chemotherapeutics. Thus, cells globally control translation initiation stringency with critical roles during the mammalian cell cycle to preserve mitotic cell physiology.

Executable Network Models of Integrated Multiomics Data

Multiomics profiling provides a holistic picture of a condition being examined and captures the complexity of signaling events, beginning from the original cause (environmental or genetic), to downstream functional changes at multiple molecular layers. Pathway enrichment analysis has been used with multiomics data sets to characterize signaling mechanisms. However, technical and biological variability between these layered data limit an integrative computational analyses. We present a Boolean network-based method, multiomics Boolean Omics Network Invariant-Time Analysis (mBONITA), to integrate omics data sets that quantify multiple molecular layers. mBONITA utilizes prior knowledge networks to perform topology-based pathway analysis. In addition, mBONITA identifies genes that are consistently modulated across molecular measurements by combining observed fold-changes and variance, with a measure of node (i.e., gene or protein) influence over signaling, and a measure of the strength of evidence for that gene across data sets. We used mBONITA to integrate multiomics data sets from RAMOS B cells treated with the immunosuppressant drug cyclosporine A under varying O₂ tensions to identify pathways involved in hypoxia-mediated chemotaxis. We compare mBONITA's performance with 6 other pathway analysis methods designed for multiomics data and show that mBONITA identifies a set of pathways with evidence of modulation across all omics layers. mBONITA is freely available at https://github.com/Thakar-Lab/mBONITA.

MDM2 inhibitors, nutlin-3a and navtemadelin, retain efficacy in human and mouse cancer cells cultured in hypoxia

Activation of p53 by small molecule MDM2 inhibitors can induce cell cycle arrest or death in p53 wildtype cancer cells. However, cancer cells exposed to hypoxia can develop resistance to other small molecules, such as chemotherapies, that activate p53. Here, we evaluated whether hypoxia could render cancer cells insensitive to two MDM2 inhibitors with different potencies, nutlin-3a and navtemadlin. Inhibitor efficacy and potency were evaluated under short-term hypoxic conditions in human and mouse cancer cells expressing different p53 genotypes (wild-type, mutant, or null). Treatment of wild-type p53 cancer cells with MDM2 inhibitors reduced cell growth by > 75% in hypoxia through activation of the p53-p21 signaling pathway; no inhibitor-induced growth reduction was observed in hypoxic mutant or null p53 cells except at very high concentrations. The concentration of inhibitors needed to induce the maximal p53 response was not significantly different in hypoxia compared to normoxia. However, inhibitor efficacy varied by species and by cell line, with stronger effects at lower concentrations observed in human cell lines than in mouse cell lines grown as 2D and 3D cultures. Together, these results indicate that MDM2 inhibitors retain efficacy in hypoxia, suggesting they could be useful for targeting acutely hypoxic cancer cells.

Ribosome profiling analysis reveals the roles of DDX41 in translational regulation

DDX41 mutation has been observed in myeloid malignancies including myelodysplastic syndromes and acute myeloid leukemia, but the underlying causative mechanisms of these diseases have not been fully elucidated. The DDX41 protein is an ATP-dependent RNA helicase with roles in RNA metabolism. We previously showed that DDX41 is involved in ribosome biogenesis by promoting the processing of newly transcribed pre-ribosomal RNA. To build on this finding, in this study, we leveraged ribosome profiling technology to investigate the involvement of DDX41 in translation. We found that DDX41 knockdown resulted in both translationally increased and decreased transcripts. Both gene set enrichment analysis and gene ontology analysis indicated that ribosome-associated genes were translationally promoted after DDX41 knockdown, in part because these transcripts had significantly shorter transcript length and higher transcriptional and translational levels. In addition, we found that transcripts with 5'-terminal oligopyrimidine motifs tended to be translationally upregulated when the DDX41 level was low. Our data suggest that a translationally regulated feedback mechanism involving DDX41 may exist for ribosome biogenesis.

Ciliogenesis requires sphingolipid-dependent membrane and axoneme interaction

Cilium formation and regeneration requires new protein synthesis, but the underlying cytosolic translational reprogramming remains largely unknown. Using ribosome footprinting, we performed global translatome profiling during cilia regeneration in Chlamydomonas and uncovered that flagellar genes undergo an early transcriptional activation but late translational repression. This pattern guided our identification of sphingolipid metabolism enzymes, including serine palmitoyltransferase (SPT), as essential regulators for ciliogenesis. Cryo-electron tomography showed that ceramide loss abnormally increased the membrane-axoneme distance and generated bulged cilia. We found that ceramides interact with intraflagellar transport (IFT) particle proteins that IFT motors transport along axoneme microtubules (MTs), suggesting that ceramide-IFT particle-IFT motor-MT interactions connect the ciliary membrane with the axoneme to form rod-shaped cilia. SPT-deficient vertebrate cells were defective in ciliogenesis, and SPT mutations from patients with hereditary sensory neuropathy disrupted cilia, which could be restored by sphingolipid supplementation. These results reveal a conserved role of sphingolipid in cilium formation and link compromised sphingolipid production with ciliopathies.

Ribosome Profiling Reveals Genome-Wide Cellular Translational Regulation in Lacticaseibacillus rhamnosus ATCC 53103 under Acid Stress

During fermentation and food processing, Lacticaseibacillus rhamnosus ATCC 53103 can encounter many adverse conditions, and acid stress is one of them. The purpose of the present study was to investigate the influence of acid stress on the global translational and transcriptional regulation of Lacticaseibacillus rhamnosus ATCC 53103. Two pH values (pH 6.0 vs. pH 5.0) were applied, the effects of which were studied via ribosome profiling and RNA sequencing assay. Under acid stress, many genes showed differential changes at the translational and transcriptional levels. A total of 10 genes showed different expression trends at the two levels. The expression of 337 genes—which mainly participated in the ABC transporters, amino acid metabolism, and ribosome functional group assembly pathways—was shown to be regulated only at the translational level. The translational efficiency of a few genes participating in the pyrimidine and amino acid metabolism pathways were upregulated. Ribosome occupancy data suggested that ribosomes accumulated remarkably in the elongation region of open reading frame regions under acid stress. This study provides new insights into Lacticaseibacillus rhamnosus ATCC 53103 gene expression under acid stress, and demonstrates that the bacterium can respond to acid stress with synergistic translational and transcriptional regulation mechanisms, improving the vitality of cells.

Uncontrolled Oxygen Levels in Cultures of Retinal Pigment Epithelium: Have We Missed the Obvious?

Retinal pigment epithelium (RPE) is the outermost layer of retina located between the photoreceptor cells and the choroid. This highly-polarized monolayer provides critical support for the functioning of the other parts of the retina, especially photoreceptors. Methods of culturing RPE have been under development since its establishment in 1920s. Despite considering various factors, oxygen (O₂) levels in RPE microenvironments during culture preparation and experimental procedure have been overlooked. O₂ is a crucial parameter in the cultures, and therefore, maintaining RPE cells at O₂ levels different from their native environment (70-90 mm Hg of O₂) could have unintended consequences. Owing to the importance of the topic, lack of sufficient discussion in the literature and to encourage future research, this paper will focus on uncontrolled O₂ level in cultures of RPE cells.

Quantitative Translation Proteomics Using mePROD

Multiplexed enhanced protein dynamic mass spectrometry (mePROD MS) enables robust quantification of translation in cell culture. Tandem mass tags (TMT) are combined with pulsed stable isotope labeling in cell culture (pSILAC) to monitor newly synthesized proteins on a proteome wide scale. While approaches combining pSILAC and TMT typically require long labeling times to reach sufficient intensity of the newly synthesized peptides in the mass spectrometer, mePROD uses a carrier signal that boosts the survey scan intensity and strongly increases identification rates. Hence, this protocol provides an easy and cost-efficient method to profile proteome-wide translatome changes at a temporal resolution of minutes.

Translational control in cell ageing: an update

Cellular ageing is one of the main drivers of organismal ageing and holds keys towards improving the longevity and quality of the extended life. Elucidating mechanisms underlying the emergence of the aged cells as well as their altered responses to the environment will help understanding the evolutionarily defined longevity preferences across species with different strategies of survival. Much is understood about the role of alterations in the DNA, including many epigenetic modifications such as methylation, in relation to the aged cell phenotype. While transcriptomes of the aged cells are beginning to be better-characterised, their translational responses remain under active investigation. Many of the translationally controlled homeostatic pathways are centred around mitigation of DNA damage, cell stress response and regulation of the proliferative potential of the cells, and thus are critical for the aged cell function. Translation profiling-type studies have boosted the opportunities in discovering the function of protein biosynthesis control and are starting to be applied to the aged cells. Here, we provide a summary of the current knowledge about translational mechanisms considered to be commonly altered in the aged cells, including the integrated stress response-, mechanistic target of Rapamycin- and elongation factor 2 kinase-mediated pathways. We enlist and discuss findings of the recent works that use broad profiling-type approaches to investigate the age-related translational pathways. We outline the limitations of the methods and the remaining unknowns in the established ageing-associated translation mechanisms, and flag translational mechanisms with high prospective importance in ageing, for future studies.

Long RNA Sequencing and Ribosome Profiling of Inflamed β-Cells Reveal an Extensive Translatome Landscape

Type 1 diabetes (T1D) is an autoimmune disease characterized by autoreactive T cell-mediated destruction of the insulin-producing pancreatic β-cells. Increasing evidence suggest that the β-cells themselves contribute to their own destruction by generating neoantigens through the production of aberrant or modified proteins that escape central tolerance. We recently demonstrated that ribosomal infidelity amplified by stress could lead to the generation of neoantigens in human β-cells, emphasizing the participation of nonconventional translation events in autoimmunity, as occurring in cancer or virus-infected tissues. Using a transcriptome-wide profiling approach to map translation initiation start sites in human β-cells under standard and inflammatory conditions, we identify a completely new set of polypeptides derived from noncanonical start sites and translation initiation within long noncoding RNA. Our data underline the extreme diversity of the β-cell translatome and may reveal new functional biomarkers for β-cell distress, disease prediction and progression, and therapeutic intervention in T1D.

Most non-canonical proteins uniquely populate the proteome or immunopeptidome

Combining RNA sequencing, ribosome profiling, and mass spectrometry, we elucidate the contribution of non-canonical translation to the proteome and major histocompatibility complex (MHC) class I immunopeptidome. Remarkably, of 14,498 proteins identified in three human B cell lymphomas, 2,503 are non-canonical proteins. Of these, 28% are novel isoforms and 72% are cryptic proteins encoded by ostensibly non-coding regions (60%) or frameshifted canonical genes (12%). Cryptic proteins are translated as efficiently as canonical proteins, have more predicted disordered residues and lower stability, and critically generate MHC-I peptides 5-fold more efficiently per translation event. Translating 5' "untranslated" regions hinders downstream translation of genes involved in transcription, translation, and antiviral responses. Novel protein isoforms show strong enrichment for signaling pathways deregulated in cancer. Only a small fraction of cryptic proteins detected in the proteome contribute to the MHC-I immunopeptidome, demonstrating the high preferential access of cryptic defective ribosomal products to the class I pathway.

Insights Into Translatomics in the Nervous System

Most neurological disorders are caused by abnormal gene translation. Generally, dysregulation of elements involved in the translational process disrupts homeostasis in neurons and neuroglia. Better understanding of how the gene translation process occurs requires detailed analysis of transcriptomic and proteomic profile data. However, a lack of strictly direct correlations between mRNA and protein levels limits translational investigation by combining transcriptomic and proteomic profiling. The much better correlation between proteins and translated mRNAs than total mRNAs in abundance and insufficiently sensitive proteomics approach promote the requirement of advances in translatomics technology. Translatomics which capture and sequence the mRNAs associated with ribosomes has been effective in identifying translational changes by genetics or projections, ribosome stalling, local translation, and transcript isoforms in the nervous system. Here, we place emphasis on the main three translatomics methods currently used to profile mRNAs attached to ribosome-nascent chain complex (RNC-mRNA). Their prominent applications in neurological diseases including glioma, neuropathic pain, depression, fragile X syndrome (FXS), neurodegenerative disorders are outlined. The content reviewed here expands our understanding on the contributions of aberrant translation to neurological disease development.

EPAS1 and VEGFA gene variants are related to the symptoms of acute mountain sickness in Chinese Han population: a cross-sectional study

BACKGROUND

More people ascend to high altitude (HA) for various activities, and some individuals are susceptible to HA illness after rapidly ascending from plains. Acute mountain sickness (AMS) is a general complaint that affects activities of daily living at HA. Although genomic association analyses suggest that single nucleotide polymorphisms (SNPs) are involved in the genesis of AMS, no major gene variants associated with AMS-related symptoms have been identified.

METHODS

In this cross-sectional study, 604 young, healthy Chinese Han men were recruited in June and July of 2012 in Chengdu, and rapidly taken to above 3700 m by plane. Basic demographic parameters were collected at sea level, and heart rate, pulse oxygen saturation (SpO₂), systolic and diastolic blood pressure and AMS-related symptoms were determined within 18-24 h after arriving in Lhasa. AMS patients were identified according to the latest Lake Louise scoring system (LLSS). Potential associations between variant genotypes and AMS/AMS-related symptoms were identified by logistic regression after adjusting for potential confounders (age, body mass index and smoking status).

RESULTS

In total, 320 subjects (53.0%) were diagnosed with AMS, with no cases of high-altitude pulmonary edema or high-altitude cerebral edema. SpO₂ was significantly lower in the AMS group than that in the non-AMS group (P = 0.003). Four SNPs in hypoxia-inducible factor-related genes were found to be associated with AMS before multiple hypothesis testing correction. The rs6756667 (EPAS1) was associated with mild gastrointestinal symptoms (P = 0.013), while rs3025039 (VEGFA) was related to mild headache (P = 0.0007). The combination of rs6756667 GG and rs3025039 CT/TT further increased the risk of developing AMS (OR = 2.70, P < 0.001).

CONCLUSIONS

Under the latest LLSS, we find that EPAS1 and VEGFA gene variants are related to AMS susceptibility through different AMS-related symptoms in the Chinese Han population; this tool might be useful for screening susceptible populations and predicting clinical symptoms leading to AMS before an individual reaches HA.

TRIAL REGISTRATION

Chinese Clinical Trial Registration, ChiCTR-RCS-12002232 . Registered 31 May 2012.

Genomic and physiological resilience in extreme environments are associated with a secure attachment style

Understanding individual capability to adjust to protracted confinement and isolation may inform adaptive plasticity and disease vulnerability/resilience, and may have long-term implications for operations requiring prolonged presence in distant and restricted environments. Individual coping depends on many different factors encompassing psychological dispositional traits, endocrine reactivity and their underlying molecular mechanisms (e.g. gene expression). A positive view of self and others (secure attachment style) has been proposed to promote individual resilience under extreme environmental conditions. Here, we tested this hypothesis and investigated the underlying molecular mechanisms in 13 healthy volunteers confined and isolated for 12 months in a research station located 1670 km away from the south geographic pole on the Antarctic Plateau at 3233 m above sea level. Study participants, stratified for attachment style, were characterised longitudinally (before, during and after confinement) for their psychological appraisal of the stressful nature of the expedition, diurnal fluctuations in endocrine stress reactivity, and gene expression profiling (transcriptomics). Predictably, a secure attachment style was associated with reduced psychological distress and endocrine vulnerability to stress. In addition, while prolonged confinement and isolation remarkably altered overall patterns of gene expression, such alteration was largely reduced in individuals characterised by a secure attachment style. Furthermore, increased resilience was associated with a reduced expression of genes involved in energy metabolism (mitochondrial function and oxidative phosphorylation). Ultimately, our data indicate that a secure attachment style may favour individual resilience in extreme environments and that such resilience can be mapped onto identifiable molecular substrates.

Functional Translatome Proteomics Reveal Converging and Dose-Dependent Regulation by mTORC1 and eIF2α

Regulation of translation is essential during stress. However, the precise sets of proteins regulated by the key translational stress responses-the integrated stress response (ISR) and mTORC1-remain elusive. We developed multiplexed enhanced protein dynamics (mePROD) proteomics, adding signal amplification to dynamic-SILAC and multiplexing, to enable measuring acute changes in protein synthesis. Treating cells with ISR/mTORC1-modulating stressors, we showed extensive translatome modulation with ∼20% of proteins synthesized at highly reduced rates. Comparing translation-deficient sub-proteomes revealed an extensive overlap demonstrating that target specificity is achieved on protein level and not by pathway activation. Titrating cap-dependent translation inhibition confirmed that synthesis of individual proteins is controlled by intrinsic properties responding to global translation attenuation. This study reports a highly sensitive method to measure relative translation at the nascent chain level and provides insight into how the ISR and mTORC1, two key cellular pathways, regulate the translatome to guide cellular survival upon stress.

Multimapping confounds ribosome profiling analysis: A case-study of the Hsp90 molecular chaperone

Ribosome profiling (Ribo-seq) can potentially provide detailed information about ribosome position on transcripts and estimates of protein translation levels in vivo. Hsp90 chaperones, which play a critical role in stress tolerance, have characteristic patterns of differential expression under nonstressed and heat shock conditions. By analyzing published Ribo-seq data for the Hsp90 chaperones in S. cerevisiae, we find wide-ranging artifacts originating from "multimapping" reads (reads that cannot be uniquely assigned to one position), which constitute ~25% of typical S. cerevisiae Ribo-seq datasets and ~80% of the reads from HEK293 cells. Estimates of Hsp90 protein production as determined by Ribo-seq are reproducible but not robust, with inferred expression levels that can change 10-fold depending on how multimapping reads are processed. The differential expression of Hsp90 chaperones under nonstressed and heat shock conditions creates artificial peaks and valleys in their ribosome profiles that give a false impression of regulated translational pausing. Indeed, we find that multimapping can even create an appearance of reproducibility to the shape of the Hsp90 ribosome profiles from biological replicates. Adding further complexity, this artificial reproducibility is dependent on the computational method used to construct the ribosome profile. Given the ubiquity of multimapping reads in Ribo-seq experiments and the complexity of artifacts associated with multimapping, we developed a publicly available computational tool to identify transcripts most at risk for multimapping artifacts. In doing so, we identify biological pathways that are enriched in multimapping transcripts, meaning that particular biological pathways will be highly susceptible to multimapping artifacts.

Translational and transcriptional responses in human primary hepatocytes under hypoxia

The liver is the primary source of a large number of plasma proteins and plays a critical role in multiple biological processes. Inadequate oxygen supply characterizing various clinical settings such as liver transplantation exposes the liver to hypoxic conditions. Studies assessing hypoxia-induced global translational changes in liver are lacking. Here, we employed a recently developed ribosome-profiling technique to assess global translational responses of human primary hepatocytes exposed to acute hypoxic stress (1% O₂) for the short term. In parallel, transcriptome profiling was performed to assess mRNA expression changes. We found that translational responses appeared earlier and were predominant over transcriptional responses. A significant decrease in translational efficiency of several ribosome genes indicated translational inhibition of new ribosome protein synthesis in hypoxia. Pathway enrichment analysis highlighted altered translational regulation of MAPK signaling, drug metabolism, oxidative phosphorylation, and nonalcoholic fatty liver disease pathways. Gene Ontology enrichment analysis revealed terms related to translation, metabolism, angiogenesis, apoptosis, and response to stress. Transcriptional induction of genes encoding heat shock proteins was observed within 30 min of hypoxia. Induction of genes encoding stress response mediators, metabolism regulators, and proangiogenic proteins was observed at 240 min. Despite the liver being the primary source of coagulation proteins and the implicated role of hypoxia in thrombosis, limited differences were observed in genes encoding coagulation-associated proteins. Overall, our study demonstrates the predominance of translational regulation over transcription and highlights differentially regulated pathways or biological processes in short-term hypoxic stress responses of human primary hepatocytes. NEW & NOTEWORTHY The novelty of this study lies in applying parallel ribosome- and transcriptome-profiling analyses to human primary hepatocytes in hypoxia. To our knowledge, this is the first study to assess global translational responses using ribosome profiling in hypoxic hepatocytes. Our results demonstrate the predominance of translational responses over transcriptional responses in early hepatic hypoxic stress responses. Furthermore, our study reveals multiple pathways and specific genes showing altered regulation in hypoxic hepatocytes.

ALFA: annotation landscape for aligned reads

Background

The last 10 years have seen the rise of countless functional genomics studies based on Next-Generation Sequencing (NGS). In the vast majority of cases, whatever the species, whatever the experiment, the two first steps of data analysis consist of a quality control of the raw reads followed by a mapping of those reads to a reference genome/transcriptome. Subsequent steps then depend on the type of study that is being made. While some tools have been proposed for investigating data quality after the mapping step, there is no commonly adopted framework that would be easy to use and broadly applicable to any NGS data type.

Results

We present ALFA, a simple but universal tool that can be used after the mapping step on any kind of NGS experiment data for any organism with available genomic annotations. In a single command line, ALFA can compute and display distribution of reads by categories (exon, intron, UTR, etc.) and biotypes (protein coding, miRNA, etc.) for a given aligned dataset with nucleotide precision. We present applications of ALFA to Ribo-Seq and RNA-Seq on Homo sapiens, CLIP-Seq on Mus musculus, RNA-Seq on Saccharomyces cerevisiae, Bisulfite sequencing on Arabidopsis thaliana and ChIP-Seq on Caenorhabditis elegans.

Conclusions

We show that ALFA provides a powerful and broadly applicable approach for post mapping quality control and to produce a global overview using common or dedicated annotations. It is made available to the community as an easy to install command line tool and from the Galaxy Tool Shed.