Mutational Bias Research Articles

Impact of the clinically approved BTK inhibitors on the conformation of full-length BTK and analysis of the development of BTK resistance mutations in chronic lymphocytic leukemia.

Inhibition of Bruton's tyrosine kinase (BTK) has proven to be highly effective in the treatment of B-cell malignancies such as chronic lymphocytic leukemia (CLL), autoimmune disorders and multiple sclerosis. Since the approval of the first BTK inhibitor (BTKi), Ibrutinib, several other inhibitors including Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib have been clinically approved. All are covalent active site inhibitors, with the exception of the reversible active site inhibitor Pirtobrutinib. The large number of available inhibitors for the BTK target creates challenges in choosing the most appropriate BTKi for treatment. Side-by-side comparisons in CLL have shown that different inhibitors may differ in their treatment efficacy. Moreover, the nature of the resistance mutations that arise in patients appears to depend on the specific BTKi administered. We have previously shown that Ibrutinib binding to the kinase active site causes unanticipated long-range effects on the global conformation of BTK (Joseph, R.E., et al., 2020, https://doi.org/10.7554/eLife.60470 ). Here we show that binding of each of the five approved BTKi to the kinase active site brings about distinct allosteric changes that alter the conformational equilibrium of full-length BTK. Additionally, we provide an explanation for the resistance mutation bias observed in CLL patients treated with different BTKi and characterize the mechanism of action of two common resistance mutations: BTK T474I and L528W.

Exploring the Relationship Between Gene Expression and Low-Frequency Somatic Mutations in Arabidopsis with Duplex Sequencing.

Intragenomic mutation rates can vary dramatically due to transcription-associated mutagenesis or transcription-coupled repair, which vary based on local epigenomic modifications that are nonuniformly distributed across genomes. One feature associated with decreased mutation is higher expression level, which depends on environmental cues. To understand the magnitude of expression-dependent mutation rate variation, we perturbed expression through a heat treatment in Arabidopsis thaliana. We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using duplex sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty. Somatic mutations in plants may be useful for understanding drivers of DNA damage and repair in the germline since plants experience late germline segregation and both somatic and germline cells share common repair machinery. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild-type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8 × 10-8 and 2.6 × 10-8, respectively), while MMR mutant frequencies were significantly elevated (1.13 × 10-6). Interestingly, in the MMR mutant lines, there was no difference in the somatic mutation frequencies between temperature treatments or between highly versus lowly expressed genes. The extremely low somatic variant frequencies in WT plants indicate that larger datasets will be needed to address fundamental evolutionary questions about whether environmental change leads to gene-specific changes in mutation rate.

AlzDiscovery: A computational tool to identify Alzheimer's disease-causing missense mutations using protein structure information.

Alzheimer's disease (AD) is one of the most common forms of dementia and neurodegenerative diseases, characterized by the formation of neuritic plaques and neurofibrillary tangles. Many different proteins participate in this complicated pathogenic mechanism, and missense mutations can alter the folding and functions of these proteins, significantly increasing the risk of AD. However, many methods to identify AD-causing variants did not consider the effect of mutations from the perspective of a protein three-dimensional environment. Here, we present a machine learning-based analysis to classify the AD-causing mutations from their benign counterparts in 21 AD-related proteins leveraging both sequence- and structure-based features. Using computational tools to estimate the effect of mutations on protein stability, we first observed a bias of the pathogenic mutations with significant destabilizing effects on family AD-related proteins. Combining this insight, we built a generic predictive model, and improved the performance by tuning the sample weights in the training process. Our final model achieved the performance on area under the receiver operating characteristic curve up to 0.95 in the blind test and 0.70 in an independent clinical validation, outperforming all the state-of-the-art methods. Feature interpretation indicated that the hydrophobic environment and polar interaction contacts were crucial to the decision on pathogenic phenotypes of missense mutations. Finally, we presented a user-friendly web server, AlzDiscovery, for researchers to browse the predicted phenotypes of all possible missense mutations on these 21 AD-related proteins. Our study will be a valuable resource for AD screening and the development of personalized treatment.

P53 mutation biases squamocolumnar junction progenitor cells towards dysplasia rather than metaplasia in Barrett’s oesophagus

BackgroundWhile p53 mutations occur early in Barrett’s oesophagus (BE) progression to oesophageal adenocarcinoma (EAC), their role in gastric cardia stem cells remains unclear.ObjectiveThis study investigates the impact of p53 mutation...

Whole genome regulatory effect of MoISW2 and consequences for the evolution of the rice plant pathogenic fungus Magnaporthe oryzae.

Isw2 proteins, ubiquitous across eukaryotes, exhibit a propensity for DNA binding and exert dynamic influences on local chromosome condensation in an ATP-dependent fashion, thereby modulating the accessibility of neighboring genes to transcriptional machinery. Here, we report the deletion of a putative MoISW2 gene, yielding substantial ramifications on plant pathogenicity. Subsequent gene complementation and chromatin immunoprecipitation sequencing (ChIP-seq) analyses were conducted to delineate binding sites. RNA sequencing (RNA-seq) assays revealed discernible impacts on global gene regulation along chromosomes in both mutant and wild-type strains, with comparative analyses against 55 external RNA-seq data sets corroborating these findings. Notably, MoIsw2-mediated binding and activities delineate genomic loci characterized by pronounced gene expression variability proximal to MoIsw2 binding sites, juxtaposed with comparatively stable expression in surrounding regions. The contingent genes influenced by MoIsw2 activity predominantly encompass niche-determinant genes, including those encoding secreted proteins, secondary metabolites, and stress-responsive elements, alongside avirulence genes. Furthermore, our investigations unveil a spatial correlation between MoIsw2 binding motifs and known transposable elements (TEs), suggesting a potential interplay wherein TE transposition at these loci could modulate the transcriptional landscape of Magnaporthe oryzae in a strain-specific manner. Collectively, these findings position MoIsw2 as a plausible master regulator orchestrating the delicate equilibrium between genes vital for biomass proliferation, akin to housekeeping genes, and niche-specific determinants crucial for ecological adaptability. Stress-induced TE transposition, in conjunction with MoIsw2 activity, emerges as a putative mechanism fostering enhanced mutagenesis and accelerated evolution of niche-determinant genes relative to housekeeping counterparts.IMPORTANCEIsw2 proteins are conserved in plants, fungi, animals, and other eukaryotes. We show that a fungal Isw2 protein in the rice pathogen Magnaporthe oryzae binds to retrotransposon (RT) DNA motifs and affects the epigenetic gene expression landscape of the fungal genome. Mainly ecological niche determinant genes close to the binding motifs are affected. RT elements occur frequently in DNA between genes in most organisms. They move place and multiply in the genome, especially under physiological stress. We further discuss the Isw2 and RT combined activities as a possible sought-after mechanism that can cause biased mutation rates and faster evolution of genes necessary for reacting to abiotic and biotic challenges. The most important biotic challenges for plant pathogens are the ones from the host plants' innate immunity. The overall result of these combined activities will be an adaptation-directed evolution of niche-determinant genes.

Abstract IA-07: RAS mutation tropism: insight into tumor initation

Abstract How we get cancer, a process termed tumor initiation, remains an enigma, which has crucial clinical implications for early detection, prevention, and interception. Backtracking to catch the moment a specific mutagenic event in a single gene from an individual cell occurred decades before manifesting as cancer is challenging. However, the RAS family of small GTPases thate are commonly mutated in cancer provide some insight. Despite over 50 possible oncogenic RAS mutations, each cancer type exhibits a predilection or ‘tropism’ towards a distinct and often unique subset of these mutations. For example, most pancreatic ductal adenocarcinomas are driven by a G12D/V mutation in KRAS. As RAS mutations are often initiating, this tropism ostensibly reflects the fundamental biology underlying tumor initiation, which becomes embedded in the genome of the resulting tumors. Fortunately, the process of RAS mutation tropism can be studied in mice exposed to the environmental carcinogen urethane. Such carcinogenesis both captures the moment of tumor initiation (the point of carcinogen exposure) and displays RAS mutation tropism (lung lesions are induced by a specific RAS mutation). By using ultra-sensitive sequencing to detect mutations immediately after urethane exposure in different mouse genetic backgrounds, we find that the mutagenic bias of this carcinogen influences the type of RAS mutations generated. However, biological selection for an optimal or ‘sweet spot’ of oncogenic signaling ultimately appears to dictate which mutations go on to initiate tumorigenesis. Analysis of sequence datasets from human lung cancer suggests a similar dominant role for biological selection in the RAS mutation tropism of this cancer. In summary, RAS mutation tropism appears to be a combination of mutational bias, and even more so, biological selection, providing critical insights into the origins of cancers like pancreatic. Citation Format: Christopher M Counter. RAS mutation tropism: insight into tumor initation [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Pancreatic Cancer Research; 2024 Sep 15-18; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(17 Suppl_2):Abstract nr IA-07.

Conservation of mutation and recombination parameters between mammals and zebra finch.

Most of our understanding of the fundamental processes of mutation and recombination stems from a handful of disparate model organisms and pedigree studies of mammals, with little known about other vertebrates. To gain a broader comparative perspective, we focused on the zebra finch (Taeniopygia castanotis), which, like other birds, differs from mammals in its karyotype (which includes many micro-chromosomes), in the mechanism by which recombination is directed to the genome, and in aspects of ontogenesis. We collected genome sequences from three generation pedigrees that provide information about 80 meioses, inferring 202 single-point de novo mutations, 1,174 crossovers, and 275 non-crossovers. On that basis, we estimated a sex-averaged mutation rate of 5.0 × 10-9 per base pair per generation, on par with mammals that have a similar generation time. Also as in mammals, we found a paternal germline mutation bias at later stages of gametogenesis (of 1.7 to 1) but no discernible difference between sexes in early development. We also examined recombination patterns, and found that the sex-averaged crossover rate on macro-chromosomes (1.05 cM/Mb) is again similar to values observed in mammals, as is the spatial distribution of crossovers, with a pronounced enrichment near telomeres. In contrast, non-crossover rates are more uniformly distributed. On micro-chromosomes, sex-averaged crossover rates are substantially higher (4.21 cM/Mb), as expected from crossover homeostasis, and both crossover and non-crossover events are more uniformly distributed. At a finer scale, recombination events overlap CpG islands more often than expected by chance, as expected in the absence of PRDM9. Despite differences in the mechanism by which recombination events are specified and the presence of many micro-chromosomes, estimates of the degree of GC-biased gene conversion (59%), the mean non-crossover conversion tract length (~23 bp), and the non-crossover to crossover ratio (6.7:1) are all comparable to those reported in primates and mice. The conservation of mutation and recombination properties from zebra finch to mammals suggest that these processes have evolved under stabilizing selection.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria.

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

Is biased mutation sufficient to save runaway sexual selection?

In the 1980s, groundbreaking theoretical studies showed that ornaments displayed during courtship can coevolve with preferences for such ornaments, leading to extreme exaggeration of both traits. Later models cast doubt on such "runaway" sexual selection, showing that even a small cost of preferences can prevent exaggerated ornaments from persisting long-term. It was subsequently shown that if mutations acting on the ornament are biased-tending to produce smaller rather than larger ornaments-then exaggeration can persist even in the presence of preference costs, seemingly vindicating the original models. Here, we unpack an implicit assumption of these "biased mutation" models: Mutations are assumed to lead, on average, to both smaller and less costly ornaments. Biased mutation consequently generates both a fitness cost (due to reduced mating success) and a fitness benefit (due to increased survival). We lift this assumption by separating an individual's investment in an ornament from their efficiency in converting such investment into ornament size. We assume that biased mutation acts only on efficiency but not on investment, and discuss the plausibility of this alternative assumption. Our model predicts that exaggerated ornaments and preferences can persist stably once they arise, but that strong initial preferences are needed to kick-start the runaway process. Consequently, biased mutation alone may not always be sufficient to save runaway sexual selection.

Rare variations within the serine/arginine-rich splicing factor PtoRSZ21 modulate stomatal size to determine drought tolerance in Populus.

Rare variants contribute significantly to the 'missing heritability' of quantitative traits. The genome-wide characteristics of rare variants and their roles in environmental adaptation of woody plants remain unexplored. Utilizing genome-wide rare variant association study (RVAS), expression quantitative trait loci (eQTL) mapping, genetic transformation, and molecular experiments, we explored the impact of rare variants on stomatal morphology and drought adaptation in Populus. Through comparative analysis of five world-wide Populus species, we observed the influence of mutational bias and adaptive selection on the distribution of rare variants. RVAS identified 75 candidate genes correlated with stomatal size (SS)/stomatal density (SD), and a rare haplotype in the promoter of serine/arginine-rich splicing factor PtoRSZ21 emerged as the foremost association signal governing SS. As a positive regulator of drought tolerance, PtoRSZ21 can recruit the core splicing factor PtoU1-70K to regulate alternative splicing (AS) of PtoATG2b (autophagy-related 2). The rare haplotype PtoRSZ21hap2 weakens binding affinity to PtoMYB61, consequently affecting PtoRSZ21 expression and SS, ultimately resulting in differential distribution of Populus accessions in arid and humid climates. This study enhances the understanding of regulatory mechanisms that underlie AS induced by rare variants and might provide targets for drought-tolerant varieties breeding in Populus.

Homopolymer switches mediate adaptive mutability in mismatch repair-deficient colorectal cancer

Mismatch repair (MMR)-deficient cancer evolves through the stepwise erosion of coding homopolymers in target genes. Curiously, the MMR genes MutS homolog 6 (MSH6) and MutS homolog 3 (MSH3) also contain coding homopolymers, and these are frequent mutational targets in MMR-deficient cancers. The impact of incremental MMR mutations on MMR-deficient cancer evolution is unknown. Here we show that microsatellite instability modulates DNA repair by toggling hypermutable mononucleotide homopolymer runs in MSH6 and MSH3 through stochastic frameshift switching. Spontaneous mutation and reversion modulate subclonal mutation rate, mutation bias and HLA and neoantigen diversity. Patient-derived organoids corroborate these observations and show that MMR homopolymer sequences drift back into reading frame in the absence of immune selection, suggesting a fitness cost of elevated mutation rates. Combined experimental and simulation studies demonstrate that subclonal immune selection favors incremental MMR mutations. Overall, our data demonstrate that MMR-deficient colorectal cancers fuel intratumor heterogeneity by adapting subclonal mutation rate and diversity to immune selection.

Transcriptomic Analysis Reveals Adaptive Evolution and Conservation Implications for the Endangered Magnolia lotungensis.

Magnolia lotungensis is an extremely endangered endemic tree in China. To elucidate the genetic basis of M. lotungensis, we performed a comprehensive transcriptome analysis using a sample integrating the plant's bark, leaves, and flowers. De novo transcriptome assembly yielded 177,046 transcripts and 42,518 coding sequences. Notably, we identified 796 species-specific genes enriched in organelle gene regulation and defense responses. A codon usage bias analysis revealed that mutation bias appears to be the primary driver of selection in shaping the species' genetic architecture. An evolutionary analysis based on dN/dS values of paralogous and orthologous gene pairs indicated a predominance of purifying selection, suggesting strong evolutionary constraints on most genes. A comparative transcriptomic analysis with Magnolia sinica identified approximately 1000 ultra-conserved genes, enriched in essential cellular processes such as transcriptional regulation, protein synthesis, and genome stability. Interestingly, only a limited number of 511 rapidly evolving genes under positive selection were detected compared to M. sinica and Magnolia kuangsiensis. These genes were enriched in metabolic processes associated with adaptation to specific environments, potentially limiting the species' ability to expand its range. Our findings contribute to understanding the genetic architecture of M. lotungensis and suggest that an insufficient number of adaptive genes contribute to its endangered status.

Comparative Study on Codon Usage Patterns across Chloroplast Genomes of Eighteen Taraxacum Species

This study investigates codon usage bias within the chloroplast genomes of 18 Taraxacum species, focusing on the base composition and various metrics including GC content, Relative Synonymous Codon Usage (RSCU), Effective Number of Codons (ENc), and GC3s. Our analysis revealed a pronounced preference for A/T-ending codons across Taraxacum species, with GC content across the first, second, and third positions of the codons (GC1, GC2, GC3) and the average GC content consistently below 50%. A detailed examination using the RSCU metric identified 29 commonly preferred A/T-ending codons, indicating a strong codon usage bias towards these endings. Specifically, the codon for leucine (UUA) emerged as highly preferred, while the codon for serine (AGC) was least favored. Through the ENc–GC3s plot analysis, we explored the forces shaping this bias, finding evidence that both mutation pressure and natural selection significantly influence codon preference, with most coding sequences showing weak bias. The PR2 plot analysis further confirmed the role of these factors by demonstrating a higher frequency of T over A and C over G at the third codon position, pointing towards a mutation bias complemented by natural selection. Collectively, our findings highlight a consistent pattern of codon usage bias in the chloroplast genomes of Taraxacum species, influenced by a combination of mutation pressure, natural selection, and possibly other unidentified factors.

The Patterns of Codon Usage between Chordates and Arthropods are Different but Co-evolving with Mutational Biases.

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.

Spontaneous mutations and mutational responses to penicillin treatment in the bacterial pathogen Streptococcus pneumoniae D39

Bacteria with functional DNA repair systems are expected to have low mutation rates due to strong natural selection for genomic stability. However, our study of the wild-type Streptococcus pneumoniae D39, a pathogen responsible for many common diseases, revealed a high spontaneous mutation rate of 0.02 per genome per cell division in mutation-accumulation (MA) lines. This rate is orders of magnitude higher than that of other non-mutator bacteria and is characterized by a high mutation bias in the A/T direction. The high mutation rate may have resulted from a reduction in the overall efficiency of selection, conferred by the tiny effective population size in nature. In line with this, S. pneumoniae D39 also exhibited the lowest DNA mismatch-repair (MMR) efficiency among bacteria. Treatment with the antibiotic penicillin did not elevate the mutation rate, as penicillin did not induce DNA damage and S. pneumoniae lacks a stress response pathway. Our findings suggested that the MA results are applicable to within-host scenarios and provide insights into pathogen evolution.

Re-examining Correlations Between Synonymous Codon Usage and Protein Bond Angles in Escherichia coli.

Rosenberg AA, Marx A, Bronstein AM (Codon-specific Ramachandran plots show amino acid backbone conformation depends on identity of the translated codon. Nat Commun. 2022:13:2815) recently found a surprising correlation between synonymous codon usage and the dihedral bond angles of the resulting amino acid. However, their analysis did not account for the strongest known correlate of codon usage: gene expression. We re-examined the relationship between bond angles and codon usage by applying the approach of Rosenberg et al. to simulated protein-coding sequences that (i) have random codon usage, (ii) codon usage determined by mutation biases, and (iii) maintain the general relationship between codon usage and gene expression via the assumption of selection-mutation-drift equilibrium. We observed correlations between dihedral bond angle and codon usage when codon usage is entirely random, indicating possible conflation of noise with differences in bond angle distributions between synonymous codons. More relevant to the general analysis of codon usage patterns, we found surprisingly good agreement between the analysis of the real sequences and the analysis of sequences simulated assuming selection-mutation-drift equilibrium, with 91% of significant synonymous codon pairs detected in the former were also detected in the latter. We believe the correlation between codon usage and dihedral bond angles resulted from the variation in codon usage across genes due to the interplay between mutation bias, natural selection for translation efficiency, and gene expression, further underscoring these factors must be controlled for when looking for novel patterns related to codon usage.

Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer.

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Correction of non-random mutational biases along a linear bacterial chromosome by the mismatch repair endonuclease NucS.

The linear chromosome of Streptomyces exhibits a highly compartmentalized structure with a conserved central region flanked by variable arms. As double strand break (DSB) repair mechanisms play a crucial role in shaping the genome plasticity of Streptomyces, we investigated the role of EndoMS/NucS, a recently characterized endonuclease involved in a non-canonical mismatch repair (MMR) mechanism in archaea and actinobacteria, that singularly corrects mismatches by creating a DSB. We showed that Streptomyces mutants lacking NucS display a marked colonial phenotype and a drastic increase in spontaneous mutation rate. In vitro biochemical assays revealed that NucS cooperates with the replication clamp to efficiently cleave G/T, G/G and T/T mismatched DNA by producing DSBs. These findings are consistent with the transition-shifted mutational spectrum observed in the mutant strains and reveal that NucS-dependent MMR specific task is to eliminate G/T mismatches generated by the DNA polymerase during replication. Interestingly, our data unveil a crescent-shaped distribution of the transition frequency from the replication origin towards the chromosomal ends, shedding light on a possible link between NucS-mediated DSBs and Streptomyces genome evolution.

Molecular insights into codon usage analysis of mitochondrial fission and fusion gene: relevance to neurodegenerative diseases.

Mitochondrial dysfunction is the leading cause of neurodegenerative disorders like Alzheimer's disease and Parkinson's disease. Mitochondria is a highly dynamic organelle continuously undergoing the process of fission and fusion for even distribution of components and maintaining proper shape, number, and bioenergetic functionality. A set of genes governs the process of fission and fusion. OPA1, Mfn1, and Mfn2 govern fusion, while Drp1, Fis1, MIEF1, and MIEF2 genes control fission. Determination of specific molecular patterns of transcripts of these genes revealed the impact of compositional constraints on selecting optimal codons. AGA and CCA codons were over-represented, and CCC, GTC, TTC, GGG, ACG were under-represented in the fusion gene set. In contrast, CTG was over-represented, and GCG, CCG, and TCG were under-represented in the fission gene set. Hydropathicity analysis revealed non-polar protein products of both fission and fusion gene set transcripts. AGA codon repeats are an integral part of translational regulation machinery and present a distinct pattern of over-representation and under-representation in different transcripts within the gene sets, suggestive of selective translational force precisely controlling the occurrence of the codon. Out of six synonymous codons, five synonymous codons encoding for leucine were used differently in both gene sets. Hence, forces regulating the occurrence of AGA and five synonymous leucine-encoding codons suggest translational selection. A correlation of mutational bias with gene expression and codon bias and GRAVY and AROMA signifies the selection pressure in both gene sets, while the correlation of compositional bias with gene expression, codon bias, protein properties, and minimum free energy signifies the presence of compositional constraints. More than 25% of codons of both gene sets showed a significant difference in codon usage. The overall analysis shed light on molecular features of gene sets involved in fission and fusion.