An Analysis of Single Nucleotide Substitution in Genetic Codons - Probabilities and Outcomes

High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2

Introduction

The propensity of nucleotide bases to form pairs, causes folding and the formation of secondary structure in the RNA. Therefore, purine (R): pyrimidine (Y) base-pairing is vital to maintain uniform lateral dimension in RNA secondary structure. Transversions or base substitutions between R and Y bases, are more detrimental to the stability of RNA secondary structure, than transitions derived from substitutions between A and G or C and T. The study of transversion and transition base substitutions is important to understand evolutionary mechanisms of RNA secondary structure in the 5'  and 3'  untranslated (UTR) regions of SARS-CoV-2. In this work, we carried out comparative analysis of transition and transversion base substitutions in the stem and loop regions of RNA secondary structure of SARS-CoV-2.

Methods

We have considered the experimentally determined and well documented stem and loop regions of 5' and 3' UTR regions of SARS-CoV-2 for base substitution analysis. The secondary structure comprising of stem and loop regions were visualized using the RNAfold web server. The GISAID repository was used to extract base sequence alignment of the UTR regions. Python scripts were developed for comparative analysis of transversion and transition frequencies in the stem and the loop regions.

Results

The results of base substitution analysis revealed a higher transition (ti) to transversion (tv) ratio (ti/tv) in the stem region of UTR of RNA secondary structure of SARS-CoV-2 reported during the early stage of the pandemic. The higher ti/tv ratio in the stem region suggested the influence of secondary structure in selecting the pattern of base substitutions. This differential pattern of ti/tv values between stem and loop regions was not observed among the Delta and Omicron variants that dominated the later stage of the pandemic. It is noteworthy that the ti/tv values in the stem and loop regions were similar among the later dominant Delta and Omicron variant strains which is to be investigated to understand the rapid evolution and global adaptation of SARS-CoV-2.

Conclusion

Our findings implicate the lower frequency of transversions than the transitions in the stem regions of UTRs of SARS-CoV-2. The RNA secondary structures are associated with replication, translation, and packaging, further investigations are needed to understand these base substitutions across different variants of SARS-CoV-2.

Uncovering the role of leucine 59 in Renilla luciferase stability and activity with error-prone PCR: Implications for protein engineering

A new variant of Renilla luciferase, named Met C-SRLuc 8, was obtained from a random mutagenesis library and expressed in Escherichia coli BL21 (DE3) plys and purified. The results of the enzyme's binding affinity, kinetic stability, and bioinformatic studies demonstrated that leucine 59, located within the hot-spot foldon in the N-terminal domain of the protein, plays a significant role in the stability and activity of Renilla luciferase. These findings may facilitate the engineering of different variants of this enzyme to achieve thermally stable versions for various biotechnological applications.

Continent-wide evolutionary trends of emerging SARS-CoV-2 variants: dynamic profiles from Alpha to Omicron

The ongoing SARS-CoV-2 evolution process has generated several variants due to its continuous mutations, making pandemics more critical. The present study illustrates SARS-CoV-2 evolution and its emerging mutations in five directions. First, the significant mutations in the genome and S-glycoprotein were analyzed in different variants. Three linear models were developed with the regression line to depict the mutational load for S-glycoprotein, total genome excluding S-glycoprotein, and whole genome. Second, the continent-wide evolution of SARS-CoV-2 and its variants with their clades and divergence were evaluated. It showed the region-wise evolution of the SARS-CoV-2 variants and their clustering event. The major clades for each variant were identified. One example is clade 21K, a major clade of the Omicron variant. Third, lineage dynamics and comparison between SARS-CoV-2 lineages across different countries are also illustrated, demonstrating dominant variants in various countries over time. Fourth, gene-wise mutation patterns and genetic variability of SARS-CoV-2 variants across various countries are illustrated. High mutation patterns were found in the ORF10, ORF6, S, and low mutation pattern E genes. Finally, emerging AA point mutations (T478K, L452R, N501Y, S477N, E484A, Q498R, and Y505H), their frequencies, and country-wise occurrence were identified, and the highest event of two mutations (T478K and L452R) was observed.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

Approximation algorithm for rearrangement distances considering repeated genes and intergenic regions

The rearrangement distance is a method to compare genomes of different species. Such distance is the number of rearrangement events necessary to transform one genome into another. Two commonly studied events are the transposition, which exchanges two consecutive blocks of the genome, and the reversal, which reverts a block of the genome. When dealing with such problems, seminal works represented genomes as sequences of genes without repetition. More realistic models started to consider gene repetition or the presence of intergenic regions, sequences of nucleotides between genes and in the extremities of the genome. This work explores the transposition and reversal events applied in a genome representation considering both gene repetition and intergenic regions. We define two problems called Minimum Common Intergenic String Partition and Reverse Minimum Common Intergenic String Partition. Using a relation with these two problems, we show a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta \left( k \right)$$\end{document}Θk-approximation for the Intergenic Transposition Distance, the Intergenic Reversal Distance, and the Intergenic Reversal and Transposition Distance problems, where k is the maximum number of copies of a gene in the genomes. Our practical experiments on simulated genomes show that the use of partitions improves the estimates for the distances.

Mutational and phylogenetic status of west siberian strains of BLV

The study is devoted of full-genome BLV sequences circulating in cattle populations of the Novosibirsk region, Russia. The phylogenetic tree shows that the West Siberian isolates are quite closely related to such previously isolated strains as AF399704 (Brazil), AP018007, AP018016, AP018019, LC007988, LC007991 (Japan) and EF065638 (Belgium) we calculations show that the number of mutations that could independently occur in parallel evolving BLV strains significantly exceeds the expected number based on the probability of corresponding substitutions. It was also found that the studied isolates have some mutations, the presence of which, at first glance, is possible only with their divergent development in different independently evolving branches. However, calculations show that the probability of an independent origin of an identical mutation is extremely small, which indicates the possibility of exchanging RNA sites between isolates circulating in West Siberian cattle populations.

CRISPR/Cas with ribonucleoprotein complexes and transiently selected telomere vectors allows highly efficient marker-free and multiple genome editing in Botrytis cinerea

CRISPR/Cas has become the state-of-the-art technology for genetic manipulation in diverse organisms, enabling targeted genetic changes to be performed with unprecedented efficiency. Here we report on the first establishment of robust CRISPR/Cas editing in the important necrotrophic plant pathogen Botrytis cinerea based on the introduction of optimized Cas9-sgRNA ribonucleoprotein complexes (RNPs) into protoplasts. Editing yields were further improved by development of a novel strategy that combines RNP delivery with cotransformation of transiently stable vectors containing telomeres, which allowed temporary selection and convenient screening for marker-free editing events. We demonstrate that this approach provides superior editing rates compared to existing CRISPR/Cas-based methods in filamentous fungi, including the model plant pathogen Magnaporthe oryzae. Genome sequencing of edited strains revealed very few additional mutations and no evidence for RNP-mediated off-targeting. The high performance of telomere vector-mediated editing was demonstrated by random mutagenesis of codon 272 of the sdhB gene, a major determinant of resistance to succinate dehydrogenase inhibitor (SDHI) fungicides by in bulk replacement of the codon 272 with codons encoding all 20 amino acids. All exchanges were found at similar frequencies in the absence of selection but SDHI selection allowed the identification of novel amino acid substitutions which conferred differential resistance levels towards different SDHI fungicides. The increased efficiency and easy handling of RNP-based cotransformation is expected to accelerate molecular research in B. cinerea and other fungi.

In this study, we describe the establishment of the CRISPR/Cas technology for genome editing in the gray mold fungus Botrytis cinerea, one of the economically most important plant pathogens worldwide. We report the development of a strategy which combines the introduction of an optimized nuclear-targeted Cas9-single guide RNA ribonucleoprotein complex (RNP) and a repair template together with unstable telomere vectors for transient selection into fungal protoplasts. A high proportion of the transformants contains the desired genetic changes, and the telomere vector is lost subsequently when selection is stopped. This system allowed introduction of changes into the genome without the requirement of selection markers. It shows superior editing efficiencies compared to existing CRISPR/Cas protocols for filamentous fungi, and leads to a very low number of additional off-target mutations. To demonstrate the performance of our protocol, we conducted for the first time a site-directed, random mutagenesis in a gene encoding an important fungicide target. This approach allows new applications such as in vivo structure-function analysis of proteins and rational fungicide resistance studies. As demonstrated with the rice blast pathogen Magnaporthe oryzae, the RNP-based CRISPR/Cas toolset with telomere vectors can be transferred to other fungi and is expected to boost their genetic manipulation.

Deep2Full: Evaluating strategies for selecting the minimal mutational experiments for optimal computational predictions of deep mutational scan outcomes

Performing a complete deep mutational scan with all single point mutations may not be practical, and may not even be required, especially if predictive computational models can be developed. Computational models are however naive to cellular response in the myriads of assay-conditions. In a realistic paradigm of assay context-aware predictive hybrid models that combine minimal experimental data from deep mutational scans with structure, sequence information and computational models, we define and evaluate different strategies for choosing this minimal set. We evaluated the trivial strategy of a systematic reduction in the number of mutational studies from 85% to 15%, along with several others about the choice of the types of mutations such as random versus site-directed with the same 15% data completeness. Interestingly, the predictive capabilities by training on a random set of mutations and using a systematic substitution of all amino acids to alanine, asparagine and histidine (ANH) were comparable. Another strategy we explored, augmenting the training data with measurements of the same mutants at multiple assay conditions, did not improve the prediction quality. For the six proteins we analyzed, the bin-wise error in prediction is optimal when 50-100 mutations per bin are used in training the computational model, suggesting that good prediction quality may be achieved with a library of 500-1000 mutations.

Biological Roles of Protein-Coding Tandem Repeats in the Yeast Candida Albicans

Tandem repeat (TR) DNA mutates faster than other DNA by insertion and deletion of repeats. Large parts of eukaryotic proteomes are encoded by ORFs containing protein-coding TRs (TR-ORFs, pcTRs) with largely unknown biological consequences. We explored these in the yeast Candida albicans, an opportunistic human pathogen. We found that almost half of C. albicans’ proteins are encoded by TR-ORFs. pcTR frequency differed only moderately between different gene (GO) categories. Bioinformatic predictions of genome-wide mutation rates and clade-specific differences in pcTR allele frequencies indicated that pcTRs (i) significantly increase the genome-wide mutation rate; (ii) significantly impact on fitness and (iii) allow the evolution of selectively advantageous clade-specific protein variants. Synonymous mutations reduced the repetitiveness of many amino acid repeat-encoding pcTRs. A survey, in 58 strains, revealed that in some pcTR regions in which repetitiveness was not significantly diminished by synonymous mutations the habitat predicted which alleles were present, suggesting roles of pcTR mutation in short-term adaptation and pathogenesis. In C. albicans pcTR mutation apparently is an important mechanism for mutational advance and possibly also rapid adaptation, with synonymous mutations providing a mechanism for adjusting mutation rates of individual pcTRs. Analyses of Arabidopsis and human pcTRs showed that the latter also occurs in other eukaryotes.