All Roads Lead to (Ch)Romosome 19!

A hepatocellular carcinoma (HCC) co-regulatory network exists between chromosome 19 microRNA cluster (C19MC) at 19q13.42, melanoma-A antigens, IFN-γ and p53, promoting an oncogenic role of C19MC that is disrupted by metal ions zinc and nickel. IFN-γ plays a co-operative role whereas IL-6 is antagonistic, each have a major bearing on the expression of HLA molecules on cancer cells. Analysis of Mesenchymal stem cells and cancer cells predicted C19MC modulation of apoptosis in induced pluripotency and tumorigenesis.

Key, differentially expressed genes in HCC included cancer-related transcription factors (TF) EGR1, FOS, and FOSB. From mRNA and miRNA expression profiles these were most enriched in the p53 signaling pathway where mRNA levels of each decreased in HCC tissues. In addition, mRNA levels of CCNB1, CCNB2, and CHEK1, key markers of the p53 signaling pathway, were all increased. miR-181a-5p regulated FOS and EGR1 to promote the invasion and progression of HCC by p53 signaling pathway and it plays an important role in maturation or impairment of natural killer (NK) cells.

A pan-cancer analysis, on microRNA-associated gene activation, produced the top 57 miRNAs that positively correlated with at least 100 genes. miR-150, at 19q13.33 was the most active, it positively correlated with 1009 different genes each covering at least 10 cancers. It is an important hematopoietic, especially B, T, and NK, cell specific miRNA.

Rapid functional impairment of NK cells following tumor entry limits anti-tumor immunity. Gene regulatory network analysis revealed downregulation of TF regulons, over pseudo-time, as NK cells transition to their impaired end state. These included AP-1 complex TF's, Fos, Fosb (19q13.32), Jun, Junb (19p13.13), which are activated during NK cell cytolytic programs and down regulated by interactions with inhibitory ligands. Other down-regulated TF's included Irf8, Klf2 (19p13.11), Myc, which support NK cell activation and proliferation. There were no significantly upregulated TF's suggesting that the tumor-retained NK state arises from the reduced activity of core transcription factors associated with promoting mature NK cell development and expansion.

Innate immune, intra-tumoral, stimulatory dendritic cells (SDCs) and NK cells cluster together and are necessary for enhanced T cell tumor responses. In human melanoma, SDC abundance is associated with intra-tumoral expression of the cytokine producing gene FLT3LG (19q13.33) that is predominantly produced by NK cells in tumors. Computed tomography exposes patients to ionizing X-irradiation. Determined trends in the expression of 24 radiation-responsive genes linked to cancer, in vivo, found that TP53 and FLT3LG expression increased linearly with CT dose. 

Undifferentiated embryonal sarcoma of the liver displays high aneuploidy with recurrent alterations of 19q13.4 that are uniformly associated with aberrantly high levels of transcriptional activity of C19MC microRNA. Further, TP53 mutation or loss was present with all samples that also display C19MC changes. The 19q13.4 locus is gene-poor with highly repetitive sequences. Given the noncoding nature and lack of an obvious oncogene, disruption of the nearby C19MC regulatory region became a target for tumorigenesis. 

The endogenous retroviral, hot-spot deletion rate at 19p13.11-19p13.12 and 19q33-19q42 occurs at double the background deletion rate. Clustered in and around these regions are many gene families including KIR, Siglec, Leukocyte immunoglobulin-like receptors and cytokines that associate important NK gene features to proximal NK genes that were overrepresented in a meta analysis of blood pressure. 

Endogenous retroviruses that invite p53 and its transcriptional network, at retroviral hot-spots, suggest that lymphocyte progenitors, such as ILC's and expanded, NK cells are synergistically responsive to transcription from this busy region including by the top differentially expressed blood pressure genes MYADM, GZMB, CD97, NKG7, CLC, PPP1R13L , GRAMD1A as well as (RAS-KKS) Kallikrein related peptidases to educate early and expanded NK cells that shape immune responses.  

Can Ancient Pathways Defeat Cancer?

It has been widely acknowledged that non-coding RNAs are master-regulators of genomic function. The association between human introns and ncRNAs has a pronounced synergistic effect with important implications for fine-tuning gene expression patterns across the entire genome. There is also strong preference of ncRNA from intronic regions particularly associated with the transcribed strand. 

Accumulating evidence demonstrates that, analogous to other small ncRNAs (e.g. miRNAs, siRNA's etc.) piRNAs have both oncogenic and tumor suppressive roles in cancer development. Functionally, piRNAs maintain genomic integrity and cell age by silencing repetitive, transposable elements, and are capable of regulating the expression of specific downstream target genes in a post-transcriptional manner. 

Unlike miRNAs and siRNAs, the precursors of piRNAs are single stranded transcripts without any prominent secondary hairpin structures. These precursors are usually generated from specific genomic locations containing repetitive elements, a process that is typically orchestrated via a Dicer-independent pathway. 

Without restraint, the ancient, L1 class of transposable elements can interrupt the genome through insertions, deletions, rearrangements, and copy number variations. L1 activity has contributed to instability and evolution of genomes, and is tightly regulated by DNA methylation, histone modifications, and piRNA. They can impact genome variation by mispairing and unequal crossing-over during meiosis due to repetitive DNA sequences. Indeed meiotic double-strand breaks are the proximal trigger for retrotransposon eruptions as highlighted in animals lacking p53.

Through a novel 28-base small piRNA of the KIR3DL1 gene, antisense transcripts mediate Killer Ig-like receptor (KIR) transcriptional silencing in immune somatic, Natural Killer (NK) cell lineage, a mechanism that may be broadly used in orchestrating immune development. Expressed on NK cells, KIR's are important determinants of NK cell function. Silencing  individual KIR genes is strongly correlated with the presence of CpG dinucleotide methylation within the promoter. 

Structural research exposed the enormous binding complexity behind KIR haplotypes and HLA allotypes. Not only via protein structures, but also plasticity and selective binding behavior's as influenced by extrinsic factors. One study links a specific recognition of HLA-C*05:01 by KIR2DS4 receptor through a peptide highly conserved among bacteria pathogenic in humans. Another demonstrated a hierarchy of functional peptide selectivity by KIR–HLA-C interactions, including cross-reactive binding, with relevance to NK cell biology and human disease associations. Additionally a p53 peptide most overlapped other high performance peptides for a HLA-C allotype C*02:02 that shares identical contact residues with C*05:01.

Ancient pathways linking p53 to attenuation of aberrant stem cell proliferation may predate the divergence between vertebrates and invertebrates. Human stem cell proliferation, as determined by p53 transposable element silencing, may also serve a NK progenitor to promote the repertoire of more than 30,000 NK cell subsets. 

A recent study showed that wild type p53 can restrain transposon mobility through interaction with PIWI-piRNA complex. Also, cellular metabolism regulates sensitivity to NK cells depending on P53 status and P53 pathway is coupled to NK cell maturation leaving open the possibility that a direct relationship exists. Further, functional interactions between KIR and HLA modify risks of basal cell carcinoma (BCC) and squamous cell carcinomas (SCC) and KIR B haplotypes provide selective pressure for altered P53 in BCC tumors. 

Anticipating p53's broader influences or responses, cells, extracted from 48 different sections of 7 tumor biopsies were sequenced and TP53 DNA computed using Codondex algorithm. Each section produced a TP53 Consensus Variant (CV), represented by its intron1, ncDNA Key Sequence's (KS). Bioinformatic correlations between each KS and cytotoxicity resulting from NK coculture with the section may predict KIR-HLA and extrinsic factor plasticity to reliably determine from KS's, optimal cell/tissue selections for NK cell education and licensing. 

Tolerating Your Non-self!

Immune cells get comfortable with cancer
Courtesy https://deepai.org

A hallmark of cancer, autoimmunity and disease is the aberrant transcription of typically silenced, repetitive genetic elements that mimic Pathogen-Associated Molecular Patterns (PAMP's) that bind Pattern Recognition Receptors (PPR's) triggering the innate immune system and inflammation. Unrestrained, this 'viral mimicry' activates a generally conserved mechanism that, under restraint, supports homeostasis. These repetitive viral DNA sequences normally act as a quality control over genomic dysregulation responding in ways that preferentially promote immune conditions for stability. If aberrantly unrestrained and the 'viral mimicry' is transcribed it may result in undesirable immune reactions that disrupt the homeostasis of cells.

Mitochondrial DNA (mtDNA) are one source of cytosolic double stranded RNA (dsRNA) that is commonly present in cells. Trp53 Mutant Embryonic Fibroblasts (MEF's) contain innate immune stimulating endogenous dsRNA, from mtDNA that mimic PAMP's. The immune response, via RIG-1 like PRR, leads to expression of type 1 interferon (IFN) and proinflammatory cytokine genes. Further, Natural Killer cells also produce a multitude of cytokines that can promote or dampen an immune response. Wild-type p53 suppresses viral repeats and contributes to innate immunity by enhancing IFN-dependent antiviral activity independent of its function as a proapoptotic and tumor suppressor gene. 

Post-translationally modified P53, located in the cytoplasm, enhances the permeability of the mitochondrial outer membrane thus stimulating apoptosis. However, treating Trp53 mutant MEF's with DNA demethylating agent caused a huge increase in the level of transcripts encoding short interspersed nuclear elements and other species of noncoding RNAs that generated a strong type 1 IFN response. This did not occur in p53 wild-type MEF's. Thus it appears that another function of p53 is to silence repeats that can accidentally induce an immune response.

This has several implications for how we understand self versus non-self discrimination. When pathogen-associated features were quantified, specific repeats in the genome not only display PAMP's capable of stimulating PRRs but, in some instances, have seemingly maintained such features under selection. For organisms with a high degree of epigenetic regulation and chromosomal organization immuno-stimulatory repeats release a danger signal, such as repeats released after p53 mutations. Here, immune stimulation may act as back-up for the failure of other p53 functions such as apoptosis or senescence due to mutation. This supports the hypothesis that specific repeats gained favor by maintaining non-self PAMPs to act as sensors for loss of heterochromatin as an epigenetic checkpoint of quality control that avoids genome instability generally. 

When P53 mutates it begins to fail its restraint of viral suppression, this enables a 'viral mimicry' and aberrant immune reactions. These may promote survival of cells that can leverage immunity, promote angiogenesis and heightened proliferation of cancers, or other diseases under modified conditions for non-self tolerance. 

Toward Customized Natural Killer Cells

An important role of Natural Killer (NK) cells is to eliminate other cells that extinguish or diminish expression of self-MHC class I molecules or Human Leukocyte Antigen (HLA), which commonly occurs as a result of viral infection or cellular transformation. This capacity arises because NK cells express stimulatory and inhibitory receptors that engage ligands on normal cells. The majority of inhibitory receptors belong to the Killer-cell immunoglobulin-like receptors (KIR) and CD94/NKG2A  families and are specific for MHC I molecules. When an NK cell encounters a normal cell, engagement of the inhibitory receptors conveys signals that counteract stimulatory signaling. Lysis occurs when inhibition is lost because the target cell lacks one or more self-MHC molecules or when target cells express high levels of stimulatory ligands that counter inhibition.

Mitochondrial DNA (MtDNA) embedded in the genomes of 66,000 humans was associated with adverse consequences including cancer. Overall tumor specific nuclear embedded MtDNA was more common on Chromosome (Chr)19, less common on Chr6 and tended to involve non-coding, repetitive elements or satellite repeats. 

The dimorphic relationship between genes on Chr6, encoding HLA and  Chr19, encoding KIRs  may elucidate how, why and when NK cells determine self restraint or attack cells infected by pathogens and disease. Chr19 has also been linked to blood pressure mechanics, immunity and checkpoints associated with P53. Cancer mutation burden is shaped by G4 DNA, cell cycle replication stress, DNA repair pathway and mitochondrial dysfunction. G4 DNA overrepresentation generally occurs in tumors with mutations in tumor suppressor gene's such as TP53. 

Whether KIR-HLA relationships are associated with p53 status of NK cells and of its target is unknown. However, it has been reported that cellular metabolism regulates a cells sensitivity to NK cells depending on its P53 status and that P53 pathway is coupled to NK cell maturation leaving open the possibility that a relationship exists. 

KIR and HLA genes are polymorphic and display significant variations, The independent segregation of these unlinked gene families produces extraordinary diversity in the number and type of KIR-HLA pairs inherited in individuals. Variation affects the KIR repertoire of NK cell clones, NK cell maturation, the capability to deliver signals, and consequently the NK cell response to human diseases.

One study suggests that functional interactions between KIR and HLA modify risks of basal cell carcinoma (BCC) and squamous cell carcinomas (SCC) and that KIR B haplotypes provide selective pressure for altered P53 in BCC tumors.

MtDNA and other insertions into nuclear DNA may have altered Chr19-Chr6 linkage relationships and KIR-HLA validity, affecting the integrity of NK missing-self surveillance. Therefore, P53 dependent metabolism and P53 coupled NK cell education may point to a required synchronicity, obtained through NK education, licensing KIR-HLA and other receptor-ligand combinations for a global NK symbiosis.

The altered landscape of cancer is often characterized by a heterogeneous mix of immunosuppressive metabolites, glucose and amino acid deprivation, hypoxia and acidity, which, in concert, prevent effective anti-tumor immunity, here NK therapies herald great potential.

NK cell co-culture with patient cells selected using precise P53 rankings for a distinct P53-coupled-NK cell education may realize a mature NK subset with P53-paired characteristics. Trojan therapy using autologous or combined allogeneic NK cells may promote licensing, through a broad synchronization including at least KIR-HLA. This ex-vivo approach may resist re-education in vivo and activate against P53-decoupled-KIR-HLA affected cells. The objective is an NK subset that, in vivo will initiate and progress a limited innate immune response and disrupt near-neighbor targets that will contribute to a broader immune response.  

Retroviral Defense And Mitochondrial Offense

Chromosomal DNA has played host to the long game of viral insertions that repeat and continue as a genetic and epigenetic symbiosis along its phosphate and pentose sugar backbone. But, the bacterial origin of mitochondria and its hosted DNA also promotes its offense. 

Research suggests that retrovirus insertions evolved from a type of transposon called a retrotransposon. The evolutionary time scales of inherited, endogenous retroviruses (ERV) and the appearance of the zinc finger gene that binds its unique sequences occur over same time scales of primate evolution. Additionaly the zinc-finger genes that inactivate transposable elements are commonly located on chromosome 19. The recurrence of independent ERV invasions can be countered by a reservoir of zinc-finger repressors that are continuously generated on copy number variant (CNV) formation hotspots.

One of the more intiguing aspects of prevalent CNV hotspots on chromosome 19 are their proximity to killer immunoglobulin receptor gene's (KIR's) and other critical gene's of the innate immune system.

Frequently occuring DNA breaks can cause genomic instability, which is a hallmark of cancer. These breaks are over represented at G4 DNA quadruplexes within, hominid-specific, SVA retrotransposons and generally occur in tumors with mutations in tumor suppressor genes, such as TP53. Cancer mutational burden is shaped by G4 DNA, replication stress and mitochondrial dysfunction, that in lung adenocarcinoma downlregulates SPATA18, a mitochondrial eating protein (MIEAP) that contributes to mitophagy. 

Genetic variations, in non-coding regions can control the activity of conserved protein-coding genes resulting in the establishment of species-specific transcriptional networks. A chromosome 19 zinc finger, ZNF558 evolved as a suppressor of LINE-1 transposons, but has since been co-opted to singly regulate SPATA18. These variations are evident from a panel of 409 human lymphoblastoid cell lines where the lengths of the ZNF558 variable number tandem repeats (VNTR) negatively correlated with its expression. 

Colon cancer cells with p53 deletion were used to analyze deregulated p53 target genes in HCT116 p53 null cells compared to HCT116-p53 +/+ cells. SPATA18 was the most upregulted gene in the differential expression providing further insight to p53 and mitophagy via SPATA18-MIEAP.

p53 response elements (p53RE) can be shaped by long terminal repeats from endogenous retroviruses, long interspersed nuclear repeats, and ALU repeats in humans and fuzzy tandem repeats in mice. Further, p53 pervasively binds to p53REs derived from retrotransposons or other mobile genetic elements and can suppress transcription of retroelements. The p53- mediated mechanisms conferring protection from retroelements is also conserved through evolution. Certainly, p53 has been shown to have other roles in DNA  context, such as playing an important role in replication restart and replication fork progression. The absence of these p53-dependent processes can lead to further genomic instability. 

The frequency of variable length, long or short nucleotide repeats and their locations within a gene may be key to the repression of DNA sequences that would otherwise cause genomic instability or protein expressions that would eat bacterial mitochondria or destroy its cell host. 

The complexity of variable length insertions is made evident when exhaustively analyzing a simple length 12 sequence for the potential frequency of each of its variable length repeats starting from a minumum variable length of 8.

Then, for TGTGGGCCCACA(12)

All possible internal variable length combinations from and including length 8:

TGTGGGCC(8)|GTGGGCCC(8)|TGTGGGCCC(9)|TGGGCCCA(8)|GTGGGCCCA(9)|TGTGGGCCCA(10|GGGCCCAC(8)|TGGGCCCAC(9)|GTGGGCCCAC(10)|TGTGGGCCCAC(11)|GGCCCACA(8)|GGGCCCACA(9)|TGGGCCCACA(10)|GTGGGCCCACA(11)|TGTGGGCCCACA(12)

For example, reviewing length (8) only:

TGTGGGCC (8) occurs 5 times

GTGGGCCC (8) occurs 8 times

TGGGCCCA (8) occurs 9 times

GGGCCCAC (8) occurs 8 times

GGCCCACA (8) occurs 5 times

Any repeat can be ranked based on its ocurrence within all possible combinations of a given sequence, known as the repeats' iScore rank. This illustrates a potential useful statistical ranking that, subject to biology may describe a repeats inherency to be more or less effective, in increments of the gene sequence. 

Repression of the most active sequences, especially in context of repeats may result in genetic variation. 

First Intron DNA - Site for a Genetic Brain?

DNA Methylation

The first intron of a gene, regardless of tissue or species is conserved as a site of downstream methylation with an inverse relationship to transcription and gene expression. Therefore, it is an informative gene feature regarding the relationship between DNA methylation and gene expression. But, expression in induced pluripotent stem cells (iPSC's) has been a major challenge to the stem cell industry, because by comparison these cells have not yet reached the state of natural pluripotent or embryonic stem cells (ESC's).

In mice two X chromosomes (XC) are active in the epiblasts of blastocysts as well as in pluripotent stem cells. One XC is inactivated triggered by Xist (non coding) RNA transcripts coating it to become silent. Designer transcription factor (dTF) repressors, binding the Xist intron 1 enhancer region caused higher H3K9me3 methylation and led to XC's opening and X-linked gene repression in MEFs. This substantially improved iPSC production and somatic cell nuclear transfer (SCNT) preimplantation embryonic development. This also correlated with much fewer abnormally expressed genes frequently associated with SCNT, even though it did not affect Xist expression. In stark contrast, the dTF activator targeting the same enhancer region drastically decreased both iPSC generation and SCNT efficiencies and induced ESC differentiation. 

A genome-wide, tissue-independent quasi-linear, inverse relationship exists between DNA methylation of the first intron and gene expression. More tissue-specific, differentially methylated regions exist in the first intron than in any other gene feature. These have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. CpGs in transcription factor binding motifs are enriched in the first intron and methylation tends to increase with distance from the first exon–first intron boundary, with a concomitant decrease in gene expression.

Since the relationship between sequence, methylation, repression and transcription is determinative in ESC differentiation it may also suggest a broader link to differential translation. Translation is required for miRNA-dependent transcript destabilization that alters levels of coding and noncoding transcripts. But, steady-state abundance and decay rates of cytosolic long non-coding RNA's (lncRNAs) are insensitive to miRNA loss. Instead lncRNAs fused to protein-coding reporter sequences become susceptible to miRNA-mediated decay. 

In this model, first intron DNA sequences that are differentially methylated, bind transcription factors that effect transcription, impact splicing, expressions of coding or non-coding transcripts and transcript destabilizations resulting in differential rates and possible variations in translation. This bottom-up, dynamic view of the classical process may elevate the first intron from 'junk' to a DNA 'brain' because it plays a more extensive role, heading the process toward translation of any gene or switching it off entirely.  

For this reason, among others Codondex uses first intron k-mers relative to the transcripts mRNA as the basis for comparing same gene transcripts in diseased cells or tissue samples. Further, p53 and BRCA1 miRNA key sequences, discovered using Codondex iScore algorithm, when transfected into HeLa cells resulted in significantly reduced proliferation that may result from this accelerated, transfected miRNA dependent decay.

 

Short Sequences of Proximally Disordered DNA

Oxford Nanopore Device Reducing Sequencing Cost

Relationships exist between short sequences of proximal DNA (SSPD) of a gene that when transcribed into RNA present stronger or weaker binding attractions to RNA binding proteins (RBP'S) that settle, edit, splice and resolve messenger RNA (mRNA). Responsive to epigenetic stimuli on Histones and DNA, mRNA are constantly transcribed in different quantity, at different times such that different mRNA strands are transported from the nucleus to cytoplasm where they are translated into and produce any of more than 30,000 different proteins.

Single nucleotide polymorphisms and DNA mutations can alter SSPD combinations in different diseased cells thus altering sequence proximity, ordering that affects transcribed RNA's attraction and optimal binding of RBP's. This may result in modified splicing of RNA, assembly of mRNA and slight or major variations in some or all translated protein derived from that gene. 

The specific effects of these DNA variations, on the multitude of proteins produced are generally unknown. However, it remains important to understand their effects in disease, diagnosis and therapy. Typically these have historically been researched by large scale analysis of RBP on RNA as opposed to the more fundamental, yet underrepresented massive array of diseased variant DNA to mRNA transitions.

Most pharmaceutical research is directed to a molecular interference targeting an aberrant protein to cure widely represented or highly impactful disease conditions of society. Economic assessments generally influence government decisions to support research based on loss of GDP contribution by a specific disease in a  patient cohort. However, in the modern multi-omics era top down research into protein-RNA activity is descending deeper into the cell to include RNA-mRNA and mRNA-DNA customizable therapies that will eventually resolve individually assessed diseases at a price that addresses much larger array of patient needs.  

SNP's and other mutations can vary considerably in cells. These variations can cause instability during division and lead to translated differences that can ultimately drive cancerous cell growth to escape patient immunity. Like a 'whack-a-mole' game, pattern variation and mechanistic persistence eventually beat the player. Without effective immune clearance these cells can replicate into tumors and contribute to microenvironments that support their existence.

Link to video on tumor microenvironment https://youtu.be/Z9H2utcnBic

We thought to analyze DNA and mRNA transcripts from cells in tumors and their microenvironments to see if we could expose the SSPD disordered combinations that may have promoted sub-optimal RBP attractions and led to sustained immune escape. Given the complexity of DNA to mRNA transcription, for any given gene many distortions in gene data sets have to be filtered. To do that we focused on p53, the most mutated gene in cancer. We designed a method to compare sequences arrays of DNA and mRNA Ensembl transcripts, from the consensus of healthy patients to multiple cell samples extracted from different sections of a patients tumor and tumor microenvironment.     

We previously identified and measured different levels of Natural Killer (NK) cell cytotoxicity, produced from cocultures with the extracted samples of each of the multiple sites of a biopsy. We will measure the different p53 transcript SSPD combinations associated with each sample and determine whether disordered SSPD's corelate with NK cytotoxicity from each coculture. We expect to identify whether biopsied tumor cells, ranked by SSPD's predict the cytotoxicity resulting from NK cell cocultures. We will narrow our research to identify the varied expressions of receptor combinations associated with degrees of cytotoxicity. We will test immune efficacy to lyse and destroy tumor cells. Finally we will test for adaptive immune response. 

Our vision is for per-patient, predictable cell co-culture pairings, for innate immune cell education based on ranking DNA-mRNA combinations to lead to multiple effective therapies. The falling cost of sequencing and sophistication of GMP laboratories presently servicing oncologists may support a successful use of this analytical approach to laboratory assisted disease management.

   

 

Non-Coding DNA Key Sequences

DNA Structural Inherency

Wind two strands of elastic, eventually it will knot, ultimately it will double up on itself. Separate the strands. From the point of unwinding, forces will be directed to different regions and the separation will approximately return to the wound state of the band. Do the same with each of 10 different bands or strings of any type, they will all behave in much the same way. For a given section of DNA being transcribed, the effect of separation will be much the same. For a given gene, there will be sequences that can tolerate force to greater or lesser degrees. For different transcripts, of a gene variation at those sequences may be crucial to the integrity of transcription machinery that separates DNA strands to initiate replication to RNA and for the outcome.

Cellular biology is enormously complex in all regards. The physics of molecular interaction, fluid dynamics, and chemistry combine in a system where cause and effect is near impossible to predict. At the most elementary level we hypothesize some non-coding DNA (ncDNA) possess structural inherencies that can be deployed to direct gene proteins and cell function for diagnosis or therapy.

Coding DNA and its regulatory, non-coding gene compliment is transcribed and spliced from a transcribed gene. Transcription to RNA, edited mRNA, spliced non-coding RNA and ultimately mRNA translation to protein can produce wide ranging, variable outcomes that may not be re-captured experimentally. 

A single nucleotide polymorphism (SNP) or SNP combinations within a gene may affect the finely tuned balance that results. Under different environmental conditions this could be material to the protein produced. Additionally other mutations of the gene could add complexity to the environment and/or the  resulting protein translation. 

At this level of cellular biology, genetic DNA stores instruction for protein assemblies to produce new protein required for the fully functional cell. However, DNA's stored mutations can lead to different functional or non-functional versions of protein depending on many different factors. Relationships between ncDNA, including mutations and the transcripts' edited, protein coding mRNA may represent unexplored inherencies that can regulate the gene's mRNA or translated protein.

We built an algorithm to elaborately compare ncDNA sequences of multiple protein coding transcripts of the same gene. For each transcript it steps through every variable length ncDNA sequence (kmer) (specifically intron1), computes a signature for each and indexes it to the constant of the transcripts' mRNA signature. For each step these signatures order the kmers for each of the transcript's. The order is represented in a vector of all the transcripts being compared.  

At millions of successive steps (depending on total intron 1 length's) transcripts mostly retain their vector ordering except, as expected at a kmer length change. Mostly transcript order in the vector does not change, occasionally a few positions change, vary rarely do all positions change. Position changes that cause another, like a domino effect are filtered out. For the rarest positions changes at a step, we look to the root causes in the kmer (sequence). We call this a Key Sequence because it is identified by the significance of changes to transcript positions in the vector compared to the vector at the next step. 

Therefore, Key Sequences cause the most position changes between transcripts being compared by the algorithm. This relative measure is step dependent and Key Sequences are discovered by comparing transcript positions in the vector at the next step location. Logically, this infers a genes structural inherency discovered through ncDNA Key Sequence relationships to mRNA, to other transcripts, error in gene alignments, sequenced reads or the algorithm. 

In assay testing we were able to predict and synthesize non-coding RNA Key Sequences that significantly reduced proliferation of HeLa cells. In our pre-clinical work, based on comparisons to transcripts of the TP53 we will be predicting the efficacy of cell and tissue selections that educate and activate Natural Killer cells.

If Key Sequences are inherent they could open a new frontier for diagnosis and therapy.