<figure><img src="image_1.png" /></figure>

Check for
updates

# Systematic Nucleotide Exchange Analysis of ESTs From the Human Cancer Genome Project Report: Origins of 347 Unknown ESTs Indicate Putative Transcription of Non-Coding Genomic Regions

OPEN ACCESS

**Edited by:**

Dapeng Wang,
University of Leeds,
United Kingdom

**Reviewed by:**

Rachita Yadav,
Harvard Medical School,
United States
Hauke Busch,
Universität zu Lübeck,
Germany

*Correspondence:*

Ganesh Warthi
g.warthi6791@gmail.com

**Specialty section:**

This article was submitted to
Bioinformatics and Computational
Biology, a section of the journal
Frontiers in Genetics

**Received:** 28 October 2019

**Accepted:** 15 January 2020

**Published:** 11 February 2020

**Citation:**

Warthi G, Fournier P-E and Seligmann H (2020) Systematic Nucleotide Exchange Analysis of ESTs From the Human Cancer Genome Project Report: Origins of 347 Unknown ESTs Indicate Putative Transcription of Non-Coding Genomic Regions. *Front. Genet.* 11:42.
doi: 10.3389/fgene.2020.00042

Ganesh Warthi<fnref n="1" />,²*, Pierre-Edouard Fournier<fnref n="1" />,² and Hervé Seligmann³,<fnref n="4" />

Expressed sequence tags (ESTs) provide an imprint of cellular RNA diversity irrespectively of sequence homology with template genomes. NCBI databases include many unknown RNAs from various normal and cancer cells. These are usually ignored assuming sequencing artefacts or contamination due to their lack of sequence homology with template DNA. Here, we report genomic origins of 347 ESTs previously assumed artefacts/unknown, from the FAPESP/LICR Human Cancer Genome Project. EST template detection uses systematic nucleotide exchange analyses called swinger transformations. Systematic nucleotide exchanges replace systematically particular nucleotides with different nucleotides. Among 347 unknown ESTs, 51 ESTs match mitogenome transcription, 17 and 2 ESTs are from nuclear chromosome non-coding regions, and uncharacterized nuclear genes. Identified ESTs mapped on 205 protein-coding genes, 10 genes had swinger RNAs in several biosamples. Whole cell transcriptome searches for 17 ESTs mapping on non-coding regions confirmed their transcription. The 10 swinger-transcribed genes identified more than once associate with cancer induction and progression, suggesting swinger transformation occurs mainly in highly transcribed genes. Swinger transformation is a unique method to identify noncanonical RNAs obtained from NGS, which identifies putative ncRNA transcribed regions. Results suggest that swinger transcription occurs in highly active genes in normal and genetically unstable cancer cells.

**Keywords:** swinger RNAs, systematic nucleotide exchanges, transcriptome, cDNA, cancer

<footnote n="1">¹ Aix Marseille Univ, IRD, APHM, SSA, VITROME, IHU-Méditerranée Infection, Marseille, France, ² IHU-Méditerranée Infection, Marseille, France, ³ The National Natural History Collections, The Hebrew University of Jerusalem, Jerusalem, Israel,</footnote>

<footnote n="4">⁴ Université Grenoble Alpes, Faculty of Medicine, Laboratory AGEIS EA 7407, Team Tools for e-Gnosis Medical & Labcom CNRS/UGA/OrangeLabs Telecoms4Health, La Tronche, France</footnote>
INTRODUCTION

Systematic nucleotide exchanges, also called swinger
transformations, systematically exchange specific nucleotides
with other specific nucleotides during DNA replication and/or
RNA transcription. DNA or RNA molecules corresponding to
systematic nucleotide changes are called swinger DNAs or swinger
RNAs. Only 23 systematic nucleotide exchanges (Figure 1) are
possible with four nucleotide bases (A, T, C and G), i.e. nine
symmetric (X ↔ Y, e.g. A ↔ G) (Seligmann, 2013a; Seligmann,
2013b) and 14 asymmetric (X → Y → Z → X, e.g. A → G → T →
A) (Seligmann, 2013b; Seligmann, 2013c). For example, in
symmetric exchange A ↔ G, all As are replaced by Gs and all
Gs by As. The 14 asymmetric exchanges are directional, e.g. A →
G → T → A: all As are replaced by Gs, Gs by Ts and Ts by As. It is
unclear whether these systematic exchanges occur during DNA/
RNA polymerizations or result from posttranscriptional editions.
Their relatively long lengths (> 100 nucleotides) favors the former.
Previous correlation analyses (Seligmann, 2013b; Seligmann,
2013c; Michel and Seligmann, 2014) show that the lengths and
abundances of swinger RNAs are approximately proportional to
rates calculated on the basis of corresponding single nucleotide
misinsertions by the human mitochondrial gamma DNA
polymerase (from Lee and Johnson, 2006). This suggests that
swinger RNAs result from polymerizations where the polymerase
is stabilized in the usually transient, unstable state that causes
regular single nucleotide misinsertions.

Systematic nucleotide exchange analyses identify genomic
origins of unknown DNA/RNA sequences. Several types of
RNAs with no sequence homology to template genomes exist.
This includes RNA-DNA differences (Blank et al., 1986; Ninio,
1991; Li et al., 2011; Bahn et al., 2012; Strathern et al., 2012;
Bar-Yaacov et al., 2013; Knippa and Peterson, 2013; Zhou

et al., 2013; Wang et al., 2016), post-transcriptional editing (Bass, 2002; Schaub and Keller, 2002; Chen and Bundschuh, 2012; Park et al., 2012; Wang IX et al., 2014a; Lee et al., 2015), post-transcriptional hyper editing (Porath et al., 2014) and polymerase template switching (Lee et al., 2007; Löytynoja and Goldman, 2017). Transcript fusion explains some noncanonical RNAs (Kumar et al., 2016; López-Nieva et al., 2019; Singh et al., 2019). Some rare DNAs and RNAs seem to result from spontaneous, template-free polymerization (Béguin et al., 2015; Seligmann and Raoult, 2018). Due to the systematic nature of swinger transformations, identified, previously unknown sequence reads overcome the widely argued possibility of by-chance alignment with randomly amplified sequencing artefacts since randomness does not produce systematic exchanges.

Mitogenomic human swinger RNA reads produced by
Illumina confirmed corresponding EST data (Seligmann,
2016b). The swinger transcriptome of the amoeban-hosted
Mimivirus was also confirmed by two different sequencing
techniques (SOLID and 454, Seligmann and Raoult, 2018).
Detected peptides matching translation of the swinger-
transformed mitogenome tend to map on detected swinger
RNAs (Seligmann, 2016a; Seligmann, 2016b; Seligmann, 2016c;
Seligmann, 2017a). These findings were further confirmed in
purified mitochondrial transcriptomes (Warthi and Seligmann,
2018), rejecting the possibility of cytoplasmic contamination.
Chimeric mitochondrial swinger RNAs also exist, partly
following swinger polymerization and partly regular
polymerization, with abrupt switches between these parts
(Seligmann, 2015a; Seligmann, 2015b). This observation is
paralleled by chimeric peptides, corresponding to translation of
adjacent regular and swinger RNA (Seligmann, 2016d). Swinger
RNA coverage associates with secondary structure formation

Systematic nucleotide exchanges

A Symmetric

1. A → C

2. A → G

3. A → T

4. C → G

5. C → T

6. G → T

7. A → C + G → T

8. A → G + C → T

9. A → T + C → G

E.g.: A↔C

ATGTACTGC

<figure><img src="image_2.png" /></figure>

CTGTCATGA

B Asymmetric

1. A → C → T → A

2. A → C → G → A

3. A → G → C → A

4. A → G → T → A

5. A → T → C → A

6. A → T → G → A

7. C → G → T → C

8. C → T → G → C

9. A → C → G → T → A

10. A → C → T → G → A

11. A → G → C → T → A

12. A → G → T → C → A

13. A → T → C → G → A

14. A → T → G → C → A

E.g. : A → C → T → A

ATGTACTGC

<figure><img src="image_3.png" /></figure>

CAGACTAGT

**FIGURE 1** | Twenty three systematic nucleotide exchanges. (A) 9 Symmetric exchanges in which one nucleotide is exchanged by another during DNA replication or transcription, e.g. A ↔ C, where all As are replaced by Cs and all Cs by As. (B) 14 Asymmetric exchanges that are directional in nature, e.g. A→C→T→A: all As are replaced by Cs, Cs by Ts and Ts by As.
(Seligmann, 2016e). Swinger DNA also occurs (Seligmann, 2014a; Seligmann, 2014b). Molecular functions and associations of swinger polymerizations with healthy or unhealthy cells remain unknown. Recently, swinger RNAs with $A \leftrightarrow T + C \leftrightarrow G$ transformation were identified in HIV mediated Non-Hodgkin's Lymphoma (Warthi et al., 2019). The identified RNAs with $A \leftrightarrow T + C \leftrightarrow G$ transformation are similar to the transcriptional product by polymerase template switching, known in retroviruses (Kandel and Nudler, 2002) and the human genome (Löytynoja and Goldman, 2017).

The Encyclopedia of DNA Elements (ENCODE) project concluded that a human genome includes approximately 20,000 protein-coding genes and predicts that 80% of non-coding regions regulate gene expression (Dunham et al., 2012). However, the transcriptomes generated by current RNA sequencing technologies do not cover the complete human genome, because analyses of sequencing reads assume only canonical transcription. Analyses assuming systematic nucleotide exchanges could identify non-canonical RNAs and their genomic origin.

Expressed sequence tags (ESTs) are short DNA sequences (100–600 bp) generated from the sequencing of cDNA libraries. The ESTs represent RNAs derived from a particular cell irrespective of sequence similarity or dissimilarity to its template genome, unlike transcriptomes, in which sequencing reads with sequence similarity are considered as true reads while ignoring non-homologous reads. This underlines a major limitation in identifying and studying non-canonical RNAs and their association with various genetic diseases like cancers. Therefore, here we apply swinger transformations to identify unknown expressed sequence tags (ESTs) occurring in the FAPESP/LICR Human Cancer Genome Project (Neto et al., 2000). We expect to identify unknown ESTs reported in cancer cells, identify and report their genomic origins and confirm that swinger-transformed RNAs occur beyond mitochondria.

# METHODS

## Identification of Unknown Ests

The FAPESP/LICR Human Cancer Genome Project corresponds to 891,011 published (Neto et al., 2000) and 55,248 unpublished ESTs in the NCBI database. These ESTs were blasted with the “Human genomic + transcript” database using the “highly similar sequence” algorithm (megablast) (Zhang et al, 2000; Morgulis et al., 2008). A total of 149,500 ESTs with no sequence similarity with the human genome were further blasted across all sequences in the NCBI database to exclude ESTs corresponding to potential contaminations.

## Swingier Transformation of Unknown Ests

The unknown 149,500 ESTs were then swinger-transformed according to the 23 bijective transformations (Figure 1) (Michel and Seligmann, 2014). These 149,500 unknown ESTs were blasted again with the “Human genomic + transcript” database using megablast, in order to detect their genomic template.

# RESULTS AND DISCUSSION

In total, 347 ESTs (0.23%) were identified using swinger transformation analyses (Tables 1–3). This rate is 4x lower than that estimated for human mitochondria (about 100/10000 ESTs, 1% (data from Seligmann, 2012; Seligmann, 2013b; Seligmann, 2013c). Swinger-transformed sequences of these identified 347 ESTs are in Supplementary Table 1. The remaining 99.77% of unidentified ESTs could be sequencing artefacts. These could also be ESTs with non-systematic post-transcriptional hyper editing (Porath et al., 2014) or resulting from systematic nucleotide deletion transcription (Warthi and Seligmann, 2019). Among 347 ESTs, 223 ESTs were symmetrically transformed swinger RNAs (219 ESTs were $A \leftrightarrow T$ transformed, two ESTs were $A \leftrightarrow C$ transformed, and one of each was $C \leftrightarrow T$ and $G \leftrightarrow T$ transformed). The average length of alignment of $A \leftrightarrow T$, $A \leftrightarrow C$, $C \leftrightarrow T$ and $G \leftrightarrow T$ transformed ESTs was 325 bp, 225 bp, 196 bp, 152 bp with an average percentage identity of 96.95, 90.52, 96.42, and 89.47% respectively. Similarly, $A \rightarrow T \rightarrow C \rightarrow G \rightarrow A$ transformed asymmetric swinger transformations account for 124 ESTs with 96.95% average percentage identity and 308 bp average alignment length. The largest aligned swinger RNA was 685 bp with 41 bp as the smallest. Unlike symmetrical transformations, ESTs identified as $A \rightarrow T \rightarrow C \rightarrow G \rightarrow A$ transformed have to be transcribed assuming $A \rightarrow G \rightarrow C \rightarrow T \rightarrow A$ from the template DNA region (Figure 2). From here on, these identified ESTs are considered swinger RNAs. Abundances of swinger RNA classes are proportional to their mean length ($r = 0.96$, two tailed $P = 0.00945$). Assuming that different swinger classes result from different polymerase states, the positive correlation indicates that the same factor that promotes switching to a given mode of swinger transcription also favors the stability of this mode. The opposite, indicated by a negative association, would mean that frequent types of switches are unstable. We suggest that systematic nucleotide transformations happen during transcription (swinger transcription). This might result from polymerase enzyme fatigue. In carcinogenesis, during malignant transformation, cancer cells produce mutated RNAs and proteins because of genetic instability (Ciocca and Calderwood, 2005; Ruckova et al., 2012). Therefore, swinger RNAs are probably mainly non- or dys-functional, and produce non- and dysfunctional proteins with probable carcinogenic effects.

## Swingier-Transformed Genes in Different Cancer Tissues

These 347 swinger-transformed RNAs matched 205 known and two uncharacterized genes (Tables 1–3). These identified swinger RNAs might be artefacts. However, among these 207 genes, swinger RNAs from 10 genes were detected in multiple cancer types and biosamples (Table 2). These genes are adenylate kinase 1 (AK1), cystatin B (CSTB), DLG associated protein (DLGAP), heat shock protein family A (Hsp70) member 8 (HSPA8), mediator complex subunits (MED), poly(A) binding protein cytoplasmic 1 (PABPC1), MT-16s rRNA, MT-ND1, MT-ND2 and MT-ND4.
<table><caption>TABLE 1 | List of swinger transformed ESTs identified using BLAST. The list includes ESTs identified using swinger transformations and mapped on genes identified only in one biosamples.</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AI902333.1</td><td>SAMN00156544</td><td>Breast</td><td>173</td><td>A &rarr; C</td><td>NM_001278462.1</td><td>12</td><td>MDM2 proto-oncogene (MDM2)</td><td>100</td><td>97.000</td><td>2.37e-39</td><td>65</td><td>163</td></tr><tr><td>AI903883.1</td><td>SAMN00156581</td><td>Breast</td><td>193</td><td>G &rarr; T</td><td>NM_080386.3</td><td>2</td><td>tubulin alpha 3d (TUBA3D)</td><td>147</td><td>93.197</td><td>4.39e-57</td><td>1,095</td><td>950</td></tr><tr><td>AI939933.1</td><td>SAMN00156800</td><td>Colon</td><td>283</td><td>A &rarr; T</td><td>NM_001185092.1</td><td>1</td><td>nitrilase 1 (NIT1)</td><td>168</td><td>94.643</td><td>6.67e-67</td><td>402</td><td>238</td></tr><tr><td>AI939934.1*</td><td>SAMN00156800</td><td>Colon</td><td>273</td><td>A &rarr; T</td><td>NM_001037333.2</td><td>5</td><td>cytoplasmic FMR1 interacting protein 2 (CYFIP2)</td><td>200</td><td>96.500</td><td>1.80e-87</td><td>738</td><td>540</td></tr><tr><td>AI939935.1*</td><td>SAMN00156800</td><td>Colon</td><td>156</td><td>A &rarr; T</td><td>NM_001037333.2</td><td>5</td><td>cytoplasmic FMR1 interacting protein 2 (CYFIP2)</td><td>142</td><td>95.775</td><td>1.45e-58</td><td>733</td><td>592</td></tr><tr><td>AI939936.1*</td><td>SAMN00156800</td><td>Colon</td><td>303</td><td>A &rarr; T</td><td>NM_001037333.2</td><td>5</td><td>cytoplasmic FMR1 interacting protein 2 (CYFIP2)</td><td>202</td><td>98.515</td><td>1.78e-97</td><td>741</td><td>540</td></tr><tr><td>AI939937.1</td><td>SAMN00156800</td><td>Colon</td><td>253</td><td>A &rarr; T</td><td>NM_001008491.2</td><td>2</td><td>septin 2 (SEPT2)</td><td>234</td><td>98.291</td><td>6.26e-112</td><td>2,935</td><td>2,703</td></tr><tr><td>AI939938.1*</td><td>SAMN00156800</td><td>Colon</td><td>266</td><td>A &rarr; T</td><td>NM_001037333.2</td><td>5</td><td>cytoplasmic FMR1 interacting protein 2 (CYFIP2)</td><td>202</td><td>94.059</td><td>8.44e-81</td><td>541</td><td>741</td></tr><tr><td>AI939939.1</td><td>SAMN00156800</td><td>Colon</td><td>377</td><td>A &rarr; T</td><td>NW_019805496.1</td><td>11</td><td>NADH:ubiquinone oxidoreductase core subunit S3 (NDUFS3)</td><td>246</td><td>98.374</td><td>1.34e-118</td><td>10,294</td><td>10,049</td></tr><tr><td>AI939940.1</td><td>SAMN00156800</td><td>Colon</td><td>361</td><td>A &rarr; T</td><td>NM_001008491.2</td><td>2</td><td>septin 2 (SEPT2)</td><td>299</td><td>94.314</td><td>2.85e-125</td><td>2,905</td><td>2,608</td></tr><tr><td>AI939941.1</td><td>SAMN00156801</td><td>Colon</td><td>303</td><td>A &rarr; T</td><td>NM_001267774.1</td><td>17</td><td>intraflagellar transport 20 (IFT20)</td><td>220</td><td>98.182</td><td>1.05e-104</td><td>609</td><td>390</td></tr><tr><td>AI939942.1</td><td>SAMN00156801</td><td>Colon</td><td>182</td><td>A &rarr; T</td><td>NC_000008.11</td><td>8</td><td>vir like m6A methyltransferase associated (VIRMA)</td><td>77</td><td>97.403</td><td>5.45e-28</td><td>94,525,717</td><td>94,525,793</td></tr><tr><td>AI939943.1</td><td>SAMN00156801</td><td>Colon</td><td>312</td><td>A &rarr; T</td><td>NM_001282538.1</td><td>1</td><td>G protein subunit beta 1 (GNB1)</td><td>301</td><td>96.013</td><td>1.01e-134</td><td>1,212</td><td>1,512</td></tr><tr><td>AI939947.1*</td><td>SAMN00156461</td><td>Colon</td><td>210</td><td>A &rarr; T</td><td>NM_001009570.2</td><td>2</td><td>chaperonin containing TCP1 subunit 7 (CCT7)</td><td>182</td><td>95.604</td><td>5.08e-78</td><td>1,127</td><td>946</td></tr><tr><td>AI939949.1*</td><td>SAMN00156461</td><td>Colon</td><td>622</td><td>A &rarr; T</td><td>NM_001009570.2</td><td>2</td><td>chaperonin containing TCP1 subunit 7 (CCT7)</td><td>572</td><td>94.231</td><td>0.0</td><td>576</td><td>1,143</td></tr><tr><td>AI939950.1*</td><td>SAMN00156461</td><td>Colon</td><td>512</td><td>A &rarr; T</td><td>NM_001009570.2</td><td>2</td><td>chaperonin containing TCP1 subunit 7 (CCT7)</td><td>484</td><td>98.347</td><td>0.0</td><td>1,183</td><td>702</td></tr><tr><td>AI939951.1*</td><td>SAMN00156802</td><td>Colon</td><td>317</td><td>A &rarr; T</td><td>NR_037849.1</td><td>2</td><td>WW domain binding protein 1 (WBP1) INO80B-WBP1 readthrough (NMD candidate) (INO80B-WBP1), long non-coding RNA</td><td>310</td><td>94.839</td><td>4.70e-133</td><td>1,103</td><td>796</td></tr><tr><td>AI939953.1</td><td>SAMN00156802</td><td>Colon</td><td>609</td><td>A &rarr; T</td><td>NM_001134671.2</td><td>8</td><td>derlin 1 (DERL1)</td><td>608</td><td>97.862</td><td>0.0</td><td>2,407</td><td>1,802</td></tr><tr><td>AI939959.1*</td><td>SAMN00156802</td><td>Colon</td><td>218</td><td>A &rarr; T</td><td>NR_037849.1</td><td>2</td><td>WW domain binding protein 1 (WBP1) INO80B-WBP1 readthrough (NMD candidate) (INO80B-WBP1), long non-coding RNA</td><td>219</td><td>97.717</td><td>6.34e-102</td><td>1,043</td><td>826</td></tr><tr><td>AI939961.1*</td><td>SAMN00156802</td><td>Colon</td><td>231</td><td>A &rarr; T</td><td>NM_033017.3</td><td>7</td><td>tripartite motif containing 4 (TRIM4)</td><td>204</td><td>95.588</td><td>2.33e-86</td><td>2,432</td><td>2,633</td></tr><tr><td>AI939962.1*</td><td>SAMN00156802</td><td>Colon</td><td>212</td><td>A &rarr; T</td><td>NM_001300821.2</td><td>12</td><td>eukaryotic translation initiation factor 4B (EIF4B)</td><td>160</td><td>96.250</td><td>6.67e-67</td><td>2,434</td><td>2,276</td></tr><tr><td>AI939963.1*</td><td>SAMN00156802</td><td>Colon</td><td>340</td><td>A &rarr; T</td><td>NM_033017.3</td><td>7</td><td>tripartite motif containing 4 (TRIM4)</td><td>291</td><td>95.189</td><td>2.85e-125</td><td>2,629</td><td>2,344</td></tr><tr><td>AI939965.1</td><td>SAMN00156802</td><td>Colon</td><td>214</td><td>A &rarr; T</td><td>NR_038854.1</td><td>12</td><td>negative regulator of antiviral response (non-protein coding) (NRAV), long non-coding RNA</td><td>220</td><td>94.545</td><td>1.39e-88</td><td>632</td><td>848</td></tr><tr><td>AI939966.1*</td><td>SAMN00156802</td><td>Colon</td><td>366</td><td>A &rarr; T</td><td>NR_037849.1</td><td>2</td><td>WW domain binding protein 1 (WBP1) INO80B-WBP1 readthrough (NMD candidate) (INO80B-WBP1), long non-coding RNA</td><td>317</td><td>94.637</td><td>2.81e-135</td><td>1,717</td><td>2,032</td></tr><tr><td>AI939968.1</td><td>SAMN00156802</td><td>Colon</td><td>270</td><td>A &rarr; T</td><td>NC_000003.12</td><td>3</td><td>CUB domain containing protein 1 (CDCP1)</td><td>270</td><td>98.889</td><td>1.31e-133</td><td>45,123,674</td><td>45,123,943</td></tr><tr><td>AI939969.1*</td><td>SAMN00156802</td><td>Colon</td><td>143</td><td>A &rarr; T</td><td>NM_001300821.2</td><td>12</td><td>eukaryotic translation initiation factor 4B (EIF4B)</td><td>123</td><td>98.374</td><td>6.63e-54</td><td>2,313</td><td>2,434</td></tr><tr><td>AI939970.1*</td><td>SAMN00156802</td><td>Colon</td><td>255</td><td>A &rarr; T</td><td>NM_033017.3</td><td>7</td><td>tripartite motif containing 4 (TRIM4)</td><td>231</td><td>94.805</td><td>1.04e-96</td><td>2,406</td><td>2,633</td></tr><tr><td>AI939972.1</td><td>SAMN00156802</td><td>Colon</td><td>168</td><td>A &rarr; T</td><td>NC_000017.11</td><td>17</td><td>target of myb1 like 2 membrane trafficking protein (TOM1L2)</td><td>166</td><td>96.386</td><td>1.39e-70</td><td>17,849,628</td><td>17,849,465</td></tr><tr><td>AI939973.1</td><td>SAMN00156802</td><td>Colon</td><td>337</td><td>A &rarr; T</td><td>NR_037616.1</td><td>6</td><td>thioredoxin domain containing 5 (TXNDC5) and BLOC1S5-TXNDC5 readthrough (NMD candidate) (BLOC1S5-TXNDC5), long non-coding RNA</td><td>340</td><td>97.059</td><td>4.42e-160</td><td>1,880</td><td>1,542</td></tr><tr><td>AI939976.1*</td><td>SAMN00156802</td><td>Colon</td><td>200</td><td>A &rarr; T</td><td>NR_037849.1</td><td>2</td><td>WW domain binding protein 1 (WBP1) INO80B-WBP1 readthrough (NMD candidate) (INO80B-WBP1), long non-coding RNA</td><td>200</td><td>100.000</td><td>1.34e-100</td><td>1,785</td><td>1,984</td></tr><tr><td>AI939980.1</td><td>SAMN00156802</td><td>Colon</td><td>155</td><td>A &rarr; T</td><td>NM_144570.2</td><td>16</td><td>Jupiter microtubule associated homolog 2 (JPT2)</td><td>156</td><td>94.872</td><td>3.04e-62</td><td>2,949</td><td>2,794</td></tr><tr><td>AI939981.1*</td><td>SAMN00156802</td><td>Colon</td><td>229</td><td>A &rarr; T</td><td>NR_037849.1</td><td>2</td><td>WW domain binding protein 1 (WBP1) INO80B-WBP1 readthrough (NMD candidate) (INO80B-WBP1), long non-coding RNA</td><td>235</td><td>96.170</td><td>2.22e-103</td><td>1,996</td><td>1,762</td></tr><tr><td>AI939988.1*</td><td>SAMN00156802</td><td>Colon</td><td>290</td><td>A &rarr; T</td><td>NR_037849.1</td><td>2</td><td>WW domain binding protein 1 (WBP1) INO80B-WBP1 read through (NMD candidate) (INO80B-WBP1), long non-coding RNA</td><td>280</td><td>96.786</td><td>7.71e-128</td><td>2,021</td><td>1,744</td></tr><tr><td>AI940584.1*</td><td>SAMN00156483</td><td>Head_neck</td><td>583</td><td>A &rarr; T</td><td>NM_001166285.1</td><td>2</td><td>chaperonin containing TCP1 subunit 7 (CCT7)</td><td>568</td><td>93.838</td><td>0.0</td><td>1,890</td><td>1,326</td></tr><tr><td>AW062762.1</td><td>SAMN00156993</td><td>Colon</td><td>433</td><td>A &rarr; T</td><td>NM_001354666.1</td><td>3</td><td>peroxisome proliferator activated receptor gamma (PPARG), transcript</td><td>409</td><td>98.778</td><td>0.0</td><td>1,137</td><td>1,545</td></tr><tr><td>AW062764.1</td><td>SAMN00156993</td><td>Colon</td><td>610</td><td>A &rarr; T</td><td>NC_000002.12</td><td>2</td><td>protein kinase C epsilon (PRKCE)</td><td>491</td><td>99.389</td><td>0.0</td><td>45,959,855</td><td>45,959,366</td></tr><tr><td>AW062765.1</td><td>SAMN00156993</td><td>Colon</td><td>360</td><td>A &rarr; T</td><td>NC_000009.12</td><td>9</td><td>sorting nexin family member 30 (SNX30)</td><td>341</td><td>99.120</td><td>2.01e-173</td><td>112,811,006</td><td>112,811,345</td></tr><tr><td>AW062769.1</td><td>SAMN00156993</td><td>Colon</td><td>595</td><td>A &rarr; T</td><td>NM_001198773.2</td><td>1</td><td>phosphatidylinositol 4-kinase beta (PI4KB)</td><td>466</td><td>98.712</td><td>0.0</td><td>3,151</td><td>3,615</td></tr><tr><td>AW062770.1</td><td>SAMN00156993</td><td>Colon</td><td>589</td><td>A &rarr; T</td><td>NM_000016.10</td><td>16</td><td>Rho GTPase activating protein 17 (ARHGAP17)</td><td>586</td><td>99.147</td><td>0.0</td><td>24,921,862</td><td>24,922,445</td></tr><tr><td>AW062771.1</td><td>SAMN00156993</td><td>Colon</td><td>643</td><td>A &rarr; T</td><td>NM_006386.4</td><td>22</td><td>DEAD-box helicase 17 (DDX17)</td><td>605</td><td>96.694</td><td>0.0</td><td>1,878</td><td>2,480</td></tr><tr><td>AW062772.1</td><td>SAMN00156993</td><td>Colon</td><td>540</td><td>A &rarr; T</td><td>NM_000381.3</td><td>X</td><td>midline 1 (MID1)</td><td>519</td><td>99.037</td><td>0.0</td><td>5,758</td><td>5,240</td></tr><tr><td>AW062773.1</td><td>SAMN00156993</td><td>Colon</td><td>378</td><td>A &rarr; T</td><td>NM_002795.3</td><td>17</td><td>proteasome subunit beta 3 (PSMB3)</td><td>368</td><td>98.913</td><td>0.0</td><td>384</td><td>751</td></tr><tr><td>AW062775.1</td><td>SAMN00156993</td><td>Colon</td><td>650</td><td>A &rarr; T</td><td>NM_000358.2</td><td>5</td><td>transforming growth factor beta induced (TGFB1)</td><td>325</td><td>97.846</td><td>2.05e-158</td><td>1,166</td><td>1,489</td></tr><tr><td>AW062776.1</td><td>SAMN00156993</td><td>Colon</td><td>722</td><td>A &rarr; T</td><td>XM_017008304.2</td><td>4</td><td>PREDICTED: Homo sapiens RNA binding motif protein 47 (RBM47)</td><td>357</td><td>97.479</td><td>9.36e-172</td><td>1,924</td><td>2,280</td></tr><tr><td>AW062777.1</td><td>SAMN00156993</td><td>Colon</td><td>449</td><td>A &rarr; T</td><td>NM_001345950.1</td><td>10</td><td>tripartite motif containing 8 (TRIM8)</td><td>110</td><td>92.727</td><td>4.07e-36</td><td>1,009</td><td>901</td></tr><tr><td>AW062778.1</td><td>SAMN00156993</td><td>Colon</td><td>491</td><td>A &rarr; T</td><td>NM_000010.11</td><td>10</td><td>C-terminal binding protein 2 (CTBP2)</td><td>474</td><td>98.734</td><td>0.0</td><td>125,121,894</td><td>125,121,421</td></tr><tr><td>AW062779.1</td><td>SAMN00156993</td><td>Colon</td><td>252</td><td>A &rarr; T</td><td>NM_001136153.1</td><td>6</td><td>activating transcription factor 6 beta (ATF6B)</td><td>234</td><td>98.718</td><td>6.09e-114</td><td>1,207</td><td>1,440</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AW062782.1</td><td>SAMN00156993</td><td>Colon</td><td>558</td><td>A ↔ T</td><td>NC_000008.11</td><td>8</td><td>ArfGAP with SH3 domain, ankyrin repeat and PH domain 1 (ASAP1)</td><td>356</td><td>95.787</td><td>3.41e-161</td><td>130,131,998</td><td>130,132,351</td></tr><tr><td>AW062784.1</td><td>SAMN00156993</td><td>Colon</td><td>542</td><td>A ↔ T</td><td>NM_015397.3</td><td>9</td><td>DDB1 and CUL4 associated factor 12 (DCAF12)</td><td>326</td><td>97.546</td><td>4.45e-155</td><td>2,795</td><td>3,116</td></tr><tr><td>AW062787.1</td><td>SAMN00156993</td><td>Colon</td><td>647</td><td>A ↔ T</td><td>NC_000012.12</td><td>12</td><td>signal transducer and activator of transcription 2</td><td>608</td><td>91.776</td><td>0.0</td><td>56,349,156</td><td>56,348,551</td></tr><tr><td>AW062788.1</td><td>SAMN00156993</td><td>Colon</td><td>458</td><td>A ↔ T</td><td>NC_000021.9</td><td>21</td><td>VPS26 endosomal protein sorting factor C (VPS26C)</td><td>293</td><td>97.611</td><td>4.54e-140</td><td>37,236,577</td><td>37,236,868</td></tr><tr><td>AW062789.1</td><td>SAMN00156993</td><td>Colon</td><td>685</td><td>A ↔ T</td><td>NC_000017.11</td><td>17</td><td>Myosin XVIIA (MYO18A)</td><td>406</td><td>98.276</td><td>0.0</td><td>29,141,358</td><td>29,141,763</td></tr><tr><td>AW062792.1</td><td>SAMN00156993</td><td>Colon</td><td>688</td><td>A ↔ T</td><td>NM_001082480.2</td><td>8</td><td>zinc finger protein 623 (ZNF623)</td><td>671</td><td>97.914</td><td>0.0</td><td>3,041</td><td>3,707</td></tr><tr><td>AW062793.1</td><td>SAMN00156993</td><td>Colon</td><td>283</td><td>A ↔ T</td><td>NC_000007.14</td><td>7</td><td>PREDICTED: Homo sapiens speedy/RINGO cell cycle regulator family member E2B (SPDYE2B)</td><td>206</td><td>84.466</td><td>1.86e-49</td><td>102,357,109</td><td>102,357,313</td></tr><tr><td>AW062794.1</td><td>SAMN00156993</td><td>Colon</td><td>115</td><td>A ↔ T</td><td>NM_001282507.1</td><td>16</td><td>PKD1P4-NPIPA8 readthrough (PKD1P4-NPIPA8), non-coding RNA+ nuclear pore complex interacting protein family member A7 (NPIPA7)</td><td>98</td><td>96.939</td><td>2.43e-38</td><td>778</td><td>875</td></tr><tr><td>AW062796.1</td><td>SAMN00156993</td><td>Colon</td><td>196</td><td>A ↔ T</td><td>NM_001256410.1</td><td>10</td><td>RAB18, member RAS oncogene family (RAB18)</td><td>183</td><td>98.361</td><td>1.37e-85</td><td>876</td><td>1,058</td></tr><tr><td>AW062799.1</td><td>SAMN00156993</td><td>Colon</td><td>316</td><td>A ↔ T</td><td>NM_022130.3</td><td>5</td><td>golgi phosphoprotein 3 (GOLPH3), transcript variant X1, mRNA</td><td>276</td><td>98.551</td><td>1.27e-135</td><td>1,152</td><td>1,427</td></tr><tr><td>AW062801.1</td><td>SAMN00156993</td><td>Colon</td><td>263</td><td>A ↔ T</td><td>NM_020700.1</td><td>12</td><td>protein phosphatase, Mg2+/Mn2+ dependent 1H (PPM1H)</td><td>245</td><td>99.592</td><td>6.00e-124</td><td>4,901</td><td>4,657</td></tr><tr><td>AW062802.1</td><td>SAMN00156993</td><td>Colon</td><td>360</td><td>A ↔ T</td><td>NM_000610.3</td><td>11</td><td>CD44 molecule (Indian blood group) (CD44),</td><td>346</td><td>97.399</td><td>4.38e-165</td><td>2,799</td><td>2,454</td></tr><tr><td>AW062803.1</td><td>SAMN00156993</td><td>Colon</td><td>642</td><td>A ↔ T</td><td>NM_005814.2</td><td>1</td><td>glycoprotein A33 (GPA33), transcript variant X1, mRNA</td><td>351</td><td>92.877</td><td>3.51e-141</td><td>825</td><td>1,171</td></tr><tr><td>AW062804.1</td><td>SAMN00156993</td><td>Colon</td><td>640</td><td>A ↔ T</td><td>NM_000701.7</td><td>1</td><td>ATPase Na+/K+ transporting subunit alpha 1 (ATP1A1</td><td>636</td><td>96.855</td><td>0.0</td><td>3,394</td><td>2,760</td></tr><tr><td>AW062805.1</td><td>SAMN00156993</td><td>Colon</td><td>622</td><td>A ↔ T</td><td>NC_000009.12</td><td>9</td><td>golgi membrane protein 1 (GOLM1)</td><td>400</td><td>94.500</td><td>1.55e-174</td><td>86,074,944</td><td>86,075,342</td></tr><tr><td>AW062807.1*</td><td>SAMN00156993</td><td>Colon</td><td>616</td><td>A ↔ T</td><td>NM_001079524.1</td><td>4</td><td>phosphoribosylaminoimidazole carboxylase and phosphoribosylaminoimidazolesuccinocarboxamide synthase (PAICS)</td><td>582</td><td>88.144</td><td>0.0</td><td>1,737</td><td>2,276</td></tr><tr><td>AW062808.1</td><td>SAMN00156993</td><td>Colon</td><td>325</td><td>A ↔ T</td><td>NM_015658.3</td><td>2</td><td>NOC2 like nucleolar associated transcriptional repressor (NOC2L)</td><td>290</td><td>98.621</td><td>7.55e-143</td><td>953</td><td>1,242</td></tr><tr><td>AW062809.1</td><td>SAMN00156993</td><td>Colon</td><td>609</td><td>A ↔ T</td><td>NM_001278699.2</td><td>3</td><td>RNA polymerase II subunit H (POLR2H)</td><td>544</td><td>97.978</td><td>0.0</td><td>1,229</td><td>1,771</td></tr><tr><td>AW062815.1*</td><td>SAMN00156993</td><td>Colon</td><td>621</td><td>A ↔ T</td><td>NW_017852933.1</td><td>16</td><td>nuclear pore complex-interacting protein family member A7-like (LOC102724993)</td><td>606</td><td>98.515</td><td>0.0</td><td>1,699,909</td><td>1,699,305</td></tr><tr><td>AW062818.1</td><td>SAMN00156993</td><td>Colon</td><td>363</td><td>A ↔ T</td><td>NC_000003.12</td><td>3</td><td>coiled-coil-helix-coiled-coil-helix domain containing 6 (CHCHD6)</td><td>349</td><td>98.281</td><td>7.23e-173</td><td>126,914,442</td><td>126,914,789</td></tr><tr><td>AW062819.1</td><td>SAMN00156993</td><td>Colon</td><td>524</td><td>A ↔ T</td><td>NM_015658.3</td><td>1</td><td>NOC2 like nucleolar associated transcriptional repressor (NOC2L)</td><td>252</td><td>99.603</td><td>7.71e-128</td><td>991</td><td>1,242</td></tr><tr><td>AW062820.1*</td><td>SAMN00156993</td><td>Colon</td><td>508</td><td>A ↔ T</td><td>NM_001293298.1</td><td>15</td><td>cell migration inducing hyaluronidase 1 (CEMIP)</td><td>388</td><td>97.423</td><td>0.0</td><td>6,478</td><td>6,092</td></tr><tr><td>AW062821.1*</td><td>SAMN00156993</td><td>Colon</td><td>336</td><td>A ↔ T</td><td>NM_002295.5</td><td>3</td><td>ribosomal protein SA (RPSA)</td><td>329</td><td>95.441</td><td>9.69e-147</td><td>182</td><td>510</td></tr><tr><td>AW062822.1</td><td>SAMN00156993</td><td>Colon</td><td>350</td><td>A ↔ T</td><td>NM_007278.1</td><td>17</td><td>GABA type A receptor-associated protein (GABARAP)</td><td>155</td><td>96.774</td><td>1.09e-66</td><td>650</td><td>497</td></tr><tr><td>AW062823.1</td><td>SAMN00156993</td><td>Colon</td><td>619</td><td>A ↔ T</td><td>NM_001300832.2</td><td>11</td><td>FERM domain containing 8 (FRMD8)</td><td>499</td><td>98.998</td><td>0.0</td><td>2,751</td><td>2,253</td></tr><tr><td>AW062824.1</td><td>SAMN00156993</td><td>Colon</td><td>294</td><td>A ↔ T</td><td>NM_001731.2</td><td>12</td><td>BTG anti-proliferation factor 1 (BTG1), mRNA</td><td>289</td><td>97.924</td><td>1.63e-139</td><td>1,410</td><td>1,122</td></tr><tr><td>AW062825.1</td><td>SAMN00156993</td><td>Colon</td><td>358</td><td>A ↔ T</td><td>NM_002086.4</td><td>17</td><td>growth factor receptor bound protein 2 (GRB2), transcr</td><td>246</td><td>98.780</td><td>1.30e-120</td><td>2,793</td><td>3,038</td></tr><tr><td>AW062826.1</td><td>SAMN00156993</td><td>Colon</td><td>311</td><td>A ↔ T</td><td>NC_000007.14</td><td>7</td><td>PREDICTED: Homo sapiens speedy/RINGO cell cycle regulator family member E2B (SPDYE2B)</td><td>206</td><td>84.466</td><td>1.86e-49</td><td>102,357,109</td><td>102,357,313</td></tr><tr><td>AW062827.1</td><td>SAMN00156993</td><td>Colon</td><td>500</td><td>A ↔ T</td><td>NM_014455.3</td><td>1</td><td>ring finger protein 115 (RNF115)</td><td>380</td><td>99.474</td><td>0.0</td><td>1,064</td><td>685</td></tr><tr><td>AW062828.1</td><td>SAMN00156993</td><td>Colon</td><td>490</td><td>A ↔ T</td><td>NM_001271969.1</td><td>6</td><td>heat shock protein 90 alpha family class B member 1 (HSP90AB1),</td><td>364</td><td>98.077</td><td>1.54e-179</td><td>2,091</td><td>1,730</td></tr><tr><td>AW062829.1</td><td>SAMN00156993</td><td>Colon</td><td>588</td><td>A ↔ T</td><td>NM_001024216.2</td><td>1</td><td>filamin binding LIM protein 1 (FBLM1)</td><td>525</td><td>99.048</td><td>0.0</td><td>2,354</td><td>2,877</td></tr><tr><td>AW062830.1</td><td>SAMN00156993</td><td>Colon</td><td>191</td><td>A ↔ T</td><td>NM_003302.2</td><td>7</td><td>thyroid hormone receptor interactor 6 (TRIP6)</td><td>173</td><td>96.532</td><td>4.98e-75</td><td>1,639</td><td>1,467</td></tr><tr><td>AW062831.1</td><td>SAMN00156993</td><td>Colon</td><td>157</td><td>A ↔ T</td><td>NM_003749.2</td><td>13</td><td>insulin receptor substrate 2 (IRS2), mRNA</td><td>144</td><td>94.444</td><td>5.12e-55</td><td>5,643</td><td>5,503</td></tr><tr><td>AW175732.1*</td><td>SAMN00156993</td><td>Colon</td><td>392</td><td>A ↔ T</td><td>NC_000015.10</td><td>15</td><td>cell migration inducing hyaluronidase 1 (CEMIP)</td><td>237</td><td>96.203</td><td>1.72e-104</td><td>80,878,828</td><td>80,879,064</td></tr><tr><td>AW175733.1</td><td>SAMN00156993</td><td>Colon</td><td>748</td><td>A ↔ T</td><td>NM_032833.4</td><td>1</td><td>protein phosphatase 1 regulatory subunit 15B (PPP1R15B)</td><td>574</td><td>95.645</td><td>0.0</td><td>71</td><td>643</td></tr><tr><td>AW175736.1</td><td>SAMN00156993</td><td>Colon</td><td>391</td><td>A ↔ T</td><td>NM_005080.3</td><td>22</td><td>X-box binding protein 1 (XBP1)</td><td>284</td><td>97.183</td><td>2.13e-133</td><td>497</td><td>779</td></tr><tr><td>AW175738.1*</td><td>SAMN00156993</td><td>Colon</td><td>195</td><td>A ↔ T</td><td>NC_000003.12</td><td>3</td><td>chromosome 3 open reading frame 14 (C3orf14)</td><td>137</td><td>95.620</td><td>6.63e-54</td><td>62,332,406</td><td>62,332,270</td></tr><tr><td>AW175740.1*</td><td>SAMN00156993</td><td>Colon</td><td>382</td><td>A ↔ T</td><td>NC_000008.11</td><td>8</td><td>pericentriolar material 1 (PCM1)</td><td>288</td><td>97.917</td><td>2.11e-138</td><td>17,968,415</td><td>17,968,700</td></tr><tr><td>AW175741.1</td><td>SAMN00156993</td><td>Colon</td><td>277</td><td>A ↔ T</td><td>NC_000002.12</td><td>2</td><td>phosphodiesterase 11A (PDE11A)</td><td>203</td><td>95.567</td><td>3.80e-86</td><td>177,849,206</td><td>177,849,004</td></tr><tr><td>AW175742.1</td><td>SAMN00156993</td><td>Colon</td><td>701</td><td>A ↔ T</td><td>NC_000001.11</td><td>1</td><td>nuclear receptor subfamily 5 (NR5A2)</td><td>190</td><td>98.421</td><td>1.75e-89</td><td>200,040,067</td><td>200,040,255</td></tr><tr><td>AW175747.1*</td><td>SAMN00156993</td><td>Colon</td><td>182</td><td>A ↔ T</td><td>NC_000001.12</td><td>12</td><td>ribosomal protein SA pseudogene 52 (RPSAP52)</td><td>122</td><td>92.623</td><td>5.23e-40</td><td>65,765,193</td><td>65,765,313</td></tr><tr><td>AW175749.1</td><td>SAMN00156993</td><td>Colon</td><td>379</td><td>A ↔ T</td><td>NC_000003.12</td><td>2</td><td>tigger transposable element also zinc finger protein 621 (ZNF621) derived 1 (TIGD1)</td><td>244</td><td>98.361</td><td>7.82e-118</td><td>3,165,299</td><td>3,165,056</td></tr><tr><td>AW175751.1</td><td>SAMN00156993</td><td>Colon</td><td>403</td><td>A ↔ T</td><td>NT_187585.1</td><td>11</td><td>potassium voltage-gated channel subfamily Q member 1 (KCNQ1)</td><td>336</td><td>94.048</td><td>4.54e-140</td><td>49,605</td><td>49,278</td></tr><tr><td>AW175753.1</td><td>SAMN00156993</td><td>Colon</td><td>800</td><td>A ↔ T</td><td>NM_000005.10</td><td>5</td><td>MCC regulator of WNT signaling pathway (MCC)</td><td>140</td><td>84.286</td><td>6.87e-29</td><td>113,144,132</td><td>113,144,267</td></tr><tr><td>AW176932.1*</td><td>SAMN00156993</td><td>Colon</td><td>655</td><td>A ↔ T</td><td>NM_001079524.1</td><td>4</td><td>phosphoribosylaminoimidazole carboxylase and phosphoribosylaminoimidazolesuccinocarboxamide synthase (PAICS)</td><td>624</td><td>97.917</td><td>0.0</td><td>2,426</td><td>1,803</td></tr><tr><td>AW176933.1*</td><td>SAMN00156993</td><td>Colon</td><td>286</td><td>A ↔ T</td><td>NC_000001.12</td><td>12</td><td>ATPase sarcoplasmic/endoplasmic reticulum Ca2+ transporting 2</td><td>283</td><td>97.173</td><td>3.56e-131</td><td>110,285,901</td><td>110,285,620</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>EST Acc. No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AW176934.1*</td><td>SAMN00156993</td><td>Colon</td><td>294</td><td>A ↔ T</td><td>NC_000012.12</td><td>12</td><td>ATPase sarcoplasmic/endoplasmic reticulum Ca2+ transporting 2</td><td>285</td><td>99.649</td><td>3.49e-146</td><td>110,285,904</td><td>110,285,620</td></tr><tr><td>AW176938.1</td><td>SAMN00156993</td><td>Colon</td><td>183</td><td>A ↔ T</td><td>NM_014280.2</td><td>1</td><td>DnaJ heat shock protein family (Hsp40) member C8 (DNAJC8)</td><td>179</td><td>96.648</td><td>6.40e-79</td><td>661</td><td>483</td></tr><tr><td>AW176940.1</td><td>SAMN00156993</td><td>Colon</td><td>334</td><td>A ↔ T</td><td>NC_000014.9</td><td>14</td><td>serine palmitoyltransferase long chain base subunit 2 (SPTLC2)</td><td>306</td><td>98.693</td><td>2.68e-152</td><td>77,544,607</td><td>77,544,912</td></tr><tr><td>AW176941.1</td><td>SAMN00156993</td><td>Colon</td><td>664</td><td>A ↔ T</td><td>NC_000002.12</td><td>2</td><td>MIR4435-2 host gene (MIR4435-2HG)</td><td>580</td><td>98.276</td><td>0.0</td><td>111,301,721</td><td>111,302,299</td></tr><tr><td>AW176942.1</td><td>SAMN00156993</td><td>Colon</td><td>337</td><td>A ↔ T</td><td>NM_001265603.1</td><td>15</td><td>mortality factor 4 like 1 (MORF4L1)</td><td>279</td><td>98.208</td><td>4.57e-135</td><td>565</td><td>843</td></tr><tr><td>AW176943.1</td><td>SAMN00156993</td><td>Colon</td><td>627</td><td>A ↔ T</td><td>NC_000016.10</td><td>16</td><td>membrane bound transcription factor peptidase, site 1 (MBTPS1)</td><td>500</td><td>98.200</td><td>0.0</td><td>84,061,142</td><td>84,061,638</td></tr><tr><td>AW176945.1</td><td>SAMN00156993</td><td>Colon</td><td>361</td><td>A ↔ T</td><td>NM_030629.2</td><td>16</td><td>c-Maf inducing protein (CMIP)</td><td>326</td><td>98.466</td><td>9.49e-162</td><td>2,775</td><td>2,450</td></tr><tr><td>AW176946.1</td><td>SAMN00156993</td><td>Colon</td><td>664</td><td>A ↔ T</td><td>NM_001873.3</td><td>4</td><td>carboxypeptidase E (CPE)</td><td>307</td><td>98.046</td><td>5.79e-149</td><td>1,193</td><td>890</td></tr><tr><td>AW176947.1</td><td>SAMN00156993</td><td>Colon</td><td>680</td><td>A ↔ T</td><td>NC_000004.12</td><td>4</td><td>huntingtin (HTT)</td><td>660</td><td>98.485</td><td>0.0</td><td>3,200,416</td><td>3,201,074</td></tr><tr><td>AW176948.1</td><td>SAMN00156993</td><td>Colon</td><td>360</td><td>A ↔ T</td><td>NM_001256799.2</td><td>12</td><td>glyceraldehyde-3-phosphate dehydrogenase (GAPDH)</td><td>311</td><td>99.357</td><td>2.05e-158</td><td>831</td><td>1,140</td></tr><tr><td>AW176949.1</td><td>SAMN00156993</td><td>Colon</td><td>304</td><td>A ↔ T</td><td>NM_005731.3</td><td>2</td><td>actin related protein 2/3 complex subunit 2 (ARPC2)</td><td>272</td><td>94.118</td><td>7.87e-113</td><td>867</td><td>597</td></tr><tr><td>AW176952.1</td><td>SAMN00156993</td><td>Colon</td><td>629</td><td>A ↔ T</td><td>NC_000002.12</td><td>2</td><td>integrin subunit alpha 6 (ITGA6)</td><td>630</td><td>99.841</td><td>0.0</td><td>172,430,001</td><td>172,429,372</td></tr><tr><td>AW176953.1</td><td>SAMN00156993</td><td>Colon</td><td>375</td><td>A ↔ T</td><td>NM_001204510.1</td><td>17</td><td>SENP3-EIF4A1 readthrough (NMD candidate) (SENP3-EIF4A1), long non-coding RNA+eukaryotic translation initiation factor 4A1 (EIF4A1)</td><td>338</td><td>99.408</td><td>5.59e-174</td><td>824</td><td>487</td></tr><tr><td>AW176954.1</td><td>SAMN00156993</td><td>Colon</td><td>282</td><td>A ↔ T</td><td>NM_006761.4</td><td>17</td><td>tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein epsilon (YWHAE), non-coding RNA</td><td>264</td><td>98.485</td><td>2.14e-128</td><td>961</td><td>699</td></tr><tr><td>AW176955.1</td><td>SAMN00156993</td><td>Colon</td><td>470</td><td>A ↔ T</td><td>NC_000007.14</td><td>10</td><td>general transcription factor III</td><td>264</td><td>100.000</td><td>3.54e-136</td><td>74,722,574</td><td>74,722,311</td></tr><tr><td>AW176958.1</td><td>SAMN00156993</td><td>Colon</td><td>255</td><td>A ↔ T</td><td>NM_001291484.2</td><td>19</td><td>carcinoembryonic antigen related cell adhesion molecule 5 (CEACAM5)</td><td>227</td><td>99.119</td><td>1.02e-111</td><td>1,130</td><td>1,356</td></tr><tr><td>AW176959.1*</td><td>SAMN00156993</td><td>Colon</td><td>682</td><td>A ↔ T</td><td>NM_001402.5</td><td>6</td><td>eukaryotic translation elongation factor 1 alpha 1 (EEF1A1)</td><td>647</td><td>92.427</td><td>0.0</td><td>923</td><td>1,566</td></tr><tr><td>AW176960.1*</td><td>SAMN00156993</td><td>Colon</td><td>622</td><td>A ↔ T</td><td>NM_001402.5</td><td>6</td><td>eukaryotic translation elongation factor 1 alpha 1 (EEF1A1)</td><td>610</td><td>99.672</td><td>0.0</td><td>891</td><td>1,499</td></tr><tr><td>AW176962.1</td><td>SAMN00156993</td><td>Colon</td><td>309</td><td>A ↔ T</td><td>NM_005514.7</td><td>6</td><td>major histocompatibility complex, class I, B (HLA-B)</td><td>276</td><td>97.464</td><td>1.28e-130</td><td>1,338</td><td>1,063</td></tr><tr><td>AW176963.1*</td><td>SAMN00156993</td><td>Colon</td><td>451</td><td>A ↔ T</td><td>NC_000001.11</td><td>1</td><td>erythrocyte membrane protein band 4 (EPB41)</td><td>381</td><td>98.163</td><td>0.0</td><td>29,114,294</td><td>29,113,915</td></tr><tr><td>AW176964.1*</td><td>SAMN00156993</td><td>Colon</td><td>392</td><td>A ↔ T</td><td>NC_000001.11</td><td>1</td><td>erythrocyte membrane protein band 4 (EPB41)</td><td>381</td><td>98.688</td><td>0.0</td><td>29,114,294</td><td>29,113,915</td></tr><tr><td>AW176965.1</td><td>SAMN00156993</td><td>1Colon</td><td>628</td><td>A ↔ T</td><td>NM_001137604.2</td><td>2</td><td>RNA polymerase I subunit B (POLR1B)</td><td>421</td><td>99.050</td><td>0.0</td><td>3,159</td><td>3,579</td></tr><tr><td>AW176968.1</td><td>SAMN00156993</td><td>Colon</td><td>633</td><td>A ↔ T</td><td>NM_001172663.1</td><td>16</td><td>RAB40C, member RAS oncogene family (RAB40C)</td><td>160</td><td>98.750</td><td>1.79e-74</td><td>1,801</td><td>1,643</td></tr><tr><td>AW176969.1</td><td>SAMN00156993</td><td>Colon</td><td>348</td><td>A ↔ T</td><td>NM_002588.3</td><td>5</td><td>protocadherin gamma subfamily B, 3 (PCDHGB3)</td><td>350</td><td>97.714</td><td>1.57e-169</td><td>4,400</td><td>4,051</td></tr><tr><td>AW176971.1*</td><td>SAMN00156993</td><td>Colon</td><td>650</td><td>A ↔ T</td><td>NM_001402.5</td><td>6</td><td>eukaryotic translation elongation factor 1 alpha 1 (EEF1A1)</td><td>307</td><td>98.697</td><td>7.44e-153</td><td>1,305</td><td>1,000</td></tr><tr><td>AW176972.1</td><td>SAMN00156993</td><td>Colon</td><td>374</td><td>A ↔ T</td><td>NM_001167671.2</td><td>3</td><td>LIM domain containing preferred translocation partner in lipoma (LPP)</td><td>369</td><td>98.103</td><td>0.0</td><td>16,859</td><td>17,225</td></tr><tr><td>AW176973.1</td><td>SAMN00156993</td><td>Colon</td><td>645</td><td>A ↔ T</td><td>NC_000016.10</td><td>16</td><td>zinc and ring finger 1 (ZNRF1)</td><td>587</td><td>98.978</td><td>0.0</td><td>75,058,197</td><td>75,058,783</td></tr><tr><td>AW176974.1</td><td>SAMN00156993</td><td>Colon</td><td>652</td><td>A ↔ T</td><td>NM_001276418.1</td><td>9</td><td>SEC16 homolog A, endoplasmic reticulum export factor (SEC16A)</td><td>585</td><td>99.316</td><td>0.0</td><td>8,695</td><td>8,111</td></tr><tr><td>AW176982.1*</td><td>SAMN00156993</td><td>Colon</td><td>560</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>cox1</td><td>513</td><td>99.220</td><td>0.0</td><td>6,537</td><td>6,025</td></tr><tr><td>AW903505.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>436</td><td>A → T → C → G → A</td><td>NR_148543.1</td><td>2</td><td>clathrin heavy chain linker domain containing 1 (CLHC1)</td><td>437</td><td>98.169</td><td>0.0</td><td>569</td><td>1,001</td></tr><tr><td>AW903508.1</td><td>SAMN00159106</td><td>Nervous normal</td><td>443</td><td>A → T → C → G → A</td><td>NM_015018.3</td><td>6</td><td>DOP1 leucine zipper like protein A (DOP1A)</td><td>362</td><td>99.448</td><td>0.0</td><td>7,625</td><td>7,264</td></tr><tr><td>AW903509.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>449</td><td>A → T → C → G → A</td><td>NM_001351724.1</td><td>4</td><td>Rap guanine nucleotide exchange factor 2 (RAPGEF2)</td><td>451</td><td>97.561</td><td>0.0</td><td>5,954</td><td>6,400</td></tr><tr><td>AW903510.1</td><td>SAMN00159106</td><td>Nervous normal</td><td>472</td><td>A → T → C → G → A</td><td>NM_004859.3</td><td>17</td><td>clathrin heavy chain (CLTC)</td><td>406</td><td>99.015</td><td>0.0</td><td>3,799</td><td>4,202</td></tr><tr><td>AW903511.1</td><td>SAMN00159106</td><td>Nervous normal</td><td>464</td><td>A → T → C → G → A</td><td>NC_000014.9</td><td>14</td><td>endoplasmic reticulum oxidoreductase 1 alpha (ERO1A)</td><td>444</td><td>97.973</td><td>0.0</td><td>52,686,677</td><td>52,687,114</td></tr><tr><td>AW903512.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>185</td><td>A → T → C → G → A</td><td>NM_024996.5</td><td>3</td><td>G elongation factor mitochondrial 1 (GFM1)</td><td>186</td><td>90.323</td><td>9.61e-61</td><td>3,446</td><td>3,264</td></tr><tr><td>AW903513.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>254</td><td>A → T → C → G → A</td><td>NM_024996.5</td><td>3</td><td>G elongation factor mitochondrial 1 (GFM1)</td><td>166</td><td>100.000</td><td>2.60e-81</td><td>3,237</td><td>3,072</td></tr><tr><td>AW903514.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>423</td><td>A → T → C → G → A</td><td>NM_001351724.1</td><td>4</td><td>Rap guanine nucleotide exchange factor 2 (RAPGEF2)</td><td>424</td><td>99.057</td><td>0.0</td><td>5,962</td><td>6,384</td></tr><tr><td>AW903515.1</td><td>SAMN00159106</td><td>Nervous normal</td><td>447</td><td>A → T → C → G → A</td><td>NC_000008.11</td><td>8</td><td>mitochondrial calcium uptake family member 3</td><td>447</td><td>100.000</td><td>0.0</td><td>17,028,849</td><td>17,028,403</td></tr><tr><td>AW903516.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>446</td><td>A → T → C → G → A</td><td>NM_007353.2</td><td>7</td><td>G protein subunit alpha 12 (GNA12)</td><td>409</td><td>99.511</td><td>0.0</td><td>3,841</td><td>4,249</td></tr><tr><td>AW903517.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>382</td><td>A → T → C → G → A</td><td>NM_004529.3</td><td>9</td><td>MLLT3, super elongation complex subunit (MLLT3)</td><td>379</td><td>99.736</td><td>0.0</td><td>4,755</td><td>4,377</td></tr><tr><td>AW903518.1*</td><td>SAMN00159106</td><td>Nervous normal</td><td>261</td><td>A → T → C → G → A</td><td>NC_012920.1</td><td>mito</td><td>cox1</td><td>133</td><td>99.248</td><td>9.61e-61</td><td>6,634</td><td>6,766</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>EST Acc. No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AW903519.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>457</td><td>A → T → C → G → A</td><td>NC_000001.11</td><td>1</td><td>membrane associated guanylate kinase (MAGI3)</td><td>443</td><td>97.065</td><td>0.0</td><td>113,403,203</td><td>113,403,645</td></tr><tr><td>AW903520.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>461</td><td>A → T → C → G → A</td><td>NR_029435.1</td><td>14</td><td>PSMA3 antisense RNA 1 (PSMA3-AS1), long non-coding RNA</td><td>447</td><td>100.000</td><td>0.0</td><td>2,426</td><td>1,980</td></tr><tr><td>AW903521.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>449</td><td>A → T → C → G → A</td><td>NC_012920.1</td><td>mito</td><td>cox1</td><td>451</td><td>98.670</td><td>0.0</td><td>6,534</td><td>6,982</td></tr><tr><td>AW903523.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>447</td><td>A → T → C → G → A</td><td>NC_000004.12</td><td>4</td><td>NEDD4 binding protein 2 (N4BP2)</td><td>445</td><td>99.775</td><td>0.0</td><td>40,063,375</td><td>40,063,375</td></tr><tr><td>AW903524.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>455</td><td>A → T → C → G → A</td><td>NM_001286691.1</td><td>9</td><td>MLLT3, super elongation complex subunit (MLLT3)</td><td>449</td><td>100.000</td><td>0.0</td><td>4,671</td><td>4,223</td></tr><tr><td>AW903525.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>431</td><td>A → T → C → G → A</td><td>NC_000003.12</td><td>3</td><td>C-X9-C motif containing 1 (269)</td><td>269</td><td>99.628</td><td>2.40e-136</td><td>28,242,970</td><td>28,242,702</td></tr><tr><td>AW903526.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>466</td><td>A → T → C → G → A</td><td>NR_029435.1</td><td>14</td><td>PSMA3 antisense RNA 1 (PSMA3-AS1) and long non-coding RNA</td><td>456</td><td>99.781</td><td>0.0</td><td>2,426</td><td>1,971</td></tr><tr><td>AW903527.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>446</td><td>A → T → C → G → A</td><td>NM_004529.3</td><td>9</td><td>MLLT3, super elongation complex subunit (MLLT3)</td><td>446</td><td>99.776</td><td>0.0</td><td>4,213</td><td>4,658</td></tr><tr><td>AW903528.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>461</td><td>A → T → C → G → A</td><td>NR_148543.1</td><td>2</td><td>clathrin heavy chain linker domain containing 1 (CLHC1)</td><td>457</td><td>99.344</td><td>0.0</td><td>1,034</td><td>578</td></tr><tr><td>AW903529.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>464</td><td>A → T → C → G → A</td><td>NR_049782.1</td><td>2</td><td>sorting nexin 17 (SNX17)</td><td>455</td><td>98.901</td><td>0.0</td><td>1,025</td><td>573</td></tr><tr><td>AW903530.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>418</td><td>A → T → C → G → A</td><td>NM_004800.2</td><td>13</td><td>transmembrane 9 superfamily member 2 (TM9SF2)</td><td>395</td><td>98.734</td><td>0.0</td><td>2,509</td><td>2,115</td></tr><tr><td>AW903531.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>408</td><td>A → T → C → G → A</td><td>NC_000022.11</td><td>22</td><td>phosphatidylinositol 4-kinase alpha (PI4KA)</td><td>407</td><td>97.543</td><td>0.0</td><td>20,799,225</td><td>20,799,629</td></tr><tr><td>AW903532.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>279</td><td>A → T → C → G → A</td><td>NC_000003.12</td><td>3</td><td>zinc finger protein 148 (ZNF148)</td><td>194</td><td>98.969</td><td>5.50e-93</td><td>125,335,569</td><td>125,335,762</td></tr><tr><td>AW903533.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>473</td><td>A → T → C → G → A</td><td>NM_001293092.1</td><td>7</td><td>G protein subunit alpha 12 (GNA12)</td><td>409</td><td>99.267</td><td>0.0</td><td>3,790</td><td>4,198</td></tr><tr><td>AW903534.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>469</td><td>A → T → C → G → A</td><td>NC_000007.14</td><td>7</td><td>ataxin 7 like 1 (ATXN7L1)</td><td>461</td><td>99.132</td><td>0.0</td><td>105,869,501</td><td>105,869,961</td></tr><tr><td>AW903537.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>448</td><td>A → T → C → G → A</td><td>NC_000016.10</td><td>16</td><td>RNA binding fox-1 homolog 1 (RBFOX1)</td><td>431</td><td>99.304</td><td>0.0</td><td>6,537,483</td><td>6,537,913</td></tr><tr><td>AW903538.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>454</td><td>A → T → C → G → A</td><td>NC_000023.11</td><td>X</td><td>euroligin 4 X-linked (NLGN4X)</td><td>454</td><td>99.559</td><td>0.0</td><td>6,115,554</td><td>6,115,102</td></tr><tr><td>AW903539.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>321</td><td>A → T → C → G → A</td><td>NC_000011.10</td><td>11</td><td>chromosome 11 open reading frame 80 (C11orf80)</td><td>248</td><td>93.952</td><td>9.08e-101</td><td>66,766,914</td><td>66,767,161</td></tr><tr><td>AW903540.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>329</td><td>A → T → C → G → A</td><td>NC_000007.14</td><td>7</td><td>family with sequence similarity 221 member A (FAM221A)</td><td>329</td><td>100.000</td><td>6.36e-172</td><td>23,695,685</td><td>23,696,013</td></tr><tr><td>AW903541.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>470</td><td>A → T → C → G → A</td><td>NC_000005.10</td><td>5</td><td>teneurin transmembrane protein 2 (TENM2)</td><td>455</td><td>99.341</td><td>0.0</td><td>167,499,052</td><td>167,499,506</td></tr><tr><td>AW903542.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>464</td><td>A → T → C → G → A</td><td>NM_024996.5</td><td>3</td><td>G elongation factor mitochondrial 1 (GFM1)</td><td>455</td><td>99.560</td><td>0.0</td><td>2,980</td><td>3,434</td></tr><tr><td>AW903543.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>470</td><td>A → T → C → G → A</td><td>NC_000001.11</td><td>1</td><td>vacuolar protein sorting 72 homolog (VPS72)</td><td>468</td><td>99.359</td><td>0.0</td><td>151,182,693</td><td>151,183,160</td></tr><tr><td>AW903544.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>458</td><td>A → T → C → G → A</td><td>NM_018947.5</td><td>7</td><td>cytochrome c, somatic (CYCS), mRNA</td><td>460</td><td>99.130</td><td>0.0</td><td>703</td><td>245</td></tr><tr><td>AW903545.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>458</td><td>A → T → C → G → A</td><td>NC_000005.10</td><td>5</td><td>teneurin transmembrane protein 2 (TENM2)</td><td>459</td><td>99.782</td><td>0.0</td><td>167,499,058</td><td>167,499,516</td></tr><tr><td>AW903546.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>236</td><td>A → T → C → G → A</td><td>NC_000011.10</td><td>11</td><td>GRB2 associated binding protein 2 (GAB2)</td><td>154</td><td>97.403</td><td>2.05e-67</td><td>78,340,987</td><td>78,341,140</td></tr><tr><td>AW903547.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>441</td><td>A → T → C → G → A</td><td>NC_000001.11</td><td>1</td><td>membrane associated guanylate kinase (MAGI3)</td><td>427</td><td>97.190</td><td>0.0</td><td>113,403,203</td><td>113,403,629</td></tr><tr><td>AW903550.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>460</td><td>A → T → C → G → A</td><td>NM_172350.2</td><td>1</td><td>CD46 molecule (CD46)</td><td>462</td><td>99.134</td><td>0.0</td><td>2,971</td><td>2,510</td></tr><tr><td>AW903551.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>215</td><td>A → T → C → G → A</td><td>NC_012920.1</td><td>mito</td><td>cox1</td><td>109</td><td>99.083</td><td>2.11e-47</td><td>7,416</td><td>7,308</td></tr><tr><td>AW903553.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>468</td><td>A → T → C → G → A</td><td>NM_001351562.1</td><td>9</td><td>TLE family member 4 (TLE4)</td><td>469</td><td>97.655</td><td>0.0</td><td>4176</td><td>4643</td></tr><tr><td>AW903554.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>357</td><td>A → T → C → G → A</td><td>NM_018947.5</td><td>7</td><td>cytochrome c, somatic (CYCS), mRNA</td><td>355</td><td>99.155</td><td>0.0</td><td>682</td><td>328</td></tr><tr><td>AW903556.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>444</td><td>A → T → C → G → A</td><td>NM_004529.3</td><td>9</td><td>MLLT3, super elongation complex subunit (MLLT3)</td><td>446</td><td>99.327</td><td>0.0</td><td>4,746</td><td>4,301</td></tr><tr><td>AW903557.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>458</td><td>A → T → C → G → A</td><td>NC_000012.12</td><td>12</td><td>synaptotagmin 1 (SYT1)</td><td>462</td><td>98.701</td><td>0.0</td><td>79,289,394</td><td>79,288,933</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>EST Acc. No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AW903558.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>476</td><td>A → T → C → G → A</td><td>NC_000012.12</td><td>12</td><td>kinesin family member 21A (KIF21A)</td><td>477</td><td>99.161</td><td>0.0</td><td>39,406,699</td><td>39,407,174</td></tr><tr><td>AW903559.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>453</td><td>A → T → C → G → A</td><td>NC_000017.11</td><td>17</td><td>dynein axonemal heavy chain 17 (DNAH17)</td><td>290</td><td>99.310</td><td>6.59e-147</td><td>78,428,317</td><td>78,428,606</td></tr><tr><td>AW903560.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>459</td><td>A → T → C → G → A</td><td>NC_000006.12</td><td>6</td><td>capping protein regulator and myosin 1 linker 1 (CARMIL1)</td><td>403</td><td>96.030</td><td>0.0</td><td>25,338,780</td><td>25,338,379</td></tr><tr><td>AW903561.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>463</td><td>A → T → C → G → A</td><td>NC_000014.9</td><td>14</td><td>PSMA3 antisense RNA 1 (PSMA3-AS1) and long non-coding RNA</td><td>449</td><td>98.886</td><td>0.0</td><td>58,265,829</td><td>58,265,381</td></tr><tr><td>AW903563.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>427</td><td>A → T → C → G → A</td><td>NM_001352634.1</td><td>8</td><td>pericentriolar material 1 (PCM1)</td><td>423</td><td>99.291</td><td>0.0</td><td>5,954</td><td>6,376</td></tr><tr><td>AW903564.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>457</td><td>A → T → C → G → A</td><td>NM_025054.4</td><td>8</td><td>valosin containing protein interacting protein 1 (VCPIP1)</td><td>374</td><td>99.198</td><td>0.0</td><td>7,603</td><td>7,976</td></tr><tr><td>AW903565.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>434</td><td>A → T → C → G → A</td><td>NC_000007.14</td><td>7</td><td>ataxin 7 like 1 (ATXN7L1)</td><td>427</td><td>99.297</td><td>0.0</td><td>105,869,531</td><td>105,869,956</td></tr><tr><td>AW903566.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>283</td><td>A → T → C → G → A</td><td>NM_015541.2</td><td>3</td><td>leucine rich repeats and immunoglobulin like domains 1 (LRIG1)</td><td>190</td><td>98.421</td><td>1.54e-88</td><td>3,447</td><td>3,258</td></tr><tr><td>AW903567.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>444</td><td>A → T → C → G → A</td><td>NC_000022.11</td><td>22</td><td>RNA binding fox-1 homolog 2 (RBFOX2)</td><td>344</td><td>92.442</td><td>3.11e-135</td><td>35,867,779</td><td>35,868,118</td></tr><tr><td>AW903568.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>428</td><td>A → T → C → G → A</td><td>NM_020755.3</td><td>6</td><td>serine incorporator 1 (SERINC1)</td><td>326</td><td>99.080</td><td>1.07e-164</td><td>2,332</td><td>2,008</td></tr><tr><td>AW903570.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>454</td><td>A → T → C → G → A</td><td>NT_187614.1</td><td>17</td><td>aminopeptidase puromycin sensitive pseudogene (NPEPPS1)</td><td>200</td><td>97.000</td><td>1.54e-88</td><td>2,236,540</td><td>2,236,343</td></tr><tr><td>AW903571.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>480</td><td>A → T → C → G → A</td><td>NM_001242866.2</td><td>21</td><td>protein arginine methyltransferase 2 (PRMT2)</td><td>297</td><td>99.663</td><td>1.82e-152</td><td>1,248</td><td>952</td></tr><tr><td>AW903573.1*</td><td>SAMN00159106</td><td>Nervous_normal</td><td>336</td><td>A → T → C → G → A</td><td>NC_012920.1</td><td>mito</td><td>cox1</td><td>236</td><td>98.729</td><td>1.15e-114</td><td>7,112</td><td>7,347</td></tr><tr><td>AW903576.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>472</td><td>A → T → C → G → A</td><td>NM_203401.1</td><td>4</td><td>stathmin 1 (STMN1)</td><td>468</td><td>98.932</td><td>0.0</td><td>361</td><td>827</td></tr><tr><td>AW903577.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>466</td><td>A → T → C → G → A</td><td>NC_000007.14</td><td>7</td><td>LUC7 like 2 (LUC7L2), pre-mRNA splicing factor</td><td>461</td><td>98.915</td><td>0.0</td><td>139,362,483</td><td>139,362,023</td></tr><tr><td>BE085177.1*</td><td>SAMN00159564</td><td>Breast</td><td>269</td><td>A → T → C → G → A</td><td>NC_000010.11</td><td>10</td><td>protein tyrosine phosphatase receptor type E (PTPRE)</td><td>242</td><td>98.347</td><td>6.87e-117</td><td>127,960,155</td><td>127,959,914</td></tr><tr><td>BE085178.1</td><td>SAMN00159564</td><td>Breast</td><td>176</td><td>A → T → C → G → A</td><td>NM_015888.4</td><td>1</td><td>hook microtubule tethering protein 1 (HOOK1)</td><td>171</td><td>92.982</td><td>2.67e-61</td><td>2,185</td><td>2,024</td></tr><tr><td>BE085183.1</td><td>SAMN00159564</td><td>Breast</td><td>293</td><td>A → T → C → G → A</td><td>NC_000010.11</td><td>10</td><td>SMC5-SMC6 complex localization factor 2 (SLF2)</td><td>92</td><td>95.652</td><td>5.99e-33</td><td>100,947,017</td><td>100,947,107</td></tr><tr><td>BE085186.1*</td><td>SAMN00159564</td><td>Breast</td><td>302</td><td>A → T → C → G → A</td><td>NM_002840.4</td><td>1</td><td>protein tyrosine phosphatase, receptor type F (PTPRF)</td><td>268</td><td>96.642</td><td>2.45e-121</td><td>2,985</td><td>2,724</td></tr><tr><td>BE085188.1</td><td>SAMN00159564</td><td>Breast</td><td>280</td><td>A → T → C → G → A</td><td>NM_007042.4</td><td>3</td><td>ribonuclease P/MRP subunit p14 (RPP14)</td><td>193</td><td>96.373</td><td>1.55e-83</td><td>328</td><td>139</td></tr><tr><td>BE085197.1</td><td>SAMN00159564</td><td>Breast</td><td>196</td><td>A → T → C → G → A</td><td>NM_001281724.2</td><td>3</td><td>biotinidase (BTD)</td><td>106</td><td>99.057</td><td>2.73e-46</td><td>2,275</td><td>2,170</td></tr><tr><td>BE085202.1</td><td>SAMN00159564</td><td>Breast</td><td>170</td><td>A → T → C → G → A</td><td>NM_021009.6</td><td>12</td><td>ubiquitin C (UBC)</td><td>102</td><td>89.216</td><td>1.01e-25</td><td>1,869</td><td>1,772</td></tr><tr><td>BE145238.1</td><td>SAMN00157869</td><td>Head_neck</td><td>333</td><td>A ↔ T</td><td>NM_001042465.2</td><td>10</td><td>prosaposin (PSAP)</td><td>334</td><td>97.904</td><td>1.50e-163</td><td>1,739</td><td>2,072</td></tr><tr><td>BE145239.1*</td><td>SAMN00157869</td><td>Head_neck</td><td>293</td><td>A ↔ T</td><td>NM_001104629.1</td><td>4</td><td>chromosome 4 open reading frame 19 (C4orf19)</td><td>292</td><td>97.260</td><td>5.59e-138</td><td>1,828</td><td>2,119</td></tr><tr><td>BE145240.1</td><td>SAMN00157869</td><td>Head_neck</td><td>169</td><td>A ↔ T</td><td>NM_001033925.1</td><td>10</td><td>TIA1 cytotoxic granule associated RNA binding protein like 1 (TIAL1)</td><td>169</td><td>96.450</td><td>4.73e-74</td><td>1,013</td><td>1,181</td></tr><tr><td>BE145242.1</td><td>SAMN00157869</td><td>Head_neck</td><td>452</td><td>A ↔ T</td><td>NM_022470.3</td><td>3</td><td>zinc finger matrin-type 3 (ZMAT3)</td><td>404</td><td>98.267</td><td>0.0</td><td>6,273</td><td>6,675</td></tr><tr><td>BE145243.1</td><td>SAMN00157869</td><td>Head_neck</td><td>110</td><td>A ↔ T</td><td>NM_001101662.1</td><td>1</td><td>nardilysin convertase (NRDC)</td><td>110</td><td>99.091</td><td>1.76e-48</td><td>1,445</td><td>1,336</td></tr><tr><td>BE145245.1</td><td>SAMN00157869</td><td>Head_neck</td><td>421</td><td>A ↔ T</td><td>NC_000005.10</td><td>5</td><td>small integral membrane protein 3 (SMIM3)</td><td>416</td><td>97.837</td><td>0.0</td><td>150,784,346</td><td>150,784,758</td></tr><tr><td>BE145246.1</td><td>SAMN00157869</td><td>Head_neck</td><td>266</td><td>A ↔ T</td><td>NM_001033047.2</td><td>4</td><td>nephronectin (NPNT)</td><td>267</td><td>92.884</td><td>9.75e-106</td><td>391</td><td>656</td></tr><tr><td>BE145247.1</td><td>SAMN00157869</td><td>Head_neck</td><td>359</td><td>A ↔ T</td><td>NW_019805496.1</td><td>11</td><td>mitochondrial carrier 2 (MTCH2)</td><td>356</td><td>96.067</td><td>5.39e-163</td><td>55,682</td><td>56,035</td></tr><tr><td>BE145248.1*</td><td>SAMN00157869</td><td>Head_neck</td><td>247</td><td>A ↔ T</td><td>NC_000014.9</td><td>14</td><td>valosin containing protein lysine methyltransferase (VCPKMT)</td><td>240</td><td>92.083</td><td>9.96e-91</td><td>50,108,378</td><td>50,108,617</td></tr><tr><td>BE145250.1*</td><td>SAMN00157869</td><td>Head_neck</td><td>289</td><td>A ↔ T</td><td>NC_000014.9</td><td>14</td><td>valosin containing protein lysine methyltransferase (VCPKMT)</td><td>290</td><td>95.517</td><td>1.58e-128</td><td>50,108,285</td><td>50,108,574</td></tr><tr><td>BE145251.1</td><td>SAMN00157869</td><td>Head_neck</td><td>329</td><td>A ↔ T</td><td>NC_000009.12</td><td>9</td><td>spectrin alpha, non-erythrocytic 1 (SPTAN1)</td><td>330</td><td>96.061</td><td>9.15e-151</td><td>128,610,482</td><td>128,610,155</td></tr><tr><td>BF332836.1</td><td>SAMN00160205</td><td>Breast</td><td>298</td><td>A → T → C → G → A</td><td>NC_000001.11</td><td>1</td><td>putative homeodomain transcription factor 1 (PHTF1)</td><td>230</td><td>95.652</td><td>4.22e-99</td><td>113,710,355</td><td>113,710,126</td></tr><tr><td>BF332837.1</td><td>SAMN00160205</td><td>Breast</td><td>376</td><td>A → T → C → G → A</td><td>NC_000002.12</td><td>2</td><td>BCL2 like 11(BCL2 L11)</td><td>241</td><td>98.340</td><td>3.20e-115</td><td>111,159,962</td><td>111,160,199</td></tr><tr><td>BF332838.1</td><td>SAMN00160205</td><td>Breast</td><td>359</td><td>A → T → C → G → A</td><td>NC_000008.11</td><td>8</td><td>integrator complex subunit 9 (INTS9)</td><td>237</td><td>97.468</td><td>1.16e-109</td><td>28,784,748</td><td>28,784,981</td></tr><tr><td>BF332840.1</td><td>SAMN00160205</td><td>Breast</td><td>300</td><td>A → T → C → G → A</td><td>NM_018364.4</td><td>1</td><td>round spermatid basic protein 1 (RSBN1)</td><td>209</td><td>94.737</td><td>1.20e-84</td><td>3,267</td><td>3,062</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>EST Acc. No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>BF332841.1</td><td>SAMN00160205</td><td>Breast</td><td>372</td><td>A → T → C → G → A</td><td>NM_002266.3</td><td>17</td><td>karyopherin subunit alpha 2 (KPNA2)</td><td>269</td><td>96.654</td><td>6.82e-122</td><td>2,022</td><td>1,757</td></tr><tr><td>BF332842.1</td><td>SAMN00160205</td><td>Breast</td><td>378</td><td>A → T → C → G → A</td><td>NM_003761.4</td><td>2</td><td>vesicle associated membrane protein 8 (VAMP8)</td><td>253</td><td>99.209</td><td>8.76e-126</td><td>296</td><td>548</td></tr><tr><td>BF332846.1</td><td>SAMN00160205</td><td>Breast</td><td>280</td><td>A → T → C → G → A</td><td>NM_017409.3</td><td>12</td><td>homeobox C10 (HOXC10)</td><td>190</td><td>95.789</td><td>4.35e-79</td><td>1,958</td><td>1,773</td></tr><tr><td>BF332847.1</td><td>SAMN00160205</td><td>Breast</td><td>405</td><td>A → T → C → G → A</td><td>NC_000007.14</td><td>7</td><td>3-hydroxyisobutyrate dehydrogenase (HIBADH)</td><td>275</td><td>97.818</td><td>2.42e-131</td><td>27,629,869</td><td>27,629,596</td></tr><tr><td>BF332849.1</td><td>SAMN00160205</td><td>Breast</td><td>308</td><td>A → T → C → G → A</td><td>NC_000014.9</td><td>14</td><td>dishevelled associated activator of morphogenesis 1 (DAAM1)</td><td>206</td><td>95.631</td><td>1.99e-87</td><td>59,268,775</td><td>59,268,570</td></tr><tr><td>BF332850.1</td><td>SAMN00160205</td><td>Breast</td><td>383</td><td>A → T → C → G → A</td><td>NC_000012.12</td><td>12</td><td>spermatogenesis associated serine rich 2 (SPATS2)</td><td>296</td><td>96.284</td><td>4.02e-134</td><td>49,458,046</td><td>49,457,752</td></tr><tr><td>BF332851.1</td><td>SAMN00160205</td><td>Breast</td><td>335</td><td>A → T → C → G → A</td><td>NC_000013.11</td><td>13</td><td>SMAD family member 9 (SMAD9)</td><td>249</td><td>98.394</td><td>8.83e-121</td><td>36,916,214</td><td>36,916,462</td></tr><tr><td>BF332854.1</td><td>SAMN00160205</td><td>Breast</td><td>293</td><td>A → T → C → G → A</td><td>NC_000003.12</td><td>3</td><td>coiled-coil domain containing 66 (CCDC66)</td><td>221</td><td>98.190</td><td>4.19e-104</td><td>56,561,670</td><td>56,561,451</td></tr><tr><td>BF332855.1</td><td>SAMN00160205</td><td>Breast</td><td>309</td><td>A → T → C → G → A</td><td>NC_000014.9</td><td>14</td><td>RAD51 paralog B (RAD51B)</td><td>222</td><td>96.396</td><td>1.97e-97</td><td>68,178,422</td><td>68,178,203</td></tr><tr><td>BF353453.1</td><td>SAMN00161174</td><td>Head_neck</td><td>395</td><td>A ↔ T</td><td>NM_002274.3</td><td>17</td><td>keratin 13 (KRT13)</td><td>53</td><td>100.000</td><td>8.30e-19</td><td>1,679</td><td>1,627</td></tr><tr><td>BF354529.1*</td><td>SAMN00160318</td><td>Head_neck</td><td>397</td><td>A → T → C → G → A</td><td>NM_015596.2</td><td>19</td><td>kallikrein related peptidase 13 (KLK13)</td><td>317</td><td>98.107</td><td>1.09e-154</td><td>1,348</td><td>1,664</td></tr><tr><td>BF354541.1</td><td>SAMN00160318</td><td>Head_neck</td><td>363</td><td>A → T → C → G → A</td><td>NM_030634.2</td><td>1</td><td>zinc finger protein 436 (ZNF436)</td><td>295</td><td>98.644</td><td>8.52e-146</td><td>3,915</td><td>4,209</td></tr><tr><td>BF354544.1</td><td>SAMN00160318</td><td>Head_neck</td><td>449</td><td>A → T → C → G → A</td><td>NC_000006.12</td><td>6</td><td>protein phosphatase 1 regulatory subunit 10 (PPP1R10)</td><td>384</td><td>98.698</td><td>0.0</td><td>30,600,715</td><td>30,601,097</td></tr><tr><td>BF354553.1</td><td>SAMN00160318</td><td>Head_neck</td><td>387</td><td>A → T → C → G → A</td><td>NW_016107297.1</td><td>5</td><td>GUSB Pseudogene (GUSBP1)</td><td>332</td><td>98.193</td><td>1.79e-162</td><td>163,285</td><td>162,954</td></tr><tr><td>BF354559.1</td><td>SAMN00160318</td><td>Head_neck</td><td>193</td><td>A → T → C → G → A</td><td>XM_011524110.3</td><td>17</td><td>small G protein signaling modulator 2 (SGSM2)</td><td>99</td><td>93.939</td><td>5.99e-33</td><td>2,457</td><td>2,360</td></tr><tr><td>BF354564.1*</td><td>SAMN00160318</td><td>Head_neck</td><td>402</td><td>A → T → C → G → A</td><td>NM_015596.2</td><td>19</td><td>kallikrein related peptidase 13 (KLK13)</td><td>306</td><td>98.366</td><td>1.09e-149</td><td>1,359</td><td>1,664</td></tr><tr><td>BF354572.1</td><td>SAMN00160318</td><td>Head_neck</td><td>389</td><td>A → T → C → G → A</td><td>XM_011514699.2</td><td>6</td><td>glucose-fructose oxidoreductase domain containing 1 (GFOD1)</td><td>346</td><td>98.844</td><td>3.80e-174</td><td>3,008</td><td>3,353</td></tr><tr><td>BF370229.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>282</td><td>A → T → C → G → A</td><td>NM_000138.4</td><td>15</td><td>fibrillin 1 (FBN1)</td><td>127</td><td>95.276</td><td>5.87e-48</td><td>8,263</td><td>8,388</td></tr><tr><td>BF370230.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>282</td><td>A → T → C → G → A</td><td>NM_006887.4</td><td>15</td><td>ZFP36 ring finger protein like 2 (ZFP36L2)</td><td>346</td><td>97.399</td><td>3.85e-164</td><td>2,544</td><td>2,201</td></tr><tr><td>BF370232.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>410</td><td>A → T → C → G → A</td><td>NC_012920.1</td><td>mito</td><td>atp6</td><td>359</td><td>98.329</td><td>6.31e-177</td><td>8,663</td><td>9,020</td></tr><tr><td>BF370233.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>437</td><td>A → T → C → G → A</td><td>NC_000023.11</td><td>x</td><td>ubiquitin specific peptidase 9 X-linked (USP9X)</td><td>373</td><td>98.928</td><td>0.0</td><td>41,186,360</td><td>41,185,989</td></tr><tr><td>BF370234.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>413</td><td>A → T → C → G → A</td><td>NC_000012.12</td><td>12</td><td>lysine demethylase 5A (KDM5A)</td><td>349</td><td>97.994</td><td>2.96e-170</td><td>326,420</td><td>326,075</td></tr><tr><td>BF370235.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>294</td><td>A → T → C → G → A</td><td>NM_001288733.1</td><td>17</td><td>testis expressed 2 (TEX2)</td><td>127</td><td>98.425</td><td>9.68e-56</td><td>4,533</td><td>4,659</td></tr><tr><td>BF370236.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>294</td><td>A → T → C → G → A</td><td>NC_000009.12</td><td>9</td><td>basonuclin 2 (BNC2)</td><td>194</td><td>98.969</td><td>5.50e-93</td><td>16,514,557</td><td>16,514,750</td></tr><tr><td>BF370237.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>423</td><td>A → T → C → G → A</td><td>NM_001146227.1</td><td>8</td><td>ribosomal protein S20 (RPS20)</td><td>304</td><td>99.013</td><td>6.54e-152</td><td>493</td><td>191</td></tr><tr><td>BF370238.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>332</td><td>A → T → C → G → A</td><td>NC_000004.12</td><td>4</td><td>sec1 family domain containing 2 (SCFD2)</td><td>243</td><td>98.354</td><td>2.47e-116</td><td>52,874,276</td><td>52,874,518</td></tr><tr><td>BF370239.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>355</td><td>A → T → C → G → A</td><td>NM_018046.4</td><td>5</td><td>angiogenic factor with G-patch and FHA domains 1 (AGGF1)</td><td>293</td><td>98.294</td><td>1.43e-143</td><td>4,058</td><td>4,350</td></tr><tr><td>BF370241.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>305</td><td>A → T → C → G → A</td><td>NC_000017.11</td><td>17</td><td>tubulin folding cofactor D (TBCD)</td><td>173</td><td>98.844</td><td>2.60e-81</td><td>82,761,458</td><td>82,761,630</td></tr><tr><td>BF370242.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>250</td><td>A → T → C → G → A</td><td>NC_000001.11</td><td>1</td><td>GNG12, DIRAS3 and WLS antisense RNA 1 (GNG12-AS1) and Wnt ligand secretion mediator (WLS)</td><td>127</td><td>98.425</td><td>9.68e-56</td><td>68,182,918</td><td>68,182,792</td></tr><tr><td>BF370243.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>413</td><td>A → T → C → G → A</td><td>NC_000005.10</td><td>5</td><td>ubiquitin conjugating enzyme E2 D2 (UBE2D2)</td><td>272</td><td>97.059</td><td>8.76e-126</td><td>139,610,960</td><td>139,610,690</td></tr><tr><td>BF798643.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>420</td><td>A ↔ T</td><td>NR_045104.1</td><td>15</td><td>integrator complex subunit 14 (INTS14)</td><td>290</td><td>92.759</td><td>5.75e-114</td><td>163</td><td>448</td></tr><tr><td>BF798647.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>415</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>12s</td><td>271</td><td>94.096</td><td>3.46e-111</td><td>915</td><td>649</td></tr><tr><td>BF798651.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>363</td><td>A ↔ T</td><td>NM_032593.2</td><td>9</td><td>histidine triad nucleotide binding protein 2 (HINT2)</td><td>57</td><td>94.737</td><td>5.12e-15</td><td>429</td><td>484</td></tr><tr><td>BF798653.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>425</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>12s</td><td>270</td><td>94.815</td><td>5.75e-114</td><td>649</td><td>915</td></tr><tr><td>BF798654.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>375</td><td>A ↔ T</td><td>NM_001291920.1</td><td>9</td><td>retinoid X receptor alpha (RXRA)</td><td>136</td><td>95.588</td><td>8.10e-53</td><td>2,040</td><td>2,170</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>EST Length</th><th>Swinger transformation</th><th>Aligned sequence AI</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>e-value</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>BF798660.1</td><td>SAMN00162090</td><td>Colon_ibs</td><td>464</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>12s</td><td>269</td><td>97.398</td><td>1.22e-125</td><td>649</td><td>915</td></tr><tr><td>BF798665.1</td><td>SAMN00162090</td><td>Colon_ibs</td><td>485</td><td>A ↔ T</td><td>NM_004069.4</td><td>19</td><td>adaptor related protein complex 2 sigma 1 subunit (AP2S1)</td><td>373</td><td>98.123</td><td>0.0</td><td>564</td><td>193</td></tr><tr><td>BF798667.1</td><td>SAMN00162090</td><td>Colon_ibs</td><td>521</td><td>A ↔ T</td><td>NM_005561.3</td><td>13</td><td>lysosomal associated membrane protein 1 (LAMP1)</td><td>481</td><td>95.842</td><td>0.0</td><td>392</td><td>872</td></tr></tbody></table>

Genes corresponding to more than one EST sequence are indicated by *.

Two A ↔ T-transformed swinger RNAs detected in the colon and head-neck cancer lines mapped on the adenylate kinase 1 (AK1) coding gene (Table 2). AK1 plays an important role in tumor suppression and it is often downregulated in cancer cells (Collavin et al., 1999; Janssen et al., 2004; Vasseur et al., 2005; Jan et al., 2019). Three swinger RNAs detected in breast (A → G → C → T → A) and two colon cancer (A ↔ T) lines mapped on cystatin B (CSTB) (Table 2). CSTB plays an important role in expression and epigenetic regulation and is downregulated in lung, gastric and colorectal cancers (Zhang et al., 2016; Ma et al., 2017). CSTB promotes cell proliferation, migration and suppresses apoptosis in gastric cancer cells (Zhang et al., 2016). Downregulation of CSTB also promotes gastric cancer (Zhang et al., 2016). Similarly, swinger RNAs for HSPA8 and PABPC1 were identified. We identified two swinger RNAs from the HSPA8 protein-coding region. The carboxy-terminus of Hsc70 interacting protein (CHIP) plays an important role in cancer initiation and progression (Hatakeyama et al., 2005; Kajiro et al., 2009; Gaude et al., 2012) and has an anti-tumor effect in many cancer types including colon and gastric cancers (Kajiro et al., 2009; Ahmed et al., 2012; Wang et al., 2013; Ying et al., 2013; Wang T et al., 2014; Wang Y et al., 2014). Thirty-three swinger RNAs were transcribed from the PABPC1 gene from colon cancer and head-neck cancer biosamples, respectively (Table 2). PABPC proteins are RNA processing proteins associated with gene expression regulation (Liu et al., 2012) and are upregulated in prostate and colorectal cancers (Eisermann et al., 2015). PABPC also has a tumor suppressor role in head and neck squamous cell carcinoma (Zeng et al., 2018). Swingertransformed RNAs produced from these identified genes will produce non-homologous and putatively nonfunctional mRNAs translating dysfunctional nonhomologous proteins, which supports the association between downregulation of these genes and cancer cell types. These swinger RNAs could be early stage factors responsible for cancer induction or result from genetic instability in later stages of malignant cancers (Ruckova et al., 2012; Ciocca and Calderwood, 2005). Results favor the former because the swinger transformed mRNA of tumor suppressor genes would result in disruption of gene function inhibiting programmed cell death. Analyses also identify two swinger RNAs mapped on DLG associated proteins (DLGAP1 and DLGAP4) coding genes in healthy nervous tissue samples and colon cancer cells. DLGAP is a protein overexpressed in the brain (Fagerberg et al., 2014) and promotes invasiveness in cancer cell lines (Li et al., 2018). Similarly, three swinger RNAs were transcribed from mediator complex subunit genes (MED) in colon and breast cancer tissues (Table 2). Mutations in MED12 are associated with tumorigenesis (Bullerdiek and Rommel, 2018; Xie et al., 2018) and cause benign breast fibroepithelial lesions (Pareja et al., 2019). These two observations suggest that swinger transformation of MED should induce carcinogenesis, whereas, swinger RNAs are transcribed from overexpressing DLGAP genes due to enzyme fatigue in cancer cells.

Among the detected swinger RNAs, 51 swinger RNAs match mitochondrial genes (Tables 1 and 2). It is widely known and
<table><caption>TABLE 2 | List of ESTs based on the genes identified in more than one biosample.</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>Length of EST</th><th>Swinger transformation</th><th>Acc. No. of aligned sequence</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>evalue</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>BE085181.1</td><td>SAMN00159564</td><td>Breast</td><td>109</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>16s</td><td>41</td><td>95.122</td><td>6.21e<br>-08</td><td>2,885</td><td>2,924</td></tr><tr><td>BE085203.1</td><td>SAMN00159564</td><td>Breast</td><td>282</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>16s</td><td>189</td><td>98.413</td><td>1.54e<br>-88</td><td>1,950</td><td>1,764</td></tr><tr><td>BE085205.1</td><td>SAMN00159564</td><td>Breast</td><td>134</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>16s</td><td>89</td><td>97.753</td><td>1.29e<br>-34</td><td>2,882</td><td>2,794</td></tr><tr><td>BE085207.1</td><td>SAMN00159564</td><td>Breast</td><td>245</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>16s</td><td>185</td><td>97.838</td><td>1.20e<br>-84</td><td>2,838</td><td>3,021</td></tr><tr><td>BF332839.1</td><td>SAMN00160205</td><td>Breast</td><td>249</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>16s</td><td>156</td><td>97.436</td><td>5.70e<br>-68</td><td>2,258</td><td>2,105</td></tr><tr><td>BF798658.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>624</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>16s+trnAL+ND1</td><td>464</td><td>98.276</td><td>0.0</td><td>3,075</td><td>3,537</td></tr><tr><td>BF798678.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>589</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>16s+trnAL+ND1</td><td>464</td><td>97.414</td><td>0.0</td><td>3,075</td><td>3,537</td></tr><tr><td>AI939948.1</td><td>SAMN00156461</td><td>Colon</td><td>529</td><td>A ↔ T</td><td>NM_000476.2</td><td>9</td><td>adenylate kinase 1 (AK1)</td><td>504</td><td>97.024</td><td>0.0</td><td>198</td><td>700</td></tr><tr><td>AI940588.1</td><td>SAMN00156483</td><td>Head_neck</td><td>650</td><td>A ↔ T</td><td>NM_001318121.1</td><td>9</td><td>adenylate kinase 1 (AK1)</td><td>574</td><td>98.606</td><td>0.0</td><td>177</td><td>749</td></tr><tr><td>BF332844.1</td><td>SAMN00160205</td><td>Breast</td><td>251</td><td>A → T → C →<br>G → A</td><td>NM_000100.3</td><td>21</td><td>cystatin B (CSTB)</td><td>170</td><td>96.471</td><td>7.32e<br>-72</td><td>452</td><td>287</td></tr><tr><td>AW176935.1</td><td>SAMN00156993</td><td>Colon</td><td>386</td><td>A ↔ T</td><td>NC_000021.9</td><td>21</td><td>cystatin b (CSTB)</td><td>349</td><td>99.140</td><td>2.58e<br>-177</td><td>43,774,237</td><td>43,774,583</td></tr><tr><td>AW176936.1</td><td>SAMN00156993</td><td>Colon</td><td>368</td><td>A ↔ T</td><td>NC_000021.9</td><td>21</td><td>cystatin b (CSTB)</td><td>348</td><td>99.425</td><td>5.55e<br>-179</td><td>43,774,237</td><td>43,774,583</td></tr><tr><td>AW903549.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>444</td><td>A → T → C →<br>G → A</td><td>NM_001308390.1</td><td>18</td><td>DLG associated protein 1 (DLGAP1)</td><td>446</td><td>99.327</td><td>0.0</td><td>4,908</td><td>4,464</td></tr><tr><td>AW176970.1</td><td>SAMN00156993</td><td>Colon</td><td>701</td><td>A ↔ T</td><td>NC_000020.11</td><td>20</td><td>DLG associated protein 4 (DLGAP4)</td><td>685</td><td>98.686</td><td>0.0</td><td>36,525,990</td><td>36,526,674</td></tr><tr><td>BF332852.1</td><td>SAMN00160205</td><td>Breast</td><td>406</td><td>A → T → C →<br>G → A</td><td>NM_153201.3</td><td>11</td><td>heat shock protein family A (Hsp70)<br>member 8 (HSPA8)</td><td>302</td><td>97.351</td><td>1.84e<br>-142</td><td>925</td><td>626</td></tr><tr><td>BF798649.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>529</td><td>A ↔ T</td><td>NM_006597.5</td><td>11</td><td>heat shock protein family A (Hsp70)<br>member 8 (HSPA8)</td><td>486</td><td>95.062</td><td>0.0</td><td>386</td><td>865</td></tr><tr><td>AW062795.1</td><td>SAMN00156993</td><td>Colon</td><td>639</td><td>A ↔ T</td><td>NM_005120.2</td><td>X</td><td>mediator complex subunit 12<br>(MED12)</td><td>502</td><td>99.801</td><td>0.0</td><td>4,475</td><td>3,974</td></tr><tr><td>BF332843.1</td><td>SAMN00160205</td><td>Breast</td><td>340</td><td>A → T → C →<br>G → A</td><td>NM_005121.2</td><td>17</td><td>mediator complex subunit 13<br>(MED13)</td><td>214</td><td>98.131</td><td>9.08e<br>-101</td><td>943</td><td>731</td></tr><tr><td>BF332848.1</td><td>SAMN00160205</td><td>Breast</td><td>372</td><td>A → T → C →<br>G → A</td><td>NM_001271811.1</td><td>12</td><td>mediator complex subunit 21<br>(MED21)</td><td>245</td><td>97.959</td><td>2.47e<br>-116</td><td>1,097</td><td>854</td></tr><tr><td>AW176957.1</td><td>SAMN00156993</td><td>Colon</td><td>517</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd1</td><td>495</td><td>99.192</td><td>0.0</td><td>3,663</td><td>4,157</td></tr><tr><td>BF798657.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>363</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd1</td><td>210</td><td>93.333</td><td>1.30e<br>-80</td><td>3,333</td><td>3,537</td></tr><tr><td>AW903507.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>215</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>214</td><td>99.533</td><td>2.51e<br>-106</td><td>4,739</td><td>4,952</td></tr><tr><td>BF332853.1</td><td>SAMN00160205</td><td>Breast</td><td>308</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>216</td><td>99.074</td><td>9.02e<br>-106</td><td>4,724</td><td>4,939</td></tr><tr><td>BF354531.1</td><td>SAMN00160318</td><td>Head_neck</td><td>373</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>242</td><td>98.760</td><td>1.48e<br>-118</td><td>4,482</td><td>4,723</td></tr><tr><td>BF354537.1</td><td>SAMN00160318</td><td>Head_neck</td><td>345</td><td>A → T → C →<br>G → A</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>232</td><td>100.000</td><td>5.31e<br>-118</td><td>4,492</td><td>4,723</td></tr><tr><td>BI032899.1</td><td>SAMN00162655</td><td>Nervous_normal</td><td>424</td><td>A ↔ C</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>351</td><td>84.046</td><td>2.80e<br>-87</td><td>5,371</td><td>5,033</td></tr><tr><td>AI939946.1</td><td>SAMN00156461</td><td>Colon</td><td>488</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>393</td><td>97.710</td><td>0.0</td><td>5,062</td><td>4,670</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>Length of EST</th><th>Swinger transformation</th><th>Acc. No. of aligned sequence</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>evalue</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AI939952.1</td><td>SAMN00156802</td><td>Colon</td><td>263</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>236</td><td>95.339</td><td>8.21e<br>-101</td><td>4,958</td><td>4,723</td></tr><tr><td>AI939954.1</td><td>SAMN00156802</td><td>Colon</td><td>252</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>225</td><td>96.444</td><td>2.95e<br>-100</td><td>4,945</td><td>4,723</td></tr><tr><td>AI939955.1</td><td>SAMN00156802</td><td>Colon</td><td>212</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>113</td><td>99.115</td><td>1.14e<br>-49</td><td>4,905</td><td>4,794</td></tr><tr><td>AI939956.1</td><td>SAMN00156802</td><td>Colon</td><td>200</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>186</td><td>98.387</td><td>6.48e<br>-87</td><td>4,930</td><td>4,746</td></tr><tr><td>AI939957.1</td><td>SAMN00156802</td><td>Colon</td><td>269</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>218</td><td>97.706</td><td>2.28e<br>-101</td><td>4,751</td><td>4,967</td></tr><tr><td>AI939958.1</td><td>SAMN00156802</td><td>Colon</td><td>223</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>198</td><td>96.465</td><td>6.48e<br>-87</td><td>4,919</td><td>4,723</td></tr><tr><td>AI939960.1</td><td>SAMN00156802</td><td>Colon</td><td>286</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>245</td><td>94.694</td><td>2.28e<br>-101</td><td>4,962</td><td>4,723</td></tr><tr><td>AI939964.1</td><td>SAMN00156802</td><td>Colon</td><td>251</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>229</td><td>96.070</td><td>2.95e<br>-100</td><td>4,949</td><td>4,723</td></tr><tr><td>AI939967.1</td><td>SAMN00156802</td><td>Colon</td><td>281</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>245</td><td>97.551</td><td>3.74e<br>-114</td><td>4,723</td><td>4,967</td></tr><tr><td>AI939971.1</td><td>SAMN00156802</td><td>Colon</td><td>251</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>240</td><td>92.500</td><td>1.36e<br>-90</td><td>4,962</td><td>4,723</td></tr><tr><td>AI939974.1</td><td>SAMN00156802</td><td>Colon</td><td>243</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>231</td><td>94.805</td><td>3.74e<br>-96</td><td>4,952</td><td>4,723</td></tr><tr><td>AI939975.1</td><td>SAMN00156802</td><td>Colon</td><td>152</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>152</td><td>96.711</td><td>5.05e<br>-65</td><td>4,898</td><td>4,748</td></tr><tr><td>AI939977.1</td><td>SAMN00156802</td><td>Colon</td><td>226</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>201</td><td>95.522</td><td>1.37e<br>-85</td><td>4,922</td><td>4,723</td></tr><tr><td>AI939978.1</td><td>SAMN00156802</td><td>Colon</td><td>212</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>188</td><td>97.340</td><td>1.77e<br>-84</td><td>4,914</td><td>4,729</td></tr><tr><td>AI939979.1</td><td>SAMN00156802</td><td>Colon</td><td>162</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>163</td><td>95.706</td><td>3.02e<br>-67</td><td>4,949</td><td>4,789</td></tr><tr><td>AI939982.1</td><td>SAMN00156802</td><td>Colon</td><td>248</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>233</td><td>95.708</td><td>3.72e<br>-101</td><td>4,955</td><td>4,723</td></tr><tr><td>AI939984.1</td><td>SAMN00156802</td><td>Colon</td><td>233</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>214</td><td>94.393</td><td>2.93e<br>-87</td><td>4,933</td><td>4,723</td></tr><tr><td>AI939985.1</td><td>SAMN00156802</td><td>Colon</td><td>145</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>122</td><td>96.721</td><td>1.86e<br>-49</td><td>4,847</td><td>4,967</td></tr><tr><td>AI939986.1</td><td>SAMN00156802</td><td>Colon</td><td>166</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>163</td><td>93.865</td><td>8.45e<br>-63</td><td>4,891</td><td>4,729</td></tr><tr><td>AI940581.1</td><td>SAMN00156483</td><td>Head_neck</td><td>533</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>392</td><td>99.490</td><td>0.0</td><td>4,671</td><td>5,062</td></tr><tr><td>AI940582.1</td><td>SAMN00156483</td><td>Head_neck</td><td>509</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>393</td><td>99.237</td><td>0.0</td><td>4,670</td><td>5,062</td></tr><tr><td>AI940585.1</td><td>SAMN00156483</td><td>Head_neck</td><td>475</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>394</td><td>96.701</td><td>0.0</td><td>4,670</td><td>5,062</td></tr><tr><td>AI940586.1</td><td>SAMN00156483</td><td>Head_neck</td><td>296</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>294</td><td>98.299</td><td>5.83e<br>-144</td><td>4,692</td><td>4,984</td></tr><tr><td>AI940587.1</td><td>SAMN00156483</td><td>Head_neck</td><td>474</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd2</td><td>393</td><td>99.237</td><td>0.0</td><td>5,062</td><td>4,670</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>Length of EST</th><th>Swinger transformation</th><th>Acc. No. of aligned sequence</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>evalue</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>BF354521.1</td><td>SAMN00160318</td><td>Head neck</td><td>383</td><td>A → T → C → G → A</td><td>NC_012920.1</td><td>mito</td><td>nd4</td><td>343</td><td>95.335</td><td>2.35e-151</td><td>11,016</td><td>11,354</td></tr><tr><td>CK327105.1</td><td>SAMN00157633</td><td>Colon</td><td>474</td><td>A ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd4</td><td>391</td><td>99.233</td><td>0.0</td><td>11,611</td><td>11,999</td></tr><tr><td>BF370011.1</td><td>SAMN00160853</td><td>Prostate normal</td><td>252</td><td>C ↔ T</td><td>NC_012920.1</td><td>mito</td><td>nd4</td><td>196</td><td>96.429</td><td>2.65e-85</td><td>11,610</td><td>11,801</td></tr><tr><td>AW062763.1</td><td>SAMN00156993</td><td>Colon</td><td>624</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>387</td><td>97.674</td><td>0.0</td><td>1,882</td><td>1,497</td></tr><tr><td>AW062774.1</td><td>SAMN00156993</td><td>Colon</td><td>673</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.753</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW062781.1</td><td>SAMN00156993</td><td>Colon</td><td>197</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>186</td><td>99.462</td><td>1.36e-90</td><td>1,707</td><td>1,523</td></tr><tr><td>AW062786.1</td><td>SAMN00156993</td><td>Colon</td><td>583</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.012</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW062790.1</td><td>SAMN00156993</td><td>Colon</td><td>685</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>95.309</td><td>0.0</td><td>1,888</td><td>1,485</td></tr><tr><td>AW062791.1</td><td>SAMN00156993</td><td>Colon</td><td>533</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.012</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW062797.1</td><td>SAMN00156993</td><td>Colon</td><td>331</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>321</td><td>99.065</td><td>2.64e-162</td><td>1,485</td><td>1,805</td></tr><tr><td>AW062798.1</td><td>SAMN00156993</td><td>Colon</td><td>312</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>303</td><td>97.690</td><td>1.25e-145</td><td>1,503</td><td>1,805</td></tr><tr><td>AW062806.1</td><td>SAMN00156993</td><td>Colon</td><td>472</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>409</td><td>99.511</td><td>0.0</td><td>1,893</td><td>1,485</td></tr><tr><td>AW062810.1</td><td>SAMN00156993</td><td>Colon</td><td>625</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>625</td><td>99.520</td><td>0.0</td><td>2,430</td><td>1,806</td></tr><tr><td>AW062811.1</td><td>SAMN00156993</td><td>Colon</td><td>426</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.753</td><td>0.0</td><td>1,485</td><td>1,889</td></tr><tr><td>AW062812.1</td><td>SAMN00156993</td><td>Colon</td><td>420</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>404</td><td>99.257</td><td>0.0</td><td>1,486</td><td>1,889</td></tr><tr><td>AW062813.1</td><td>SAMN00156993</td><td>Colon</td><td>688</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>100.000</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW062814.1</td><td>SAMN00156993</td><td>Colon</td><td>443</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.012</td><td>0.0</td><td>1,485</td><td>1,889</td></tr><tr><td>AW062816.1</td><td>SAMN00156993</td><td>Colon</td><td>536</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.753</td><td>0.0</td><td>1,485</td><td>1,889</td></tr><tr><td>AW062817.1</td><td>SAMN00156993</td><td>Colon</td><td>250</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>221</td><td>93.665</td><td>4.91e-85</td><td>1,574</td><td>1,781</td></tr><tr><td>AW176937.1</td><td>SAMN00156993</td><td>Colon</td><td>360</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>346</td><td>97.977</td><td>7.29e-168</td><td>1,485</td><td>1,825</td></tr><tr><td>AW176944.1</td><td>SAMN00156993</td><td>Colon</td><td>603</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.012</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW176950.1</td><td>SAMN00156993</td><td>Colon</td><td>419</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>86</td><td>100.000</td><td>3.15e-37</td><td>1,804</td><td>1,889</td></tr><tr><td>AW176951.1</td><td>SAMN00156993</td><td>Colon</td><td>192</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>142</td><td>100.000</td><td>2.33e-68</td><td>1,748</td><td>1,889</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>Length of EST</th><th>Swinger transformation</th><th>Acc. No. of aligned sequence</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>evalue</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AW176956.1</td><td>SAMN00156993</td><td>Colon</td><td>348</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>319</td><td>99.060</td><td>3.41e-161</td><td>1,485</td><td>1,803</td></tr><tr><td>AW176966.1</td><td>SAMN00156993</td><td>Colon</td><td>497</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.506</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW176967.1</td><td>SAMN00156993</td><td>Colon</td><td>158</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>159</td><td>97.484</td><td>1.39e-70</td><td>1,917</td><td>2,074</td></tr><tr><td>AW176975.1</td><td>SAMN00156993</td><td>Colon</td><td>481</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>373</td><td>97.855</td><td>0.0</td><td>1,485</td><td>1,854</td></tr><tr><td>AW176976.1</td><td>SAMN00156993</td><td>Colon</td><td>605</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>354</td><td>88.701</td><td>2.17e-118</td><td>1,834</td><td>1,485</td></tr><tr><td>AW176977.1</td><td>SAMN00156993</td><td>Colon</td><td>639</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.259</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>AW176978.1</td><td>SAMN00156993</td><td>Colon</td><td>434</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>407</td><td>99.017</td><td>0.0</td><td>1,485</td><td>1,889</td></tr><tr><td>AW176979.1</td><td>SAMN00156993</td><td>Colon</td><td>599</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>406</td><td>99.507</td><td>0.0</td><td>1,485</td><td>1,889</td></tr><tr><td>AW176980.1</td><td>SAMN00156993</td><td>Colon</td><td>426</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.012</td><td>0.0</td><td>1,485</td><td>1,889</td></tr><tr><td>AW176981.1</td><td>SAMN00156993</td><td>Colon</td><td>599</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>405</td><td>99.012</td><td>0.0</td><td>1,889</td><td>1,485</td></tr><tr><td>BF333938.1</td><td>SAMN00157869</td><td>Head_neck</td><td>665</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>79</td><td>100.000</td><td>2.92e-33</td><td>1,563</td><td>1,485</td></tr><tr><td>BF333939.1</td><td>SAMN00157869</td><td>Head_neck</td><td>611</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>319</td><td>98.746</td><td>1.89e-159</td><td>1,803</td><td>1,485</td></tr><tr><td>BF333940.1</td><td>SAMN00157869</td><td>Head_neck</td><td>509</td><td>A ↔ T</td><td>NM_002568.3</td><td>8</td><td>poly(A) binding protein cytoplasmic 1 (PABPC1)</td><td>311</td><td>91.318</td><td>4.33e-116</td><td>1,792</td><td>1,485</td></tr></tbody></table>
<table><caption>TABLE 3 | List of ESTs mapped on non-coding genomic regions and uncharacterized genes.</caption><thead><tr><th>EST Acc No.</th><th>Biosample</th><th>Cancer tissue</th><th>Length of EST</th><th>Swinger transformation</th><th>Acc No of aligned sequence</th><th>Chromosome</th><th>Gene</th><th>Aligned Seq. length</th><th>%ID</th><th>evalue</th><th>5'</th><th>3'</th></tr></thead><tbody><tr><td>AI939983.1</td><td>SAMN00156802</td><td>Colon</td><td>379</td><td>A ↔ T</td><td>NC_000009.12</td><td>9</td><td>Non-coding region</td><td>360</td><td>97.778</td><td>4.32e-175</td><td>120,855,482</td><td>1,21E+08</td></tr><tr><td>AI939987.1</td><td>SAMN00156802</td><td>Colon</td><td>367</td><td>A ↔ T</td><td>NC_000009.12</td><td>9</td><td>Non-coding region</td><td>345</td><td>97.971</td><td>7.29e-168</td><td>120,855,464</td><td>1,21E+08</td></tr><tr><td>AW062766.1</td><td>SAMN00156993</td><td>Colon</td><td>457</td><td>A ↔ T</td><td>NC_000016.10</td><td>16</td><td>Non-coding region</td><td>313</td><td>99.681</td><td>9.49e-162</td><td>74,318,755</td><td>74319067</td></tr><tr><td>AW175734.1</td><td>SAMN00156993</td><td>Colon</td><td>328</td><td>A ↔ T</td><td>NC_000004.12</td><td>4</td><td>Non-coding region</td><td>136</td><td>90.441</td><td>1.12e-41</td><td>131,268,224</td><td>1,31E+08</td></tr><tr><td>AW175735.1</td><td>SAMN00156993</td><td>Colon</td><td>620</td><td>A ↔ T</td><td>NC_000004.12</td><td>4</td><td>Non-coding region</td><td>191</td><td>95.288</td><td>1.78e-79</td><td>131,268,167</td><td>1,31E+08</td></tr><tr><td>AW175739.1</td><td>SAMN00156993</td><td>Colon</td><td>222</td><td>A ↔ T</td><td>NC_000008.11</td><td>8</td><td>Non-coding region</td><td>139</td><td>97.122</td><td>1.83e-59</td><td>102,600,346</td><td>1,03E+08</td></tr><tr><td>AW175744.1</td><td>SAMN00156993</td><td>Colon</td><td>506</td><td>A ↔ T</td><td>NC_000008.11</td><td>8</td><td>Non-coding region</td><td>313</td><td>97.764</td><td>3.46e-151</td><td>88,009,959</td><td>88,009,648</td></tr><tr><td>AW903506.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>422</td><td>A → T → C → G → A</td><td>NC_000006.12</td><td>6</td><td>Non-coding region</td><td>350</td><td>97.143</td><td>8.28e-166</td><td>132,009,130</td><td>1,32E+08</td></tr><tr><td>AW903535.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>440</td><td>A → T → C → G → A</td><td>NC_000013.11</td><td>13</td><td>Non-coding region</td><td>442</td><td>99.321</td><td>0.0</td><td>95,650,104</td><td>95,650,545</td></tr><tr><td>AW903555.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>434</td><td>A → T → C → G → A</td><td>NC_000010.11</td><td>10</td><td>Non-coding region</td><td>436</td><td>99.312</td><td>0.0</td><td>20,554,187</td><td>20,553,752</td></tr><tr><td>AW903569.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>430</td><td>A → T → C → G → A</td><td>NC_000004.12</td><td>4</td><td>Non-coding region</td><td>357</td><td>98.599</td><td>1.36e-178</td><td>17,483,414</td><td>17,483,770</td></tr><tr><td>AW903572.1</td><td>SAMN00159106</td><td>Nervous_normal</td><td>248</td><td>A → T → C → G → A</td><td>NC_000013.11</td><td>13</td><td>Non-coding region</td><td>56</td><td>100.000</td><td>3.66e-20</td><td>95,650,461</td><td>95,650,406</td></tr><tr><td>BE085219.1</td><td>SAMN00159564</td><td>Breast</td><td>314</td><td>A → T → C → G → A</td><td>NC_000003.12</td><td>3</td><td>Non-coding region</td><td>311</td><td>97.749</td><td>3.94e-149</td><td>169,791,191</td><td>1,7E+08</td></tr><tr><td>BF332859.1</td><td>SAMN00160205</td><td>Breast</td><td>211</td><td>A → T → C → G → A</td><td>NC_000020.11</td><td>20</td><td>Non-coding region</td><td>143</td><td>90.210</td><td>2.12e-42</td><td>59,862,273</td><td>59,862,140</td></tr><tr><td>BF332861.1</td><td>SAMN00160205</td><td>Breast</td><td>364</td><td>A → T → C → G → A</td><td>NC_000010.11</td><td>10</td><td>Non-coding region</td><td>256</td><td>97.266</td><td>4.11e-119</td><td>8,810,014</td><td>8,810,267</td></tr><tr><td>BF354533.1</td><td>SAMN00160318</td><td>Head neck</td><td>388</td><td>A → T → C → G → A</td><td>NC_000005.10</td><td>5</td><td>Non-coding region</td><td>331</td><td>99.698</td><td>2.29e-171</td><td>69,932,968</td><td>69,933,298</td></tr><tr><td>BF370228.1</td><td>SAMN00160862</td><td>Prostate_normal</td><td>264</td><td>A → T → C → G → A</td><td>NC_000023.11</td><td>x</td><td>Non-coding region</td><td>195</td><td>95.897</td><td>2.01e-82</td><td>136,928,879</td><td>1,37E+08</td></tr><tr><td>AW062783.1</td><td>SAMN00156993</td><td>Colon</td><td>640</td><td>A ↔ T</td><td>NC_000009.12</td><td>9</td><td>Uncharacterized LOC101927086</td><td>557</td><td>75.404</td><td>5.09e-60</td><td>70,473,455</td><td>70,472,931</td></tr><tr><td>BF798648.1</td><td>SAMN00162090</td><td>Colon_ins</td><td>344</td><td>A ↔ T</td><td>NC_000003.12</td><td>3</td><td>Uncharacterized LOC102723512</td><td>148</td><td>96.622</td><td>2.87e-62</td><td>4,006,872</td><td>4,006,726</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Relationship between A → T → C → G → A and A → G → C → T → A asymmetric swinger transformation.</figcaption></figure>

proven that mitochondrial genes are overexpressed in various cancer types to meet up the metabolic requirements of cancer cells and are strongly associated with cancer (Tarantul et al., 2001; Modica-Napolitano and Singh, 2004; Księżakowska-Łakoma et al., 2014; Lin et al., 2018). Among 51 swinger RNAs, 7, 1, 30 and 3 swinger RNAs are transcribed from mitochondrial 16s rRNA, ND1, ND2 and ND4, respectively. These were identified in several biosamples (Table 2). Sixty-six percent (33 ESTs) of the identified swinger RNAs were transcribed from the MT-ND region. A previous study showed bias for swinger transformation of MT-ND genes (Warthi and Seligmann, 2018).

Interestingly, detection of swinger transformed DLGAP1 mRNAs (overexpressed ~15× in neurons compared to other tissue cells; Fagerberg et al., 2014), in normal nervous tissue samples suggests that the swinger transformation is not directly associated with cancer but perhaps associated with highly expressed genes. Carcinogenesis could be the outcome of swinger transformed dysfunctional mRNAs. Previously identified swinger RNAs from highly active cancer-associated genes (Warthi et al., 2019) also support this finding. Indeed, even within mitochondrial genes, rRNAs are more expressed than other genes, and swinger rRNAs are the most frequently observed mitochondrial swinger RNAs in previous publications (Seligmann, 2013b; Seligmann, 2013c; Seligmann, 2014b; Seligmann, 2015b), matching the pattern of positive association between regular and swinger transcriptions. Hence, the positive association between expression and swinger transformation occurs independently for human nuclear and mitochondrial genes. Hence, at least at this qualitative level of analysis, observations on these highly expressed nuclear and mitochondrial genes support that swinger transformations associate with highly active genomic regions that result in cancer induction and progression due to nonfunctional transcripts.

In order to test whether swinger transformations are biased towards some regions of genes, genes for which swinger RNAs were detected in more than one biosample (Table 2) are regionally compartmentalized in three equal regions i.e. 5' region, mid region and 3' regions, each region spanning 1/3 of the gene. Mitochondrial genes like 16S rRNA (five of seven ESTs), ND4 (two of three ESTs), and nuclear gene Cystatin B (three of three ESTs) mapped on the 3' extremity of their respective genes.

Among the 30 ESTs mapped on ND2 gene across multiple biosamples, 28 mapped on the central region (5' to 3'; 4,692 to 4,984 bp) of the mitochondrial genome (NC_012920.1), while two mapped on the 5' terminal of ND2 (4,482–4723 bp) gene. Interestingly, for MT-ND4, among the three ESTs identified across three different biosamples for three different swinger transformations, two ESTs mapped on the same 3' terminal ND4 region i.e. 11,610–11,801 bp on the mitochondrial genome (NC_012920.1). For Nuclear gene poly(A) binding protein cytoplasmic 1 (PABPC1), all the 33 identified ESTs preferentially mapped on the mid-region (5' 1,485 bp to 3' 1,803 bp on NM_002568.3) of PABPC1 gene.

Considering mitochondrial and nuclear genes separately, 38 among 43 (mitochondrial, 88.4%, one tailed sign test P = 0.0011) and 33 among 33 (nuclear, 100%, one tailed sign test P = 1.6 × 10⁻⁶) swinger RNAs map either on the mid or the 3' regions of the gene. These tests assume that the probability of mapping randomly on these regions is 2/3. One tailed sign tests are justified by the working hypothesis that polymerase enzyme fatigue (occurring more downstream from transcription initiation) causes swinger transcription. These observations indicate that the mid and 3' regions of highly expressed genes are more prone to produce swinger transformed RNAs than the 5' region, for each nuclear and mitochondrial genes. Similarly, swinger transformed ESTs (from only one biosample) mapping on identified genes, preferentially mapped on the same region (Table 1 shows the 5' and 3' positions of such ESTs).

Among 275 ESTs from nuclear protein coding genes, 196 ESTs mapped on exons of genes and 76 swinger transformed ESTs mapped on gene introns. Three ESTs mapped partly on gene exons and partly on introns, suggesting that swinger transformations occur before post transcriptional modification, and support our working hypothesis of swinger polymerization during replication/transcription.

The 347 ESTs correspond to 203 protein coding genes and 19 RNAs from non-coding regions. The median and mean size of a protein coding gene in human genome are ~26,288 bp and ~66,577 bp (Piovesan et al., 2016). The mean and median length of identified 203 swinger transformed protein coding genes (including mitochondrial genes) are 124,427 bp and 55,291 bp which is almost twice the mean and median lengths of human protein coding genes. These genes for which swinger transcripts were detected also included the largest protein coding gene in the human genome i.e. RNA binding fox-1 homolog 1 (RBFOX1). These observations indicate that swinger transformation probably due to polymerase fatigue is not only associated to highly active genes but could also associate with large genes. This could also explain the biased mapping of swinger RNAs on the mid and 3' regions of genes, vs the first third of genes.

## Swingier RNAs From Non-Coding and Uncharacterized Genomic Regions

Seventeen and two identified swinger RNAs were transcribed from 15 non-coding genomic regions and two uncharacterized genomic regions, respectively (Table 3). The 15 non-coding
swinger RNAs were not from 5' or 3' UTR or intronic regions of protein-coding genes. To test if these 15 identified genomic regions on various human chromosomes are transcriptionally active, we did a transcriptome search in seventy-one samples (SRX768406-SRX768476, Garzon et al., 2014) with the identified swinger RNAs transcribed from these non-coding genomic regions. Good and complete alignments of the reads on the searched sequences confirmed transcriptions of these 15 non-coding regions (Supplementary Table 2). This result suggests that these 15 regions are transcriptionally active, with potential roles in nerve cells and possible associations with carcinogenesis. However, these specific cancer tissue samples are unavailable, preventing further *in-vitro* and *in-vivo* tests and analyses.

The identified swinger RNAs were mapped on human nuclear chromosomes. No compartmentalizations on any chromosomal arm, region or band of a chromosome were detected in relation to swinger expressed regions. Swinger RNAs produced by transformations that are more conservative at the amino acid level are also more abundant than swinger RNAs produced by transformations that cause more amino acid changes (Seligmann, 2018). This effect was also observed when comparing transformations involving the same number of nucleotides in the transformation (2, 3 or 4). For the distribution of swinger RNAs detected in the present study, too few classes of transformations were detected to enable such subdivisions. However, when considering all transformations, conservation increases with swinger RNA abundances (r = 0.4187, one tailed P = 0.0234). C ↔ G is the least conservative among the transformations involving only two nucleotides, and is one among two transformations for which no swinger RNA was detected in the currently examined data.

## CONCLUSION

We report genomic origins of 347 previously unidentified ESTs generated by the FAPESP/LICR Human Cancer Genome Project (Neto et al., 2000). These represent 0.23% of the 149,500 unidentified ESTs from that project. Note that other types of systematic transformations apparently produce non-canonical RNAs, such as systematic deletions, which produce so-called delRNAs (Seligmann, 2015c; Seligmann, 2016f; El Houmami and Seligmann, 2017; Seligmann, 2017b; Warthi and Seligmann, 2019) and corresponding peptides (Seligmann, 2015c; Seligmann, 2016a; Seligmann, 2016b; Seligmann, 2016c; Seligmann, 2016d; Seligmann, 2016e), including chimeric peptides (Seligmann and Warthi, 2019). This underlines a general approach to identify unknown sequences generated by various sequencing methods. This study identifies swinger RNAs transcribed from multiple cancer-associated genes (Tables 1 and 2), suggesting highly active genes produce swinger transcripts possibly due to enzyme fatigue and promoting cancer progression. The identified sequences

might be sequencing artefacts. However, random artefacts generated by sequencing equipment could not produce systematic exchanges. Systematic sequencing errors should produce swinger RNAs for all sequenced ESTs, however, biosamples where swinger RNAs are detected include regular canonical RNAs. Analyses also identify transcriptionally active non-coding regions on human chromosomes discovering putative ncRNA transcribing regions with potentially significant roles in normal and cancer cells. We provide a unique method to study and identify unknown sequencing reads, reducing loss of important genetic information in raw sequence data. Systematic editing of RNA might contribute to solve the dark DNA conundrum.

## DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the article/ Supplementary Material.

## AUTHOR CONTRIBUTIONS

GW conducted design, analysis, observation and writing. P-EF contributed to result assessments. HS designed the systematic nucleotide exchange framework. P-EF and HS contributed to writing and supervising this study.

## FUNDING

This work was supported by the A\*MIDEX project (no ANR-11-IDEX-0001-02) funded by the « Investissements d'Avenir » French Government program, managed by the French National Research Agency (ANR).

## ACKNOWLEDGMENTS

We thank the reviewers for helpful and constructive comments on earlier draft of the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00042/full#supplementary-material

**SUPPLEMENTARY TABLE 1** | Swinger transformed version of 347 ESTs mapped on human genome.

**SUPPLEMENTARY TABLE 2** | Whole cell transcriptome BLAST result for swinger transformed ESTs mapped on non coding genomic region.
REFERENCES

Ahmed, S. F., Deb, S., Paul, I., Chatterjee, A., Mandal, T., Chatterjee, U., et al. (2012). The chaperone-assisted E3 ligase C terminus of Hsc70-interacting protein (CHIP) targets PTEN for proteasomal degradation. *J. Biol. Chem.* 287, 15996-16006. doi: 10.1074/jbc.M111.321083

Béguin, P., Gill, S., Charpin, N., and Forterre, P. (2015). Synergistic template-free synthesis of dsDNA by Thermococcus nautilus primase PolrTN2, DNA polymerase PolB, and pTN2 helicase. *Extremophiles* 19 (1), 69–76. doi: 10.1007/s00792-014-0706-1

Bahn, J. H., Lee, J. H., Li, G., Greer, C., Peng, G., and Xiao, X. (2012). Accurate identification of A-to-I RNA editing in human by transcriptome sequencing. *Genome Res.* 22, 142–150. doi: 10.1101/gr.124107.111

Bar-Yaacov, D., Avital, G., Levin, L., Richards, A. L., Hachen, N., Rebolledo Jaramillo, B., et al. (2013). RNA-DNA differences in human mitochondria restore ancestral form of 16S ribosomal RNA. *Genome Res.* 23, 1789–1796. doi: 10.1101/gr.161265.113

Bass, B. L. (2002). RNA editing by adenosine deaminases that act on RNA. *Annu. Rev. Biochem.* 71, 817–846. doi: 10.1146/annurev.biochem.71.110601.135501

Blank, A., Gallant, J. A., Burgess, R. R., and Loeb, L. A. (1986). An RNA polymerase mutant with reduced accuracy of chain elongation. *Biochemistry* 25, 5920–5928. doi: 10.1021/bi00368a013

Bullerdiek, J., and Rommel, B. (2018). Factors targeting MED12 to drive tumorigenesis? *F1000Research* 7, 359. doi: 10.12688/f1000research.14227.1

Chen, C., and Bundschuh, R. (2012). Systematic investigation of insertional and deletional RNA-DNA differences in the human transcriptome. *BMC Genomics* 13, 616. doi: 10.1186/1471-2164-13-616

Ciocca, D. R., and Calderwood, S. K. (2005). Heat shock proteins in cancer: diagnostic, prognostic, predictive, and treatment implications. *Cell Stress Chaperones* 10 (2), 86–103. doi: 10.1379/CSC-99r.1

Collavin, L., Lazarević, D., Utrera, R., Marzinotto, S., Monte, M., and Schneider, C. (1999). wt p53 dependent expression of a membrane-associated isoform of adenylate kinase. *Oncogene* 18, 5879–5888. doi: 10.1038/sj.onc.1202970

Dunham, I., Kundaje, A., Aldred, S. F., Collins, P. J., Davis, C. A., Doyle, F., et al. (2012). An integrated encyclopedia of DNA elements in the human genome. *Nature* 489 (7414), 57–74. doi: 10.1038/nature11247

Eisermann, K., Dar, J. A., Dong, J., Wang, D., Masoodi, K. Z., and Wang, Z. (2015). Poly (A) binding protein cytoplasmic 1 is a novel co-regulator of the androgen receptor. *PloS One* 10 (7), e0128495. doi: 10.1371/journal.pone.0128495

El Houmami, N., and Seligmann, H. (2017). Evolution of nucleotide punctuation marks: from structural to linear signals. *Front. Genet.* 27, 8. doi: 10.3389/fgene.2017.00036

Fagerberg, L., Hallstrom, B. M., Oksvold, P., Kampf, C., Djureinovic, D., Odeberg, J., et al. (2014). Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. *Mol. Cell. Proteomics* 13(2), 397–406. doi: 10.1074/mcp.M113.035600

Garzon, R., Volinia, S., Papaioannou, D., Nicolet, D., Kohlschmidt, J., Yan, P. S., et al. (2014). Expression and prognostic impact of lncRNAs in acute myeloid leukemia. *Proc. Natl. Acad. Sci. U. S. A.* 111 (52), 18679–18684. doi: 10.1073/pnas.1422050112

Gaude, H., Aznar, N., Delay, A., Bres, A., Buchet-Poyau, K., Caillat, C., et al. (2012). Molecular chaperone complexes with antagonizing activities regulate stability and activity of the tumor suppressor LKB1. *Oncogene* 31, 1582–1591. doi: 10.1038/onc.2011.342

Hatakeyama, S., Watanabe, M., Fujii, Y., and Nakayama, K. I. (2005). Targeted destruction of c-Myc by an engineered ubiquitin ligase suppresses cell transformation and tumor formation. *Cancer Res.* 65, 7874–7879. doi: 10.1158/0008-5472.CAN-05-1581

Jan, Y. H., Lai, T. C., Yang, C. J., Huang, M. S., and Hsiao, M. (2019). A co-expressed gene status of adenylate kinase 1/4 reveals prognostic gene signature associated with prognosis and sensitivity to EGFR targeted therapy in lung adenocarcinoma. *Sci. Rep.* 9 (1), 12329. doi: 10.1038/s41598-019-48243-9

Janssen, E., Kuiper, J., Hodgson, D., Zingman, L. V., Alekseev, A. E., Terzic, A., et al. (2004). Two structurally distinct and spatially compartmentalized adenylate kinases are expressed from the AK1 gene in mouse brain. *Mol. Cell. Biochem.* 256-257, 59–72. doi: 10.1023/B:MCBI.000000985.15267.db

Kajiro, M., Hirota, R., Nakajima, Y., Kawanowa, K., So-Ma, K., Ito, I., et al. (2009). The ubiquitin ligase CHIP acts as an upstream regulator of oncogenic pathways. *Nat. Cell Biol.* 11, 312–319. doi: 10.1038/ncb1839

Kandel, E. S., and Nudler, E. (2002). Template switching by RNA polymerase II *in vivo*: evidence and implications from a retroviral system. *Mol. Cell.* 10 (6), 1495–1502. doi: 10.1016/S1097-2765(02)00777-3

Knippa, K., and Peterson, D. O. (2013). Fidelity of RNA polymerase II transcription: role of Rbp9 in error detection and proofreading. *Biochemistry* 52, 7807–7817. doi: 10.1021/bi4009566

Księżakowska-Lakoma, K., Żyla, M., and Wilczyński, J. R. (2014). Mitochondrial dysfunction in cancer. *Przeglad menopauzalny = Menopause Rev.* 13 (2), 136–144. doi: 10.5114/pm.2014.42717

Kumar, S., Razzaq, S. K., Vo, A. D., Gautam, M., and Li, H. (2016). Identifying fusion transcripts using next generation sequencing. *Wiley Interdiscip. Rev. RNA* 7 (6), 811–823. doi: 10.1002/wrna.1382

López-Nieva, P., Fernández-Navarro, P., Graña-Castro, O., Andrés-León, E., Santos, J., Villa-Morales, M., et al. (2019). Detection of novel fusion-transcripts by RNA-Seq in T-cell lymphoblastic lymphoma. *Sci. Rep.* 9, 5179. doi: 10.1038/s41598-019-41675-3

Löytynoja, A., and Goldman, N. (2017). Short template switch events explain mutation clusters in the human genome. *Genome Res.* 27 (6), 1039–1049. doi: 10.1101/gr.214973.116

Lee, H. R., and Johnson, K. A. (2006). Fidelity of the human mitochondrial DNA polymerase. *J. Biol. Chem.* 281, 36236–36240. doi: 10.1074/jbc.M607964200

Lee, J. A., Carvalho, C. M. B., and Lupski, J. R. (2007). A DNA replication mechanism for generating nonrecurrent rearrangements associated with genomic disorders. *Cell.* 131(7), 1235–1247. doi: 10.1016/j.cell.2007.11.037

Lee, S. Y., Joung, J. G., Park, C. H., Park, J. H., and Kim, J. H. (2015). RCARE: RNA Sequence Comparison and Annotation for RNA Editing. *BMC Med. Genomics* 8 (2), S8. doi: 10.1186/1755-8794-8-S2-S8

Li, M., Wang, I. X., Li, Y., Bruzel, A., Richards, A. L., Toung, J. M., et al. (2011). Widespread RNA and DNA sequence differences in the human transcriptome. *Science* 333 (6038), 53–58. doi: 10.1126/science.1207018

Li, L., Zeng, Q., Bhutkar, A., Galván, J. A., Karamitopoulou, E., Noordermeer, D., et al. (2018). GKAP acts as a genetic modulator of NMDAR signaling to govern invasive tumor growth. *Cancer Cell* 33 (4), 736–751. doi: 10.1016/j.ccell.2018.02.011

Lin, C. S., Liu, L. T., Ou, L. H., Pan, S. C., Lin, C. I., and Wei, Y. H. (2018). Role of mitochondrial function in the invasiveness of human colon cancer cells. *Oncol. Rep.* 39 (1), 316–330. doi: 10.3892/or.2017.6087

Liu, D., Yin, B., Wang, Q., Ju, W., Chen, Y., Qiu, H., et al. (2012). Cytoplasmic Poly(A) binding protein 4 is highly expressed in human colorectal cancer and correlates with better prognosis. *J. Genet. Genomics.* 39 (8), 369–374 doi: 10.1016/j.jgg.2012.05.007

Ma, Y., Chen, Y., and Petersen, I. (2017). Expression and epigenetic regulation of cystatin B in lung cancer and colorectal cancer. *Pathol. Res. Pract.* 213 (12), 1568–1574. doi: 10.1016/j.prp.2017.06.007

Michel, C. J., and Seligmann, H. (2014). Bijective transformation circular codes and nucleotide exchanging RNA transcription. *BioSystems.* 118 (1), 39–50. doi: 10.1016/j.biosystems.2014.02.002

Modica-Napolitano, J. S., and Singh, K. K. (2004). Mitochondrial dysfunction in cancer. *Mitochondrion.* 4 (5–6 SPEC. ISS.), 755–762. doi: 10.1016/j.mito.2004.07.027

Morgulis, A., Coulouris, G., Raytselis, Y., Madden, T. L., Agarwala, R., and Schäffer, A. A. (2008). Database indexing for production MegaBLAST searches. *Bioinformatics* 24 (16), 1757–1764. doi: 10.1093/bioinformatics/btn322

Neto, E. D., Correa, R. G., Verjovski-Almeida, S., Briones, M. R. S., Nagai, M. A., Da Silva, W., et al. (2000). Shotgun sequencing of the human transcriptome with ORF expressed sequence tags. *Proc. Natl. Acad. Sci. U. S. A.* 97, 3491–3496. doi: 10.1073/pnas.97.7.3491

Ninio, J. (1991). Connections between translation, transcription and replication error-rates. *Biochimie.* 73, 1517–1523. doi: 10.1016/0300-9084(91)90186-5

Pareja, F., Da Cruz Paula, A., Murray, M. P., Hoang, T., Gularte-Mérida, R., Brown, D., et al. (2019). Recurrent MED12 exon 2 mutations in benign breast fibroepithelial lesions in adolescents and young adults. *J. Clin. Pathol.* 72 (3), 258–262. doi: 10.1136/jclinpath-2018-205570
Park, E., Williams, B., Wold, B. J., and Mortazavi, A. (2012). RNA editing in the human ENCODE RNA-seq data. *Genome Res.* 22, 1626–1633. doi: 10.1101/ gr.134957.111

Piovesan, A., Caracausi, M., Antonaros, F., Pelleri, M. C., and Vitale, L. (2016). GeneBase 1.1: a tool to summarize data from NCBI gene datasets and its application to an update of human gene statistics. *Database J. Biol. Database Curation.* 2016, 1–13. doi: 10.1093/database/baw153

Porath, H. T., Carmi, S., and Levanon, E. Y. (2014). A genome-wide map of hyper-edited RNA reveals numerous new sites. *Nat. Commun.* 5, 4726. doi: 10.1038/ ncomms5726

Ruckova, E., Muller, P., Nenutil, R., and Vojtesek, B. (2012). Alterations of the Hsp70/Hsp90 chaperone and the HOP/CHIP co-chaperone system in cancer. *Cell. Mol. Biol. Lett.* 17, 446–458. doi: 10.2478/s11658-012-0021-8

Schaub, M., and Keller, W. (2002). RNA editing by adenosine deaminases generates RNA and protein diversity. *Biochimie.* 84 (8), 791–803. doi: 10.1016/S0300-9084 (02)01446-3

Seligmann, H., and Raoult, D. (2018). Stem-loop RNA hairpins in giant viruses: Invading rRNA-like repeats and a template free RNA. *Front. Microbiol.* 9, 101. doi: 10.3389/fmicb.2018.00101

Seligmann, H., and Warthi, G. (2019). Chimeric translation for mitochondrial peptides: regular and expanded codons. *Comput. Struct. Biotech.* J. 17, 1195–1202. doi: 10.1016/j.csbj.2019.08.006

Seligmann, H. (2012). Overlapping genes coded in the 3'-to-5'-direction in mitochondrial genes and 3'-to-5' polymerization of non-complementary RNA by an 'invertase'. *J. Theor. Biol.* 315, 38–52. doi: 10.1016/j.jtbi. 2012.08.044

Seligmann, H. (2013a). "Replicational Mutation Gradients, Dipole Moments, Nearest Neighbour Effects and DNA Polymerase Gamma Fidelity in Human Mitochondrial Genomes," in *The Mechanisms of DNA Replication*. Ed. D. Stuart (London: InTech, chpater10), 257–286. doi: 10.5772/51245

Seligmann, H. (2013b). Polymerization of non-complementary RNA: Systematic symmetric nucleotide exchanges mainly involving uracil produce mitochondrial RNA transcripts coding for cryptic overlapping genes. *BioSystems* 111, 156–174. doi: 10.1016/j.biosystems.2013.01.011

Seligmann, H. (2013c). Systematic asymmetric nucleotide exchanges produce human mitochondrial RNAs cryptically encoding for overlapping protein coding genes. *J. Theor. Biol.* 324, 1–20. doi: 10.1016/j.jtbi.2013.01.024

Seligmann, H. (2014a). Species radiation by DNA replication that systematically exchanges nucleotides? *J. Theor. Biol.* 363, 216–222. doi: 10.1016/ j.jtbi.2014.08.036

Seligmann, H. (2014b). Mitochondrial swimger replication: DNA replication systematically exchanging nucleotides and short 16S ribosomal DNA swimger inserts. *BioSystems* 125, 22–31. doi: 10.1016/j.biosystems.2014.09.012

Seligmann, H. (2015a). Swinger RNAs with sharp switches between regular transcription and transcription systematically exchanging ribonucleotides: Case studies. *BioSystems* 135, 1–8. doi: 10.1016/j.biosystems.2015.07.003

Seligmann, H. (2015b). Sharp switches between regular and swimger mitochondrial replication: 16S rDNA systematically exchanging nucleotides A<-> T+C<-> G in the mitogenome of *Kamimuria wangi*. *Mitochondrial DNA* 27 (4), 2440–2446. doi: 10.3109/19401736.2015.1033691

Seligmann, H. (2015c). Codon expansion and systematic transcriptional deletions produce tetra-, pentacoded mitochondrial peptides. *J. Theor. Biol.* 387, 154–165. doi: 10.1016/j.jtbi.2015.09.030

Seligmann, H. (2016a). Natural chymotrypsin-like-cleaved human mitochondrial peptides confirm tetra-, pentacodon, non-canonical RNA translations. *BioSystems* 147, 78–93. doi: 10.1016/j.biosystems.2016.07.010

Seligmann, H. (2016b). Translation of mitochondrial swimger RNAs according to tri-, tetra- and pentacodons. *BioSystems* 140, 38–48. doi: 10.1016/ j.biosystems.2015.11.009

Seligmann, H. (2016c). Unbiased mitoproteome analyses confirm non-canonical RNA, expanded codon translations. *Comput. Struct. Biotechnol.* J. 14, 391–403. doi: 10.1016/j.csbj.2016.09.004

Seligmann, H. (2016d). Chimeric mitochondrial peptides from contiguous regular and swimger RNA. *Comput. Struct. Biotechnol.* J. 14, 283–297. doi: 10.1016/ j.csbj.2016.06.005

Seligmann, H. (2016e). Swinger RNA self-hybridization and mitochondrial non-canonical swimger transcription, transcription systematically exchanging nucleotides. *J. Theor. Biol.* 399, 84–91. doi: 10.1016/j.jtbi.2016.04.007

Seligmann, H. (2016f). Systematically frameshifting by deletion of every 4th or 4th and 5th nucleotides during mitochondrial transcription: RNA self-hybridization regulates delRNA expression. *Biosystems* 142-143, 43–51. doi: 10.1016/j.biosystems.2016.03.009

Seligmann, H. (2017a). Natural mitochondrial proteolysis confirms transcription systematically exchanging/deleting nucleotides, peptides coded by expanded codons. *J. Theor. Biol.* 414, 76–90. doi: 10.1016/j.jtbi.2016.11.021

Seligmann, H. (2017b). Reviewing evidence for systematic transcriptional deletions, nucleotide exchanges, and expanded codons, and peptide clusters in human mitochondria. *BioSystems* 160, 10–24. doi: 10.1016/ j.biosystems.2017.08.002

Seligmann, H. (2018). Bijective codon transformations show genetic code symmetries centered on cytosine's coding properties. *Theory Biosci* 137, 17–31. doi: 10.1007/s12064-017-0258-x

Singh, A., Zahra, S., Das, D., and Kumar, S. (2019). "AtFusionDB: a database of fusion transcripts in Arabidopsis thaliana. *Database* 2019, 1–9. doi: 10.1093/ database/bay135

Strathern, J. N., Jin, D. J., Court, D. L., and Kashlev, M. (2012). Isolation and characterization of transcription fidelity mutants. *Biochim. Biophys. Acta - Gene Regul. Mech.* 1819, 694–699. doi: 10.1016/j.bbagrm.2012.02.005

Tarantul, V., Nikolaev, A., Hannig, H., Kalmyrazaev, B., Muchoyan, I., Maximov, V., et al. (2001). Detection of abundantly transcribed genes and gene translocation in human immunodeficiency virus - Associated non-hodgkin's lymphoma. *Neoplasia* 3 (2), 132–142. doi: 10.1038/ sj.neo.7900137

Vasseur, S., Malicet, C., Calvo, E. L., Dagorn, J. C., and Iovanna, J. L. (2005). Gene expression profiling of tumours derived from rasv12/E1A-transformed mouse embryonic fibroblasts to identify genes required for tumour development. *Mol. Cancer* 4 (1), 4. doi: 10.1186/1476-4598-4-4

Wang, S., Wu, X., Zhang, J., Chen, Y., Xu, J., Xia, X., et al. (2013). CHIP functions as a novel suppressor of tumour angiogenesis with prognostic significance in human gastric cancer. *Gut* 62, 496–508. doi: 10.1136/gutjnl- 2011-301522

Wang, I. X., Core, L. J., Kwak, H., Brady, L., Bruzel, A., McDaniel, L., et al. (2014). RNA-DNA differences are generated in human cells within seconds after RNA exits polymerase II. *Cell Rep.* 6, 906–915. doi: 10.1016/ j.celrep.2014.01.037

Wang, T., Yang, J., Xu, J., Li, J., Cao, Z., Zhou, L., et al. (2014a). CHIP is a novel tumor suppressor in pancreatic cancer through targeting EGFR. *Oncotarget.* 5, 1969–1986. doi: 10.18632/oncotarget.1890

Wang, Y., Ren, F., Wang, Y., Feng, Y., Wang, D., Jia, B., et al. (2014b). CHIP/Stub1 functions as a tumor suppressor and represses NF-κB-mediated signaling in colorectal cancer. *Carcinogenesis* 35, 983–991. doi: 10.1093/carcin/bgt393

Wang, I. X., Grunseich, C., Chung, Y. G., Kwak, H., Ramrattan, G., Zhu, Z., et al. (2016). RNA-DNA sequence differences in Saccharomyces cerevisiae. *Genome Res.* 26, 1544–1554. doi: 10.1101/gr.207878.116

Warthi, G., and Seligmann, H. (2018). "Swinger RNAs in the Human Mitochondrial Transcriptome," in *In mitochondrial DNA-new insights*. Eds. H. Seligmann and G. Warthi, (London) 79–92. Chapter 4. doi: 10.5772/ intechopen.80805

Warthi, G., and Seligmann, H. (2019). Transcripts with systematic nucleotide deletion of 1-12 nucleotide in human mitochondrion suggest potential non-canonical transcription. *PloS One* 14, e0217356. doi: 10.1371/journal.pone.0217356

Warthi, G., Fournier, PE., and Seligmann, H. (2019). Identification of noncanonical transcripts produced by systematic nucleotide exchanges in HIV-associated centroblastic lymphoma. *DNA Cell Biol.* doi: 10.1089/ dna.2019.5066 [Epub ahead of print]

Xie, J., Ubango, J., Ban, Y., Chakravarti, D., Kim, J. J., and Wei, J. J. (2018). Comparative analysis of AKT and the related biomarkers in uterine leiomyomas with MED12, HMGA2, and FH mutations. *Genes Chromosom. Cancer* 57 (10), 485–494. doi: 10.1002/gcc.22643

Ying, Z., Haiyan, G., and Haidong, G. (2013). BAG5 regulates PTEN stability in MCF-7 cell line. *BMB Rep.* 46, 490–494. doi: 10.5483/BMBRep.2013.46. 10.268

Zeng, M., Li, F., Wang, L., Chen, C., Huang, X., Wu, X., et al. (2018). Downregulated cytoplasmic polyadenylation element-binding protein-4 is associated with the carcinogenesis of head and neck squamous cell carcinoma. *Oncol. Lett.* 15 (3), 3226–3232. doi: 10.3892/ol.2017.7661
Zhang, Z., Schwartz, S., Wagner, L., and Miller, W. (2000). A greedy algorithm for aligning DNA sequences. *J. Comput. Biol.* 7 (1–2), 203–214. doi: 10.1089/10665270050081478

Zhang, J., Shi, Z. F., Huang, J. X., and Zou, X. G. (2016). CSTB downregulation promotes cell proliferation and migration and suppresses apoptosis in gastric cancer SGC-7901 cell line. *Oncol. Res.* 24 (6), 487–494. doi: 10.3727/096504016X14685034103752

Zhou, Y. N., Lubkowska, L., Hui, M., Court, C., Chen, S., Court, D. L., et al. (2013). Isolation and characterization of RNA polymerase rpoB mutations that alter transcription slippage during elongation in Escherichia coli. *J. Biol. Chem.* 288, 2700–2710. doi: 10.1074/jbc.M112.429464

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Warthi, Fournier and Seligmann. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.