It Is Imperative to Establish a Pellucid Definition of Chimeric RNA and to Clear Up a Lot of Confusion in the Relevant Research

Updating mRNA variants of the human RSK4 gene and their expression in different stressed situations

We determined RNA spectrum of the human RSK4 (hRSK4) gene (also called RPS6KA6) and identified 29 novel mRNA variants derived from alternative splicing, which, plus the NCBI-documented ones and the five we reported previously, totaled 50 hRSK4 RNAs that, by our bioinformatics analyses, encode 35 hRSK4 protein isoforms of 35-762 amino acids. Many of the mRNAs are bicistronic or tricistronic for hRSK4. The NCBI-normalized NM_014496.5 and the protein it encodes are designated herein as the Wt-1 mRNA and protein, respectively, whereas the NM_001330512.1 and the long protein it encodes are designated as the Wt-2 mRNA and protein, respectively. Many of the mRNA variants responded differently to different situations of stress, including serum starvation, a febrile temperature, treatment with ethanol or ethanol-extracted clove buds (an herbal medicine), whereas the same stressed situation often caused quite different alterations among different mRNA variants in different cell lines. Mosifloxacin, an antibiotics and also a functional inhibitor of hRSK4, could inhibit the expression of certain hRSK4 mRNA variants. The hRSK4 gene likely uses alternative splicing as a handy tool to adapt to different stressed situations, and the mRNA and protein multiplicities may partly explain the incongruous literature on its expression and comports.

Recent advances in cancer fusion transcript detection

Extensive investigation of gene fusions in cancer has led to the discovery of novel biomarkers and therapeutic targets. To date, most studies have neglected chromosomal rearrangement-independent fusion transcripts and complex fusion structures such as double or triple-hop fusions, and fusion-circRNAs. In this review, we untangle fusion-related terminology and propose a classification system involving both gene and transcript fusions. We highlight the importance of RNA-level fusions and how long-read sequencing approaches can improve detection and characterization. Moreover, we discuss novel bioinformatic tools to identify fusions in long-read sequencing data and strategies to experimentally validate and functionally characterize fusion transcripts.

RNA expression profiling from the liquid fraction of synovial fluid in knee joint osteoarthritis patients

OBJECTIVE

To investigate the RNA profile of synovial fluid (SF) from osteoarthritis (OA) patients and carry out cluster analysis of OA-related genes.

METHODS

RNA of SF from OA patients was isolated using RNA-specific Trizol. A cDNA library was built and subjected to the second-generation sequencing using HisSeq4000 with a data size of 8G. The sequencing reads were aligned to the UCSC human reference genome (hg38) using Tophat with default parameters. Gene function enrichment was generated using DAVID.

RESULTS

The minimum weight 0.096 µg RNA of SF sample was used for sequencing analysis, which produced 66,154,562 clean reads, 91.28% of which were matched to the reference with 2,682 genes identified. Some of the unmatchable reads matched RNAs of bacteria, mainly Pseudomonas. The detected human RNAs in samples fell into different categories of genes, including protein-coding ones, processed and unprocessed pseudogenes, and long noncoding, antisense and miscellaneous RNAs that mediate various biological functions. Interestingly, 80% of the expressed genes belonged to the mitochondrial genome.

CONCLUSION

These results suggest that less than 0.1 µg RNA is sufficient for establishing a cDNA library and deep sequencing, and that the liquid fraction of SF contains a whole RNA repertoire that may reflect a history of previous microorganism infections.

Establishment of New Genetic Markers and Methods for Sex Determination of Mouse and Human Cells using Polymerase Chain Reactions and Crude DNA Samples

Background: The currently available methods for sexing human or mouse cells have weaknesses. Therefore, it is necessary to establish new methods. Methods: We used bioinformatics approach to identify genes that have alleles on both the X and Y chromosomes of mouse and human genomes and have a region showing a significant difference between the X and Y alleles. We then used polymerase chain reactions (PCR) followed by visualization of the PCR amplicons in agarose gels to establish these genomic regions as genetic sex markers. Results: Our bioinformatics analyses identified eight mouse sex markers and 56 human sex markers that are new, i.e. are previously unreported. Six of the eight mouse markers and 14 of the 56 human markers were verified using PCR and ensuing visualization of the PCR amplicons in agarose gels. Most of the tested and untested sex markers possess significant differences in the molecular weight between the X- and Y-derived PCR amplicons and are thus much better than most, if not all, previously-reported genetic sex markers. We also established several simple and essentially cost-free methods for extraction of crude genomic DNA from cultured cells, blood samples, and tissues that could be used as template for PCR amplification. Conclusion: We have established new sex genetic markers and methods for extracting genomic DNA and for sexing human and mouse cells. Our work may also lend some methodological strategies to the identification of new genetic sex markers for other organismal species.

At elevated temperatures, heat shock protein genes show altered ratios of different RNAs and expression of new RNAs, including several novel HSPB1 mRNAs encoding HSP27 protein isoforms

Heat shock proteins (HSP) serve as chaperones to maintain the physiological conformation and function of numerous cellular proteins when the ambient temperature is increased. To determine how accurate the general assumption that HSP gene expression is increased in febrile situations is, the RNA levels of the HSF1 (heat shock transcription factor 1) gene and certain HSP genes were determined in three cell lines cultured at 37˚C or 39˚C for three days. At 39˚C, the expression of HSF1, HSPB1, HSP90AA1 and HSP70A1L genes demonstrated complex changes in the ratios of expression levels between different RNA variants of the same gene. Several older versions of the RNAs of certain HSP genes that have been replaced by a newer version in the National Center for Biotechnology Information database were also detected, indicating that the older versions are actually RNA variants of these genes. The present study cloned four new RNA variants of the HSP27-encoding HSPB1 gene, which together encode three short HSP27 peptides. Reanalysis of the proteomics data from our previous studies also demonstrated that proteins from certain HSP genes could be detected simultaneously at multiple positions using SDS-PAGE, suggesting that these genes may engender multiple protein isoforms. These results collectively suggested that, besides increasing their expression, certain HSP and associated genes also use alternative transcription start sites to produce multiple RNA transcripts and use alternative splicing of a transcript to produce multiple mature RNAs, as important mechanisms for responding to an increased ambient temperature in vitro.

Computational Detection of Known Pathogenic Gene Fusions in a Normal Tissue Database and Implications for Genetic Disease Research

Several recent studies have demonstrated the utility of RNA-Seq in the diagnosis of rare inherited disease. Diagnostic rates 35% higher than those previously achievable with DNA-Seq alone have been attained. These studies have primarily profiled gene expression and splicing defects, however, some have also shown that fusion transcripts are diagnostic or phenotypically relevant in patients with constitutional disorders. Fusion transcripts have traditionally been studied as oncogenic phenomena, with relevance only to cancer testing. Consequently, fusion detection algorithms were biased toward the detection of well-known oncogenic fusions, hindering their application to rare Mendelian genetic disease studies. A recent methodology published by the authors successfully tailored a traditional algorithm to the detection of pathogenic fusion events in inherited disease. A key mechanism of decreasing false positive or biologically benign events was comparison to a database of events detected in normal tissues. This approach is akin to population frequency-based filtering of genetic variants. It is predicated on the idea that pathogenic fusion transcripts are absent from normal tissue. We report on an analysis of RNA-Seq data from the genotype-tissue expression (GTEx) project in which known pathogenic fusions are computationally detected at low levels in normal tissues unassociated with the disease phenotype. Examples include archetypal cancer fusion transcripts, as well as fusions responsible for rare inherited disease. We consider potential explanations for the detectability of such transcripts and discuss the bearing such results have on the future profiling of genetic disease patients for pathogenic gene fusions.

About three-fourths of mouse proteins unexpectedly appear at a low position of SDS-PAGE, often as additional isoforms, questioning whether all protein isoforms have been eliminated in gene-knockout cells or organisms

Most genes in evolutionarily complex genomes are expressed to multiple protein isoforms, but there is not yet any simple high-throughput approach to identify these isoforms. Using an oversimplified top-down LC-MS/MS strategy, we detected, around the 26-kD position of SDS-PAGE, proteins produced from 782 genes in a Cdk4-/- mouse embryonic fibroblast cell line. Interestingly, only 213 (27.24%, about one-fourth) of these 782 genes have their proteins with a theoretical molecular mass (TMM) 10% smaller or larger than 26 kD, that is, between 23 and 29 kD, the range set as allowed variation in SDS-PAGE. These 213 proteins are considered as the wild type (WT). The remaining three-fourths includes proteins from 66 (9.44%) genes with a TMM smaller than 23 kD and proteins from 503 (64.32%, nearly two-thirds) genes with a TMM larger than 29 kD; these proteins are categorized into a larger-group or a smaller-group, respectively, for their appearance at a higher or lower position of SDS-PAGE. For instance, at this 26-kD position we detected proteins from the Rps27a, Snrpf, Hist1h4a, and Rps25 genes whose proteins' TMM is 8.6, 9.7, 11.4, and 13.7 kD, respectively, and detected proteins from the Plelc1 and Prkdc genes, whose largest isoform is 533.9 and 471.1 kD, respectively. We extrapolate that many of those proteins migrating unexpectedly in SDS-PAGE may be isoforms besides the WT protein. Moreover, we also detected a Cdk4 protein in this Cdk4-/- cell line, thus wondering whether some of other gene-knockout cells or organisms show similar incompleteness of the knockout.

Fusion Transcripts of Adjacent Genes: New Insights into the World of Human Complex Transcripts in Cancer

The awareness of genome complexity brought a radical approach to the study of transcriptome, opening eyes to single RNAs generated from two or more adjacent genes according to the present consensus. This kind of transcript was thought to originate only from chromosomal rearrangements, but the discovery of readthrough transcription opens the doors to a new world of fusion RNAs. In the last years many possible intergenic cis-splicing mechanisms have been proposed, unveiling the origins of transcripts that contain some exons of both the upstream and downstream genes. In some cases, alternative mechanisms, such as trans-splicing and transcriptional slippage, have been proposed. Five databases, containing validated and predicted Fusion Transcripts of Adjacent Genes (FuTAGs), are available for the scientific community. A comparative analysis revealed that two of them contain the majority of the results. A complete analysis of the more widely characterized FuTAGs is provided in this review, including their expression pattern in normal tissues and in cancer. Gene structure, intergenic splicing patterns and exon junction sequences have been determined and here reported for well-characterized FuTAGs. The available functional data and the possible roles in cancer progression are discussed.

A tailored approach to fusion transcript identification increases diagnosis of rare inherited disease

BACKGROUND

RNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data.

METHODS

We describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization.

RESULTS

We demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance.

CONCLUSIONS

The approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.

Transcriptional-Readthrough RNAs Reflect the Phenomenon of "A Gene Contains Gene(s)" or "Gene(s) within a Gene" in the Human Genome, and Thus Are Not Chimeric RNAs

Tens of thousands of chimeric RNAs, i.e., RNAs with sequences of two genes, have been identified in human cells. Most of them are formed by two neighboring genes on the same chromosome and are considered to be derived via transcriptional readthrough, but a true readthrough event still awaits more evidence and trans-splicing that joins two transcripts together remains as a possible mechanism. We regard those genomic loci that are transcriptionally read through as unannotated genes, because their transcriptional and posttranscriptional regulations are the same as those of already-annotated genes, including fusion genes formed due to genetic alterations. Therefore, readthrough RNAs and fusion-gene-derived RNAs are not chimeras. Only those two-gene RNAs formed at the RNA level, likely via trans-splicing, without corresponding genes as genomic parents, should be regarded as authentic chimeric RNAs. However, since in human cells, procedural and mechanistic details of trans-splicing have never been disclosed, we doubt the existence of trans-splicing. Therefore, there are probably no authentic chimeras in humans, after readthrough and fusion-gene derived RNAs are all put back into the group of ordinary RNAs. Therefore, it should be further determined whether in human cells all two-neighboring-gene RNAs are derived from transcriptional readthrough and whether trans-splicing truly exists.

While it is not deliberate, much of today's biomedical research contains logical and technical flaws, showing a need for corrective action

Biomedical research has advanced swiftly in recent decades, largely due to progress in biotechnology. However, this rapid spread of new, and not always-fully understood, technology has also created a lot of false or irreproducible data and artifacts, which sometimes have led to erroneous conclusions. When describing various scientific issues, scientists have developed a habit of saying "on one hand… but on the other hand…", because discrepant data and conclusions have become omnipresent. One reason for this problematic situation is that we are not always thoughtful enough in study design, and sometimes lack enough philosophical contemplation. Another major reason is that we are too rushed in introducing new technology into our research without assimilating technical details. In this essay, we provide examples in different research realms to justify our points. To help readers test their own weaknesses, we raise questions on technical details of RNA reverse transcription, polymerase chain reactions, western blotting and immunohistochemical staining, as these methods are basic and are the base for other modern biotechnologies. Hopefully, after contemplation and reflection on these questions, readers will agree that we indeed know too little about these basic techniques, especially about the artifacts they may create, and thus many conclusions drawn from the studies using those ever-more-sophisticated techniques may be even more problematic.

Chimeric RNAs and their implications in cancer

Chimeric RNAs have been believed to be solely produced by gene fusions resulting from chromosomal rearrangement, thus unique features of cancer. Detected chimeric RNAs have also been viewed as surrogates for the presence of gene fusions. However, more and more research has demonstrated that chimeric RNAs in general are not a hallmark of cancer, but rather widely present in non-cancerous cells and tissues. At the same time, they may be produced by other mechanisms other than chromosomal rearrangement. The field of non-canonical chimeric RNAs is still in its infancy, with many challenges ahead, including the lack of a unified terminology. However, we believe that these non-canonical chimeric RNAs will have significant impacts in cancer detection and treatment.

Probably less than one-tenth of the genes produce only the wild type protein without at least one additional protein isoform in some human cancer cell lines

To estimate how many genes produce multiple protein isoforms, we electrophoresed proteins from MCF7 and MDA-MB231 (MB231) human breast cancer cells in SDS-PAGE and excised narrow stripes of the gel at the 48kD, 55kD and 72kD. Proteins in these stripes were identified using liquid chromatography and tandem mass spectrometry. A total of 765, 750 and 679 proteins from MB231 cells, as well as 470, 390 and 490 proteins from MCF7 cells, were identified from the 48kD, 55kD and 72kD stripes, respectively. We arbitrarily allowed a 10% technical variation from the proteins' theoretical molecular mass (TMM) and considered those proteins with their TMMs within the 43-53 kD, 49-61 kD and 65-79 kD ranges as the wild type (WT) expected from the corresponding stripe, whereas those with a TMM above or below this range as a smaller- or larger-group, respectively. Only 263 (34.4%), 269 (35.9%) and 151 (22.2%) proteins from MB231 cells and 117 (24.9%), 135 (34.6%) and 130 (26.5%) proteins from MCF7 cells from the 48kD, 55kD and 72kD stripes, respectively, belonged to the WT, while the remaining majority belonged to the smaller- or larger-groups. Only about 3-16%, on average about 10% regardless of the stripe and cell line, of the proteins appeared in only one stripe and within the WT range, while the remaining preponderance appeared also in additional stripe(s) or had a larger or smaller TMM. We conclude that few (fewer than 10%) of the human genes produce only the WT protein without additional isoform(s).