Codon usage bias: causative factors, quantification methods and genome-wide patterns: with emphasis on insect genomes

The Patterns of Codon Usage between Chordates and Arthropods are Different but Co-evolving with Mutational Biases

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.

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

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

Investigation of CaMV-host co-evolution through synonymous codon pattern

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

Fitness difference between two synonymous mutations of Phytophthora infestans ATP6 gene

BACKGROUND

Sequence variation produced by mutation provides the ultimate source of natural selection for species adaptation. Unlike nonsynonymous mutation, synonymous mutations are generally considered to be selectively neutral but accumulating evidence suggests they also contribute to species adaptation by regulating the flow of genetic information and the development of functional traits. In this study, we analysed sequence characteristics of ATP6, a housekeeping gene from 139 Phytophthora infestans isolates, and compared the fitness components including metabolic rate, temperature sensitivity, aggressiveness, and fungicide tolerance among synonymous mutations.

RESULTS

We found that the housekeeping gene exhibited low genetic variation and was represented by two major synonymous mutants at similar frequency (0.496 and 0.468, respectively). The two synonymous mutants were generated by a single nucleotide substitution but differed significantly in fitness as well as temperature-mediated spatial distribution and expression. The synonymous mutant ending in AT was more common in cold regions and was more expressed at lower experimental temperature than the synonymous mutant ending in GC and vice versa.

CONCLUSION

Our results are consistent with the argument that synonymous mutations can modulate the adaptive evolution of species including pathogens and have important implications for sustainable disease management, especially under climate change.

Comparative Analysis and Phylogenetic Study of Dawkinsia filamentosa and Pethia nigrofasciata Mitochondrial Genomes

Smiliogastrinae are recognized for their high nutritional and ornamental value. In this study, we employed high-throughput sequencing technology to acquire the complete mitochondrial genome sequences of Dawkinsia filamentosa and Pethia nigrofasciata. The gene composition and arrangement order in these species were similar to those of typical vertebrates, comprising 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 non-coding region. The mitochondrial genomes of D. filamentosa and P. nigrofasciata measure 16,598 and 16,948 bp, respectively. Both D. filamentosa and P. nigrofasciata exhibit a significant preference for AT bases and an anti-G bias. Notably, the AT and GC skew values of the ND6 gene fluctuated markedly, suggesting that the selection and mutation pressures on this gene may differ from those affecting other genes. Phylogenetic analysis, based on the complete mitochondrial genomes of 23 Cyprinidae fishes, revealed that D. filamentosa is closely related to the sister group comprising Dawkinsia denisonii and Sahyadria chalakkudiensis. Similarly, P. nigrofasciata forms a sister group with Pethia ticto and Pethia stoliczkana.

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

Splicing neoantigen discovery with SNAF reveals shared targets for cancer immunotherapy

Immunotherapy has emerged as a crucial strategy to combat cancer by "reprogramming" a patient's own immune system. Although immunotherapy is typically reserved for patients with a high mutational burden, neoantigens produced from posttranscriptional regulation may provide an untapped reservoir of common immunogenic targets for new targeted therapies. To comprehensively define tumor-specific and likely immunogenic neoantigens from patient RNA-Seq, we developed Splicing Neo Antigen Finder (SNAF), an easy-to-use and open-source computational workflow to predict splicing-derived immunogenic MHC-bound peptides (T cell antigen) and unannotated transmembrane proteins with altered extracellular epitopes (B cell antigen). This workflow uses a highly accurate deep learning strategy for immunogenicity prediction (DeepImmuno) in conjunction with new algorithms to rank the tumor specificity of neoantigens (BayesTS) and to predict regulators of mis-splicing (RNA-SPRINT). T cell antigens from SNAF were frequently evidenced as HLA-presented peptides from mass spectrometry (MS) and predict response to immunotherapy in melanoma. Splicing neoantigen burden was attributed to coordinated splicing factor dysregulation. Shared splicing neoantigens were found in up to 90% of patients with melanoma, correlated to overall survival in multiple cancer cohorts, induced T cell reactivity, and were characterized by distinct cells of origin and amino acid preferences. In addition to T cell neoantigens, our B cell focused pipeline (SNAF-B) identified a new class of tumor-specific extracellular neoepitopes, which we termed ExNeoEpitopes. ExNeoEpitope full-length mRNA predictions were tumor specific and were validated using long-read isoform sequencing and in vitro transmembrane localization assays. Therefore, our systematic identification of splicing neoantigens revealed potential shared targets for therapy in heterogeneous cancers.

Genome-wide phylogenetic and genetic evolutionary analyses of mitochondria in Hypoderma bovis and H. sinense on the Qinghai-Tibetan Plateau

Hypoderma bovis (H. bovis) and Hypoderma sinense (H. sinense) are insects that cause hypodermosis in yaks and Bos taurus. Hypodermosis is a severe skin condition that not only impairs the development of local animal husbandry but also poses threats to human health as a zoonosis. The Qinghai-Tibetan Plateau (QTP) is known as the "Roof of the World." Its unique geographical environment and climate conditions have supported the growth of a wide range of mammals, providing favorable conditions for Hypoderma spp. to complete their life cycles. In this study, the whole mitochondrial genomes of H. bovis and H. sinense collected from the QTP were sequenced and phylogenetically analyzed. We found that the whole genomes of H. bovis and H. sinense are 16,283 bp and 16,300 bp in length, respectively. Both the H. bovis and H. sinense genomes have 37 mitochondrial genes, which include two rRNA genes (16S rRNA and 12S rRNA), 22 tRNA genes, the control region (D-loop region), the light chain replication initiation region, and 13 protein-coding genes (PCGs). The phylogenetic tree generated based on the 13 PCGs revealed close phylogenetic relationships between H. sinense, H. bovis, and Hypoderma lineatum. A similar result was also found in our phylogenetic analysis based on 18S rRNA and 28S rRNA. However, analysis of cytochrome oxidase subunit I (COI) showed cluster of H. bovis, H. sinense, and Cuterebra spp. on the same branch, all belonging to Oestridae. The differentiation time generated based on 13 PCGs indicates that H. bovis and H. sinense differentiated and formed ~4.69 million years ago (Mya) and ~4.06 Mya, respectively. This timing coincides with the differentiation and appearance of yak and Bos taurus in the Pliocene (~4.7 Mya), indicating that the parasites and mammals diverged in close temporal proximity. Of note, this period also witnessed a rapid uplift of the QTP, causing significant climate and environmental changes. Thus, we conjecture that the differentiation of Hypoderma spp. is potentially related to the differentiation of their host species, as well as climate changes caused by the uplift of the QTP. Overall, our study can provide valuable data to support further studies on the phylogeny and differentiation of Hypoderma spp. on the QTP.

Adaptive evidence of mitochondrial genes in Pteromalidae and Eulophidae (Hymenoptera: Chalcidoidea)

Pteromalidae and Eulophidae are predominant and abundant taxa within Chalcidoidea (Hymenoptera: Apocrita). These taxa are found in diverse ecosystems, ranging from basin deserts (200 m) to alpine grasslands (4500 m). Mitochondria, cellular powerhouses responsible for energy production via oxidative phosphorylation, are sensitive to various environmental factors such as extreme cold, hypoxia, and intense ultraviolet radiation characteristic of alpine regions. Whether the molecular evolution of mitochondrial genes in these parasitoids corresponds to changes in the energy requirements and alpine environmental adaptations remains unknown. In this study, we performed a comparative analysis of mitochondrial protein-coding genes from 11 alpine species of Pteromalidae and Eulophidae, along with 18 lowland relatives, including 16 newly sequenced species. We further examined the codon usage preferences (RSCU, ENC-GC3s, neutrality, and PR2 bias plot) in these mitochondrial protein-coding sequences and conducted positive selection analysis based on their Bayesian phylogenetic relationships, and identified positive selection sites in the ATP6, ATP8, COX1, COX3, and CYTB genes, emphasizing the crucial role of mitochondrial gene adaptive evolution in the adaptation of Pteromalidae and Eulophidae to alpine environments. The phylogenetically independent contrast (PIC) analysis results verified the ω ratio of 13 PCGs from Pteromalidae and Eulophidae increased with elevation, and results from generalized linear model confirm that ATP6, ATP8, COX3, and ND1 are closely correlated with temperature-related environmental factors. This research not only enriched the molecular data of endemic alpine species but also underscores the significance of mitochondrial genes in facilitating the adaptation of these minor parasitoids to plateau habitats.

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

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

Risk assessment and dissemination mechanism of antibiotic resistance genes in compost

In recent years, the excessive of antibiotics in livestock and poultry husbandry, stemming from extensive industry experience, has resulted in the accumulation of residual antibiotics and antibiotic resistance genes (ARGs) in livestock manure. Composting, as a crucial approach for the utilization of manure resources, has the potential to reduce the levels of antibiotics and ARGs in manure, although complete elimination is challenging. Previous studies have primarily focused on the diversity and abundance of ARGs in compost or have solely examined the correlation between ARGs and their carriers, potentially leading to a misjudgment of the actual risk associated with ARGs in compost. To address this gap, this study investigated the transfer potential of ARGs in compost and their co-occurrence with opportunistic pathogenic bacteria by extensively analyzing metagenomic sequencing data of compost worldwide. The results demonstrated that the potential risk of ARGs in compost was significantly lower than in manure, suggesting that composting effectively reduces the risk of ARGs. Further analysis showed that the microbes shifted their life history strategy in manure and compost due to antibiotic pressure and formed metabolic interactions dominated by antibiotic-resistant microbes, increasing ARG dissemination frequency. Therefore, husbandry practice without antibiotic addition was recommended to control ARG evolution, dissemination, and abatement both at the source and throughout processing.

Description and phylogenetic analysis of the complete mitochondrial genome in Eulaelaps silvestris provides new insights into the molecular classification of the family Haemogamasidae

In this study, the mitochondrial genome of Eulaelaps silvestris, which parasitizes Apodemus chevrieri, was sequenced and assembled to fill the gap in understanding the molecular evolution of the genus Eulaelaps. The E. silvestris mitochondrial genome is a double-stranded DNA molecule with a length of 14 882 bp, with a distinct AT preference for base composition and a notably higher AT content than GC content. The arrangement between genes is relatively compact, with a total of 10 gene intergenic regions and 12 gene overlap regions. All protein-coding genes had a typical ATN initiation codon, and only 2 protein-coding genes had an incomplete termination codon T. Out of the 13 protein-coding genes, the 5 most frequently used codons ended in A/U, with only 1 codon ending in G/C had an relative synonymous codon usage value >1. Except for trnS1 and trnS2, which lacked the D arm, all other tRNAs were able to form a typical cloverleaf structure; and there were a total of 38 mismatches in the folding process of tRNA genes. Unlike the gene arrangement order of the arthropod hypothetical ancestor, the E. silvestris mitochondrial genome underwent fewer rearrangements, mainly near tRNA genes and control regions. Both the maximum likelihood tree and the Bayesian tree showed that the family Haemogamasidae is most closely related to the family Dermanyssidae. The results not only provide a theoretical basis for studying the phylogenetic relationships of the genus Eulaelaps, but also provide molecular evidence that the family Haemogamasidae does not belong to the subfamily Laelapidae.

Misclassified: identification of zoonotic transition biomarker candidates for influenza A viruses using deep neural network

Introduction: Zoonotic transition of Influenza A viruses is the cause of epidemics with high rates of morbidity and mortality. Predicting which viral strains are likely to transition from their genetic sequence could help in the prevention and response against these zoonotic strains. We hypothesized that features predictive of viral hosts could be leveraged to identify biomarkers of zoonotic viral transition. Methods: We trained deep learning models to predict viral hosts based on the virus mRNA or protein sequences. Our multi-host dataset contained 848,630 unique nucleotide sequences obtained from the NCBI Influenza Virus and Influenza Research Databases. Each sequence, representing one gene from one viral strain, was classified into one of the three host categories: Avian, Human, and Swine. Trained models were analyzed using various neural network interpretation methods to identify interesting candidates for zoonotic transition biomarkers. Results: Using mRNA sequences as input led to higher prediction accuracies than amino acids, suggesting that the codon sequence contains information relevant to viral hosts that is lost during protein translation. UMAP visualization of the latent space of our classifiers showed that viral sequences clustered according to their host of origin. Interestingly, sequences from pandemic zoonotic viral strains localized at the margins between hosts, while zoonotic sequences incapable of Human-to-Human transmission localized with non-zoonotic viruses from the same host. In addition, host prediction for pandemic zoonotic sequences had low prediction accuracy, which was not the case for the other zoonotic strains. This supports our hypothesis that ambiguously predicted viral sequences bear features associated with cross-species infectivity. Finally, we compared misclassified sequences to well-classified ones to extract interesting candidates for zoonotic transition biomarkers. While features varied significantly between pairs of species and viral genes, several codons were conserved in Swine-to-Human and Avian-to-Human misclassified sequences, and in particular in the NA, HA, and NP genes, suggesting their importance for zoonosis in Humans. Discussion: Analysis of viral sequences using neural network interpretation approaches revealed important genetic differences between zoonotic viruses with pandemic potential, compared to non-zoonotic viral strains or zoonotic viruses incapable of Human-to-Human transmission.

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.

Uncovering the first complete chloroplast genomics, comparative analysis, and phylogenetic relationships of the medicinal plants Rhamnus cathartica and Frangula alnus (Rhamnaceae)

Rhamnus cathartica and Frangula alnus are economically valuable medicinal plants from the Rhamnaceae family. However, their chloroplast genome structure, phylogenetic position, relationships, and evolution remain poorly understood. Herein, the complete chloroplast genome resources of R. cathartica and F. alnus have been added. The first comparative analysis of the Rhamnus and Frangula species based on complete chloroplast genomes was provided. The chloroplast genomes of R. cathartica and F. alnus exhibited a quadripartite structure, with total lengths of 161,149 bp and 161,255 bp, respectively. The lack of the infA and psbL genes does not negatively impact the normal functioning of Rhamnus and Frangula species. The rpl20 and rpl33 genes are undergoing rapid evolution. Rhamnus and Frangula species prefer amino acids with A/U-terminal codons. There were between 100 and 126 simple sequence repeats and between 38 and 100 long repeats. Several highly divergent intergenic regions (trnK-UUU-trnQ-UUG, atpH-atpI, trnY-GUA-trnE-UUC, trnG-GCC-trnfM-CAU, trnT-UGU-trnF-GAA, rpl20-rps12, and rpl22-rps19) and highly divergent genes (ycf3, ndhA, rpl32, and ycf1) were identified, which could serve as potential phylogenetic markers due to their variability. We reconstructed the phylogenetic relationships among Rhamnus species and F. alnus using complete chloroplast genomes. There is no significant correlation between the medicinal value of the species analyzed and their phylogenetic relationships. These results provide valuable insights for understanding the phylogenetic relationship and evolution of Rhamnus and Frangula species. These findings could serve as a foundation for future studies on the Rhamnaceae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01331-7.

Analysis of synonymous codon usage patterns in mitochondrial genomes of nine Amanita species

Introduction

Codon basis is a common and complex natural phenomenon observed in many kinds of organisms.

Methods

In the present study, we analyzed the base bias of 12 mitochondrial core protein-coding genes (PCGs) shared by nine Amanita species.

Results

The results showed that the codons of all Amanita species tended to end in A/T, demonstrating the preference of mitochondrial codons of Amanita species for a preference for this codon. In addition, we detected the correlation between codon base composition and the codon adaptation index (CAI), codon bias index (CBI), and frequency of optimal codons (FOP) indices, indicating the influence of base composition on codon bias. The average effective number of codons (ENC) of mitochondrial core PCGs of Amanita is 30.81, which is <35, demonstrating the strong codon preference of mitochondrial core PCGs of Amanita. The neutrality plot analysis and PR2-Bias plot analysis further demonstrated that natural selection plays an important role in Amanita codon bias. In addition, we obtained 5–10 optimal codons (ΔRSCU > 0.08 and RSCU > 1) in nine Amanita species, and GCA and AUU were the most widely used optimal codons. Based on the combined mitochondrial sequence and RSCU value, we deduced the genetic relationship between different Amanita species and found large variations between them.

Discussion

This study promoted the understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Folding of heterologous proteins in bacterial cell factories: Cellular mechanisms and engineering strategies

The expression of correctly folded and functional heterologous proteins is important in many biotechnological production processes, whether it is enzymes, biopharmaceuticals or biosynthetic pathways for production of sustainable chemicals. For industrial applications, bacterial platform organisms, such as E. coli, are still broadly used due to the availability of tools and proven suitability at industrial scale. However, expression of heterologous proteins in these organisms can result in protein aggregation and low amounts of functional protein. This review provides an overview of the cellular mechanisms that can influence protein folding and expression, such as co-translational folding and assembly, chaperone binding, as well as protein quality control, across different model organisms. The knowledge of these mechanisms is then linked to different experimental methods that have been applied in order to improve functional heterologous protein folding, such as codon optimization, fusion tagging, chaperone co-production, as well as strain and protein engineering strategies.

A comprehensive analysis of Usutu virus (USUV) genomes revealed lineage-specific codon usage patterns and host adaptations

The Usutu virus (USUV) is an emerging arbovirus virus maintained in the environment of Afro-Eurasia via a bird-mosquito-bird enzootic cycle and sporadically infected other vertebrates. Despite primarily asymptomatic or mild symptoms, humans infected by USUV can develop severe neurological diseases such as meningoencephalitis. However, no detailed study has yet been conducted to investigate its evolution from the perspective of codon usage patterns. Codon usage choice of viruses reflects the genetic variations that enable them to reconcile their viability and fitness toward the external environment and new hosts. This study performed a comprehensive evolution and codon usage analysis of USUVs. Our reconstructed phylogenetic tree confirmed that the circulation viruses belong to eight distinct lineages, reaffirmed by principal component analysis based on codon usage patterns. We also found a relatively small codon usage bias and that natural selection, mutation pressure, dinucleotide abundance, and evolutionary processes collectively shaped the codon usage of the USUV, with natural selection predominating over the others. Additionally, a complex interaction of codon usage between the USUV and its host was observed. This process could have enabled USUV to adapt to various hosts and vectors, including humans. Therefore, the USUV may possess a potential risk of cross-species transmission and subsequent outbreaks. In this respect, further epidemiologic surveys, diversity monitoring, and pathogenetic research are warranted.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

The identification of conserved sequence features of co-translationally decayed mRNAs and upstream open reading frames in angiosperm transcriptomes

RNA turnover is essential in maintaining messenger RNA (mRNA) homeostasis during various developmental stages and stress responses. Co-translational mRNA decay (CTRD), a process in which mRNAs are degraded while still associated with translating ribosomes, has recently been discovered to function in yeast and three angiosperm transcriptomes. However, it is still unclear how prevalent CTRD across the plant lineage. Moreover, the sequence features of co-translationally decayed mRNAs have not been well-studied. Here, utilizing a collection of publicly available degradome sequencing datasets for another seven angiosperm transcriptomes, we have confirmed that CTRD is functioning in at least 10 angiosperms and likely throughout the plant lineage. Additionally, we have identified sequence features shared by the co-translationally decayed mRNAs in these species, implying a possible conserved triggering mechanism for this pathway. Given that degradome sequencing datasets can also be used to identify actively translating upstream open reading frames (uORFs), which are quite understudied in plants, we have identified numerous actively translating uORFs in the same 10 angiosperms. These findings reveal that actively translating uORFs are prevalent in plant transcriptomes, some of which are conserved across this lineage. We have also observed conserved sequence features in the regions flanking these uORFs' stop codons that might contribute to ribosome stalling at these sequences. Finally, we discovered that there were very few overlaps between the mRNAs harboring actively translating uORFs and those sorted into the co-translational decay pathway in the majority of the studied angiosperms, suggesting that these two processes might be nearly mutually exclusive in those species. In total, our findings provide the identification of CTRD and actively translating uORFs across a broad collection of plants and provide novel insights into the important sequence features associated with these collections of mRNAs and regulatory elements, respectively.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

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

Methods

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

Results

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

Discussion

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

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

First Description of the Mitogenome Features of Neofoleyellides Genus (Nematoda: Onchocercidae) Isolated from a Wild Bird (Pyrrhocorax pyrrhocorax)

Simple Summary

Filarioidea, a superfamily of nematodes, presently includes 42 species divided into six genera, mainly in the family Onchocercidae, which have been reported to infect a wide range of hosts, including reptiles, birds, and mammals. Current limitations in molecular characterization methods and species identification are the main obstacles to a better understanding of the biology of Onchocercidae species, particularly in wildlife. Thus, the objective of the present study was to sequence and analyze the complete mt genome of Neofoleyellides sp. isolated from a wild bird (Pyrrhocorax pyrrhocorax) and to assess its phylogenetic position in the Onchocercidae family. The evaluated Neofoleyellides sp. mt genome was consistent with the molecular pattern of the Onchocercidae family: 36 subunits consisting of 12 PCGs, 2 rRNAs, and 22 tRNAs. Phylogenetic analyses based on the 18S rRNA gene, cox1 gene, and 12 PCGs showed consistent results, which strongly supported monophyly of the genus Neofoleyellides. These findings enriched the gene database and improved our knowledge of the molecular characteristics of the Onchocercidae family, which provide useful genetic markers to study the population genetics, molecular biology, and phylogenetics of these Onchocercidae nematodes.

Abstract

The Onchocercidae family is composed of more than 30 valid nematode species with notable zoonotic potential. Current limitations in molecular characterization methods and species identification are the main obstacles to a better understanding of the biology of Onchocercidae species, particularly in wildlife. This study describes for the first time the complete mitochondrial (mt) genome sequence of Neofoleyellides sp. isolated from a wild bird (Pyrrhocorax pyrrhocorax) and belonging to the Neofoleyellides genus (Nematoda: Onchocercidae). The mt genome of Neofoleyellides sp. (GenBank accession number: ON641583) was a typical circular DNA molecule of 13,628 bp in size with an AT content of 76.69%. The complete mt genome comprised 36 functional subunits, including 12 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. The most common start codon was ATT/ATG except for nad2 with TTG, and TAA was the termination codon for all protein-coding genes (PCGs). Phylogenetic analysis of the concatenated and aligned amino acid sequences of the 12 PCGs showed that the trees generated using different methods (Bayesian inference and maximum likelihood) with different partition schemes shared similar topologies. The isolated Neofoleyellides sp. was placed in the Onchocercidae family and formed a sister branch with the genera Onchocerca and Dirofilaria. The entire mt genome of Neofoleyellides sp. presented in this study could provide useful data for studying the population genetics and phylogenetic relationships of Onchocercidae species.

Analysis of the Compositional Features and Codon Usage Pattern of Genes Involved in Human Autophagy

Autophagy plays an intricate role in paradigmatic human pathologies such as cancer, and neurodegenerative, cardiovascular, and autoimmune disorders. Autophagy regulation is performed by a set of autophagy-related (ATG) genes, first recognized in yeast genome and subsequently identified in other species, including humans. Several other genes have been identified to be involved in the process of autophagy either directly or indirectly. Studying the codon usage bias (CUB) of genes is crucial for understanding their genome biology and molecular evolution. Here, we examined the usage pattern of nucleotide and synonymous codons and the influence of evolutionary forces in genes involved in human autophagy. The coding sequences (CDS) of the protein coding human autophagy genes were retrieved from the NCBI nucleotide database and analyzed using various web tools and software to understand their nucleotide composition and codon usage pattern. The effective number of codons (ENC) in all genes involved in human autophagy ranges between 33.26 and 54.6 with a mean value of 45.05, indicating an overall low CUB. The nucleotide composition analysis of the autophagy genes revealed that the genes were marginally rich in GC content that significantly influenced the codon usage pattern. The relative synonymous codon usage (RSCU) revealed 3 over-represented and 10 under-represented codons. Both natural selection and mutational pressure were the key forces influencing the codon usage pattern of the genes involved in human autophagy.

First description of the mitogenome and phylogeny:Aedes vexansand Ochlerotatus caspius of the Tribe Aedini (Diptera: Culicidae)

Culicidae, the mosquito family, includes more than 3600 species subdivided into the subfamilies Anophelinae and Culicinae. One-third of mosquitoes belong to the Aedini tribe, which is subordinate to the subfamily Culicinae, which comprises common vectors of viral zoonoses. The tribe of Aedini is extremely diverse in morphology and geographical distribution and has high ecological and medical significance. However, knowledge about the systematics of the Aedini tribe is still limited owing to its large population and the similar morphological characteristics of its species. This study provides the first description of the complete mitochondrial (mt) genome sequence of Aedes vexans and Ochlerotatus caspius belonging to the Aedini tribe. The mt genomes of A. vexans and O. caspius are circular molecules that are 15,861 bp and 15,954 bp in size, with AT contents of 78.54% and 79.36%, respectively. Both the circular mt genomes comprise 37 functional subunits, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a control region (also known as the AT-rich region). The most common start codons are ATT/ATG, apart from cox1 (TCG) and nad5 (GTG), while TAA is the termination codon for all PCGs. All tRNAs have a typical clover leaf structure, except tRNA Ser1. Phylogenetic analysis of the concatenated, aligned amino acid sequences of the 13 PCGs showed that A. vexans gathered with Aedes sp. in a sister taxon, and O. caspius gathered with Ochlerotatus sp. in a sister taxon. The findings from the present study support the concept of monophyly of all groups, ratify the current taxonomic classification, and provide vital molecular marker resources for further studies of the taxonomy, population genetics, and systematics of the Aedini tribe.

Using codon usage analysis to speculate potential animal hosts of hepatitis E virus: An exploratory study

There has been an increase in the reported number of animals worldwide that carry the hepatitis E virus (HEV). This study aimed to explore potential animal hosts for HEV through codon usage analysis. Full-length HEV sequences of six genotypes as well as codon usage of potential animal hosts were collected. Moreover, nucleotide composition and codon usage bias were compared across HEV genotypes and animal hosts. Based on the analysis for human HEV-1 and humans, the results were basically consistent with epidemiology evidence. Among 17 potential animal hosts, all HEV genotypes exhibited a preference for guanine/cytosine in the third position of synonymous codons. Furthermore, non-human primates and humans have large high-frequency codons identical to HEV in addition to a high correlation of codon fraction with HEV. Some animals in close contact with humans showed high preference for HEV, including cattle, dogs, and rats with HEV-A, cats, dogs, and swine with HEV-C1. Codon usage bias has limited efficiency in determining the hosts for HEV, but it may provide indicative clues for potential animal hosts when combined with experimental and epidemiological evidence.

The genetic evolution and codon usage pattern of severe fever with thrombocytopenia syndrome virus

Severe fever with thrombocytopenia syndrome (SFTS) is a newly emerging zoonotic infectious disease caused by the SFTS virus (SFTSV), which has been continuously circulating in Eastern Asia in recent years. Although the evolution of SFTSV has been investigated, the evolutionary changes associated with codon usage have not been reported. Thus, a comprehensive genetic and codon usage bias analysis of SFTSV was conducted to elucidate the genetic diversity and evolutionary relationships in a novel perspective. The study amplified and sequenced fifteen SFTSV strains from a prefecture of Zhejiang Province, Eastern China in 2020, where SFTS cases have been continuously reported in the past decade. Phylogenetic analysis was conducted based on the complete coding sequences of SFTSV segments. It suggested that all SFTSV strains circulating in Zhejiang were clustered with Japanese and Korean strains, which belonged to two different genotypes. Meanwhile, thirty-nine genetic reassortants classified into nineteen different reassortment forms were identified, while 45 recombination events in 41 SFTSV strains were found. Codon usage patterns were further analyzed to understand the evolutionary changes in relation to genotype and host. And it revealed that codon usage bias was mainly driven by natural selection rather than mutation pressure. In addition, the codon adaptation index (CAI) analysis demonstrated the strong adaptability of SFTSV to Gallus gallus and Homo sapiens. Similarity index (SiD) analysis indicated that Haemaphysalis longicornis posed a strong selection pressure to SFTSV. In conclusion, this study revealed that the genetic diversity of SFTSV is gradually increasing. The codon usage analysis suggested that codon usage bias of SFTSV was mainly driven by natural selection, and SFTSV has evolved host-specific codon usage patterns. This contributes to the development of control measures against SFTSV.

Calculating and comparing codon usage values in rare disease genes highlights codon clustering with disease-and tissue- specific hierarchy

We designed a novel strategy to define codon usage bias (CUB) in 6 specific small cohorts of human genes. We calculated codon usage (CU) values in 29 non-disease-causing (NDC) and 31 disease-causing (DC) human genes which are highly expressed in 3 distinct tissues, kidney, muscle, and skin. We applied our strategy to the same selected genes annotated in 15 mammalian species. We obtained CUB hierarchical clusters for each gene cohort which showed tissue-specific and disease-specific CUB fingerprints. We showed that DC genes (especially those expressed in muscle) display a low CUB, well recognizable in codon hierarchical clustering. We defined the extremely biased codons as "zero codons" and found that their number is significantly higher in all DC genes, all tissues, and that this trend is conserved across mammals. Based on this calculation in different gene cohorts, we identified 5 codons which are more differentially used across genes and mammals, underlining that some genes have favorite synonymous codons in use. Since of the muscle genes clear clusters, and, among these, dystrophin gene surprisingly does not show any "zero codon" we adopted a novel approach to study CUB, we called "mapping-on-codons". We positioned 2828 dystrophin missense and nonsense pathogenic variations on their respective codon, highlighting that its frequency and occurrence is not dependent on the CU values. We conclude our strategy consents to identify a hierarchical clustering of CU values in a gene cohort-specific fingerprints, with recognizable trend across mammals. In DC muscle genes also a disease-related fingerprint can be observed, allowing discrimination between DC and NDC genes. We propose that using our strategy which studies CU in specific gene cohorts, as rare disease genes, and tissue specific genes, may provide novel information about the CUB role in human and medical genetics, with implications on synonymous variations interpretation and codon optimization algorithms.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

Unlocking the Complete Chloroplast Genome of a Native Tree Species from the Amazon Basin, Capirona (Calycophyllum Spruceanum, Rubiaceae), and Its Comparative Analysis with Other Ixoroideae Species

Capirona (Calycophyllum spruceanum Benth.) belongs to subfamily Ixoroideae, one of the major lineages in the Rubiaceae family, and is an important timber tree. It originated in the Amazon Basin and has widespread distribution in Bolivia, Peru, Colombia, and Brazil. In this study, we obtained the first complete chloroplast (cp) genome of capirona from the department of Madre de Dios located in the Peruvian Amazon. High-quality genomic DNA was used to construct libraries. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of C. spruceanum has a 154,480 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (84,813 bp) and a small single-copy (SSC) region (18,101 bp), separated by two inverted repeat (IR) regions (25,783 bp). The annotation of C. spruceanum cp genome predicted 87 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and one pseudogene. A total of 41 simple sequence repeats (SSR) of this cp genome were divided into mononucleotides (29), dinucleotides (5), trinucleotides (3), and tetranucleotides (4). Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other six Ixoroideae species revealed that the small single copy and large single copy regions showed more divergence than inverted regions. Finally, phylogenetic analyses resolved that C. spruceanum is a sister species to Emmenopterys henryi and confirms its position within the subfamily Ixoroideae. This study reports for the first time the genome organization, gene content, and structural features of the chloroplast genome of C. spruceanum, providing valuable information for genetic and evolutionary studies in the genus Calycophyllum and beyond.

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.

In-depth analysis of amino acid and nucleotide sequences of Hsp60: how conserved is this protein?

Chaperonin Hsp60, as a protein found in all organisms, is of great interest in medicine, since it is present in many tissues and can be used both as a drug and as an object of targeted therapy. Hence, Hsp60 deserves a fundamental comparative analysis to assess its evolutionary characteristics. It was found that the percent identity of Hsp60 amino acid sequences both within and between phyla was not high enough to identify Hsp60s as highly conserved proteins. However, their ATP binding sites are largely conserved. The amino acid composition of Hsp60s remained relatively constant. At the same time, the analysis of the nucleotide sequences showed that GC content in the Hsp60 genes was comparable to or greater than the genomic values, which may indicate a high resistance to mutations due to tight control of the nucleotide composition by DNA repair systems. Natural selection plays a dominant role in the evolution of Hsp60 genes. The degree of mutational pressure affecting the Hsp60 genes is quite low, and its direction does not depend on taxonomy. Interestingly, for the Hsp60 genes from Chordata, Arthropoda, and Proteobacteria the exact direction of mutational pressure could not be determined. However, upon further division into classes, it was found that the direction of the mutational pressure for Hsp60 genes from Fish differs from that for other chordates. The direction of the mutational pressure affects the synonymous codon usage bias. The number of high and low represented codons increases with increasing GC content, which can improve codon usage. Special server has been created for bioinformatics analysis of Hsp60: http://oka.protres.ru:4202/.

Compositional features and pattern of codon usage for mitochondrial CO genes among reptiles

The phenomenon of non-random occurrence of synonymous nucleotide triplets (codons) in the coding sequences of genes is the codon usage bias (CUB). In this study, we used bioinformatic tool kit to analyze the compositional pattern and CUB of mitogenes namely COI, COII and COIII across different orders of reptiles. Estimation of overall base composition in the protein-coding sequences of COI, COII and COIII genes of the reptilian orders revealed an uneven usage of nucleotides. The overall count of A nucleotide was found to be the highest while the overall count of G nucleotide was the least. The CO genes across the three reptilian orders were prominently AT biased. Comparison of the GC proportion at each codon position displayed that GC1 percentage ranked the highest in all the three CO genes of the reptilian orders. SCUO values indicated weaker CUB, while considerable variation of SCUO values existed in the three CO genes across the studied reptiles. Relative synonymous codon usage (RSCU) values indicated that mostly the A ending codons were preferred. Based on the parameters namely neutrality plot, mutational responsive index and translational selection, we could conclude that natural selection was the major evolutionary force in COI, COII and COIII genes in the studied reptilian orders. However, correspondence analysis, parity plot and correlation studies indicated the existence of mutation pressure as well on the CO genes.

Poorly Conserved P15 Proteins of Cileviruses Retain Elements of Common Ancestry and Putative Functionality: A Theoretical Assessment on the Evolution of Cilevirus Genomes

The genus Cilevirus groups enveloped single-stranded (+) RNA virus members of the family Kitaviridae, order Martellivirales. Proteins P15, scarcely conserved polypeptides encoded by cileviruses, have no apparent homologs in public databases. Accordingly, the open reading frames (ORFs) p15, located at the 5'-end of the viral RNA2 molecules, are considered orphan genes (ORFans). In this study, we have delved into ORFs p15 and the relatively poorly understood biochemical properties of the proteins P15 to posit their importance for viruses across the genus and theorize on their origin. We detected that the ORFs p15 are under purifying selection and that, in some viral strains, the use of synonymous codons is biased, which might be a sign of adaptation to their plant hosts. Despite the high amino acid sequence divergence, proteins P15 show the conserved motif [FY]-L-x(3)-[FL]-H-x-x-[LIV]-S-C-x-C-x(2)-C-x-G-x-C, which occurs exclusively in members of this protein family. Proteins P15 also show a common predicted 3D structure that resembles the helical scaffold of the protein ORF49 encoded by radinoviruses and the phosphoprotein C-terminal domain of mononegavirids. Based on the 3D structural similarities of P15, we suggest elements of common ancestry, conserved functionality, and relevant amino acid residues. We conclude by postulating a plausible evolutionary trajectory of ORFans p15 and the 5'-end of the RNA2 of cileviruses considering both protein fold superpositions and comparative genomic analyses with the closest kitaviruses, negeviruses, nege/kita-like viruses, and unrelated viruses that share the ecological niches of cileviruses.

Codon Optimization of Insect Odorant Receptor Genes May Increase Their Stable Expression for Functional Characterization in HEK293 Cells

Insect odorant receptor (OR) genes are routinely expressed in Human Embryonic Kidney (HEK) 293 cells for functional characterization ("de-orphanization") using transient or stable expression. However, progress in this research field has been hampered because some insect ORs are not functional in this system, which may be due to insufficient protein levels. We investigated whether codon optimization of insect OR sequences for expression in human cells could facilitate their functional characterization in HEK293 cells with stable and inducible expression. We tested the olfactory receptor co-receptor (Orco) proteins from the bark beetles Ips typographus ("Ityp") and Dendroctonus ponderosae ("Dpon"), and six ItypORs previously characterized in Xenopus laevis oocytes and/or HEK cells. Western blot analysis indicated that codon optimization yielded increased cellular protein levels for seven of the eight receptors. Our experimental assays demonstrated that codon optimization enabled functional characterization of two ORs (ItypOR25 and ItypOR29) which are unresponsive when expressed from wildtype (non-codon optimized) genes. Similar to previous Xenopus oocyte recordings, ItypOR25 responded primarily to the host/conifer monoterpene (+)-3-carene. ItypOR29 responded primarily to (+)-isopinochamphone and similar ketones produced by fungal symbionts and trees. Codon optimization also resulted in significantly increased responses in ItypOR49 to its pheromone ligand (R)-(-)-ipsdienol, and improved responses to the Orco agonist VUAA1 in ItypOrco. However, codon optimization did not result in functional expression of DponOrco, ItypOR23, ItypOR27, and ItypOR28 despite higher protein levels as indicated by Western blots. We conclude that codon optimization may enable or improve the functional characterization of insect ORs in HEK cells, although this method is not sufficient for all ORs that are not functionally expressed from wildtype genes.

Edging on Mutational Bias, Induced Natural Selection From Host and Natural Reservoirs Predominates Codon Usage Evolution in Hantaan Virus

The molecular evolutionary dynamics that shape hantaviruses’ evolution are poorly understood even now, besides the contribution of virus-host interaction to their evolution remains an open question. Our study aimed to investigate these two aspects in Hantaan virus (HTNV)—the prototype of hantaviruses and an emerging zoonotic pathogen that infects humans, causing hemorrhagic fever with renal syndrome (HFRS): endemic in Far East Russia, China, and South Korea—via a comprehensive, phylogenetic-dependent codon usage analysis. We found that host- and natural reservoir-induced natural selection is the primary determinant of its biased codon choices, exceeding the mutational bias effect. The phylogenetic analysis of HTNV strains resulted in three distinct clades: South Korean, Russian, and Chinese. An effective number of codon (ENC) analysis showed a slightly biased codon usage in HTNV genomes. Nucleotide composition and RSCU analyses revealed a significant bias toward A/U nucleotides and A/U-ended codons, indicating the potential influence of mutational bias on the codon usage patterns of HTNV. Via ENC-plot, Parity Rule 2 (PR2), and neutrality plot analyses, we would conclude the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of HTNV; however, natural selection is the dominant factor influencing its codon usage bias. Codon adaptation index (CAI), Relative codon deoptimization index (RCDI), and Similarity Index (SiD) analyses uncovered the intense selection pressure from the host (Human) and natural reservoirs (Striped field mouse and Chinese white-bellied rat) in shaping HTNV biased codon choices. Our study clearly revealed the evolutionary processes in HTNV and the role of virus-host interaction in its evolution. Moreover, it opens the door for a more comprehensive codon usage analysis for all hantaviruses species to determine their molecular evolutionary dynamics and adaptability to several hosts and environments. We believe that our research will help in a better and deep understanding of HTNV evolution that will serve its future basic research and aid live attenuated vaccines design.

Genetic changes and evolutionary analysis of canine circovirus

Canine circovirus (canineCV) has been found to be associated with vasculitis, hemorrhage, hemorrhagic enteritis, and diarrhea of canines. CanineCV, like other circoviruses, may also be associated with lymphoid depletion and immunosuppression. This circovirus has been detected worldwide in different countries and species. Recombination and mutation events in the canineCV genome have been described, indicating that the virus is continuing to evolve. However, the origin, codon usage patterns, and host adaptation of canineCV remain to be studied. Here, the coding sequences of 93 canineCV sequences available in the GenBank database were used for analysis. The results showed that canineCV sequences could be classified into five genotypes, as confirmed by phylogenetic and principal component analysis (PCA). Maximum clade credibility (MCC) and maximum-likelihood (ML) trees suggested that canineCV originated from bat circovirus. G/T and A/C nucleotide biases were observed in ORF1 and ORF2, respectively, and a low codon usage bias (CUB) was found in canineCV using an effective number of codon (ENC) analysis. Correlation analysis, ENC plot analysis and neutrality plot analysis indicated that the codon usage pattern was mainly shaped by natural selection. Codon adaptation index (CAI) analysis, relative codon deoptimization index (RCDI) analysis, and similarity index (SiD) analysis revealed a better adaption to Vulpes vulpes than to Canis familiaris. Furthermore, a cross-species transmission hypothesis that canineCV may have evolved from bats (origin analysis) and subsequently adapted to wolves, arctic foxes, dogs, and red foxes, was proposed. This study contributes to our understanding of the factors related to canineCV evolution and host adaption.

Homopeptide and homocodon levels across fungi are coupled to GC/AT-bias and intrinsic disorder, with unique behaviours for some amino acids

Homopeptides (runs of one amino-acid type) are evolutionarily important since they are prone to expand/contract during DNA replication, recombination and repair. To gain insight into the genomic/proteomic traits driving their variation, we analyzed how homopeptides and homocodons (which are pure codon repeats) vary across 405 Dikarya, and probed their linkage to genome GC/AT bias and other factors. We find that amino-acid homopeptide frequencies vary diversely between clades, with the AT-rich Saccharomycotina trending distinctly. As organisms evolve, homocodon and homopeptide numbers are majorly coupled to GC/AT-bias, exhibiting a bi-furcated correlation with degree of AT- or GC-bias. Mid-GC/AT genomes tend to have markedly fewer simply because they are mid-GC/AT. Despite these trends, homopeptides tend to be GC-biased relative to other parts of coding sequences, even in AT-rich organisms, indicating they absorb AT bias less or are inherently more GC-rich. The most frequent and most variable homopeptide amino acids favour intrinsic disorder, and there are an opposing correlation and anti-correlation versus homopeptide levels for intrinsic disorder and structured-domain content respectively. Specific homopeptides show unique behaviours that we suggest are linked to inherent slippage probabilities during DNA replication and recombination, such as poly-glutamine, which is an evolutionarily very variable homopeptide with a codon repertoire unbiased for GC/AT, and poly-lysine whose homocodons are overwhelmingly made from the codon AAG.

Codon-based indices for modeling gene expression and transcript evolution

Codon usage bias (CUB) refers to the phenomena that synonymous codons are used in different frequencies in most genes and organisms. The general assumption is that codon biases reflect a balance between mutational biases and natural selection. Today we understand that the codon content is related and can affect all gene expression steps. Starting from the 1980s, codon-based indices have been used for answering different questions in all biomedical fields, including systems biology, agriculture, medicine, and biotechnology. In general, codon usage bias indices weigh each codon or a small set of codons to estimate the fitting of a certain coding sequence to a certain phenomenon (e.g., bias in codons, adaptation to the tRNA pool, frequencies of certain codons, transcription elongation speed, etc.) and are usually easy to implement. Today there are dozens of such indices; thus, this paper aims to review and compare the different codon usage bias indices, their applications, and advantages. In addition, we perform analysis that demonstrates that most indices tend to correlate even though they aim to capture different aspects. Due to the centrality of codon usage bias on different gene expression steps, it is important to keep developing new indices that can capture additional aspects that are not modeled with the current indices.

Co-evolution between codon usage and protein-protein interaction in bacteria

We study the correlation between the codon usage bias of genetic sequences and the network features of protein-protein interaction (PPI) in bacterial species. We use PCA techniques in the space of codon bias indices to show that genes with similar patterns of codon usage have a significantly higher probability that their encoded proteins are functionally connected and interacting. Importantly, this signal emerges when multiple aspects of codon bias are taken into account at the same time. The present study extends our previous observations on E. coli over a wide set of 34 bacteria. These findings could allow for future investigations on the possible effects of codon bias on the topology of the PPI network, with the aim of improving existing bioinformatics methods for predicting protein interactions.

Analysis of synonymous codon usage of transcriptome database in Rheum palmatum

Background

Rheum palmatum is an endangered and important medicinal plant in Asian countries, especially in China. However, there is little knowledge about the codon usage bias for R. palmatum CDSs. In this project, codon usage bias was determined based on the R. palmatum 2,626 predicted CDSs from R. palmatum transcriptome.

Methods

In this study, all codon usage bias parameters and nucleotide compositions were calculated by Python script, Codon W, DNA Star, CUSP of EMBOSS.

Results

The average GC and GC3 content are 46.57% and 46.6%, respectively, the results suggested that there exists a little more AT than GC in the R. palmatum genes, and the codon bias of R. palmatum genes preferred to end with A/T. We concluded that the codon bias in R. palmatum was affect by nucleotide composition, mutation pressure, natural selection, gene expression levels, and the mutation pressure is the prominent factor. In addition, we figured out 28 optimal codons and most of them ended with A or U. The project here can offer important information for further studies on enhancing the gene expression using codon optimization in heterogeneous expression system, predicting the genetic and evolutionary mechanisms in R. palmatum.

The nucleotide usages significantly impact synonymous codon usage in Mycoplasma hyorhinis

Five annotated genomes of Mycoplasma hyorhinis were analyzed for clarifying evolutionary dynamics driving the overall codon usage pattern. Information entropy used for estimating nucleotide usage pattern at the gene level indicates that multiple evolutionary dynamics participate in forcing nucleotide usage bias at every codon position. Moreover, nucleotide usage bias directly contributes to synonymous codon usage biases with two different extremes. The overrepresented codons tended to have A/T in the third codon position, and the underrepresented codons strongly used G/C in the third position. Furthermore, correspondence analysis and neutrality plot reflect an obvious interplay between mutation pressure and natural selection mediating codon usage in M. hyorhinis genome. Due to significant bias in usages between A/T and G/C at the gene level, different selective forces have been proposed to contribute to codon usage preference in M. hyorhinis genome, including nucleotide composition constraint derived from mutation pressure, translational selection involved in natural selection, and strand-specific mutational bias represented by different nucleotide skew index. The systemic analyses of codon usage for M. hyorhinis can enable us to better understand the mechanisms of evolution in this species.

A versatile toolkit for CRISPR-Cas13-based RNA manipulation in Drosophila

Advances in CRISPR technology have immensely improved our ability to manipulate nucleic acids, and the recent discovery of the RNA-targeting endonuclease Cas13 adds even further functionality. Here, we show that Cas13 works efficiently in Drosophila, both ex vivo and in vivo. We test 44 different Cas13 variants to identify enzymes with the best overall performance and show that Cas13 could target endogenous Drosophila transcripts in vivo with high efficiency and specificity. We also develop Cas13 applications to edit mRNAs and target mitochondrial transcripts. Our vector collection represents a versatile tool collection to manipulate gene expression at the post-transcriptional level.

Conservation analysis of SARS-CoV-2 spike suggests complicated viral adaptation history from bat to human

Background

The current coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome (SARS)-CoV-2, has become the most devastating public health emergency in the 21st century and one of the most influential plagues in history. Studies on the origin of SARS-CoV-2 have generally agreed that the virus probably comes from bat, closely related to a bat CoV named BCoV-RaTG13 taken from horseshoe bat (Rhinolophus affinis), with Malayan pangolin (Manis javanica) being a plausible intermediate host. However, due to the relatively low number of SARS-CoV-2-related strains available in public domain, the evolutionary history remains unclear.

Methodology

Nine hundred ninety-five coronavirus sequences from NCBI Genbank and GISAID were obtained and multiple sequence alignment was carried out to categorize SARS-CoV-2 related groups. Spike sequences were analyzed using similarity analysis and conservation analyses. Mutation analysis was used to identify variations within receptor-binding domain (RBD) in spike for SARS-CoV-2-related strains.

Results

We identified a family of SARS-CoV-2-related strains, including the closest relatives, bat CoV RaTG13 and pangolin CoV strains. Sequence similarity analysis and conservation analysis on spike sequence identified that N-terminal domain, RBD and S2 subunit display different degrees of conservation with several coronavirus strains. Mutation analysis on contact sites in SARS-CoV-2 RBD reveals that human-susceptibility probably emerges in pangolin.

Conclusion and implication

We conclude that the spike sequence of SARS-CoV-2 is the result of multiple recombination events during its transmission from bat to human, and we propose a framework of evolutionary history that resolve the relationship of BCoV-RaTG13 and pangolin coronaviruses with SARS-CoV-2.

Lay Summary

This study analyses whole-genome and spike sequences of coronavirus from NCBI using phylogenetic and conservation analyses to reconstruct the evolutionary history of severe acute respiratory syndrome (SARS)-CoV-2 and proposes an evolutionary history of spike in the progenitors of SARS-CoV-2 from bat to human through mammal hosts before they recombine into the current form.

Large-scale in silico mutagenesis experiments reveal optimization of genetic code and codon usage for protein mutational robustness

Background

How, and the extent to which, evolution acts on DNA and protein sequences to ensure mutational robustness and evolvability is a long-standing open question in the field of molecular evolution. We addressed this issue through the first structurome-scale computational investigation, in which we estimated the change in folding free energy upon all possible single-site mutations introduced in more than 20,000 protein structures, as well as through available experimental stability and fitness data.

Results

At the amino acid level, we found the protein surface to be more robust against random mutations than the core, this difference being stronger for small proteins. The destabilizing and neutral mutations are more numerous in the core and on the surface, respectively, whereas the stabilizing mutations are about 4% in both regions. At the genetic code level, we observed smallest destabilization for mutations that are due to substitutions of base III in the codon, followed by base I, bases I+III, base II, and other multiple base substitutions. This ranking highly anticorrelates with the codon-anticodon mispairing frequency in the translation process. This suggests that the standard genetic code is optimized to limit the impact of random mutations, but even more so to limit translation errors. At the codon level, both the codon usage and the usage bias appear to optimize mutational robustness and translation accuracy, especially for surface residues.

Conclusion

Our results highlight the non-universality of mutational robustness and its multiscale dependence on protein features, the structure of the genetic code, and the codon usage. Our analyses and approach are strongly supported by available experimental mutagenesis data.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Detection of Recombinant Rousettus Bat Coronavirus GCCDC1 in Lesser Dawn Bats (Eonycteris spelaea) in Singapore

Rousettus bat coronavirus GCCDC1 (RoBat-CoV GCCDC1) is a cross-family recombinant coronavirus that has previously only been reported in wild-caught bats in Yúnnan, China. We report the persistence of a related strain in a captive colony of lesser dawn bats captured in Singapore. Genomic evidence of the virus was detected using targeted enrichment sequencing, and further investigated using deeper, unbiased high throughput sequencing. RoBat-CoV GCCDC1 Singapore shared 96.52% similarity with RoBat-CoV GCCDC1 356 (NC_030886) at the nucleotide level, and had a high prevalence in the captive bat colony. It was detected at five out of six sampling time points across the course of 18 months. A partial segment 1 from an ancestral Pteropine orthoreovirus, p10, makes up the recombinant portion of the virus, which shares high similarity with previously reported RoBat-CoV GCCDC1 strains that were detected in Yúnnan, China. RoBat-CoV GCCDC1 is an intriguing, cross-family recombinant virus, with a geographical range that expands farther than was previously known. The discovery of RoBat-CoV GCCDC1 in Singapore indicates that this recombinant coronavirus exists in a broad geographical range, and can persist in bat colonies long-term.

Comparative Analysis of the Complete Plastid Genome of Five Bupleurum Species and New Insights into DNA Barcoding and Phylogenetic Relationship

Bupleurum L. (Apiaceae) is a perennial and herbal genus, most species of which have high medicinal value. However, few studies have been performed using plastome data in this genus, and the phylogenetic relationships have always been controversial. In this study, the plastid genomes of Bupleurum chinense and Bupleurum commelynoideum were sequenced, and their gene content, order, and structure were counted and analyzed. The only three published Bupleurum species (B. boissieuanum, B. falcatum, and B. latissimum) and other fifteen allied species were selected to conduct a series of comparative and phylogenetic analyses. The genomes of B. chinense and B. commelynoideum were 155,869 and 155,629 bp in length, respectively, both of which had a typical quadripartite structure. The genome length, structure, guanine and cytosine (GC) content, and gene distribution were highly similar to the other three Bupleurum species. The five Bupleurum species had nearly the same codon usages, and eight regions (petN-psbM, rbcL-accD, ccsA-ndhD, trnK(UUU)-rps16, rpl32-trnL(UAG)-ccsA, petA-psbJ, ndhF-rpl32, and trnP(UGG)-psaJ-rpl33) were found to possess relatively higher nucleotide diversity, which may be the promising DNA barcodes in Bupleurum. Phylogenetic analysis revealed that all Bupleurum species clustered into a monophyletic clade with high bootstrap support and diverged after the Chamaesium clade. Overall, our study provides new insights into DNA barcoding and phylogenetic relationship between Bupleurum and its related genera, and will facilitate the population genomics, conservation genetics, and phylogenetics of Bupleurum in Apiaceae.

Genetic Adaptations, Biases, and Evolutionary Analysis of Canine Distemper Virus Asia-4 Lineage in a Fatal Outbreak of Wild-Caught Civets in Thailand

Canine morbillivirus (CDV) is a serious pathogen that can cause fatal systemic disease in a wide range of domestic and wildlife carnivores. Outbreaks of CDV in wildlife species lead to questions regarding the dispersal of the CDV origin. In the present study, we identified a fatal CDV outbreak in caged wild-caught civets in Thailand. Full-length genetic analysis revealed that CDV from the Asia-4 lineage served as the likely causative agent, which was supported by the viral localization in tissues. Evolutionary analysis based on the CDV hemagglutinin (H) gene revealed that the present civet CDV has co-evolved with CDV strains in dogs in Thailand since about 2014. The codon usage pattern of the CDV H gene revealed that the CDV genome has a selective bias of an A/U-ended codon preference. Furthermore, the codon usage pattern of the CDV Asia-4 strain from potential hosts revealed that the usage pattern was related more to the codon usage of civets than of dogs. This finding may indicate the possibility that the discovered CDV had initially adapted its virulence to infect civets. Therefore, the CDV Asia-4 strain might pose a potential risk to civets. Further epidemiological, evolutionary, and codon usage pattern analyses of other CDV-susceptible hosts are required.