A new and updated resource for codon usage tables

Assembly and comparative analysis of the complete mitochondrial and chloroplast genome of Cyperus stoloniferus (Cyperaceae), a coastal plant possessing saline-alkali tolerance

BACKGROUND

Cyperus stoloniferus is an important species in coastal ecosystems and possesses economic and ecological value. To elucidate the structural characteristics, variation, and evolution of the organelle genome of C. stoloniferus, we sequenced, assembled, and compared its mitochondrial and chloroplast genomes.

RESULTS

We assembled the mitochondrial and chloroplast genomes of C. stoloniferus. The total length of the mitochondrial genome (mtDNA) was 927,413 bp, with a GC content of 40.59%. It consists of two circular DNAs, including 37 protein-coding genes (PCGs), 22 tRNAs, and five rRNAs. The length of the chloroplast genome (cpDNA) was 186,204 bp, containing 93 PCGs, 40 tRNAs, and 8 rRNAs. The mtDNA and cpDNA contained 81 and 129 tandem repeats, respectively, and 346 and 1,170 dispersed repeats, respectively, both of which have 270 simple sequence repeats. The third high-frequency codon (RSCU > 1) in the organellar genome tended to end at A or U, whereas the low-frequency codon (RSCU < 1) tended to end at G or C. The RNA editing sites of the PCGs were relatively few, with only 9 and 23 sites in the mtDNA and cpDNA, respectively. A total of 28 mitochondrial plastid DNAs (MTPTs) in the mtDNA were derived from cpDNA, including three complete trnT-GGU, trnH-GUG, and trnS-GCU. Phylogeny and collinearity indicated that the relationship between C. stoloniferus and C. rotundus are closest. The mitochondrial rns gene exhibited the greatest nucleotide variability, whereas the chloroplast gene with the greatest nucleotide variability was infA. Most PCGs in the organellar genome are negatively selected and highly evolutionarily conserved. Only six mitochondrial genes and two chloroplast genes exhibited Ka/Ks > 1; in particular, atp9, atp6, and rps7 may have undergone potential positive selection.

CONCLUSION

We assembled and validated the mtDNA of C. stoloniferus, which contains a 15,034 bp reverse complementary sequence. The organelle genome sequence of C. stoloniferus provides valuable genomic resources for species identification, evolution, and comparative genomic research in Cyperaceae.

Repetitive Sequences, Codon Usage Bias and Phylogenetic Analysis of the Plastome of Miliusa glochidioides

Annonaceae is the largest family in Magnoliales, exhibiting the greatest diversity among and within genera. In this study, we conducted an analysis of repetitive sequences and codon usage bias in the previously acquired plastome of Miliusa glochidioides. Using a concatenated dataset of shared genes, we constructed the phylogenetic relationships among 27 Annonaceae species. The results showed that the size of the plastomes in the Annonaceae ranged from 159 to 202 kb, with the size of the inverted repeat region ranging from 40 to 65 kb. Within the plastome of M. glochidioides, we identified 42 SSRs, 36 tandem repeats, and 9 dispersed repeats. These SSRs consist of three nucleotide types and eight motif types, with a preference for A/T bases, primarily located in the large single-copy regions and intergenic spacers. Tandem and dispersed repeat sequences were predominantly detected in the IR region. Through codon usage bias analysis, we identified 30 high-frequency codons and 11 optimal codons. The plastome of M. glochidioides demonstrated relatively weak codon usage bias, favoring codons with A/T endings, primarily influenced by natural selection. Phylogenetic analysis revealed that all four subfamilies formed monophyletic groups, with Cananga odorata (Ambavioideae) and Anaxagorea javanica (Anaxagoreoideae) successively nested outside Annonoideae + Malmeoideae. These findings improve our understanding of the plastome of M. glochidioides and provide additional insights for studying plastome evolution in Annonaceae.

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

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

Chloroplast Genomes of Vitis flexuosa and Vitis amurensis: Molecular Structure, Phylogenetic, and Comparative Analyses for Wild Plant Conservation

The chloroplast genome plays a crucial role in elucidating genetic diversity and phylogenetic relationships. Vitis vinifera L. (grapevine) is an economically important species, prompting exploration of wild genetic resources to enhance stress resilience. We meticulously assembled the chloroplast genomes of two Korean Vitis L. species, V. flexuosa Thunb. and V. amurensis Rupr., contributing valuable data to the Korea Crop Wild Relatives inventory. Through exhaustive specimen collection spanning diverse ecological niches across South Korea, we ensured comprehensive representation of genetic diversity. Our analysis, which included rigorous codon usage bias assessment and repeat analysis, provides valuable insights into amino acid preferences and facilitates the identification of potential molecular markers. The assembled chloroplast genomes were subjected to meticulous annotation, revealing divergence hotspots enriched with nucleotide diversity, thereby presenting promising candidates for DNA barcodes. Additionally, phylogenetic analysis reaffirmed intra-genus relationships and identified related crops, shedding light on evolutionary patterns within the genus. Comparative examination with chloroplast genomes of other crops uncovered conserved sequences and variable regions, offering critical insights into genetic evolution and adaptation. Our study advances the understanding of chloroplast genomes, genetic diversity, and phylogenetic relationships within Vitis species, thereby laying a foundation for enhancing grapevine genetic diversity and resilience to environmental challenges.

Recombinant multiepitope proteins expressed in Escherichia coli cells and their potential for immunodiagnosis

Recombinant multiepitope proteins (RMPs) are a promising alternative for application in diagnostic tests and, given their wide application in the most diverse diseases, this review article aims to survey the use of these antigens for diagnosis, as well as discuss the main points surrounding these antigens. RMPs usually consisting of linear, immunodominant, and phylogenetically conserved epitopes, has been applied in the experimental diagnosis of various human and animal diseases, such as leishmaniasis, brucellosis, cysticercosis, Chagas disease, hepatitis, leptospirosis, leprosy, filariasis, schistosomiasis, dengue, and COVID-19. The synthetic genes for these epitopes are joined to code a single RMP, either with spacers or fused, with different biochemical properties. The epitopes' high density within the RMPs contributes to a high degree of sensitivity and specificity. The RMPs can also sidestep the need for multiple peptide synthesis or multiple recombinant proteins, reducing costs and enhancing the standardization conditions for immunoassays. Methods such as bioinformatics and circular dichroism have been widely applied in the development of new RMPs, helping to guide their construction and better understand their structure. Several RMPs have been expressed, mainly using the Escherichia coli expression system, highlighting the importance of these cells in the biotechnological field. In fact, technological advances in this area, offering a wide range of different strains to be used, make these cells the most widely used expression platform. RMPs have been experimentally used to diagnose a broad range of illnesses in the laboratory, suggesting they could also be useful for accurate diagnoses commercially. On this point, the RMP method offers a tempting substitute for the production of promising antigens used to assemble commercial diagnostic kits.

Population-specific variations in KCNH2 predispose patients to delayed ventricular repolarization upon dihydroartemisinin-piperaquine therapy

Dihydroartemisinin-piperaquine is efficacious for the treatment of uncomplicated malaria and its use is increasing globally. Despite the positive results in fighting malaria, inhibition of the Kv11.1 channel (hERG; encoded by the KCNH2 gene) by piperaquine has raised concerns about cardiac safety. Whether genetic factors could modulate the risk of piperaquine-mediated QT prolongations remained unclear. Here, we first profiled the genetic landscape of KCNH2 variability using data from 141,614 individuals. Overall, we found 1,007 exonic variants distributed over the entire gene body, 555 of which were missense. By optimizing the gene-specific parametrization of 16 partly orthogonal computational algorithms, we developed a KCNH2-specific ensemble classifier that identified a total of 116 putatively deleterious missense variations. To evaluate the clinical relevance of KCNH2 variability, we then sequenced 293 Malian patients with uncomplicated malaria and identified 13 variations within the voltage sensing and pore domains of Kv11.1 that directly interact with channel blockers. Cross-referencing of genetic and electrocardiographic data before and after piperaquine exposure revealed that carriers of two common variants, rs1805121 and rs41314375, experienced significantly higher QT prolongations (ΔQTc of 41.8 ms and 61 ms, respectively, vs 14.4 ms in controls) with more than 50% of carriers having increases in QTc >30 ms. Furthermore, we identified three carriers of rare population-specific variations who experienced clinically relevant delayed ventricular repolarization. Combined, our results map population-scale genetic variability of KCNH2 and identify genetic biomarkers for piperaquine-induced QT prolongation that could help to flag at-risk patients and optimize efficacy and adherence to antimalarial therapy.

The complete chloroplast genome and phylogentic results support the species position of Swertia banzragczii and Swertia marginata (Gentianaceae) in Mongolia

BACKGROUND

Swertia banzragczii and S. marginata are important medicinal species in Mongolia. However, their taxonomic positions and genetic backgrounds remain unknown. In this study, we explored the complete chloroplast genomes and DNA barcoding of these species and compared them with those of closely related species within the subgenus to determine their taxonomic positions and phylogenetic relationships.

RESULT

The chloroplast genomes of S. banzragczii and S. marginata encoded 114 genes, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Among them, 16 genes contained a single intron, and 2 genes had two introns. Closely related species had a conserved genome structure and gene content. Only differences in genome length were noticed, which were caused by the expansion and contraction of the inverted repeat (IR) region and loss of exons in some genes. The trnH-GUG-psbA and trnD-GUC-trnY-GUA intergenic regions had high genetic diversity within Swertia plastomes. Overall, S. banzragczii and S. marginata are true species and belong to the subgenus Swertia.

CONCLUSIONS

These results provide valuable genetic and morphological information on rare and subendemic Swertia species in Mongolia, which can be used for further advanced studies on the Swertia genus.

Codon-optimization in gene therapy: promises, prospects and challenges

Codon optimization has evolved to enhance protein expression efficiency by exploiting the genetic code's redundancy, allowing for multiple codon options for a single amino acid. Initially observed in E. coli, optimal codon usage correlates with high gene expression, which has propelled applications expanding from basic research to biopharmaceuticals and vaccine development. The method is especially valuable for adjusting immune responses in gene therapies and has the potenial to create tissue-specific therapies. However, challenges persist, such as the risk of unintended effects on protein function and the complexity of evaluating optimization effectiveness. Despite these issues, codon optimization is crucial in advancing gene therapeutics. This study provides a comprehensive review of the current metrics for codon-optimization, and its practical usage in research and clinical applications, in the context of gene therapy.

DNA-Based Nonviral Gene Therapy─Challenging but Promising

Over the past decades, significant progress has been made in utilizing nucleic acids, including DNA and RNA molecules, for therapeutic purposes. For DNA molecules, although various DNA delivery systems have been established, viral vector systems are the go-to choice for large-scale commercial applications. However, viral systems have certain disadvantages such as immune response, limited payload capacity, insertional mutagenesis and pre-existing immunity. In contrast, nonviral systems are less immunogenic, not size limited, safer, and easier for manufacturing compared with viral systems. What's more, nonviral DNA vectors have demonstrated their capacity to mediate specific protein expression in vivo for diverse therapeutic objectives containing a wide range of diseases such as cancer, rare diseases, neurodegenerative diseases, and infectious diseases, yielding promising therapeutic outcomes. However, exogenous plasmid DNA is prone to degrade and has poor immunogenicity in vivo. Thus, various strategies have been developed: (i) designing novel plasmids with special structures, (ii) optimizing plasmid sequences for higher expression, and (iii) developing more efficient nonviral DNA delivery systems. Based on these strategies, many interesting clinical results have been reported. This Review discusses the development of DNA-based nonviral gene therapy, including novel plasmids, nonviral delivery systems, clinical advances, and prospects. These developments hold great potential for enhancing the efficacy and safety of nonviral gene therapy and expanding its applications in the treatment of various diseases.

Antiviral responses versus virus-induced cellular shutoff: a game of thrones between influenza A virus NS1 and SARS-CoV-2 Nsp1

Following virus recognition of host cell receptors and viral particle/genome internalization, viruses replicate in the host via hijacking essential host cell machinery components to evade the provoked antiviral innate immunity against the invading pathogen. Respiratory viral infections are usually acute with the ability to activate pattern recognition receptors (PRRs) in/on host cells, resulting in the production and release of interferons (IFNs), proinflammatory cytokines, chemokines, and IFN-stimulated genes (ISGs) to reduce virus fitness and mitigate infection. Nevertheless, the game between viruses and the host is a complicated and dynamic process, in which they restrict each other via specific factors to maintain their own advantages and win this game. The primary role of the non-structural protein 1 (NS1 and Nsp1) of influenza A viruses (IAV) and the pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), respectively, is to control antiviral host-induced innate immune responses. This review provides a comprehensive overview of the genesis, spatial structure, viral and cellular interactors, and the mechanisms underlying the unique biological functions of IAV NS1 and SARS-CoV-2 Nsp1 in infected host cells. We also highlight the role of both non-structural proteins in modulating viral replication and pathogenicity. Eventually, and because of their important role during viral infection, we also describe their promising potential as targets for antiviral therapy and the development of live attenuated vaccines (LAV). Conclusively, both IAV NS1 and SARS-CoV-2 Nsp1 play an important role in virus-host interactions, viral replication, and pathogenesis, and pave the way to develop novel prophylactic and/or therapeutic interventions for the treatment of these important human respiratory viral pathogens.

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

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

Bioengineering of the Marine Diatom Phaeodactylum tricornutum with Cannabis Genes Enables the Production of the Cannabinoid Precursor, Olivetolic Acid

The increasing demand for novel natural compounds has prompted the exploration of innovative approaches in bioengineering. This study investigates the bioengineering potential of the marine diatom Phaeodactylum tricornutum through the introduction of cannabis genes, specifically, tetraketide synthase (TKS), and olivetolic acid cyclase (OAC), for the production of the cannabinoid precursor, olivetolic acid (OA). P. tricornutum is a promising biotechnological platform due to its fast growth rate, amenability to genetic manipulation, and ability to produce valuable compounds. Through genetic engineering techniques, we successfully integrated the cannabis genes TKS and OAC into the diatom. P. tricornutum transconjugants expressing these genes showed the production of the recombinant TKS and OAC enzymes, detected via Western blot analysis, and the production of cannabinoids precursor (OA) detected using the HPLC/UV spectrum when compared to the wild-type strain. Quantitative analysis revealed significant olivetolic acid accumulation (0.6-2.6 mg/L), demonstrating the successful integration and functionality of the heterologous genes. Furthermore, the introduction of TKS and OAC genes led to the synthesis of novel molecules, potentially expanding the repertoire of bioactive compounds accessible through diatom-based biotechnology. This study demonstrates the successful bioengineering of P. tricornutum with cannabis genes, enabling the production of OA as a precursor for cannabinoid production and the synthesis of novel molecules with potential pharmaceutical applications.

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.

An exploratory in silico comparison of open-source codon harmonization tools

BACKGROUND

Not changing the native constitution of genes prior to their expression by a heterologous host can affect the amount of proteins synthesized as well as their folding, hampering their activity and even cell viability. Over the past decades, several strategies have been developed to optimize the translation of heterologous genes by accommodating the difference in codon usage between species. While there have been a handful of studies assessing various codon optimization strategies, to the best of our knowledge, no research has been performed towards the evaluation and comparison of codon harmonization algorithms. To highlight their importance and encourage meaningful discussion, we compared different open-source codon harmonization tools pertaining to their in silico performance, and we investigated the influence of different gene-specific factors.

RESULTS

In total, 27 genes were harmonized with four tools toward two different heterologous hosts. The difference in %MinMax values between the harmonized and the original sequences was calculated (ΔMinMax), and statistical analysis of the obtained results was carried out. It became clear that not all tools perform similarly, and the choice of tool should depend on the intended application. Almost all biological factors under investigation (GC content, RNA secondary structures and choice of heterologous host) had a significant influence on the harmonization results and thus must be taken into account. These findings were substantiated using a validation dataset consisting of 8 strategically chosen genes.

CONCLUSIONS

Due to the size of the dataset, no complex models could be developed. However, this initial study showcases significant differences between the results of various codon harmonization tools. Although more elaborate investigation is needed, it is clear that biological factors such as GC content, RNA secondary structures and heterologous hosts must be taken into account when selecting the codon harmonization tool.

AAV-based in vivo gene therapy for neurological disorders

Recent advancements in gene supplementation therapy are expanding the options for the treatment of neurological disorders. Among the available delivery vehicles, adeno-associated virus (AAV) is often the favoured vector. However, the results have been variable, with some trials dramatically altering the course of disease whereas others have shown negligible efficacy or even unforeseen toxicity. Unlike traditional drug development with small molecules, therapeutic profiles of AAV gene therapies are dependent on both the AAV capsid and the therapeutic transgene. In this rapidly evolving field, numerous clinical trials of gene supplementation for neurological disorders are ongoing. Knowledge is growing about factors that impact the translation of preclinical studies to humans, including the administration route, timing of treatment, immune responses and limitations of available model systems. The field is also developing potential solutions to mitigate adverse effects, including AAV capsid engineering and designs to regulate transgene expression. At the same time, preclinical research is addressing new frontiers of gene supplementation for neurological disorders, with a focus on mitochondrial and neurodevelopmental disorders. In this Review, we describe the current state of AAV-mediated neurological gene supplementation therapy, including critical factors for optimizing the safety and efficacy of treatments, as well as unmet needs in this field.

Time-calibrated phylogeny and full mitogenome sequence of the Galapagos sea lion (Zalophus wollebaeki) from scat DNA

Background

The Galapagos sea lion, Zalophus wollebaeki, is an endemic and endangered otariid, which is considered as a sentinel species of ecosystem dynamics in the Galapagos archipelago. Mitochondrial DNA is an important tool in phylogenetic and population genetic inference. In this work we use Illumina sequencing to complement the mitogenomic resources for Zalophus genus-the other two species employed Sanger sequencing-by a complete mitochondrial genome and a molecular clock of this species, which is not present in any case.

Materials and Methods

We used DNA obtained from a fresh scat sample of a Galapagos sea lion and shotgun-sequenced it on the Illumina NextSeq platform. The obtained raw reads were processed using the GetOrganelle software to filter the mitochondrial Zalophus DNA reads (∼16% survive the filtration), assemble them, and set up a molecular clock.

Results

From the obtained 3,511,116 raw reads, we were able to assemble a full mitogenome of a length of 16,676 bp, consisting of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNA), and two ribosomal RNAs (rRNA). A time-calibrated phylogeny confirmed the phylogenetic position of Z. wollebaeki in a clade with Z. californianus, and Z. japonicus, and sister to Z. californianus; as well as establishing the divergence time for Z. wollebaeki 0.65 million years ago. Our study illustrates the possibility of seamlessly sequencing full mitochondrial genomes from fresh scat samples of marine mammals.

Human Schlafen 11 inhibits influenza A virus production

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

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.

Schlafen 12 restricts HIV-1 latency reversal by a codon-usage dependent post-transcriptional block in CD4+ T cells

Latency is a major barrier towards virus elimination in HIV-1-infected individuals. Yet, the mechanisms that contribute to the maintenance of HIV-1 latency are incompletely understood. Here we describe the Schlafen 12 protein (SLFN12) as an HIV-1 restriction factor that establishes a post-transcriptional block in HIV-1-infected cells and thereby inhibits HIV-1 replication and virus reactivation from latently infected cells. The inhibitory activity is dependent on the HIV-1 codon usage and on the SLFN12 RNase active sites. Within HIV-1-infected individuals, SLFN12 expression in PBMCs correlated with HIV-1 plasma viral loads and proviral loads suggesting a link with the general activation of the immune system. Using an RNA FISH-Flow HIV-1 reactivation assay, we demonstrate that SLFN12 expression is enriched in infected cells positive for HIV-1 transcripts but negative for HIV-1 proteins. Thus, codon-usage dependent translation inhibition of HIV-1 proteins participates in HIV-1 latency and can restrict the amount of virus release after latency reversal.

Ribosomes in RNA Granules Are Stalled on mRNA Sequences That Are Consensus Sites for FMRP Association

Local translation in neurons is partly mediated by the reactivation of stalled polysomes. Stalled polysomes may be enriched within the granule fraction, defined as the pellet of sucrose gradients used to separate polysomes from monosomes. The mechanism of how elongating ribosomes are reversibly stalled and unstalled on mRNAs is still unclear. In the present study, we characterize the ribosomes in the granule fraction using immunoblotting, cryogenic electron microscopy (cryo-EM), and ribosome profiling. We find that this fraction, isolated from 5-d-old rat brains of both sexes, is enriched in proteins implicated in stalled polysome function, such as the fragile X mental retardation protein (FMRP) and Up-frameshift mutation 1 homologue. Cryo-EM analysis of ribosomes in this fraction indicates they are stalled, mainly in the hybrid state. Ribosome profiling of this fraction reveals (1) an enrichment for footprint reads of mRNAs that interact with FMRPs and are associated with stalled polysomes, (2) an abundance of footprint reads derived from mRNAs of cytoskeletal proteins implicated in neuronal development, and (3) increased ribosome occupancy on mRNAs encoding RNA binding proteins. Compared with those usually found in ribosome profiling studies, the footprint reads were longer and were mapped to reproducible peaks in the mRNAs. These peaks were enriched in motifs previously associated with mRNAs cross-linked to FMRP in vivo, independently linking the ribosomes in the granule fraction to the ribosomes associated with FMRP in the cell. The data supports a model in which specific sequences in mRNAs act to stall ribosomes during translation elongation in neurons.SIGNIFICANCE STATEMENT Neurons send mRNAs to synapses in RNA granules, where they are not translated until an appropriate stimulus is given. Here, we characterize a granule fraction obtained from sucrose gradients and show that polysomes in this fraction are stalled on consensus sequences in a specific state of translational arrest with extended ribosome-protected fragments. This finding greatly increases our understanding of how neurons use specialized mechanisms to regulate translation and suggests that many studies on neuronal translation may need to be re-evaluated to include the large fraction of neuronal polysomes found in the pellet of sucrose gradients used to isolate polysomes.

Detection of unintended on-target effects in CRISPR genome editing by DNA donors carrying diagnostic substitutions

CRISPR nucleases can introduce double-stranded DNA breaks in genomes at positions specified by guide RNAs. When repaired by the cell, this may result in the introduction of insertions and deletions or nucleotide substitutions provided by exogenous DNA donors. However, cellular repair can also result in unintended on-target effects, primarily larger deletions and loss of heterozygosity due to gene conversion. Here we present a strategy that allows easy and reliable detection of unintended on-target effects as well as the generation of control cells that carry wild-type alleles but have demonstratively undergone genome editing at the target site. Our 'sequence-ascertained favorable editing' (SAFE) donor approach relies on the use of DNA donor mixtures containing the desired nucleotide substitutions or the wild-type alleles together with combinations of additional 'diagnostic' substitutions unlikely to have any effects. Sequencing of the target sites then results in that two different sequences are seen when both chromosomes are edited with 'SAFE' donors containing different sets of substitutions, while a single sequence indicates unintended effects such as deletions or gene conversion. We analyzed more than 850 human embryonic stem cell clones edited with 'SAFE' donors and detect all copy number changes and almost all clones with gene conversion.

Analysis of 3.5 million SARS-CoV-2 sequences reveals unique mutational trends with consistent nucleotide and codon frequencies

BACKGROUND

Since the onset of the SARS-CoV-2 pandemic, bioinformatic analyses have been performed to understand the nucleotide and synonymous codon usage features and mutational patterns of the virus. However, comparatively few have attempted to perform such analyses on a considerably large cohort of viral genomes while organizing the plethora of available sequence data for a month-by-month analysis to observe changes over time. Here, we aimed to perform sequence composition and mutation analysis of SARS-CoV-2, separating sequences by gene, clade, and timepoints, and contrast the mutational profile of SARS-CoV-2 to other comparable RNA viruses.

METHODS

Using a cleaned, filtered, and pre-aligned dataset of over 3.5 million sequences downloaded from the GISAID database, we computed nucleotide and codon usage statistics, including calculation of relative synonymous codon usage values. We then calculated codon adaptation index (CAI) changes and a nonsynonymous/synonymous mutation ratio (dN/dS) over time for our dataset. Finally, we compiled information on the types of mutations occurring for SARS-CoV-2 and other comparable RNA viruses, and generated heatmaps showing codon and nucleotide composition at high entropy positions along the Spike sequence.

RESULTS

We show that nucleotide and codon usage metrics remain relatively consistent over the 32-month span, though there are significant differences between clades within each gene at various timepoints. CAI and dN/dS values vary substantially between different timepoints and different genes, with Spike gene on average showing both the highest CAI and dN/dS values. Mutational analysis showed that SARS-CoV-2 Spike has a higher proportion of nonsynonymous mutations than analogous genes in other RNA viruses, with nonsynonymous mutations outnumbering synonymous ones by up to 20:1. However, at several specific positions, synonymous mutations were overwhelmingly predominant.

CONCLUSIONS

Our multifaceted analysis covering both the composition and mutation signature of SARS-CoV-2 gives valuable insight into the nucleotide frequency and codon usage heterogeneity of SARS-CoV-2 over time, and its unique mutational profile compared to other RNA viruses.

In silico investigation of uncoupling protein function in avian genomes

Introduction

The uncoupling proteins (UCPs) are involved in lipid metabolism and belong to a family of mitochondrial anionic transporters. In poultry, only one UCP homologue has been identified and experimentally shown to be associated with growth, feed conversion ratio, and abdominal fat according to its predominant expression in bird muscles. In endotherm birds, cell metabolic efficiency can be tuned by the rate of mitochondrial coupling. Thus, avUCP may be a key contributor to controlling metabolic rate during particular environmental changes.

Methods

This study aimed to perform a set of in-silico investigations primarily focused on the structural, biological, and biomimetic functions of avUCP. Thereby, using in silico genome analyses among 8 avian species (chicken, turkey, swallow, manakin, sparrow, wagtail, pigeon, and mallard) and a series of bioinformatic approaches, we provide phylogenetic inference and comparative genomics of avUCPs and investigate whether sequence variation can alter coding sequence characteristics, the protein structure, and its biological features. Complementarily, a combination of literature mining and prediction approaches was also applied to predict the gene networks of avUCP to identify genes, pathways, and biological crosstalk associated with avUCP function.

Results

The results showed the evolutionary alteration of UCP proteins in different avian species. Uncoupling proteins in avian species are highly conserved trans membrane proteins as seen by sequence alignment, physio-chemical parameters, and predicted protein structures. Taken together, avUCP has the potential to be considered a functional marker for the identification of cell metabolic state, thermogenesis, and oxidative stress caused by cold, heat, fasting, transfer, and other chemical stimuli stresses in birds. It can also be deduced that avUCP, in migrant or domestic birds, may increase heat stress resistance by reducing fatty acid transport/b-oxidation and thermoregulation alongside antioxidant defense mechanisms. The predicted gene network for avUCP highlighted a cluster of 21 genes involved in response to stress and 28 genes related to lipid metabolism and the proton buffering system. Finally, among 11 enriched pathways, crosstalk of 5 signaling pathways including MAPK, adipocytokine, mTOR, insulin, ErbB, and GnRH was predicted, indicating a possible combination of positive or negative feedback among pathways to regulate avUCP functions.

Discussion

Genetic selection for fast-growing commercial poultry has unintentionally increased susceptibility to many kinds of oxidative stress, and so avUCP could be considered as a potential candidate gene for balancing energy expenditure and reactive oxygen species production, especially in breeding programs. In conclusion, avUCP can be introduced as a pleiotropic gene that requires the contribution of regulatory genes, hormones, pathways, and genetic crosstalk to allow its finely-tuned function.

CoDe: a web-based tool for codon deoptimization

Summary

We have developed a web-based tool, CoDe (Codon Deoptimization) that deoptimizes genetic sequences based on different codon usage bias, ultimately reducing expression of the corresponding protein. The tool could also deoptimize the sequence for a specific region and/or selected amino acid(s). Moreover, CoDe can highlight sites targeted by restriction enzymes in the wild-type and codon-deoptimized sequences. Importantly, our web-based tool has a user-friendly interface with flexible options to download results.

Availability and implementation

The web-based tool CoDe is freely available at https://web.iitm.ac.in/bioinfo2/codeop/landing_page.html.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Comparison of Boraginales Plastomes: Insights into Codon Usage Bias, Adaptive Evolution, and Phylogenetic Relationships

The Boraginales (Boraginaceae a.l.) comprise more than 2450 species worldwide. However, little knowledge exists of the characteristics of the complete plastid genome. In this study, three new sequences representing the first pt genome of Heliotropiaceae and Cordiaceae were assembled and compared with other Boraginales species. The pt genome sizes of Cordia dichotoma, Heliotropium arborescens, and Tournefortia montana were 151,990 bp, 156,243 bp, and 155,891 bp, respectively. Multiple optimal codons were identified, which may provide meaningful information for enhancing the gene expression of Boraginales species. Furthermore, codon usage bias analyses revealed that natural selection and other factors may dominate codon usage patterns in the Boraginales species. The boundaries of the IR/LSC and IR/SSC regions were significantly different, and we also found a signal of obvious IR region expansion in the pt genome of Nonea vesicaria and Arnebia euchroma. Genes with high nucleic acid diversity (pi) values were also calculated, which may be used as potential DNA barcodes to investigate the phylogenetic relationships in Boraginales. psaI, rpl33, rpl36, and rps19 were found to be under positive selection, and these genes play an important role in our understanding of the adaptive evolution of the Boraginales species. Phylogenetic analyses implied that Boraginales can be divided into two groups. The existence of two subfamilies (Lithospermeae and Boragineae) in Boraginaceae is also strongly supported. Our study provides valuable information on pt genome evolution and phylogenetic relationships in the Boraginales species.

The complete chloroplast genome sequences of three Broussonetia species and comparative analysis within the Moraceae

Background

Species of Broussonetia (family Moraceae) are commonly used to make textiles and high-grade paper. The distribution of Broussonetia papyrifera L. is considered to be related to the spread and location of humans. The complete chloroplast (cp) genomes of B. papyrifera, Broussonetia kazinoki Sieb., and Broussonetia kaempferi Sieb. were analyzed to better understand the status and evolutionary biology of the genus Broussonetia.

Methods

The cp genomes were assembled and characterized using SOAPdenovo2 and DOGMA. Phylogenetic and molecular dating analysis were performed using the concatenated nucleotide sequences of 35 species in the Moraceae family and were based on 66 protein-coding genes (PCGs). An analysis of the sequence divergence (pi) of each PCG among the 35 cp genomes was conducted using DnaSP v6. Codon usage indices were calculated using the CodonW program.

Results

All three cp genomes had the typical land plant quadripartite structure, ranging in size from 160,239 bp to 160,841 bp. The ribosomal protein L22 gene (RPL22) was either incomplete or missing in all three Broussonetia species. Phylogenetic analysis revealed two clades. Clade 1 included Morus and Artocarpus, whereas clade 2 included the other seven genera. Malaisia scandens Lour. was clustered within the genus Broussonetia. The differentiation of Broussonetia was estimated to have taken place 26 million years ago. The PCGs' pi values ranged from 0.0005 to 0.0419, indicating small differences within the Moraceae family. The distribution of most of the genes in the effective number of codons plot (ENc-plot) fell on or near the trend line; the slopes of the trend line of neutrality plots were within the range of 0.0363-0.171. These results will facilitate the identification, taxonomy, and utilization of the Broussonetia species and further the evolutionary studies of the Moraceae family.

Composition and Codon Usage Pattern Results in Divergence of the Zinc Binuclear Cluster (Zn(II)2Cys6) Sequences among Ascomycetes Plant Pathogenic Fungi

Zinc binuclear cluster proteins (ZBC; Zn(II)₂Cys₆) are unique to the fungi kingdom and associated with a series of functions, viz., the utilization of macromolecules, stress tolerance, and most importantly, host-pathogen interactions by imparting virulence to the pathogen. Codon usage bias (CUB) is the phenomenon of using synonymous codons in a non-uniform fashion during the translation event, which has arisen because of interactions among evolutionary forces. The Zn(II)₂Cys₆ coding sequences from nine Ascomycetes plant pathogenic species and model system yeast were analysed for compositional and codon usage bias patterns. The clustering analysis diverged the Ascomycetes fungi into two clusters. The nucleotide compositional and relative synonymous codon usage (RSCU) analysis indicated GC biasness toward Ascomycetes fungi compared with the model system S. cerevisiae, which tends to be AT-rich. Further, plant pathogenic Ascomycetes fungi belonging to cluster-2 showed a higher number of GC-rich high-frequency codons than cluster-1 and was exclusively AT-rich in S. cerevisiae. The current investigation also showed the mutual effect of the two evolutionary forces, viz. natural selection and compositional constraints, on the CUB of Zn(II)₂Cys₆ genes. The perseverance of GC-rich codons of Zn(II)₂Cys₆ in Ascomycetes could facilitate the invasion process. The findings of the current investigation show the role of CUB and nucleotide composition in the evolutionary divergence of Ascomycetes plant pathogens and paves the way to target specific codons and sequences to modulate host-pathogen interactions through genome editing and functional genomics tools.

Potassium viroporins as model systems for understanding eukaryotic ion channel behaviour

Ion channels are membrane proteins essential for a plethora of cellular functions including maintaining cell shape, ion homeostasis, cardiac rhythm and action potential in neurons. The complexity and often extensive structure of eukaryotic membrane proteins makes it difficult to understand their basic biological regulation. Therefore, this article suggests, viroporins - the miniature versions of eukaryotic protein homologs from viruses - might serve as model systems to provide insights into behaviour of eukaryotic ion channels in general. The structural requirements for correct assembly of the channel along with the basic functional properties of a K+ channel exist in the minimal design of the viral K+ channels from two viruses, Chlorella virus (Kcv) and Ectocarpus siliculosus virus (Kesv). These small viral proteins readily assemble into tetramers and they sort in cells to distinct target membranes. When these viruses-encoded channels are expressed into the mammalian cells, they utilise their protein machinery and hence can serve as excellent tools to study the cells protein sorting machinery. This combination of small size and robust function makes viral K+ channels a valuable model system for detection of basic structure-function correlations. It is believed that molecular and physiochemical analyses of these viroporins may serve as basis for the development of inhibitors or modulators to ion channel activity for targeting ion channel diseases - so called channelopathies. Therefore, it may provide a potential different scope for molecular pharmacology studies aiming at novel and innovative therapeutics associated with channel related diseases. This article reviews the structural and functional properties of Kcv and Kesv upon expression in mammalian cells and Xenopus oocytes. The mechanisms behind differential protein sorting in Kcv and Kesv are also thoroughly discussed.

CHIKV infection reprograms codon optimality to favor viral RNA translation by altering the tRNA epitranscriptome

Ample evidence indicates that codon usage bias regulates gene expression. How viruses, such as the emerging mosquito-borne Chikungunya virus (CHIKV), express their genomes at high levels despite an enrichment in rare codons remains a puzzling question. Using ribosome footprinting, we analyze translational changes that occur upon CHIKV infection. We show that CHIKV infection induces codon-specific reprogramming of the host translation machinery to favor the translation of viral RNA genomes over host mRNAs with an otherwise optimal codon usage. This reprogramming was mostly apparent at the endoplasmic reticulum, where CHIKV RNAs show high ribosome occupancy. Mechanistically, it involves CHIKV-induced overexpression of KIAA1456, an enzyme that modifies the wobble U34 position in the anticodon of tRNAs, which is required for proper decoding of codons that are highly enriched in CHIKV RNAs. Our findings demonstrate an unprecedented interplay of viruses with the host tRNA epitranscriptome to adapt the host translation machinery to viral production.

Viruses completely depend on the host translational machinery, but their genomes are often enriched in rare codons and therefore should be translated with poor efficiency. Here, Jungfleisch et al. apply Ribo-Seq and RNASeq to provide a global view on the translational changes occurring during Chikungunya virus (CHIKV) infection. CHIKV infection induces a codon-specific reprogramming of the host translation machinery to favor the translation of viral RNA genomes over host mRNAs via tRNA modification.

The Codon Statistics Database: a Database of Codon Usage Bias

We present the Codon Statistics Database, an online database that contains codon usage statistics for all the species with reference or representative genomes in RefSeq (over 15,000). The user can search for any species and access two sets of tables. One set lists, for each codon, the frequency, the Relative Synonymous Codon Usage, and whether the codon is preferred. Another set of tables lists, for each gene, its GC content, Effective Number of Codons, Codon Adaptation Index, and frequency of optimal codons. Equivalent tables can be accessed for (1) all nuclear genes, (2) nuclear genes encoding ribosomal proteins, (3) mitochondrial genes, and (4) chloroplast genes (if available in the relevant assembly). The user can also search for any taxonomic group (e.g., “primates”) and obtain a table comparing all the species in the group. The database is free to access without registration at http://codonstatsdb.unr.edu.

Enhanced effective codon numbers to understand codon usage bias

Codon usage bias is a well recognized phenomenon but the relative influence of its major causes: G+C content, mutational biases, and selection, are often difficult to disentangle. This paper presents methods to calculate modified effective codon numbers that allow the investigation of the sources of codon bias and how genes or organisms have their codon biases shaped. In particular, it demonstrates that variation in codon usage bias across organisms is likely driven more by likely mutational forces while the variation in codon usage bias within genomes is likely driven by codon selectional forces.

Structural, functional, and immunogenicity implications of F9 gene recoding

Hemophilia B is a blood clotting disorder caused by deficient activity of coagulation factor IX (FIX). Multiple recombinant FIX proteins are currently approved to treat hemophilia B, and several gene therapy products are currently being developed. Codon optimization is a frequently used technique in the pharmaceutical industry to improve recombinant protein expression by recoding a coding sequence using multiple synonymous codon substitutions. The underlying assumption of this gene recoding is that synonymous substitutions do not alter protein characteristics because the primary sequence of the protein remains unchanged. However, a critical body of evidence shows that synonymous variants can affect cotranslational folding and protein function. Gene recoding could potentially alter the structure, function, and in vivo immunogenicity of recoded therapeutic proteins. Here, we evaluated multiple recoded variants of F9 designed to further explore the effects of codon usage bias on protein properties. The detailed evaluation of these constructs showed altered conformations, and assessment of translation kinetics by ribosome profiling revealed differences in local translation kinetics. Assessment of wild-type and recoded constructs using a major histocompatibility complex (MHC)-associated peptide proteomics assay showed distinct presentation of FIX-derived peptides bound to MHC class II molecules, suggesting that despite identical amino acid sequence, recoded proteins could exhibit different immunogenicity risks. Posttranslational modification analysis indicated that overexpression from gene recoding results in suboptimal posttranslational processing. Overall, our results highlight potential functional and immunogenicity concerns associated with gene-recoded F9 products. These findings have general applicability and implications for other gene-recoded recombinant proteins.

Activation and execution of the hepatic integrated stress response by dietary essential amino acid deprivation is amino acid specific

Dietary removal of an essential amino acid (EAA) triggers the integrated stress response (ISR) in liver. Herein, we explored the mechanisms that activate the ISR and execute changes in transcription and translation according to the missing EAA. Wild‐type mice and mice lacking general control nonderepressible 2 (Gcn2) were fed an amino acid complete diet or a diet devoid of either leucine or sulfur amino acids (methionine and cysteine). Serum and liver leucine concentrations were significantly reduced within the first 6 h of feeding a diet lacking leucine, corresponding with modest, GCN2‐dependent increases in Atf4 mRNA translation and induction of selected ISR target genes (Fgf21, Slc7a5, Slc7a11). In contrast, dietary removal of the sulfur amino acids lowered serum methionine, but not intracellular methionine, and yet hepatic mRNA abundance of Atf4, Fgf21, Slc7a5, Slc7a11 substantially increased regardless of GCN2 status. Liver tRNA charging levels did not correlate with intracellular EAA concentrations or GCN2 status and remained similar to mice fed a complete diet. Furthermore, loss of Gcn2 increased the occurrence of ribosome collisions in liver and derepressed mechanistic target of rapamycin complex 1 signal transduction, but these changes did not influence execution of the ISR. We conclude that ISR activation is directed by intracellular EAA concentrations, but ISR execution is not. Furthermore, a diet devoid of sulfur amino acids does not require GCN2 for the ISR to execute changes to the transcriptome.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Evaluation of Methylotrophic Yeast Ogataea thermomethanolica TBRC 656 as a Heterologous Host for Production of an Animal Vaccine Candidate

Multiple yeast strains have been developed into versatile heterologous protein expression platforms. Earlier works showed that Ogataea thermomethanolica TBRC 656 (OT), a thermotolerant methylotrophic yeast, can efficiently produce several industrial enzymes. In this work, we demonstrated the potential of this platform for biopharmaceutical manufacturing. Using a swine vaccine candidate as a model, we showed that OT can be optimized to express and secrete the antigen based on porcine circovirus type 2d capsid protein at a respectable yield. Crucial steps for yield improvement include codon optimization and reduction of OT protease activities. The antigen produced in this system could be purified efficiently and induce robust antibody response in test animals. Improvements in this platform, especially more efficient secretion and reduced extracellular proteases, would extend its potential as a competitive platform for biopharmaceutical industries.

Large-Scale Analysis of SARS-CoV-2 Synonymous Mutations Reveals the Adaptation to the Human Codon Usage During the Virus Evolution

Many large national and transnational studies have been dedicated to the analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome, most of which focused on missense and nonsense mutations. However, approximately 30 per cent of the SARS-CoV-2 variants are synonymous, therefore changing the target codon without affecting the corresponding protein sequence.

By performing a large-scale analysis of sequencing data generated from almost 400,000 SARS-CoV-2 samples, we show that silent mutations increasing the similarity of viral codons to the human ones tend to fixate in the viral genome overtime. This indicates that SARS-CoV-2 codon usage is adapting to the human host, likely improving its effectiveness in using the human aminoacyl-tRNA set through the accumulation of deceitfully neutral silent mutations.

One-Sentence Summary. Synonymous SARS-CoV-2 mutations related to the activity of different mutational processes may positively impact viral evolution by increasing its adaptation to the human codon usage.

Codon usage studies and epitope-based peptide vaccine prediction against Tropheryma whipplei

Background

The Tropheryma whipplei causes acute gastroenteritis to neuronal damages in Homo sapiens. Genomics and codon adaptation studies would be helpful advancements of disease evolution prediction, prevention, and treatment of disease. The codon usage data and codon usage measurement tools were deployed to detect the rare, very rare codons, and also synonymous codons usage. The higher effective number of codon usage values indicates the low codon usage bias in T. whipplei and also in the 23S and 16S ribosomal RNA genes.

Results

In T. whipplei, it was found to hold low codon biasness in genomic sets. The synonymous codons possess the base content in 3rd position that was calculated as A3S% (24.47 and 22.88), C3S% (20.99 and 22.88), T3S% (21.47 and 19.53), and G3S% (33.08 and 34.71) for 23s and 16s rRNA, respectively.

Conclusion

Amino acids like valine, aspartate, leucine, and phenylalanine hold high codon usage frequency and also found to be present in epitopes KPSYLSALSAHLNDK and FKSFNYNVAIGVRQP that were screened from proteins excinuclease ABC subunit UvrC and 3-oxoacyl-ACP reductase FabG, respectively. This method opens novel ways to determine epitope-based peptide vaccines against different pathogenic organisms.

GeneTargeter: Automated In Silico Design for Genome Editing in the Malaria Parasite, Plasmodium falciparum

Functional characterization of the multitude of poorly described proteins in the human malarial pathogen, Plasmodium falciparum, requires tools to enable genome-scale perturbation studies. Here, we present GeneTargeter (genetargeter.mit.edu), a software tool for automating the design of homology-directed repair donor vectors to achieve gene knockouts, conditional knockdowns, and epitope tagging of P. falciparum genes. We demonstrate GeneTargeter-facilitated genome-scale design of six different types of knockout and conditional knockdown constructs for the P. falciparum genome and validate the computational design process experimentally with successful donor vector assembly and transfection. The software's modular nature accommodates arbitrary destination vectors and allows customizable designs that extend the genome manipulation outcomes attainable in Plasmodium and other organisms.

A detailed comparative analysis of codon usage bias in Alongshan virus

Alongshan virus (ALSV) is an emerging tick-borne pathogen that infects humans, causing febrile disease. ALSV uses Ixodes Persulcatus ticks to infect humans with a wide range of signs, from asymptomatic to encephalitis-like syndrome. There is an increasing public health concern about the ALSV infection. To get insight into the impacts of viral relations with their hosts on viral ability, survival, and evasion from hosts immune systems remain unknown. The codon usage is a driving force in viral genome evolution; therefore, we enrolled 41 ALSV strains in codon usage analysis to elucidate the molecular evolutionary dynamics of ALSV. The results indicate that the overall codon usage among ALSV isolates is relatively similar and slightly biased. Base compositions for the cds were in order of G >A >C >U and in the third position of codons G3 >A3 >C3 >T3. The RSCU values revealed that the more frequently used codons were mostly GC ended. Different codon preferences in ALSV genes in relation to codon usage of H. sapiens and Ixodes Persulcatus genes were found. Neutrality plot was determined to reveal the superiority of natural selection over directional mutation pressure in causing CUB based on GC12 versus GC3 contents. The results of these studies suggest that the emergence of ALSV in China, Russia and Finland may also be reflected in ALSV codon usage. Altogether, the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of ALSV.

Evidence from Drosophila Supports Higher Duplicability of Faster Evolving Genes

The faster rate of evolution of duplicated genes relative to singletons has been well documented in multiple lineages. This observation has generally been attributed to a presumed release from constraint following creation of a redundant, duplicate copy. However, it is not obvious that the relationship operates in this direction. An alternative possibility—that the faster rate of evolution predates the duplication event and the observed differences result from a higher propensity to duplicate in fast-evolving genes—has been tested in primates and in insects. However, these studies arrived at different conclusions and clarity is needed on whether these contrasting results relate to differences in methodology or legitimate biological differences between the lineages selected. Here, we test whether duplicable genes are faster evolving independent of duplication in the Drosophila lineage and find that our results support the conclusion that faster evolving genes are more likely to duplicate, in agreement with previous work in primates. Our findings indicate that this characteristic of gene duplication is not restricted to a single lineage and has broad implications for the interpretation of the impact of gene duplication. We identify a subset of “singletons” which defy the general trends and appear to be faster evolving. Further investigation implicates homology detection failure and suggests that these may be duplicable genes with unidentifiable paralogs.

Comprehensive Analysis of Codon Usage Patterns in Chinese Porcine Circoviruses Based on Their Major Protein-Coding Sequences

Porcine circoviruses (PCVs) are distributed in swine herds worldwide and represent a threat to the health of domestic pigs and the profits of the swine industry. Currently, four PCV species, including PCV-1, PCV-2, PCV-3 and PCV-4, have been identified in China. Considering the ubiquitous characteristic of PCVs, the new emerged PCV-4 and the large scale of swine breeding in China, an overall analysis on codon usage bias for Chinese PCV sequences was performed by using the major proteins coding sequences (ORF1 and ORF2) to better understand the relationship of these viruses with their host. The data from genome nucleotide frequency composition and relative synonymous codon usage (RSCU) analysis revealed an overrepresentation of AT pair and the existence of a certain codon usage bias in all PCVs. However, the values of an effective number of codons (ENC) revealed that the bias was of low magnitude. Principal component analysis, ENC-plot, parity rule two analysis and correlation analysis suggested that natural selection and mutation pressure were both involved in the shaping of the codon usage patterns of PCVs. However, a neutrality plot revealed a stronger effect of natural selection than mutation pressure on codon usage patterns. Good host adaptation was also shown by the codon adaptation index analysis for all these viruses. Interestingly, obtained data suggest that PCV-4 might be more adapted to its host compared to other PCVs. The present study obtained insights into the codon usage pattern of PCVs based on ORF1 and ORF2, which further helps the understanding the molecular evolution of these swine viruses.

Exploiting the Gal4/UAS System as Plant Orthogonal Molecular Toolbox to Control Reporter Expression in Arabidopsis Protoplasts

The ability of protein domains to fold independently from the rest of the polypeptide is the principle governing the generation of fusion proteins with customized functions. A clear example is the split transcription factor system based on the yeast GAL4 protein and its cognate UAS enhancer. The rare occurrence of the UAS element in the transcriptionally sensitive regions of the Arabidopsis genome makes this transcription factor an ideal orthogonal platform to control reporter induction. Moreover, heterodimeric transcriptional complexes can be generated by exploiting posttranslational modifications hampering or promoting the interaction between GAL4-fused transcriptional partners, whenever this leads to the reconstitution of a fully functional GAL4 factor.The assembly of multiple engineered proteins into a synthetic transcriptional complex requires preliminary testing, before its components can be stably introduced into the plant genome. Mesophyll protoplast transformation represents a fast and reliable technique to test and optimize synthetic regulatory modules. Remarkable properties are the possibility to transform different combinations of plasmids (co-transformation) and the physiological resemblance of these isolated cells with the original tissue.Here we describe an extensive protocol to produce and exploit Arabidopsis mesophyll protoplasts to investigate the transcriptional output of GAL4/UAS-based complexes that are sensitive to posttranslational protein modifications.

In silico proof of principle of machine learning-based antibody design at unconstrained scale

Generative machine learning (ML) has been postulated to become a major driver in the computational design of antigen-specific monoclonal antibodies (mAb). However, efforts to confirm this hypothesis have been hindered by the infeasibility of testing arbitrarily large numbers of antibody sequences for their most critical design parameters: paratope, epitope, affinity, and developability. To address this challenge, we leveraged a lattice-based antibody-antigen binding simulation framework, which incorporates a wide range of physiological antibody-binding parameters. The simulation framework enables the computation of synthetic antibody-antigen 3D-structures, and it functions as an oracle for unrestricted prospective evaluation and benchmarking of antibody design parameters of ML-generated antibody sequences. We found that a deep generative model, trained exclusively on antibody sequence (one dimensional: 1D) data can be used to design conformational (three dimensional: 3D) epitope-specific antibodies, matching, or exceeding the training dataset in affinity and developability parameter value variety. Furthermore, we established a lower threshold of sequence diversity necessary for high-accuracy generative antibody ML and demonstrated that this lower threshold also holds on experimental real-world data. Finally, we show that transfer learning enables the generation of high-affinity antibody sequences from low-N training data. Our work establishes a priori feasibility and the theoretical foundation of high-throughput ML-based mAb design.

KIF3B gene silent variant leading to sperm morphology and motility defects and male infertility

Sperm structural and functional defects are leading causes of male infertility. Patients with immotile sperm disorders suffer from axoneme failure and show a significant reduction in sperm count. The kinesin family member 3B (KIF3B) is one of the genes involved in the proper formation of sperm with a critical role in intraflagellar and intramanchette transport. A part of exon 2 and exons 3–5 of the KIF3B encodes a protein coiled-coil domain that interacts with IFT20 from the IFT protein complex. In the present study, the coding region of KIF3B coiled-coil domain was assessed in 88 oligoasthenoteratozoospermic patients, and the protein expression was evaluated in the mature spermatozoa of the case and control groups using immunocytochemistry and western blotting. According to the results, there was no genetic variation in the exons 3–5 of the KIF3B, but a new A &gt; T variant was identified within the exon 2 in 30 patients, where nothing was detected in the control group. In contrast to healthy individuals, significantly reduced protein expression was observable in oligoasthenoteratozoospermic (OAT) patients carrying variation where protein organization was disarranged, especially in the principal piece and midpiece of the sperm tail. Besides, the protein expression level was lower in the patients’ samples compared to that of the control group. According to the results of the present study the NM_004798.3:c.1032A &gt; T, p.Pro344 = variant; which has been recently submitted to the Clinvar database; although synonymous, causes alterations in the transcription factor binding site, exon skipping, and also exonic splicing enhancer-binding site. Therefore, KIF3B can play an important role in spermatogenesis and the related protein reduction can cause male infertility.

The N-Acetylglucosamine Kinase from Yarrowia lipolytica Is a Moonlighting Protein

In Yarrowia lipolytica, expression of the genes encoding the enzymes of the N-acetylglucosamine (NAGA) utilization pathway (NAG genes) becomes independent of the presence of NAGA in a Ylnag5 mutant lacking NAGA kinase. We addressed the question of whether the altered transcription was due to a lack of kinase activity or to a moonlighting role of this protein. Glucosamine-6-phosphate deaminase (Nag1) activity was measured as a reporter of NAG genes expression. The NGT1 gene encoding the NAGA transporter was deleted, creating a Ylnag5 ngt1 strain. In glucose cultures of this strain, Nag1 activity was similar to that of the Ylnag5 strain, ruling out the possibility that NAGA derived from cell wall turnover could trigger the derepression. Heterologous NAGA kinases were expressed in a Ylnag5 strain. Among them, the protein from Arabidopsis thaliana did not restore kinase activity but lowered Nag1 activity 4-fold with respect to a control. Expression in the Ylnag5 strain of YlNag5 variants F320S or D214V with low kinase activity caused a repression similar to that of the wild-type protein. Together, these results indicate that YlNag5 behaves as a moonlighting protein. An RNA-seq analysis revealed that the Ylnag5 mutation had a limited transcriptomic effect besides derepression of the NAG genes.