CAIcal: a combined set of tools to assess codon usage adaptation

High-level expression of codon-optimized Taq DNA polymerase under the control of rhaBAD promoter

A DNA polymerase from Thermus aquaticus remains the most popular among DNA polymerases. It was widely applied in various fields involving the application of polymerase chain reaction (PCR), implying the high commercial value of this enzyme. For this reason, an attempt to obtain a high yield of Taq DNA polymerase is continuously conducted. In this study, the l-rhamnose-inducible promoter rhaBAD was utilized due to its ability to produce recombinant protein under tight control in E. coli expression system. Instead of full-length Taq polymerase, an N-terminal deletion of Taq polymerase was selected. To obtain a high-level expression, we attempted to optimize the codon by reducing the rare codon and GC content, and in a second attempt, we optimized the culture conditions for protein expression. The production of Taq polymerase using the optimum culture condition improved the level of expression by up to 3-fold. This approach further proved that a high level of recombinant protein expression could be achieved by yielding a purified Taq polymerase of about 8.5 mg/L of culture. This is the first research publication on the production of Taq polymerase with N-terminal deletion in E. coli with the control of the rhaBAD promoter system.

Combinative workflow for mRNA vaccine development

Recently, mRNA has gained a lot of attention in the field of vaccines, gene therapy, and protein replacement therapies. Herein, we are demonstrating a comprehensive approach to designing, cloning, and characterizing an antigenic cassette for the development of mRNA vaccine for COVID-19. The gene encoding the antigenic spike protein of the SARS-CoV-2 Omicron variant (B.1.1.529) was designed using the databases, characterized by in-silico tools, and assembled using overlapping oligonucleotide-based assembly by PCR. Next, the gene was cloned, mRNA was synthesized, and characterized using orthogonal approaches (Capillary electrophoresis, Sanger DNA sequencing, Next-generation sequencing, HPLC, qPCR, etc.). Furthermore, the antigen expression was monitored in-vitro using an animal cell model by western blot, flow cytometer, and surface plasmon resonance. The demonstrated approach has also been followed for developing the mRNA vaccines for various other indications such as Malaria, Herpes, Dengue, HPV, etc.

Determination of codon pattern and evolutionary forces acting on genes linked to inflammatory bowel disease

Inflammatory bowel disease (IBD) is an inflammatory disorder of the gastrointestinal tract. The present study attempted to understand the codon usage preferences in genes associated with IBD progression. Compositional analysis, codon usage bias (CUB), Relative synonymous codon usage (RSCU), RNA structure, and expression analysis were performed to obtain a comprehensive picture of codon usage in IBD genes. Compositional analysis of 62 IBD-associated genes revealed that G and T are the most and least abundant nucleotides, respectively. ApG, CpA, and TpG dinucleotides were overrepresented or randomly used, while ApC, CpG, GpT, and TpA dinucleotides were either underrepresented or randomly used in genes related to IBD. The codons influencing the codon usage the most in IBD genes were CGC and AGG. A comparison of codon usage between IBD, and pancreatitis (non-IBD inflammatory disease) indicated that only codon CTG codon usage was significantly different between IBD and pancreatitis. At the same time, there were codons ATA, ACA, CGT, CAA, GTA, CCT, ATT, GCT, CGG, TTG, and CAG for whom codon usage was significantly different for IBD and housekeeping gene sets. The results suggest similar codon usage in at least two inflammatory disorders, IBD and pancreatitis. The analysis helps understand the codon biology, factors affecting gene expression of IBD-associated genes, and the evolution of these genes. The study helps reveal the molecular patterns associated with IBD.

Revealing Molecular Patterns of Alzheimer's Disease Risk Gene Expression Signatures in COVID-19 Brains

Background

Various virus infections are known to predispose to Alzheimer's disease (AD), and a linkage between COVID-19 and AD has been established. COVID-19 infection modulates the gene expression of the genes implicated in progression of AD.

Objective

Determination of molecular patterns and codon usage and context analysis for the genes that are modulated during COVID-19 infection and are implicated in AD was the target of the study.

Methods

Our study employed a comprehensive array of research methods, including relative synonymous codon usage, Codon adaptation index analysis, Neutrality and parity analysis, Rare codon analyses, and codon context analysis. This meticulous approach was crucial in determining the molecular patterns present in genes up or downregulated during COVID-19 infection.

Results

G/C ending codons were preferred in upregulated genes while not in downregulated genes, and in both gene sets, longer genes have high expressivity. Similarly, T over A nucleotide was preferred, and selection was the major evolutionary force in shaping codon usage in both gene sets. Apart from stops codons, codons CGU - Arg, AUA - Ile, UUA - Leu, UCG - Ser, GUA - Val, and CGA - Arg in upregulated genes, while CUA - Leu, UCG - Ser, and UUA - Leu in downregulated genes were present below the 0.5%. Glutamine-initiated codon pairs have high residual values in upregulated genes. Identical codon pairs GAG-GAG and GUG-GUG were preferred in both gene sets.

Conclusions

The shared and unique molecular features in the up- and downregulated gene sets provide insights into the complex interplay between COVID-19 infection and AD. Further studies are required to elucidate the relationship of these molecular patterns with AD pathology.

An immuno-informatics approach for annotation of hypothetical proteins and multi-epitope vaccine designed against the Mpox virus

A worrying new outbreak of Monkeypox (Mpox) in humans is caused by the Mpox virus (MpoxV). The pathogen has roughly 28 hypothetical proteins of unknown structure, function, and pathogenicity. Using reliable bioinformatics tools, we attempted to analyze the MpoxV genome, identify the role of hypothetical proteins (HPs), and design a potential candidate vaccine. Out of 28, we identified seven hypothetical proteins using multi-server validation with high confidence for the occurrence of conserved domains. Their physical, chemical, and functional characterizations, including molecular weight, theoretical isoelectric point, 3D structures, GRAVY value, subcellular localization, functional motifs, antigenicity, and virulence factors, were performed. We predicted possible cytotoxic T cell (CTL), helper T cell (HTL) and linear and conformational B cell epitopes, which were combined in a 219 amino acid multiepitope vaccine with human β defensin as a linker. This multi-epitopic vaccine was structurally modelled and docked with toll-like receptor-3 (TLR-3). The dynamical stability of the vaccine-TLR-3 docked complexes exhibited stable interactions based on RMSD and RMSF tests. Additionally, the modelled vaccine was cloned in-silico in an E. coli host to check the appropriate expression of the final vaccine built. Our results might conform to an immunogenic and safe vaccine, which would require further experimental validation.Communicated by Ramaswamy H. Sarma.

Codon usage bias of goose circovirus and its adaptation to host

Goose circovirus (GoCV), a potential immunosuppressive virus possessing a circular single-stranded DNA genome, is widely distributed in both domesticated and wild geese. This virus infection causes significant economic losses in the waterfowl industry. The codon usage patterns of viruses reflect the evolutionary history and genetic architecture, allowing them to adapt quickly to changes in the external environment, particularly to their hosts. In this study, we retrieved the coding sequences (Rep and Cap) and the genome of GoCV from GenBank, conducting comprehensive research to explore the codon usage patterns in 144 GoCV strains. The overall codon usage of the GoCV strains was relatively similar and exhibited a slight bias. The effective number of codons (ENC) indicated a low overall extent of codon usage bias (CUB) in GoCV. Combined with the base composition and relative synonymous codon usage (RSCU) analysis, the results revealed a bias toward A- and G-ending codons in the overall codon usage. Analysis of the ENC-GC3s plot and neutrality plot suggested that natural selection plays an important role in shaping the codon usage pattern of GoCV, with mutation pressure having a minor influence. Furthermore, the correlations between ENC and relative indices, as well as correspondence analysis (COA), showed that hydrophobicity and geographical distribution also contribute to codon usage variation in GoCV, suggesting the possible involvement of natural selection. In conclusion, GoCV exhibits comparatively slight CUB, with natural selection being the major factor shaping the codon usage pattern of GoCV. Our research contributes to a deeper understanding of GoCV evolution and its host adaptation, providing valuable insights for future basic studies and vaccine design related to GoCV.

An efficient plasmid-based system for the recovery of recombinant vesicular stomatitis virus encoding foreign glycoproteins

Viral glycoproteins mediate entry into host cells, thereby dictating host range and pathogenesis. In addition, they constitute the principal target of neutralizing antibody responses, making them important antigens in vaccine development. Recombinant vesicular stomatitis virus (VSV) encoding foreign glycoproteins can provide a convenient and safe surrogate system to interrogate the function, evolution, and antigenicity of viral glycoproteins from viruses that are difficult to manipulate or those requiring high biosafety level containment. However, the production of recombinant VSV can be technically challenging. In this work, we present an efficient and robust plasmid-based system for the production of recombinant VSV encoding foreign glycoproteins. We validate the system using glycoproteins from different viral families, including arenaviruses, coronaviruses, and hantaviruses, as well as highlight their utility for studying the effects of mutations on viral fitness. Overall, the methods described herein can facilitate the study of both native and recombinant VSV encoding foreign glycoproteins and can serve as the basis for the production of VSV-based vaccines.

Orb2 enables rare-codon-enriched mRNA expression during Drosophila neuron differentiation

Regulation of codon optimality is an increasingly appreciated layer of cell- and tissue-specific protein expression control. Here, we use codon-modified reporters to show that differentiation of Drosophila neural stem cells into neurons enables protein expression from rare-codon-enriched genes. From a candidate screen, we identify the cytoplasmic polyadenylation element binding (CPEB) protein Orb2 as a positive regulator of rare-codon-dependent mRNA stability in neurons. Using RNA sequencing, we reveal that Orb2-upregulated mRNAs in the brain with abundant Orb2 binding sites have a rare-codon bias. From these Orb2-regulated mRNAs, we demonstrate that rare-codon enrichment is important for mRNA stability and social behavior function of the metabotropic glutamate receptor (mGluR). Our findings reveal a molecular mechanism by which neural stem cell differentiation shifts genetic code regulation to enable critical mRNA stability and protein expression.

Identification and Genomic Characterization of Two Novel Hepatoviruses in Shrews from Yunnan Province, China

Hepatitis A virus (HAV), a member of the genus Hepatovirus (Picornaviridae HepV), remains a significant viral pathogen, frequently causing enterically transmitted hepatitis worldwide. In this study, we conducted an epidemiological survey of HepVs carried by small terrestrial mammals in the wild in Yunnan Province, China. Utilizing HepV-specific broad-spectrum RT-PCR, next-generation sequencing (NGS), and QNome nanopore sequencing (QNS) techniques, we identified and characterized two novel HepVs provisionally named EpMa-HAV and EpLe-HAV, discovered in the long-tailed mountain shrew (Episoriculus macrurus) and long-tailed brown-toothed shrew (Episoriculus leucops), respectively. Our sequence and phylogenetic analyses of EpMa-HAV and EpLe-HAV indicated that they belong to the species Hepatovirus I (HepV-I) clade II, also known as the Chinese shrew HepV clade. Notably, the codon usage bias pattern of novel shrew HepVs is consistent with that of previously identified Chinese shrew HepV. Furthermore, our structural analysis demonstrated that shrew HepVs differ from other mammalian HepVs in RNA secondary structure and exhibit variances in key protein sites. Overall, the discovery of two novel HepVs in shrews expands the host range of HepV and underscores the existence of genetically diverse animal homologs of human HAV within the genus HepV.

Complete mitogenome and intra-family comparative mitogenomics showed distinct position of Pama Croaker Otolithoides pama

The Pama Croaker, Otolithoides pama, is an economically important fish species in Bangladesh. Intra-family similarities in morphology and typical barcode sequences of cox1 create ambiguities in its identification. Therefore, morphology and the complete mitochondrial genome of O. pama, and comparative mitogenomics within the family Sciaenidae have been studied. Extracted genomic DNA was subjected to Illumina-based short read sequencing for De-Novo mitogenome assembly. The complete mitogenome of O. pama (Accession: OQ784575.1) was 16,513 bp, with strong AC biasness and strand asymmetry. Relative synonymous codon usage (RSCU) among 13 protein-coding genes (PCGs) of O. pama was also analyzed. The studied mitogenomes including O. pama exhibited consistent sizes and gene orders, except for the genus Johnius which possessed notably longer mitogenomes with unique gene rearrangements. Different genetic distance metrics across 30 species of Sciaenidae family demonstrated 12S rRNA and the control region (CR) as the most conserved and variable regions, respectively, while most of the PCGs undergone a purifying selection. Different phylogenetic trees were congruent with one another, where O. pama was distinctly placed. This study would contribute to distinguishing closely related fish species of Sciaenidae family and can be instrumental in conserving the genetic diversity of O. pama.

Mitogenomic Insights into the Evolution, Divergence Time, and Ancestral Ranges of Coturnix Quails

The Old-World quails, Coturnix coturnix (common quail) and Coturnix japonica (Japanese quail), are morphologically similar yet occupy distinct geographic ranges. This study aimed to elucidate their evolutionary trajectory and ancestral distribution patterns through a thorough analysis of their mitochondrial genomes. Mitogenomic analysis revealed high structural conservation, identical translational mechanisms, and similar evolutionary pressures in both species. Selection analysis revealed significant evidence of positive selection across the Coturnix lineage for the nad4 gene tree owing to environmental changes and acclimatization requirements during its evolutionary history. Divergence time estimations imply that diversification among Coturnix species occurred in the mid-Miocene (13.89 Ma), and their current distributions were primarily shaped by dispersal rather than global vicariance events. Phylogenetic analysis indicates a close relationship between C. coturnix and C. japonica, with divergence estimated at 2.25 Ma during the Pleistocene epoch. Ancestral range reconstructions indicate that the ancestors of the Coturnix clade were distributed over the Oriental region. C. coturnix subsequently dispersed to Eurasia and Africa, and C. japonica to eastern Asia. We hypothesize that the current geographic distributions of C. coturnix and C. japonica result from their unique dispersal strategies, developed to evade interspecific territoriality and influenced by the Tibetan Plateau's geographic constraints. This study advances our understanding of the biogeographic and evolutionary processes leading to the diversification of C. coturnix and C. japonica, laying important groundwork for further research on this genus.

Phylogeny and evolution of SARS-CoV-2 during Delta and Omicron variant waves in India

SARS-CoV-2 evolution has continued to generate variants, responsible for new pandemic waves locally and globally. Varying disease presentation and severity has been ascribed to inherent variant characteristics and vaccine immunity. This study analyzed genomic data from 305 whole genome sequences from SARS-CoV-2 patients before and through the third wave in India. Delta variant was reported in patients without comorbidity (97%), while Omicron BA.2 was reported in patients with comorbidity (77%). Tissue adaptation studies brought forth higher propensity of Omicron variants to bronchial tissue than lung, contrary to observation in Delta variants from Delhi. Study of codon usage pattern distinguished the prevalent variants, clustering them separately, Omicron BA.2 isolated in February grouped away from December strains, and all BA.2 after December acquired a new mutation S959P in ORF1b (44.3% of BA.2 in the study) indicating ongoing evolution. Loss of critical spike mutations in Omicron BA.2 and gain of immune evasion mutations including G142D, reported in Delta but absent in BA.1, and S371F instead of S371L in BA.1 could explain very brief period of BA.1 in December 2021, followed by complete replacement by BA.2. Higher propensity of Omicron variants to bronchial tissue, probably ensured increased transmission while Omicron BA.2 became the prevalent variant possibly due to evolutionary trade-off. Virus evolution continues to shape the epidemic and its culmination.Communicated by Ramaswamy H. Sarma.

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Biodegradation of phthalate acid esters and whole-genome analysis of a novel Streptomyces sp. FZ201 isolated from natural habitats

Di-n-butyl phthalate (DBP) is one of the most extensively used phthalic acid esters (PAEs) and is considered to be an emerging, globally concerning pollutant. The genus Streptomyces holds promise as a degrader of various organic pollutants, but PAE biodegradation mechanisms by Streptomyces species remain unsolved. In this study, a novel PAE-degrading Streptomyces sp. FZ201 isolated from natural habitats efficiently degraded various PAEs. FZ201 had strong resilience against DBP and exhibited immediate degradation, with kinetics adhering to a first-order model. The comprehensive biodegradation of DBP involves de-esterification, β-oxidation, trans-esterification, and aromatic ring cleavage. FZ201 contains numerous catabolic genes that potentially facilitate PAE biodegradation. The DBP metabolic pathway was reconstructed by genome annotation and intermediate identification. Streptomyces species have an open pangenome with substantial genome expansion events during the evolutionary process, enabling extensive genetic diversity and highly plastic genomes within the Streptomyces genus. FZ201 had a diverse array of highly expressed genes associated with the degradation of PAEs, potentially contributing significantly to its adaptive advantage and efficiency of PAE degradation. Thus, FZ201 is a promising candidate for remediating highly PAE-contaminated environments. These findings enhance our preliminary understanding of the molecular mechanisms employed by Streptomyces for the removal of PAEs.

Analysis of codon usage bias of exonuclease genes in invertebrate iridescent viruses

Invertebrate iridescent viruses (IIVs) are double-stranded DNA viruses that belong to the Iridoviridae family. IIVs result diseases that vary in severity from subclinical to lethal in invertebrate hosts. Codon usage bias (CUB) analysis is a versatile method for comprehending the genetic and evolutionary aspects of species. In this study, we analyzed the CUB in 10 invertebrate iridescent viruses exonuclease genes by calculating and comparing the nucleotide contents, effective number of codons (ENC), codon adaptation index (CAI), relative synonymous codon usage (RSCU), and others. The results revealed that IIVs exonuclease genes are rich in A/T. The ENC analysis displayed a low codon usage bias in IIVs exonuclease genes. ENC-plot, neutrality plot, and parity rule 2 plot demonstrated that besides mutational pressure, other factors like natural selection, dinucleotide content, and aromaticity also contributed to CUB. The findings could enhance our understanding of the evolution of IIVs exonuclease genes.

Genetic background of adaptation of Crimean-Congo haemorrhagic fever virus to the different tick hosts

Crimean-Congo haemorrhagic fever orthonairovirus (CCHFV) is a negative-sense, single-stranded RNA virus with a segmented genome and the causative agent of a severe Crimean-Congo haemorrhagic fever (CCHF) disease. The virus is transmitted mainly by tick species in Hyalomma genus but other ticks such as representatives of genera Dermacentor and Rhipicephalus may also be involved in virus life cycle. To improve our understanding of CCHFV adaptation to its tick species, we compared nucleotide composition and codon usage patterns among the all CCHFV strains i) which sequences and other metadata as locality of collection and date of isolation are available in GenBank and ii) which were isolated from in-field collected tick species. These criteria fulfilled 70 sequences (24 coding for S, 23 for M, and 23 for L segment) of virus isolates originating from different representatives of Hyalomma and Rhipicephalus genera. Phylogenetic analyses confirmed that Hyalomma- and Rhipicephalus-originating CCHFV isolates belong to phylogenetically distinct CCHFV clades. Analyses of nucleotide composition among the Hyalomma- and Rhipicephalus-originating CCHFV isolates also showed significant differences, mainly in nucleotides located at the 3rd codon positions indicating changes in codon usage among these lineages. Analyses of codon adaptation index (CAI), effective number of codons (ENC), and other codon usage statistics revealed significant differences between Hyalomma- and Rhipicephalus-isolated CCHFV strains. Despite both sets of strains displayed a higher adaptation to use codons that are preferred by Hyalomma ticks than Rhipicephalus ticks, there were distinct codon usage preferences observed between the two tick species. These findings suggest that over the course of its long co-evolution with tick vectors, CCHFV has optimized its codon usage to efficiently utilize translational resources of Hyalomma species.

Comprehensive Genomics Investigation of Neboviruses Reveals Distinct Codon Usage Patterns and Host Specificity

Neboviruses (NeVs) from the Caliciviridae family have been linked to enteric diseases in bovines and have been detected worldwide. As viruses rely entirely on the cellular machinery of the host for replication, their ability to thrive in a specific host is greatly impacted by the specific codon usage preferences. Here, we systematically analyzed the codon usage bias in NeVs to explore the genetic and evolutionary patterns. Relative Synonymous Codon Usage and Effective Number of Codon analyses indicated a marginally lower codon usage bias in NeVs, predominantly influenced by the nucleotide compositional constraints. Nonetheless, NeVs showed a higher codon usage bias for codons containing G/C at the third codon position. The neutrality plot analysis revealed natural selection as the primary factor that shaped the codon usage bias in both the VP1 (82%) and VP2 (57%) genes of NeVs. Furthermore, the NeVs showed a highly comparable codon usage pattern to bovines, as reflected through Codon Adaptation Index and Relative Codon Deoptimization Index analyses. Notably, yak NeVs showed considerably different nucleotide compositional constraints and mutational pressure compared to bovine NeVs, which appear to be predominantly host-driven. This study sheds light on the genetic mechanism driving NeVs' adaptability, evolution, and fitness to their host species.

Schlafen14 Impairs HIV-1 Expression in a Codon Usage-Dependent Manner

Schlafen (SLFN) is a family of proteins upregulated by type I interferons with a regulatory role in translation. Intriguingly, SLFN14 associates with the ribosome and can degrade rRNA, tRNA, and mRNA in vitro, but a role in translation is still unknown. Ribosomes are important regulatory hubs during translation elongation of mRNAs rich in rare codons. Therefore, we evaluated the potential role of SLFN14 in the expression of mRNAs enriched in rare codons, using HIV-1 genes as a model. We found that, in a variety of cell types, including primary immune cells, SLFN14 regulates the expression of HIV-1 and non-viral genes based on their codon adaptation index, a measurement of the synonymous codon usage bias; consequently, SLFN14 inhibits the replication of HIV-1. The potent inhibitory effect of SLFN14 on the expression of the rare codon-rich transcript HIV-1 Gag was minimized by codon optimization. Mechanistically, we found that the endoribonuclease activity of SLFN14 is required, and that ribosomal RNA degradation is involved. Therefore, we propose that SLFN14 impairs the expression of HIV-1 transcripts rich in rare codons, in a catalytic-dependent manner.

Compositional Features and Codon Usage Pattern of Genes Associated with Parkinson's Disease

Codon usage bias (CUB) is the phenomenon of non-uniform usage of synonymous codons in which some codons are more used than others and it helps in understanding the molecular organization of genome. Bioinformatic approach was used to analyze the protein-coding sequences of genes associated with Parkinson's disease (PD) to explore compositional features and codon usage pattern as no details work was reported yet. The average improved effective number of codons (Nc) and Nc prime were 42.74 and 44.26 respectively, indicated that CUB was low in these genes. In most of the genes, the overall GC content was almost 50% and GC content at the 1st codon position was the highest while GC content at the 2nd codon position was lowest. Relative synonymous codon usage (RSCU) analysis elucidated over-represented (p > 1.6) and under-represented codons (p < 0.6). The GTG (Val) is the only codon over-represented in all genes. Over-represented codons except (GTG) were A or T ending while under-represented codons (except ACT) were G or C ending. The codons namely TTA (Leu), CTA (Leu), ATC (Ile), ATA (Ile), AGT (Ser), AAC (Asn), TGT (Cys), TGC (Cys), CGC (Arg), AGA (Arg), and AGG (Arg) were absent in SNCA1 to SNCA8 genes. The codon TCG (Ser) was absent in all genes except UCHL1 and PINK1. Correspondence analysis (COA) revealed that the pattern of codon usage differs among genes associated with PD. Neutrality plot analysis indicated some of the points are diagonal distribution suggested that mutation pressure influenced the CUB in genes associated with PD.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

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.

Synthetic biology approach revealed enhancement in haeme oxygenase-1 gene expression by codon pair optimization while reduction by codon deoptimization

Haem oxygenase-1 (HO-1) is a ubiquitously expressed gene involved in cellular homoeostasis, and its imbalance in expression results in various disorders. To alleviate such disorders, HO-1 gene expression needs to be modulated. Codon usage bias results from evolutionary forces acting on any nucleotide sequence and determines the gene expression. Like codon usage bias, codon pair bias also exists, playing a role in gene expression. In the present study, HO-1 gene was recoded by manipulating codon and codon pair bias, and four such constructs were made through codon/codon pair deoptimization and codon/codon pair optimization to reduce and enhance the HO-1 gene expression. Codon usage analysis was done for these constructs for four tissues brain, heart, pancreas and liver. Based on codon usage in different tissues, gene expression of these tissues was determined in terms of the codon adaptation index. Based on the codon adaptation index, minimum free energy, and translation efficiency, constructs were evaluated for enhanced or decreased HO-1 expression. The analysis revealed that for enhancing gene expression, codon pair optimization, while for reducing gene expression, codon deoptimization is efficacious. The recoded constructs developed in the study could be used in gene therapy regimens to cure HO-1 over or underexpression-associated disorders.

Relative synonymous codon usage and codon pair analysis of depression associated genes

Depression negatively impacts mood, behavior, and mental and physical health. It is the third leading cause of suicides worldwide and leads to decreased quality of life. We examined 18 genes available at the genetic testing registry (GTR) from the National Center for Biotechnological Information to investigate molecular patterns present in depression-associated genes. Different genotypes and differential expression of the genes are responsible for ensuing depression. The present study, investigated codon pattern analysis, which might play imperative roles in modulating gene expression of depression-associated genes. Of the 18 genes, seven and two genes tended to up- and down-regulate, respectively, and, for the remaining genes, different genotypes, an outcome of SNPs were responsible alone or in combination with differential expression for different conditions associated with depression. Codon context analysis revealed the abundance of identical GTG-GTG and CTG-CTG pairs, and the rarity of methionine-initiated codon pairs. Information based on codon usage, preferred codons, rare, and codon context might be used in constructing a deliverable synthetic construct to correct the gene expression level of the human body, which is altered in the depressive state. Other molecular signatures also revealed the role of evolutionary forces in shaping codon usage.

The genome of Salmacisia buchloëana, the parasitic puppet master pulling strings of sexual phenotypic monstrosities in buffalograss

To complete its parasitic lifecycle, Salmacisia buchloëana, a biotrophic fungus, manipulates reproductive organ development, meristem determinacy, and resource allocation in its dioecious plant host, buffalograss (Bouteloua dactyloides; Poaceae). To gain insight into S. buchloëana's ability to manipulate its host, we sequenced and assembled the 20.1 Mb genome of S. buchloëana into 22 chromosome-level pseudomolecules. Phylogenetic analysis suggests that S. buchloëana is nested within the genus Tilletia and diverged from Tilletia caries and Tilletia walkeri ∼40 MYA. We find that S. buchloëana contains a novel chromosome arm with no syntenic relationship to other publicly available Tilletia genomes, and that genes on the novel arm are upregulated upon infection, suggesting that this unique chromosomal segment may have played a critical role in S. buchloëana's evolution and host specificity. Salmacisia buchloëana has one of the largest fractions of serine peptidases (1.53% of the proteome) and one of the highest GC contents (62.3%) in all classified fungi. Analysis of codon base composition indicated that GC content is controlled more by selective constraints than directional mutation, and that S. buchloëana has a unique bias for the serine codon UCG. Finally, we identify 3 inteins within the S. buchloëana genome, 2 of which are located in a gene often used in fungal taxonomy. The genomic and transcriptomic resources generated here will aid plant pathologists and breeders by providing insight into the extracellular components contributing to sex determination in dioecious grasses.

Tobacco Plant: A Novel and Promising Heterologous Bioreactor for the Production of Recombinant Bovine Chymosin

The natural source of chymosin, a key enzyme in the dairy industry, is insufficient for rapidly growing cheese industries. Large-scale production of recombinant proteins in heterologous hosts provides an efficient alternative solution. Here, the codon-optimized synthetic prochymosin gene, which has a CAI index of 0.926, was subcloned from a cloning vector (pUC57-bCYM) into the pBI121 vector, resulting in the construct named pBI121-bCYM. CAI ranges from 0 to 1 and higher CAI improves gene expression in heterologous hosts. The overexpression of the prochymosin gene was under the control of constitutive CaMV 35S promoter and NOS terminator and was transferred into the tobacco via A. tumefaciens strain LBA4404. Explant type, regeneration method, inoculation temperature, cell density (OD600) of Agrobacterium for inoculation, and acetosyringone concentration were leaf explants, direct somatic embryogenesis, 19 °C, 0.1, and 100 µM, respectively. The successful integration and expression of the prochymosin gene, along with the bioactivity of recombinant chymosin, were confirmed by PCR, RT-PCR, and milk coagulation assay, respectively. Overall, this study reports the first successful overexpression of the codon-optimized prochymosin form of the bovine chymosin enzyme in the tobacco via indirect transformation. Production of recombinant bovine chymosin in plants can be an easy-to-scale-up, safe, and inexpensive platform.

Improved prediction of bacterial CRISPRi guide efficiency from depletion screens through mixed-effect machine learning and data integration

CRISPR interference (CRISPRi) is the leading technique to silence gene expression in bacteria; however, design rules remain poorly defined. We develop a best-in-class prediction algorithm for guide silencing efficiency by systematically investigating factors influencing guide depletion in genome-wide essentiality screens, with the surprising discovery that gene-specific features substantially impact prediction. We develop a mixed-effect random forest regression model that provides better estimates of guide efficiency. We further apply methods from explainable AI to extract interpretable design rules from the model. This study provides a blueprint for predictive models for CRISPR technologies where only indirect measurements of guide activity are available.

Circular replication-associated protein-encoding single-stranded DNA virus with risk of spillover is widely prevalent in domestic animals in China

Circular replication-associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses are highly diverse and have a broad range of hosts. In this study, we report the detection of Bo-Circo-like virus AH20-1 in the feces of diarrheal cattle. The virus has a circular genome of 3,912 nucleotides, three major putative open reading frames, and encodes a Rep gene of 310 amino acids. We found that the virus is closely related to the Bo-Circo-like virus CH strain, which belongs to the novel Kirkoviridae family. Furthermore, we conducted a nationwide surveillance program and found that the virus is prevalent in China (23.6%, 205/868), with the BCLa subtype being the predominant strain. Our findings suggest that the virus can infect sheep, highlighting the potential for cross-species transmission. Our pressure analysis indicates that the CRESS-DNA Kirkoviridae family has broad host adaptation, and that selection pressure played an important role in the evolution of its Rep genes. Our study underscores the need for continued epidemiological surveillance of this virus due to its widespread prevalence in our ruminant population and potential for cross-species transmission.

Codon Pattern and Context Analysis in Genes Triggering Alzheimer's Disease and Latent Tau Protein Aggregation Post-Anesthesia Exhibited Unique Molecular Patterns Associated with Functional Aspects

Background

Previous reports have demonstrated post-operative dementia and Alzheimer's disease (AD), and increased amyloid-β levels and tau hyperphosphorylation have been observed in animal models post-anesthesia.

Objective

After surgical interventions, loss in memory has been observed that has been found linked with genes modulated after anesthesia. Present study aimed to study molecular pattern present in genes modulated post anesthesia and involved in characters progressing towards AD.

Methods

In the present study, 17 transcript variants belonging to eight genes, which have been found to modulate post-anesthesia and contribute to AD progression, were envisaged for their compositional features, molecular patterns, and codon and codon context-associated studies.

Results

The sequences' composition was G/C rich, influencing dinucleotide preference, codon preference, codon usage, and codon context. The G/C nucleotides being highly occurring nucleotides, CpGdinucleotides were also preferred; however, CpG was highly disfavored at p3-1 at the codon junction. The nucleotide composition of Cytosine exhibited a unique feature, and unlike other nucleotides, it did not correlate with codon bias. Contrarily, it correlated with the sequence lengths. The sequences were leucine-rich, and multiple leucine repeats were present, exhibiting the functional role of neuroprotection from neuroinflammation post-anesthesia.

Conclusions

The analysis pave the way to elucidate unique molecular patterns in genes modulated during anesthetic treatment and might help ameliorate the ill effects of anesthetics in the future.

CodonU: A Python Package for Codon Usage Analysis

Codon Usage Analysis (CUA) has been accompanied by several web servers and independent programs written in several programming languages. Also this diversity speaks for the need of a reusable software that can be helpful in reading, manipulating and acting as a pipeline for such data and file formats. This kind of analyses use multiple tools to address the multifaceted aspects of CUA. So, we propose CodonU, a package written in Python language to integrate all aspects. It is compatible with existing file formats and can be used solely or with a group of other such packages. The proposed package incorporates various statistical measures necessary for codon usage analysis. The measures vary with nature of the sequences, viz. for nucleotide, codon adaptation index (CAI), codon bias index (CBI), tRNA adaptation index (tAI) etc. and for protein sequences Gravy score etc. Users can also perform the correspondence analysis (COA). This package also provides the liberty to generate graphics to users, and also develop phylogenetic tree. Capabilities of the proposed package were checked thoroughly on a genomic set of Staphylococcus aureus.

In-silico Codon Context and Synonymous Usage Analysis of Genes for Molecular Mechanisms Inducing Autophagy and Apoptosis with Reference to Neurodegenerative Disorders

Background

Autophagy and apoptosis are cellular processes that maintain cellular homeostasis and remove damaged or aged organelles or aggregated and misfolded proteins. Stress factors initiate the signaling pathways common to autophagy and apoptosis. An imbalance in the autophagy and apoptosis, led by cascade of molecular mechanism prior to both processes culminate into neurodegeneration.

Objective

In present study, we urge to investigate the codon usage pattern of genes which are common before initiating autophagy and apoptosis.

Methods

In the present study, we took up eleven genes (DAPK1, BECN1, PIK3C3 (VPS34), BCL2, MAPK8, BNIP3 L (NIX), PMAIP1, BAD, BID, BBC3, MCL1) that are part of molecular signaling mechanism prior to autophagy and apoptosis. We analyzed dinucleotide odds ratio, codon bias, usage, context, and rare codon analysis.

Results

CpC and GpG dinucleotides were abundant, with the dominance of G/C ending codons as preferred codons. Clustering analysis revealed that MAPK8 had a distinct codon usage pattern compared to other envisaged genes. Both positive and negative contexts were observed, and GAG-GAG followed by CTG-GCC was the most abundant codon pair. Of the six synonymous arginine codons, two codons CGT and CGA were the rarest.

Conclusions

The information presented in the study may be used to manipulate the process of autophagy and apoptosis and to check the pathophysiology associated with their dysregulation.

Codon usage bias analysis of the spike protein of human coronavirus 229E and its host adaptability

Human coronavirus 229E (HCoV-229E) represents one of the known coronaviruses capable of infecting humans and causes mild respiratory symptoms. It is also considered to have a zoonotic source, originating from animals and being transmitted the humans. In this study, a comprehensive phylogenetic and codon usage analysis of the spike (S) gene of HCoV-229E was conducted. Utilizing phylogenetic analysis and principal component analysis, HCoV-229E was categorized into four distinct clusters, each demonstrating unique host affiliations. Furthermore, it was observed that the codon usage bias within the S gene of HCoV-229E is relatively low, primarily influenced by natural selection patterns, with contributions from mutation pressure and dinucleotide abundance. Comparative analysis involving Codon Adaptation Index (CAI) and Relative Codon Deoptimization Index (RCDI) revealed that the codon usage pattern of HCoV-229E mirrors more closely that of camels, as opposed to alpacas and humans. The elucidation of the codon usage pattern within HCoV-229E, which we have meticulously examined, offers valuable insights for a more comprehensive comprehension of viral features, history, and evolutionary trajectory.

Differential Synonymous Codon Selection in the B56 Gene Family of PP2A Regulatory Subunits

Protein phosphatase 2A (PP2A) functions as a tumor suppressor and consists of a scaffolding, catalytic, and regulatory subunit. The B56 gene family of regulatory subunits impart distinct functions onto PP2A. Codon usage bias (CUB) involves the selection of synonymous codons, which can affect gene expression by modulating processes such as transcription and translation. CUB can vary along the length of a gene, and differential use of synonymous codons can be important in the divergence of gene families. The N-termini of the gene product encoded by B56α possessed high CUB, high GC content at the third codon position (GC3), and high rare codon content. In addition, differential CUB was found in the sequence encoding two B56γ N-terminal splice forms. The sequence encoding the N-termini of B56γ/γ, relative to B56δ/γ, displayed CUB, utilized more frequent codons, and had higher GC3 content. B56α mRNA had stronger than predicted secondary structure at their 5' end, and the B56δ/γ splice variants had long regions of weaker than predicted secondary structure at their 5' end. The data suggest that B56α is expressed at relatively low levels as compared to the other B56 isoforms and that the B56δ/γ splice variant is expressed more highly than B56γ/γ.

A small stretch of poor codon usage at the beginning of dengue virus open reading frame may act as a translational checkpoint

OBJECTIVE

Rare codons were previously shown to be enriched at the beginning of the dengue virus (DENV) open reading frame. However, the role of rare codons in regulating translation efficiency and replication of DENV remains unclear. The present study aims to clarify the significance of rare codon usage at the beginning of DENV transcripts using the codon adaptation index (CAI).

METHODOLOGY

CAIs of the whole starting regions of DENV transcripts as well as 18-codon sliding windows of the regions were analyzed.

RESULTS

One of the intriguing findings is that those rare codons do not typically result in uniformly low CAI in the starting region with rare codons. However, it shows a notable local drop in CAI around the 50th codon in all dengue serotypes. This suggests that there may be a translational checkpoint at this site and that the rare codon usage upstream to this checkpoint may not be related to translational control.

Codon optimization of a gene encoding DNA polymerase from Pyrococcus furiosus and its expression in Escherichia coli

BACKGROUND

DNA polymerase is an essential component in PCR assay for DNA synthesis. Improving DNA polymerase with characteristics indispensable for a powerful assay is crucial because it can be used in wide-range applications. Derived from Pyrococcus furiosus, Pfu DNA polymerase (Pfu pol) is one of the excellent polymerases due to its high fidelity. Therefore, we aimed to develop Pfu pol from a synthetic gene with codon optimization to increase its protein yield in Escherichia coli.

RESULTS

Recombinant Pfu pol was successfully expressed and purified with a two-step purification process using nickel affinity chromatography, followed by anion exchange chromatography. Subsequently, the purified Pfu pol was confirmed by Western blot analysis, resulting in a molecular weight of approximately 90 kDa. In the final purification process, we successfully obtained a large amount of purified enzyme (26.8 mg/L). Furthermore, the purified Pfu pol showed its functionality and efficiency when tested for DNA amplification using the standard PCR.

CONCLUSIONS

Overall, a high-level expression of recombinant Pfu pol was achieved by employing our approach in the present study. In the future, our findings will be useful for studies on synthesizing recombinant DNA polymerase in E. coli expression system.

Base composition, adaptation, and evolution of goose astroviruses: codon-based investigation

Goose astroviruses (GoAstVs) are causative agents that account for fatal infection of goslings characterized by visceral urate deposition, resulting in severe economic losses in major goose-producing regions in China since 2017. In this study, we sought to unravel the intrinsic properties associated with adaptation and evolution in the host environment of GoAstVs. Consistent results from phylogenetic analysis and correspondence analysis performed on the codon usage patterns (CUPs) reveal 2 clusters of GoAstVs, namely, GoAstV-1 and GoAstV-2. However, multiple similar compositional characteristics were found, despite the high divergence between GoAstV-1 and GoAstV-2. Studies on the base composition of GoAstVs reveal an A/U bias, indicating a compositional constraint, while natural selection prevailed in determining the CUPs in the virus genome based on our neutrality plot analysis, reflecting high adaptive pressure to fit the host environment. Codon adaptation index (CAI) analysis revealed a higher degree of fitness to the CUPs of the corresponding host for GoAstVs than avian influenza virus and betacoronaviruses, which may be a favorable factor contributing to the high pathogenicity and wide distribution of GoAstVs in goslings. In addition, GoAstVs were less adapted to ducks and chickens, with significantly lower CAI values than to geese, which may be a reason for the different prevalence of GoAstVs among these species. Extensive investigations on dinucleotide distribution revealed a significant suppression of the CpG and UpA motifs in the virus genome, which may facilitate adaptation to the host's innate immune system by evading surveillance. In addition, our study reported the trends of increasing fitness to the host's microenvironment for GoAstVs through increasing adaptation to host CUPs and ongoing reduction of CpG motifs in the virus genome. The present analysis deepens our understanding of the basic biology, pathogenesis, adaptation and evolutionary pattern of GoAstVs, and contributes to the development of novel antiviral strategies.

Codon usage bias analysis of the gene encoding NAD+-dependent DNA ligase protein of Invertebrate iridescent virus 6

The genome of Invertebrate iridescent virus 6 (IIV6) contains a sequence that shows similarity to eubacterial NAD+-dependent DNA ligases. The 615-amino acid open reading frame (ORF 205R) consists of several domains, including an N-terminal domain Ia, followed by an adenylation domain, an OB-fold domain, a helix-hairpin-helix (HhH) domain, and a BRCT domain. Notably, the zinc finger domain, typically present in NAD+-dependent DNA ligases, is absent in ORF 205R. Since the protein encoded by ORF 205R (IIV6 DNA ligase gene) is involved in critical functions such as DNA replication, modification, and repair, it is crucial to comprehend the codon usage associated with this gene. In this paper, the codon usage bias (CUB) in DNA ligase gene of IIV6 and 11 reference iridoviruses was analyzed by comparing the nucleotide contents, relative synonymous codon usage (RSCU), effective number of codons (ENC), codon adaptation index (CAI), relative abundance of dinucleotides and other indices. Both the base content and the RCSU analysis indicated that the A- and T-ending codons were mostly favored in the DNA ligase gene of IIV6. The ENC value of 35.64 implied a high CUB in the IIV6 DNA ligase gene. The ENC plot, neutrality plot, parity rule 2 plot, correspondence analysis revealed that mutation pressure and natural selection had an impact on the CUB of the IIVs DNA ligase genes. Additionally, the analysis of codon adaptation index demonstrated that the IIV6 DNA ligase gene is strongly adapted to its host. These findings will improve our comprehension of the CUB of IIV6 DNA ligase and reference genes, which may provide the required information for a fundamental evolutionary analysis of these genes.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

Mutation and codon bias analysis of the spike protein of Omicron, the recent variant of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant is a heavily mutated virus and designated as a variant of concern. To investigate the codon usage pattern of this new variant, we performed mutation and codon bias analysis for Omicron as well as for its sub-lineages BA.1 and BA.2 and compared them with the original SARS-CoV-2 and the Delta variant sequences obtained in this study. Our results indicate that the sub-lineage BA.1 and BA.2 have up to 23 sites of difference on the spike protein, which have minimal impact on function. The Omicron variant and its sub-lineages have similar codon usage patterns and A/U ending codons appear to be preferred over G/C ending codons. The Omicron has a lower degree of codon usage bias in spite of evidence that natural selection, mutation pressure and dinucleotide abundance shape the codon usage bias of Omicron, with natural selection being more significant on BA.2 than the other sub-lineages of Omicron. The codon usage pattern of Omicron variant that we explored provides valid information for a clearer understanding of Omicron and its sub-lineages, which could find application in vaccine development and optimization.

Analysis of codon usage patterns of porcine enteric alphacoronavirus and its host adaptability

Porcine enteric alphacoronavirus (PEAV) is a newly emerging swine enteropathogen that poses a threat to the swine industry. To understand the PEAV genome evolution, we performed a comprehensive analysis of the codon usage patterns in fifty-nine PEAV strains currently available. Phylogenetic analysis showed that PEAV can be divided into six lineages. Effective number of codons analysis demonstrated that the PEAV genome exhibits a low codon usage bias (CUB). Nucleotide composition analysis indicated that the PEAV genome has the most abundant nucleotide U content, with GC content (39.37% ± 0.08%) much lower than AU content (60.63% ± 0.08%). Neutrality and effective number of codons plot analyses suggested that natural selection rather than mutation pressure dominates the CUB of PEAV. Host adaptation analysis revealed that PEAV fits the codon usage pattern of non-human primates, humans and mice better than that of pigs. Our data enriches information on PEAV evolution, host adaptability, and cross-species transmission.

Non-adaptive evolution in codon usage of human-origin monkeypox virus

Monkeypox virus (Mpox) is a zoonotic infectious disease that threatens human and animal health, with a global outbreak of the low-pathogenic Mpox beginning from 2022. In this study, we analyzed the codon usage of Mpox between two clades, Clade-I and Clade-IIb-B, to understand changes in host adaptation. Clade-IIb-B of the Mpox genome underwent non-adaptive evolution making it less adapted to its host than Clade-I. The analysis of individual genes revealed that 48 genes exhibited non-adaptive mutation, while 38 genes underwent adaptive mutations. Genes involved in replication, transcription, and host-modulation exhibited a mix of adaptive and non-adaptive evolutionary patterns. This study also found that the mutations of Mpox led to changes in non-adaptative genes in different organs. Additionally, we found that codon usage of Mpox was less similar to that of up-regulated host genes and more similar to that of down-regulated host genes post-infection, indicating that codon usage affects host gene expression. Overall, the study highlights the non-adaptive changes in codon usage as a potential cause of differences in Mpox virulence and provides insights into the evolutionary and adaptive mechanisms of Mpox and its potential impact on pathogenicity and host adaptation.

Human Schlafen 11 inhibits influenza A virus production

Schlafen (SLFN) proteins are a subset of interferon-stimulated early response genes with antiviral properties. An antiviral mechanism of SLFN11 was previously demonstrated in human immunodeficiency virus type 1 (HIV-1)-infected cells, and it was shown that SLFN11 inhibited HIV-1 virus production in a codon usage-specific manner. The codon usage patterns of many viruses are vastly different from those of their hosts. The codon usage-specific inhibition of HIV-1 expression by SLFN11 suggests that SLFN11 may be able to inhibit other viruses with a suboptimal codon usage pattern. However, the effect of SLFN11 on the replication of influenza A virus (IAV) has never been reported. The induction of SLFN11 expression was observed upon IAV infection. The reduction of SLFN11 expression also promotes influenza virus replication. Moreover, we found that overexpression of SLFN11 could reduce the expression of a reporter gene with a viral codon usage pattern, and the inhibition of viral hemagglutinin (HA) gene was codon-specific as the expression of codon optimized HA was not affected. These results indicate that SLFN11 inhibits the influenza A virus in a codon-specific manner and that SLFN11 may contribute to innate defense against influenza A viruses.

Genetic Engineering of Acidithiobacillus ferridurans Using CRISPR Systems To Mitigate Toxic Release in Biomining

Biomining processes utilize microorganisms, such as Acidithiobacillus, to extract valuable metals by producing sulfuric acid and ferric ions that dissolve sulfidic minerals. However, excessive production of these compounds can result in metal structure corrosion and groundwater contamination. Synthetic biology offers a promising solution to improve Acidithiobacillus strains for sustainable, eco-friendly, and cost-effective biomining, but genetic engineering of these slow-growing microorganisms is challenging with current inefficient and time-consuming methods. To address this, we established a CRISPR-dCas9 system for gene knockdown in A. ferridurans JAGS, successfully downregulating the transcriptional levels of two genes involved in sulfur oxidation. More importantly, we constructed an all-in-one CRISPR-Cas9 system for fast and efficient genome editing in A. ferridurans JAGS, achieving seamless gene deletion (HdrB3), promoter substitution (Prus to Ptac), and exogenous gene insertion (GFP). Additionally, we created a HdrB-Rus double-edited strain and performed biomining experiments to extract Ni from pyrrhotite tailings. The engineered strain demonstrated a similar Ni recovery rate to wild-type A. ferridurans JAGS but with significantly lower production of iron ions and sulfuric acid in leachate. These high-efficient CRISPR systems provide a powerful tool for studying gene functions and creating useful recombinants for synthetic biology-assisted biomining applications in the future.

Phylogenetic Analysis and Codon Usage Bias Reveal the Base of Feline and Canine Chaphamaparvovirus for Cross-Species Transmission

Chaphamaparvoviruses (ChPVs) are ancient viruses that have been detected in a variety of hosts. In this study, through a phylogenetic analysis and the adaptability of ChPV to multiple hosts, we evaluated the basis for the ability of feline (FeChPV) and canine ChPV (CaChPV) for cross-species transmission. Phylogenetic analysis showed that FeChPV and CaChPV were closely related. Notably, two strains of ChPVs isolated from domestic cats and two from dogs clustered together with CaChPVs and FeChPVs, respectively, suggesting that the stringent boundaries between canine and feline ChPV may be broken. Further analysis revealed that CaChPV and FeChPV were more adapted to dogs than to cats. Mutation analysis identified several shared mutations in cross-species-transmissible strains. Furthermore, the VP structures of FeChPV and CaChPV exhibited a high degree of similarity across both cross-species-transmissible and non-cross-species-transmissible strains. However, it is crucial to note that these results are largely computational, and limitations exist in terms of the number and diversity of samples analyzed; the capacity for cross-species transmission should be approached with caution and elucidated in further studies.

Orb2 enables rare-codon-enriched mRNA expression during Drosophila neuron differentiation

Regulation of codon optimality is an increasingly appreciated layer of cell- and tissue-specific protein expression control. Here, we use codon-modified reporters to show that differentiation of Drosophila neural stem cells into neurons enables protein expression from rare-codon-enriched genes. From a candidate screen, we identify the cytoplasmic polyadenylation element binding (CPEB) protein Orb2 as a positive regulator of rare-codon-dependent expression in neurons. Using RNA sequencing, we reveal that Orb2-upregulated mRNAs in the brain with abundant Orb2 binding sites have a rare-codon bias. From these Orb2-regulated mRNAs, we demonstrate that rare-codon enrichment is important for expression control and social behavior function of the metabotropic glutamate receptor (mGluR). Our findings reveal a molecular mechanism by which neural stem cell differentiation shifts genetic code regulation to enable critical mRNA and protein expression.

A bacterial-like Pictet-Spenglerase drives the evolution of fungi to produce β-carboline glycosides together with separate genes

Diverse β-carboline (βC) alkaloids are produced by microbes, plants, and animals with myriad bioactivities and drug potentials. However, the biosynthetic mechanism of βCs remains largely elusive, especially regarding the hydroxyl and glucosyl modifications of βCs. Here, we report the presence of the bacterial-like Pictet-Spenglerase gene Fcs1 in the entomopathogenic Beauveria fungi that can catalyze the biosynthesis of the βC skeleton. The overexpression of Fcs1 in Beauveria bassiana led to the identification of six βC methyl glycosides, termed bassicarbosides (BCSs) A-F. We verified that the cytochrome P450 (CYP) genes adjacent to Fcs1 cannot oxidize βCs. Alternatively, the separated CYP684B2 family gene Fcs2 was identified to catalyze βC hydroxylation together with its cofactor gene Fcs3. The functional homologue of Fcs2 is only present in the Fcs1-containing fungi and highly similar to the Fcs1-connected yet nonfunctional CYP. Both evolved quicker than those from fungi without Fcs1 homologues. Finally, the paired methyl/glucosyl transferase genes were verified to mediate the production of BCSs from hydroxy-βCs. All these functionally verified genes are located on different chromosomes of Beauveria, which is in contrast to the typical content-clustered feature of fungal biosynthetic gene clusters (BGCs). We also found that the production of BCSs selectively contributed to fungal infection of different insect species. Our findings shed light on the biosynthetic mechanism of βC glycosides, including the identification of a βC hydroxylase. The results of this study also propose an evolving process of fungal BGC formation following the horizontal transfer of a bacterial gene to fungi.

Structural and biochemical characterization of Schlafen11 N-terminal domain

Schlafen11 (SLFN11) is one of the most studied Schlafen proteins that plays vital roles in cancer therapy and virus-host interactions. Herein, we determined the crystal structure of the Sus scrofa SLFN11 N-terminal domain (NTD) to 2.69 Å resolution. sSLFN11-NTD is a pincer-shaped molecule that shares an overall fold with other SLFN-NTDs but exhibits distinct biochemical characteristics. sSLFN11-NTD is a potent RNase cleaving type I and II tRNAs and rRNAs, and with preference to type II tRNAs. Consistent with the codon usage-based translation suppression activity of SLFN11, sSLFN11-NTD cleaves synonymous serine and leucine tRNAs with different efficiencies in vitro. Mutational analysis revealed key determinates of sSLFN11-NTD nucleolytic activity, including the Connection-loop, active site, and key residues essential for substrate recognition, among which E42 constrains sSLFN11-NTD RNase activity, and all nonconservative mutations of E42 stimulated RNase activities. sSLFN11 inhibited the translation of proteins with a low codon adaptation index in cells, which mainly dependent on the RNase activity of the NTD because E42A enhanced the inhibitory effect, but E209A abolished inhibition. Our findings provide structural characterization of an important SLFN11 protein and expand our understanding of the Schlafen family.

Genome-resolved metagenomics provides insights into the ecological roles of the keystone taxa in heavy-metal-contaminated soils

Microorganisms that exhibit resistance to environmental stressors, particularly heavy metals, have the potential to be used in bioremediation strategies. This study aimed to explore and identify microorganisms that are resistant to heavy metals in soil environments as potential candidates for bioremediation. Metagenomic analysis was conducted using microbiome metagenomes obtained from the rhizosphere of soil contaminated with heavy metals and mineral-affected soil. The analysis resulted in the recovery of a total of 175 metagenome-assembled genomes (MAGs), 73 of which were potentially representing novel taxonomic levels beyond the genus level. The constructed ecological network revealed the presence of keystone taxa, including Rhizobiaceae, Xanthobacteraceae, Burkholderiaceae, and Actinomycetia. Among the recovered MAGs, 50 were associated with these keystone taxa. Notably, these MAGs displayed an abundance of genes conferring resistance to heavy metals and other abiotic stresses, particularly those affiliated with the keystone taxa. These genes were found to combat excessive accumulation of zinc/manganese, arsenate/arsenite, chromate, nickel/cobalt, copper, and tellurite. Furthermore, the keystone taxa were found to utilize both organic and inorganic energy sources, such as sulfur, arsenic, and carbon dioxide. Additionally, these keystone taxa exhibited the ability to promote vegetation development in re-vegetated mining areas through phosphorus solubilization and metabolite secretion. In summary, our study highlights the metabolic adaptability and ecological significance of microbial keystone taxa in mineral-affected soils. The MAGs associated with keystone taxa exhibited a markedly higher number of genes related to abiotic stress resistance and plant growth promotion compared to non-keystone taxa MAGs.

Identification and characterization of codon usage pattern and influencing factors in HFRS-causing hantaviruses

Hemorrhagic fever with renal syndrome (HFRS) is an acute viral zoonosis carried and transmitted by infected rodents through urine, droppings, or saliva. The etiology of HFRS is complex due to the involvement of viral factors and host immune and genetic factors which hinder the development of potential therapeutic solutions for HFRS. Hantaan virus (HTNV), Dobrava-Belgrade virus (DOBV), Seoul virus (SEOV), and Puumala virus (PUUV) are predominantly found in hantaviral species that cause HFRS in patients. Despite ongoing prevention and control efforts, HFRS remains a serious economic burden worldwide. Furthermore, recent studies reported that the hantavirus nucleocapsid protein is a multi-functional protein and plays a major role in the replication cycle of the hantavirus. However, the precise mechanism of the nucleoproteins in viral pathogenesis is not completely understood. In the framework of the current study, various in silico approaches were employed to identify the factors influencing the codon usage pattern of hantaviral nucleoproteins. Based on the relative synonymous codon usage (RSCU) values, a comparative analysis was performed between HFRS-causing hantavirus and their hosts, suggesting that HTNV, DOBV, SEOV, and PUUV, were inclined to evolve their codon usage patterns that were comparable to those of their hosts. The results indicated that most of the overrepresented codons had AU-endings, which revealed that mutational pressure is the major force shaping codon usage patterns. However, the influence of natural selection and geographical factors cannot be ignored on viral codon usage bias. Further analysis also demonstrated that HFRS causing hantaviruses adapted host-specific codon usage patterns to sustain successful replication and transmission chains within hosts. To our knowledge, no study to date reported the factors influencing the codon usage pattern within hantaviral nucleoproteins. Thus, the proposed computational scheme can help in understanding the underlying mechanism of codon usage patterns in HFRS-causing hantaviruses which lend a helping hand in designing effective anti-HFRS treatments in future. This study, although comprehensive, relies on in silico methods and thus necessitates experimental validation for more solid outcomes. Beyond the identified factors influencing viral behavior, there could be other yet undiscovered influences. These potential factors should be targets for further research to improve HFRS therapeutic strategies.

Application of codon usage and context analysis in genes up- or down-regulated in neurodegeneration and cancer to combat comorbidities

Introduction

Neurodegeneration and cancer present in comorbidities with inverse effects due to the expression of genes and pathways acting in opposition. Identifying and studying the genes simultaneously up or downregulated during morbidities helps curb both ailments together.

Methods

This study examines four genes. Three of these (Amyloid Beta Precursor Protein (APP), Cyclin D1 (CCND1), and Cyclin E2 (CCNE2) are upregulated, and one protein phosphatase 2 phosphatase activator (PTPA) is simultaneously downregulated in both disorders. We investigated molecular patterns, codon usage, codon usage bias, nucleotide bias in the third codon position, preferred codons, preferred codon pairs, rare codons, and codon context.

Results

Parity analysis revealed that T is preferred over A, and G is preferred over C in the third codon position, suggesting composition plays no role in nucleotide bias in both the upregulated and downregulated gene sets and that mutational forces are stronger in upregulated gene sets than in downregulated ones. Transcript length influenced the overall %A composition and codon bias, and the codon AGG exerted the strongest influence on codon usage in both the upregulated and downregulated gene sets. Codons ending in G/C were preferred for 16 amino acids, and glutamic acid-, aspartic acid-, leucine-, valine-, and phenylalanine-initiated codon pairs were preferred in all genes. Codons CTA (Leu), GTA (Val), CAA (Gln), and CGT (Arg) were underrepresented in all examined genes.

Discussion

Using advanced gene editing tools such as CRISPR/Cas or any other gene augmentation technique, these recoded genes may be introduced into the human body to optimize gene expression levels to augment neurodegeneration and cancer therapeutic regimens simultaneously.

In silico and experimental methods for designing a potent anticancer arazyme-herceptin fusion protein in HER2-positive breast cancer

CONTEXT

Breast cancer is the most prevalent type of malignancies among women worldwide and is associated with serious physical and mental consequences. Current chemotherapies may lack successful outcomes; thus, the development of targeted recombinant immunotoxins is plausible. The predicted B cell and T cell epitopes of arazyme of the fusion protein are able to elicit immune response. The results of codon adaptation tool of herceptin-arazyme have improved from 0.4 to 1. The in silico immune simulation results showed significant response for immune cells. In conclusion, our findings show that the known multi-epitope fusion protein may activate humoral and cellular immune responses and maybe a possible candidate for breast cancer treatment.

METHODS

In this study, the selected monoclonal antibody constituting herceptin and the bacterial metalloprotease, arazyme, was used with different peptide linkers to design a novel fusion protein to predict different B cell and T cell epitopes by the means of the relevant databases. Modeler 10.1 and I-TASSER online server were used to predict and validate the 3D structure and then docked to HER2-receptor using HADDOCK2.4 web server. The molecular dynamics (MD) simulations of the arazyme-linker-herceptin-HER2 complex were performed by GROMACS 2019.6 software. The sequence of arazyme-herceptin was optimized for the expression in prokaryotic host using online servers and cloned into pET-28a plasmid. The recombinant pET28a was transferred into the Escherichia coli BL21DE3. Expression and binding affinity of arazyme-herceptin and arazyme to human breast cancer cell lines (SK-BR-3/HER2 + and MDA-MB-468/HER2 -) were validated by the SDS-PAGE and cell‑ELISA, respectively.