Base composition and translational selection are insufficient to explain codon usage bias in plant viruses

Genetic analysis of tomato brown rugose fruit virus reveals evolutionary adaptation and codon usage bias patterns

Tomato brown rugose fruit virus (ToBRFV) poses a significant threat to tomato production worldwide, prompting extensive research into its genetic diversity, evolutionary dynamics, and adaptive strategies. In this study, we conducted a comprehensive analysis of ToBRFV at the codon level, focusing on codon usage bias, selection pressures, and evolutionary patterns across multiple genes. Our analysis revealed distinct patterns of codon usage bias and selection pressures within the ToBRFV genome, with varying levels of genetic diversity and evolutionary constraints among different genes. We observed a transition/transversion bias of 2.07 across the entire ToBRFV genome, with the movement protein (MP) gene exhibiting the highest transition/transversion bias and SNP density, suggesting potential evolutionary pressures or a higher mutation rate in this gene. Furthermore, our study identified episodic positive selection primarily in the MP gene, highlighting specific codons subject to adaptive changes in response to host immune pressures or environmental factors. Comparative analysis of codon usage bias in the coat protein (CP) and RNA-dependent RNA polymerase (RdRp) genes revealed gene-specific patterns reflecting functional constraints and adaptation to the host's translational machinery. Our findings provide valuable insights into the molecular mechanisms driving ToBRFV evolution and adaptation, with implications for understanding viral pathogenesis, host-virus interactions, and the development of control strategies. Future research directions include further elucidating the functional significance of codon usage biases, exploring the role of episodic positive selection in viral adaptation, and leveraging these insights to inform the development of effective antiviral strategies and crop protection measures.

Circular Single-Stranded DNA: Discovery, Biological Effects, and Applications

The field of nucleic acid therapeutics has witnessed a significant surge in recent times, as evidenced by the increasing number of approved genetic drugs. However, current platform technologies containing plasmids, lipid nanoparticle-mRNAs, and adeno-associated virus vectors encounter various limitations and challenges. Thus, we are devoted to finding a novel nucleic acid vector and have directed our efforts toward investigating circular single-stranded DNA (CssDNA), an ancient form of nucleic acid. CssDNAs are ubiquitous, but generally ignored. Accumulating evidence suggests that CssDNAs possess exceptional properties as nucleic acid vectors, exhibiting great potential for clinical applications in genetic disorders, gene editing, and immune cell therapy. Here, we comprehensively review the discovery and biological effects of CssDNAs as well as their applications in the field of biomedical research for the first time. Undoubtedly, as an ancient form of DNA, CssDNA holds immense potential and promises novel insights for biomedical research.

Compositional biases and evolution of the largest plant RNA virus order Patatavirales

Patatavirales is the largest order of plant RNA viruses and exclusively contains the family Potyviridae, accounting for 30 % of all known plant viruses. The composition bias of animal RNA viruses and several plant RNA viruses has been determined. However, the comprehensive nucleic acid composition, codon pair usage patterns, dinucleotide preference and codon pair preference of plant RNA viruses have not been investigated to date. In this study, integrated analysis and discussion of the nucleic acid composition, codon usage patterns, dinucleotide composition and codon pair bias of potyvirids were performed using 3732 complete genome coding sequences. The nucleic acid composition of potyvirids was significantly enriched in A/U. Interestingly, the A/U-rich nucleotide composition of Patatavirales is essential for determining the preferred A-ended and U-ended codons and the overexpression of UpG and CpA dinucleotides. The codon usage patterns and codon pair bias of potyvirids were significantly correlated with their nucleic acid composition. Additionally, the codon usage pattern, dinucleotide composition and codon-pair bias of potyvirids are more dependent on the classification of the virus compared with their hosts. Our analysis provides a better understanding of future research on the origin and evolution patterns of the order Patatavirales.

Evolution and host adaptability of plant RNA viruses: Research insights on compositional biases

During recent decades, many new emerging or re-emerging RNA viruses have been found in plants through the development of deep-sequencing technology and big data analysis. These findings largely changed our understanding of the origin, evolution and host range of plant RNA viruses. There is evidence that their genetic composition originates from viruses, and host populations play a key role in the evolution and host adaptability of plant RNA viruses. In this mini-review, we describe the state of our understanding of the evolution of plant RNA viruses in view of compositional biases and explore how they adapt to the host. It appears that adenine rich (A-rich) coding sequences, low CpG and UpA dinucleotide frequencies and lower codon usage patterns were found in the vast majority of plant RNA viruses. The codon usage pattern of plant RNA viruses was influenced by both natural selection and mutation pressure, and natural selection mostly from hosts was the dominant factor. The codon adaptation analyses support that plant RNA viruses probably evolved a dynamic balance between codon adaptation and deoptimization to maintain efficient replication cycles in multiple hosts with various codon usage patterns. In the future, additional combinations of computational and experimental analyses of the nucleotide composition and codon usage of plant RNA viruses should be addressed.

Synonymous Codon Usage Analysis of Three Narcissus Potyviruses

Narcissus degeneration virus (NDV), narcissus late season yellows virus (NLSYV) and narcissus yellow stripe virus (NYSV), which belong to the genus Potyvirus of the family Potyviridae, cause significant losses in the ornamental value and quality of narcissus. Several previous studies have explored the genetic diversity and evolution rate of narcissus viruses, but the analysis of the synonymous codons of the narcissus viruses is still unclear. Herein, the coat protein (CP) of three viruses is used to analyze the viruses’ phylogeny and codon usage pattern. Phylogenetic analysis showed that NYSV, NDV and NLSYV isolates were divided into five, three and five clusters, respectively, and these clusters seemed to reflect the geographic distribution. The effective number of codon (ENC) values indicated a weak codon usage bias in the CP coding region of the three narcissus viruses. ENC-plot and neutrality analysis showed that the codon usage bias of the three narcissus viruses is all mainly influenced by natural selection compared with the mutation pressure. The three narcissus viruses shared the same best optimal codon (CCA) and the synonymous codon prefers to use codons ending with A/U, compared to C/G. Our study shows the codon analysis of different viruses on the same host for the first time, which indicates the importance of the evolutionary-based design to control these viruses.

A Predicted Stem Loop in Coat Protein-Coding Sequence of Tobacco Vein Banding Mosaic Virus Is Required for Efficient Replication

Potyviral coat protein (CP) is involved in the replication and movement of potyviruses. However, little information is available on the roles of CP-coding sequence in potyviral infection. Here, we introduced synonymous substitutions to the codon C574G575C576 coding conserved residue arginine at position 192 (R192) of tobacco vein banding mosaic virus (TVBMV) CP. Substitution of the codon C574G575C576 to A574G575A576 or A574G575G576, but not C574G575A576, C574G575T576, or C574G575G576, reduced the replication, cell-to-cell movement, and accumulation of TVBMV in Nicotiana benthamiana plants, suggesting that C574 was critical for replication of TVBMV. Nucleotides 531 to 576 of the TVBMV CP-coding sequence were predicted to form a stem-loop structure, in which four consecutive C-G base pairs (C576-G531, C532-G575, C574-G533, and C534-G573) were located at the stem. Synonymous substitutions of R178-codon C532G533C534 to A532G533A534 and A532G533G534, but not C532G533A534, C532G533T534, or C532G533G534, reduced the replication levels, cell-to-cell, and systemic movement of TVBMV, suggesting that C532 was critical for TVBMV replication. Synonymous substitutions disrupting base pairs C576-G531 and C534-G573 did not affect viral accumulation. After three serial-passage inoculations, the accumulation of spontaneous mutant viruses was restored, and codons A532G533A534, A532G533G534, A574G575A576, or A574G575G576 of mutants were each separately changed to C532G533A534, C532G533G534, C574G575A576, or C574G575G576. Synonymous mutation of R178 and R192 also reduced viral accumulation in N. tabacum plants. Therefore, we concluded that the two consecutive C532-G575 and C574-G533 base pairs played critical roles in TVBMV replication via maintaining the stability of the stem-loop structures formed by nucleotides 531 to 576 of the CP-coding sequence.

Novel Fusari- and Toti-like Viruses, with Probable Different Origins, in the Plant Pathogenic Oomycete Globisporangiumultimum

To further classify the oomycete viruses that have been discovered in recent years, we investigated virus infection in the plant-parasitic oomycete Globisporangium ultimum in Japan. Double-stranded RNA detection, high-throughput sequencing, and RT-PCR revealed that the G. ultimum isolate UOP226 contained two viruses related to fusarivirus and totivirus, named Pythium ultimum RNA virus 1 (PuRV1) and Pythium ultimum RNA virus 2 (PuRV2), respectively. Phylogenetic analysis of the deduced amino acid sequence of the RNA-dependent RNA polymerase (RdRp) showed that fusari-like PuRV1 belonged to a different phylogenetic group than Plasmopara viticola lesion-associated fusari virus (PvlaFV) 1-3 from oomycete Plasmopara viticola. Codon usage bias of the PuRV1 RdRp gene was more similar to those of fungi than Globisporangium and Phytophthora, suggesting that the PuRV1 ancestor horizontally transmitted to G. ultimum ancestor from fungi. Phylogenetic analysis of the deduced amino acid sequence of the RdRp of toti-like PuRV2 showed a monophyletic group with the other toti-like oomycete viruses from Globisporangium, Phytophthora, and Pl. viticola. However, the nucleotide sequences of toti-like oomycete viruses were not so homologous, suggesting the possibility of convergent evolution of toti-like oomycete viruses.

Biased Mutation and Selection in RNA Viruses

RNA viruses are responsible for some of the worst pandemics known to mankind, including outbreaks of Influenza, Ebola, and COVID-19. One major challenge in tackling RNA viruses is the fact they are extremely genetically diverse. Nevertheless, they share common features that include their dependence on host cells for replication, and high mutation rates. We set out to search for shared evolutionary characteristics that may aid in gaining a broader understanding of RNA virus evolution, and constructed a phylogeny-based data set spanning thousands of sequences from diverse single-stranded RNA viruses of animals. Strikingly, we found that the vast majority of these viruses have a skewed nucleotide composition, manifested as adenine rich (A-rich) coding sequences. In order to test whether A-richness is driven by selection or by biased mutation processes, we harnessed the effects of incomplete purifying selection at the tips of virus phylogenies. Our results revealed consistent mutational biases toward U rather than A in genomes of all viruses. In +ssRNA viruses, we found that this bias is compensated by selection against U and selection for A, which leads to A-rich genomes. In -ssRNA viruses, the genomic mutational bias toward U on the negative strand manifests as A-rich coding sequences, on the positive strand. We investigated possible reasons for the advantage of A-rich sequences including weakened RNA secondary structures, codon usage bias, and selection for a particular amino acid composition, and conclude that host immune pressures may have led to similar biases in coding sequence composition across very divergent RNA viruses.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

Analysis of codon usage bias in potato virus Y non-recombinant strains

Potato virus Y (PVY) is a member of the genus Potyvirus, family Potyviridae, is considered one of the most devastating pest affecting economically important crops, such as potato, tobacco, tomato and pepper, representing a serious threat due to high incidence and worldwide distribution. Its economic significance as well as it biological and molecular complexities have aroused great attention, thus several studies have explore it genetic characteristics. However, little is known about PVY codon usage. To shed light on the relation of codon usage among viruses and their hosts is extremely important to understand virus survival, fitness and evolution. In this study, we performed a comprehensive analysis of codon usage and composition of PVY non-recombinant strains (PVY^N-NA, PVY^Eu-N, PVY^O, PVY^O5, PVY^C) based on 130 complete open reading frame sequences extracted from public databases. Furthermore, similarities between the synonymous codon usage of PVY and its main hosts were investigated. The results obtained in the current study suggest that the overall codon usage among PVY genotypes is similar and slightly biased. PVY codon usage is strongly influenced by mutational bias, but also by G + C compositional constraint and dinucleotide composition. Furthermore, similarities among codon usage preferences between PVY strains and analyzed hosts were observed.

Genome-wide codon usage pattern analysis reveals the correlation between codon usage bias and gene expression in Cuscuta australis

The protein-coding genes and pseudogenes of Cuscuta australis had the diverse contribution to the formation and evolution of parasitism. The codon usage pattern analysis of these two type genes could be used to understand the gene transcription and translation. In this study, we systematically analyzed the codon usage patterns of protein-coding sequences and pseudogenes sequences in C. australis. The results showed that the high frequency codons of protein coding sequences and pseudogenes had the same A/U bias in the third position. However, these two sequences had converse bias at the third base in optimal codons: the protein coding sequences preferred G/C-ending codons while pseudogene sequences preferred A/U-ending codons. Neutrality plot and effective number of codons plot revealed that natural selection played a more important role than mutation pressure in two sequences codon usage bias. Furthermore, the gene expression level had a significant positive correlation with codon usage bias in C. australis. Highly-expressed protein coding genes exhibited a higher codon bias than lowly-expressed genes. Meanwhile, the high-expression genes tended to use G/C-ending synonymous codons. This result further verified the optimal codons usage bias and its correlation with the gene expression in C. australis.

Molecular characteristics of segment 5, a unique fragment encoding two partially overlapping ORFs in the genome of rice black-streaked dwarf virus

Rice black-streaked dwarf virus (RBSDV), a ds-RNA virus in Fijivirus genus with family Reoviridae, which is transmitted by the small brown planthopper, is responsible for incidence of maize rough dwarf disease (MRDD) and rice black-streaked dwarf disease (RBSDD). To understand the variation and evolution of S5, a unique fragment in the genome of RBSDV which encodes two partially overlapping ORFs (ORF5-1 and ORF5-2), we analyzed 127 sequences from maize and rice exhibiting symptoms of dwarfism. The nucleotide diversity of both ORF5-1 (π = 0.039) and ORF5-2 (π = 0.027) was higher than that of the overlapping region (π = 0.011) (P < 0.05). ORF5-2 was under the greatest selection pressure based on codon bias analysis, and its activation was possibly influenced by the overlapping region. The recombinant fragments of three recombinant events (14NM23, 14BM20, and 14NM17) cross the overlapping region. Based on neighbor-joining tree analysis, the overlapping region could represent the evolutionary basis of the full-length S5, which was classified into three main groups. RBSDV populations were expanding and haplotype diversity resulted mainly from the overlapping region. The genetic differentiation of combinations (T127-B35, T127-J34, A58-B35, A58-J34, and B35-J34) reached significant or extremely significant levels. Gene flow was most frequent between subpopulations A58 and B35, with the smallest |Fst| (0.02930). We investigated interactions between 13 RBSDV proteins by two-hybrid screening assays and identified interactions between P5-1/P6, P6/P9-1, and P3/P6. We also observed self-interactive effects of P3, P6, P7-1, and P10. In short, we have proven that RBSDV populations were expanding and the overlapping region plays an important role in the genetic variation and evolution of RBSDV S5. Our results enable ongoing research into the evolutionary history of RBSDV-S5 with two partly overlapping ORFs.

Genome-Wide Variation in Potyviruses

Potyviruses (family Potyviridae, genus Potyvirus) are the result of an initial radiation event that occurred 6,600 years ago. The genus currently consists of 167 species that infect monocots or dicots, including domesticated and wild plants. Potyviruses are transmitted in a non-persistent way by more than 200 species of aphids. As indicated by their wide host range, worldwide distribution, and diversity of their vectors, potyviruses have an outstanding capacity to adapt to new hosts and environments. However, factors that confer adaptability are poorly understood. Viral RNA-dependent RNA polymerases introduce nucleotide substitutions that generate genetic diversity. We hypothesized that selection imposed by hosts and vectors creates a footprint in areas of the genome involved in host adaptation. Here, we profiled genomic and polyprotein variation in all species in the genus Potyvirus. Results showed that the potyviral genome is under strong negative selection. Accordingly, the genome and polyprotein sequence are remarkably stable. However, nucleotide and amino acid substitutions across the potyviral genome are not randomly distributed and are not determined by codon usage. Instead, substitutions preferentially accumulate in hypervariable areas at homologous locations across potyviruses. At a frequency that is higher than that of the rest of the genome, hypervariable areas accumulate non-synonymous nucleotide substitutions and sites under positive selection. Our results show, for the first time, that there is correlation between host range and the frequency of sites under positive selection. Hypervariable areas map to the N terminal part of protein P1, N and C terminal parts of helper component proteinase (HC-Pro), the C terminal part of protein P3, VPg, the C terminal part of NIb (RNA-dependent RNA polymerase), and the N terminal part of the coat protein (CP). Additionally, a hypervariable area at the NIb-CP junction showed that there is variability in the sequence of the NIa protease cleavage sites. Structural alignment showed that the hypervariable area in the CP maps to the N terminal flexible loop and includes the motif required for aphid transmission. Collectively, results described here show that potyviruses contain fixed hypervariable areas in key parts of the genome which provide mutational robustness and are potentially involved in host adaptation.

Eukaryotic Circular Rep-Encoding Single-Stranded DNA (CRESS DNA) Viruses: Ubiquitous Viruses With Small Genomes and a Diverse Host Range

While single-stranded DNA (ssDNA) was once thought to be a relatively rare genomic architecture for viruses, modern metagenomics sequencing has revealed circular ssDNA viruses in most environments and in association with diverse hosts. In particular, circular ssDNA viruses encoding a homologous replication-associated protein (Rep) have been identified in the majority of eukaryotic supergroups, generating interest in the ecological effects and evolutionary history of circular Rep-encoding ssDNA viruses (CRESS DNA) viruses. This review surveys the explosion of sequence diversity and expansion of eukaryotic CRESS DNA taxonomic groups over the last decade, highlights similarities between the well-studied geminiviruses and circoviruses with newly identified groups known only through their genome sequences, discusses the ecology and evolution of eukaryotic CRESS DNA viruses, and speculates on future research horizons.

Large-Scale Synonymous Substitutions in Cucumber Mosaic Virus RNA 3 Facilitate Amino Acid Mutations in the Coat Protein

The effect of large-scale synonymous substitutions in a small icosahedral, single-stranded RNA viral genome on virulence, viral titer, and protein evolution were analyzed. The coat protein (CP) gene of the Fny stain of cucumber mosaic virus (CMV) was modified. We created four CP mutants in which all the codons of nine amino acids in the 5' or 3' half of the CP gene were replaced by either the most frequently or the least frequently used synonymous codons in monocot plants. When the dicot host (Nicotiana benthamiana) was inoculated with these four CP mutants, viral RNA titers in uninoculated symptomatic leaves decreased, while all mutants eventually showed mosaic symptoms similar to those for the wild type. The codon adaptation index of these four CP mutants against dicot genes was similar to those of the wild-type CP gene, indicating that the reduction of viral RNA titer was due to deleterious changes of the secondary structure of RNAs 3 and 4. When two 5' mutants were serially passaged in N. benthamiana, viral RNA titers were rapidly restored but competitive fitness remained decreased. Although no nucleic acid changes were observed in the passaged wild-type CMV, one to three amino acid changes were observed in the synonymously mutated CP of each passaged virus, which were involved in recovery of viral RNA titer of 5' mutants. Thus, we demonstrated that deleterious effects of the large-scale synonymous substitutions in the RNA viral genome facilitated the rapid amino acid mutation(s) in the CP to restore the viral RNA titer.IMPORTANCE Recently, it has been known that synonymous substitutions in RNA virus genes affect viral pathogenicity and competitive fitness by alteration of global or local RNA secondary structure of the viral genome. We confirmed that large-scale synonymous substitutions in the CP gene of CMV resulted in decreased viral RNA titer. Importantly, when viral evolution was stimulated by serial-passage inoculation, viral RNA titer was rapidly restored, concurrent with a few amino acid changes in the CP. This novel finding indicates that the deleterious effects of large-scale nucleic acid mutations on viral RNA secondary structure are readily tolerated by structural changes in the CP, demonstrating a novel part of the adaptive evolution of an RNA viral genome. In addition, our experimental system for serial inoculation of large-scale synonymous mutants could uncover a role for new amino acid residues in the viral protein that have not been observed in the wild-type virus strains.

The Evolution of Molecular Compatibility between Bacteriophage ΦX174 and its Host

Viruses rely upon their hosts for biosynthesis of viral RNA, DNA and protein. This dependency frequently engenders strong selection for virus genome compatibility with potential hosts, appropriate gene regulation and expression necessary for a successful infection. While bioinformatic studies have shown strong correlations between codon usage in viral and host genomes, the selective factors by which this compatibility evolves remain a matter of conjecture. Engineered to include codons with a lesser usage and/or tRNA abundance within the host, three different attenuated strains of the bacterial virus ФX174 were created and propagated via serial transfers. Molecular sequence data indicate that biosynthetic compatibility was recovered rapidly. Extensive computational simulations were performed to assess the role of mutational biases as well as selection for translational efficiency in the engineered phage. Using bacteriophage as a model system, we can begin to unravel the evolutionary processes shaping codon compatibility between viruses and their host.

Tomato Spotted Wilt Virus NSs Protein Supports Infection and Systemic Movement of a Potyvirus and Is a Symptom Determinant

Plant viruses are inducers and targets of antiviral RNA silencing. To condition susceptibility, most plant viruses encode silencing suppressor proteins that interfere with antiviral RNA silencing. The NSs protein is an RNA silencing suppressor in orthotospoviruses, such as the tomato spotted wilt virus (TSWV). The mechanism of RNA silencing suppression by NSs and its role in virus infection and movement are poorly understood. Here, we cloned and tagged TSWV NSs and expressed it from a GFP-tagged turnip mosaic virus (TuMV-GFP) carrying either a wild-type or suppressor-deficient (AS9) helper component proteinase (HC-Pro). When expressed in cis, NSs restored pathogenicity and promoted systemic infection of suppressor-deficient TuMV-AS9-GFP in Nicotiana benthamiana and Arabidopsis thaliana. Inactivating mutations were introduced in NSs RNA-binding domain one. A genetic analysis with active and suppressor-deficient NSs, in combination with wild-type and mutant plants lacking essential components of the RNA silencing machinery, showed that the NSs insert is stable when expressed from a potyvirus. NSs can functionally replace potyviral HC-Pro, condition virus susceptibility, and promote systemic infection and symptom development by suppressing antiviral RNA silencing through a mechanism that partially overlaps that of potyviral HC-Pro. The results presented provide new insight into the mechanism of silencing suppression by NSs and its effect on virus infection.

Analysis of synonymous codon usage bias and phylogeny of coat protein gene in banana bract mosaic virus isolates

The banana is one of the world's most important livelihood crops. Banana plants are principally infected by four virus species, Banana bunchy top virus (genus Babuvirus), Cucumber mosaic virus (genus Cucumovirus), Banana streak virus (genus Badnavirus) and Banana bract mosaic virus (genus Potyvirus). The objective of this study is to understand the codon usage pattern and phylogeny of coat protein gene in different banana bract mosaic virus (BBrMV) isolates. The BBrMV Coat Protein (CP) gene was amplified from BBrMV infected banana plant samples collected from different districts of Tamil Nadu and Karnataka, India. Six new BBrMV isolates were submitted to National Center for Biotechnology Information. Phylogenetic analysis and codon usage indices were studied along with other isolates of BBrMV. Phylogenetic analysis of CP genes shows that most of BBrMV isolates are closely related to each other except KF385484.1 and KF385478.1. Relative codon usage patterns among different BBrMV isolates were calculated by software CodonW version 1.4.2. In BBrMV, codons with A-ended or U ended are the most preferential except the Leu and Gln whose optimized codons are CAG and UUG ending by G. The codon usage patterns of BBrMV isolates are principally influenced by mutational bias; however, compositional constraints along with mutational bias also play a major role.

Characterization of a Sorghum mosaic virus (SrMV) isolate in China

Sorghum mosaic virus (SrMV), a causal agent of the destructive sugarcane mosaic disease, has a global presence. An isolate of SrMV infecting a commercially-grown sugarcane plant was recovered from the Hainan province of China. The virions were visualized by an electron microscope, and the coat proteins (CPs) were sequenced by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and tandem mass spectrometry. Discrepancies between the CP predicted and actual amino acid sequences were noted, which confounded the phylogenetic assignment of the isolate. The apparent variations may have physiological effects on the pathogenicity and virulence of SrMV.

Phylogeny and synonymous codon usage pattern of Papaya ringspot virus coat protein gene in the sub-Himalayan region of north-east India

Sub-Himalayan West Bengal is favorable for the production of several fruits and vegetables. Papaya is one of the common plants cultivated in the area. Most of the papaya plants of the area are susceptible to Papaya ringspot virus (PRSV). Coat protein genes of 6 PRSV isolates of the area were sequenced following RT–PCR. Phylogenetic study of the PRSV isolates showed about 80%–90% similarity with Cuban isolates. The codon usage pattern of our isolates was also analyzed, along with several other isolates. PRSV isolates of our study showed a preference for 8 putative optimal codons. Correspondence analysis of the genes of different isolates along the first 2 major axes were done, as the first 2 axes contributed more in shaping codon usage pattern. In the phylogenetic tree constructed by the neighbour-joining method, our isolates clustered together with the east Indian, north Indian, and Bangladeshi isolates. The diversity and codon usage pattern of the PRSV isolates of different regions were studied, and it was observed that the codon usage pattern of PRSV isolates is probably influenced by translational selection along with mutational bias.

Molecular Genetic Analysis and Evolution of Segment 7 in Rice Black-Streaked Dwarf Virus in China

Rice black-streaked dwarf virus (RBSDV) causes maize rough dwarf disease or rice black-streaked dwarf disease and can lead to severe yield losses in maize and rice. To analyse RBSDV evolution, codon usage bias and genetic structure were investigated in 111 maize and rice RBSDV isolates from eight geographic locations in 2013 and 2014. The linear dsRNA S7 is A+U rich, with overall codon usage biased toward codons ending with A (A3s, S7-1: 32.64%, S7-2: 29.95%) or U (U3s, S7-1: 44.18%, S7-2: 46.06%). Effective number of codons (Nc) values of 45.63 in S7-1 (the first open reading frame of S7) and 39.96 in S7-2 (the second open reading frame of S7) indicate low degrees of RBSDV-S7 codon usage bias, likely driven by mutational bias regardless of year, host, or geographical origin. Twelve optimal codons were detected in S7. The nucleotide diversity (π) of S7 sequences in 2013 isolates (0.0307) was significantly higher than in 2014 isolates (0.0244, P = 0.0226). The nucleotide diversity (π) of S7 sequences in isolates from Jinan (0.0391) was higher than that from the other seven locations (P < 0.01). Only one S7 recombinant was detected in Baoding. RBSDV isolates could be phylogenetically classified into two groups according to S7 sequences, and further classified into two subgroups. S7-1 and S7-2 were under negative and purifying selection, with respective Ka/Ks ratios of 0.0179 and 0.0537. These RBSDV populations were expanding (P < 0.01) as indicated by negative values for Tajima's D, Fu and Li's D, and Fu and Li's F. Genetic differentiation was detected in six RBSDV subpopulations (P < 0.05). Absolute Fst (0.0790) and Nm (65.12) between 2013 and 2014, absolute Fst (0.1720) and Nm (38.49) between maize and rice, and absolute Fst values of 0.0085-0.3069 and Nm values of 0.56-29.61 among these eight geographic locations revealed frequent gene flow between subpopulations. Gene flow between 2013 and 2014 was the most frequent.

Recurrent positive selection and heterogeneous codon usage bias events leading to coexistence of divergent pigeon circoviruses

The capsid genes from 14 pigeon circovirus (PiCV) sequences, collected from Taiwan between 2009 and 2010, were sequenced and compared with 14 PiCV capsid gene sequences from GenBank. Based on pairwise comparison, PiCV strains from Taiwan shared 73.9-100% nucleotide identity and 72-100% amino acid identity with those of the 14 reported PiCV sequences. Phylogenetic analyses revealed that Taiwanese PiCV isolates can be grouped into two clades: clade 1 comprising isolates from Belgium, Australia, USA, Italy and China, and clade 2 showing close relation to isolates from Germany and France. Recurrent positive selection was detected in clade 1 PiCV lineages, which may contribute to the diversification of predominant PiCV sequences in Taiwan. Further observations suggest that synonymous codon usage variations between PiCV clade 1 and clade 2 may reflect the adaptive divergence on translation efficiency of capsid genes in infectious hosts. Variation in selective pressures acting on the evolutionary divergence and codon usage bias of both clades explains the regional coexistence of virus sequences congeners prevented from competitive exclusion within an island such as Taiwan. Our genotyping results also provide insight into the aetiological agents of PiCV outbreak in Taiwan and we present a comparative analysis of the central coding region of PiCV genome. From the sequence comparison results of 28 PiCVs which differs in regard to the geographical origin and columbid species, we identified conserved regions within the capsid gene that are likely to be suitable for primer selection and vaccine development.

System analysis of synonymous codon usage biases in archaeal virus genomes

Recent studies of geothermally heated aquatic ecosystems have found widely divergent viruses with unusual morphotypes. Archaeal viruses isolated from these hot habitats usually have double-stranded DNA genomes, linear or circular, and can infect members of the Archaea domain. In this study, the synonymous codon usage bias (SCUB) and dinucleotide composition in the available complete archaeal virus genome sequences have been investigated. It was found that there is a significant variation in SCUB among different Archaeal virus species, which is mainly determined by the base composition. The outcome of correspondence analysis (COA) and Spearman׳s rank correlation analysis shows that codon usage of selected archaeal virus genes depends mainly on GC richness of genome, and the gene׳s function, albeit with smaller effects, also contributes to codon usage in this virus. Furthermore, this investigation reveals that aromaticity of each protein is also critical in affecting SCUB of these viral genes although it was less important than that of the mutational bias. Especially, mutational pressure may influence SCUB in SIRV1, SIRV2, ARV1, AFV1, and PhiCh1 viruses, whereas translational selection could play a leading role in HRPV1׳s SCUB. These conclusions not only can offer an insight into the codon usage biases of archaeal virus and subsequently the possible relationship between archaeal viruses and their host, but also may help in understanding the evolution of archaeal viruses and their gene classification, and more helpful to explore the origin of life and the evolution of biology.

•

The SCUB of archaeal virus genes depends mainly on GC richness of genome.

•

The mutational pressure is the main factor that influences SCUB.

•

The aromaticity of each protein is also critical in affecting SCUB.

•

The translational selection could play a leading role in HRPV1׳s SCUB.

•

The mode is helpful to explore the origin of life and the evolution of biology.

Structure and assembly of filamentous bacteriophages

Filamentous bacteriophages are interesting paradigms in structural molecular biology, in part because of the unusual mechanism of filamentous phage assembly. During assembly, several thousand copies of an intracellular DNA-binding protein bind to each copy of the replicating phage DNA, and are then displaced by membrane-spanning phage coat proteins as the nascent phage is extruded through the bacterial plasma membrane. This complicated process takes place without killing the host bacterium. The bacteriophage is a semi-flexible worm-like nucleoprotein filament. The virion comprises a tube of several thousand identical major coat protein subunits around a core of single-stranded circular DNA. Each protein subunit is a polymer of about 50 amino-acid residues, largely arranged in an α-helix. The subunits assemble into a helical sheath, with each subunit oriented at a small angle to the virion axis and interdigitated with neighbouring subunits. A few copies of "minor" phage proteins necessary for infection and/or extrusion of the virion are located at each end of the completed virion. Here we review both the structure of the virion and aspects of its function, such as the way the virion enters the host, multiplies, and exits to prey on further hosts. In particular we focus on our understanding of the way the components of the virion come together during assembly at the membrane. We try to follow a basic rule of empirical science, that one should chose the simplest theoretical explanation for experiments, but be prepared to modify or even abandon this explanation as new experiments add more detail.

Evidence of pervasive biologically functional secondary structures within the genomes of eukaryotic single-stranded DNA viruses

Single-stranded DNA (ssDNA) viruses have genomes that are potentially capable of forming complex secondary structures through Watson-Crick base pairing between their constituent nucleotides. A few of the structural elements formed by such base pairings are, in fact, known to have important functions during the replication of many ssDNA viruses. Unknown, however, are (i) whether numerous additional ssDNA virus genomic structural elements predicted to exist by computational DNA folding methods actually exist and (ii) whether those structures that do exist have any biological relevance. We therefore computationally inferred lists of the most evolutionarily conserved structures within a diverse selection of animal- and plant-infecting ssDNA viruses drawn from the families Circoviridae, Anelloviridae, Parvoviridae, Nanoviridae, and Geminiviridae and analyzed these for evidence of natural selection favoring the maintenance of these structures. While we find evidence that is consistent with purifying selection being stronger at nucleotide sites that are predicted to be base paired than at sites predicted to be unpaired, we also find strong associations between sites that are predicted to pair with one another and site pairs that are apparently coevolving in a complementary fashion. Collectively, these results indicate that natural selection actively preserves much of the pervasive secondary structure that is evident within eukaryote-infecting ssDNA virus genomes and, therefore, that much of this structure is biologically functional. Lastly, we provide examples of various highly conserved but completely uncharacterized structural elements that likely have important functions within some of the ssDNA virus genomes analyzed here.

The evolution of codon usage in structural and non-structural viral genes: the case of Avian coronavirus and its natural host Gallus gallus

•

Avian coronavirus codon usage is closer to chicken for non-structural proteins.

•

Spike glycoprotein gene shares less preferred codons with chicken.

•

Avian coronavirus balances replication and immune camouflage for its codon usage.

To assess the codon evolution in virus–host systems, Avian coronavirus and its natural host Gallus gallus were used as a model. Codon usage (CU) was measured for the viral spike (S), nucleocapsid (N), nonstructural protein 2 (NSP2) and papain-like protease (PL^pro) genes from a diverse set of A. coronavirus lineages and for G. gallus genes (lung surfactant protein A, intestinal cholecystokinin, oviduct ovomucin alpha subunit, kidney vitamin D receptor and the ubiquitary beta-actin) for different A. coronavirus replicating sites. Relative synonymous codon usage (RSCU) trees accommodating all virus and host genes in a single topology showed a higher proximity of A. coronavirus CU to the respiratory tract for all genes. The codon adaptation index (CAI) showed a lower adaptation of S to G. gallus compared to NSP2, PL^pro and N. The effective number of codons (Nc) and GC_3% revealed that natural selection and genetic drift are the evolutionary forces driving the codon usage evolution of both A. coronavirus and G. gallus regardless of the gene being considered. The spike gene showed only one 100% conserved amino acid position coded by an A. coronavirus preferred codon, a significantly low number when compared to the three other genes (p < 0.0001). Virus CU evolves independently for each gene in a manner predicted by the protein function, with a balance between natural selection and mutation pressure, giving further molecular basis for the viruses’ ability to exploit the host's cellular environment in a concerted virus–host molecular evolution.