A longer time to relapse is associated with a larger increase in differences between paired primary and recurrent IDH wild-type glioblastomas at both the transcriptomic and genomic levels.

Glioblastoma (GBM) is the most common malignant brain tumor in adults, which remains incurable and often recurs rapidly after initial therapy. While large efforts have been dedicated to uncover genomic/transcriptomic alternations associated with the recurrence of GBMs, the evolutionary trajectories of matched pairs of primary and recurrent (P-R) GBMs remain largely elusive. It remains challenging to identify genes associated with time to relapse (TTR) and construct a stable and effective prognostic model for predicting TTR of primary GBM patients. By integrating RNA-sequencing and genomic data from multiple datasets of patient-matched longitudinal GBMs of isocitrate dehydrogenase wild-type (IDH-wt), here we examined the associations of TTR with heterogeneities between paired P-R GBMs in gene expression profiles, tumor mutation burden (TMB), and microenvironment. Our results revealed a positive correlation between TTR and transcriptomic/genomic differences between paired P-R GBMs, higher percentages of non-mesenchymal-to-mesenchymal transition and mesenchymal subtype for patients with a short TTR than for those with a long TTR, a high correlation between paired P-R GBMs in gene expression profiles and TMB, and a negative correlation between the fitting level of such a paired P-R GBM correlation and TTR. According to these observations, we identified 55 TTR-associated genes and thereby constructed a seven-gene (ZSCAN10, SIGLEC14, GHRHR, TBX15, TAS2R1, CDKL1, and CD101) prognostic model for predicting TTR of primary IDH-wt GBM patients using univariate/multivariate Cox regression analyses. The risk scores estimated by the model were significantly negatively correlated with TTR in the training set and two independent testing sets. The model also segregated IDH-wt GBM patients into two groups with significantly divergent progression-free survival outcomes and showed promising performance for predicting 1-, 2-, and 3-year progression-free survival rates in all training and testing sets. Our findings provide new insights into the molecular understanding of GBM progression at recurrence and potential targets for therapeutic treatments.

Long-term hematopoietic transfer of the anti-cancer and lifespan-extending capabilities of a genetically engineered blood system by transplantation of bone marrow mononuclear cells.

A causal relationship exists among the aging process, organ decay and disfunction, and the occurrence of various diseases including cancer. A genetically engineered mouse model, termed Klf1K74R/K74R or Klf1(K74R), carrying mutation on the well-conserved sumoylation site of the hematopoietic transcription factor KLF1/EKLF has been generated that possesses extended lifespan and healthy characteristics, including cancer resistance. We show that the healthy longevity characteristics of the Klf1(K74R) mice, as exemplified by their higher anti-cancer capability, are likely gender-, age-, and genetic background-independent. Significantly, the anti-cancer capability, in particular that against melanoma as well as hepatocellular carcinoma, and lifespan-extending property of Klf1(K74R) mice, could be transferred to wild-type mice via transplantation of their bone marrow mononuclear cells at a young age of the latter. Furthermore, NK(K74R) cells carry higher in vitro cancer cell-killing ability than wild-type NK cells. Targeted/global gene expression profiling analysis has identified changes in the expression of specific proteins, including the immune checkpoint factors PDCD and CD274, and cellular pathways in the leukocytes of the Klf1(K74R) that are in the directions of anti-cancer and/or anti-aging. This study demonstrates the feasibility of developing a transferable hematopoietic/blood system for long-term anti-cancer and, potentially, for anti-aging.

Detecting intragenic trans-splicing events from non-co-linearly spliced junctions by hybrid sequencing.

Trans-spliced RNAs (ts-RNAs) are a type of non-co-linear (NCL) transcripts that consist of exons in an order topologically inconsistent with the corresponding DNA template. Detecting ts-RNAs is often interfered by experimental artifacts, circular RNAs (circRNAs) and genetic rearrangements. Particularly, intragenic ts-RNAs, which are derived from separate precursor mRNA molecules of the same gene, are often mistaken for circRNAs through analyses of RNA-seq data. Here we developed a bioinformatics pipeline (NCLscan-hybrid), which integrated short and long RNA-seq reads to minimize false positives and proposed out-of-circle and rolling-circle long reads to distinguish between intragenic ts-RNAs and circRNAs. Combining NCLscan-hybrid screening and multiple experimental validation steps successfully confirmed that four NCL events, which were previously regarded as circRNAs in databases, originated from trans-splicing. CRISPR-based endogenous genome modification experiments further showed that flanking intronic complementary sequences can significantly contribute to ts-RNA formation, providing an efficient/specific method to deplete ts-RNAs. We also experimentally validated that one ts-RNA (ts-ARFGEF1) played an important role for p53-mediated apoptosis through affecting the PERK/eIF2a/ATF4/CHOP signaling pathway in breast cancer cells. This study thus described both bioinformatics procedures and experimental validation steps for rigorous characterization of ts-RNAs, expanding future studies for identification, biogenesis, and function of these important but understudied transcripts.

A-to-I RNA editing contributes to the persistence of predicted damaging mutations in populations.

Adenosine-to-inosine (A-to-I) RNA editing is a very common co-/posttranscriptional modification that can lead to A-to-G changes at the RNA level and compensate for G-to-A genomic changes to a certain extent. It has been shown that each healthy individual can carry dozens of missense variants predicted to be severely deleterious. Why strongly detrimental variants are preserved in a population and not eliminated by negative natural selection remains mostly unclear. Here, we ask if RNA editing correlates with the burden of deleterious A/G polymorphisms in a population. Integrating genome and transcriptome sequencing data from 447 human lymphoblastoid cell lines, we show that nonsynonymous editing activities (prevalence/level) are negatively correlated with the deleteriousness of A-to-G genomic changes and positively correlated with that of G-to-A genomic changes within the population. We find a significantly negative correlation between nonsynonymous editing activities and allele frequency of A within the population. This negative editing-allele frequency correlation is particularly strong when editing sites are located in highly important genes/loci. Examinations of deleterious missense variants from the 1000 Genomes Project further show a significantly higher proportion of rare missense mutations for G-to-A changes than for other types of changes. The proportion for G-to-A changes increases with increasing deleterious effects of the changes. Moreover, the deleteriousness of G-to-A changes is significantly positively correlated with the percentage of editing enzyme binding motifs at the variants. Overall, we show that nonsynonymous editing is associated with the increased burden of G-to-A missense mutations in healthy individuals, expanding RNA editing in pathogenomics studies.

An Evolutionary Landscape of A-to-I RNA Editome across Metazoan Species.

Adenosine-to-inosine (A-to-I) editing is widespread across the kingdom Metazoa. However, for the lack of comprehensive analysis in nonmodel animals, the evolutionary history of A-to-I editing remains largely unexplored. Here, we detect high-confidence editing sites using clustering and conservation strategies based on RNA sequencing data alone, without using single-nucleotide polymorphism information or genome sequencing data from the same sample. We thereby unveil the first evolutionary landscape of A-to-I editing maps across 20 metazoan species (from worm to human), providing unprecedented evidence on how the editing mechanism gradually expands its territory and increases its influence along the history of evolution. Our result revealed that highly clustered and conserved editing sites tended to have a higher editing level and a higher magnitude of the ADAR motif. The ratio of the frequencies of nonsynonymous editing to that of synonymous editing remarkably increased with increasing the conservation level of A-to-I editing. These results thus suggest potentially functional benefit of highly clustered and conserved editing sites. In addition, spatiotemporal dynamics analyses reveal a conserved enrichment of editing and ADAR expression in the central nervous system throughout more than 300 Myr of divergent evolution in complex animals and the comparability of editing patterns between invertebrates and between vertebrates during development. This study provides evolutionary and dynamic aspects of A-to-I editome across metazoan species, expanding this important but understudied class of nongenomically encoded events for comprehensive characterization.

Integrative transcriptome sequencing identifies trans-splicing events with important roles in human embryonic stem cell pluripotency.

Trans-splicing is a post-transcriptional event that joins exons from separate pre-mRNAs. Detection of trans-splicing is usually severely hampered by experimental artifacts and genetic rearrangements. Here, we develop a new computational pipeline, TSscan, which integrates different types of high-throughput long-/short-read transcriptome sequencing of different human embryonic stem cell (hESC) lines to effectively minimize false positives while detecting trans-splicing. Combining TSscan screening with multiple experimental validation steps revealed that most chimeric RNA products were platform-dependent experimental artifacts of RNA sequencing. We successfully identified and confirmed four trans-spliced RNAs, including the first reported trans-spliced large intergenic noncoding RNA ("tsRMST"). We showed that these trans-spliced RNAs were all highly expressed in human pluripotent stem cells and differentially expressed during hESC differentiation. Our results further indicated that tsRMST can contribute to pluripotency maintenance of hESCs by suppressing lineage-specific gene expression through the recruitment of NANOG and the PRC2 complex factor, SUZ12. Taken together, our findings provide important insights into the role of trans-splicing in pluripotency maintenance of hESCs and help to facilitate future studies into trans-splicing, opening up this important but understudied class of post-transcriptional events for comprehensive characterization.

Evolutionary expansion of SPOP and associated TD/POZ gene family: impact of evolutionary route on gene expression pattern.

Evolutionary expansion of a gene family may occur at both the DNA and RNA levels. The rat testis-specific Rtdpoz-T2 and -T1 (rT2 and rT1) retrogenes are members of the TD/POZ gene family which also includes the well-characterized SPOP gene. In this study, rT2/rT1 transcriptional activation in cancer cells is demonstrated; the cancer rT2/rT1 transcripts are structurally similar to the embryonic transcripts reported previously in frequent exonization of transposed elements. On database interrogation, we have identified an uncharacterized rT2/rT1-like SPOP paralog, designated as SPOP-like (SPOPL), in the human and rodent genomes. Ka/Ks analysis indicates that the SPOPL genes are under functional constraints implicating biological functions. Phylogenetic analyses further suggest that segmental duplication and retrotransposition events had occurred giving rise to new gene members or retrogenes in the human-rodent ancestors during the evolution of the TD/POZ gene family. Based on this and previous works, a model is proposed to map the routes of evolutionary expansion of the TD/POZ gene family. More importantly, different gene expression patterns of members of the family are depicted: intron-harboring members are ubiquitously expressed whereas retrogenes are expressed in tissue-specific and developmentally regulated manner, and are fortuitously re-activated in cancer cells involving exonization of transposed elements.