Codon Usage Preference Research Articles

Genetic Insights and Phylogenetic Relationships of Lycorma delicatula (Hemiptera: Fulgoridae, 1845) Revealed Through Mitochondrial Genome Sequence Analysis

Lycorma delicatula, a phytophagous agricultural pest, is invading various regions across the globe. We employed the mitochondrial genome and mitochondrial marker NADH dehydrogenase subunit 2 (ND2) to elucidate the mitochondrial genome characterization and phylogenetic relationships. The mitochondrial genome revealed a circular mitochondrial structure with a length of 16 138 bp encompassing 37 genes and a GC% of 23.41%. Codon usage preferences exhibited significant variation, with A or U being the most preferred bases in the third codon position. Most tRNAs displayed a typical cloverleaf structure, while a few exhibited base mismatches, primarily in G-U pairs. Phylogenetic analysis using the ND2 genes indicated that the population in northwest China could be tightly near together, with a genetic distance ranging from 0 to 0.01. Speculations arising from the mitochondrial genome analysis suggest that species from northern China were closer relatives to invaders from other regions, Lycorma meliae emerged as a close relative of Lycorma delicatula among the Hemipterans. Collectively, our findings enhance the understanding of phylogenetic relationships and complement the mitochondrial genomic information of Lycorma delicatula, thereby contributing valuable data for subsequent studies on this pest species.

Analysis of Preferred Codon Usage in the ZIKA Virus Genome and Their Implications for Genome Evolution and Vaccine Design

Analysis of Preferred Codon Usage in the ZIKA Virus Genome and Their Implications for Genome Evolution and Vaccine Design

The complete mitochondrial genome of Castanopsis carlesii and Castanea henryi reveals the rearrangement and size differences of mitochondrial DNA molecules

BackgroundCastanopsis carlesii is a dominant tree species in subtropical evergreen broad-leaved forests and holds significant ecological value. It serves as an excellent timber tree species and raw material for cultivating edible fungi. Henry Chinquapin (Castanea henryi) wood is known for its hardness and resistance to water and moisture, making it an exceptional timber species. Additionally, its fruit has a sweet and fruity taste, making it a valuable food source. However, the mitogenomes of these species have not been previously reported. To gain a better understanding of them, this study successfully assembled high-quality mitogenomes of C. carlesii and Ca. henryi for the first time.ResultsOur research reveals that the mitochondrial DNA (mtDNA) of C. carlesii exhibits a unique multi-branched conformation, while Ca. henryi primarily exists in the form of two independent molecules that can be further divided into three independent molecules through one pair of long repetitive sequences. The size of the mitogenomes of C. carlesii and Ca. henryi are 592,702 bp and 379,929 bp respectively, which are currently the largest and smallest Fagaceae mitogenomes recorded thus far. The primary factor influencing mitogenome size is dispersed repeats. Comparison with published mitogenomes from closely related species highlights differences in size, gene loss patterns, codon usage preferences, repetitive sequences, as well as mitochondrial plastid DNA segments (MTPTs).ConclusionsOur study enhances the understanding of mitogenome structure and evolution in Fagaceae, laying a crucial foundation for future research on cell respiration, disease resistance, and other traits in this family.

Expression analysis of the apple HSP70 gene family in abiotic stress and phytohormones and expression validation of candidate MdHSP70 genes

Heat shock protein 70 (HSP70) is one kind of molecular chaperones which are widely found in organisms, and its members are highly conserved among each other, with important roles in plant growth and development. In this study, 56 HSP70 genes were identified from the apple genome database. Analysis of gene duplication events showed that tandem and segmental duplication events play an important role in promoting the amplification of the MdHSP70 gene family. Collinearity analysis showed that HSP70 family members of apple were more closely related to HSP70 family members of Arabidopsis, tomato and soybean. The promoter region of the apple HSP70 genes contains a large number of cis-acting elements in response to hormones and stress. Tissue-specific expression analysis showed that some of the genes were associated with various stages of the apple growth process. Codon preference analysis showed small differences between codon bases 1 and 3 in the apple HSP70 genome, and the codon base composition had a small effect on codon usage preference. The multiple expression patterns of the MdHSP70 gene suggested that MdHSP70 gene members play important roles in growth and development and in response to hormonal and abiotic stresses. The yeast two-hybrid (Y2H) demonstrated that MdHSP70-53 interacts with MdDVH24_032563. The qRT-PCR analysis showed that most MdHSP70 members’ hormonal and abiotic stresses (MdHSP70-6, MdHSP70-26 and MdHSP70-45) appeared to be highly expressed. To further elucidate the function of MdHSP70 (6, 26, 45), we introduced them into tobacco to confirm subcellular locations and noted that these genes are located in the cytoplasm and cell membrane. This study serves as a theoretical basis for further studies of the MdHSP70 gene and helps to further investigate the functional characterization of MdHSP70 gene.

GenRCA: a user-friendly rare codon analysis tool for comprehensive evaluation of codon usage preferences based on coding sequences in genomes

BackgroundThe study of codon usage bias is important for understanding gene expression, evolution and gene design, providing critical insights into the molecular processes that govern the function and regulation of genes. Codon Usage Bias (CUB) indices are valuable metrics for understanding codon usage patterns across different organisms without extensive experiments. Considering that there is no one-fits-all index for all species, a comprehensive platform supporting the calculation and analysis of multiple CUB indices for codon optimization is greatly needed.ResultsHere, we release GenRCA, an updated version of our previous Rare Codon Analysis Tool, as a free and user-friendly website for all-inclusive evaluation of codon usage preferences of coding sequences. In this study, we manually reviewed and implemented up to 31 codon preference indices, with 65 expression host organisms covered and batch processing of multiple gene sequences supported, aiming to improve the user experience and provide more comprehensive and efficient analysis.ConclusionsOur website fills a gap in the availability of comprehensive tools for species-specific CUB calculations, enabling researchers to thoroughly assess the protein expression level based on a comprehensive list of 31 indices and further guide the codon optimization.

Comparative analysis of codon usage bias in the chloroplast genomes of eighteen Ampelopsideae species (Vitaceae)

BackgroundThe tribe Ampelopsideae plants are important garden plants with both medicinal and ornamental values. The study of codon usage bias (CUB) facilitates a deeper comprehension of the molecular genetic evolution of species and their adaptive strategies. The joint analysis of CUB in chloroplast genomes (cpDNA) offers valuable insights for in-depth research on molecular genetic evolution, biological resource conservation, and elite breeding within this plant family.ResultsThe base composition and codon usage preferences of the eighteen chloroplast genomes were highly similar, with the GC content of bases at all positions of their codons being less than 50%. This indicates that they preferred A/T bases. Their effective codon numbers were all in the range of 35–61, which indicates that the codon preferences of the chloroplast genomes of the 18 Ampelopsideae plants were relatively weak. A series of analyses indicated that the codon preference of the chloroplast genomes of the 18 Ampelopsideae plants was influenced by a combination of multiple factors, with natural selection being the primary influence. The clustering tree generated based on the relative usage of synonymous codons is consistent with some of the results obtained from the phylogenetic tree of chloroplast genomes, which indicates that the clustering tree based on the relative usage of synonymous codons can be an important supplement to the results of the sequence-based phylogenetic analysis. Eventually, 10 shared best codons were screened on the basis of the chloroplast genomes of 18 species.ConclusionThe codon preferences of the chloroplast genome in Ampelopsideae plants are relatively weak and are primarily influenced by natural selection. The codon composition of the chloroplast genomes of the eighteen Ampelopsideae plants and their usage preferences were sufficiently similar to demonstrate that the chloroplast genomes of Ampelopsideae plants are highly conserved. This study provides a scientific basis for the genetic evolution of chloroplast genes in Ampelopsideae species and their suitable strategies.

Genomic factors shaping codon usage across the Saccharomycotina subphylum.

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. It has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. We characterized codon usage across 1,154 strains from 1,051 species from the fungal subphylum Saccharomycotina to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns. We found a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is distinct between the 12 orders to such a degree that yeasts can be classified with an accuracy >90% using a machine learning algorithm. We also characterized the degree to which codon usage bias is impacted by translational selection. We found it was influenced by a combination of features, including the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that avoidance of the CGN codons is associated with a decline in arginine tRNA function. Consistent with previous findings, codon usage bias within the Saccharomycotina is shaped by genomic features and GC bias. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Codon Usage Analysis Reveals Distinct Evolutionary Patterns and Host Adaptation Strategies in Duck Hepatitis Virus 1 (DHV-1) Phylogroups.

Duck hepatitis virus 1 (DHV-1) is a major threat to the global poultry industry, causing significant economic losses due to high mortality rates in young ducklings. To better understand the evolution and host adaptation strategies of DHV-1, we conducted a comprehensive codon usage analysis of DHV-1 genomes. Our phylogenetic analysis revealed three well-supported DHV-1 phylogroups (Ia, Ib, and II) with distinct genetic diversity patterns. Comparative analyses of the codon usage bias and dinucleotide abundance uncovered a strong preference for A/U-ended codons and a biased pattern of dinucleotide usage in the DHV-1 genome, with CG dinucleotides being extremely underrepresented. Effective number of codons (ENC) analysis indicated a low codon usage bias in the DHV-1 ORF sequences, suggesting adaptation to host codon usage preferences. PR2 bias, ENC plot, and neutrality analyses revealed that both mutation pressure and natural selection influence the codon usage patterns of DHV-1. Notably, the three DHV-1 phylogroups exhibited distinct evolutionary trends, with phylogroups Ia and Ib showing evidence of neutral evolution accompanied by selective pressure, while the phylogroup II evolution was primarily driven by random genetic drift. Comparative analysis of the codon usage indices (CAI, RCDI, and SiD) among the phylogroups highlighted significant differences between subgroups Ia and Ib, suggesting distinct evolutionary pressures or adaptations influencing their codon usage. These findings contribute to our understanding of DHV-1 evolution and host adaptation, with potential implications for the development of effective control measures and vaccines.

Comparative genome wise analysis of codon usage of Staphylococcus Genus.

The genus Staphylococcus encompasses a diverse array of bacteria with significant implications for human health, including disreputable pathogens such as Staphylococcus aureus and Staphylococcus epidermidis. Understanding the genetic composition and codon usage patterns of Staphylococcus species is crucial for unraveling their evolutionary dynamics, adaptive strategies, and pathogenic potential. In this study, we conducted a comprehensive analysis of codon usage patterns across 48 species within the Staphylococcus genus. Our findings uncovered variations in genomic G-C content across Staphylococcus species, impacting codon usage preferences, with a notable preference for A/T-rich codons observed in pathogenic strains. This preference for A/T-rich codons suggests an energy-saving strategy in pathogenic organisms. Analysis of dinucleotide pair expression patterns unveiled insights into genomic dynamics, with overrepresented codon pairs reflecting trends in dinucleotide expression across genomes. Additionally, a significant correlation between CAI and genomic G-C content underscored the intricate relationship between codon usage patterns and gene expression strategies. Amino acid usage analysis highlighted preferences for energetically cheaper amino acids, suggesting adaptive strategies promoting energy efficiency. This comprehensive analysis sheds light on the evolutionary dynamics and adaptive mechanisms employed by Staphylococcus species, providing valuable insights into their pathogenic potential and clinical implications. Understanding these genomic features is crucial for devising strategies to combat staphylococcal infections and improve public health outcomes.

Integration of Illumina and PacBio HiFi Sequencing Reveals a Three-Linear-Molecule Mitogenome with RNA-Editing Sites and Phylogeny in Arrow Bamboo (Fargesia qinlingensis)

Arrow bamboo (Fargesia qinlingensis) is endemic to the Qinling Mountains and has remarkable adaptive resilience to changing climates. However, its complete mitogenome remains unknown. Using the Illumina and PacBio HiFi sequencing platforms, we found that the mitogenome assembly of the F. qinlingensis has a multi-branched skeleton comprising three linear molecules (M1, M2, and M3), with a length of 442,368 bp and a GC content of 44.05%. Thirty-five unique PCGs were identified in the complete mitogenome, including twenty-four core structural genes, eleven noncore structural genes, three rRNAs, and sixteen tRNAs. The GCU for alanine and CAA for glutamine represented the most significant frequency (RSCU = 1.55) in the codon usage preference. A total of 51, 28, and 14 SSRs were determined on M1, M2, and M3, respectively. The mitogenome contained 149 pairs of dispersed repeats with lengths greater than 30 bp, the most abundant of which were 82 forward and 67 palindromic repeats. A long repeat sequence (14,342 bp) was characterized in mediating mitogenome recombination. DNA transfer analyses suggested that 44 MTPTs (30,943 bp, 6.99%) originated from the plastome. Among the 482 potential C-U/T RNA-editing sites predicted in 35 PCGs, ccmFn (38 times) and ccmC (36 times) shoed the highest frequency. Collinearity and phylogenetic trees revealed the close relationship between F. qinlingensis and Bambusa oldhamii. The primary features of the mitogenome of F. qinlingensis will help decipher the functional mitochondrial traits related to growth performance and climate resilience. Moreover, our findings provide insights into the evolution, environmental adaptation, and sustainable use of subalpine bamboo resources in the Qinling Mountains.

Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae.

The mitochondrial genomes of three species of Yangiella were sequenced, annotated, and analyzed. The genome length of the three species of the genus is 15,070-15,202 bp, with a typical gene number, including a control region, 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs). It was found that the mitochondrial genome of Yangiella had AT bias. Except for the lack of a DHU arm of the trnS1 gene, the other tRNAs had a typical cloverleaf structure, and the codon usage preferences of the three species exhibited high similarity. In addition, tRNA gene rearrangements were observed among the three subfamilies of Aradidae (Mezirinae, Calisiinae, Aradinae), and it was found that codon usage preferences appeared to be less affected by base mutation and more by natural selection. The Pi and Ka/Ks values indicated that cox1 was the most conserved gene in the mitochondrial genome of Aradidae, while atp8 and nad6 were rapidly evolved genes. Substitution saturation level analysis showed that the nucleic acid sequence of mitochondrial protein-coding genes in Aradidae did not reach saturation, suggesting the rationality of the phylogenetic analysis data. Bayesian and maximum likelihood methods were used to analyze the phylogeny of 16 species of Hemiptera insects, which supported the monophyly of Aneurinae, Carventinae, and Mezirinae, as well as the monophyly of Yangiella. Based on fossils and previous studies, the differentiation time was inferred, indicating that Yangiella diverged about 57 million years ago.

Patterns in Genome-Wide Codon Usage Bias in Representative Species of Lycophytes and Ferns.

The latest research shows that ferns and lycophytes have distinct evolutionary lineages. The codon usage patterns of lycophytes and ferns have not yet been documented. To investigate the gene expression profiles across various plant lineages with respect to codon usage, analyze the disparities and determinants of gene evolution in primitive plant species, and identify appropriate exogenous gene expression platforms, the whole-genome sequences of four distinct species were retrieved from the NCBI database. The findings indicated that Ceratopteris richardii, Adiantum capillus-veneris, and Selaginella moellendorffii exhibited an elevated A/U content in their codon base composition and a tendency to end with A/U. Additionally, S. capillus-veneris had more C/G in its codons and a tendency to end with C/G. The ENC values derived from both ENC-plot and ENC-ratio analyses deviated significantly from the standard curves, suggesting that the codon usage preferences of these four species were primarily influenced by genetic mutations and natural selection, with natural selection exerting a more prominent influence. This finding was further supported by PR2-Plot, neutrality plot analysis, and COA. A combination of RSCU and ENC values was used as a reference criterion to rank the codons and further identify the optimal codons. The study identified 24 high-frequency codons in C. richardii, A. capillus-veneris, and Diphasiastrum complanatum, with no shared optimal codons among the four species. Arabidopsis thaliana and Ginkgo biloba exhibited similar codon preferences to the three species, except for S. moellendorffii. This research offers a theoretical framework at the genomic codon level for investigating the phylogenetic relationships between lycophytes and ferns, shedding light on gene codon optimization and its implications for genetic engineering in breeding.

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

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

Mitogenome Assembly Reveals Gene Migration and RNA Editing Events in Plateau Hongliu (Myricaria elegans Royle.)

The Plateau Hongliu (Myricaria elegans Royle.) is a woody shrub halophyte that thrives in arid areas of western Tibet, in the Himalayan Mountains. It is acclaimed as superior in saline stress acclimation and as a critical pharmaceutical resource of the Tibetan traditional herb. Nevertheless, the mitogenome in the genus Myricaria remains unknown. Here, using the Illumina and PacBio sequencing assays, the first complete mitogenome of the M. elegans revealed a multi-branched skeleton with a total length of 416,354 bp and GC content of 44.33%, comprising two circular molecules (M1 and 2). The complete mitogenome annotates 31 unique protein-encoding genes (PEGs), fifteen tRNAs, and three rRNA genes. The UAA exhibits the most prominent codon usage preference as a termination, followed by UUA codons for leucine. The mitogenome contains 99 simple sequence repeats and 353 pairs of dispersed repeats, displaying the most frequent in palindromic repeats. Gene transfer analyses identified 8438 bp of 18 homologous fragments from the plastome, accounting for 2.03% of the total length. Using the PREP suite, 350 C-U RNA editing sites were predicted, of which nad4 and ccmB were on the top frequency. Syntenic and phylogenetic analyses suggested weakly conserved patterns of M. elegans in Caryophyllales owing to the genome rearrangement. In summary, the deciphered unique features and complexities of the mitogenome in M. elegans provide novel insights into understanding the evolution and biological conservation underlying climate resilience in halophytes.

Genetic background of adaptation of Crimean-Congo haemorrhagic fever virus to the different tick hosts.

Crimean-Congo haemorrhagic fever orthonairovirus (CCHFV) is a negative-sense, single-stranded RNA virus with a segmented genome and the causative agent of a severe Crimean-Congo haemorrhagic fever (CCHF) disease. The virus is transmitted mainly by tick species in Hyalomma genus but other ticks such as representatives of genera Dermacentor and Rhipicephalus may also be involved in virus life cycle. To improve our understanding of CCHFV adaptation to its tick species, we compared nucleotide composition and codon usage patterns among the all CCHFV strains i) which sequences and other metadata as locality of collection and date of isolation are available in GenBank and ii) which were isolated from in-field collected tick species. These criteria fulfilled 70 sequences (24 coding for S, 23 for M, and 23 for L segment) of virus isolates originating from different representatives of Hyalomma and Rhipicephalus genera. Phylogenetic analyses confirmed that Hyalomma- and Rhipicephalus-originating CCHFV isolates belong to phylogenetically distinct CCHFV clades. Analyses of nucleotide composition among the Hyalomma- and Rhipicephalus-originating CCHFV isolates also showed significant differences, mainly in nucleotides located at the 3rd codon positions indicating changes in codon usage among these lineages. Analyses of codon adaptation index (CAI), effective number of codons (ENC), and other codon usage statistics revealed significant differences between Hyalomma- and Rhipicephalus-isolated CCHFV strains. Despite both sets of strains displayed a higher adaptation to use codons that are preferred by Hyalomma ticks than Rhipicephalus ticks, there were distinct codon usage preferences observed between the two tick species. These findings suggest that over the course of its long co-evolution with tick vectors, CCHFV has optimized its codon usage to efficiently utilize translational resources of Hyalomma species.

Comprehensive Genomics Investigation of Neboviruses Reveals Distinct Codon Usage Patterns and Host Specificity.

Neboviruses (NeVs) from the Caliciviridae family have been linked to enteric diseases in bovines and have been detected worldwide. As viruses rely entirely on the cellular machinery of the host for replication, their ability to thrive in a specific host is greatly impacted by the specific codon usage preferences. Here, we systematically analyzed the codon usage bias in NeVs to explore the genetic and evolutionary patterns. Relative Synonymous Codon Usage and Effective Number of Codon analyses indicated a marginally lower codon usage bias in NeVs, predominantly influenced by the nucleotide compositional constraints. Nonetheless, NeVs showed a higher codon usage bias for codons containing G/C at the third codon position. The neutrality plot analysis revealed natural selection as the primary factor that shaped the codon usage bias in both the VP1 (82%) and VP2 (57%) genes of NeVs. Furthermore, the NeVs showed a highly comparable codon usage pattern to bovines, as reflected through Codon Adaptation Index and Relative Codon Deoptimization Index analyses. Notably, yak NeVs showed considerably different nucleotide compositional constraints and mutational pressure compared to bovine NeVs, which appear to be predominantly host-driven. This study sheds light on the genetic mechanism driving NeVs' adaptability, evolution, and fitness to their host species.

Codon usage bias and phylogenetic analysis of chloroplast genome in 36 gracilariaceae species.

Gracilariaceae is a group of marine large red algae and main source of agar with important economic and ecological value. The codon usage patterns of chloroplast genomes in 36 species from Graciliaceae show that GC range from 0.284 to 0.335, the average GC3 range from 0.135 to 0.243 and the value of ENC range from 35.098 to 42.327, which indicates these genomes are rich in AT and prefer to use codons ending with AT in these species. Nc plot, PR2 plot, neutrality plot analyses and correlation analysis indicate that these biases may be caused by multiple factors, such as natural selection and mutation pressure, but prolonged natural selection is the main driving force influencing codon usage preference. The cluster analysis and phylogenetic analysis show that the differentiation relationship of them is different and indicate that codons with weak or unbiased preferences may also play an irreplaceable role in these species' evolution. In addition, we identified 26 common high-frequency codons and 8-18 optimal codons all ending in A/U in these 36 species. Our results will not only contribute to carrying out transgenic work in Gracilariaceae species to maximize the protein yield in the future, but also lay a theoretical foundation for further exploring systematic classification of them.

Codon usage bias of human papillomavirus type 33 and 58: A comprehensive analysis.

Cervical cancer is closely linked to specific strains of human papillomavirus (HPV), notably HPV-33 and HPV-58, which exhibit a significant prevalence among women in China. Nevertheless, the codon usage bias in HPV-33 and HPV-58 is not well comprehended. The objective of this research is to analyze the codon usage patterns HPV-33 and HPV-58, pinpoint the primary factors that influence codon preference. The overall preference for codon usage in two HPV genotypes is not significant. Both HPV genotypes exhibit a preference for codons that end with A/U. The GC3 content for HPV-33 is 25.43% ± 0.35%, and for HPV-58, it is 29.44% ± 0.57%. Out of the 26 favored codons in HPV-33 and HPV-58 (relative synonymous codon usage (RSCU) > 1), 25 conclude with A/U. Principal component analysis (PCA) shows a tight clustering of the entire genome sequences of HPV-33 and HPV-58, suggesting a similarity in their RSCU preferences. Moreover, an examination of dinucleotide abundance indicated that translation selection influenced the development of a distinctive dinucleotide usage pattern in HPV-33 and HPV-58. Additionally, a combined analysis involving an effective number of codons plot, parity rule 2, and neutrality analysis demonstrated that, for HPV-33 and HPV-58, the primary determinant influencing codon usage preference is natural selection. HPV-33 and HPV-58 exhibit a restricted set of favored codons in common with humans, potentially mitigating competition for translation resources. Our discoveries could provide valuable perspectives on the evolutionary patterns and codon usage preferences of HPV-33 and HPV-58 viruses, contributing to the development and application of relevant HPV subtype vaccines.

Optimized transgene expression in the red alga Porphyridium purpureum and efficient recombinant protein secretion into the culture medium

Microalgae represent a promising but yet underexplored production platform for biotechnology. The vast majority of studies on recombinant protein expression in algae have been conducted in a single species, the green alga Chlamydomonas reinhardtii. However, due to epigenetic silencing, transgene expression in Chlamydomonas is often inefficient. Here we have investigated parameters that govern efficient transgene expression in the red microalga Porphyridium purpureum. Porphyridium is unique in that the introduced transformation vectors are episomally maintained as autonomously replicating plasmids in the nucleus. We show that full codon optimization to the preferred codon usage in the Porphyridium genome confers superior transgene expression, not only at the level of protein accumulation, but also at the level of mRNA accumulation, indicating that high translation rates increase mRNA stability. Our optimized expression constructs resulted in YFP accumulation to unprecedented levels of up to 5% of the total soluble protein. We also designed expression cassettes that target foreign proteins to the secretory pathway and lead to efficient protein secretion into the culture medium, thus simplifying recombinant protein harvest and purification. Our study paves the way to the exploration of red microalgae as expression hosts in molecular farming for recombinant proteins and metabolites.

Codon Optimization Enables the Geneticin Resistance Gene to Be Applied Efficiently to the Genetic Manipulation of the Plant Pathogenic Fungus Botrytis cinerea.

Botrytis cinerea can infect almost all of the important horticultural crops and cause severe economic losses globally every year. Modifying candidate genes and studying the phenotypic changes are among the most effective ways to unravel the pathogenic mechanism of this crop killer. However, few effective positive selection markers are used for B. cinerea genetic transformation, which limits multiple modifications to the genome, especially genes involving redundant functions. Here, we optimized a geneticin resistance gene, BcNPTII, based on the codon usage preference of B. cinerea. We found that BcNPTII can greatly increase the transformation efficiency of B. cinerea under G418 selection, with approximately 30 times higher efficiency than that of NPTII, which is applied efficiently to transform Magnaporthe oryzae. Using the gene replacement method, we successfully knocked out the second gene BOT2, with BcNPTII as the selection marker, from the mutant ΔoahA, in which OAHA was first replaced by the hygromycin resistance gene HPH in a field strain. We obtained the double knockout mutant ΔoahA Δbot2. Our data show that the codon-optimized BcNPTII is an efficient positive selection marker for B. cinerea transformation and can be used for various genetic manipulations in B. cinerea, including field wild-type strains.