RNA Populations Research Articles

RNA-Seq Experiment and Data Analysis.

With the ability to obtain several millions of reads per sample, high-throughput RNA sequencing (RNA-Seq) enables investigation of any transcriptome at a fine resolution. Not just the messenger RNA (mRNA), but a wide variety of different RNA populations (e.g., total RNA, microRNA, long ncRNA, pre-mRNA) can also be investigated using RNA-Seq. While facilitating accurate quantification of gene expression, RNA-Seq offers the opportunity to estimate abundance of isoforms and find novel transcripts and allele-specific transcripts. In this chapter, we describe a protocol to construct an RNA-Seq library for sequencing on Illumina NGS platforms and a computational pipeline to perform RNA-Seq data analysis. The protocols described in this chapter can be applied to the analysis of differential gene expression in control versus 17β-estradiol treatment of in vivo or in vitro systems.

Extending and Running the Mosquito Small RNA Genomics Resource Pipeline.

The Mosquito Small RNA Genomics (MSRG) resource is a repository of analyses on the small RNA transcriptomes of mosquito cell cultures and somatic and gonadal tissues. This resource allows for comparing the regulation dynamics of small RNAs generated from transposons and viruses across mosquito species. This chapter covers the procedures to set up the MSRG resource pipeline as a new installation by detailing the necessary collection of genome reference and annotation files and lists of microRNAs (miRNAs) hairpin sequences, transposon repeats consensus sequences, and virus genome sequences. Proper execution of the MSRG resource pipeline yields outputs amenable to biologists to further analyze with desktop and spreadsheet software to gain insights into the balance between arthropod endogenous small RNA populations and the proportions of virus-derived small RNAs that include Piwi-interacting RNAs (piRNAs) and endogenous small interfering RNAs (siRNAs).

Analysis of the circRNA and T-UCR populations identifies convergent pathways in mouse and human models of Rett syndrome

Noncoding RNAs play regulatory roles in physiopathology, but their involvement in neurodevelopmental diseases is poorly understood. Rett syndrome is a severe, progressive neurodevelopmental disorder linked to loss-of-function mutations of the MeCP2 gene for which no cure is yet available. Analysis of the noncoding RNA profile corresponding to the brain-abundant circular RNA (circRNA) and transcribed-ultraconserved region (T-UCR) populations in a mouse model of the disease reveals widespread dysregulation and enrichment in glutamatergic excitatory signaling and microtubule cytoskeleton pathways of the corresponding host genes. Proteomic analysis of hippocampal samples from affected individuals confirms abnormal levels of several cytoskeleton-related proteins together with key alterations in neurotransmission. Importantly, the glutamate receptor GRIA3 gene displays altered biogenesis in affected individuals and in vitro human cells and is influenced by expression of two ultraconserved RNAs. We also describe post-transcriptional regulation of SIRT2 by circRNAs, which modulates acetylation and total protein levels of GluR-1. As a consequence, both regulatory mechanisms converge on the biogenesis of AMPA receptors, with an effect on neuronal differentiation. In both cases, the noncoding RNAs antagonize MeCP2-directed regulation. Our findings indicate that noncoding transcripts may contribute to key alterations in Rett syndrome and are not only useful tools for revealing dysregulated processes but also molecules of biomarker value.

Nanopore-Based Direct RNA-Sequencing Reveals a High-Resolution Transcriptional Landscape of Porcine Reproductive and Respiratory Syndrome Virus.

The TRS-mediated discontinuous transcription process is a hallmark of Arteriviruses. Precise assessment of the intricate subgenomic RNA (sg mRNA) populations is required to understand the kinetics of viral transcription. It is difficult to reconstruct and comprehensively quantify splicing events using short-read sequencing, making the identification of transcription-regulatory sequences (TRS) particularly problematic. Here, we applied long-read direct RNA sequencing to characterize the recombined RNA molecules produced in porcine alveolar macrophages during early passage infection of porcine reproductive and respiratory syndrome virus (PRRSV). Based on sequencing two PRRSV isolates, namely XM-2020 and GD, we revealed a high-resolution and diverse transcriptional landscape in PRRSV. The data revealed intriguing differences in subgenomic recombination types between the two PRRSVs while also demonstrating TRS-independent heterogeneous subpopulation not previously observed in Arteriviruses. We find that TRS usage is a regulated process and share the common preferred TRS in both strains. This study also identified a substantial number of TRS-mediated transcript variants, including alternative-sg mRNAs encoding the same annotated ORF, as well as putative sg mRNAs encoded nested internal ORFs, implying that the genetic information encoded in PRRSV may be more intensively expressed. Epigenetic modifications have emerged as an essential regulatory layer in gene expression. Here, we gained a deeper understanding of m5C modification in poly(A) RNA, elucidating a potential link between methylation and transcriptional regulation. Collectively, our findings provided meaningful insights for redefining the transcriptome complexity of PRRSV. This will assist in filling the research gaps and developing strategies for better control of the PRRS.

Significance and Relevance of Spermatozoal RNAs to Male Fertility in Livestock.

Livestock production contributes to a significant part of the economy in developing countries. Although artificial insemination techniques brought substantial improvements in reproductive efficiency, male infertility remains a leading challenge in livestock. Current strategies for the diagnosis of male infertility largely depend on the evaluation of semen parameters and fail to diagnose idiopathic infertility in most cases. Recent evidences show that spermatozoa contains a suit of RNA population whose profile differs between fertile and infertile males. Studies have also demonstrated the crucial roles of spermatozoal RNA (spRNA) in spermatogenesis, fertilization, and early embryonic development. Thus, the spRNA profile may serve as unique molecular signatures of fertile sperm and may play pivotal roles in the diagnosis and treatment of male fertility. This manuscript provides an update on various spRNA populations, including protein-coding and non-coding RNAs, in livestock species and their potential role in semen quality, particularly sperm motility, freezability, and fertility. The contribution of seminal plasma to the spRNA population is also discussed. Furthermore, we discussed the significance of rare non-coding RNAs (ncRNAs) such as long ncRNAs (lncRNAs) and circular RNAs (circRNAs) in spermatogenic events.

Analysis of subcellular transcriptomes by RNA proximity labeling with Halo-seq.

Thousands of RNA species display nonuniform distribution within cells. However, quantification of the spatial patterns adopted by individual RNAs remains difficult, in part by a lack of quantitative tools for subcellular transcriptome analysis. In this study, we describe an RNA proximity labeling method that facilitates the quantification of subcellular RNA populations with high spatial specificity. This method, termed Halo-seq, pairs a light-activatable, radical generating small molecule with highly efficient Click chemistry to efficiently label and purify spatially defined RNA samples. We compared Halo-seq with previously reported similar methods and found that Halo-seq displayed a higher efficiency of RNA labeling, indicating that it is well suited to the investigation of small, precisely localized RNA populations. We then used Halo-seq to quantify nuclear, nucleolar and cytoplasmic transcriptomes, characterize their dynamic nature following perturbation, and identify RNA sequence features associated with their composition. Specifically, we found that RNAs containing AU-rich elements are relatively enriched in the nucleus. This enrichment becomes stronger upon treatment with the nuclear export inhibitor leptomycin B, both expanding the role of HuR in RNA export and generating a comprehensive set of transcripts whose export from the nucleus depends on HuR.

Aedes aegypti Piwi4 Structural Features Are Necessary for RNA Binding and Nuclear Localization.

The PIWI-interacting RNA (piRNA) pathway provides an RNA interference (RNAi) mechanism known from Drosophila studies to maintain the integrity of the germline genome by silencing transposable elements (TE). Aedes aegypti mosquitoes, which are the key vectors of several arthropod-borne viruses, exhibit an expanded repertoire of Piwi proteins involved in the piRNA pathway, suggesting functional divergence. Here, we investigate RNA-binding dynamics and subcellular localization of A. aegypti Piwi4 (AePiwi4), a Piwi protein involved in antiviral immunity and embryonic development, to better understand its function. We found that AePiwi4 PAZ (Piwi/Argonaute/Zwille), the domain that binds the 3′ ends of piRNAs, bound to mature (3′ 2′ O-methylated) and unmethylated RNAs with similar micromolar affinities (KD = 1.7 ± 0.8 μM and KD of 5.0 ± 2.2 μM, respectively; p = 0.05) in a sequence independent manner. Through site-directed mutagenesis studies, we identified highly conserved residues involved in RNA binding and found that subtle changes in the amino acids flanking the binding pocket across PAZ proteins have significant impacts on binding behaviors, likely by impacting the protein secondary structure. We also analyzed AePiwi4 subcellular localization in mosquito tissues. We found that the protein is both cytoplasmic and nuclear, and we identified an AePiwi4 nuclear localization signal (NLS) in the N-terminal region of the protein. Taken together, these studies provide insights on the dynamic role of AePiwi4 in RNAi and pave the way for future studies aimed at understanding Piwi interactions with diverse RNA populations.

CD63+ and MHC Class I+ Subsets of Extracellular Vesicles Produced by Wild-Type and CD47-Deficient Jurkat T Cells Have Divergent Functional Effects on Endothelial Cell Gene Expression.

T cells and endothelial cells engage in bidirectional communication that regulates angiogenesis and T cell transmigration. Extracellular vesicles (EVs) mediate intercellular communication by the transfer of bioactive molecules including RNAs. EVs produced by a given cell type are heterogeneous in their RNA content, but it is unclear how specific EV surface markers relate to their functional effects on target cells. Our previous work established that Jurkat T cell EVs bearing CD63, MHC-I, or CD47 surface markers contain distinct noncoding RNA populations. The present study reveals that CD63+ and MHC-I+ EVs from CD47-deficient Jurkat T cells are enriched in small non-coding RNAs relative to EVs from wild-type Jurkat T cells. CD47-deficient Jurkat T cells secrete more CD63+ and MHC-I+ EVs, but MHC-I+ EVs are selectively taken up more by human umbilical vein endothelial cells. Transcriptomics analysis of endothelial cells treated with CD63+ or MHC-I+ EVs showed surface marker- and CD47-dependent changes in gene expression in the target cells. Gene set enrichment analysis identified CD47-dependent, and surface marker-dependent effects of T cell EVs on VEGF and inflammatory signaling, cell cycle, and lipid and cholesterol metabolism. Thus, subsets of T cell EVs differentially regulate endothelial cell metabolism and inflammatory and angiogenic responses.

All you need to know about sperm RNAs.

Spermatogenesis generates a small and highly specialised type of cell that is apparently incapable of transcription and translation. For many years, this dogma was supported by the assumption that (i) the compact sperm nucleus, resulting from the substitution of histones by protamine during spermatogenesis, renders the genome inaccessible to the transcriptional machinery; and (ii) the loss of most organelles, including endoplasmic reticulum and ribosomes, limits or prevents translational activity. Despite these observations, several types of coding and non-coding RNAs have been identified in human sperm. Their functional roles, particularly during fertilisation and embryonic development, are only now becoming apparent. This review aimed to summarise current knowledge of the origin, types and functional roles of sperm RNAs, and to evaluate the clinical benefits of employing these transcripts as biomarkers of male fertility and reproductive outcomes. The possible contribution of sperm RNAs to intergenerational or transgenerational phenotypic inheritance is also addressed. A comprehensive literature search on PubMed was conducted using the search terms 'sperm' AND 'RNA'. Searches focussed upon articles written in English and published prior to August 2020. The development of more sensitive and accurate RNA technologies, including RNA sequencing, has enabled the identification and characterisation of numerous transcripts in human sperm. Though a majority of these RNAs likely arise during spermatogenesis, other data support an epididymal origin of RNA transmitted to maturing sperm by extracellular vesicles. A minority may also be synthesised by de novo transcription in mature sperm, since a small portion of the sperm genome remains packed by histones. This complex RNA population has important roles in paternal chromatin packaging, sperm maturation and capacitation, fertilisation, early embryogenesis and developmental maintenance. In recent years, additional lines of evidence from animal models support a role for sperm RNAs in intergenerational or transgenerational inheritance, modulating both the genotype and phenotype of progeny. Importantly, several reports indicate that the sperm RNA content of fertile and infertile men differs considerably and is strongly modulated by the environment, lifestyle and pathological states. Transcriptional profiling has considerable potential for the discovery of fertility biomarkers. Understanding the role of sperm transcripts and comparing the sperm RNA fingerprint of fertile and infertile men could help to elucidate the regulatory pathways contributing to male factor infertility. Such data might also provide a molecular explanation for several causes of idiopathic male fertility. Ultimately, transcriptional profiling may be employed to optimise ART procedures and overcome some of the underlying causes of male infertility, ensuring the birth of healthy children.

Mutually Orthogonal Bioconjugation of Vinyl Nucleosides for RNA Metabolic Labeling.

We report a strategy for the orthogonal conjugation of the vinyl nucleosides, 5-vinyluridine (5-VU) and 2-vinyladenosine (2-VA), via selective reactivity with maleimide and tris(2-carboxyethyl)phosphine (TCEP), respectively. The orthogonality was investigated using density functional theory (DFT) and confirmed by reactions with vinyl nucleosides. Further, these chemistries were used to modify RNA for fluorescent cell imaging. These reactions allow for the expanded use of RNA metabolic labeling to study nascent RNA expression within different RNA populations.

SARS-CoV-2 Quasispecies Provides an Advantage Mutation Pool for the Epidemic Variants.

ABSTRACTThe dynamics of quasispecies afford RNA viruses a great fitness on cell tropism and host range. To study the quasispecies features and the intra-host evolution of SARS-CoV-2, we collected nine confirmed patients and sequenced the haplotypes of spike gene using a single-molecule real-time platform. Fourteen samples were extracted from sputum, nasopharyngeal swabs, or stool, which in total produced 283,655 high-quality circular consensus sequences. We observed a stable quasispecies structure that one master mutant (mean abundance ∼0.70), followed by numerous minor mutants (mean abundance ∼1.21 × 10−3). Under high selective pressure, minor mutants may obtain a fitness advantage and become the master ones. The later predominant substitution D614G existed in the minor mutants of more than one early patient. An epidemic variant had a possibility to be independently originated from multiple hosts. The mutant spectrums covered ∼85% amino acid variations of public genomes (GISAID; frequency ≥ 0.1) and likely provided an advantage mutation pool for the current/future epidemic variants. Notably, 32 of 35 collected antibody escape substitutions were preexistent in the early quasispecies. Virus populations in different tissues/organs revealed potentially independent replications. The quasispecies complexity of sputum samples was significantly lower than that of nasopharyngeal swabs (P = 0.02). Evolution analysis revealed that three continuous S2 domains (HR1, CH, and CD) had undergone a positive selection. Cell fusion-related domains may play a crucial role in adapting to the intrahost immune system. Our findings suggested that future epidemiologic investigations and clinical interventions should consider the quasispecies information that has missed by routine single consensus genome.IMPORTANCE RNA virus population in a host does not consist of a consensus single haplotype but rather an ensemble of related sequences termed quasispecies. The dynamics of quasispecies afford SARS-CoV-2 a great ability on genetic fitness during intrahost evolution. The process is likely achieved by changing the genetic characteristics of key functional genes, such as the spike glycoprotein. Previous studies have applied the next-generation sequencing (NGS) technology to evaluate the quasispecies of SARS-CoV-2, and results indicated a low genetic diversity of the spike gene. However, the NGS platform cannot directly obtain the full haplotypes without assembling, and it is also difficult to predict the extremely low-frequency variations. Therefore, we introduced a single-molecule real-time technology to directly obtain the haplotypes of the RNA population and further study the quasispecies features and intrahost evolution of the spike gene.

An essential role for the piRNA pathway in regulating the ribosomal RNA pool in C.elegans.

Piwi-interacting RNAs (piRNAs) are RNA effectors with key roles in maintaining genome integrity and promoting fertility in metazoans. In Caenorhabditis elegans loss of piRNAs leads to a transgenerational sterility phenotype. The plethora of piRNAs and their ability to silence transcripts with imperfect complementarity have raised several (non-exclusive) models for the underlying drivers of sterility. Here, we report the extranuclear and transferable nature of the sterility driver, its suppression via mutations disrupting the endogenous RNAi and poly-uridylation machinery, and copy-number amplification at the ribosomal DNA locus. In piRNA-deficient animals, several small interfering RNA (siRNA) populations become increasingly overabundant in the generations preceding loss of germline function, including ribosomal siRNAs (risiRNAs). A concomitant increase in uridylated sense rRNA fragments suggests that poly-uridylation may potentiate RNAi-mediated gene silencing of rRNAs. We conclude that loss of the piRNA machinery allows for unchecked amplification of siRNA populations, originating from abundant highly structured RNAs, to deleterious levels.

Locked-nucleotide antagonists to varicella zoster virus small non-coding RNA block viral growth and have potential as an anti-viral therapy

Herpes zoster (HZ) remains a significant health burden with millions of cases in North America and Europe annually. HZ is frequently followed by long-term pain or post-herpetic neuralgia (PHN). Although effective vaccines for HZ are available, currently used nucleotide analogues often have limited effectiveness against HZ and especially PHN, so there remains a need for additional antiviral therapies for HZ. We recently identified a population of small non-coding RNA (sncRNA) encoded by Varicella Zoster Virus (VZV) and showed that single locked-nucleic acid antagonists (LNAA) to some sncRNA can modulate VZV replication in cell culture. In this work, we explored the antiviral effects of combinations of LNAA oligonucleotides targeting VZVsncRNA. Combinations of LNAA targeting three VZVsncRNA encoded in and near a critical viral regulatory gene were additive, achieving 96 % reduction in virus growth in a cell line. VZV growth was also inhibited by more than 90 % in primary human skin fibroblast cultures by individual and combinations of LNAA to VZVsncRNA. The inhibition by VZVsncRNA was specific and not a consequence of innate immune responses since LNAA to a different VZVsncRNA enhanced VZV growth. Targeted VZVsncRNA lack homologous sequences in the human transcriptome suggesting that LNAA to them would have reduced cytotoxicity if used as therapeutics. These results support further development of oligonucleotides targeting VZVsncRNA as a novel treatment for HZ.

Clinical Evidence That Deficiency of the Methylcytosine Dioxygenase TET2 Gene Promotes the Clonal Expansion and Therapeutic Efficacy of CD19-Directed T Cells

Adoptive transfer of autologous chimeric antigen receptor (CAR) T cells that target tumor antigens can induce remissions in patients with hematologic malignancies. A major obstacle to successful clinical responses in some patient populations has been infusion with ineffective cell populations that are terminally differentiated, exhausted, poorly proliferative and short-lived. Understanding and manipulating the mechanistic pathways that sustain T cell potency will potentially unlock durable responses to tumors. Here we present an unusual case that helps clarify determinants of CAR T cell persistence. A 78-year-old individual with advanced chronic lymphocytic leukemia (CLL) was infused with CAR T cells that target CD19 (CTL019). CAR T cell expansion was delayed and not evident until 50 days following infusion and anti-tumor activity was observed in the peripheral blood, lymph nodes as well as bone marrow, accompanied by complete remission that has been sustained for over 4.5 years. Retrospective analyses revealed that at the peak of in vivo CAR T cell activity, over 94% percent of CD8+ CTL019 cells originated from a single clone in which lentiviral vector-mediated insertion of the CAR transgene led to monoallelic loss of TET2 . Accordingly, examination of polyadenylated TET2 RNA populations showed the appearance of new chimeric RNAs that spliced from TET2 exon 9 into the vector and terminated, truncating the encoded protein. Further analysis revealed that this patient had a germline polymorphism in the other allele of TET2 . This was a missense variant at amino acid 1879 (exon 11), converting the wildtype residue, glutamic acid, to glutamine (protein change: pE1879Q; cDNA change: c.5635G>C; allele frequency of approximately 50%). TET2 deficient CTL019 cells exhibited reduced cytosine hydroxymethylation and a chromatin accessibility profile measured by ATAC-Seq consistent with altered T cell differentiation. Furthermore, at the peak of in vivo engraftment, CAR T cells from this subject showed an early (central) memory phenotype which differs from CTL019 cells of other complete responding patients that are instead characterized by late (effector) memory differentiation at the height of their response. Experimental knock-down of TET2 in healthy donor CAR T cells altered their differentiation and augmented CTL019 effector function in a manner that was consistent with the effects of TET2 deficiency observed in this subject. Our findings suggest that the profound clonal expansion of a single TET2 deficient CAR T cell resulted in a complete and durable remission of CLL, highlighting the potential for epigenetic modulation of ineffective patient T cells to generate curative T cell products. DisclosuresFraietta:Novartis: Patents & Royalties: patents, Research Funding. Morrissette:Novartis: Consultancy. Liu:Novartis: Patents & Royalties. Young:Novartis: Research Funding. Levine:Novartis Pharmaceuticals Corporation: Patents & Royalties, Research Funding; Tmunity Therapeutics: Equity Ownership, Research Funding; GE Healthcare: Consultancy; Brammer Bio: Consultancy. Zhao:Novartis: Patents & Royalties. Kalos:Novartis: Research Funding. Porter:Incyte: Honoraria; Servier: Honoraria, Other: Travel reimbursement; Genentech/Roche: Employment, Other: Family member employment, stock ownship - family member; Immunovative Therapies: Other: Member DSMB; Novartis: Honoraria, Patents & Royalties, Research Funding. Lacey:Novartis: Research Funding; Genentech: Honoraria. June:Tmunity Therapeutics: Equity Ownership, Research Funding; WIRB/Copernicus Group: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celldex: Honoraria, Membership on an entity's Board of Directors or advisory committees; Immune Design: Equity Ownership, Membership on an entity's Board of Directors or advisory committees; Novartis: Patents & Royalties, Research Funding. Melenhorst:Novartis: Research Funding.

Global and Targeted RNA-Sequencing Assessment of Cancer Testis Antigen Re-Expression Following Azacitidine in Myelodysplastic Syndromes

Background and Aim: Azanucleoside DNA-hypomethylating agents have remarkable clinical activity in myelodysplastic sindromes (MDS) and acute myeloid leukemia (AML), particularly at low, non-cytotoxic doses favoring hypomethylation over cytotoxicity. Cancer/testis antigens (CTAs) encoding immunogenic proteins are normally expressed in the testicles and trophoblastic cells of the ovary and in some malignant tissues: these antigens are expressed to variable degrees particularly in melanoma, multiple myeloma, lung and breast cancer, but expression is quite infrequent in myeloid malignancies. However, they have been shown to undergo derepression after the use of demethylating agents in cell line models. Due to their tumor-restricted expression and their potential immunogenicity, we aimed to comprehensively study CTAs re-expression following azacitidine (AZA) treatment in MDS patients.Methods and Patients: We enrolled 32 patients with MDS (n =22) and chronic myelomonocytic leukemia (n = 10). Deep RNA-sequencing (RNA-Seq) was performed on RNA extracted from blood at day 0 and at day +28 of the first AZA cycle, from 12 patients and compared with 9 healthy donors. RNA-Seq (mean 106 million reads per sample) was sequenced on a HiSeq4000. After TopHat alignment and Cufflinks assembly, expression profile was calculated as FPKM. Targeted RNA-Seq (tRNA-Seq) designed to capture 214 CTA genes was performed on a validation cohort of 20 patients (40 samples). Coverage analysis as well as mapping the reads and alignment was done using the Ion Torrent Browser Suite. We used LIMMA package or empirical analysis of digital gene expression data in R (edgeR) to identify differential expressed genes between day 0 and d+28 samples, and between responder and non-responder groups. A log2 fold change (FC) > 10 and a p<0.1 were considered as significant.Results: The median age of the cohort was 69 (range 48-81 years). The cohort consisted of MDS (75%) and MDS/MPN (25%) patients. MDS Patients were stratified based on IPSS as low (8%), intermediate-1 (23%), intermediate-2 (46%), and high (23%) risk groups and MDS/MPN patients based on CPSS as low (20%), intermediate-1 (20%) and intermediate-2 (60%) . A homogeneous treatment schedule with AZA (75 mg/m2/day x 7 days every 4 weeks) was received by each patient. Seventeen patients were classified as responders (achieving complete remission (CR), marrow CR, partial remission, or stable disease with haematological improvement according to the 2006 International Working Group (IWG) criteria for evaluating MDS. Five patients progressed within the 6 first cycles; otherwise response evaluation was performed after the sixth cycle.In the screening RNA-Seq experiment, MAGEB3, AKAP4, POTEB2 and TMPRSS12 were non-expressed in control or MDS pretreatment, but re-expressed after 1 cycle of Aza; ADAM29, and XAGE1E were non expressed in controls, but expressed in MDS samples: an upregulation which increased after the first cycle.Regarding the tRNA-Seq cohort, on average, we achieved 3.888.308 mapped reads per sample (p25-p75, 1.2x106 - 5.9x106) which represents a sufficient depth for digital gene expression profiling of 214 genes. CTAs re-expressed after one AZA cycle in tRNA-Seq validation cohort: CT45A5, BRDT, TFDP3, TSPY3, TSPY2 showed an FC >10 and a p<0.05 and ADAM29. SLCO6A1, PAGE5 and RQCD1 showed a p<0.05 but a FC<10. CTAs with an increased expression in MDS after Aza: FMR1NB, ADAM29, and XAGE1E were not expressed in controls, re-expressed in MDS before treatment and showed an upregulation after AZA. Responders vs non responders: strikingly CT45A4, and SPAG17 were markedly overexpressed in non-responders than in responders (FC >10, = 0.006; p=0.01), and viceversa for PRM1 and TFDP3 (FC>10; p=0.01, p 0.006).Conclusions: This is the first study to assess, in a global and sentitive design, the targetable re-expression of CTAs following therapeutical demethylation in MDS. Strikingly, in addition to finding attractive targets for cancer immunotherapy, some CTAs were exclusively and markedly re-expressed in non-responders. We also show that targeted RNA-Seq can overcome the test posed by the wide range of the RNA population by sequencing targeted genes providing a vast enrichment of weakly expressed transcripts. DisclosuresNo relevant conflicts of interest to declare.

Molecular insights on potato yellow vein crinivirus infections in the highlands of Colombia.

Potato yellow vein virus (PYVV) was detected in potatoes grown in the Central highlands, north of Bogotá (~3000 m altitude), Colombia. At this altitude viral whitefly vectors are largely absent, but infection persists because of the use of uncertified tubers. Plants with typical PYVV-induced yellowing symptoms, as well as with atypical yellowing or non-symptomatic symptoms were sampled at three separate geographical locations. PYVV presence was assessed by RT-PCR, and several plants were subjected to high-throughput sequencing (HTS) of their small RNA (sRNA) populations. Complete or almost complete sequences of four PYVV isolates were thus reconstructed, all from symptomatic plants. Three viral isolates infected plants singly, while the fourth co-infected the plant together with a potyvirus. Relative proportions of sRNAs to each of the three crinivirus genomic RNAs were found to remain comparable among the four infections. Genomic regions were identified as hotspots of sRNA formation, or as regions that poorly induced sRNAs. Furthermore, PYVV titres in the mixed versus single infections remained comparable, indicating an absence of synergistic/antagonistic effects of the potyvirus on the accumulation of PYVV. Daughter plants raised in the greenhouse from tubers of the infected, field-sampled plants displayed mild PYVV infection symptoms that disappeared with time, demonstrating the occurrence of recovery and asymptomatic infection phenotypes in this pathosystem.

Participation of Multifunctional RNA in Replication, Recombination and Regulation of Endogenous Plant Pararetroviruses (EPRVs)

Pararetroviruses, taxon Caulimoviridae, are typical of retroelements with reverse transcriptase and share a common origin with retroviruses and LTR retrotransposons, presumably dating back 1.6 billion years and illustrating the transition from an RNA to a DNA world. After transcription of the viral genome in the host nucleus, viral DNA synthesis occurs in the cytoplasm on the generated terminally redundant RNA including inter- and intra-molecule recombination steps rather than relying on nuclear DNA replication. RNA recombination events between an ancestral genomic retroelement with exogenous RNA viruses were seminal in pararetrovirus evolution resulting in horizontal transmission and episomal replication. Instead of active integration, pararetroviruses use the host DNA repair machinery to prevail in genomes of angiosperms, gymnosperms and ferns. Pararetrovirus integration – leading to Endogenous ParaRetroViruses, EPRVs – by illegitimate recombination can happen if their sequences instead of homologous host genomic sequences on the sister chromatid (during mitosis) or homologous chromosome (during meiosis) are used as template. Multiple layers of RNA interference exist regulating episomal and chromosomal forms of the pararetrovirus. Pararetroviruses have evolved suppressors against this plant defense in the arms race during co-evolution which can result in deregulation of plant genes. Small RNAs serve as signaling molecules for Transcriptional and Post-Transcriptional Gene Silencing (TGS, PTGS) pathways. Different populations of small RNAs comprising 21–24 nt and 18–30 nt in length have been reported for Citrus, Fritillaria, Musa, Petunia, Solanum and Beta. Recombination and RNA interference are driving forces for evolution and regulation of EPRVs.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus.

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Dynamics of small non-coding RNAs in bovine scNT embryos through the maternal-to-embryonic transition†.

The efficiency of somatic cell nuclear transfer (scNT) for production of viable offspring is relatively low as compared to in vitro fertilization (IVF), presumably due to deficiencies in epigenetic reprogramming of the donor cell genome. Such defects may also involve the population of small non-coding RNAs (sncRNAs), which are important during early embryonic development. The objective of this study was to examine dynamic changes in relative abundance of sncRNAs during the maternal-to-embryonic transition (MET) in bovine embryos produced by scNT as compared to IVF by using RNA sequencing. When comparing populations of miRNA in scNT versus IVF embryos, only miR-2340, miR-345, and miR34a were differentially expressed in morulae, though many more miRNAs were differentially expressed when comparing across developmental stages. Also of interest, distinct populations of piwi-interacting like RNAs (pilRNAs) were identified in bovine embryos prior to and during embryonic genome activation (EGA) as compared bovine embryos post-EGA and differentiated cells. Overall, sncRNA sequencing analysis of preimplantation embryos revealed largely similar profiles of sncRNAs for IVF and scNT embryos at the 2-cell, 8-cell, morula, and blastocyst stages of development. However, these sncRNA profiles, including miRNA, piRNA, and tRNA fragments, were notably distinct prior to and after completion of the MET.

Multiple Viral Infections Detected in Phytophthora condilina by Total and Small RNA Sequencing.

Marine oomycetes have recently been shown to be concurrently infected by (−)ssRNA viruses of the order Bunyavirales. In this work, even higher virus variability was found in a single isolate of Phytophthora condilina, a recently described member of Phytophthora phylogenetic Clade 6a, which was isolated from brackish estuarine waters in southern Portugal. Using total and small RNA-seq the full RdRp of 13 different potential novel bunya-like viruses and two complete toti-like viruses were detected. All these viruses were successfully confirmed by reverse transcription polymerase chain reaction (RT-PCR) using total RNA as template, but complementarily one of the toti-like and five of the bunya-like viruses were confirmed when dsRNA was purified for RT-PCR. In our study, total RNA-seq was by far more efficient for de novo assembling of the virus sequencing but small RNA-seq showed higher read numbers for most viruses. Two main populations of small RNAs (21 nts and 25 nts-long) were identified, which were in accordance with other Phytophthora species. To the best of our knowledge, this is the first study using small RNA sequencing to identify viruses in Phytophthora spp.