Codon usage is an important determinant of gene expression levels largely through its effects on transcription

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.

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

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

A novel de novo synonymous variant in GREB1L impacts the mRNA splicing associated with aplasia of the urogenital system

GREB1-like retinoic acid receptor coactivator (GREB1L) gene is associated with autosomal dominant renal hypodysplasia/aplasia 3 (RHDA3) and deafness, autosomal dominant 80 (DFNA80). Among the GREB1L variants reported, most of them are missense or frameshift, while no pathogenic synonymous variants have been recorded. Classical theory paid little attention to synonymous variants and classified it as nonpathogenic; however, recent studies suggest that the variants might be equally important. Here, we report a 7-year-old girl with new symptoms of clitoromegaly, uterovaginal, and ovarian agenesis as well as right kidney missing. A novel de novo GREB1L synonymous variant (NM_001142966: c.4731C>T, p.G1577=) was identified via whole exome sequencing. The variant was predicted to be disease-causing through in silico analysis and was classified as likely pathogenic. Minigene splicing assays confirmed a 6 bp deletion in mutant cDNA comparing with the wild type, leading to two amino acids lost in GREB1L protein. Secondary and tertiary structure modeling showed alterations in protein structure. Our finding reveals a novel GREB1L variant with a new phenotype of urogenital system and is the first to report a pathogenic synonymous variant in GREB1L which affects mRNA splicing, suggesting synonymous variants cannot be ignored in prenatal diagnosis and genetic counseling.

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.

smAMPsTK: a toolkit to unravel the smORFome encoding AMPs of plant species

The pervasive repertoire of plant molecules with the potential to serve as a substitute for conventional antibiotics has led to obtaining better insights into plant-derived antimicrobial peptides (AMPs). The massive distribution of Small Open Reading Frames (smORFs) throughout eukaryotic genomes with proven extensive biological functions reflects their practicality as antimicrobials. Here, we have developed a pipeline named smAMPsTK to unveil the underlying hidden smORFs encoding AMPs for plant species. By applying this pipeline, we have elicited AMPs of various functional activity of lengths ranging from 5 to 100 aa by employing publicly available transcriptome data of five different angiosperms. Later, we studied the coding potential of AMPs-smORFs, the inclusion of diverse translation initiation start codons, and amino acid frequency. Codon usage study signifies no such codon usage biases for smORFs encoding AMPs. Majorly three start codons are prominent in generating AMPs. The evolutionary and conservational study proclaimed the widespread distribution of AMPs encoding genes throughout the plant kingdom. Domain analysis revealed that nearly all AMPs have chitin-binding ability, establishing their role as antifungal agents. The current study includes a developed methodology to characterize smORFs encoding AMPs, and their implications as antimicrobial, antibacterial, antifungal, or antiviral provided by SVM score and prediction status calculated by machine learning-based prediction models. The pipeline, complete package, and the results derived for five angiosperms are freely available at https://github.com/skbinfo/smAMPsTK.Communicated by Ramaswamy H. Sarma.

Analysis of Codon Usage Bias in Chloroplast Genomes of Dryas octopetala var. asiatica (Rosaceae)

Dryas octopetala var. asiatica, a dwarf shrub belonging to the Rosaceae family and native to Asia, exhibits notable plasticity in photosynthesis in response to temperature variations. However, the codon usage patterns and factors influencing them in the chloroplast genome of this species have not yet been documented. This study sequenced and assembled the complete genome of D. octopetala var. asiatica. The annotated genes in the chloroplast genome were analyzed for codon composition through multivariate statistical methods including a neutrality plot, a parity rule 2 (PR2) bias plot, and an effective number of codons (ENC) plot using CodonW 1.4.2 software. The results indicated that the mean GC content of 53 CDSs was 38.08%, with the average GC content at the third codon base position being 27.80%, suggesting a preference for A/U(T) at the third codon position in chloroplast genes. Additionally, the chloroplast genes exhibited a weak overall codon usage bias (CUB) based on ENC values and other indicators. Correlation analysis showed a significant negative correlation between ENC value and GC2, an extremely positive correlation with GC3, but no correlation with GC1 content. These findings highlight the importance of the codon composition at the third position in influencing codon usage bias. Furthermore, our analysis indicated that the CUB of the chloroplast genome of D. octopetala var. asiatica was primarily influenced by natural selection and other factors. Finally, this study identified UCA, CCU, GCU, AAU, GAU, and GGU as the optimal codons. These results offer a foundational understanding for genetic modification and evolutionary dynamics of the chloroplast genome of D. octopetala var. asiatica.

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.

Selection on synonymous sites: the unwanted transcript hypothesis

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Alter codon bias of the P. pastoris genome to overcome a bottleneck in codon optimization strategy development and improve protein expression

Apart from its role in translation, codon bias is also an important mechanism to regulate mRNA levels. The traditional frequency-based codon optimization strategy is rather efficient in organisms such as N. crassa, but much less in yeast P. pastoris which is a popular host for heterologous protein expression. This is because that unlike N. crassa, the preferred codons of P. pastoris are actually AU-rich and hence codon optimization for extremely low GC content comes with issues of pre-mature transcriptional termination or low RNA stability in spite of translational advantages. To overcome this bottleneck, we focused on three reporter genes in P. pastoris first and confirmed the great advantage of GC-prone codon optimization on mRNA levels. Then we altered the codon bias profile of P. pastoris by introducing additional rare tRNA gene copies. Prior to that we constructed IPTG-regulated tRNA species to enable chassis cells to switch between different codon bias status. As demonstrated again with reporter genes, protein yield of luc and 0788 was successfully increased by 4-5 folds in chassis cells. In summary, here we provide an alternative codon optimization strategy for genes with unsatisfactory performance under traditional codon frequency-based optimization.

Chloroplast genome analysis and evolutionary insights in the versatile medicinal plant Calendula officinalis L

Calendula officinalis L.is a versatile medicinal plant with numerous applications in various fields. However, its chloroplast genome structure, features, phylogeny, and patterns of evolution and mutation remain largely unexplored. This study examines the chloroplast genome, phylogeny, codon usage bias, and divergence time of C. officinalis, enhancing our understanding of its evolution and adaptation. The chloroplast genome of C. officinalis is a 150,465 bp circular molecule with a G + C content of 37.75% and comprises 131 genes. Phylogenetic analysis revealed a close relationship between C. officinalis, C. arvensis, and Osteospermum ecklonis. A key finding is the similarity in codon usage bias among these species, which, coupled with the divergence time analysis, supports their close phylogenetic proximity. This similarity in codon preference and divergence times underscores a parallel evolutionary adaptation journey for these species, highlighting the intricate interplay between genetic evolution and environmental adaptation in the Asteraceae family. Moreover unique evolutionary features in C. officinalis, possibly associated with certain genes were identified, laying a foundation for future research into the genetic diversity and medicinal value of C. officinalis.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Designing of a multi-epitopes based vaccine against Haemophilius parainfluenzae and its validation through integrated computational approaches

Haemophilus parainfluenzae is a Gram-negative opportunist pathogen within the mucus of the nose and mouth without significant symptoms and has an ability to cause various infections ranging from ear, eye, and sinus to pneumonia. A concerning development is the increasing resistance of H. parainfluenzae to beta-lactam antibiotics, with the potential to cause dental infections or abscesses. The principal objective of this investigation is to utilize bioinformatics and immuno-informatic methodologies in the development of a candidate multi-epitope Vaccine. The investigation focuses on identifying potential epitopes for both B cells (B lymphocytes) and T cells (helper T lymphocytes and cytotoxic T lymphocytes) based on high non-toxic and non-allergenic characteristics. The selection process involves identifying human leukocyte antigen alleles demonstrating strong associations with recognized antigenic and overlapping epitopes. Notably, the chosen alleles aim to provide coverage for 90% of the global population. Multi-epitope constructs were designed by using suitable linker sequences. To enhance the immunological potential, an adjuvant sequence was incorporated using the EAAAK linker. The final vaccine construct, comprising 344 amino acids, was achieved after the addition of adjuvants and linkers. This multi-epitope Vaccine demonstrates notable antigenicity and possesses favorable physiochemical characteristics. The three-dimensional conformation underwent modeling and refinement, validated through in-silico methods. Additionally, a protein-protein molecular docking analysis was conducted to predict effective binding poses between the multi-epitope Vaccine and the Toll-like receptor 4 protein. The Molecular Dynamics (MD) investigation of the docked TLR4-vaccine complex demonstrated consistent stability over the simulation period, primarily attributed to electrostatic energy. The docked complex displayed minimal deformation and enhanced rigidity in the motion of residues during the dynamic simulation. Furthermore, codon translational optimization and computational cloning was performed to ensure the reliability and proper expression of the multi-Epitope Vaccine. It is crucial to emphasize that despite these computational validations, experimental research in the laboratory is imperative to demonstrate the immunogenicity and protective efficacy of the developed vaccine. This would involve practical assessments to ascertain the real-world effectiveness of the multi-epitope Vaccine.

The complete Chloroplast genome of Stachys geobombycis and comparative analysis with related Stachys species

Herb genomics, at the forefront of traditional Chinese medicine research, combines genomics with traditional practices, facilitating the scientific validation of ancient remedies. This integration enhances public understanding of traditional Chinese medicine's efficacy and broadens its scope in modern healthcare. Stachys species encompass annual or perennial herbs or small shrubs, exhibiting simple petiolate or sessile leaves. Despite their wide-ranging applications across various fields, molecular data have been lacking, hindering the precise identification and taxonomic elucidation of Stachys species. To address this gap, we assembled the complete chloroplast (CP) genome of Stachys geobombycis and conducted reannotation and comparative analysis of seven additional species within the Stachys genus. The findings demonstrate that the CP genomes of these species exhibit quadripartite structures, with lengths ranging from 14,523 to 150,599 bp. Overall, the genome structure remains relatively conserved, hosting 131 annotated genes, including 87 protein coding genes, 36 tRNA genes, and 8 rRNA genes. Additionally, 78 to 98 SSRs and long repeat sequences were detected , and notably, 6 highly variable regions were identified as potential molecular markers in the CP genome through sequence alignment. Phylogenetic analysis based on Bayesian inference and maximum likelihood methods strongly supported the phylogenetic position of the genus Stachys as a member of Stachydeae tribe. Overall, this comprehensive bioinformatics study of Stachys CP genomes lays the groundwork for phylogenetic classification, plant identification, genetic engineering, evolutionary studies, and breeding research concerning medicinal plants within the Stachys genus.

Plastid genome data provide new insights into the dynamic evolution of the tribe Ampelopsideae (Vitaceae)

BACKGROUND

Ampelopsideae J. Wen & Z.L. Nie is a small-sized tribe of Vitaceae Juss., including ca. 47 species from four genera showing a disjunct distribution worldwide across all the continents except Antarctica. There are numerous species from the tribe that are commonly used as medicinal plants with immune-modulating, antimicrobial, and anti-hypertensive properties. The tribe is usually recognized into three clades, i.e., Ampelopsis Michx., Nekemias Raf., and the Southern Hemisphere clade. However, the relationships of the three clades differ greatly between the nuclear and the plastid topologies. There has been limited exploration of the chloroplast phylogenetic relationships within Ampelopsideae, and studies on the chloroplast genome structure of this tribe are only available for a few individuals. In this study, we aimed to investigate the evolutionary characteristics of plastid genomes of the tribe, including their genome structure and evolutionary insights.

RESULTS

We sequenced, assembled, and annotated plastid genomes of 36 species from the tribe and related taxa in the family. Three main clades were recognized within Ampelopsideae, corresponding to Ampelopsis, Nekemias, and the Southern Hemisphere lineage, respectively, and all with 100% bootstrap supports. The genome sequences and content of the tribe are highly conserved. However, comparative analyses suggested that the plastomes of Nekemias demonstrate a contraction in the large single copy region and an expansion in the inverted repeat region, and possess a high number of forward and palindromic repeat sequences distinct from both Ampelopsis and the Southern Hemisphere taxa.

CONCLUSIONS

Our results highlighted plastome variations in genome length, expansion or contraction of the inverted repeat region, codon usage bias, and repeat sequences, are corresponding to the three lineages of the tribe, which probably faced with different environmental selection pressures and evolutionary history. This study provides valuable insights into understanding the evolutionary patterns of plastid genomes within the Ampelopsideae of Vitaceae.

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

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

Specific codons control cellular resources and fitness

As cellular engineering progresses from simply overexpressing proteins to imparting complex phenotypes through multigene expression, judicious appropriation of cellular resources is essential. Since codon use is degenerate and biased, codons may control cellular resources at a translational level. We investigate how partitioning transfer RNA (tRNA) resources by incorporating dissimilar codon usage can drastically alter interdependence of expression level and burden on the host. By isolating the effect of individual codons' use during translation elongation while eliminating confounding factors, we show that codon choice can trans-regulate fitness of the host and expression of other heterologous or native genes. We correlate specific codon usage patterns with host fitness and derive a coding scheme for multigene expression called the Codon Health Index (CHI, χ). This empirically derived coding scheme (χ) enables the design of multigene expression systems that avoid catastrophic cellular burden and is robust across several proteins and conditions.

The Genome of Plasmodium gonderi: Insights into the Evolution of Human Malaria Parasites

Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Can the Revolution in mRNA-Based Vaccine Technologies Solve the Intractable Health Issues of Current Ruminant Production Systems?

To achieve the World Health Organization's global Sustainable Development Goals, increased production of high-quality protein for human consumption is required while minimizing, ideally reducing, environmental impacts. One way to achieve these goals is to address losses within current livestock production systems. Infectious diseases are key limiters of edible protein production, affecting both quantity and quality. In addition, some of these diseases are zoonotic threats and potential contributors to the emergence of antimicrobial resistance. Vaccination has proven to be highly successful in controlling and even eliminating several livestock diseases of economic importance. However, many livestock diseases, both existing and emerging, have proven to be recalcitrant targets for conventional vaccination technologies. The threat posed by the COVID-19 pandemic resulted in unprecedented global investment in vaccine technologies to accelerate the development of safe and efficacious vaccines. While several vaccination platforms emerged as front runners to meet this challenge, the clear winner is mRNA-based vaccination. The challenge now is for livestock industries and relevant stakeholders to harness these rapid advances in vaccination to address key diseases affecting livestock production. This review examines the key features of mRNA vaccines, as this technology has the potential to control infectious diseases of importance to livestock production that have proven otherwise difficult to control using conventional approaches. This review focuses on the challenging diseases of ruminants due to their importance in global protein production. Overall, the current literature suggests that, while mRNA vaccines have the potential to address challenges in veterinary medicine, further developments are likely to be required for this promise to be realized for ruminant and other livestock species.

Molecular Mechanisms and the Significance of Synonymous Mutations

Synonymous mutations result from the degeneracy of the genetic code. Most amino acids are encoded by two or more codons, and mutations that change a codon to another synonymous codon do not change the amino acid in the gene product. Historically, such mutations have been considered silent because they were assumed to have no to very little impact. However, research in the last few decades has produced several examples where synonymous mutations play important roles. These include optimizing expression by enhancing translation initiation and accelerating or decelerating translation elongation via codon usage and mRNA secondary structures, stabilizing mRNA molecules and preventing their breakdown before translation, and faulty protein folding or increased degradation due to enhanced ubiquitination and suboptimal secretion of proteins into the appropriate cell compartments. Some consequences of synonymous mutations, such as mRNA stability, can lead to different outcomes in prokaryotes and eukaryotes. Despite these examples, the significance of synonymous mutations in evolution and in causing disease in comparison to nonsynonymous mutations that do change amino acid residues in proteins remains controversial. Whether the molecular mechanisms described by which synonymous mutations affect organisms can be generalized remains poorly understood and warrants future research in this area.

A codon-based live-cell biomonitoring system for assessing intracellular phenylalanine bioavailability in cyanobacteria

Phenylalanine, as an essential aromatic amino acid, is not only needed for protein and vital molecules such as neurotransmitter and hormone synthesis but also a substrate for the biosynthesis of phenylpropanoids and various bioactive compounds. The metabolism of phenylalanine is dynamic and transitory, which would otherwise inhibit cell growth. Therefore, it is challenging and imperative to monitor intracellular phenylalanine bioavailability in real time, which has great significance for evaluating the effectiveness of introducing pathway-specific genetic modifications to enhance phenylalanine generation. In this study, we proposed a live-cell biomonitoring system to assess phenylalanine bioavailability in real time in cyanobacteria based on codon degeneracy and species-specific usage bias. The biomonitoring system was generated through genetic modification of phenylalanine codons in the chloramphenicol antibiotic resistance gene to wholly preferred and rare codons, in combination with an orthogonal constitutive promoter Trc to express these genes. Cyanobacterial cells equipped with a preferred codon-based gene showed a significant growth advantage over those with rare codons under antibiotic pressure, while the delayed growth caused by rare codon-based genes could be rescued by supplementing phenylalanine in the cultivation medium. Increasing intracellular phenylalanine bioavailability could promote rare codon-based gene containing cell growth to a similar level as wild-type strains harboring preferred codon-based gene, providing a live-cell visualized screening method to relatively define phenylalanine content from either random mutation libraries or pathway-specific engineering cyanobacterial chassis before conducting labor-intensive quantitative measurements.

Expanding the synthetic biology toolbox of Cupriavidus necator for establishing fatty acid production

The Gram-negative betaproteobacterium Cupriavidus necator is a chemolithotroph that can convert carbon dioxide into biomass. Cupriavidus necator has been engineered to produce a variety of high-value chemicals in the past. However, there is still a lack of a well-characterized toolbox for gene expression and genome engineering. Development and optimization of biosynthetic pathways in metabolically engineered microorganisms necessitates control of gene expression via functional genetic elements such as promoters, ribosome binding sites (RBSs), and codon optimization. In this work, a set of inducible and constitutive promoters were validated and characterized in C. necator, and a library of RBSs was designed and tested to show a 50-fold range of expression for green fluorescent protein (gfp). The effect of codon optimization on gene expression in C. necator was studied by expressing gfp and mCherry genes with varied codon-adaptation indices and was validated by expressing codon-optimized variants of a C12-specific fatty acid thioesterase to produce dodecanoic acid. We discuss further hurdles that will need to be overcome for C. necator to be widely used for biosynthetic processes.

Transcriptional landscape of the pMP7017 megaplasmid and its impact on the Bifidobacterium breve UCC2003 transcriptome

The 190 kb megaplasmid pMP7017 of Bifidobacterium breve JCM7017 represents the first conjugative and largest plasmid characterised within this genus to date. In the current study, we adopted an integrated approach combining transcriptomics, whole genome comparative analysis and metagenomic data mining to understand the biology of pMP7017 and related megaplasmids, and to assess the impact of plasmid-carriage on the host strain. The data generated revealed variations within basic features of promoter elements which correlate with a high level of transcription on the plasmid and highlight the transcriptional activity of genes encoding both offensive and defensive adaptations, including a Type IIL restriction-modification system, an anti-restriction system and four Type II toxin-antitoxin systems. Furthermore, a highly transcribed tmRNA, which likely provides translational support to the host strain, was identified, making pMP7017 the first plasmid of the Bifidobacterium genus and the smallest plasmid known to express a tmRNA. Analyses of synteny and variability among pMP7017 and related plasmids indicate substantial diversity in gene organisation and accessory gene cargo highlighting diverse (co-)evolution and potential host-specific rearrangements and adaptations. Systematic analysis of the codon usage profile of transcriptionally active pMP7017-encoded genes suggests that pMP7017 originated from (sub)species of Bifidobacterium longum. Furthermore, mining of metagenomic data suggests the presence of pMP7017-homologues in ~10% of microbiome samples, mostly infants and/or mothers from various geographical locations. Comparative transcriptome analysis of the B. breve UCC2003 chromosome in the presence or absence of pMP7017 revealed differential expression of genes representing 8% of the total gene pool. Genes involved in genetic information processing were exclusively upregulated, while altered expression of genes involved in biofilm production and polysaccharide biosynthesis was also observed.

Synonymous codon usage regulates translation initiation

Nonoptimal synonymous codons repress gene expression, but the underlying mechanisms are poorly understood. We and others have previously shown that nonoptimal codons slow translation elongation speeds and thereby trigger messenger RNA (mRNA) degradation. Nevertheless, transcript levels are often insufficient to explain protein levels, suggesting additional mechanisms by which codon usage regulates gene expression. Using reporters in human and Drosophila cells, we find that transcript levels account for less than half of the variation in protein abundance due to codon usage. This discrepancy is explained by translational differences whereby nonoptimal codons repress translation initiation. Nonoptimal transcripts are also less bound by the translation initiation factors eIF4E and eIF4G1, providing a mechanistic explanation for their reduced initiation rates. Importantly, translational repression can occur without mRNA decay and deadenylation, and it does not depend on the known nonoptimality sensor, CNOT3. Our results reveal a potent mechanism of regulation by codon usage where nonoptimal codons repress further rounds of translation.

Decoding stoichiometric protein synthesis in E. coli through translation rate parameters

E. coli is one of the most widely used organisms for understanding the principles of cellular and molecular genetics. However, we are yet to understand the origin of several experimental observations related to the regulation of gene expression in E. coli. One of the prominent examples in this context is the proportional synthesis in multiprotein complexes where all of their obligate subunits are produced in proportion to their stoichiometry. In this work, by combining the next-generation sequencing data with the stochastic simulations of protein synthesis, we explain the origin of proportional protein synthesis in multicomponent complexes. We find that the estimated initiation rates for the translation of all subunits in those complexes are proportional to their stoichiometry. This constraint on protein synthesis kinetics enforces proportional protein synthesis without requiring any feedback mechanism. We also find that the translation initiation rates in E. coli are influenced by the coding sequence length and the enrichment of A and C nucleotides near the start codon. Thus, this study rationalizes the role of conserved and nonrandom features of genes in regulating the translation kinetics and unravels a key principle of the regulation of protein synthesis.

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

OBJECTIVE

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

METHODOLOGY

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

RESULTS

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

Proteome based analysis of circulating SARS-CoV-2 variants: approach to a universal vaccine candidate

The discovery of the first infectious variant in Wuhan, China, in December 2019, has posed concerns over global health due to the spread of COVID-19 and subsequent variants. While the majority of patients experience flu-like symptoms such as cold and fever, a small percentage, particularly those with compromised immune systems, progress from mild illness to fatality. COVID-19 is caused by a RNA virus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our approach involved utilizing immunoinformatic to identify vaccine candidates with multiple epitopes and ligand-binding regions in reported SARS-CoV-2 variants. Through analysis of the spike glycoprotein, we identified dominant epitopes for T-cells and B-cells, resulting in a vaccine construct containing two helper T-cell epitopes, six cytotoxic T-cell epitopes, and four linear B-cell epitopes. Prior to conjugation with adjuvants and linkers, all epitopes were evaluated for antigenicity, toxicity, and allergenicity. Additionally, we assessed the vaccine Toll-Like Receptors complex (2, 3, and 4). The vaccine construct demonstrated antigenicity, non-toxicity, and non-allergenicity, thereby enabling the host to generate antibodies with favorable physicochemical characteristics. Furthermore, the 3D structure of the B-cell construct exhibited a ProSA-web z-score plot with a value of -1.71, indicating the reliability of the designed structure. The Ramachandran plot analysis revealed that 99.6% of the amino acid residues in the vaccine subunit were located in the high favored observation region, further establishing its strong candidacy as a vaccination option.

Degradation of polypropylene by fungi Coniochaeta hoffmannii and Pleurostoma richardsiae

The urgent need for better disposal and recycling of plastics has motivated a search for microbes with the ability to degrade synthetic polymers. While microbes capable of metabolizing polyurethane and polyethylene terephthalate have been discovered and even leveraged in enzymatic recycling approaches, microbial degradation of additive-free polypropylene (PP) remains elusive. Here we report the isolation and characterization of two fungal strains with the potential to degrade pure PP. Twenty-seven fungal strains, many isolated from hydrocarbon contaminated sites, were screened for degradation of commercially used textile plastic. Of the candidate strains, two identified as Coniochaeta hoffmannii and Pleurostoma richardsiae were found to colonize the plastic fibers using scanning electron microscopy (SEM). Further experiments probing degradation of pure PP films were performed using C. hoffmannii and P. richardsiae and analyzed using SEM, Raman spectroscopy and Fourier transform infrared spectroscopy with attenuated total reflectance (FTIR-ATR). The results showed that the selected fungi were active against pure PP, with distinct differences in the bonds targeted and the degree to which each was altered. Whole genome and transcriptome sequencing was conducted for both strains and the abundance of carbohydrate active enzymes, GC content, and codon usage bias were analyzed in predicted proteomes for each. Enzymatic assays were conducted to assess each strain's ability to degrade naturally occurring compounds as well as synthetic polymers. These investigations revealed potential adaptations to hydrocarbon-rich environments and provide a foundation for further investigation of PP degrading activity in C. hoffmannii and P. richardsiae.

Maximizing Heterologous Expression of Engineered Type I Polyketide Synthases: Investigating Codon Optimization Strategies

Type I polyketide synthases (T1PKSs) hold enormous potential as a rational production platform for the biosynthesis of specialty chemicals. However, despite great progress in this field, the heterologous expression of PKSs remains a major challenge. One of the first measures to improve heterologous gene expression can be codon optimization. Although controversial, choosing the wrong codon optimization strategy can have detrimental effects on the protein and product levels. In this study, we analyzed 11 different codon variants of an engineered T1PKS and investigated in a systematic approach their influence on heterologous expression in Corynebacterium glutamicum, Escherichia coli, and Pseudomonas putida. Our best performing codon variants exhibited a minimum 50-fold increase in PKS protein levels, which also enabled the production of an unnatural polyketide in each of these hosts. Furthermore, we developed a free online tool (https://basebuddy.lbl.gov) that offers transparent and highly customizable codon optimization with up-to-date codon usage tables. In this work, we not only highlight the significance of codon optimization but also establish the groundwork for the high-throughput assembly and characterization of PKS pathways in alternative hosts.

Genetic engineering strategy for generating a stable dsRNA virus vector using a virus-like codon-modified transgene

IMPORTANCE

The stabilities of transgenes in RNA virus vectors differ between the genes of interest, but the molecular mechanisms determining genetic stability remain unknown. This study demonstrated that the stability of a transgene was affected by the nucleotide composition, and altering the codon usage of transgenes to resemble that of the viral genome significantly increased transgene stability in double-stranded RNA virus vectors. The virus-like codon modification strategy enabled generation of stable rotavirus and mammalian orthoreovirus vectors, which could be developed as machinery for gene delivery to the intestines and/or respiratory organs. This technology has further potential to be expanded to other RNA viruses.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Prevalence and genomic analysis of t203-like G9 (G9-VI) rotaviruses circulating in children with gastroenteritis in Beijing, China

Our previous surveillance revealed that t203-like G9 (tentatively designated subtype G9-VI) rotaviruses re-emerged in 2010 in Beijing and rapidly prevailed over the G9-III subtype (the most common G9 subtype globally) and previously predominant G genotypes over the following two years. G9-VI belongs to the VP7 evolutionary lineage VI, which includes unusual and sporadic human rotaviruses from China (t203) and Japan. To obtain insight into the epidemiology, evolution, and transmission advantages of G9-VI rotavirus, we performed follow-up surveillance (2014-2017) and whole-genome analysis of 12 representative G9 strains. The results showed that the G9 genotype was predominant (77.4%), with a marked increase in prevalence (previously 43.5%). Within the G9 genotype, subtype G9-VI accounted for the majority (98.3%) of cases. The most prevalent P-genotype was P[8] (93.7%), within which subtype P[8]b was rare (0.7%). Phylogenetically, the G9-VI subtype strains in this study clustered closely with contemporary emerging human rotaviruses from many other countries in VP7 lineage VI, indicating that this subtype is capable of spreading globally. These currently emerging G9-VI rotaviruses formed a distinct monophyletic subcluster when compared to early G9-VI rotaviruses. Furthermore, four specific amino acid substitutions and synonymous codon substitutions were observed in the VP7 genes between the current G9-VI and globally common G9-III rotaviruses. The remaining nine genes of all of the analyzed representative G9 strains, whether G9-VI or G9-III, combined with the P[8]a, P[8]b, or P[6] genotype and exhibited the same Wa-like backbone constellation.

Compositional shifts associated with major evolutionary transitions in plants

Heterogeneity in gene trees, morphological characters, and composition has been associated with several major plant clades. Here, we examine heterogeneity in composition across a large transcriptomic dataset of plants to better understand whether locations of shifts in composition are shared across gene regions and whether directions of shifts within clades are shared across gene regions. We estimate mixed models of composition for both nucleotide and amino acids across a recent large-scale transcriptomic dataset for plants. We find shifts in composition across both nucleotide and amino acid datasets, with more shifts detected in nucleotides. We find that Chlorophytes and lineages within experience the most shifts. However, many shifts occur at the origins of land, vascular, and seed plants. While genes in these clades do not typically share the same composition, they tend to shift in the same direction. We discuss potential causes of these patterns. Compositional heterogeneity has been highlighted as a potential problem for phylogenetic analysis, but the variation presented here highlights the need to further investigate these patterns for the signal of biological processes.

Growth-dependent Gene Expression Variation Influences the Strength of Codon Usage Biases

The most highly expressed genes in microbial genomes tend to use a limited set of synonymous codons, often referred to as "preferred codons." The existence of preferred codons is commonly attributed to selection pressures on various aspects of protein translation including accuracy and/or speed. However, gene expression is condition-dependent and even within single-celled organisms transcript and protein abundances can vary depending on a variety of environmental and other factors. Here, we show that growth rate-dependent expression variation is an important constraint that significantly influences the evolution of gene sequences. Using large-scale transcriptomic and proteomic data sets in Escherichia coli and Saccharomyces cerevisiae, we confirm that codon usage biases are strongly associated with gene expression but highlight that this relationship is most pronounced when gene expression measurements are taken during rapid growth conditions. Specifically, genes whose relative expression increases during periods of rapid growth have stronger codon usage biases than comparably expressed genes whose expression decreases during rapid growth conditions. These findings highlight that gene expression measured in any particular condition tells only part of the story regarding the forces shaping the evolution of microbial gene sequences. More generally, our results imply that microbial physiology during rapid growth is critical for explaining long-term translational constraints.

Antimicrobial peptides in Dendrobium officinale: Genomic parameters, peptide structures, and gene expression patterns

A weak codon usage bias was found in Dendrobium catenatum (D. officiale) antimicrobial peptides (AMPs), after the analysis of relative synonymous codon usage, GC contents, and the effective number of codons. The codon usage preference was mainly influenced by natural selection pressure. The self-optimized prediction method and SWISS-MODEL were applied for peptide structural and domain analyses, and some typical antimicrobial domains were found in D. officinale AMP amino sequences, such as knot1 domain, gibberellins-stimulated domain, cupin_1 domain, defensin_like domain, and SLR1-BP (S locus-related glycoprotein 1 binding pollen coat protein) domain. To investigate the AMPs gene expression pattern, abiotic stresses, such as salt stress, drought stress, salicylic acid (SA), and methyl jasmonate (JA), were applied and the gene expression levels were detected by the real-time fluorescent quantitative polymerase chain reaction. Results showed that, even though the basic AMPs gene expressions were low, some AMPs can still be induced by salt dress, while the drought dress did not show the same impact. The SA and JA signaling pathways might be involved in most of the AMPs expressions. The natural selection of the D. officinale AMPs and thus forming diverse types of AMPs enhanced the plant's innate immunity and disease resistance capability, which would lead to a better understanding of the molecular mechanism for D. officinale adapting to the environment. The finding that salt stress, SA, and JA signaling pathways can induce AMP expression lays a foundation for the further development and functional verification of D. officinale AMPs.

Comparative analysis of the medicinal plant Polygonatum kingianum (Asparagaceae) with related verticillate leaf types of the Polygonatum species based on chloroplast genomes

Background

Polygonatum kingianum has been widely used as a traditional Chinese medicine as well as a healthy food. Because of its highly variable morphology, this medicinal plant is often difficult to distinguish from other related verticillate leaf types of the Polygonatum species. The contaminants in P. kingianum products not only decrease the products' quality but also threaten consumer safety, seriously inhibiting the industrial application of P. kingianum.

Methods

Nine complete chloroplast (cp) genomes of six verticillate leaf types of the Polygonatum species were de novo assembled and systematically analyzed.

Results

The total lengths of newly sequenced cp genomes ranged from 155,437 to 155,977 bp, including 86/87 protein-coding, 38 tRNA, and 8 rRNA genes, which all exhibited well-conserved genomic structures and gene orders. The differences in the IR/SC (inverted repeats/single-copy) boundary regions and simple sequence repeats were detected among the verticillate leaf types of the Polygonatum cp genomes. Comparative cp genomes analyses revealed that a higher similarity was conserved in the IR regions than in the SC regions. In addition, 11 divergent hotspot regions were selected, providing potential molecular markers for the identification of the Polygonatum species with verticillate leaf types. Phylogenetic analysis indicated that, as a super barcode, plastids realized a fast and efficient identification that clearly characterized the relationships within the verticillate leaf types of the Polygonatum species. In brief, our results not only enrich the data on the cp genomes of the genus Polygonatum but also provide references for the P. kingianum germplasm resource protection, herbal cultivation, and drug production.

Conclusion

This study not only accurately identifies P. kingianum species, but also provides valuable information for the development of molecular markers and phylogenetic analyses of the Polygonatum species with verticillate leaf types.

Optimising the zebrafish Cre/Lox toolbox. Codon improved iCre, new gateway tools, Cre protein and guidelines

We recently introduced the Cre/Lox technology in our laboratory for both transient (mRNA injections) and stable/transgenic experiments. We experienced significant issues such as silencing, mosaicism, and partial recombination using both approaches. Reviewing the literature gave us the impression that these issues are common among the zebrafish community using the Cre/Lox system. While some researchers took advantage of these problems for specific applications, such as cell and lineage tracing using the Zebrabow construct, we tried here to improve the efficiency and reliability of this system by constituting and testing a new set of tools for zebrafish genetics. First, we implemented a codon-improved Cre version (iCre) designed for rodent studies to counteract some of the aforementioned problems. This eukaryotic-like iCre version was engineered to i) reduce silencing, ii) increase mRNA stability, iii) enhance translational efficiency, and iv) improve nuclear translocation. Second, we established a new set of tol2-kit compatible vectors to facilitate the generation of either iCre-mRNA or iCre-transgenes for transient and transgenic experiments, respectively. We then validated the use of this material and are providing tips for users. Interestingly, during the validation steps, we found that maternal iCRE-mRNA and/or protein deposition from female transgenics systematically led to complete/homogeneous conversion of all tested Lox-responder-transgenes, as opposed to some residual imperfect conversion when using males-drivers or mRNA injections. Considering that we did not find any evidence of Cre-protein soaking and injections in the literature as it is usually conducted with cells, we tested these approaches. While soaking of cell-permeant CRE-protein did not lead to any detectable Lox-conversion, 1ng-10 ng protein injections led to robust and homogeneous Lox-recombination, suggesting that the use of protein could be a robust option for exogenous delivery. This approach may be particularly useful to manipulate housekeeping genes involved in development, sex determination and reproduction which are difficult to investigate with traditional knockout approaches. All in all, we are providing here a new set of tools that should be useful in the field.

Base Composition, Codon Usage, and Patterns of Gene Sequence Evolution in Butterflies

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage, and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC, and selection on codon usage.

Single nucleotide polymorphisms associated with wine fermentation and adaptation to nitrogen limitation in wild and domesticated yeast strains

For more than 20 years, Saccharomyces cerevisiae has served as a model organism for genetic studies and molecular biology, as well as a platform for biotechnology (e.g., wine production). One of the important ecological niches of this yeast that has been extensively studied is wine fermentation, a complex microbiological process in which S. cerevisiae faces various stresses such as limited availability of nitrogen. Nitrogen deficiencies in grape juice impair fermentation rate and yeast biomass production, leading to sluggish or stuck fermentations, resulting in considerable economic losses for the wine industry. In the present work, we took advantage of the "1002 Yeast Genomes Project" population, the most complete catalogue of the genetic variation in the species and a powerful resource for genotype-phenotype correlations, to study the adaptation to nitrogen limitation in wild and domesticated yeast strains in the context of wine fermentation. We found that wild and domesticated yeast strains have different adaptations to nitrogen limitation, corroborating their different evolutionary trajectories. Using a combination of state-of-the-art bioinformatic (GWAS) and molecular biology (CRISPR-Cas9) methodologies, we validated that PNP1, RRT5 and PDR12 are implicated in wine fermentation, where RRT5 and PDR12 are also involved in yeast adaptation to nitrogen limitation. In addition, we validated SNPs in these genes leading to differences in fermentative capacities and adaptation to nitrogen limitation. Altogether, the mapped genetic variants have potential applications for the genetic improvement of industrial yeast strains.

Growth-dependent gene expression variation influences the strength of codon usage biases

The most highly expressed genes in microbial genomes tend to use a limited set of synonymous codons, often referred to as "preferred codons." The existence of preferred codons is commonly attributed to selection pressures on various aspects of protein translation including accuracy and/or speed. However, gene expression is condition-dependent and even within single-celled organisms transcript and protein abundances can vary depending on a variety of environmental and other factors. Here, we show that growth rate-dependent expression variation is an important constraint that significantly influences the evolution of gene sequences. Using large-scale transcriptomic and proteomic data sets in Escherichia coli and Saccharomyces cerevisiae, we confirm that codon usage biases are strongly associated with gene expression but highlight that this relationship is most pronounced when gene expression measurements are taken during rapid growth conditions. Specifically, genes whose relative expression increases during periods of rapid growth have stronger codon usage biases than comparably expressed genes whose expression decreases during rapid growth conditions. These findings highlight that gene expression measured in any particular condition tells only part of the story regarding the forces shaping the evolution of microbial gene sequences. More generally, our results imply that microbial physiology during rapid growth is critical for explaining long-term translational constraints.

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

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

Development of a Streptomyces-based system for facile thioholgamide library generation and analysis

Derivatizing natural products (NPs) is essential in structure-activity relationship (SAR) studies, compound optimization, and drug development. Ribosomally synthesized and post-translationally modified peptides (RiPPs) represent one of the major classes of natural products. Thioholgamide represents thioamitide - a recently emerged family of RiPPs with unique structures and great potential in anticancer drug development. Although the method for generating the RiPP library by codon substitutions in the precursor peptide gene is straightforward, the techniques to perform RiPP derivatization in Actinobacteria remain limited and time-consuming. Here, we report a facile system for producing a library of randomized thioholgamide derivatives utilizing an optimized Streptomyces host. This technique enabled us to access all possible amino acid substitutions of the thioholgamide molecule, one position at a time. Out of 152 potential derivatives, 85 were successfully detected, revealing the impact of amino acid substitutions on thioholgamide post-translational modifications (PTMs). Moreover, new PTMs were observed among thioholgamide derivatives: thiazoline heterocycles, which have not yet been reported for thioamitides, and S-methylmethionine, which is very rare in nature. The obtained library was subsequently used for thioholgamide SAR studies and stability assays.

Directed stepwise tracing of polysynaptic neuronal circuits with replication-deficient pseudorabies virus

Brain functions are accomplished by polysynaptic circuits formed by neurons wired together through multiple orders of synaptic connections. Polysynaptic connectivity has been difficult to examine due to a lack of methods of continuously tracing the pathways in a controlled manner. Here, we demonstrate directed, stepwise retrograde polysynaptic tracing by inducible reconstitution of replication-deficient trans-neuronal pseudorabies virus (PRVΔIE) in the brain. Furthermore, PRVΔIE replication can be temporally restricted to minimize its neurotoxicity. With this tool, we delineate a wiring diagram between the hippocampus and striatum-two major brain systems for learning, memory, and navigation-that consists of projections from specific hippocampal domains to specific striatal areas via distinct intermediate brain regions. Therefore, this inducible PRVΔIE system provides a tool for dissecting polysynaptic circuits underlying complex brain functions.

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.

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

Introduction

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

Methods

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

Results

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

Discussion

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

Expression of blue pigment synthetase a from Streptomyces lavenduale reveals insights on the effects of refactoring biosynthetic megasynthases for heterologous expression in Escherichia coli

High GC bacteria from the genus Streptomyces harbor expansive secondary metabolism. The expression of biosynthetic proteins and the characterization and identification of biological "parts" for synthetic biology purposes from such pathways are of interest. However, the high GC content of proteins from actinomycetes in addition to the large size and multi-domain architecture of many biosynthetic proteins (such as non-ribosomal peptide synthetases; NRPSs, and polyketide synthases; PKSs often called "megasynthases") often presents issues with full-length translation and folding. Here we evaluate a non-ribosomal peptide synthetase (NRPS) from Streptomyces lavenduale, a multidomain "megasynthase" gene that comes from a high GC (72.5%) genome. While a preliminary step in revealing differences, to our knowledge this presents the first head-to-head comparison of codon-optimized sequences versus a native sequence of proteins of streptomycete origin heterologously expressed in E. coli. We found that any disruption in co-translational folding from codon mismatch that reduces the titer of indigoidine is explainable via the formation of more inclusion bodies as opposed to compromising folding or posttranslational modification in the soluble fraction. This result supports that one could apply any refactoring strategies that improve soluble expression in E. coli without concern that the protein that reaches the soluble fraction is differentially folded.

Codon Optimization Improves the Prediction of Xylose Metabolism from Gene Content in Budding Yeasts

Xylose is the second most abundant monomeric sugar in plant biomass. Consequently, xylose catabolism is an ecologically important trait for saprotrophic organisms, as well as a fundamentally important trait for industries that hope to convert plant mass to renewable fuels and other bioproducts using microbial metabolism. Although common across fungi, xylose catabolism is rare within Saccharomycotina, the subphylum that contains most industrially relevant fermentative yeast species. The genomes of several yeasts unable to consume xylose have been previously reported to contain the full set of genes in the XYL pathway, suggesting the absence of a gene-trait correlation for xylose metabolism. Here, we measured growth on xylose and systematically identified XYL pathway orthologs across the genomes of 332 budding yeast species. Although the XYL pathway coevolved with xylose metabolism, we found that pathway presence only predicted xylose catabolism about half of the time, demonstrating that a complete XYL pathway is necessary, but not sufficient, for xylose catabolism. We also found that XYL1 copy number was positively correlated, after phylogenetic correction, with xylose utilization. We then quantified codon usage bias of XYL genes and found that XYL3 codon optimization was significantly higher, after phylogenetic correction, in species able to consume xylose. Finally, we showed that codon optimization of XYL2 was positively correlated, after phylogenetic correction, with growth rates in xylose medium. We conclude that gene content alone is a weak predictor of xylose metabolism and that using codon optimization enhances the prediction of xylose metabolism from yeast genome sequence data.

Comparative Analysis of Codon Usage Patterns in Nuclear and Chloroplast Genome of Dalbergia (Fabaceae)

The Dalbergia plants are widely distributed across more than 130 tropical and subtropical countries and have significant economic and medicinal value. Codon usage bias (CUB) is a critical feature for studying gene function and evolution, which can provide a better understanding of biological gene regulation. In this study, we comprehensively analyzed the CUB patterns of the nuclear genome, chloroplast genome, and gene expression, as well as systematic evolution of Dalbergia species. Our results showed that the synonymous and optimal codons in the coding regions of both nuclear and chloroplast genome of Dalbergia preferred ending with A/U at the third codon base. Natural selection was the primary factor affecting the CUB features. Furthermore, in highly expressed genes of Dalbergia odorifera, we found that genes with stronger CUB exhibited higher expression levels, and these highly expressed genes tended to favor the use of G/C-ending codons. In addition, the branching patterns of the protein-coding sequences and the chloroplast genome sequences were very similar in the systematic tree, and different with the cluster from the CUB of the chloroplast genome. This study highlights the CUB patterns and features of Dalbergia species in different genomes, explores the correlation between CUB preferences and gene expression, and further investigates the systematic evolution of Dalbergia, providing new insights into codon biology and the evolution of Dalbergia plants.

Nascent MSKIK Peptide Cancels Ribosomal Stalling by Arrest Peptides in Escherichia coli

The insertion of the DNA sequence encoding SKIK peptide adjacent to the M start codon of a difficult-to-express protein enhances protein production in Escherichia coli. In this report, we reveal that the increased production of the SKIK-tagged protein is not due to codon usage of the SKIK sequence. Furthermore, we found that insertion of SKIK or MSKIK just before the SecM arrest peptide (FSTPVWISQAQGIRAGP), which causes ribosomal stalling on mRNA, greatly increased the production of the protein containing the SecM arrest peptide in the E. coli reconstituted cell-free protein synthesis system (PURE system). A similar translation enhancement phenomenon by MSKIK was observed for the CmlA leader peptide, a ribosome arrest peptide, whose arrest is induced by chloramphenicol. These results strongly suggest that the nascent MSKIK peptide prevents or releases ribosomal stalling immediately following its generation during the translation process, resulting in an increase of protein production.

Synonymous Codon Variant Analysis for Autophagic Genes Dysregulated in Neurodegeneration

Neurodegenerative disorders are often a culmination of the accumulation of abnormally folded proteins and defective organelles. Autophagy is a process of removing these defective proteins, organelles, and harmful substances from the body, and it works to maintain homeostasis. If autophagic removal of defective proteins has interfered, it affects neuronal health. Some of the autophagic genes are specifically found to be associated with neurodegenerative phenotypes. Non-functional, mutated, or gene copies having silent mutations, often termed synonymous variants, might explain this. However, these synonymous variant which codes for exactly similar proteins have different translation rates, stability, and gene expression profiling. Hence, it would be interesting to study the pattern of synonymous variant usage. In the study, synonymous variant usage in various transcripts of autophagic genes ATG5, ATG7, ATG8A, ATG16, and ATG17/FIP200 reported to cause neurodegeneration (if dysregulated) is studied. These genes were analyzed for their synonymous variant usage; nucleotide composition; any possible nucleotide skew in a gene; physical properties of autophagic protein including GRAVY and AROMA; hydropathicity; instability index; and frequency of acidic, basic, neutral amino acids; and gene expression level. The study will help understand various evolutionary forces acting on these genes and the possible augmentation of a gene if showing unusual behavior.

The efficacy and safety of SARS-CoV-2 vaccines mRNA1273 and BNT162b2 might be complicated by rampant C-to-U RNA editing

The SARS-CoV-2 RNA vaccines are smartly designed to increase the synonymous codon usage by introducing multiple U-to-C mutations. This design would elevate the translation efficiency of vaccine RNAs. However, we found evidence to reason that the designed cytidines might be converted to uridines again by C-to-U RNA deamination in host cells. This C-to-U mechanism might be a main factor that affects the efficacy and safety of RNA vaccines.