Internal oligo(dT) priming in bulk and single cell RNA sequencing

Significant advances in RNA sequencing have been recently made possible by the use of oligo(dT) primers for simultaneous mRNA enrichment and reverse transcription priming. The associated increase in efficiency has enabled more economical bulk RNA sequencing methods as well as the advent of high throughput single cell RNA sequencing, now already one of the most widely adopted new methods in the study of transcriptomics. However, the effects of off-target oligo(dT) priming on gene expression quantification have not been fully appreciated. In the present study, we describe the extent, the possible causes, and the consequences of internal oligo(dT) priming across multiple publicly available datasets obtained from a variety of bulk and single cell RNA sequencing platforms. In order to explore and address this issue, we developed a computational algorithm for identification of sequencing read alignments that likely resulted from internal oligo(dT) priming and their subsequent removal from the data. Directly comparing filtered datasets to those obtained by an alternative method reveals significant improvements in gene expression measurement. Finally, we infer a list of genes whose expression quantification is most likely to be affected by internal oligo(dT) priming.

Arg-Vasotocin Directly Activates Isotocin Receptors and Induces COX2 Expression in Ovoviviparous Guppies

Oxytocin (OT) is a crucial regulator of reproductive behaviors, including parturition in mammals. Arg-vasopressin (AVP) is a nonapeptide homologous to Arg-vasotocin (AVT) in teleosts that has comparable affinity for the OT receptor. In the present study, ovoviviparous guppies (Poecilia reticulata) were used to study the effect of AVT on delivery mediated by the activation of prostaglandin (PG) biosynthesis via isotocin (IT) receptors (ITRs). One copy each of it and avt and two copies of itrs were identified in guppies. The results of the affinity assay showed that various concentrations of AVT and IT (10−6, 10−7, and 10−8 mol/L) significantly activated itr1 (P < 0.05). In vitro experiments revealed significant upregulation (P < 0.05) of cyclooxygenase 2 (cox2), which is the rate-limiting enzyme involved in PG biosynthesis, and itr1 by AVT and IT. Furthermore, dual in situ hybridization detected positive signals for itr1 and cox2 at the same site, implying that ITR1 may regulate cox2 gene expression. Measurement of prostaglandin F2a (PGF2a) concentrations showed that AVT induced PGF2a synthesis (P < 0.05) and that the effect of IT was not significant. Finally, intraperitoneal administration of PGF2a significantly induced premature parturition of guppies. This study is the first to identify and characterize AVT and ITRs in guppies. The findings suggest that AVT promotes PG biosynthesis via ITR and that PGF2a induces delivery behavior in ovoviviparous guppies.

GPR55 Receptor Activation by the N-Acyl Dopamine Family Lipids Induces Apoptosis in Cancer Cells via the Nitric Oxide Synthase (nNOS) Over-Stimulation

GPR55 is a GPCR of the non-CB1/CB2 cannabinoid receptor family, which is activated by lysophosphatidylinositol (LPI) and stimulates the proliferation of cancer cells. Anandamide, a bioactive lipid endocannabinoid, acts as a biased agonist of GPR55 and induces cancer cell death, but is unstable and psychoactive. We hypothesized that other endocannabinoids and structurally similar compounds, which are more hydrolytically stable, could also induce cancer cell death via GPR55 activation. We chemically synthesized and tested a set of fatty acid amides and esters for cell death induction via GPR55 activation. The most active compounds appeared to be N-acyl dopamines, especially N-docosahexaenoyl dopamine (DHA-DA). Using a panel of cancer cell lines and a set of receptor and intracellular signal transduction machinery inhibitors together with cell viability, Ca2+, NO, ROS (reactive oxygen species) and gene expression measurement, we showed for the first time that for these compounds, the mechanism of cell death induction differed from that published for anandamide and included neuronal nitric oxide synthase (nNOS) overstimulation with concomitant oxidative stress induction. The combination of DHA-DA with LPI, which normally stimulates cancer proliferation and is increased in cancer setting, had an increased cytotoxicity for the cancer cells indicating a therapeutic potential.

Transposase-assisted tagmentation of RNA/DNA hybrid duplexes

Tn5-mediated transposition of double-strand DNA has been widely utilized in various high-throughput sequencing applications. Here, we report that the Tn5 transposase is also capable of direct tagmentation of RNA/DNA hybrids in vitro. As a proof-of-concept application, we utilized this activity to replace the traditional library construction procedure of RNA sequencing, which contains many laborious and time-consuming processes. Results of Transposase-assisted RNA/DNA hybrids Co-tagmEntation (termed ‘TRACE-seq’) are compared to traditional RNA-seq methods in terms of detected gene number, gene body coverage, gene expression measurement, library complexity, and differential expression analysis. At the meantime, TRACE-seq enables a cost-effective one-tube library construction protocol and hence is more rapid (within 6 hr) and convenient. We expect this tagmentation activity on RNA/DNA hybrids to have broad potentials on RNA biology and chromatin research.

Breast tumor stiffness instructs bone metastasis via mechanical memory

The mechanical microenvironment of primary breast tumors plays a substantial role in promoting tumor progression. While the transitory response of cancer cells to pathological stiffness in their native microenvironment has been well described, it is still unclear how mechanical stimuli in the primary tumor influence distant, late-stage metastatic phenotypes across time and space in absentia. Here, we show that primary tumor stiffness promotes stable, non-genetically heritable phenotypes in breast cancer cells. This mechanical memory instructs cancer cells to adopt and maintain increased cytoskeletal dynamics, traction force, and 3D invasion in vitro, in addition to promoting osteolytic bone metastasis in vivo. Furthermore, we established a mechanical conditioning (MeCo) score comprised of mechanically-regulated genes as a global gene expression measurement of tumor stiffness response. Clinically, we show that a high MeCo score is strongly associated with bone metastasis in patients. Using a discovery approach, we mechanistically traced mechanical memory in part to ERK-mediated mechanotransductive activation of RUNX2, an osteogenic gene bookmarker and bone metastasis driver. The combination of these RUNX2 traits permits the stable transactivation of osteolytic target genes that remain upregulated after cancer cells disseminate from their activating microenvironment in order to modify a distant microenvironment. Using genetic, epigenetic, and functional approaches, we were able to simulate, repress, select and extend RUNX2-mediated mechanical memory and alter cancer cell behavior accordingly. In concert with previous studies detailing the influence of biochemical properties of the primary tumor stroma on distinct metastatic phenotypes, our findings detailing the influence of biomechanical properties support a generalized model of cancer progression in which the integrated properties of the primary tumor microenvironment govern the secondary tumor microenvironment, i.e., soil instructs soil.

Not Only Mutations Matter: Molecular Picture of Acute Myeloid Leukemia Emerging from Transcriptome Studies

The last two decades of genome-scale research revealed a complex molecular picture of acute myeloid leukemia (AML). On the one hand, a number of mutations were discovered and associated with AML diagnosis and prognosis; some of them were introduced into diagnostic tests. On the other hand, transcriptome studies, which preceded AML exome and genome sequencing, remained poorly translated into clinics. Nevertheless, gene expression studies significantly contributed to the elucidation of AML pathogenesis and indicated potential therapeutic directions. The power of transcriptomic approach lies in its comprehensiveness; we can observe how genome manifests its function in a particular type of cells and follow many genes in one test. Moreover, gene expression measurement can be combined with mutation detection, as high-impact mutations are often present in transcripts. This review sums up 20 years of transcriptome research devoted to AML. Gene expression profiling (GEP) revealed signatures distinctive for selected AML subtypes and uncovered the additional within-subtype heterogeneity. The results were particularly valuable in the case of AML with normal karyotype which concerns up to 50% of AML cases. With the use of GEP, new classes of the disease were identified and prognostic predictors were proposed. A plenty of genes were detected as overexpressed in AML when compared to healthy control, includingKIT,BAALC,ERG,MN1,CDX2,WT1,PRAME,andHOXgenes. High expression of these genes constitutes usually an unfavorable prognostic factor. Upregulation ofFLT3andNPM1genes, independent on their mutation status, was also reported in AML and correlated with poor outcome. However, transcriptome is not limited to the protein-coding genes; other types of RNA molecules exist in a cell and regulate genome function. It was shown that microRNA (miRNA) profiles differentiated AML groups and predicted outcome not worse than protein-coding gene profiles. For example, upregulation ofmiR-10a,miR-10b, andmiR-196band downregulation ofmiR-192were found as typical of AML withNPM1mutation whereas overexpression ofmiR-155was associated withFLT3-internal tandem duplication (FLT3-ITD). Development of high-throughput technologies and microarray replacement by next generation sequencing (RNA-seq) enabled uncovering a real variety of leukemic cell transcriptomes, reflected by gene fusions, chimeric RNAs, alternatively spliced transcripts, miRNAs, piRNAs, long noncoding RNAs (lncRNAs), and their special type, circular RNAs. Many of them can be considered as AML biomarkers and potential therapeutic targets. The relations between particular RNA puzzles and other components of leukemic cells and their microenvironment, such as exosomes, are now under investigation. Hopefully, the results of this research will shed the light on these aspects of AML pathogenesis which are still not completely understood.

Examination of Independent Prognostic Power of Gene Expressions and Histopathological Imaging Features in Cancer

Cancer prognosis is of essential interest, and extensive research has been conducted searching for biomarkers with prognostic power. Recent studies have shown that both omics profiles and histopathological imaging features have prognostic power. There are also studies exploring integrating the two types of measurements for prognosis modeling. However, there is a lack of study rigorously examining whether omics measurements have independent prognostic power conditional on histopathological imaging features, and vice versa. In this article, we adopt a rigorous statistical testing framework and test whether an individual gene expression measurement can improve prognosis modeling conditional on high-dimensional imaging features, and a parallel analysis is conducted reversing the roles of gene expressions and imaging features. In the analysis of The Cancer Genome Atlas (TCGA) lung adenocarcinoma and liver hepatocellular carcinoma data, it is found that multiple individual genes, conditional on imaging features, can lead to significant improvement in prognosis modeling; however, individual imaging features, conditional on gene expressions, only offer limited prognostic power. Being among the first to examine the independent prognostic power, this study may assist better understanding the “connectedness” between omics profiles and histopathological imaging features and provide important insights for data integration in cancer modeling.