Interaction of Proteins with Inverted Repeats and Cruciform Structures in Nucleic Acids

A comprehensive study of Z-DNA density and its evolutionary implications in birds

BACKGROUND

Z-DNA, a left-handed helical form of DNA, plays a significant role in genomic stability and gene regulation. Its formation, associated with high GC content and repetitive sequences, is linked to genomic instability, potentially leading to large-scale deletions and contributing to phenotypic diversity and evolutionary adaptation.

RESULTS

In this study, we analyzed the density of Z-DNA-prone motifs of 154 avian genomes using the non-B DNA Motif Search Tool (nBMST). Our findings indicate a higher prevalence of Z-DNA motifs in promoter regions across all avian species compared to other genomic regions. A negative correlation was observed between Z-DNA density and developmental time in birds, suggesting that species with shorter developmental periods tend to have higher Z-DNA densities. This relationship implies that Z-DNA may influence the timing and regulation of development in avian species. Furthermore, Z-DNA density showed associations with traits such as body mass, egg mass, and genome size, highlighting the complex interactions between genome architecture and phenotypic characteristics. Gene Ontology (GO) analysis revealed that Z-DNA motifs are enriched in genes involved in nucleic acid binding, kinase activity, and translation regulation, suggesting a role in fine-tuning gene expression essential for cellular functions and responses to environmental changes. Additionally, the potential of Z-DNA to drive genomic instability and facilitate adaptive evolution underscores its importance in shaping phenotypic diversity.

CONCLUSIONS

This study emphasizes the role of Z-DNA as a dynamic genomic element contributing to gene regulation, genomic stability, and phenotypic diversity in avian species. Future research should experimentally validate these associations and explore the molecular mechanisms by which Z-DNA influences avian biology.

invertiaDB: A Database of Inverted Repeats Across Organismal Genomes

Inverted repeats are repetitive elements that can form hairpin and cruciform structures. They are linked to genomic instability, however they also have various biological functions. Their distribution differs markedly across taxonomic groups in the tree of life, and they exhibit high polymorphism due to their inherent genomic instability. Advances in sequencing technologies and declined costs have enabled the generation of an ever-growing number of complete genomes for organisms across taxonomic groups in the tree of life. However, a comprehensive database encompassing inverted repeats across diverse organismal genomes has been lacking. We present InvertiaDB, the first comprehensive database of inverted repeats spanning multiple taxa, featuring repeats identified in the genomes of 118,070 organisms across all major taxonomic groups. The database currently hosts 30,067,666 inverted repeat sequences, serving as a centralized, user-friendly repository to perform searches, interactive visualization, and download existing inverted repeat data for independent analysis. invertiaDB is implemented as a web portal for browsing, analyzing and downloading inverted repeat data. invertiaDB is publicly available at https://invertiadb.netlify.app/homepage.html.

A machine learning enhanced EMS mutagenesis probability map for efficient identification of causal mutations in Caenorhabditis elegans

Chemical mutagenesis-driven forward genetic screens are pivotal in unveiling gene functions, yet identifying causal mutations behind phenotypes remains laborious, hindering their high-throughput application. Here, we reveal a non-uniform mutation rate caused by Ethyl Methane Sulfonate (EMS) mutagenesis in the C. elegans genome, indicating that mutation frequency is influenced by proximate sequence context and chromatin status. Leveraging these factors, we developed a machine learning enhanced pipeline to create a comprehensive EMS mutagenesis probability map for the C. elegans genome. This map operates on the principle that causative mutations are enriched in genetic screens targeting specific phenotypes among random mutations. Applying this map to Whole Genome Sequencing (WGS) data of genetic suppressors that rescue a C. elegans ciliary kinesin mutant, we successfully pinpointed causal mutations without generating recombinant inbred lines. This method can be adapted in other species, offering a scalable approach for identifying causal genes and revitalizing the effectiveness of forward genetic screens.

Extraribosomal Functions of Bacterial Ribosomal Proteins-An Update, 2023

Ribosomal proteins (r-proteins) are abundant, highly conserved, and multifaceted cellular proteins in all domains of life. Most r-proteins have RNA-binding properties and can form protein-protein contacts. Bacterial r-proteins govern the co-transcriptional rRNA folding during ribosome assembly and participate in the formation of the ribosome functional sites, such as the mRNA-binding site, tRNA-binding sites, the peptidyl transferase center, and the protein exit tunnel. In addition to their primary role in a cell as integral components of the protein synthesis machinery, many r-proteins can function beyond the ribosome (the phenomenon known as moonlighting), acting either as individual regulatory proteins or in complexes with various cellular components. The extraribosomal activities of r-proteins have been studied over the decades. In the past decade, our understanding of r-protein functions has advanced significantly due to intensive studies on ribosomes and gene expression mechanisms not only in model bacteria like Escherichia coli or Bacillus subtilis but also in little-explored bacterial species from various phyla. The aim of this review is to update information on the multiple functions of r-proteins in bacteria.

Diversity, Distribution, and Chromosomal Rearrangements of TRIP1 Repeat Sequences in Escherichia coli

The bacterial genome contains numerous repeated sequences that greatly affect its genomic plasticity. The Escherichia coli K-12 genome contains three copies of the TRIP1 repeat sequence (TRIP1a, TRIP1b, and TRIP1c). However, the diversity, distribution, and role of the TRIP1 repeat sequence in the E. coli genome are still unclear. In this study, after screening 6725 E. coli genomes, the TRIP1 repeat was found in the majority of E. coli strains (96%: 6454/6725). The copy number and direction of the TRIP1 repeat sequence varied in each genome. Overall, 2449 genomes (36%: 2449/6725) had three copies of TRIP1 (TRIP1a, TRIP1b, and TRIP1c), which is the same as E. coli K-12. Five types of TRIP1 repeats, including two new types (TRIP1d and TRIP1e), are identified in E. coli genomes, located in 4703, 3529, 5741, 1565, and 232 genomes, respectively. Each type of TRIP1 repeat is localized to a specific locus on the chromosome. TRIP1 repeats can cause intra-chromosomal rearrangements. A total of 156 rearrangement events were identified, of which 88% (137/156) were between TRIP1a and TRIP1c. These findings have important implications for future research on TRIP1 repeats.

Pithoviruses Are Invaded by Repeats That Contribute to Their Evolution and Divergence from Cedratviruses

Pithoviridae are amoeba-infecting giant viruses possessing the largest viral particles known so far. Since the discovery of Pithovirus sibericum, recovered from a 30,000-yr-old permafrost sample, other pithoviruses, and related cedratviruses, were isolated from various terrestrial and aquatic samples. Here, we report the isolation and genome sequencing of 2 Pithoviridae from soil samples, in addition to 3 other recent isolates. Using the 12 available genome sequences, we conducted a thorough comparative genomic study of the Pithoviridae family to decipher the organization and evolution of their genomes. Our study reveals a nonuniform genome organization in 2 main regions: 1 concentrating core genes and another gene duplications. We also found that Pithoviridae genomes are more conservative than other families of giant viruses, with a low and stable proportion (5% to 7%) of genes originating from horizontal transfers. Genome size variation within the family is mainly due to variations in gene duplication rates (from 14% to 28%) and massive invasion by inverted repeats. While these repeated elements are absent from cedratviruses, repeat-rich regions cover as much as a quarter of the pithoviruses genomes. These regions, identified using a dedicated pipeline, are hotspots of mutations, gene capture events, and genomic rearrangements that contribute to their evolution.

Molecular mechanism controlling anthocyanin composition and content in radish plants with different root colors

Radish (Raphanus sativus) roots exhibit various colors that reflect their anthocyanin compositions and contents. However, the details of the mechanism linking the expression of anthocyanin biosynthesis and their transcriptional regulators to anthocyanin composition in radish roots remained unknown. Here, we characterized the role of the anthocyanin biosynthetic enzyme flavonoid 3'-hydroxylase (RsF3'H), together with the R2R3 MYB transcription factor (TF) RsMYB1 and the basic helix-loop-helix (bHLH) TF TRANSPARENT TESTA 8 (RsTT8), in four radish plants with different root colors: white (W), deep red (DR), dark purple (DP), and dark greyish purple (DGP). The DR plant contained heterozygous for RsF3'H with low expression level and accumulated a large amount of pelargonidin, resulting in deep red color. While, the DP and DGP plants accumulated the cyanidin due to the higher expression level of functional RsF3'H. Notably, RsMYB1 and RsTT8 transcripts were abundant in all pigmented roots, but not in white roots. To investigate the differential expression of RsMYB1 and RsTT8, we compared the sequences of their promoter regions among the four radish plants, revealing variations in the numbers of cis-elements and in promoter architecture. Promoter activation assays demonstrated that variation in the RsMYB1 and RsTT8 promoters may contribute to the expression level of these genes, and RsMYB1 can activate its own expression as well as promote the RsTT8 expression. These results suggested that RsF3'H plays a vital role in anthocyanin composition and the expression level of both RsMYB1 and RsTT8 are crucial determinants for anthocyanin content in radish roots. Overall, these findings provide insight into the molecular basis of anthocyanin composition and level in radish roots.

Variability of Inverted Repeats in All Available Genomes of Bacteria

Noncanonical secondary structures in nucleic acids have been studied intensively in recent years. Important biological roles of cruciform structures formed by inverted repeats (IRs) have been demonstrated in diverse organisms, including humans. Using Palindrome analyser, we analyzed IRs in all accessible bacterial genome sequences to determine their frequencies, lengths, and localizations. IR sequences were identified in all species, but their frequencies differed significantly across various evolutionary groups. We detected 242,373,717 IRs in all 1,565 bacterial genomes. The highest mean IR frequency was detected in the Tenericutes (61.89 IRs/kbp) and the lowest mean frequency was found in the Alphaproteobacteria (27.08 IRs/kbp). IRs were abundant near genes and around regulatory, tRNA, transfer-messenger RNA (tmRNA), and rRNA regions, pointing to the importance of IRs in such basic cellular processes as genome maintenance, DNA replication, and transcription. Moreover, we found that organisms with high IR frequencies were more likely to be endosymbiotic, antibiotic producing, or pathogenic. On the other hand, those with low IR frequencies were far more likely to be thermophilic. This first comprehensive analysis of IRs in all available bacterial genomes demonstrates their genomic ubiquity, nonrandom distribution, and enrichment in genomic regulatory regions. IMPORTANCE Our manuscript reports for the first time a complete analysis of inverted repeats in all fully sequenced bacterial genomes. Thanks to the availability of unique computational resources, we were able to statistically evaluate the presence and localization of these important regulatory sequences in bacterial genomes. This work revealed a strong abundance of these sequences in regulatory regions and provides researchers with a valuable tool for their manipulation.

Concentration of inverted repeats along human DNA

This work aims to describe the observed enrichment of inverted repeats in the human genome; and to identify and describe, with detailed length profiles, the regions with significant and relevant enriched occurrence of inverted repeats. The enrichment is assessed and tested with a recently proposed measure (z-scores based measure). We simulate a genome using an order 7 Markov model trained with the data from the real genome. The simulated genome is used to establish the critical values which are used as decision thresholds to identify the regions with significant enriched concentrations. Several human genome regions are highly enriched in the occurrence of inverted repeats. This is observed in all the human chromosomes. The distribution of inverted repeat lengths varies along the genome. The majority of the regions with severely exaggerated enrichment contain mainly short length inverted repeats. There are also regions with regular peaks along the inverted repeats lengths distribution (periodic regularities) and other regions with exaggerated enrichment for long lengths (less frequent). However, adjacent regions tend to have similar distributions.

Widely spaced and divergent inverted repeats become a potent source of chromosomal rearrangements in long single-stranded DNA regions

DNA inverted repeats (IRs) are widespread across many eukaryotic genomes. Their ability to form stable hairpin/cruciform secondary structures is causative in triggering chromosome instability leading to several human diseases. Distance and sequence divergence between IRs are inversely correlated with their ability to induce gross chromosomal rearrangements (GCRs) because of a lesser probability of secondary structure formation and chromosomal breakage. In this study, we demonstrate that structural parameters that normally constrain the instability of IRs are overcome when the repeats interact in single-stranded DNA (ssDNA). We established a system in budding yeast whereby >73 kb of ssDNA can be formed in cdc13-707fs mutants. We found that in ssDNA, 12 bp or 30 kb spaced Alu-IRs show similarly high levels of GCRs, while heterology only beyond 25% suppresses IR-induced instability. Mechanistically, rearrangements arise after cis-interaction of IRs leading to a DNA fold-back and the formation of a dicentric chromosome, which requires Rad52/Rad59 for IR annealing as well as Rad1-Rad10, Slx4, Msh2/Msh3 and Saw1 proteins for nonhomologous tail removal. Importantly, using structural characteristics rendering IRs permissive to DNA fold-back in yeast, we found that ssDNA regions mapped in cancer genomes contain a substantial number of potentially interacting and unstable IRs.

Impacts of Molecular Structure on Nucleic Acid-Protein Interactions

Interactions between nucleic acids and proteins are some of the most important interactions in biology because they are the cornerstones for fundamental biological processes, such as replication, transcription, and recombination [...].

DNA Motifs and an Accessory CRISPR Factor Determine Cas1 Binding and Integration Activity in Sulfolobus islandicus

CRISPR-Cas systems empower prokaryotes with adaptive immunity against invasive mobile genetic elements. At the first step of CRISPR immunity adaptation, short DNA fragments from the invaders are integrated into CRISPR arrays at the leader-proximal end. To date, the mechanism of recognition of the leader-proximal end remains largely unknown. Here, in the Sulfolobus islandicus subtype I-A system, we show that mutations destroying the proximal region reduce CRISPR adaptation in vivo. We identify that a stem-loop structure is present on the leader-proximal end, and we demonstrate that Cas1 preferentially binds the stem-loop structure in vitro. Moreover, we demonstrate that the integrase activity of Cas1 is modulated by interacting with a CRISPR-associated factor Csa3a. When translocated to the CRISPR array, the Csa3a-Cas1 complex is separated by Csa3a binding to the leader-distal motif and Cas1 binding to the leader-proximal end. Mutation at the leader-distal motif reduces CRISPR adaptation efficiency, further confirming the in vivo function of leader-distal motif. Together, our results suggest a general model for binding of Cas1 protein to a leader motif and modulation of integrase activity by an accessory factor.