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.

Human parainfluenza virus 3 vaccine candidates attenuated by codon-pair deoptimization are immunogenic and protective in hamsters.

Human parainfluenza virus type 3 (HPIV3) is a major pediatric respiratory pathogen lacking available vaccines or antiviral drugs. We generated live-attenuated HPIV3 vaccine candidates by codon-pair deoptimization (CPD). HPIV3 open reading frames (ORFs) encoding the nucleoprotein (N), phosphoprotein (P), matrix (M), fusion (F), hemagglutinin-neuraminidase (HN), and polymerase (L) were modified singly or in combination to generate 12 viruses designated Min-N, Min-P, Min-M, Min-FHN, Min-L, Min-NP, Min-NPM, Min-NPL, Min-PM, Min-PFHN, Min-MFHN, and Min-PMFHN. CPD of N or L severely reduced growth in vitro and was not further evaluated. CPD of P or M was associated with increased and decreased interferon (IFN) response in vitro, respectively, but had little effect on virus replication. In Vero cells, CPD of F and HN delayed virus replication, but final titers were comparable to wild-type (wt) HPIV3. In human lung epithelial A549 cells, CPD F and HN induced a stronger IFN response, viral titers were reduced 100-fold, and the expression of F and HN proteins was significantly reduced without affecting N or P or the relative packaging of proteins into virions. Following intranasal infection in hamsters, replication in the nasal turbinates and lungs tended to be the most reduced for viruses bearing CPD F and HN, with maximum reductions of approximately 10-fold. Despite decreased in vivo replication (and lower expression of CPD F and HN in vitro), all viruses induced titers of serum HPIV3-neutralizing antibodies similar to wt and provided complete protection against HPIV3 challenge. In summary, CPD of HPIV3 yielded promising vaccine candidates suitable for further development.

Molecular Evolution of Protein Sequences and Codon Usage in Monkeypox Viruses.

The monkeypox virus (mpox virus, MPXV) epidemic in 2022 has posed a significant public health risk. Yet, the evolutionary principles of MPXV remain largely unknown. Here, we examined the evolutionary patterns of protein sequences and codon usage in MPXV. We first demonstrated the signal of positive selection in OPG027, specifically in the Clade I lineage of MPXV. Subsequently, we discovered accelerated protein sequence evolution over time in the variants responsible for the 2022 outbreak. Furthermore, we showed strong epistasis between amino acid substitutions located in different genes. The codon adaptation index (CAI) analysis revealed that MPXV genes tended to use more non-preferred codons compared to human genes, and the CAI decreased over time and diverged between clades, with Clade I > IIa and IIb-A > IIb-B. While the decrease in fatality rate among the three groups aligned with the CAI pattern, it remains unclear whether this correlation was coincidental or if the deoptimization of codon usage in MPXV led to a reduction in fatality rates. This study sheds new light on the mechanisms that govern the evolution of MPXV in human populations.

The systematic codon usage bias has an important effect on genetic adaption in native species.

Random mutations increase genetic variety and natural selection enhances adaption over generations. Codon usage biases (CUB) provide clues about the genome adaptation mechanisms of native species and extremophile species. Significant numbers of gene (CDS) of nine classes of endangered, native species, including extremophiles and mesophiles were utilised to compute CUB. Codon usage patterns differ among the lineages of endangered and extremophiles with native species. Polymorphic usage of nucleotides with codon burial suggests parallelism of native species within relatively confined taxonomic groups. Utilizing the deviation pattern of CUB of endangered and native species, I present a calculation parameter to estimate the extinction risk of endangered species. Species diversity and extinction risk are both positively associated with the propensity of random mutation in CDS (Coding DNA sequence). Codon bias tenet profoundly selected and it governs to adaptive evolution of native species.

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection.

BACKGROUND: Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures. RESULTS: Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection. CONCLUSIONS: Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

Profiling expression strategies for a type III polyketide synthase in a lysate-based, cell-free system.

Some of the most metabolically diverse species of bacteria (e.g., Actinobacteria) have higher GC content in their DNA, differ substantially in codon usage, and have distinct protein folding environments compared to tractable expression hosts like Escherichia coli. Consequentially, expressing biosynthetic gene clusters (BGCs) from these bacteria in E. coli often results in a myriad of unpredictable issues with regard to protein expression and folding, delaying the biochemical characterization of new natural products. Current strategies to achieve soluble, active expression of these enzymes in tractable hosts can be a lengthy trial-and-error process. Cell-free expression (CFE) has emerged as a valuable expression platform as a testbed for rapid prototyping expression parameters. Here, we use a type III polyketide synthase from Streptomyces griseus, RppA, which catalyzes the formation of the red pigment flaviolin, as a reporter to investigate BGC refactoring techniques. We applied a library of constructs with different combinations of promoters and rppA coding sequences to investigate the synergies between promoter and codon usage. Subsequently, we assess the utility of cell-free systems for prototyping these refactoring tactics prior to their implementation in cells. Overall, codon harmonization improves natural product synthesis more than traditional codon optimization across cell-free and cellular environments. More importantly, the choice of coding sequences and promoters impact protein expression synergistically, which should be considered for future efforts to use CFE for high-yield protein expression. The promoter strategy when applied to RppA was not completely correlated with that observed with GFP, indicating that different promoter strategies should be applied for different proteins. In vivo experiments suggest that there is correlation, but not complete alignment between expressing in cell free and in vivo. Refactoring promoters and/or coding sequences via CFE can be a valuable strategy to rapidly screen for catalytically functional production of enzymes from BCGs, which advances CFE as a tool for natural product research.

GNN Codon Adjacency Tunes Protein Translation.

The central dogma treats the ribosome as a molecular machine that reads one mRNA codon at a time as it adds each amino acid to its growing peptide chain. However, this and previous studies suggest that ribosomes actually perceive pairs of adjacent codons as they take three-nucleotide steps along the mRNA. We examined GNN codons, which we find are surprisingly overrepresented in eukaryote protein-coding open reading frames (ORFs), especially immediately after NNU codons. Ribosome profiling experiments in yeast revealed that ribosomes with NNU at their aminoacyl (A) site have particularly elevated densities when NNU is immediately followed (3') by a GNN codon, indicating slower mRNA threading of the NNU codon from the ribosome's A to peptidyl (P) sites. Moreover, if the assessment was limited to ribosomes that have only recently arrived at the next codon, by examining 21-nucleotide ribosome footprints (21-nt RFPs), elevated densities were observed for multiple codon classes when followed by GNN. This striking translation slowdown at adjacent 5'-NNN GNN codon pairs is likely mediated, in part, by the ribosome's CAR surface, which acts as an extension of the A-site tRNA anticodon during ribosome translocation and interacts through hydrogen bonding and pi stacking with the GNN codon. The functional consequences of 5'-NNN GNN codon adjacency are expected to influence the evolution of protein coding sequences.

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica.

MAIN CONCLUSION: In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

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.

mRNA context and translation factors determine decoding in alternative nuclear genetic codes.

The genetic code is a set of instructions that determine how the information in our genetic material is translated into amino acids. In general, it is universal for all organisms, from viruses and bacteria to humans. However, in the last few decades, exceptions to this rule have been identified both in pro- and eukaryotes. In this review, we discuss the 16 described alternative eukaryotic nuclear genetic codes and observe theories of their appearance in evolution. We consider possible molecular mechanisms that allow codon reassignment. Most reassignments in nuclear genetic codes are observed for stop codons. Moreover, in several organisms, stop codons can simultaneously encode amino acids and serve as termination signals. In this case, the meaning of the codon is determined by the additional factors besides the triplets. A comprehensive review of various non-standard coding events in the nuclear genomes provides a new insight into the translation mechanism in eukaryotes.

Analysis of chloroplast genome characteristics and codon usage bias in 14 species of Annonaceae.

For the study of species evolution, chloroplast gene expression, and transformation, the chloroplast genome is an invaluable resource. Codon usage bias (CUB) analysis is a tool that is utilized to improve gene expression and investigate evolutionary connections in genetic transformation. In this study, we analysed chloroplast genome differences, codon usage patterns and the sources of variation on CUB in 14 Annonaceae species using bioinformatics tools. The study showed that there was a significant variation in both gene sizes and numbers between the 14 species, but conservation was still maintained. It's worth noting that there were noticeable differences in the IR/SC sector boundary and the types of SSRs among the 14 species. The mono-nucleotide repeat type was the most common, with A/T repeats being more prevalent than G/C repeats. Among the different types of repeats, forward and palindromic repeats were the most abundant, followed by reverse repeats, and complement repeats were relatively rare. Codon composition analysis revealed that all 14 species had a frequency of GC lower than 50%. Additionally, it was observed that the proteins in-coding sequences of chloroplast genes tend to end with A/T at the third codon position. Among these species, 21 codons exhibited bias (RSCU > 1), and there were 8 high-frequency (HF) codons and 5 optimal codons that were identical across the species. According to the ENC-plot and Neutrality plot analysis, natural selection had less impact on the CUB of A. muricate and A. reticulata. Based on the PR2-plot, it was evident that base G had a higher frequency than C, and T had a higher frequency A. The correspondence analysis (COA) revealed that codon usage patterns different in Annonaceae.

An important consideration when expressing mAbs in Escherichiacoli.

Monoclonal antibodies (mAbs) are a driving force in the biopharmaceutical industry. Therapeutic mAbs are usually produced in mammalian cells, but there has been a push towards the use of alternative production hosts, such as Escherichia coli. When the genes encoding for a mAb heavy and light chains are codon-optimized for E. coli expression, a truncated form of the heavy chain can form along with the full-length product. In this work, the role of codon optimization in the formation of a truncated product was investigated. This study used the amino acid sequences of several therapeutic mAbs and multiple optimization algorithms. It was found that several algorithms incorporate sequences that lead to a truncated product. Approaches to avoid this truncated form are discussed.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Generation of genetic codes with 2-adic codon algebra and adaptive dynamics.

This is a brief review on modeling genetic codes with the aid of 2-adic dynamical systems. In this model amino acids are encoded by the attractors of such dynamical systems. Each genetic code is coupled to the special class of 2-adic dynamics. We consider the discrete dynamical systems, These are the iterations of a function F:Z2→Z2, where Z2 is the ring of 2-adic numbers (2-adic tree). A genetic code is characterized by the set of attractors of a function belonging to the code generating functional class. The main mathematical problem is to reduce degeneration of dynamic representation and select the optimal generating function. Here optimality can be treated in many ways. One possibility is to consider the Lipschitz functions playing the crucial role in general theory of iterations. Then we minimize the Lip-constant. The main issue is to find the proper biological interpretation of code-functions. One can speculate that the evolution of the genetic codes can be described in information space of the nucleotide-strings endowed with ultrametric (treelike) geometry. A code-function is a fitness function; the solutions of the genetic code optimization problem are attractors of the code-function. We illustrate this approach by generation of the standard nuclear and (vertebrate) mitochondrial genetics codes.

Decoding Evolution of Rubioideae: Plastomes Reveal Sweet Secrets of Codon Usage, Diagnostides, and Superbarcoding.

Galium genus belongs to the Rubiaceae family, which consists of approximately 14,000 species. In comparison to its well-known relatives, the plastomes of the Galium genus have not been explored so far. The plastomes of this genus have a typical, quadripartite structure, but differ in gene content, since the infA gene is missing in Galium palustre and Galium trfidum. An evaluation of the effectiveness of using entire chloroplast genome sequences as superbarcodes for accurate plant species identification revealed the high potential of this method for molecular delimitation within the genus and tribe. The trnE-UUC-psbD region showed the biggest number of diagnostides (diagnostic nucleotides) which might be new potential barcodes, not only in Galium, but also in other closely related genera. Relative synonymous codon usage (RSCU) appeared to be connected with the phylogeny of the Rubiaceae family, showing that during evolution, plants started preferring specific codons over others.

Comprehensive analysis of two hotspot codons in the TUBB4B gene and associated phenotypes.

Our purpose was to elucidate the genotype and ophthalmological and audiological phenotype in TUBB4B-associated inherited retinal dystrophy (IRD) and sensorineural hearing loss (SNHL), and to model the effects of all possible amino acid substitutions at the hotspot codons Arg390 and Arg391. Six patients from five families with heterozygous missense variants in TUBB4B were included in this observational study. Ophthalmological testing included best-corrected visual acuity, fundus examination, optical coherence tomography, fundus autofluorescence imaging, and full-field electroretinography (ERG). Audiological examination included pure-tone and speech audiometry in adult patients and auditory brainstem response testing in a child. Genetic testing was performed by disease gene panel analysis based on genome sequencing. The molecular consequences of the substitutions of residues 390 and 391 on TUBB4B and its interaction with α-tubulin were predicted in silico on its three-dimensional structure obtained by homology modelling. Two independent patients had amino acid exchanges at position 391 (p.(Arg391His) or p.(Arg391Cys)) of the TUBB4B protein. Both had a distinct IRD phenotype with peripheral round yellowish lesions with pigmented spots and mild or moderate SNHL, respectively. Yet the phenotype was milder with a sectorial pattern of bone spicules in one patient, likely due to a genetically confirmed mosaicism for p.(Arg391His). Three patients were heterozygous for an amino acid exchange at position 390 (p.(Arg390Gln) or p.(Arg390Trp)) and presented with another distinct retinal phenotype with well demarcated pericentral retinitis pigmentosa. All showed SNHL ranging from mild to severe. One additional patient showed a variant distinct from codon 390 or 391 (p.(Tyr310His)), and presented with congenital profound hearing loss and reduced responses in ERG. Variants at codon positions 390 and 391 were predicted to decrease the structural stability of TUBB4B and its complex with α-tubulin, as well as the complex affinity. In conclusion, the twofold larger reduction in heterodimer affinity exhibited by Arg391 substitutions suggested an association with the more severe retinal phenotype, compared to the substitution at Arg390.

Codon Optimization-based Whole-gene Scanning Identifies Hidden Nucleotides Essential for Bombyx mori Nucleopolyhedrovirus polyhedrin Hyperexpression.

During the late stage of infection, alphabaculoviruses produce many occlusion bodies (OBs) in the nuclei of the insect host's cells through the hyperexpression of polyhedrin (POLH), a major OB component encoded by polh. The strong polh promoter has been used to develop a baculovirus expression vector system for recombinant protein expression in cultured insect cells and larvae. However, the relationship between POLH accumulation and the polh coding sequence remains largely unelucidated. This study aimed to assess the importance of polh codon usage and/or nucleotide sequences in POLH accumulation by generating a baculovirus Bombyx mori nucleopolyhedrovirus (BmNPV) expressing mutant polh (co-polh) optimized according to the codon preference of its host insect. Although the deduced amino acid sequence of CO-POLH was the same as that of wild-type POLH, POLH accumulation was significantly lower in cells infected with the co-polh mutant. This reduction was due to decreased polh mRNA levels rather than translational repression. Analysis of mutant viruses with chimeric polh revealed that a 30 base-pair (bp) 5' proximal polh coding region was necessary for maintaining high polh mRNA levels. Sequence comparison of wild-type polh and co-polh identified five nucleotide differences in this region, indicating that these nucleotides were critical for polh hyperexpression. Furthermore, luciferase reporter assays showed that the 30 bp 5' coding region was sufficient for maintaining the polh promoter-driven high level of polh mRNA. Thus, our whole-gene scanning by codon optimization identified important hidden nucleotides for polh hyperexpression in alphabaculoviruses.

DNA damage induces p53-independent apoptosis through ribosome stalling.

In response to excessive DNA damage, human cells can activate p53 to induce apoptosis. Cells lacking p53 can still undergo apoptosis upon DNA damage, yet the responsible pathways are unknown. We observed that p53-independent apoptosis in response to DNA damage coincided with translation inhibition, which was characterized by ribosome stalling on rare leucine-encoding UUA codons and globally curtailed translation initiation. A genetic screen identified the transfer RNAse SLFN11 and the kinase GCN2 as factors required for UUA stalling and global translation inhibition, respectively. Stalled ribosomes activated a ribotoxic stress signal conveyed by the ribosome sensor ZAKα to the apoptosis machinery. These results provide an explanation for the frequent inactivation of SLFN11 in chemotherapy-unresponsive tumors and highlight ribosome stalling as a signaling event affecting cell fate in response to DNA damage.

Massively parallel CRISPR-assisted homologous recombination enables saturation editing of full-length endogenous genes in yeast.

Performing saturation editing of chromosomal genes will enable the study of genetic variants in situ and facilitate protein and cell engineering. However, current in vivo editing of endogenous genes either lacks flexibility or is limited to discrete codons and short gene fragments, preventing a comprehensive exploration of genotype-phenotype relationships. To enable facile saturation editing of full-length genes, we used a protospacer adjacent motif-relaxed Cas9 variant and homology-directed repair to achieve above 60% user-defined codon replacement efficiencies in Saccharomyces cerevisiae genome. Coupled with massively parallel DNA design and synthesis, we developed a saturation gene editing method termed CRISPR-Cas9- and homology-directed repair-assisted saturation editing (CHASE) and achieved highly saturated codon swapping of long genomic regions. By applying CHASE to massively edit a well-studied global transcription factor gene, we found known and unreported genetic variants affecting an industrially relevant microbial trait. The user-defined codon editing capability and wide targeting windows of CHASE substantially expand the scope of saturation gene editing.

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.