Polyadenylation of ribosomal RNA in human cells

Brainstem transcriptomic changes in male Wistar rats after acute stress, comparing the use of duplex specific nuclease (DSN)

In this work, we have analyzed the transcriptomic changes in the brainstem of male Wistar rats 2 h after an acute stress exposure. We performed duplex-specific nuclease normalization of cDNA libraries and compared the results back-to-back for the first time. Based on our RNAseq data, we selected reference genes for RT-qPCR that are best suited for acute stress experiments. Most genes were upregulated. We detected a massive shift in neuropeptide Crh, Trh,Cga, Tshb, Uts2b, Tac4, Lep and neuropeptide receptor Hcrtr1, Sstr5, Bdkrb2, Crhr2 signaling, as well as glutamate Grin3b, Grm2 and GABA Gpr156, acetylcholine Chrm4,Chrne, adrenergic Adra2b receptors expression. A strong increase in the expression of intermediate filaments Krt83/Krt86/Krt80/Krt84/Krt87/Krt4/Krt76 and motor proteins Myo7a, Klc3 was detected. Remarkably, in the absence of astrocyte activation, we also observed signs of microglial activation at this time point. Both expression of anti-inflammatory cytokines Il13, Ccl24 and pro-inflammatory cytokine receptors Il9r, Il12rb1, Tnfrsf14, Tnfrsf13c, Tnfrsf25, Tnfrsf1b were increased. In the Wnt signaling pathway, we observed increased expression of ligands-receptors Wnt1, Wnt11, Ror2 and also negative regulators Notum, Sfrp5, Sost. RNAseq results after DSN treatment correlated at a high level with RNAseq results without DSN, but there was a proportion of genes that shifted their logFC values. They are mostly rare transcripts TPM 1-10 with higher 0.5-0.9 GC content.

RNA tailing machinery drives amyloidogenic phase transition

The RNA tailing machinery adds nucleotides to the 3'-end of RNA molecules that are implicated in various biochemical functions, including protein synthesis and RNA stability. Here, we report a role for the RNA tailing machinery as enzymatic modifiers of intracellular amyloidogenesis. A targeted RNA interference screen identified Terminal Nucleotidyl-transferase 4b (TENT4b/Papd5) as an essential participant in the amyloidogenic phase transition of nucleoli into solid-like Amyloid bodies. Full-length-and-mRNA sequencing uncovered starRNA, a class of unusually long untemplated RNA molecules synthesized by TENT4b. StarRNA consists of short rRNA fragments linked to long, linear mixed tails that operate as polyanionic stimulators of amyloidogenesis in cells and in vitro. Ribosomal intergenic spacer noncoding RNA (rIGSRNA) recruit TENT4b in intranucleolar foci to coordinate starRNA synthesis driving their amyloidogenic phase transition. The exoribonuclease RNA Exosome degrades starRNA and functions as a general suppressor of cellular amyloidogenesis. We propose that amyloidogenic phase transition is under tight enzymatic control by the RNA tailing and exosome axis.

Disruption of mitochondrial energy metabolism is a putative pathogenesis of Diamond-Blackfan anemia

Energy metabolism in the context of erythropoiesis and related diseases remains largely unexplored. Here, we developed a primary cell model by differentiating hematopoietic stem progenitor cells toward the erythroid lineage and suppressing the mitochondrial oxidative phosphorylation (OXPHOS) pathway. OXPHOS suppression led to differentiation failure of erythroid progenitors and defects in ribosome biogenesis. Ran GTPase-activating protein 1 (RanGAP1) was identified as a target of mitochondrial OXPHOS for ribosomal defects during erythropoiesis. Overexpression of RanGAP1 largely alleviated erythroid defects resulting from OXPHOS suppression. Coenzyme Q10, an activator of OXPHOS, largely rescued erythroid defects and increased RanGAP1 expression. Patients with Diamond-Blackfan anemia (DBA) exhibited OXPHOS suppression and a concomitant suppression of ribosome biogenesis. RNA-seq analysis implied that the substantial mutation (approximately 10%) in OXPHOS genes accounts for OXPHOS suppression in these patients. Conclusively, OXPHOS disruption and the associated disruptive mitochondrial energy metabolism are linked to the pathogenesis of DBA.

Counter-regulation of RNA stability by UPF1 and TDP43

RNA quality control is crucial for proper regulation of gene expression. Disruption of nonsense mediated mRNA decay (NMD), the primary RNA decay pathway responsible for the degradation of transcripts containing premature termination codons (PTCs), can disrupt development and lead to multiple diseases in humans and other animals. Similarly, therapies targeting NMD may have applications in hematological, neoplastic and neurological disorders. As such, tools capable of accurately quantifying NMD status could be invaluable for investigations of disease pathogenesis and biomarker identification. Toward this end, we assemble, validate, and apply a next-generation sequencing approach (NMDq) for identifying and measuring the abundance of PTC-containing transcripts. After validating NMDq performance and confirming its utility for tracking RNA surveillance, we apply it to determine pathway activity in two neurodegenerative diseases, amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) characterized by RNA misprocessing and abnormal RNA stability. Despite the genetic and pathologic evidence implicating dysfunctional RNA metabolism, and NMD in particular, in these conditions, we detected no significant differences in PTC-encoding transcripts in ALS models or disease. Contrary to expectations, overexpression of the master NMD regulator UPF1 had little effect on the clearance of transcripts with PTCs, but rather restored RNA homeostasis through differential use and decay of alternatively poly-adenylated isoforms. Together, these data suggest that canonical NMD is not a significant contributor to ALS/FTD pathogenesis, and that UPF1 promotes neuronal survival by regulating transcripts with abnormally long 3'UTRs.

Profiling the polyadenylated transcriptome of extracellular vesicles with long-read nanopore sequencing

BACKGROUND

While numerous studies have described the transcriptomes of extracellular vesicles (EVs) in different cellular contexts, these efforts have typically relied on sequencing methods requiring RNA fragmentation, which limits interpretations on the integrity and isoform diversity of EV-targeted RNA populations. It has been assumed that mRNA signatures in EVs are likely to be fragmentation products of the cellular mRNA material, and the extent to which full-length mRNAs are present within EVs remains to be clarified.

RESULTS

Using long-read nanopore RNA sequencing, we sought to characterize the full-length polyadenylated (poly-A) transcriptome of EVs released by human chronic myelogenous leukemia K562 cells. We detected 443 and 280 RNAs that were respectively enriched or depleted in EVs. EV-enriched poly-A transcripts consist of a variety of biotypes, including mRNAs, long non-coding RNAs, and pseudogenes. Our analysis revealed that 10.58% of all EV reads, and 18.67% of all cellular (WC) reads, corresponded to known full-length transcripts, with mRNAs representing the largest biotype for each group (EV = 58.13%, WC = 43.93%). We also observed that for many well-represented coding and non-coding genes, diverse full-length transcript isoforms were present in EV specimens, and these isoforms were reflective-of but often in different ratio compared to cellular samples.

CONCLUSION

This work provides novel insights into the compositional diversity of poly-A transcript isoforms enriched within EVs, while also underscoring the potential usefulness of nanopore sequencing to interrogate secreted RNA transcriptomes.

Non-coding RNAs in immunoregulation and autoimmunity: Technological advances and critical limitations

Immune cell function is critically dependent on precise control over transcriptional output from the genome. In this respect, integration of environmental signals that regulate gene expression, specifically by transcription factors, enhancer DNA elements, genome topography and non-coding RNAs (ncRNAs), are key components. The first three have been extensively investigated. Even though non-coding RNAs represent the vast majority of cellular RNA species, this class of RNA remains historically understudied. This is partly because of a lag in technological and bioinformatic innovations specifically capable of identifying and accurately measuring their expression. Nevertheless, recent progress in this domain has enabled a profusion of publications identifying novel sub-types of ncRNAs and studies directly addressing the function of ncRNAs in human health and disease. Many ncRNAs, including circular and enhancer RNAs, have now been demonstrated to play key functions in the regulation of immune cells and to show associations with immune-mediated diseases. Some ncRNAs may function as biomarkers of disease, aiding in diagnostics and in estimating response to treatment, while others may play a direct role in the pathogenesis of disease. Importantly, some are relatively stable and are amenable to therapeutic targeting, for example through gene therapy. Here, we provide an overview of ncRNAs and review technological advances that enable their study and hold substantial promise for the future. We provide context-specific examples by examining the associations of ncRNAs with four prototypical human autoimmune diseases, specifically rheumatoid arthritis, psoriasis, inflammatory bowel disease and multiple sclerosis. We anticipate that the utility and mechanistic roles of these ncRNAs in autoimmunity will be further elucidated in the near future.

Single-cell RNA sequencing in the context of neuropathic pain: progress, challenges, and prospects

Neuropathic pain, characterized by persistent or intermittent spontaneous pain as well as some unpleasant abnormal sensations, is one of the most prevalent health problems in the world. Ectopic nerve activity, central and peripheral nociceptive sensitization and many other potential mechanisms may participate in neuropathic pain. The complexity and ambiguity of neuropathic pain mechanisms result in difficulties in pain management, and existing treatment plans provide less-than-satisfactory relief. In recent years, single-cell RNA sequencing (scRNA-seq) has been increasingly applied and has become a powerful means for biological researchers to explore the complexity of neurobiology. This technique can be used to perform unbiased, high-throughput and high-resolution transcriptional analyses of neuropathic pain-associated cells, improving the understanding of neuropathic pain mechanisms and enabling individualized pain management. To date, scRNA-seq has been preliminarily used in neuropathic pain research for applications such as compiling a dorsal root ganglion atlas, identifying new cell types and discovering gene regulatory networks associated with neuropathic pain. Although scRNA-seq is a relatively new technique in the neuropathic pain field, there have been several studies based on animal models. However, because of the various differences between animals and humans, more attention should be given to translational medicine research. With the aid of scRNA-seq, researchers can further explore the mechanism of neuropathic pain to improve the clinical understanding of the diagnosis, treatment and management of neuropathic pain.

3’UTR-Seq analysis of chicken abdominal adipose tissue reveals widespread intron retention in 3’UTR and provides insight into molecular basis of feed efficiency

Feed efficiency (FE) is an important trait in the broiler industry due to its direct correlation to efficient muscle growth instead of fat deposition. The present study characterized and compared gene expression profiles in abdominal fat from broiler chickens of different FE levels to enhance the understanding of FE biology. Specifically, traditional whole-transcript RNA-sequencing (RNA-seq) and 3’ UTR-sequencing (3’ UTR-seq) were applied to 22 and 61 samples, respectively. Overall, these two sequencing techniques shared a high correlation (0.76) between normalized counts, although 3’ UTR-seq showed a higher variance in sequencing and mapping performance statistics across samples and a lower rate of uniquely mapped reads. A higher percentage of 3’ UTR-seq reads mapped to introns suggested the frequent presence of cleavage sites in introns, thus warranting future research to study its regulatory function. Differential expression analysis identified 1198 differentially expressed genes (DEGs) between high FE (HFE) and intermediate FE (IFE) chickens with False Discovery Rate < 0.05 and fold change > 1.2. The processes that were significantly enriched by the DEGs included extracellular matrix remodeling and mechanisms impacting gene expression at the transcriptional and translational levels. Gene ontology enrichment analysis suggested that the divergence in fat deposition and FE in broiler chickens could be associated with peroxisome and lipid metabolism possibly regulated by G0/G1 switch gene 2 (G0S2).

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein.

Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences.

Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data.

Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Dynamics and thermal sensitivity of ribosomal RNA maturation paths in plants

Ribosome biogenesis is a constitutive fundamental process for cellular function. Its rate of production depends on the rate of maturation of precursor ribosomal RNA (pre-rRNA). The rRNA maturation paths are marked by four dominant rate-limiting intermediates with cell-type variation of the processivity rate. We have identified that high temperature stress in plants, while halting the existing pre-rRNA maturation schemes, also transiently triggers an atypical pathway for 35S pre-rRNA processing. This pathway leads to production of an aberrant precursor rRNA, reminiscent of yeast 24S, encompassing 18S and 5.8S rRNA that do not normally co-occur together at sub-unit levels; this response is elicited specifically by high and not low temperatures. We show this response to be conserved in two other model crop plant species (Rice and Tomato). This pathway persists even after returning to normal growth conditions for 1 hour and is reset between 1-6 hours after stress treatment, likely, due to resumption of normal 35S pre-rRNA synthesis and processing. The heat-induced ITS2 cleavage-derived precursors and stalled P-A2-like precursors were heterogeneous in nature with a fraction containing polymeric (A) tails. Furthermore, high temperature treatment and subsequent fractionation resulted in polysome and precursor rRNA depletion.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein. Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences. Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data. Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein. Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences. Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data. Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Comparison of urinary extracellular vesicle isolation methods for transcriptomic biomarker research in diabetic kidney disease

Urinary Extracellular Vesicles (uEV) have emerged as a source for biomarkers of kidney damage, holding potential to replace the conventional invasive techniques including kidney biopsy. However, comprehensive studies characterizing uEV isolation methods with patient samples are rare. Here we compared performance of three established uEV isolation workflows for their subsequent use in transcriptomics analysis for biomarker discovery in diabetic kidney disease. We collected urine samples from individuals with type 1 diabetes with macroalbuminuria and healthy controls. We isolated uEV by Hydrostatic Filtration Dialysis (HFD), ultracentrifugation (UC), and a commercial kit- based isolation method (NG), each with different established urine clearing steps. Purified EVs were analysed by electron microscopy, nanoparticle tracking analysis, and Western blotting. Isolated RNAs were subjected to miRNA and RNA sequencing. HFD and UC samples showed close similarities based on mRNA sequencing data. NG samples had a lower number of reads and different mRNA content compared to HFD or UC. For miRNA sequencing data, satisfactory miRNA counts were obtained by all methods, but miRNA contents differed slightly. This suggests that the isolation workflows enrich specific subpopulations of miRNA-rich uEV preparation components. Our data shows that HFD,UC and the kit-based method are suitable methods to isolate uEV for miRNA-seq. However, only HFD and UC were suitable for mRNA-seq in our settings.

SMALL ORGAN4 Is a Ribosome Biogenesis Factor Involved in 5.8S Ribosomal RNA Maturation

Ribosome biogenesis is crucial for cellular metabolism and has important implications for disease and aging. Human (Homo sapiens) glioma tumor-suppressor candidate region gene2 (GLTSCR2) and yeast (Saccharomyces cerevisiae) Nucleolar protein53 (Nop53) are orthologous proteins with demonstrated roles as ribosome biogenesis factors; knockdown of GLTSCR2 impairs maturation of 18S and 5.8S ribosomal RNAs (rRNAs), and Nop53 is required for maturation of 5.8S and 25S rRNAs. Here, we characterized SMALL ORGAN4 (SMO4), the most likely ortholog of human GLTSCR2 and yeast Nop53 in Arabidopsis (Arabidopsis thaliana). Loss of function of SMO4 results in a mild morphological phenotype; however, we found that smo4 mutants exhibit strong cytological and molecular phenotypes: nucleolar hypertrophy and disorganization, overaccumulation of 5.8S and 18S rRNA precursors, and an imbalanced 40S:60S ribosome subunit ratio. Like yeast Nop53 and human GLTSCR2, Arabidopsis SMO4 participates in 5.8S rRNA maturation. In yeast, Nop53 cooperates with mRNA transport4 (Mtr4) for 5.8S rRNA maturation. In Arabidopsis, we found that SMO4 plays similar roles in the 5.8S rRNA maturation pathway than those described for MTR4. However, SMO4 seems not to participate in the degradation of by-products derived from the 5'-external transcribed spacer (ETS) of 45S pre-rRNA, as MTR4 does.

Gapmer Antisense Oligonucleotides Targeting 5S Ribosomal RNA Can Reduce Mature 5S Ribosomal RNA by Two Mechanisms

In this study, we demonstrate that 5S ribosomal RNA (rRNA), a highly structured and protein-bound RNA, is quite difficult to reduce with antisense oligonucleotides (ASOs). However, we found a single accessible site that was targetable with a high-affinity complementary ASO. The ASO appeared to bind to the site, recruit RNaseH1, and cause degradation of the 5S RNA. Intriguingly, we also observed that the same ASO induced an accumulation of pre-5S RNA, which may contribute to reduced levels of mature 5S rRNA. As expected, ASO mediated reduction of 5S RNA, and modest inhibition of processing of pre-5S RNA resulted in nucleolar toxicity. However, the toxicity induced was minimal compared with actinomycin D, consistent with its modest effects on pre-5S rRNA. Mechanistically, we show that the accumulation of pre-5S rRNA required ASO hybridization to the cognate rRNA sequence but was independent of RNaseH1 activity. We found that Ro60 and La, proteins known to bind misprocessed RNAs, likely sequester the ASO-pre-5S rRNA species and block RNaseH1 activity, thus identifying another example of competitive mechanisms mediated by proteins that compete with RNaseH1 for binding to ASO-RNA heteroduplexes.

A tale of non-canonical tails: gene regulation by post-transcriptional RNA tailing

RNA tailing, or the addition of non-templated nucleotides to the 3' end of RNA, is the most frequent and conserved type of RNA modification. The addition of tails and their composition reflect RNA maturation stages and have important roles in determining the fate of the modified RNAs. Apart from canonical poly(A) polymerases, which add poly(A) tails to mRNAs in a transcription-coupled manner, a family of terminal nucleotidyltransferases (TENTs), including terminal uridylyltransferases (TUTs), modify RNAs post-transcriptionally to control RNA stability and activity. The human genome encodes 11 different TENTs with distinct substrate specificity, intracellular localization and tissue distribution. In this Review, we discuss recent advances in our understanding of non-canonical RNA tails, with a focus on the functions of human TENTs, which include uridylation, mixed tailing and post-transcriptional polyadenylation of mRNAs, microRNAs and other types of non-coding RNA.

Biallelic loss of FAM46C triggers tumor growth with concomitant activation of Akt signaling in multiple myeloma cells

Loss of heterozygosity or mutation of the family with sequence similarity 46, member C (FAM46C) gene on chromosome band 1p12 is associated with shorter overall survival of patients with multiple myeloma (MM). In this study, using human MM cell lines (KMS‐11, OCI‐My5, and ANBL‐6), we generated FAM46C^−/− cell clones and examined the effect of disruption of FAM46C on cell survival and cellular signaling. Cell proliferation assays showed increased clonogenicity of FAM46C^−/− KMS‐11 cells compared to WT cells. Xenograft experiments showed significantly shorter overall survival of mice harboring the FAM46C^−/− cell‐derived tumors than mice with the FAM46C^WT cell‐derived tumors. Notably, levels of phosphorylated Akt and its substrates increased both in vitro and in vivo in the FAM46C^−/− cells compared to WT cells. In addition, caspase activities decreased in the FAM46C^−/− cells. Results of gene set enrichment analysis showed that loss of FAM46C significantly activated serum‐responsive genes while inactivating phosphatase and tensin homolog (PTEN)‐related genes. Mechanistically, loss of FAM46C decreased the PTEN activity, number of apoptotic cells, and caspase activities. PF‐04691502, a selective PI3K inhibitor, suppressed the augmented phosphorylation of Akt and its substrate FoxO3a. Treatment with afuresertib (a specific Akt inhibitor) in combination with bortezomib additively decreased FAM46C^−/− MM cell survival. Collectively, this study is the first to report that loss of FAM46C triggers the concomitant activation of the PI3K‐Akt signaling pathway, which might be a therapeutic target for MM with abnormalities in the FAM46C gene.

Consistency of Targeted Metatranscriptomics and Morphological Characterization of Phytoplankton Communities

The composition of phytoplankton community is the basis for environmental monitoring and assessment of the ecological status of aquatic ecosystems. Community composition studies of phytoplankton have been based on time-consuming and expertise-demanding light microscopy analyses. Molecular methods have the potential to replace microscopy, but the high copy number variation of ribosomal genes and the lack of universal primers for simultaneous amplification of prokaryotic and eukaryotic genes complicate data interpretation. In this study, we used our previously developed directional primer-independent high-throughput sequencing (HTS) approach to analyze 16S and 18S rRNA community structures. Comparison of 83 boreal lake samples showed that the relative abundances of eukaryotic phytoplankton at class level and prokaryotic cyanobacteria at order level were consistent between HTS and microscopy results. At the genus level, the results had low correspondence, mainly due to lack of sequences in the reference library. HTS was superior to identify genera that are extensively represented in the reference databases but lack specific morphological characteristics. Targeted metatranscriptomics proved to be a feasible method to complement the microscopy analysis. The metatranscriptomics can also be applied without linking the sequences to taxonomy. However, direct indexing of the sequences to their environmental indicator values needs collections of more comprehensive sample sets, as long as the coverage of molecular barcodes of eukaryotic species remains insufficient.

Single cell transcriptomics: moving towards multi-omics

Single-cell multi-omics analysis helps characterize multiple layers of molecular features at a single-cell scale to provide insights into cellular processes and functions.

The Implication of mRNA Degradation Disorders on Human DISease: Focus on DIS3 and DIS3-Like Enzymes

RNA degradation is considered a critical posttranscriptional regulatory checkpoint, maintaining the correct functioning of organisms. When a specific RNA transcript is no longer required in the cell, it is signaled for degradation through a number of highly regulated steps. Ribonucleases (or simply RNases) are key enzymes involved in the control of RNA stability. These enzymes can perform the RNA degradation alone or cooperate with other proteins in RNA degradation complexes. Important findings over the last years have shed light into eukaryotic RNA degradation by members of the RNase II/RNB family of enzymes. DIS3 enzyme belongs to this family and represents one of the catalytic subunits of the multiprotein complex exosome. This RNase has a diverse range of functions, mainly within nuclear RNA metabolism. Humans encode two other DIS3-like enzymes: DIS3L (DIS3L1) and DIS3L2. DIS3L1 also acts in association with the exosome but is strictly cytoplasmic. In contrast, DIS3L2 acts independently of the exosome and shows a distinctive preference for uridylated RNAs. These enzymes have been shown to be involved in important cellular processes, such as mitotic control, and associated with human disorders like cancer. This review shows how the impairment of function of each of these enzymes is implicated in human disease.

Virulence- and signaling-associated genes display a preference for long 3'UTRs during rice infection and metabolic stress in the rice blast fungus

Generation of mRNA isoforms by alternative polyadenylation (APA) and their involvement in regulation of fungal cellular processes, including virulence, remains elusive. Here, we investigated genome-wide polyadenylation site (PAS) selection in the rice blast fungus to understand how APA regulates pathogenicity. More than half of Magnaporthe oryzae transcripts undergo APA and show novel motifs in their PAS region. Transcripts with shorter 3'UTRs are more stable and abundant in polysomal fractions, suggesting they are being translated more efficiently. Importantly, rice colonization increases the use of distal PASs of pathogenicity genes, especially those participating in signalling pathways like 14-3-3B, whose long 3'UTR is required for infection. Cleavage factor I (CFI) Rbp35 regulates expression and distal PAS selection of virulence and signalling-associated genes, tRNAs and transposable elements, pointing its potential to drive genomic rearrangements and pathogen evolution. We propose a noncanonical PAS selection mechanism for Rbp35 that recognizes UGUAH, unlike humans, without CFI25. Our results showed that APA controls turnover and translation of transcripts involved in fungal growth and environmental adaptation. Furthermore, these data provide useful information for enhancing genome annotations and for cross-species comparisons of PASs and PAS usage within the fungal kingdom and the tree of life.

Nucleotide Modifications Decrease Innate Immune Response Induced by Synthetic Analogs of snRNAs and snoRNAs

Short nuclear regulatory RNAs play a key role in the main stages of maturation of the precursors of the major RNA species. Small nuclear RNAs (snRNAs) form the core of the spliceosome and are responsible for the splicing of pre-mRNA molecules. Small nucleolar RNAs (snoRNAs) direct post-transcriptional modification of pre-rRNAs. A promising strategy for the development of non-coding RNA (ncRNAs) mimicking molecules is the introduction of modified nucleotides, which are normally present in natural ncRNAs, into the structure of synthetic RNAs. We have created a set of snoRNAs and snRNA analogs and studied the effect of base modifications, specifically, pseudouridine (Ψ) and 5-methylcytidine (m⁵C), on the immune-stimulating and cytotoxic properties of these RNAs. Here, we performed a whole-transcriptome study of the influence of synthetic snoRNA analogs with various modifications on gene expression in human cells. Moreover, we confirmed the role of PKR in the recognition of snoRNA and snRNA analogs using the short hairpin RNA (shRNA) technique. We believe that the data obtained will contribute to the understanding of the role of nucleotide modification in ncRNA functions, and can be useful for creating the agents for gene regulation based on the structure of natural snoRNAs and snRNAs.

Nucleolar Dominance in a Tetraploidy Hybrid Lineage Derived From Carassius auratus red var. () × Megalobrama amblycephala ()

Nucleolar dominance is related to the expression of 45S rRNA genes inherited from one progenitor due to the silencing of the other progenitor’s rRNA genes. To investigate nucleolar dominance associated with tetraploidization, we analyzed the changes regarding the genetic traits and expression of 45S rRNA genes in tetraploidy hybrid lineage including F₁ allotetraploids (4n = 148) and F₂ autotetraploids (4n = 200) derived from the distant hybridization of Carassius auratus red var. (2n = 100) () ×Megalobrama amblycephala (2n = 48) (). Results showed that nucleolar dominance from the females was established in F₁ hybrids and it was inherited in F₂ hybrids, suggesting that tetraploidization can lead to rapid establishment of nucleolar dominance in the hybrid origin’s tetraploid lineage. These results extend the knowledge of nucleolar dominance in polyploidy hybrid animals, which are of significance for the evolution of hybrids in vertebrates.

Tales around the clock: Poly(A) tails in circadian gene expression

Circadian rhythms are ubiquitous time-keeping processes in eukaryotes with a period of ~24 hr. Light is perhaps the main environmental cue (zeitgeber) that affects several aspects of physiology and behaviour, such as sleep/wake cycles, orientation of birds and bees, and leaf movements in plants. Temperature can serve as the main zeitgeber in the absence of light cycles, even though it does not lead to rhythmicity through the same mechanism as light. Additional cues include feeding patterns, humidity, and social rhythms. At the molecular level, a master oscillator orchestrates circadian rhythms and organizes molecular clocks located in most cells. The generation of the 24 hr molecular clock is based on transcriptional regulation, as it drives intrinsic rhythmic changes based on interlocked transcription/translation feedback loops that synchronize expression of genes. Thus, processes and factors that determine rhythmic gene expression are important to understand circadian rhythms. Among these, the poly(A) tails of RNAs play key roles in their stability, translational efficiency and degradation. In this article, we summarize current knowledge and discuss perspectives on the role and significance of poly(A) tails and associating factors in the context of the circadian clock. This article is categorized under: RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > 3' End Processing.

Transcriptome-wide discovery of coding and noncoding RNA-binding proteins

Transcriptome-wide identification of RNA-binding proteins (RBPs) is a prerequisite for understanding the posttranscriptional gene regulation networks. However, proteomic profiling of RBPs has been mostly limited to polyadenylated mRNA-binding proteins, leaving RBPs on nonpoly(A) RNAs, including most noncoding RNAs (ncRNAs) and pre-mRNAs, largely undiscovered. Here we present a click chemistry-assisted RNA interactome capture (CARIC) strategy, which enables unbiased identification of RBPs, independent of the polyadenylation state of RNAs. CARIC combines metabolic labeling of RNAs with an alkynyl uridine analog and in vivo RNA-protein photocross-linking, followed by click reaction with azide-biotin, affinity enrichment, and proteomic analysis. Applying CARIC, we identified 597 RBPs in HeLa cells, including 130 previously unknown RBPs. These newly discovered RBPs can likely bind ncRNAs, thus uncovering potential involvement of ncRNAs in processes previously unknown to be ncRNA-related, such as proteasome function and intermediary metabolism. The CARIC strategy should be broadly applicable across various organisms to complete the census of RBPs.

Single-cell RNA sequencing: Technical advancements and biological applications

Cells are the basic building blocks of organisms and each cell is unique. Single-cell RNA sequencing has emerged as an indispensable tool to dissect the cellular heterogeneity and decompose tissues into cell types and/or cell states, which offers enormous potential for de novo discovery. Single-cell transcriptomic atlases provide unprecedented resolution to reveal complex cellular events and deepen our understanding of biological systems. In this review, we summarize and compare single-cell RNA sequencing technologies, that were developed since 2009, to facilitate a well-informed choice of method. The applications of these methods in different biological contexts are also discussed. We anticipate an ever-increasing role of single-cell RNA sequencing in biology with further improvement in providing spatial information and coupling to other cellular modalities. In the future, such biological findings will greatly benefit medical research.

Diurnal Oscillations in Liver Mass and Cell Size Accompany Ribosome Assembly Cycles

The liver plays a pivotal role in metabolism and xenobiotic detoxification, processes that must be particularly efficient when animals are active and feed. A major question is how the liver adapts to these diurnal changes in physiology. Here, we show that, in mice, liver mass, hepatocyte size, and protein levels follow a daily rhythm, whose amplitude depends on both feeding-fasting and light-dark cycles. Correlative evidence suggests that the daily oscillation in global protein accumulation depends on a similar fluctuation in ribosome number. Whereas rRNA genes are transcribed at similar rates throughout the day, some newly synthesized rRNAs are polyadenylated and degraded in the nucleus in a robustly diurnal fashion with a phase opposite to that of ribosomal protein synthesis. Based on studies with cultured fibroblasts, we propose that rRNAs not packaged into complete ribosomal subunits are polyadenylated by the poly(A) polymerase PAPD5 and degraded by the nuclear exosome.

Recent advances in dynamic m6A RNA modification

The identification of m⁶A demethylases and high-throughput sequencing analysis of methylated transcriptome corroborated m⁶A RNA epigenetic modification as a dynamic regulation process, and reignited its investigation in the past few years. Many basic concepts of cytogenetics have been revolutionized by the growing understanding of the fundamental role of m⁶A in RNA splicing, degradation and translation. In this review, we summarize typical features of methylated transcriptome in mammals, and highlight the ‘writers’, ‘erasers’ and ‘readers’ of m⁶A RNA modification. Moreover, we emphasize recent advances of biological functions of m⁶A and conceive the possible roles of m⁶A in the regulation of immune response and related diseases.

Mitochondrial poly(A) polymerase is involved in tRNA repair

Transcription of the mitochondrial genome results in polycistronic precursors, which are processed mainly by the release of tRNAs interspersed between rRNAs and mRNAs. In many metazoan mitochondrial genomes some tRNA genes overlap with downstream genes; in the case of human mitochondria the genes for tRNA^Tyr and tRNA^Cys overlap by one nucleotide. It has previously been shown that processing of the common precursor releases an incomplete tRNA^Tyr lacking the 3′-adenosine. The 3′-terminal adenosine has to be added before addition of the CCA end and subsequent aminoacylation. We show that the mitochondrial poly(A) polymerase (mtPAP) is responsible for this A addition. In vitro, a tRNA^Tyr lacking the discriminator is a substrate for mtPAP. In vivo, an altered mtPAP protein level affected tRNA^Tyr maturation, as shown by sequencing the 3′ ends of mitochondrial tRNAs. Complete repair could be reconstituted in vitro with three enzymes: mtPAP frequently added more than one A to the 3′ end of the truncated tRNA, and either the mitochondrial deadenylase PDE12 or the endonuclease RNase Z trimmed the oligo(A) tail to a single A before CCA addition. An enzyme machinery that evolved primarily for other purposes thus allows to tolerate the frequent evolutionary occurrence of gene overlaps.

Microarray experiments and factors which affect their reliability

Oligonucleotide microarrays belong to the basic tools of molecular biology and allow for simultaneous assessment of the expression level of thousands of genes. Analysis of microarray data is however very complex, requiring sophisticated methods to control for various factors that are inherent to the procedures used. In this article we describe the individual steps of a microarray experiment, highlighting important elements and factors that may affect the processes involved and that influence the interpretation of the results. Additionally, we describe methods that can be used to estimate the influence of these factors, and to control the way in which they affect the expression estimates. A comprehensive understanding of the experimental protocol used in a microarray experiment aids the interpretation of the obtained results. By describing known factors which affect expression estimates this article provides guidelines for appropriate quality control and pre-processing of the data, additionally applicable to other transcriptome analysis methods that utilize similar sample handling protocols.

Mitochondrial gene order variation in the brachiopod Lingula anatina and its implications for mitochondrial evolution in lophotrochozoans

Vertebrate mitochondrial (mt) genomes display highly conserved gene order and relatively low evolutionary rates. However, these features are variable in marine invertebrates. Here we present the mt genome of the lingulid brachiopod, Lingula anatina, from Amami Island, Japan, as part of the nuclear genome project. We obtain ~2000-fold coverage of the 17.9-kb mt genome using Illumina sequencing, and we identify hypervariable regions within the same individual. Transcriptome analyses show that mt transcripts are polycistronic and expressed differentially. Unexpectedly, we find that the mt gene order of Amami Lingula is completely shuffled compared to that of a specimen from Yanagawa, suggesting that there may be cryptic species. Using breakpoint distance analyses with 101 metazoan mt genomes, we show that the evolutionary history of mt gene order among lophotrochozoans is unique. Analyses of non-synonymous substitution rates reveal that mt protein-coding genes of Lingula have experienced rapid evolution comparable to that expected for interspecific comparisons. Whole genome phylogenetic analyses suggest that mt genomes have limited value for inferring the phylogenetic positions of lophotrochozoans because of their high evolutionary rates in brachiopods and bivalves.

Single-cell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos

Circular RNAs (circRNAs) are a new class of non-polyadenylated non-coding RNAs that may play important roles in many biological processes. Here we develop a single-cell universal poly(A)-independent RNA sequencing (SUPeR-seq) method to sequence both polyadenylated and non-polyadenylated RNAs from individual cells. This method exhibits robust sensitivity, precision and accuracy. We discover 2891 circRNAs and 913 novel linear transcripts in mouse preimplantation embryos and further analyze the abundance of circRNAs along development, the function of enriched genes, and sequence features of circRNAs. Our work is key to deciphering regulation mechanisms of circRNAs during mammalian early embryonic development.

Adapted approach to profile genes while reconciling Vegf-a mRNA expression in the developing and injured lung

During lung development and injury, messenger RNA (mRNA) transcript levels of genes fluctuate over both space and time. Quantitative PCR (qPCR) is a highly sensitive, widely used technique to measure the mRNA levels. The sensitivity of this technique can be disadvantageous and errors amplified when each qPCR assay is not validated. In contrast to other organs, lungs have high RNase activity, resulting in less than optimal RNA integrity. We implemented a strategy to address these limitations in developing and injured lungs. Parameters were established and a filter designed that optimized amplicon length and included or excluded samples based on RNA integrity. This approach was illustrated and validated by measuring mRNA levels including Vegf-a in newborn mouse lungs that were injured by 85% oxygen (hyperoxia) for 12 days and compared with control (normoxia). We demonstrate that, in contrast to contradictory Vegf-a expression when normalized to the least suitable housekeeping genes, application of this filter and normalization to most suitable three housekeeping genes, Hprt, Eef2, and Rpl13a, gave reproducible Vegf-a expression, thus corroborating the sample filter. Accordingly, both short amplicon length and proper normalization to ranked, evaluated genes minimized erroneous fluctuation and qPCR amplification issues associated with nonideal RNA integrity in injured and developing lungs. Furthermore, our work uncovers how RNA integrity, purity, amplicon length, and discovery of stable candidate reference genes enhance precision of qPCR results and utilizes the advantages of qPCR in developmental studies.

Ribonucleic acid (RNA) biosynthesis in human cancer

In many respects, the most remarkable chemical substances within the genome of eukaryotic cells are remarkable proteins which are the critical structural and functional units of living cells. The specifications for everything that goes in the cell are natural digital-to-digital decoding process in an archive sequence by deoxyribonucleic acid (DNA) and an articulate construction by ribonucleic acid (RNA). The products of DNA transcription are long polymers of ribonucleotides rather than deoxyribonucleotides and are termed ribonucleic acids. Certain deoxyribonucleotide sequences, or genes, give rise to transfer RNA (tRNA) and other ribosomal RNA (rRNA) when transcribed. The ribonucleotide sequences fold extensively and rRNA is associated with specific proteins to yield the essential cell components, ribosomes. Transcription of other special sequences yields messenger RNAs (mRNAs) that contain ribonucleotide sequences that will be ultimately translated into new types of amino acid sequences of functional cellular protein molecules. This switch to a different variety of cellular molecular sequences is complex, but each sequence of the three ribonucleotides specifies the insertion of one particular amino acid into the polypeptide chain under production. Whilst mRNA is considered the vehicle by which genetic information is transmitted from the genome and allocated in the appropriate cytoplasmic sites for translation into protein via cap-dependent mechanism, the actual translation depends also on the presence of other so-called household and luxury protein molecules. Recent evidence suggests RNA species are required at initiation, because treatment of cells with antibiotics or drugs that inhibit RNA synthesis cause a decrease in protein synthesis. The rRNA is necessary as a structural constituent of the ribosomes upon which translation takes place, whereas tRNA is necessary as an adaptor in amino acid activation and elongation protein chains to ribosomes. In this article, we review malignant tumor, with stem like properties, and recent technical advances into the phenomenon of micro-particles and micro-vesicles containing cell-free nucleic acids that circulate plasma. New areas of research have been opened into screening tumor telomerase progression, prognosis of aptamers targeting cell surface, monitoring the efficacy of anticancer therapies, oncogenic transformation of host cell, and RNA polymerases role in the cell cycle progression and differentiation.

Identification of a recently active mammalian SINE derived from ribosomal RNA

Complex eukaryotic genomes are riddled with repeated sequences whose derivation does not coincide with phylogenetic history and thus is often unknown. Among such sequences, the capacity for transcriptional activity coupled with the adaptive use of reverse transcription can lead to a diverse group of genomic elements across taxa, otherwise known as selfish elements or mobile elements. Short interspersed nuclear elements (SINEs) are nonautonomous mobile elements found in eukaryotic genomes, typically derived from cellular RNAs such as tRNAs, 7SL or 5S rRNA. Here, we identify and characterize a previously unknown SINE derived from the 3'-end of the large ribosomal subunit (LSU or 28S rDNA) and transcribed via RNA polymerase III. This new element, SINE28, is represented in low-copy numbers in the human reference genome assembly, wherein we have identified 27 discrete loci. Phylogenetic analysis indicates these elements have been transpositionally active within primate lineages as recently as 6 MYA while modern humans still carry transcriptionally active copies. Moreover, we have identified SINE28s in all currently available assembled mammalian genome sequences. Phylogenetic comparisons indicate that these elements are frequently rederived from the highly conserved LSU rRNA sequences in a lineage-specific manner. We propose that this element has not been previously recognized as a SINE given its high identity to the canonical LSU, and that SINE28 likely represents one of possibly many unidentified, active transposable elements within mammalian genomes.

The dynamic epitranscriptome: N6-methyladenosine and gene expression control

N(6)-methyladenosine (m(6)A) is a modified base that has long been known to be present in non-coding RNAs, ribosomal RNA, polyadenylated RNA and at least one mammalian mRNA. However, our understanding of the prevalence of this modification has been fundamentally redefined by transcriptome-wide m(6)A mapping studies, which have shown that m(6)A is present in a large subset of the transcriptome in specific regions of mRNA. This suggests that mRNA may undergo post-transcriptional methylation to regulate its fate and function, which is analogous to methyl modifications in DNA. Thus, the pattern of methylation constitutes an mRNA 'epitranscriptome'. The identification of adenosine methyltransferases ('writers'), m(6)A demethylating enzymes ('erasers') and m(6)A-binding proteins ('readers') is helping to define cellular pathways for the post-transcriptional regulation of mRNAs.

Presence of poly(A) and poly(A)-rich tails in a positive-strand RNA virus known to lack 3׳ poly(A) tails

Here we show that Tobacco mosaic virus (TMV), a positive-strand RNA virus known to end with 3׳ tRNA-like structures, does possess a small fraction of gRNA bearing polyadenylate tails. Particularly, many tails are at sites corresponding to the 3׳ end of near full length gRNA, and are composed of poly(A)-rich sequences containing the other nucleotides in addition to adenosine, resembling the degradation-stimulating poly(A) tails observed in all biological kingdoms. Further investigations demonstrate that these polyadenylated RNA species are not enriched in chloroplasts. Silencing of cpPNPase, a chloroplast-localized polynucleotide polymerase known to not only polymerize the poly(A)-rich tails but act as a 3׳ to 5׳ exoribonuclease, does not change the profile of polyadenylate tails associated with TMV RNA. Nevertheless, because similar tails were also detected in other phylogenetically distinct positive-strand RNA viruses lacking poly(A) tails, such kind of polyadenylation may reflect a common but as-yet-unknown interface between hosts and viruses.

Circularized RT-PCR (cRT-PCR): analysis of the 5' ends, 3' ends, and poly(A) tails of RNA

Many techniques that are used to characterize individual RNA molecules can potentially alter the original transcript sequence or its posttranscriptional modifications, such as polyadenylation. Methods that are designed to define the ends of an RNA molecule, for example, oligonucleotide ligation, avoid altering the transcript sequence but can usually fulfill only one objective per experiment (e.g., define the 5' or the 3' end). In contrast, not only does circularized reverse transcription coupled with PCR (cRT-PCR) preserve the original 5' and 3' ends of the transcript and posttranscriptionally added extensions, but also the material from one experimental procedure can be utilized in order to characterize both the 5' and 3' ends. Furthermore, if suitable oligonucleotide primers are designed, cRT-PCR can be used to isolate truncated, adenylated (and nonadenylated) molecules that are intermediates in RNA decay (Slomovic and Schuster, 2008).

Oligo(dT)-primed RT-PCR isolation of polyadenylated RNA degradation intermediates

The posttranscriptional modification of RNA by polyadenylation serves various purposes, among them to assist in RNA degradation (see an alternative protocol for measuring RNA degradation on Method for measuring mRNA decay rate in Saccharomyces cerevisiae). This function, once thought to occur in prokaryotic or organellar systems alone, is now known to operate in the nuclei and cytoplasm of eukaryotes as well (Slomovic et al., 2008; Slomovic et al., 2010; Houseley and Tollervey, 2009; Deutscher, 2006). Poly(A)-assisted RNA decay begins with the endonucleolytic cleavage of the transcript. Following this, a poly(A) or oligo(A) tail is added to the 3' end of the cleavage product. This tag serves as a 'landing pad' for 3'-5' exoribonucleases that then begin to digest the RNA fragment. Truncated RNA molecules that have undergone tail addition but have yet to be degraded are called degradation intermediates. The detection of such intermediates is considered a tell-tale sign that poly(A)-assisted RNA decay occurs in the organism being studied. Determination of the tail nucleotide composition by DNA sequencing often aids the researcher in identifying the enzyme responsible for tail synthesis since tails can be either homopolymeric (exclusively A residues) or heteropolymeric (A-rich tails that may include other nucleotides). The following protocol, based on oligo(dT)-primed reverse transcription, describes the step-by-step detection and isolation of adenylated degradation intermediates in the study of poly(A)-assisted RNA decay.

Co-expressed mitochondrial genomes: recently masculinized, recombinant mitochondrial genome is co-expressed with the female-transmitted mtDNA genome in a male Mytilus trossulus mussel from the Baltic Sea

BACKGROUND

Few exceptions have been described from strict maternal inheritance of mitochondrial DNA in animals, including sea mussels (Mytilidae), clams (Donacidae, Veneridae and Solenidae) and freshwater mussels (Unionoidae) order. In these bivalves mitochondria and their DNA are transferred through two separate routes. The females inherit only the maternal mitochondrial DNA whereas the males inherit maternal as well as paternal mitochondrial DNA, which is usually present only in gonads and sperm. The mechanism controlling this phenomenon is unclear but leads to the existence of two separate mitochondrial DNA lineages in a single species. The lineages are usually well differentiated: up to 20-50% divergence in nucleotide sequence. Occasionally, a maternal mitochondrial DNA can invade the paternal transmission route, eventually replacing the diverged M-type and lowering the divergence. Such role reversal (masculinization) event has happened recently in the Mytilus population of the Baltic Sea which consists of M. edulis × M. trossulus hybrids, but the functional status of the resulting mitochondrial genome was unknown.

RESULTS

In this paper we sequenced transcripts from one specimen that was identified as male carrying both the female mitochondrial genome and a recently masculinized mitochondrial genome. Additionally, the analysis of the control region has showed that the recently masculinized, recombinant genome, not only has an M-type control region and all coding regions derived from the F-type, but also is transcriptionally active along side the maternally inherited F-type genome. In the comparative analysis, the two genomes exhibit different substitution patterns, typical for the M vs. F genome comparisons. The genetic distances and ratios of non-synonymous substitutions also suggest that one of the genomes is transitioning from the maternal to the paternal inheritance mode, consistent with its recent masculinization.

CONCLUSION

We have shown, for the first time, that the recently masculinized mitochondrial genome is active and that it accumulates excess of non-synonymous substitutions across its coding sequence. This suggests, that, under certain cytonuclear incompatibility conditions, masculinization may serve to restore the endangered functionality of the paternally inherited genome. This is also another example of a mitochondrial genome in which the recombination in the control region predated its transition from paternal to maternal transmission route.

Polyadenylation in bacteria and organelles

Polyadenylation is a posttranscriptional modification present throughout all the kingdoms of life with important roles in regulation of RNA stability, translation, and quality control. Functions of polyadenylation in prokaryotic and organellar RNA metabolism are still not fully characterized, and poly(A) tails appear to play contrasting roles in different systems. Here we present a general overview of the polyadenylation process and the factors involved in its regulation, with an emphasis on the diverse functions of 3' end modification in the control of gene expression in different biological systems.

RNA polyadenylation sites on the genomes of microorganisms, animals, and plants

Pre–messenger RNA (mRNA) 3′-end cleavage and subsequent polyadenylation strongly regulate gene expression. In comparison with the upstream or downstream motifs, relatively little is known about the feature differences of polyadenylation [poly(A)] sites among major kingdoms. We suspect that the precise poly(A) sites are very selective, and we therefore mapped mRNA poly(A) sites on complete and nearly complete genomes using mRNA sequences available in the National Center for Biotechnology Information (NCBI) Nucleotide database. In this paper, we describe the mRNA nucleotide [i.e., the poly(A) tail attachment position] that is directly in attachment with the poly(A) tail and the pre-mRNA nucleotide [i.e., the poly(A) tail starting position] that corresponds to the first adenosine of the poly(A) tail in the 29 most-mapped species (2 fungi, 2 protists, 18 animals, and 7 plants). The most representative pre-mRNA dinucleotides covering these two positions were UA, CA, and GA in 17, 10, and 2 of the species, respectively. The pre-mRNA nucleotide at the poly(A) tail starting position was typically an adenosine [i.e., A-type poly(A) sites], sometimes a uridine, and occasionally a cytidine or guanosine. The order was U>C>G at the attachment position but A>>U>C≥G at the starting position. However, in comparison with the mRNA nucleotide composition (base composition), the poly(A) tail attachment position selected C over U in plants and both C and G over U in animals, in both A-type and non-A-type poly(A) sites. Animals, dicot plants, and monocot plants had clear differences in C/G ratios at the poly(A) tail attachment position of the non-A-type poly(A) sites. This study of poly(A) site evolution indicated that the two positions within poly(A) sites had distinct nucleotide compositions and were different among kingdoms.

The human nuclear poly(a)-binding protein promotes RNA hyperadenylation and decay

Control of nuclear RNA stability is essential for proper gene expression, but the mechanisms governing RNA degradation in mammalian nuclei are poorly defined. In this study, we uncover a mammalian RNA decay pathway that depends on the nuclear poly(A)-binding protein (PABPN1), the poly(A) polymerases (PAPs), PAPα and PAPγ, and the exosome subunits RRP6 and DIS3. Using a targeted knockdown approach and nuclear RNA reporters, we show that PABPN1 and PAPα, redundantly with PAPγ, generate hyperadenylated decay substrates that are recognized by the exosome and degraded. Poly(A) tail extension appears to be necessary for decay, as cordycepin treatment or point mutations in the PAP-stimulating domain of PABPN1 leads to the accumulation of stable transcripts with shorter poly(A) tails than controls. Mechanistically, these data suggest that PABPN1-dependent promotion of PAP activity can stimulate nuclear RNA decay. Importantly, efficiently exported RNAs are unaffected by this decay pathway, supporting an mRNA quality control function for this pathway. Finally, analyses of both bulk poly(A) tails and specific endogenous transcripts reveals that a subset of nuclear RNAs are hyperadenylated in a PABPN1-dependent fashion, and this hyperadenylation can be either uncoupled or coupled with decay. Our results highlight a complex relationship between PABPN1, PAPα/γ, and nuclear RNA decay, and we suggest that these activities may play broader roles in the regulation of human gene expression.

In eukaryotes, mRNAs include a stretch of adenosine nucleotides at their 3′ end termed the poly(A) tail. In the cytoplasm, the poly(A) tail stimulates translation of the mRNA into protein, and protects the transcript from degradation. Evidence suggests that poly(A) tails may play distinct roles in RNA metabolism in the nucleus, but little is known about these functions and mechanisms. We show here that poly(A) tails can stimulate transcript decay in the nucleus, a function mediated by the ubiquitous nuclear poly(A) binding protein PABPN1. We find that PABPN1 is required for the degradation of a viral nuclear noncoding RNA as well as an inefficiently exported human mRNA. Importantly, the targeting of RNAs to this decay pathway requires the PABPN1 and poly(A) polymerase-dependent extension of the poly(A) tail. Nuclear transcripts with longer poly(A) tails are then selectively degraded by components of the nuclear exosome. These studies elucidate mechanisms that mammalian cells use to ensure proper mRNA “quality control” and may be important to regulate the expression of nuclear noncoding RNAs. Furthermore, our results suggest that the poly(A) tail has diverse and context-specific roles in gene expression.

Polyadenylation of 18S rRNA in algae(1)

Polyadenylation is best known for occurring to mRNA of eukaryotes transcribed by RNA polymerase II to stabilize mRNA molecules and promote their translation. rRNAs transcribed by RNA polymerase I or III are typically believed not to be polyadenylated. However, there is increasing evidence that polyadenylation occurs to nucleus-encoded rRNAs as part of the RNA degradation pathway. To examine whether the same polyadenylation-assisted degradation pathway occurs in algae, we surveyed representative species of algae including diatoms, chlorophytes, dinoflagellates and pelagophytes using oligo (dT)-primed reversed transcription PCR (RT-PCR). In all the algal species examined, truncated 18S rRNA or its precursor molecules with homo- or hetero-polymeric poly(A) tails were detected. Mining existing algal expressed sequence tag (EST) data revealed polyadenylated truncated 18S rRNA in four additional phyla of algae. rRNA polyadenylation occurred at various internal positions along the 18S rRNA and its precursor sequences. Moreover, putative homologs of noncanonical poly(A) polymerase (ncPAP) Trf4p, which is responsible for polyadenylating nuclear-encoded RNA and targeting it for degradation, were detected from the genomes and transcriptomes of five phyla of algae. Our results suggest that polyadenylation-assisted RNA degradation mechanism widely exists in algae, particularly for the nucleus-encoded rRNA and its precursors.

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5'ETS molecule using three distinct methods and located the acceptor site between two known 5'ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5'ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5'ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5'ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

Extensive alternative polyadenylation during zebrafish development

The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3' untranslated regions (3' UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-seq reads substantially increased and improved existing 3' UTR annotations, resulting in confidently identified 3' UTRs for >79% of the annotated protein-coding genes in zebrafish. mRNAs from most zebrafish genes undergo alternative CPA, with those from more than a thousand genes using different dominant 3' UTRs at different stages. These included one of the poly(A) polymerase genes, for which alternative CPA reinforces its repression in the ovary. 3' UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3' UTRs are highly expressed in the ovary, yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3' UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At 2 h post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3' UTRs provide a resource for studying gene regulation during vertebrate development.

Extensive alternative polyadenylation during zebrafish development

The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3′ untranslated regions (3′ UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-seq reads substantially increased and improved existing 3′ UTR annotations, resulting in confidently identified 3′ UTRs for >79% of the annotated protein-coding genes in zebrafish. mRNAs from most zebrafish genes undergo alternative CPA, with those from more than a thousand genes using different dominant 3′ UTRs at different stages. These included one of the poly(A) polymerase genes, for which alternative CPA reinforces its repression in the ovary. 3′ UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3′ UTRs are highly expressed in the ovary, yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3′ UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At 2 h post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3′ UTRs provide a resource for studying gene regulation during vertebrate development.

Gut pathology and responses to the microsporidium Nosema ceranae in the honey bee Apis mellifera

The microsporidium Nosema ceranae is a newly prevalent parasite of the European honey bee (Apis mellifera). Although this parasite is presently spreading across the world into its novel host, the mechanisms by it which affects the bees and how bees respond are not well understood. We therefore performed an extensive characterization of the parasite effects at the molecular level by using genetic and biochemical tools. The transcriptome modifications at the midgut level were characterized seven days post-infection with tiling microarrays. Then we tested the bee midgut response to infection by measuring activity of antioxidant and detoxification enzymes (superoxide dismutases, glutathione peroxidases, glutathione reductase, and glutathione-S-transferase). At the gene-expression level, the bee midgut responded to N. ceranae infection by an increase in oxidative stress concurrent with the generation of antioxidant enzymes, defense and protective response specifically observed in the gut of mammals and insects. However, at the enzymatic level, the protective response was not confirmed, with only glutathione-S-transferase exhibiting a higher activity in infected bees. The oxidative stress was associated with a higher transcription of sugar transporter in the gut. Finally, a dramatic effect of the microsporidia infection was the inhibition of genes involved in the homeostasis and renewal of intestinal tissues (Wnt signaling pathway), a phenomenon that was confirmed at the histological level. This tissue degeneration and prevention of gut epithelium renewal may explain early bee death. In conclusion, our integrated approach not only gives new insights into the pathological effects of N. ceranae and the bee gut response, but also demonstrate that the honey bee gut is an interesting model system for studying host defense responses.

New kid on the ID block: neural functions of the Nab2/ZC3H14 class of Cys₃His tandem zinc-finger polyadenosine RNA binding proteins

Polyadenosine RNA binding proteins (Pabs) play critical roles in regulating the polyadenylation, nuclear export, stability, and translation of cellular RNAs. Although most Pabs are ubiquitously expressed and are thought to play general roles in post-transcriptional regulation, mutations in genes encoding these factors have been linked to tissue-specific diseases including muscular dystrophy and now intellectual disability (ID). Our recent work defined this connection to ID, as we showed that mutations in the gene encoding the ubiquitously expressed Cys3His tandem zinc-finger (ZnF) Pab, ZC3H14 (Zinc finger protein, CCCH-type, number 14) are associated with non-syndromic autosomal recessive intellectual disability (NS-ARID). This study provided a first link between defects in Pab function and a brain disorder, suggesting that ZC3H14 plays a required role in regulating RNAs in nervous system cells. Here we highlight key questions raised by our study of ZC3H14 and its ortholog in the fruit fly Drosophila melanogaster, dNab2, and comment on future approaches that could provide insights into the cellular and molecular roles of this class of zinc finger-containing Pabs. We propose a summary model depicting how ZC3H14-type Pabs might play particularly important roles in neuronal RNA metabolism.

Identification of a tubulin-α gene specifically expressed in testis and adductor muscle during stable reference gene selection in the hermaphrodite gonad of the lion's paw scallop Nodipecten subnodosus

For non-model species, as many used for aquaculture, with minimal or no genomic information, relative quantification of gene expression studies requires preliminary research including the isolation of potential reference genes and the identification of those stably expressed under the biological conditions of interest. Here we report on the isolation of five partial gene sequences from gonad tissue cDNA in the functional hermaphrodite scallop Nodipecten subnodosus to be evaluated as reference genes: 18S-rRNA, riboprotein l8 (rp-l8), actin-β (act-β), elongation factor 1α (ef-1α) and alpha-tubulin-α (tub-α). We found that 18S-rRNA was stably expressed independently of the priming method used to reverse transcribe RNA to cDNA, oligo-dT or random hexamer. Stability analysis for the five putative reference genes with geNorm and NormFinder indicated that 18S together with rp-l8 were the most stable genes for normalization of gene expression during gonad development in both, male and female sexual regions of the hermaphrodite N. subnodosus. The least stable gene was tub-α, showing a biased expression profile between sexual regions of the gonad, therefore this gene was analyzed thereafter as a target gene together with vitellogenin (vit) and a DEAD-box RNA helicase (dbx) gene. Relative expression, estimated by normalization with the combination of 18S and rp-l8 as reference genes, indicated that as gonad development advanced two of the target genes were up-regulated, tub-α in the male region and vit in the female region. Whereas an increased expression was expected during development for vit for its known role in vitellogenesis, the increased expression of tub-α in the male sexual region was unexpected, and pointed toward this gene being a testis-specific α-tubulin isotype. Further analyses of gene expression among tissues indicated that tub-α is specifically and highly expressed in the male gonad, although expression in adductor muscle was also observed at significantly lower levels. The existence of testis specific α- and β-tubulins has been previously reported in other taxa, relating their function to sperm axoneme formation. Tissue-specific tubulin genes, particularly their promoters, have recently found an application as native promoters for transgene tissue-specific expression in research and reproductive control of insect plagues. The third target gene, a putative member of the DEAD-box RNA helicase family (dbx), showed no changes in expression during gonad development or between sexual regions, therefore it was chosen to discuss the different statistical inferences resulting from the arbitrary use of 'randomly chosen' reference genes when normalizing gene expression.