Simple Repeats as Building Blocks for Genetic Computers

Flipons enable genomes to learn by intermediating the exchange of energy for information

Recent findings have confirmed the long-held belief that alternative DNA conformations encoded by genetic elements called flipons have important biological roles. Many of these alternative structures are formed by sequences originally spread throughout the human genome by endogenous retroelements (ERE) that captured 50% of the territory before being disarmed. Only 2.6% of the remaining DNA codes for proteins. Other organisms have instead streamlined their genomes by eliminating invasive retroelements and other repeat elements. The question arises, why retain any ERE at all? A new synthesis suggests that flipons enable genomes to learn and programme the context-specific readout of information by altering the transcripts produced. The exchange of energy for information is mediated through changes in DNA topology. Here I provide a formulation for how genomes learn and describe the underlying p-bit algorithm through which flipons are tuned. The framework suggests new strategies for the therapeutic reprogramming of cells.

Structural insights into the recognition of purine-pyrimidine dinucleotide repeats by zinc finger protein ZBTB43

Purine-pyrimidine repeats (PPRs) can form left-handed Z-form DNA and induce DNA double-strand breaks (DSBs), posing a risk for genomic rearrangements and cancer. The zinc finger (ZF) and BTB domain-containing protein 43 (ZBTB43) is a transcription factor containing two Cys2-His2 (C2H2) and one C3H1 zinc fingers and plays a crucial role in maintaining genomic and epigenomic integrity by converting mutagenic Z-form PPRs to the B-form in prospermatogonia. Despite its importance, the molecular mechanism underlying the recognition of PPRs by ZBTB43 remains elusive. In this study, we determined the X-ray crystal structure of the ZBTB43 ZF1-3 in complex with the B-form DNA containing the CA repeats sequence. The structure reveals that ZF1 and ZF2 primarily recognize the CACA sequence through specific hydrogen-bonding and van der Waals contacts via a quadruple center involving Arg389, Met411, His413, and His414. These interactions were further validated by fluorescence-based DNA-binding assays using mutated ZBTB43 variants. Our structural investigation provides valuable insights into the recognition mechanism of PPRs by ZBTB43 and suggests a potential role for ZBTB43 in the transformation of Z-DNA to B-DNA, contributing to the maintenance of genomic stability.

Z-DNA formation in promoters conserved between human and mouse are associated with increased transcription reinitiation rates

A long-standing question concerns the role of Z-DNA in transcription. Here we use a deep learning approach DeepZ that predicts Z-flipons based on DNA sequence, structural properties of nucleotides and omics data. We examined Z-flipons that are conserved between human and mouse genomes after generating whole-genome Z-flipon maps and then validated them by orthogonal approaches based on high resolution chemical mapping of Z-DNA and the transformer algorithm Z-DNABERT. For human and mouse, we revealed similar pattern of transcription factors, chromatin remodelers, and histone marks associated with conserved Z-flipons. We found significant enrichment of Z-flipons in alternative and bidirectional promoters associated with neurogenesis genes. We show that conserved Z-flipons are associated with increased experimentally determined transcription reinitiation rates compared to promoters without Z-flipons, but without affecting elongation or pausing. Our findings support a model where Z-flipons engage Transcription Factor E and impact phenotype by enabling the reset of preinitiation complexes when active, and the suppression of gene expression when engaged by repressive chromatin complexes.

Expanding horizons of tandem repeats in biology and medicine: Why 'genomic dark matter' matters

Approximately half of the human genome includes repetitive sequences, and these DNA sequences (as well as their transcribed repetitive RNA and translated amino-acid repeat sequences) are known as the repeatome. Within this repeatome there are a couple of million tandem repeats, dispersed throughout the genome. These tandem repeats have been estimated to constitute ∼8% of the entire human genome. These tandem repeats can be located throughout exons, introns and intergenic regions, thus potentially affecting the structure and function of tandemly repetitive DNA, RNA and protein sequences. Over more than three decades, more than 60 monogenic human disorders have been found to be caused by tandem-repeat mutations. These monogenic tandem-repeat disorders include Huntington's disease, a variety of ataxias, amyotrophic lateral sclerosis and frontotemporal dementia, as well as many other neurodegenerative diseases. Furthermore, tandem-repeat disorders can include fragile X syndrome, related fragile X disorders, as well as other neurological and psychiatric disorders. However, these monogenic tandem-repeat disorders, which were discovered via their dominant or recessive modes of inheritance, may represent the 'tip of the iceberg' with respect to tandem-repeat contributions to human disorders. A previous proposal that tandem repeats may contribute to the 'missing heritability' of various common polygenic human disorders has recently been supported by a variety of new evidence. This includes genome-wide studies that associate tandem-repeat mutations with autism, schizophrenia, Parkinson's disease and various types of cancers. In this article, I will discuss how tandem-repeat mutations and polymorphisms could contribute to a wide range of common disorders, along with some of the many major challenges of tandem-repeat biology and medicine. Finally, I will discuss the potential of tandem repeats to be therapeutically targeted, so as to prevent and treat an expanding range of human disorders.

The Intransitive Logic of Directed Cycles and Flipons Enhances the Evolution of Molecular Computers by Augmenting the Kolmogorov Complexity of Genomes

Cell responses are usually viewed as transitive events with fixed inputs and outputs that are regulated by feedback loops. In contrast, directed cycles (DCs) have all nodes connected, and the flow is in a single direction. Consequently, DCs can regenerate themselves and implement intransitive logic. DCs are able to couple unrelated chemical reactions to each edge. The output depends upon which node is used as input. DCs can also undergo selection to minimize the loss of thermodynamic entropy while maximizing the gain of information entropy. The intransitive logic underlying DCs enhances their programmability and impacts their evolution. The natural selection of DCs favors the persistence, adaptability, and self-awareness of living organisms and does not depend solely on changes to coding sequences. Rather, the process can be RNA-directed. I use flipons, nucleic acid sequences that change conformation under physiological conditions, as a simple example and then describe more complex DCs. Flipons are often encoded by repeats and greatly increase the Kolmogorov complexity of genomes by adopting alternative structures. Other DCs allow cells to regenerate, recalibrate, reset, repair, and rewrite themselves, going far beyond the capabilities of current computational devices. Unlike Turing machines, cells are not designed to halt but rather to regenerate.

Flipons and small RNAs accentuate the asymmetries of pervasive transcription by the reset and sequence-specific microcoding of promoter conformation

The role of alternate DNA conformations such as Z-DNA in the regulation of transcription is currently underappreciated. These structures are encoded by sequences called flipons, many of which are enriched in promoter and enhancer regions. Through a change in their conformation, flipons provide a tunable mechanism to mechanically reset promoters for the next round of transcription. They act as actuators that capture and release energy to ensure that the turnover of the proteins at promoters is optimized to cell state. Likewise, the single-stranded DNA formed as flipons cycle facilitates the docking of RNAs that are able to microcode promoter conformations and canalize the pervasive transcription commonly observed in metazoan genomes. The strand-specific nature of the interaction between RNA and DNA likely accounts for the known asymmetry of epigenetic marks present on the histone tetramers that pair to form nucleosomes. The role of these supercoil-dependent processes in promoter choice and transcriptional interference is reviewed. The evolutionary implications are examined: the resilience and canalization of flipon-dependent gene regulation is contrasted with the rapid adaptation enabled by the spread of flipon repeats throughout the genome. Overall, the current findings underscore the important role of flipons in modulating the readout of genetic information and how little we know about their biology.

Z-flipon variants reveal the many roles of Z-DNA and Z-RNA in health and disease

Identifying roles for Z-DNA remains challenging given their dynamic nature. Here, we perform genome-wide interrogation with the DNABERT transformer algorithm trained on experimentally identified Z-DNA forming sequences (Z-flipons). The algorithm yields large performance enhancements (F1 = 0.83) over existing approaches and implements computational mutagenesis to assess the effects of base substitution on Z-DNA formation. We show Z-flipons are enriched in promoters and telomeres, overlapping quantitative trait loci for RNA expression, RNA editing, splicing, and disease-associated variants. We cross-validate across a number of orthogonal databases and define BZ junction motifs. Surprisingly, many effects we delineate are likely mediated through Z-RNA formation. A shared Z-RNA motif is identified in SCARF2, SMAD1, and CACNA1 transcripts, whereas other motifs are present in noncoding RNAs. We provide evidence for a Z-RNA fold that promotes adaptive immunity through alternative splicing of KRAB domain zinc finger proteins. An analysis of OMIM and presumptive gnomAD loss-of-function datasets reveals an overlap of Z-flipons with disease-causing variants in 8.6% and 2.9% of Mendelian disease genes, respectively, greatly extending the range of phenotypes mapped to Z-flipons.

HnRNP K reduces viral gene expression by targeting cytosine-rich sequences in porcine reproductive and respiratory syndrome virus-2 genome to dampen the viral growth

HnRNP K is a well-known member of HnRNP family proteins that has been implicated in the regulation of protein expression. Currently, the impact of HnRNP K on the reproduction cycle of a broad range of virus were reported, while the precise function for PRRSV was lacking. In this study, we determined that both PRRSV infection and ectopic expression of N protein induced an enrichment of HnRNP K in the cytoplasm. Using RNA pulldown and RNA immunoprecipitation, we described the interactions between the KH2 domain of HnRNP K and cytosine-rich sequences (CRS) in PRRSV genomic RNA corresponding to Nsp7α coding region. Meanwhile, overexpression of HnRNP K inhibited viral gene expression and PRRSV replication, while silencing of HnRNP K resulted in an increased in virus yield. Taken together, this study assists in the understanding of PRRSV-host interactions, and the development of vaccines based on viral genome engineering.

Global abundance of short tandem repeats is non-random in rodents and primates

Background

While of predominant abundance across vertebrate genomes and significant biological implications, the relevance of short tandem repeats (STRs) (also known as microsatellites) to speciation remains largely elusive and attributed to random coincidence for the most part. Here we collected data on the whole-genome abundance of mono-, di-, and trinucleotide STRs in nine species, encompassing rodents and primates, including rat, mouse, olive baboon, gelada, macaque, gorilla, chimpanzee, bonobo, and human. The collected data were used to analyze hierarchical clustering of the STR abundances in the selected species.

Results

We found massive differential STR abundances between the rodent and primate orders. In addition, while numerous STRs had random abundance across the nine selected species, the global abundance conformed to three consistent < clusters>, as follows: <rat, mouse>, <gelada, macaque, olive baboon>, and <gorilla, chimpanzee, bonobo, human>, which coincided with the phylogenetic distances of the selected species (p < 4E-05). Exceptionally, in the trinucleotide STR compartment, human was significantly distant from all other species.

Conclusion

Based on hierarchical clustering, we propose that the global abundance of STRs is non-random in rodents and primates, and probably had a determining impact on the speciation of the two orders. We also propose the STRs and STR lengths, which predominantly conformed to the phylogeny of the selected species, exemplified by (t)10, (ct)6, and (taa4). Phylogenetic and experimental platforms are warranted to further examine the observed patterns and the biological mechanisms associated with those STRs.

Mononucleotide A-repeats may Play a Regulatory Role in Endothermic Housekeeping Genes

Background:

Coding and non-coding short tandem repeats (STRs) facilitate a great diversity of phenotypic traits. The imbalance of mononucleotide A-repeats around transcription start sites (TSSs) was found in 3 mammals: H. sapiens, M. musculus, and R. norvegicus.

Principal Findings:

We found that the imbalance pattern originated in some vertebrates. A similar pattern was observed in mammals and birds, but not in amphibians and reptiles. We proposed that the enriched A-repeats upstream of TSSs is a novel hallmark of endotherms or warm-blooded animals. Gene ontology analysis indicates that the primary function of upstream A-repeats involves metabolism, cellular transportation, and sensory perception (smell and chemical stimulus) through housekeeping genes.

Conclusions:

Upstream A-repeats may play a regulatory role in the metabolic process of endothermic animals.

Interaction of HnRNP F with the guanine-rich segments in viral antigenomic RNA enhances porcine reproductive and respiratory syndrome virus-2 replication

Background

Heterogeneous nuclear ribonucleoprotein (HnRNP) F is a member of HnRNP family proteins that participate in splicing of cellular newly synthesized mRNAs by specifically recognizing tandem guanine-tracts (G-tracts) RNA sequences. Whether HnRNP F could recognize viral-derived tandem G-tracts and affect virus replication remain poorly defined.

Methods

The effect of HnRNP F on porcine reproductive and respiratory syndrome virus (PRRSV) propagation was evaluated by real-time PCR, western blotting, and plaque-forming unit assay. The association between HnRNP F and PRRSV guanine-rich segments (GRS) were analyzed by RNA pulldown and RNA immunoprecipitation. The expression pattern of HnRNP F was investigated by western blotting and nuclear and cytoplasmic fractionation.

Results

Knockdown of endogenous HnRNP F effectively blocks the synthesis of viral RNA and nucleocapsid (N) protein. Conversely, overexpression of porcine HnRNP F has the opposite effect. Moreover, RNA pulldown and RNA immunoprecipitation assays reveal that the qRMM1 and qRRM2 domains of HnRNP F recognize the GRS in PRRSV antigenomic RNA. Finally, HnRNP F is redistributed into the cytoplasm and forms a complex with guanine-quadruplex (G4) helicase DHX36 during PRRSV infection.

Conclusions

These findings elucidate the potential functions of HnRNP F in regulating the proliferation of PRRSV and contribute to a better molecular understanding of host-PRRSV interactions.

Mono a Mano: ZBP1's Love-Hate Relationship with the Kissing Virus

Z-DNA binding protein (ZBP1) very much represents the nuclear option. By initiating inflammatory cell death (ICD), ZBP1 activates host defenses to destroy infectious threats. ZBP1 is also able to induce noninflammatory regulated cell death via apoptosis (RCD). ZBP1 senses the presence of left-handed Z-DNA and Z-RNA (ZNA), including that formed by expression of endogenous retroelements. Viruses such as the Epstein-Barr "kissing virus" inhibit ICD, RCD and other cell death signaling pathways to produce persistent infection. EBV undergoes lytic replication in plasma cells, which maintain detectable levels of basal ZBP1 expression, leading us to suggest a new role for ZBP1 in maintaining EBV latency, one of benefit for both host and virus. We provide an overview of the pathways that are involved in establishing latent infection, including those regulated by MYC and NF-κB. We describe and provide a synthesis of the evidence supporting a role for ZNA in these pathways, highlighting the positive and negative selection of ZNA forming sequences in the EBV genome that underscores the coadaptation of host and virus. Instead of a fight to the death, a state of détente now exists where persistent infection by the virus is tolerated by the host, while disease outcomes such as death, autoimmunity and cancer are minimized. Based on these new insights, we propose actionable therapeutic approaches to unhost EBV.

Predominant monomorphism of the RIT2 and GPM6B exceptionally long GA blocks in human and enriched divergent alleles in the disease compartment

Across human protein-coding genes, the human neuron-specific genes, RIT2 and GPM6B, contain the two longest GA short tandem repeats (STRs) of 11 and 9-repeats, respectively, the length ranges of which are functional, and result in gene expression alteration. Here we sequenced the RIT2 and GPM6B STRs in 600 human subjects, consisting of late-onset neurocognitive disorder (n = 200), multiple sclerosis (n = 200), and controls (n = 200). Furthermore, we selected two large human databases, including the general-population-based gnomAD ( https://gnomad.broadinstitute.org ) and a mainly disease-phenotype-archiving database, TOPMed ( https://www.nhlbiwgs.org ), to compare allele frequencies in the general populations vs. the disease compartment. The RIT2 and GPM6B GA-repeats were monomorphic in the human subjects studied, at lengths of 11 and 9-repeats, respectively, and were predominantly human-specific in formula. Exception included a 9/11 genotype of the RIT2 GA-STR in an isolate case of female multiple sclerosis. Exceedingly rare alleles of the two GA repeats were significantly enriched in TOPMed vs. the gnomAD. We report prime instances of predominant monomorphism for specific lengths of STRs in human, and possible enrichment of rare divergent alleles in the disease phenotype compartment. While STRs are most attended because of their high polymorphic nature, STR monomorphism is an underappreciated feature, which may have a link with natural selection and disease.

The Simple Biology of Flipons and Condensates Enhances the Evolution of Complexity

The classical genetic code maps nucleotide triplets to amino acids. The associated sequence composition is complex, representing many elaborations during evolution of form and function. Other genomic elements code for the expression and processing of RNA transcripts. However, over 50% of the human genome consists of widely dispersed repetitive sequences. Among these are simple sequence repeats (SSRs), representing a class of flipons, that under physiological conditions, form alternative nucleic acid conformations such as Z-DNA, G4 quartets, I-motifs, and triplexes. Proteins that bind in a structure-specific manner enable the seeding of condensates with the potential to regulate a wide range of biological processes. SSRs also encode the low complexity peptide repeats to patch condensates together, increasing the number of combinations possible. In situations where SSRs are transcribed, SSR-specific, single-stranded binding proteins may further impact condensate formation. Jointly, flipons and patches speed evolution by enhancing the functionality of condensates. Here, the focus is on the selection of SSR flipons and peptide patches that solve for survival under a wide range of environmental contexts, generating complexity with simple parts.

Heterogeneous nuclear ribonucleoprotein E1 binds polycytosine DNA and monitors genome integrity

hnRNP E1 binds polycytosine tracts of DNA and monitors genome integrity.

Heterogeneous nuclear ribonucleoprotein E1 (hnRNP E1) is a tumor suppressor protein that binds site- and structure-specifically to RNA sequences to regulate mRNA stability, facilitate alternative splicing, and suppress protein translation on several metastasis-associated mRNAs. Here, we show that hnRNP E1 binds polycytosine-rich DNA tracts present throughout the genome, including those at promoters of several oncogenes and telomeres and monitors genome integrity. It binds DNA in a site- and structure-specific manner. hnRNP E1-knockdown cells displayed increased DNA damage signals including γ-H2AX at its binding sites and also showed increased mutations. UV and hydroxyurea treatment of hnRNP E1-knockdown cells exacerbated the basal DNA damage signals with increased cell cycle arrest, activation of checkpoint proteins, and monoubiquitination of proliferating cell nuclear antigen despite no changes in deubiquitinating enzymes. DNA damage caused by genotoxin treatment localized to hnRNP E1 binding sites. Our work suggests that hnRNP E1 facilitates functions of DNA integrity proteins at polycytosine tracts and monitors DNA integrity at these sites.

To "Z" or not to "Z": Z-RNA, self-recognition, and the MDA5 helicase

Double-stranded RNA (dsRNA) is produced both by virus and host. Its recognition by the melanoma differentiation-associated gene 5 (MDA5) initiates type I interferon responses. How can a host distinguish self-transcripts from nonself to ensure that responses are targeted correctly? Here, I discuss a role for MDA5 helicase in inducing Z-RNA formation by Alu inverted repeat (AIR) elements. These retroelements have highly conserved sequences that favor Z-formation, creating a site for the dsRNA-specific deaminase enzyme ADAR1 to dock. The subsequent editing destabilizes the dsRNA, ending further interaction with MDA5 and terminating innate immune responses directed against self. By enabling self-recognition, Alu retrotransposons, once invaders, now are genetic elements that keep immune responses in check. I also discuss the possible but less characterized roles of the other helicases in modulating innate immune responses, focusing on DExH-box helicase 9 (DHX9) and Mov10 RISC complex RNA helicase (MOV10). DHX9 and MOV10 function differently from MDA5, but still use nucleic acid structure, rather than nucleotide sequence, to define self. Those genetic elements encoding the alternative conformations involved, referred to as flipons, enable helicases to dynamically shape a cell's repertoire of responses. In the case of MDA5, Alu flipons switch off the dsRNA-dependent responses against self. I suggest a number of genetic systems in which to study interactions between flipons and helicases further.