Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Deep mutational scanning of CYP2C19 in human cells reveals a substrate specificity-abundance tradeoff

The cytochrome P450s enzyme family metabolizes ∼80% of small molecule drugs. Variants in cytochrome P450s can substantially alter drug metabolism, leading to improper dosing and severe adverse drug reactions. Due to low sequence conservation, predicting variant effects across cytochrome P450s is challenging. Even closely related cytochrome P450s like CYP2C9 and CYP2C19, which share 92% amino acid sequence identity, display distinct phenotypic properties. Using variant abundance by massively parallel sequencing, we measured the steady-state protein abundance of 7,660 single amino acid variants in CYP2C19 expressed in cultured human cells. Our findings confirmed critical positions and structural features essential for cytochrome P450 function, and revealed how variants at conserved positions influence abundance. We jointly analyzed 4,670 variants whose abundance was measured in both CYP2C19 and CYP2C9, finding that the homologs have different variant abundances in substrate recognition sites within the hydrophobic core. We also measured the abundance of all single and some multiple wild type amino acid exchanges between CYP2C19 and CYP2C9. While most exchanges had no effect, substitutions in substrate recognition site 4 reduced abundance in CYP2C19. Double and triple mutants showed distinct interactions, highlighting a region that points to differing thermodynamic properties between the 2 homologs. These positions are known contributors to substrate specificity, suggesting an evolutionary tradeoff between stability and enzymatic function. Finally, we analyzed 368 previously unannotated human variants, finding that 43% had decreased abundance. By comparing variant effects between these homologs, we uncovered regions underlying their functional differences, advancing our understanding of this versatile family of enzymes.

A novel m.5906G > a variant in MT-CO1 causes MELAS/Leigh overlap syndrome

The MELAS/Leigh overlap syndrome manifests with a blend of clinical and radiographic traits from both MELAS and LS. However, the association of MELAS/Leigh overlap syndrome with MT-CO1 gene variants has not been previously reported. In this study, we report a patient diagnosed with MELAS/Leigh overlap syndrome harboring the m.5906G > A variant in MT-CO1, with biochemical evidence supporting the pathogenicity of the variant. The variant m.5906G > A that led to a synonymous variant in the start codon of MT-CO1 was filtered as the candidate disease-causing variant of the patient. Patient-derived fibroblasts were used to generate a series of monoclonal cells carrying different m.5906G > A variant loads for further functional assays. The oxygen consumption rate, ATP production, mitochondrial membrane potential and lactate assay indicated an impairment of cellular bioenergetics due to the m.5906G > A variant. Blue native PAGE analysis revealed that the m.5906G > A variant caused a deficiency in the content of mitochondrial oxidative phosphorylation complexes. Furthermore, molecular biology assays performed for the pathogenesis, mtDNA copy number, mtDNA-encoded subunits, and recovery capacity of mtDNA were all deficient due to the m.5906G > A variant, which might be caused by mtDNA replication deficiency. Overall, our findings demonstrated the pathogenicity of m.5906G > A variant and proposed a potential pathogenic mechanism, thereby expanding the genetic spectrum of MELAS/Leigh overlap syndrome.

Codon Bias of the DDR1 Gene and Transcription Factor EHF in Multiple Species

Milk production is an essential economic trait in cattle, and understanding the genetic regulation of this trait can enhance breeding strategies. The discoidin domain receptor 1 (DDR1) gene has been identified as a key candidate gene that influences milk production, and ETS homologous factor (EHF) is recognized as a critical transcription factor that regulates DDR1 expression. Codon usage bias, which affects gene expression and protein function, has not been fully explored in cattle. This study aims to examine the codon usage bias of DDR1 and EHF transcription factors to understand their roles in dairy production traits. Data from 24 species revealed that both DDR1 and EHF predominantly used G/C-ending codons, with the GC3 content averaging 75.49% for DDR1 and 61.72% for EHF. Synonymous codon usage analysis identified high-frequency codons for both DDR1 and EHF, with 17 codons common to both genes. Correlation analysis indicated a negative relationship between the effective number of codons and codon adaptation index for both DDR1 and EHF. Phylogenetic and clustering analyses revealed similar codon usage patterns among closely related species. These findings suggest that EHF plays a crucial role in regulating DDR1 expression, offering new insights into genetically regulating milk production in cattle.

Characterisation of Three Ovine KRTAP13 Family Genes and Their Association with Wool Traits in Chinese Tan Sheep

Understanding the genetic basis of wool traits is crucial for improving wool production. In this study, we investigated the ovine KAP13 gene family, which in humans contains multiple members, while only one member has been identified to date in sheep. Three ovine KRTAP13 genes, likely representing KRTAP13-1, KRTAP13-2, and KRTAP13-4, were identified through sequence analysis and phylogenetic comparisons. These genes are positioned on chromosome 1, between KRTAP15-1 and KRTAP27-1, in a pattern that is like the arrangement in humans but not identical. Analyses revealed multiple sequence variants of each gene in 356 sheep from a variety of wool, meat, and dual-purpose breeds. The effect of these genes on four fibre traits: mean fibre curvature (MFC), mean fibre diameter (MFD), coefficient of variation of fibre diameter (CVFD), and fibre diameter standard deviation (FDSD), was assessed in 240 lambs of the Chinese Tan sheep breed. An allele of KRTAP13-2 was revealed to be associated with a decrease in FDSD and CVFD in heterotypic fibres. No associations were found between KRTAP13-4 variation and wool traits, and an association analysis for KRTAP13-1 was not conducted because no variation was found in this gene in the Chinese Tan sheep studied. These findings suggest a potential role for KRTAP13-2 in regulating wool traits, particularly fibre diameter uniformity in larger heterotypic hair fibres, and suggest its potential use as a marker for improving wool traits.

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

CbCyp51-Mediated Demethylation Inhibitor Resistance Is Modulated by Codon Bias

Cercospora leaf spot, caused by the fungus Cercospora beticola, is the most destructive foliar disease of sugarbeet worldwide. Resistance to the sterol demethylation inhibitor (DMI) fungicide tetraconazole has been previously correlated with synonymous and nonsynonymous mutations in CbCyp51. Here, we extend these analyses to the DMI fungicides prothioconazole, difenoconazole, and mefentrifluconazole in addition to tetraconazole to confirm whether the synonymous and nonsynonymous mutations at amino acid positions 144 and 170 are associated with resistance to these fungicides. Nearly half of the 593 isolates of C. beticola collected in the Red River Valley of North Dakota and Minnesota in 2021 were resistant to all four DMIs. Another 20% were resistant to tetraconazole and prothioconazole but sensitive to difenoconazole and mefentrifluconazole. A total of 13% of isolates were sensitive to all DMIs tested. We found five CbCyp51 haplotypes and associated them with phenotypes to the four DMIs. The most predominant haplotype (E170_A/L144F_C) correlated with resistance to all four DMIs with up to 97.6% accuracy. The second most common haplotype (E170_A/L144) consisted of isolates associated with resistance phenotypes to tetraconazole and prothioconazole while also exhibiting sensitive phenotypes to difenoconazole and mefentrifluconazole with up to 98.4% accuracy. Quantitative PCR did not identify differences in CbCyp51 expression between haplotypes. This study offers an understanding of the importance of codon usage in fungicide resistance and provides crop management acuity for fungicide application decision-making.

Unlocking the potential of Extensin Signal peptide and Elastin-like polypeptide tag fused to Shigella dysenteriae's IpaDSTxB to improve protein expression and purification in Nicotiana tabacum and Medicagosativa

Plants are often seen as a potent tool in the recombinant protein production industry. However, unlike bacterial expression, it is not a popular method due to the low yield and difficulty of protein extraction and purification. Therefore, developing a new high efficient and easy to purify platform is crucial. One of the best approaches to make extraction easier is to utilize the Extensin Signal peptide (EXT) to translocate the recombinant protein to the outside of the cell, along with incorporating an Elastin-like polypeptide tag (ELP) to enhance purification and accumulation rates. In this research, we transiently expressed Shigella dysenteriae's IpaDSTxB fused to both NtEXT and ELP in both Nicotiana tabacum and Medicago sativa. Our results demonstrated that N. tabacum, with an average yield of 6.39 ng/μg TSP, outperforms M. sativa, which had an average yield of 3.58 ng/μg TSP. On the other hand, analyzing NtEXT signal peptide indicated that merging EXT to the constructs facilitates translocation of IpaDSTxB to the apoplast by 78.4% and 65.9% in N. tabacum and M. sativa, respectively. Conversely, the mean level for constructs without EXT was below 25% for both plants. Furthermore, investigation into the orientation of ELP showed that merging it to the C-terminal of IpaDSTxB leads to a higher accumulation rate in both N. tabacum and M. sativa by 1.39 and 1.28 times, respectively. It also facilitates purification rate by over 70% in comparison to 20% of the 6His tag. The results show a highly efficient and easy to purify platform for the expression of heterologous proteins in plant.

Characterisation of Ovine KRTAP19-3 and Its Impact on Wool Traits in Chinese Tan Sheep

Wool, a natural fibre derived from sheep, can present a challenge to wool processing and manufacturing industries because of the variation in fibre traits. Genetic improvement offers one solution to this challenge, and having a better understanding of the genes that affect wool fibre traits is therefore important. Here, we describe ovine KRTAP19-3, a new member of the KAP19 gene family. Phylogenetic analysis revealed its relationship to other known KRTAP19 gene sequences, and an analysis of the nucleotide sequence variation in KRTAP19-3 from 288 sheep of a variety of breeds revealed six unique variant sequences. Among these variants, eleven single nucleotide polymorphisms (SNPs) were detected, with six located in the coding region. Three of these coding region SNPs were non-synonymous and would result in amino acid changes. Associations were observed between the presence of specific sequence variants in Chinese Tan sheep and wool trait variation, particularly an increase in fibre diameter variability in the heterotypic hair fibres. These findings enhance our understanding of the genes that encode sheep wool proteins.

A synonymous mutation of rs1137070 cause the mice Maoa gene transcription and translation to decrease

The Monoamine Oxidase-A (MAOA) EcoRV polymorphism (rs1137070) is a unique synonymous mutation (c.1409 T > C) within the MAOA gene, which plays a crucial role in Maoa gene expression and function. This study aimed to explore the relationship between the mouse Maoa rs1137070 genotype and differences in MAOA gene expression. Mice carrying the CC genotype of rs1137070 exhibited a significantly lower Maoa expression level, with an odds ratio of 2.44 compared to the T carriers. Moreover, the wild-type TT genotype of MAOA demonstrated elevated mRNA expression and a longer half-life. We also delved into the significant expression and structural disparities among genotypes. Furthermore, it was evident that different aspartic acid synonymous codons within Maoa influenced both MAOA expression and enzyme activity, highlighting the association between rs1137070 and MAOA. To substantiate these findings, a dual-luciferase reporter assay confirmed that GAC was more efficient than GAT binding. Conversely, the synonymous mutation altered Maoa gene expression in individual mice. An RNA pull-down assay suggested that this alteration could impact the interaction with RNA-binding proteins. In summary, our results illustrate that synonymous mutations can indeed regulate the downregulation of gene expression, leading to changes in MAOA function and their potential association with neurological-related diseases.

Newly developed mRNA vaccines induce immune responses in Litopenaeus vannamei shrimps during primary vaccination

White spot syndrome virus (WSSV) causes highly destructive infection in crustacean aquaculture, often resulting in 100% mortality within a week. However, there is lack of studies addressing the safety issues of WSSV vaccines in shrimps. In this study, WSSV VP28 mRNA vaccines were developed using codon deoptimization approach. These vaccines were administered to Litopenaeus vannamei shrimps at various dosages to access their safety and the shrimps' immune responses using quantification PCR (qPCR). The findings of this study indicate that the expression level of codon deoptimized VP28 mRNA vaccines are lower compared to the wild type VP28 vaccines, as observed through a comparison of bioinformatic predictions and experimental results. Additionally, the total haemocyte count (THC) in shrimps injected with codon deoptimized VP28 vaccine was higher than those injected with wild type VP28 vaccines. Furthermore, the expression of immune-related genes differed between codon deoptimized and wild type VP28 vaccines. In summary, the results suggest that 0.01 μg codon deoptimized VP28-D1 mRNA vaccine is the most promising WSSV mRNA vaccine, displaying low pathogenicity and expression in shrimps. To the best of our knowledge, this research represents the first attempt to attenuate WSSV using codon deoptimization method and development of a potential mRNA vaccine for shrimp purpose. The study addresses an important gap in shrimp vaccine research, offering potential solutions for WSSV control in shrimps.

Assessing Mechanisms of Potential Local Adaptation Through a Seascape Genomic Approach in a Marine Gastropod, Littoraria flava

Understanding the combined effects of environmental heterogeneity and evolutionary processes on marine populations is a primary goal of seascape genomic approaches. Here, we utilized genomic approaches to identify local adaptation signatures in Littoraria flava, a widely distributed marine gastropod in the tropical West Atlantic population. We also performed molecular evolution analyses to investigate potential selective signals across the genome. After obtaining 6,298 and 16,137 single nucleotide polymorphisms derived from genotyping-by-sequencing and RNA sequencing, respectively, 69 from genotyping-by-sequencing (85 specimens) and four from RNA sequencing (40 specimens) candidate single nucleotide polymorphisms were selected and further evaluated. The correlation analyses support different evolutionary pressures over transcribed and non-transcribed regions. Thus, single nucleotide polymorphisms within transcribed regions could account for the genotypic and possibly phenotypic divergences in periwinkles. Our molecular evolution tests based on synonymous and non-synonymous ratio (kN/kS) showed that genotype divergences containing putative adaptive single nucleotide polymorphisms arose mainly from synonymous and/or UTR substitutions rather than polymorphic proteins. The distribution of genotypes across different localities seems to be influenced by marine currents, pH, and temperature variations, suggesting that these factors may impact the species dispersion. The combination of RNA sequencing and genotyping-by-sequencing derived datasets provides a deeper understanding of the molecular mechanisms underlying selective forces responses on distinct genomic regions and could guide further investigations on seascape genomics for non-model species.

Synonymous codon substitutions modulate transcription and translation of a divergent upstream gene by modulating antisense RNA production

Synonymous codons were originally viewed as interchangeable, with no phenotypic consequences. However, substantial evidence has now demonstrated that synonymous substitutions can perturb a variety of gene expression and protein homeostasis mechanisms, including translational efficiency, translational fidelity, and cotranslational folding of the encoded protein. To date, most studies of synonymous codon-derived perturbations have focused on effects within a single gene. Here, we show that synonymous codon substitutions made far within the coding sequence of Escherichia coli plasmid-encoded chloramphenicol acetyltransferase (cat) can significantly increase expression of the divergent upstream tetracycline resistance gene, tetR. In four out of nine synonymously recoded cat sequences tested, expression of the upstream tetR gene was significantly elevated due to transcription of a long antisense RNA (asRNA) originating from a transcription start site within cat. Surprisingly, transcription of this asRNA readily bypassed the native tet transcriptional repression mechanism. Even more surprisingly, accumulation of the TetR protein correlated with the level of asRNA, rather than total tetR RNA. These effects of synonymous codon substitutions on transcription and translation of a neighboring gene suggest that synonymous codon usage in bacteria may be under selection to both preserve the amino acid sequence of the encoded gene and avoid DNA sequence elements that can significantly perturb expression of neighboring genes. Avoiding such sequences may be especially important in plasmids and prokaryotic genomes, where genes and regulatory elements are often densely packed. Similar considerations may apply to the design of genetic circuits for synthetic biology applications.

Synonymous rpsH variants: the common denominator in Escherichia coli adapting to ionizing radiation

Ionizing radiation (IR) in high doses is generally lethal to most organisms. Investigating mechanisms of radiation resistance is crucial for gaining insights into the underlying cellular responses and understanding the damaging effects of IR. In this study, we conducted a comprehensive analysis of sequencing data from an evolutionary experiment aimed at understanding the genetic adaptations to ionizing radiation in Escherichia coli. By including previously neglected synonymous mutations, we identified the rpsH c.294T > G variant, which emerged in all 17 examined isolates across four subpopulations. The identified variant is a synonymous mutation affecting the 30S ribosomal protein S8, and consistently exhibited high detection and low allele frequencies in all subpopulations. This variant, along with two additional rpsH variants, potentially influences translational control of the ribosomal spc operon. The early emergence and stability of these variants suggest their role in adapting to environmental stress, possibly contributing to radiation resistance. Our findings shed light on the dynamics of ribosomal variants during the evolutionary process and their potential role in stress adaptation, providing valuable implications for understanding clinical radiation sensitivity and improving radiotherapy.

Advancements and Future Prospects of CRISPR-Cas-Based Population Replacement Strategies in Insect Pest Management

Population replacement refers to the process by which a wild-type population of insect pests is replaced by a population possessing modified traits or abilities. Effective population replacement necessitates a gene drive system capable of spreading desired genes within natural populations, operating under principles akin to super-Mendelian inheritance. Consequently, releasing a small number of genetically edited insects could potentially achieve population control objectives. Currently, several gene drive approaches are under exploration, including the newly adapted CRISPR-Cas genome editing system. Multiple studies are investigating methods to engineer pests that are incapable of causing crop damage or transmitting vector-borne diseases, with several notable successful examples documented. This review summarizes the recent advancements of the CRISPR-Cas system in the realm of population replacement and provides insights into research methodologies, testing protocols, and implementation strategies for gene drive techniques. The review also discusses emerging trends and prospects for establishing genetic tools in pest management.

The Development of Epitope-Based Recombinant Protein Vaccines against SARS-CoV-2

The COVID-19 pandemic continues to cause infections and deaths, which are attributable to the SARS-CoV-2 Omicron variant of concern (VOC). Moderna's response to the declining protective efficacies of current SARS-CoV-2 vaccines against Omicron was to develop a bivalent booster vaccine based on the Spike (S) protein from the Wuhan and Omicron BA.4/BA.5 strains. This approach, while commendable, is unfeasible in light of rapidly emerging mutated viral strains. PubMed and Google Scholar were systematically reviewed for peer-reviewed papers up to January 2024. Articles included focused on specific themes such as the clinical history of recombinant protein vaccine development against different diseases, including COVID-19, the production of recombinant protein vaccines using different host expression systems, aspects to consider in recombinant protein vaccine development, and overcoming problems associated with large-scale recombinant protein vaccine production. In silico approaches to identify conserved and immunogenic epitopes could provide broad protection against SARS-CoV-2 VOCs but require validation in animal models. The recombinant protein vaccine development platform has shown a successful history in clinical development. Recombinant protein vaccines incorporating conserved epitopes may utilize a number of expression systems, such as yeast (Saccharomyces cerevisiae), baculovirus-insect cells (Sf9 cells), and Escherichia coli (E. coli). Current multi-epitope subunit vaccines against SARS-CoV-2 utilizing synthetic peptides are unfeasible for large-scale immunizations. Recombinant protein vaccines based on conserved and immunogenic proteins produced using E. coli offer high production yields, convenient purification, and cost-effective production of large-scale vaccine quantities capable of protecting against the SARS-CoV-2 D614G strain and its VOCs.

Comparative Analysis of Codon Usage Bias in Six Eimeria Genomes

The codon usage bias (CUB) of genes encoded by different species' genomes varies greatly. The analysis of codon usage patterns enriches our comprehension of genetic and evolutionary characteristics across diverse species. In this study, we performed a genome-wide analysis of CUB and its influencing factors in six sequenced Eimeria species that cause coccidiosis in poultry: Eimeria acervulina, Eimeria necatrix, Eimeria brunetti, Eimeria tenella, Eimeria praecox, and Eimeria maxima. The GC content of protein-coding genes varies between 52.67% and 58.24% among the six Eimeria species. The distribution trend of GC content at different codon positions follows GC1 > GC3 > GC2. Most high-frequency codons tend to end with C/G, except in E. maxima. Additionally, there is a positive correlation between GC3 content and GC3s/C3s, but a significantly negative correlation with A3s. Analysis of the ENC-Plot, neutrality plot, and PR2-bias plot suggests that selection pressure has a stronger influence than mutational pressure on CUB in the six Eimeria genomes. Finally, we identified from 11 to 15 optimal codons, with GCA, CAG, and AGC being the most commonly used optimal codons across these species. This study offers a thorough exploration of the relationships between CUB and selection pressures within the protein-coding genes of Eimeria species. Genetic evolution in these species appears to be influenced by mutations and selection pressures. Additionally, the findings shed light on unique characteristics and evolutionary traits specific to the six Eimeria species.

The Essential Functions of Molecular Chaperones and Folding Enzymes in Maintaining Endoplasmic Reticulum Homeostasis

It has been estimated that up to one-third of the proteins encoded by the human genome enter the endoplasmic reticulum (ER) as extended polypeptide chains where they undergo covalent modifications, fold into their native structures, and assemble into oligomeric protein complexes. The fidelity of these processes is critical to support organellar, cellular, and organismal health, and is perhaps best underscored by the growing number of disease-causing mutations that reduce the fidelity of protein biogenesis in the ER. To meet demands encountered by the diverse protein clientele that mature in the ER, this organelle is populated with a cadre of molecular chaperones that prevent protein aggregation, facilitate protein disulfide isomerization, and lower the activation energy barrier of cis-trans prolyl isomerization. Components of the lectin (glycan-binding) chaperone system also reside within the ER and play numerous roles during protein biogenesis. In addition, the ER houses multiple homologs of select chaperones that can recognize and act upon diverse peptide signatures. Moreover, redundancy helps ensure that folding-compromised substrates are unable to overwhelm essential ER-resident chaperones and enzymes. In contrast, the ER in higher eukaryotic cells possesses a single member of the Hsp70, Hsp90, and Hsp110 chaperone families, even though several homologs of these molecules reside in the cytoplasm. In this review, we discuss specific functions of the many factors that maintain ER quality control, highlight some of their interactions, and describe the vulnerabilities that arise from the absence of multiple members of some chaperone families.

Translation Rates and Protein Folding

The mRNA coding sequence defines not only the amino acid sequence of the protein, but also the speed at which the ribosomes move along the mRNA while making the protein. The non-uniform local kinetics - denoted as translational rhythm - is similar among mRNAs coding for related protein folds. Deviations from this conserved rhythm can result in protein misfolding. In this review we summarize the experimental evidence demonstrating how local translation rates affect cotranslational protein folding, with the focus on the synonymous codons and patches of charged residues in the nascent peptide as best-studied examples. Alterations in nascent protein conformations due to disturbed translational rhythm can persist off the ribosome, as demonstrated by the effects of synonymous codon variants of several disease-related proteins. Charged amino acid patches in nascent chains also modulate translation and cotranslational protein folding, and can abrogate translation when placed at the N-terminus of the nascent peptide. During cotranslational folding, incomplete nascent chains navigate through a unique conformational landscape in which earlier intermediate states become inaccessible as the nascent peptide grows. Precisely tuned local translation rates, as well as interactions with the ribosome, guide the folding pathway towards the native structure, whereas deviations from the natural translation rhythm may favor pathways leading to trapped misfolded states. Deciphering the 'folding code' of the mRNA will contribute to understanding the diseases caused by protein misfolding and to rational protein design.

The Effects of Codon Usage on Protein Structure and Folding

The rate of protein synthesis is slower than many folding reactions and varies depending on the synonymous codons encoding the protein sequence. Synonymous codon substitutions thus have the potential to regulate cotranslational protein folding mechanisms, and a growing number of proteins have been identified with folding mechanisms sensitive to codon usage. Typically, these proteins have complex folding pathways and kinetically stable native structures. Kinetically stable proteins may fold only once over their lifetime, and thus, codon-mediated regulation of the pioneer round of protein folding can have a lasting impact. Supporting an important role for codon usage in folding, conserved patterns of codon usage appear in homologous gene families, hinting at selection. Despite these exciting developments, there remains few experimental methods capable of quantifying translation elongation rates and cotranslational folding mechanisms in the cell, which challenges the development of a predictive understanding of how biology uses codons to regulate protein folding.

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.

Navigating the complexities of multi-domain protein folding

Proteome complexity has expanded tremendously over evolutionary time, enabling biological diversification. Much of this complexity is achieved by combining a limited set of structural units into long polypeptides. This widely used evolutionary strategy poses challenges for folding of the resulting multi-domain proteins. As a consequence, their folding differs from that of small single-domain proteins, which generally fold quickly and reversibly. Co-translational processes and chaperone interactions are important aspects of multi-domain protein folding. In this review, we discuss some of the recent experimental progress toward understanding these processes.

Optimizing enzyme properties to enhance dihydroxyacetone production via methylglyoxal biosensor development

BACKGROUND

Dihydroxyacetone (DHA) stands as a crucial chemical material extensively utilized in the cosmetics industry. DHA production through the dephosphorylation of dihydroxyacetone phosphate, an intermediate product of the glycolysis pathway in Escherichia coli, presents a prospective alternative for industrial production. However, insights into the pivotal enzyme, dihydroxyacetone phosphate dephosphorylase (HdpA), remain limited for informed engineering. Consequently, the development of an efficient tool for high-throughput screening of HdpA hypermutants becomes imperative.

RESULTS

This study introduces a methylglyoxal biosensor, based on the formaldehyde-responding regulator FrmR, for the selection of HdpA. Initial modifications involved the insertion of the FrmR binding site upstream of the -35 region and into the spacer region between the -10 and -35 regions of the constitutive promoter J23110. Although the hybrid promoter retained constitutive expression, expression of FrmR led to complete repression. The addition of 350 μM methylglyoxal promptly alleviated FrmR inhibition, enhancing promoter activity by more than 40-fold. The methylglyoxal biosensor system exhibited a gradual increase in fluorescence intensity with methylglyoxal concentrations ranging from 10 to 500 μM. Notably, the biosensor system responded to methylglyoxal spontaneously converted from added DHA, facilitating the separation of DHA producing and non-producing strains through flow cytometry sorting. Subsequently, the methylglyoxal biosensor was successfully applied to screen a library of HdpA mutants, identifying two strains harboring specific mutants 267G > T and D110G/G151C that showed improved DHA production by 68% and 114%, respectively. Expressing of these two HdpA mutants directly in a DHA-producing strain also increased DHA production from 1.45 to 1.92 and 2.29 g/L, respectively, demonstrating the enhanced enzyme properties of the HdpA mutants.

CONCLUSIONS

The methylglyoxal biosensor offers a novel strategy for constructing genetically encoded biosensors and serves as a robust platform for indirectly determining DHA levels by responding to methylglyoxal. This property enables efficiently screening of HdpA hypermutants to enhance DHA production.

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

Translational adaptation in breast cancer metastasis and emerging therapeutic opportunities

Breast cancer's tendency to metastasize poses a critical barrier to effective treatment, making it a leading cause of mortality among women worldwide. A growing body of evidence is showing that translational adaptation is emerging as a key mechanism enabling cancer cells to thrive in the dynamic tumor microenvironment (TME). Here, we systematically summarize how breast cancer cells utilize translational adaptation to drive metastasis, highlighting the intricate regulation by specific translation machinery and mRNA attributes such as sequences and structures, along with the involvement of tRNAs and other trans-acting RNAs. We provide an overview of the latest findings and emerging concepts in this area, discussing their potential implications for therapeutic strategies in breast cancer.

"Metabolic burden" explained: stress symptoms and its related responses induced by (over)expression of (heterologous) proteins in Escherichia coli

BACKGROUND

Engineering bacterial strains to redirect the metabolism towards the production of a specific product has enabled the development of industrial biotechnology. However, rewiring the metabolism can have severe implications for a microorganism, rendering cells with stress symptoms such as a decreased growth rate, impaired protein synthesis, genetic instability and an aberrant cell size. On an industrial scale, this is reflected in processes that are not economically viable.

MAIN TEXT

In literature, most stress symptoms are attributed to "metabolic burden", however the actual triggers and stress mechanisms involved are poorly understood. Therefore, in this literature review, we aimed to get a better insight in how metabolic engineering affects Escherichia coli and link the observed stress symptoms to its cause. Understanding the possible implications that chosen engineering strategies have, will help to guide the reader towards optimising the envisioned process more efficiently.

CONCLUSION

This review addresses the gap in literature and discusses the triggers and effects of stress mechanisms that can be activated when (over)expressing (heterologous) proteins in Escherichia coli. It uncovers that the activation of the different stress mechanisms is complex and that many are interconnected. The reader is shown that care has to be taken when (over)expressing (heterologous) proteins as the cell's metabolism is tightly regulated.

Non-Coding RNAs: Regulators of Stress, Ageing, and Developmental Decisions in Yeast?

Cells must change their properties in order to adapt to a constantly changing environment. Most of the cellular sensing and regulatory mechanisms described so far are based on proteins that serve as sensors, signal transducers, and effectors of signalling pathways, resulting in altered cell physiology. In recent years, however, remarkable examples of the critical role of non-coding RNAs in some of these regulatory pathways have been described in various organisms. In this review, we focus on all classes of non-coding RNAs that play regulatory roles during stress response, starvation, and ageing in different yeast species as well as in structured yeast populations. Such regulation can occur, for example, by modulating the amount and functional state of tRNAs, rRNAs, or snRNAs that are directly involved in the processes of translation and splicing. In addition, long non-coding RNAs and microRNA-like molecules are bona fide regulators of the expression of their target genes. Non-coding RNAs thus represent an additional level of cellular regulation that is gradually being uncovered.

DNA flexibility can shape the preferential hypermutation of antibody genes

Antibody-coding genes accumulate somatic mutations to achieve antibody affinity maturation. Genetic dissection using various mouse models has shown that intrinsic hypermutations occur preferentially and are predisposed in the DNA region encoding antigen-contacting residues. The molecular basis of nonrandom/preferential mutations is a long-sought question in the field. Here, we summarize recent findings on how single-strand (ss)DNA flexibility facilitates activation-induced cytidine deaminase (AID) activity and fine-tunes the mutation rates at a mesoscale within the antibody variable domain exon. We propose that antibody coding sequences are selected based on mutability during the evolution of adaptive immunity and that DNA mechanics play a noncoding role in the genome. The mechanics code may also determine other cellular DNA metabolism processes, which awaits future investigation.

Unravelling molecular dynamics in living cells: Fluorescent protein biosensors for cell biology

Genetically encoded, fluorescent protein (FP)-based Förster resonance energy transfer (FRET) biosensors are microscopy imaging tools tailored for the precise monitoring and detection of molecular dynamics within subcellular microenvironments. They are characterised by their ability to provide an outstanding combination of spatial and temporal resolutions in live-cell microscopy. In this review, we begin by tracing back on the historical development of genetically encoded FP labelling for detection in live cells, which lead us to the development of early biosensors and finally to the engineering of single-chain FRET-based biosensors that have become the state-of-the-art today. Ultimately, this review delves into the fundamental principles of FRET and the design strategies underpinning FRET-based biosensors, discusses their diverse applications and addresses the distinct challenges associated with their implementation. We place particular emphasis on single-chain FRET biosensors for the Rho family of guanosine triphosphate hydrolases (GTPases), pointing to their historical role in driving our understanding of the molecular dynamics of this important class of signalling proteins and revealing the intricate relationships and regulatory mechanisms that comprise Rho GTPase biology in living cells.

Co-translational Assembly Pathways of Nuclear Multiprotein Complexes Involved in the Regulation of Gene Transcription

Most factors that regulate gene transcription in eukaryotic cells are multimeric, often large, protein complexes. The understanding of the biogenesis pathways of such large and heterogeneous protein assemblies, as well as the dimerization partner choice among transcription factors, is crucial to interpret and control gene expression programs and consequent cell fate decisions. Co-translational assembly (Co-TA) is thought to play key roles in the biogenesis of protein complexes by directing complex formation during protein synthesis. In this review we discuss the principles of Co-TA with a special focus for the assembly of transcription regulatory complexes. We outline the expected molecular advantages of establishing co-translational interactions, pointing at the available, or missing, evidence for each of them. We hypothesize different molecular mechanisms based on Co-TA to explain the allocation "dilemma" of paralog proteins and subunits shared by different transcription complexes. By taking as a paradigm the different assembly pathways employed by three related transcription regulatory complexes (TFIID, SAGA and ATAC), we discuss alternative Co-TA strategies for nuclear multiprotein complexes and the widespread - yet specific - use of Co-TA for the formation of nuclear complexes involved in gene transcription. Ultimately, we outlined a series of open questions which demand well-defined lines of research to investigate the principles of gene regulation that rely on the coordinated assembly of protein complexes.

Systematic analysis of nonprogrammed frameshift suppression in E. coli via translational tiling proteomics

The synthesis of proteins as encoded in the genome depends critically on translational fidelity. Nevertheless, errors inevitably occur, and those that result in reading frame shifts are particularly consequential because the resulting polypeptides are typically nonfunctional. Despite the generally maladaptive impact of such errors, the proper decoding of certain mRNAs, including many viral mRNAs, depends on a process known as programmed ribosomal frameshifting. The fact that these programmed events, commonly involving a shift to the -1 frame, occur at specific evolutionarily optimized "slippery" sites has facilitated mechanistic investigation. By contrast, less is known about the scope and nature of error (i.e., nonprogrammed) frameshifting. Here, we examine error frameshifting by monitoring spontaneous frameshift events that suppress the effects of single base pair deletions affecting two unrelated test proteins. To map the precise sites of frameshifting, we developed a targeted mass spectrometry-based method called "translational tiling proteomics" for interrogating the full set of possible -1 slippage events that could produce the observed frameshift suppression. Surprisingly, such events occur at many sites along the transcripts, involving up to one half of the available codons. Only a subset of these resembled canonical "slippery" sites, implicating alternative mechanisms potentially involving noncognate mispairing events. Additionally, the aggregate frequency of these events (ranging from 1 to 10% in our test cases) was higher than we might have anticipated. Our findings point to an unexpected degree of mechanistic diversity among ribosomal frameshifting events and suggest that frameshifted products may contribute more significantly to the proteome than generally assumed.

Novel Plastid Genome Characteristics in Fugacium kawagutii and the Trend of Accelerated Evolution of Plastid Proteins in Dinoflagellates

Typical (peridinin-containing) dinoflagellates possess plastid genomes composed of small plasmids named "minicircles". Despite the ecological importance of dinoflagellate photosynthesis in corals and marine ecosystems, the structural characteristics, replication dynamics, and evolutionary forcing of dinoflagellate plastid genomes remain poorly understood. Here, we sequenced the plastid genome of the symbiodiniacean species Fugacium kawagutii and conducted comparative analyses. We identified psbT-coding minicircles, features previously not found in Symbiodiniaceae. The copy number of F. kawagutii minicircles showed a strong diel dynamics, changing between 3.89 and 34.3 copies/cell and peaking in mid-light period. We found that F. kawagutii minicircles are the shortest among all dinoflagellates examined to date. Besides, the core regions of the minicircles are highly conserved within genus in Symbiodiniaceae. Furthermore, the codon usage bias of the plastid genomes in Heterocapsaceae, Amphidiniaceae, and Prorocentraceae species are greatly influenced by selection pressure, and in Pyrocystaceae, Symbiodiniaceae, Peridiniaceae, and Ceratiaceae species are influenced by both natural selection pressure and mutation pressure, indicating a family-level distinction in codon usage evolution in dinoflagellates. Phylogenetic analysis using 12 plastid-encoded proteins and five nucleus-encoded plastid proteins revealed accelerated evolution trend of both plastid- and nucleus-encoded plastid proteins in peridinin- and fucoxanthin-dinoflagellate plastids compared to plastid proteins of nondinoflagellate algae. These findings shed new light on the structure and evolution of plastid genomes in dinoflagellates, which will facilitate further studies on the evolutionary forcing and function of the diverse dinoflagellate plastids. The accelerated evolution documented here suggests plastid-encoded sequences are potentially useful for resolving closely related dinoflagellates.

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.

A founder COL4A3 pathogenic variant resulting in Alport syndrome and thin basement membrane disease: a case report series

Alport syndrome is a rare genetic condition characterized by kidney disease, hearing impairment, and ocular abnormalities. It exhibits various inheritance patterns involving pathogenic variants in COL4A3, COL4A4, and COL4A5 genes. The phenotypes can range from isolated hematuria with a non-progressive or very slowly progressive course to progressive kidney disease with extrarenal abnormalities. Timely diagnosis of Alport syndrome facilitates the early and effective implementation of treatment, as well as genetic counseling. Here, we report the COL4A3 c.765G > A, p.((=)) mutation in three ethnically Azerbaijani, apparently unrelated, consanguineous families from the village of Algeti in the Marneuli region of Georgia. We speculate that this variant could represent a founder mutation within this population and recommend offering genetic testing to Algeti village residents with persistent hematuria.

A SNP in the ovine cathepsin K (CTSK) gene is associated with yearling growth performance in a crossbred sheep population

Cathepsin K (CTSK) is a lysosomal protease existent in the skeletal muscles which is involved in biochemical processes related to obesity. Several studies have reported the effects of CTSK gene on body weight and fat deposition in human, mice and pigs. However, information about its structure and functions in sheep is very limited. Thus, this study was performed to evaluate the association between CTSK gene variants and yearling growth performance in Afshari × Booroola-Merino crossbred sheep. A fragment of 500 bp in exon 6 and partial of intron 5 of CTSK gene was amplified with polymerase chain reaction (PCR). All animals were genotyped by single-stranded conformation polymorphism (SSCP) and further confirmed by sequencing. Association analysis using a fixed linear model indicated that g.106510225G > A SNP was significantly related to average daily weight gain (ADWG) per year, fat-tail weight to carcass weight ratio (FW/CW), muscle thickness (MT) and muscle cross-sectional area (MCSA) of animals (p < 0.05). Due to the low polymorphic information content (PIC <0.25) for targeted locus in studied population, more association studies are needed to confirm the CTSK gene effects on growth traits in sheep.

Genome mining yields putative disease-associated ROMK variants with distinct defects

Bartter syndrome is a group of rare genetic disorders that compromise kidney function by impairing electrolyte reabsorption. Left untreated, the resulting hyponatremia, hypokalemia, and dehydration can be fatal, and there is currently no cure. Bartter syndrome type II specifically arises from mutations in KCNJ1, which encodes the renal outer medullary potassium channel, ROMK. Over 40 Bartter syndrome-associated mutations in KCNJ1 have been identified, yet their molecular defects are mostly uncharacterized. Nevertheless, a subset of disease-linked mutations compromise ROMK folding in the endoplasmic reticulum (ER), which in turn results in premature degradation via the ER associated degradation (ERAD) pathway. To identify uncharacterized human variants that might similarly lead to premature degradation and thus disease, we mined three genomic databases. First, phenotypic data in the UK Biobank were analyzed using a recently developed computational platform to identify individuals carrying KCNJ1 variants with clinical features consistent with Bartter syndrome type II. In parallel, we examined genomic data in both the NIH TOPMed and ClinVar databases with the aid of Rhapsody, a verified computational algorithm that predicts mutation pathogenicity and disease severity. Subsequent phenotypic studies using a yeast screen to assess ROMK function-and analyses of ROMK biogenesis in yeast and human cells-identified four previously uncharacterized mutations. Among these, one mutation uncovered from the two parallel approaches (G228E) destabilized ROMK and targeted it for ERAD, resulting in reduced cell surface expression. Another mutation (T300R) was ERAD-resistant, but defects in channel activity were apparent based on two-electrode voltage clamp measurements in X. laevis oocytes. Together, our results outline a new computational and experimental pipeline that can be applied to identify disease-associated alleles linked to a range of other potassium channels, and further our understanding of the ROMK structure-function relationship that may aid future therapeutic strategies to advance precision medicine.

Analysis for codon usage bias in membrane anchor of nonstructural protein 5A from BVDV

The nonstructural protein 5A (NS5A) of the bovine viral diarrhea virus (BVDV) is a monotopic membrane protein. This protein can anchor to the cell membrane by an in-plane amphipathic ⍺-helix, which participates in the viral replication complex. In this study, the effects of synonymous codon usage pattern of NS5A and the overall transfer RNA (tRNA) abundance in cells on the formation of the in-plane membrane anchor of NS5A were analyzed, based on NS5A coding sequences of different BVDV genotypes. BVDV NS5A coding sequences represent the most potential for BVDV genotyping. Moreover, the nucleotide usage of BVDV NS5A dominates the genotype-specific pattern of synonymous codon usage. There is an obvious relationship between synonymous codon usage bias and the spatial conformation of the in-plane membrane anchor. Furthermore, the overall tRNA abundance profiling displays that codon positions with a high level of tRNA abundance are more than ones with a low level of tRNA abundance in the in-plane membrane anchor, implying that high translation speed probably acts on the spatial conformation of in-plane membrane anchor of BVDV NS5A. These results give a new opinion on the effect of codon usage bias in the formation of the in-plane membrane anchor of BVDV NS5A.

Single-molecule imaging reveals distinct elongation and frameshifting dynamics between frames of expanded RNA repeats in C9ORF72-ALS/FTD

C9ORF72 hexanucleotide repeat expansion is the most common genetic cause of both amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). One pathogenic mechanism is the accumulation of toxic dipeptide repeat (DPR) proteins like poly-GA, GP and GR, produced by the noncanonical translation of the expanded RNA repeats. However, how different DPRs are synthesized remains elusive. Here, we use single-molecule imaging techniques to directly measure the translation dynamics of different DPRs. Besides initiation, translation elongation rates vary drastically between different frames, with GP slower than GA and GR the slowest. We directly visualize frameshift events using a two-color single-molecule translation assay. The repeat expansion enhances frameshifting, but the overall frequency is low. There is a higher chance of GR-to-GA shift than in the reversed direction. Finally, the ribosome-associated protein quality control (RQC) factors ZNF598 and Pelota modulate the translation dynamics, and the repeat RNA sequence is important for invoking the RQC pathway. This study reveals that multiple translation steps modulate the final DPR production. Understanding repeat RNA translation is critically important to decipher the DPR-mediated pathogenesis and identify potential therapeutic targets in C9ORF72-ALS/FTD.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Base editor screens for in situ mutational scanning at scale

A fundamental challenge in biology is understanding the molecular details of protein function. How mutations alter protein activity, regulation, and response to drugs is of critical importance to human health. Recent years have seen the emergence of pooled base editor screens for in situ mutational scanning: the interrogation of protein sequence-function relationships by directly perturbing endogenous proteins in live cells. These studies have revealed the effects of disease-associated mutations, discovered novel drug resistance mechanisms, and generated biochemical insights into protein function. Here, we discuss how this "base editor scanning" approach has been applied to diverse biological questions, compare it with alternative techniques, and describe the emerging challenges that must be addressed to maximize its utility. Given its broad applicability toward profiling mutations across the proteome, base editor scanning promises to revolutionize the investigation of proteins in their native contexts.

Synonymous Variants of Uncertain Silence

Synonymous variants, traditionally regarded as silent mutations due to their lack of impact on protein sequence, structure and function, have been the subject of increasing scrutiny. This commentary explores the emerging evidence challenging the notion of synonymous variants as functionally inert. Analysis of the activity of 70 synonymous variants in the HIV Tat transcription factor revealed that 50% of the variants exhibited significant deviations from wild-type activity. Our analysis supports previous work and raises important questions about the broader impact of non-silent synonymous variants in human genes. Considering the potential functional implications, the authors propose classifying such variants as "synonymous variants of uncertain silence" (sVUS), highlighting the need for cautious interpretation and further investigations in clinical and genetic testing settings.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Sustainable and high-level microbial production of plant hemoglobin in Corynebacterium glutamicum

BACKGROUND

Plant hemoglobin shows great potential as a food additive to circumvent the controversy of using animal materials. Microbial fermentation with engineered microorganisms is considered as a promising strategy for sustainable production of hemoglobin. As an endotoxin-free and GRAS (generally regarded as safe) bacterium, Corynebacterium glutamicum is an attractive host for hemoglobin biosynthesis.

RESULTS

Herein, C. glutamicum was engineered to efficiently produce plant hemoglobin. Hemoglobin genes from different sources including soybean and maize were selected and subjected to codon optimization. Interestingly, some candidates optimized for the codon usage bias of Escherichia coli outperformed those for C. glutamicum regarding the heterologous expression in C. glutamicum. Then, saturated synonymous mutation of the N-terminal coding sequences of hemoglobin genes and fluorescence-based high-throughput screening produced variants with 1.66- to 3.45-fold increase in hemoglobin expression level. To avoid the use of toxic inducers, such as isopropyl-β-D-thiogalactopyranoside, two native inducible expression systems based on food additives propionate and gluconate were developed. Promoter engineering improved the hemoglobin expression level by 2.2- to 12.2-fold. Combination of these strategies and plasmid copy number modification allowed intracellular production of hemoglobin up to approximately 20% of total protein. Transcriptome and proteome analyses of the hemoglobin-producing strain revealed the cellular response to excess hemoglobin accumulation. Several genes were identified as potential targets for further enhancing hemoglobin production.

CONCLUSIONS

In this study, production of plant hemoglobin in C. glutamicum was systematically engineered by combining codon optimization, promoter engineering, plasmid copy number modification, and multi-omics-guided novel target discovery. This study offers useful design principles to genetically engineer C. glutamicum for the production of hemoglobin and other recombinant proteins.

KATP channel mutations in congenital hyperinsulinism: Progress and challenges towards mechanism-based therapies

Congenital hyperinsulinism (CHI) is the most common cause of persistent hypoglycemia in infancy/childhood and is a serious condition associated with severe recurrent attacks of hypoglycemia due to dysregulated insulin secretion. Timely diagnosis and effective treatment are crucial to prevent severe hypoglycemia that may lead to life-long neurological complications. In pancreatic β-cells, adenosine triphosphate (ATP)-sensitive K+ (KATP) channels are a central regulator of insulin secretion vital for glucose homeostasis. Genetic defects that lead to loss of expression or function of KATP channels are the most common cause of HI (KATP-HI). Much progress has been made in our understanding of the molecular genetics and pathophysiology of KATP-HI in the past decades; however, treatment remains challenging, in particular for patients with diffuse disease who do not respond to the KATP channel activator diazoxide. In this review, we discuss current approaches and limitations on the diagnosis and treatment of KATP-HI, and offer perspectives on alternative therapeutic strategies.

De Novo Assembly and Characterization of the Transcriptome of an Omnivorous Camel Cricket (Tachycines meditationis)

Tachycines meditationis (Orthoptera: Rhaphidophoridae: Tachycines) is a widely distributed insect in eastern Asia. This species is common in urban environments, and its unique omnivorous diet may contribute to its success in various habitats. However, molecular studies on the species are scarce. Here, we obtained the first transcriptome sequence of T. meditationis and performed preliminary analyses to test whether the evolution of coding sequences fits the expectations based on the species' ecology. We retrieved 476,495 effective transcripts and annotated 46,593 coding sequences (CDS). We analysed the codon usage and found that directional mutation pressure was the leading cause of codon usage bias in this species. This genome-wide relaxed codon usage pattern in T. meditationis is surprising, given the potentially large population size of this species. Moreover, despite the omnivorous diet, the chemosensory genes of this species do not exhibit codon usage deviating significantly from the genome-level pattern. They also do not seem to experience more gene family expansion than other cave cricket species do. A thorough search for rapidly evolved genes using the dN/dS value showed that genes associated with substance synthesis and metabolic pathways, such as retinol metabolism, aminoacyl-tRNA biosynthesis, and fatty acid metabolism, underwent species-specific positive selection. While some results seem to contradict the species ecology, our transcriptome assembly provides a valuable molecular resource for future studies on camel cricket evolution and molecular genetics for feeding ecology in insects, in general.

Insights into Synonymous Codon Usage Bias in Hepatitis C Virus and Its Adaptation to Hosts

Hepatitis C virus (HCV) is enveloped RNA virus, encoding for a polyprotein that is processed by cellular proteases. The virus is responsible for liver cirrhosis, allograft rejection, and human hepatocellular carcinoma. Based on studies including compositional analysis, odds ratio analysis, parity analysis, skew analysis, relative synonymous codon usage, codon bias, and protein properties, it was evident that codon usage bias in HCV is dependent upon the nucleotide composition. Codon context analysis revealed CTC-CTG as a preferred codon pair. While CGA and CGT codons were rare, none of the codons were rare in HCV-like viruses envisaged in the present study. Many of the preferred codon pairs were valine amino acid-initiated, which possibly infers viral infectivity; hence the role of selection forces appears to act on the HCV genome, which was further validated by neutrality analysis where selection accounted for 87.28%, while mutation accounted for 12.72% force shaping codon usage. Furthermore, codon usage was correlated with the length of the genome. HCV viruses prefer valine-initiated codon pairs, while HCV-like viruses prefer alanine-initiated codon pairs. The HCV host range is very narrow and is confined to only humans and chimpanzees. Based on indices including codon usage correlation analysis, similarity index, and relative codon deoptimization index, it is evident in the study that the chimpanzee is the primary host of the virus. The present study helped elucidate the preferred host for HCV. The information presented in the study paved the way for generating an attenuated vaccine candidate through viral recoding, with finely tuned nucleotide composition and a perfect balance of preferred and rare codons.

Implementing computational methods in tandem with synonymous gene recoding for therapeutic development

Synonymous gene recoding, the substitution of synonymous variants into the genetic sequence, has been used to overcome many production limitations in therapeutic development. However, the safety and efficacy of recoded therapeutics can be difficult to evaluate because synonymous codon substitutions can result in subtle, yet impactful changes in protein features and require sensitive methods for detection. Given that computational approaches have made significant leaps in recent years, we propose that machine-learning (ML) tools may be leveraged to assess gene-recoded therapeutics and foresee an opportunity to adapt codon contexts to enhance some powerful existing tools. Here, we examine how synonymous gene recoding has been used to address challenges in therapeutic development, explain the biological mechanisms underlying its effects, and explore the application of computational platforms to improve the surveillance of functional variants in therapeutic design.

Deep mutational scanning of essential bacterial proteins can guide antibiotic development

Deep mutational scanning is a powerful approach to investigate a wide variety of research questions including protein function and stability. Here, we perform deep mutational scanning on three essential E. coli proteins (FabZ, LpxC and MurA) involved in cell envelope synthesis using high-throughput CRISPR genome editing, and study the effect of the mutations in their original genomic context. We use more than 17,000 variants of the proteins to interrogate protein function and the importance of individual amino acids in supporting viability. Additionally, we exploit these libraries to study resistance development against antimicrobial compounds that target the selected proteins. Among the three proteins studied, MurA seems to be the superior antimicrobial target due to its low mutational flexibility, which decreases the chance of acquiring resistance-conferring mutations that simultaneously preserve MurA function. Additionally, we rank anti-LpxC lead compounds for further development, guided by the number of resistance-conferring mutations against each compound. Our results show that deep mutational scanning studies can be used to guide drug development, which we hope will contribute towards the development of novel antimicrobial therapies.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown.

Methods

In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains.

Results

The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species.

Discussion

This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

YMLA: A comparative platform to carry out functional enrichment analysis for multiple gene lists in yeast

Comparative analysis among multiple gene lists on their functional features is now a routine task due to the advancement of high-throughput experiments. Several enrichment analysis tools were developed in the past. However, these tools mainly focus on one gene list and contain only gene ontology or interaction features. What makes it worse, comparative investigation and customized feature set reanalysis are still unavailable. Therefore, we constructed the YMLA (Yeast Multiple List Analyzer) platform in this research. YMLA includes 39 yeast features and facilitates comparative analysis among multiple gene lists via tabular views, heatmaps, and network plots. Moreover, the customized feature set reanalysis function was implemented in YMLA to help form mechanism hypotheses based on a selected enriched feature subset. We demonstrated the biological applicability of YMLA via example lists consisting of genes with top/bottom translation efficiency values. The analysis results provided by YMLA reveal novel facts consistent with previous experiments. YMLA is available at https://cosbi7.ee.ncku.edu.tw/YMLA/.