G4-forming Sequences Research Articles

G-Quadruplex Structures as Epigenetic Regulatory Elements in Priming of Defense Genes upon Short-Term Trichoderma atroviride Inoculation in Maize.

Symbiosis establishment between Trichoderma atroviride and plant roots triggers the priming of defense responses, among other effects. Currently, there is no clear evidence regarding the molecular mechanisms that allow the plant to remain alert to future stimulus, either by pathogen attack or any other abiotic stress. Epigenetic modifications have emerged as a strategy to explain the increased defense response of plants in a priming state conferred by Trichoderma. Recently, various non-canonical structures of nucleic acids, especially G-quadruplex structures (G-quadruplexes or G4s), have been identified as potential targets during the establishment or maintenance of plant signals. In the present study, we developed a screening test for the identification of putative G4-forming sequences (PQSs) in previously identified Z. mays priming genes. Bioinformatic analysis revealed the presence of PQSs in the promoter region of five essential genes playing a critical role in priming in maize. Biophysical and spectroscopy studies showed the formation of G4s by these PQSs in vitro, and ChIP assays demonstrate their formation in vivo. Therefore, G4 formation could play a role as an epigenetic regulatory mechanism involved in the long-lasting primed state in maize plants.

MMG4: Recognition of G4-Forming Sequences Based on Markov Model.

G-quadruplexes (G4s) are special nucleic acid structures with various important biological functions. Existing tools and technologies for G4-forming sequences recognition are limited to time-consuming and costly methods such as circular dichroism and nuclear magnetic resonance. Developing a fast and accurate model for G4-forming sequences recognition has far-reaching significance. In this study, MMG4, a novel model to recognize G4-forming sequences based on Markov model (MM), was developed and the phenomenon of high recognition accuracy in the central region of the sequence and low accuracy in the two end regions was discovered. It was further found that the differences in base transfer probabilities, ratio distribution, and G4-motif structural content in different regions may be the causes of this phenomenon. The study also explored the impact of sequence length on recognition accuracy and found the optimal recognition interval to be [910-1049], with the highest recognition accuracy reaching 85.95%. By extracting sequence features, the study constructed three types of machine learning models: random forest (RF), support vector machine, and back-propagation neural network. It was found that recognition performance of MM was significantly better than that of the other three machine learning models, proving that the recognition method based on MM can effectively capture the correlation information between adjacent nucleotides of G4. By combining MM with the three machine learning models, the predictive performance of MMG4 improved. Among them, the RF model combined with MM has the best performance, achieving an area under the receiver operating characteristic curve value of 0.93 and an area under the precision-recall curve value of 0.9. Finally, the study validated the model robustness and generalization ability through independent testing dataset.

Complete analysis of G-quadruplex forming sequences in the gapless assembly of human chromosome Y

Recent advancements have finally delivered a complete human genome assembly, including the elusive Y chromosome. This accomplishment closes a significant knowledge gap. Prior efforts were hampered by challenges in sequencing repetitive DNA structures such as direct and inverted repeats. We used the G4Hunter algorithm to analyze the presence of G-quadruplex forming sequences (G4s) within the current human reference genome (GRCh38) and the new telomere-to-telomere (T2T) Y chromosome assemblies. This analysis served a dual purpose: identifying the location of potential G4s within the genomes and exploring their association with functionally annotated sequences. Compared to GRCh38, the T2T assembly exhibited a significantly higher prevalence of G-quadruplex forming sequences. Notably, these repeats were abundantly located around precursor RNA, exons, genes, and within protein binding sites. This remarkable co-occurrence of G4-forming sequences with these critical regulatory regions suggests their role in fundamental DNA regulation processes. Our findings indicate that the current human reference genome significantly underestimated the number of G4s, potentially overlooking their functional importance.

ZNF91 is an endogenous repressor of the molecular phenotype associated with X-linked dystonia–parkinsonism (XDP)

X-linked dystonia-parkinsonism (XDP) is a severe neurodegenerative disorder resulting from an inherited intronic SINE-Alu-VNTR (SVA) retrotransposon in the TAF1 gene that causes dysregulation of TAF1 transcription. The specific mechanism underlying this dysregulation remains unclear, but it is hypothesized to involve the formation of G-quadruplexes (G4) structures within the XDP-SVA that impede transcription. In this study, we show that ZNF91, a critical repressor of SVA retrotransposons, specifically binds to G4-forming DNA sequences. Further, we found that genetic deletion of ZNF91 exacerbates the molecular phenotype associated with the XDP-SVA insertion in patient cells, while no difference was observed when ZNF91 was deleted from isogenic control cells. Additionally, we observed a significant age-related reduction in ZNF91 expression in whole blood and brain, indicating a progressive loss of repression of the XDP-SVA in XDP. These findings indicate that ZNF91 plays a crucial role in controlling the molecular phenotype associated with XDP. Since ZNF91 binds to G4-forming DNA sequences in SVAs, this suggests that interactions between ZNF91 and G4-forming sequences in the XDP-SVA minimize the severity of the molecular phenotype. Our results showing that ZNF91 expression levels significantly decrease with age provide a potential explanation for the age-related progressive neurodegenerative character of XDP. Collectively, our study provides important insights into the protective role of ZNF91 in XDP pathogenesis and suggests that restoring ZNF91 expression, destabilization of G4s, or targeted repression of the XDP-SVA could be future therapeutic strategies to prevent or treat XDP.

An Upstream G-Quadruplex DNA Structure Can Stimulate Gene Transcription.

Four-stranded G-quadruplexes (G4s) are DNA secondary structures that can form in the human genome. G4 structures have been detected in gene promoters and are associated with transcriptionally active chromatin and the recruitment of transcription factors and chromatin remodelers. We adopted a controlled, synthetic biology approach to understand how G4s can influence transcription. We stably integrated G4-forming sequences into the promoter of a synthetic reporter gene and inserted these into the genome of human cells. The integrated G4 sequences were shown to fold into a G4 structure within a cellular genomic context. We demonstrate that G4 structure formation within a gene promoter stimulates transcription compared to the corresponding G4-negative control promoter in a way that is not dependent on primary sequence or inherent G-richness. Systematic variation in the stability of folded G4s showed that in this system, transcriptional levels increased with higher stability of the G4 structure. By creating and manipulating a chromosomally integrated synthetic promoter, we have shown that G4 structure formation in a defined gene promoter can cause gene transcription to increase, which aligns with earlier observational correlations reported in the literature linking G4s to active transcription.

Genetic variations in G-quadruplex forming sequences affect the transcription of human disease-related genes.

Guanine-rich DNA strands can fold into non-canonical four-stranded secondary structures named G-quadruplexes (G4s). G4s folded in proximal promoter regions (PPR) are associated either with positive or negative transcriptional regulation. Given that single nucleotide variants (SNVs) affecting G4 folding (G4-Vars) may alter gene transcription, and that SNVs are associated with the human diseases' onset, we undertook a novel comprehensive study of the G4-Vars genome-wide (G4-variome) to find disease-associated G4-Vars located into PPRs. We developed a bioinformatics strategy to find disease-related SNVs located into PPRs simultaneously overlapping with putative G4-forming sequences (PQSs). We studied five G4-Vars disturbing in vitro the folding and stability of the G4s located into PPRs, which had been formerly associated with sporadic Alzheimer's disease (GRIN2B), a severe familiar coagulopathy (F7), atopic dermatitis (CSF2), myocardial infarction (SIRT1)and deafness (LHFPL5). Results obtained in cultured cells for these five G4-Vars suggest that the changes in the G4s affect the transcription, potentially contributing to the development of the mentioned diseases. Collectively, data reinforce the general idea that G4-Vars may impact on the different susceptibilities to human genetic diseases' onset, and could be novel targets for diagnosis and drug design in precision medicine.

G-quadruplex forming sequences in the genes coding for cytochrome P450 enzymes and their potential roles in drug metabolism

The majority of drugs are metabolized by cytochrome P450 (CYP) enzymes, primarily belonging to the CYP1, CYP2 and CYP3 families. Genetic variations are the main cause of inter-individual differences in drug response, which constitutes a major concern in pharmacotherapy. G-quadruplexes (G4s), are non-canonical DNA and RNA secondary structures formed by guanine-rich sequences. G4s have been implicated in cancer and gene regulation. In this study, we investigated putative G4-forming sequences (PQSs) in the CYP genes. Our findings reveal a high density of PQSs in the full genes of CYP family 2. Moreover, we observe an increased density of PQSs in the promoters of CYP family 1 genes compared to non-CYP450 genes. Importantly, stable PQSs were also identified in all studied CYP genes. Subsequently, we assessed the impact of the most frequently reported genetic mutations in the selected genes and the possible effect of these mutations on G4 formation as well as on the thermodynamic stability of predicted G4s. We found that 4 SNPs overlap G4 sequences and lead to mutated DNA and RNA G4 forming sequences in their context. Notably, the mutation in the CYP2C9 gene, which is associated with impaired (S)-warfarin metabolism in patients, alters a G4 sequence. We then demonstrated that at least 10 of the 13 chosen cytochrome P450 G4 candidates form G-quadruplex structures in vitro, using a combination of spectroscopic methods. In conclusion, our findings indicate the potential role of G-quadruplexes in the regulation of cytochrome genes, and emphasize the importance of G-quadruplexes in drug metabolism.

Genomic Instability of G-Quadruplex Sequences in Escherichia coli: Roles of DinG, RecG, and RecQ Helicases.

Guanine-rich DNA can fold into highly stable four-stranded DNA structures called G-quadruplexes (G4). Originally identified in sequences from telomeres and oncogene promoters, they can alter DNA metabolism. Indeed, G4-forming sequences represent obstacles for the DNA polymerase, with important consequences for cell life as they may lead to genomic instability. To understand their role in bacterial genomic instability, different G-quadruplex-forming repeats were cloned into an Escherichia coli genetic system that reports frameshifts and complete or partial deletions of the repeat when the G-tract comprises either the leading or lagging template strand during replication. These repeats formed stable G-quadruplexes in single-stranded DNA but not naturally supercoiled double-stranded DNA. Nevertheless, transcription promoted G-quadruplex formation in the resulting R-loop for (G3T)4 and (G3T)8 repeats. Depending on genetic background and sequence propensity for structure formation, mutation rates varied by five orders of magnitude. Furthermore, while in vitro approaches have shown that bacterial helicases can resolve G4, it is still unclear whether G4 unwinding is important in vivo. Here, we show that a mutation in recG decreased mutation rates, while deficiencies in the structure-specific helicases DinG and RecQ increased mutation rates. These results suggest that G-quadruplex formation promotes genetic instability in bacteria and that helicases play an important role in controlling this process in vivo.

In Silico Identification of Potential Quadruplex Forming Sequences in LncRNAs of Cervical Cancer.

Long non-coding RNAs (lncRNAs) have emerged as auxiliary regulators of gene expression influencing tumor microenvironment, metastasis and radio-resistance in cancer. The presence of lncRNA in extracellular fluids makes them promising diagnostic markers. LncRNAs deploy higher-order structures to facilitate a complex range of functions. Among such structures, G-quadruplexes (G4s) can be detected or targeted by small molecular probes to drive theranostic applications. The in vitro identification of G4 formation in lncRNAs can be a tedious and expensive proposition. Bioinformatics-driven strategies can provide comprehensive and economic alternatives in conjunction with suitable experimental validation. We propose a pipeline to identify G4-forming sequences, protein partners and biological functions associated with dysregulated lncRNAs in cervical cancer. We identified 17 lncRNA clusters which possess transcripts that can fold into a G4 structure. We confirmed in vitro G4 formation in the four biologically active isoforms of SNHG20, MEG3, CRNDE and LINP1 by Circular Dichroism spectroscopy and Thioflavin-T-assisted fluorescence spectroscopy and reverse-transcriptase stop assay. Gene expression data demonstrated that these four lncRNAs can be potential prognostic biomarkers of cervical cancer. Two approaches were employed for identifying G4 specific protein partners for these lncRNAs and FMR2 was a potential interacting partner for all four clusters. We report a detailed investigation of G4 formation in lncRNAs that are dysregulated in cervical cancer. LncRNAs MEG3, CRNDE, LINP1 and SNHG20 are shown to influence cervical cancer progression and we report G4 specific protein partners for these lncRNAs. The protein partners and G4s predicted in lncRNAs can be exploited for theranostic objectives.

Stabilization of RNA G-quadruplexes in the SARS-CoV-2 genome inhibits viral infection via translational suppression.

The G-quadruplex (G4) formed in single-stranded DNAs or RNAs plays a key role in diverse biological processes and is considered as a potential antiviral target. In the genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), 25 putative G4-forming sequences are predicted; however, the effects of G4-binding ligands on SARS-CoV-2 replication have not been studied in the context of viral infection. In this study, we investigated whether G4-ligands suppressed SARS-CoV-2 replication and whether their antiviral activity involved stabilization of viral RNA G4s and suppression of viral gene expression. We found that pyridostatin (PDS) suppressed viral gene expression and genome replication as effectively as the RNA polymerase inhibitor remdesivir. Biophysical analyses revealed that the 25 predicted G4s in the SARS-CoV-2 genome formed a parallel G4 structure. In particular, G4-644 and G4-3467 located in the 5' region of ORF1a, formed a G4 structure that could be effectively stabilized by PDS. We also showed that PDS significantly suppressed translation of the reporter genes containing these G4s. Taken together, our results demonstrate that stabilization of RNA G4s by PDS in the SARS-CoV-2 genome inhibits viral infection via translational suppression, highlighting the therapeutic potential of G4-ligands in SARS-CoV-2 infection.

Conserved G-Quadruplex-Forming Sequences in Mammalian TERT Promoters and Their Effect on Mutation Frequency.

Somatic mutations in the promoter region of the human telomerase reverse transcriptase (hTERT) gene have been identified in many types of cancer. The hTERT promoter is known to be enriched with sequences that enable the formation of G-quadruplex (G4) structures, whose presence is associated with elevated mutagenicity and genome instability. Here, we used a bioinformatics tool (QGRS mapper) to search for G4-forming sequences (G4 motifs) in the 1000 bp TERT promoter regions of 141 mammalian species belonging to 20 orders, 5 of which, including primates and predators, contain more than 10 species. Groups of conserved G4 motifs and single-nucleotide variants within these groups were discovered using a block alignment approach (based on the Nucleotide PanGenome explorer). It has been shown that: (i) G4 motifs are predominantly located in the region proximal to the transcription start site (up to 400 bp) and are over-represented on the non-coding strand of the TERT promoters, (ii) 11 to 22% of the G4 motifs found are evolutionarily conserved across the related organisms, and (iii) a statistically significant higher frequency of nucleotide substitutions in the conserved G4 motifs compared to the surrounding regions was confirmed only for the order Primates. These data support the assumption that G4s can interfere with the DNA repair process and affect the evolutionary adaptation of organisms and species.

Presence, Location and Conservation of Putative G-Quadruplex Forming Sequences in Arboviruses Infecting Humans.

Guanine quadruplexes (G4s) are non-canonical nucleic acid structures formed by guanine (G)-rich tracts that assemble into a core of stacked planar tetrads. G4s are found in the human genome and in the genomes of human pathogens, where they are involved in the regulation of gene expression and genome replication. G4s have been proposed as novel pharmacological targets in humans and their exploitation for antiviral therapy is an emerging research topic. Here, we report on the presence, conservation and localization of putative G4-forming sequences (PQSs) in human arboviruses. The prediction of PQSs was performed on more than twelve thousand viral genomes, belonging to forty different arboviruses that infect humans, and revealed that the abundance of PQSs in arboviruses is not related to the genomic GC content, but depends on the type of nucleic acid that constitutes the viral genome. Positive-strand ssRNA arboviruses, especially Flaviviruses, are significantly enriched in highly conserved PQSs, located in coding sequences (CDSs) or untranslated regions (UTRs). In contrast, negative-strand ssRNA and dsRNA arboviruses contain few conserved PQSs. Our analyses also revealed the presence of bulged PQSs, accounting for 17-26% of the total predicted PQSs. The data presented highlight the presence of highly conserved PQS in human arboviruses and present non-canonical nucleic acid-structures as promising therapeutic targets in arbovirus infections.

Phase Separation Modulates the Formation and Stabilities of DNA Guanine Quadruplex.

In the presence of monovalent alkali metal ions, G-rich DNA sequences containing four runs of contiguous guanines can fold into G-quadruplex (G4) structures. Recent studies showed that these structures are located in critical regions of the human genome and assume important functions in many essential DNA metabolic processes, including replication, transcription, and repair. However, not all potential G4-forming sequences are actually folded into G4 structures in cells, where G4 structures are known to be dynamic and modulated by G4-binding proteins as well as helicases. It remains unclear whether there are other factors influencing the formation and stability of G4 structures in cells. Herein, we showed that DNA G4s can undergo phase separation in vitro. In addition, immunofluorescence microscopy and ChIP-seq experiments with the use of BG4, a G4 structure-specific antibody, revealed that disruption of phase separation could result in global destabilization of G4 structures in cells. Together, our work revealed phase separation as a new determinant in modulating the formation and stability of G4 structures in human cells.

A new G-quadruplex-specific photosensitizer inducing genome instability in cancer cells by triggering oxidative DNA damage and impeding replication fork progression.

Photodynamic therapy (PDT) ideally relies on the administration, selective accumulation and photoactivation of a photosensitizer (PS) into diseased tissues. In this context, we report a new heavy-atom-free fluorescent G-quadruplex (G4) DNA-binding PS, named DBI. We reveal by fluorescence microscopy that DBI preferentially localizes in intraluminal vesicles (ILVs), precursors of exosomes, which are key components of cancer cell proliferation. Moreover, purified exosomal DNA was recognized by a G4-specific antibody, thus highlighting the presence of such G4-forming sequences in the vesicles. Despite the absence of fluorescence signal from DBI in nuclei, light-irradiated DBI-treated cells generated reactive oxygen species (ROS), triggering a 3-fold increase of nuclear G4 foci, slowing fork progression and elevated levels of both DNA base damage, 8-oxoguanine, and double-stranded DNA breaks. Consequently, DBI was found to exert significant phototoxic effects (at nanomolar scale) toward cancer cell lines and tumor organoids. Furthermore, in vivo testing reveals that photoactivation of DBI induces not only G4 formation and DNA damage but also apoptosis in zebrafish, specifically in the area where DBI had accumulated. Collectively, this approach shows significant promise for image-guided PDT.

Predicting the Three-Dimensional Structure of the c-KIT Proto-Oncogene Promoter and the Dynamics of Its Strongly Coupled Guanine Quadruplexes.

Guanine quadruplexes (G4s) play essential protective and regulatory roles within cells, influencing gene expression. In several gene-promoter regions, multiple G4-forming sequences are in close proximity and may form three-dimensional arrangements. We analyze the interplay among the three neighboring G4s in the c-KIT proto-oncogene promoter (WK1, WSP, and WK2). We highlight that the three G4s are structurally linked and their cross-talk favors the formation of a parallel structure for WSP. Relying on all-atom molecular dynamic simulations exceeding the μs time scale and using enhanced sampling methods, we provide the first computationally resolved structure of a well-organized G4 cluster in the promoter of a crucial gene involved in cancer development. Our results indicate that neighboring G4s influence their mutual three-dimensional arrangement and provide a powerful tool to predict and interpret complex DNA structures that can ultimately be used as a starting point for drug discovery.

Stable bulged G-quadruplexes in the human genome: identification, experimental validation and functionalization.

DNA sequence composition determines the topology and stability of G-quadruplexes (G4s). Bulged G-quadruplex structures (G4-Bs) are a subset of G4s characterized by 3D conformations with bulges. Current search algorithms fail to capture stable G4-B, making their genome-wide study infeasible. Here, we introduced a large family of computationally defined and experimentally verified potential G4-B forming sequences (pG4-BS). We found 478263 pG4-BS regions that do not overlap 'canonical' G4-forming sequences in the human genome and are preferentially localized in transcription regulatory regions including R-loops and open chromatin. Over 90% of protein-coding genes contain pG4-BS in their promoter or gene body. We observed generally higher pG4-BS content inR-loops and their flanks,longer genes that are associated with brain tissue, immune and developmental processes. Also, the presence of pG4-BS on both template and non-template strands in promoters is associated with oncogenesis, cardiovascular disease and stemness. Our G4-BS models predicted G4-forming ability in vitro with 91.5% accuracy. Analysis of G4-seq and CUT&Tag data strongly supports the existence of G4-BS conformations genome-wide. We reconstructed a novel G4-B 3D structure located in the E2F8 promoter. This study defines a large family of G4-like sequences, offering new insights into the essential biological functions and potential future therapeutic uses of G4-B.

Asymmetric triplex metallohelices stabilise DNA G-quadruplexes in promoter oncogene sequences and efficiently reduce their expression in cancer cells

Some metallo-supramolecular helical assemblies with size, shape, charge and amphipathic architectures similar to short cationic α-helical peptides have been shown to target and stabilise DNA G-quadruplexes (G4s) in vitro and downregulate the expression of G4-regulated genes in human cells. To expand the library of metallohelical structures that can act as efficient DNA G4 binders and downregulate genes containing G4-forming sequences in their promoter regions, we investigated the interaction of the two enantiomeric pairs of asymmetric Fe(II) triplex metallohelices with a series of five different DNA G4s formed by the human telomeric sequence (hTelo) and in the promoter regions of c-MYC, c-KIT, and k-RAS oncogenes. The metallohelices display preferential binding to G4s over duplex DNA in all investigated G4-forming sequences and induced arrest of DNA polymerase on template strands containing G4-forming sequences. Moreover, the investigated metallohelices suppressed the expression of c-MYC and k-RAS genes at mRNA and protein levels in HCT116 human cancer cells, as revealed by RT-qPCR analysis and western blotting.

Characterization of RNA G-quadruplexes in porcine epidemic diarrhea virus genome and the antiviral activity of G-quadruplex ligands

Porcine epidemic diarrhea virus (PEDV), an enteropathogenic coronavirus, has catastrophic impacts on the global pig industry. However, there are still no anti-PEDV drugs with accurate targets. G-quadruplexes (G4s) are non-canonical secondary structures formed within guanine-rich regions of DNA or RNA, and have attracted great attention as potential targets for antiviral strategy. In this study, we reported two putative G4-forming sequences (PQS) in S and Nsp5 genes of PEDV genome based on bioinformatic analysis, and identified that S-PQS and Nsp5-PQS were enabled to fold into G4 structure by using circular dichroism spectroscopy and fluorescence turn-on assay. Furthermore, we verified that both S-PQS and Nsp5-PQS PQS could form G4 structure in live cells by immunofluorescence microscopy. In addition, G4-specific compounds, such as TMPyP4 and PDS, could significantly inhibit transcription, translation and proliferation of PEDV in vitro. Importantly, these compounds exert antiviral activity at the post-entry step of PEDV infection cycle, by inhibiting viral genome replication and protein expression. Lastly, we demonstrated that TMPyP4 can inhibit reporter gene expression by targeting G4 structure in Nsp5. Taken together, these findings not only reinforce the presence of viral G-quadruplex sequences in PEDV genome but also provide new insights into developing novel antiviral drugs targeting PEDV RNA G-quadruplexes.

G-quadruplexes formed by Varicella-Zoster virus reiteration sequences suppress expression of glycoprotein C and regulate viral cell-to-cell spread.

G-quadruplex (G4) formed by repetitive guanosine-rich sequences plays important roles in diverse cellular processes; however, its roles in viral infection are not fully understood. In this study, we investigated the genome-wide distribution of G4-forming sequences (G4 motifs) in Varicella-Zoster virus (VZV) and found that G4 motifs are enriched in the internal repeat short and the terminal repeat short regions flanking the unique short region and also in some reiteration (R) sequence regions. A high density of G4 motifs in the R2 region was found on the template strand of ORF14, which encodes glycoprotein C (gC), a virulent factor for viral growth in skin. Analyses such as circular dichroism spectroscopy, thermal difference spectra, and native polyacrylamide gel electrophoresis with oligodeoxynucleotides demonstrated that several G4 motifs in ORF14 form stable G4 structures. In transfection assays, gC expression from the G4-disrupted ORF14 gene was increased at the transcriptional level and became more resistant to suppression by G4-ligand treatment. The recombinant virus containing the G4-disrupted ORF14 gene expressed a higher level of gC mRNA, while it showed a slightly reduced growth. This G4-disrupted ORF14 virus produced smaller plaques than the wild-type virus. Our results demonstrate that G4 formation via reiteration sequences suppresses gC expression during VZV infection and regulates viral cell-to-cell spread.

Genomic Analysis of Non-B Nucleic Acids Structures in SARS-CoV-2: Potential Key Roles for These Structures in Mutability, Translation, and Replication?

Non-B nucleic acids structures have arisen as key contributors to genetic variation in SARS-CoV-2. Herein, we investigated the presence of defining spike protein mutations falling within inverted repeats (IRs) for 18 SARS-CoV-2 variants, discussed the potential roles of G-quadruplexes (G4s) in SARS-CoV-2 biology, and identified potential pseudoknots within the SARS-CoV-2 genome. Surprisingly, there was a large variation in the number of defining spike protein mutations arising within IRs between variants and these were more likely to occur in the stem region of the predicted hairpin stem-loop secondary structure. Notably, mutations implicated in ACE2 binding and propagation (e.g., ΔH69/V70, N501Y, and D614G) were likely to occur within IRs, whilst mutations involved in antibody neutralization and reduced vaccine efficacy (e.g., T19R, ΔE156, ΔF157, R158G, and G446S) were rarely found within IRs. We also predicted that RNA pseudoknots could predominantly be found within, or next to, 29 mutations found in the SARS-CoV-2 spike protein. Finally, the Omicron variants BA.2, BA.4, BA.5, BA.2.12.1, and BA.2.75 appear to have lost two of the predicted G4-forming sequences found in other variants. These were found in nsp2 and the sequence complementary to the conserved stem-loop II-like motif (S2M) in the 3' untranslated region (UTR). Taken together, non-B nucleic acids structures likely play an integral role in SARS-CoV-2 evolution and genetic diversity.