When Immunity Fails Programmed Cell Death

DNA Damage Response

Telomeric repeat (TR) sequences are responsible for genome integrity, where instability is a primary factor that leads to activation of p53. Introduction of a TR into cells leads to stabilization of p53, specific to TRs and not observed in plasmids containing non-TR sequences. TR-activated p53 exhibited enhanced transcriptional activity and induced p53-dependent growth suppression, measured as a reduction in colony formation. Sub-telomeric p53 binding prevents accumulation of DNA damage at human telomeres.  

Healthy cells experience thousands of DNA lesions per day. Micronuclei, containing broken fragments of DNA or chromosomes, that have become isolated, are recognized as one mediator of DNA damage response (DDR)-associated immune recognition. Like micronuclear DNA, mitochondrial DNA (mtDNA) is recognized by cGAS to drive STING-mediated inflammatory signaling. Mitochondrial damage can intersect DNA repair and inflammatory cascades with programmed cell death, through p53. In human fibroblasts and conditionally immortalized vascular smooth muscle cells p53 mediates CD54 (ICAM-1) overexpression in senescence.

Replicative senescence, an autophagy dependent program and crisis are anti-proliferative barriers that human cells must evade to gain immortality. Telomere-to-mitochondria signaling by ZBP1 mediates replicative crisis. Dysfunctional telomeres activate innate immune responses (IFN) through mitochondrial TR RNA (TERRA)–ZBP1 complexes. Senescence occurs when shortened telomeres elicit a p53 and RB dependent DNA-damage response. A crisis-associated isoform of ZBP1(innate immune sensor) is induced by the cGAS–STING DNA-sensing pathway, but reaches full activation only when associated with TERRA transcripts from dysfunctional telomeres. p53 utilizes the cGAS/STING innate immune system pathway for both cell intrinsic and cell extrinsic tumor suppressor activities. cGAS-STING activation induces the production of IFN-b and increases CD54 expression in  human cerebral microvascular endothelial cells.

In melanoma patients there is a significant correlation between cGAS expression levels and survival and between NK cell receptor expression levels and survival. Loss of cGAS expression by tumor cells could permit the tumor cell to circumvent senescence or prevent immunostimulatory NKG2D ligands expression. Loss of p53 and gain of oncogenic RAS exacerbated pro-malignant paracrine signaling activities of senescence-associated secretory phenotypes. Results imply that heterogeneity in cGAS activity, across tumors, could be an important predictor of cancer prognosis and response to treatment and suggest that NK cells could play an important role in mediating anti-tumor effects. Coculture of wild-type p53-induced human tumor cells with primary human NK cells enhanced NKG2D-dependent degranulation and IFN-γ production by NK cells. 

When p53 consensus sequences are modified and DNA damage response is compromised, replicative crisis ensues, mitochondrial membranes misfunction, mtDNA expression is downregulated and IFN signaling upregulates. A cell may then express activating immune ligands that bind NK receptors signaling non-self and cytolytic death or inhibitory receptors that signal self and immortality. 

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. 

Expanding Treatment Horizons

An unrecognized link between p53 function and the immunosurveillance of cancer and infection led to an understanding how p53 influences the expression of MHC molecules at the cell surface via binding interaction with endoplasmic reticulum ERAP1.

Targeted mutations in multiple cancers revealed TP53 gene expression ranged between the 89th and 100th percentile of all expressed transcripts, and raised the possibility that p53 peptides arising from these common mutations might be immunogenic in these patients.

Select KIR-HLA composition favoring antitumor activity could be a promising immunotherapeutic strategy against breast cancer using autologous activated Natural Killer (NK) cell clones. Coexistence of inhibitory and activating killer-cell immunoglobulin-like receptors (KIR) to the same cognate HLA-C2 and HLA-Bw4 ligands conferred breast cancer risk. Inhibitory KIR(iKIR)-HLA pairs without their activating KIR (aKIR)-HLA counterparts were significantly higher in normal controls. Contrarily and adding complexity this suggests NK cells expressing iKIR, to cognate HLA-ligands in the absence of specific aKIR counterparts are instrumental in antitumor response. 

Identification and characterization of the peptides presented by HLA-C, G and E molecules has been lacking behind the more abundant HLA-A and HLA-B gene products. The peptide specificities of these HLA molecules were elucidated using a comprehensive analysis of naturally presented peptides. The 15 most frequently expressed HLA-C alleles as well as HLA-E*01:01 and HLA-G*01:01 were transfected into lymphoblastoid C1R B-cells expressing low endogenous HLA. 

The results (above) include allotype C*02:02 for p53 presentation and indicate the overlap of HLA source protein and top 500 peptides demonstrating the enormous complexity for multivariate analysis of immune response. However,  C*02:02 and C*05:01 have identical contact residues for p8 and p9, the residues of the bound peptide that influences HLA-C interaction with KIR. This suggests peptide effects could contribute to the broader and stronger binding reactions of these two HLA-C allotypes. Interestingly SART3 and MAGEA3 proteins both interact through the p53 pathway and are reported in the peptide study (above) in addition to TP53 to present ligands on C*02:02 and C*05:01. 

Moreover, in vitro  models demonstrated that p53 is required for upregulation of NK ligands. Further, there was a strong association between the KIR B haplotype and p53 alteration in Basal Cell Carcinoma (BCC), with a higher likelihood that KIR B carriers harbor abnormal p53 (p<0.004). Together the data suggests functional interactions between KIR and HLA modify risks of BCC and Squamous Cell Carcinoma and that KIR encoded by the B genes provide selective pressure for altered p53 in BCC tumors.

Notwithstanding the enormous complexity between iKIR, aKIR - HLA interactions, immunoterapy must address the highly specific characteristics of autologous precision and discover methods to sensitively educate NK cells so that minimally invasive treatments can be extended to patients who fall outside the patient cohort for strictly regulated treatments. 

Of course, its never that simple...

Blood Pressure, Immunity and p53 Checkpoint.

Background

A few chromosome 19 curiosities developed into a deep-dive after looking into the primordial immune complex, the origins of MHC Class I and antigen receptors as revealed by comparative genomics. And the plot thickened because repressors (of endogenous retroviruses) that gained their binding affinity to retrovirus sequences at the same time their targets invaded the human lineage are preferentially located on chromosome 19. Further, the deletion rate in Zinc Finger clusters (ZNF) located around 19p.12 and 19q13.42, particularly between 51,012,739 and 55,620,741 are about twofold higher than the background deletion rate. A lot going on at this very active location which motivated this article.

At 19q13.42 kallikrein related peptidase (KLK’s), leukocyte immunoglobulin-like receptors (LILR’s) including killer-cell immunoglobulin-like receptor (KIR’s) as well MYADM, an important blood pressure related gene may also provide some clues to immunity variables that originate from or are influenced by this volatile region.

The retrotransposon bombardment of 19q13.42 and double background deletion rate is a significant remnant. However, after evolutionary MHC changed chromosomes ZNF, and within its range the chromosome 19 miRNA cluster (C19MC - 53,671,968 and 54,264,387) were still subjected to the deleterious effect of transposons. Regardless, suppression mechanics have kept epigenetic, regulatory and transcription processes, across gene’s far and wide on the move at a relatively stable rates. For example, reverse-transcribed SARS-CoV-2 RNA can integrate into the genome of cultured human cells and can be expressed in patient-derived tissues, but the effects of suppression may be sufficient to illicit a more permanent natural defense. In any event insertions and DNA damage are closely related and associated with loss of p53 that results in centrosome amplification. 

As cells pass through epithelial to mesenchymal transition (EMT), DNA damage prevents the normal reduction of p53 levels diverting the transcriptional program toward mesoderm without induction of an apoptotic response. In contrast, TP53-deficient cells differentiate to endoderm with high efficiency after DNA damage, suggesting that p53 enforces a “differentiation checkpoint” in early endoderm differentiation that alters cell fate in response to DNA damage.

Reproduction, Blood Pressure and NK

In reproduction, some of the 59 known miRNAs from primate-specific C19MC are highly expressed in human placentas and in the serum of pregnant women. They are also packaged into extracellular vesicles of diverse sizes, including exosomes and endow non-trophoblast cells with resistance to a variety of viruses. At least miR-517a-3p (a C19MC from fetal placenta) was incorporated into maternal NK cells in the third trimester, and it was rapidly cleared after delivery. miRNA's regulate the migration of human trophoblasts and suppress EMT genes critical for maintaining the epithelial cytotrophoblasts stem cell phenotype. 

Maternal uterine or decidual Natural Killer cells (dNK) express AT1, AT2, ANP, proteins of Renin Angiotensin System (RAS) suggesting dNK have the potential to contribute to changes in blood pressure that occur between days 5 and 12 of pregnancy in mice. And, pressure related mechanical stretch on endothelial cells interconnects innate and adaptive immune response in hypertension.

Pressure variables in cells and tissues may result from infection, inflammation and membrane stretch, including inner mitochondrial membrane that affects electron transport chain, endoplasmic reticulum, antigen production, presentation and exosome bound p53 / miRNA release.  ANP colocalization to dNK’s suggests that dNK RAS, at day 12 infers a localized RAS related responsiveness. STAT3 in monocytes was activated by increased endothelial stretch and is involved in driving almost all of the pathways that control NK cytolytic activity as well as the reciprocal regulatory interactions between NK cells and other components of the immune system. The crosstalk between STAT3 and p53/RAS signaling controls cancer cell metastasis and cisplatin resistance via the Slug/MAPK/PI3K/AKT-mediated regulation of EMT and autophagy.

Educating NK Subsets 

Looking into some of the ~15 genes scattered among C19MC (~sixty miRNA's) between 53,671,968 and 54,264,387;

1. MYADM was one of two blood pressure signature genes (copper uptake protein the other) differentially expressed for systolic, diastolic blood pressure and hypertension. Of the ~35 identified genes, several more strongly related to immune cell functions including PRF1, GNLY, TAGAP, IL2RB, GZMB and CD97, NKG7, CLC that are located on chromosome 19. The endothelium maintains a barrier between blood and tissue that becomes more permeable during inflammation. MYADM controls endothelial barrier function through ezrin, radixin, and moesin dependent regulation of ICAM-1 expression an essential receptor for NK interaction.

2. PRPF31 is recruited to introns following the attachment of U4 and U6 (spliceosome) RNA’s. Experiments using PRPF31 determined p53 activation is a general consequence of interfering with the spliceosome. 

3. At 54,617,158 LILRB1 receptor is expressed on immune cells where it binds to MHC class I molecules on antigen-presenting cells and transduces a negative signal that inhibits stimulation of an immune response. LILRB1 has a polymorphic regulatory region that enhances transcription in NK Cells and recruits zinc finger protein YY1 that inhibits p53. It also educates expanded human NK cells and defines a unique antitumor NK cell subset with potent antibody-dependent cellular cytotoxicity.

Monocyte/macrophage immunoglobulin-like receptors (MIR) genes are closely linked to the KIR gene family and the gene for FcαR at 19q13.4. The linkage was discovered in 1997 when a mouse sequence related to MIR mapped to a region on chromosome 7 syntenic with human 19q13.4. In 2012 a cluster of genetic loci, from multiple mouse strains and across anatomical sites was found to jointly contribute to the development of both thymic and splenic invariant natural killer T-cell NKT-cell levels. The dominant cluster was on mouse chromosome 7 and included almost all the non-C19MC genes located within the human C19MC region:– MYADM, CACNG7, VSTM1, TARM1, PRKCC(G), TFPT, NDUFA3, CNOT3, LENG1, TSEN34, RPS9. 

Four of nineteen knockout genes, that enhanced NK cell function were on chromosome 19 including GSK3 that phosphorylates Mdm2 to regulate p53 abundance, which would contribute to NK enhancement. 

A study of MHC disassortative mating in humans found Israeli’s were more gene similar, but MHC dissimilar than Europeans who were gene dissimilar and MHC dissimilar . Now, a recent study in American Indians found remarkably low KIR and HLA diversity in Amerindians that revealed signatures of strong purifying selection shaping the centromeric KIR region. This narrows to the importance of LILR-KIR region on chromosome19 that codes for the strongest NK cell educator receptors.

p53 regulates exosomes and miRNA’s directly influence NK responsiveness including regulation of dNK during pregnancy. Exosomes regulated by p53 also transfer it and can suppress growth and proliferation of p53 negative cells. Further, miRNA’s, induced by p53 can directly target ULBP2 mRNA and reduce its cell-surface expression.

Disease highlights

 

rs78378222 polymorphism in the 3'-untranslated region of TP53 contributes to development of age-associated cataracts by modifying miRNA-125b-induced apoptosis of lens epithelial cells. miRNA-125b is a novel negative regulator of p53. Deleting PRPF31 activates the p53 pathway and triggers retinal progenitor cells apoptosis. The members of the miR-125 family (miR-125a on chromosome 19q13.4 and miR-125b on chromosome 21q21.1) reside in two distinct human miRNA clusters with the let-7 and miR-99 families and these miRNAs are thus likely co-transcribed.

  

More succinctly, NK cells are alerted to induction of p53 in cancer cells by upregulation of the NKG2D ligands ULBP1 and ULBP2. p53 also induces expression of miR-34a and miR-34c, which target ULBP2 mRNA for destabilization. Observations suggest two possibly contrasting roles for p53 in NKG2DL expression and requires more investigation into how the regulation is fine-tuned. Extending this model to human populations would suggest that p53 must be inactivated among those with a robust NK response (those with B haplotypes). 

Taken together, our data suggest functional interactions between KIR and HLA modify risks of basal cell carcinoma (BCC) and squamous cell carcinoma, and that KIR encoded by the B genes provide selective pressure for altered p53 in BCC tumors. 

Conclusion

The convergence of several important cellular mechanisms that point back to a 19q13.42 address may illustrate ancient and conserved elements that perpetuate and function as integrated biological units effecting blood pressure, reproduction and immunity. Many of these impart education to innate immunity.

Immunity keeping p53 in check!

In a 2012 study on the topology of the human and mouse m6A RNA methylomes, Gene Ontology (GO) analysis of differentially expressed genes (DEG's) indicated a noteworthy enrichment of the p53 signaling pathway: 22/23 genes had differentially expressed splice variants, of which 18 were methylated. Moreover, 15 other members of the signaling pathway, which were not significant DEG's, exhibited significant differential isoform expressions. For example, isoforms of MDM4, needed for p53 inactivation were downregulated. Similar pro-apoptotic effects were observed in other pathway genes including MDM2, FAS and BAX. Higher apoptosis rate in HaCaT-T cells resulted with knockdown of m6A subunit METTL3, which also reversed a significant decrease in p53 activity. Modulation of p53 signaling through splicing may be relevant to induction of apoptosis by silencing of METTL3. 

Then, in 2019 a study of arsenite-induced human keratinocyte transformation demonstrated that knockdown of METTL3 significantly decreased m6A level, restored p53 activation and inhibited cellular transformation phenotypes in the-transformed cells. Further, m6A downregulated the expression of the positive p53 regulator, PRDM2, through the YTHDF2-promoted decay of PRDM2 mRNAs. m6A also upregulated expression of negative p53 regulator, YY1 and MDM2 through YTHDF1-stimulated translation of YY1 and MDM2 mRNA. Taken together, the study revealed the novel role of m6A in mediating human keratinocyte transformation by suppressing p53 activation and sheds light on the mechanisms of arsenic carcinogenesis via RNA epigenetics.

Finally in 2021 a discovery that YTHDF2 is upregulated in NK cells upon activation by cytokines, tumors, and cytomegalovirus infection. Ythdf2 deficiency in NK cells impaired its anti-tumor and anti-viral activity in vivo. YTHDF2 maintains NK cell homeostasis and terminal maturation, correlating with modulating NK cell trafficking and regulating Eomes, respectively. It promotes NK cell effector function and is required for IL-15-mediated NK cell survival and proliferation by forming a STAT5-YTHDF2 positive feedback loop. Analysis showed significant enrichment in cell cycle, division, and division-related processes, including mitotic cytokinesis, chromosome segregation, spindle, nucleosome, midbody, and chromosome. This data supports roles of YTHDF2 in regulating NK proliferation, survival, and effector functions. Transcriptome-wide screening identified Tardbp (TDP-43) to be involved in cell proliferation or survival as a YTHDF2-binding target in NK cells.

Downregulation of METTL3, which in spinal cord contributes with YTHDF2 to modulate inflammatory pain may upregulate differentially expressed p53 network splice variants that oppose YTHDF2 induced downregulation of p53, via PRDM2 leading to apoptotic or diseased cells. In diseased environments cytokines may upregulate YTHDF2 in NK cells leading to downregulation of p53 and cytoskeletal transformation that may be sufficient, at an immune synapse to advance cytolysis.

p53 signals may inform selections of cells and tissue that prime NK cells for advanced, personalized immune therapy. 

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.

 

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.