Analysis of the codon usage pattern of the RdRP gene of mycovirus infecting Aspergillus spp.
Virology Journal volume 16, Article number: 10 (2019)
Mycoviruses that infect fungi generally do not have a significant effect on the host and, instead, reduce the toxicity of the fungi. However, recent studies have shown that polymycovirus-1, a mycovirus that infects Aspergillus species known to cause disease in humans, is related to increased virulence of the fungus.
Comparative analysis was performed of RdRP gene codon usage patterns of Aspergillus fumigatus polymycovirus-1 (AfuPmV-1) and other mycoviruses known to infect Aspergillus spp. to examine the genetic characteristics of AfuPmV-1. In addition, codon usage analysis was performed to determine whether the nucleotide composition and codon usage characteristics of AfuPmV-1 were also present in other polymycoviruses and hypervirulence-related mycoviruses. Phylogenetic analysis was also performed to investigate their evolutionary relationship.
Analysis of nucleotide composition indicated that AfuPmV-1 had the highest GC content among analyzed mycoviruses and relative synonymous codon usage analysis indicated that all of the codons preferred by AfuPmV-1 ended with C or G, while codons ending with A or U were not observed. Moreover, the effective number of codons, the codon adaptation index, and correspondence analysis showed that AfuPmV-1 had greater codon preference compared with other mycoviruses and that AfuPmV-1 had relatively high adaptability to humans and fungi. These results were generally similar among polymycoviruses.
The codon usage pattern of AfuPmV-1 differs from other mycoviruses that infect Aspergillus spp. This difference may be related to the hypervirulence effect of AfuPmV-1. Analysis of AfuPmV-1 codon usage patterns could contribute to the identification and prediction of virulence effects of mycoviruses with similar genetic characteristics.
Mycoviruses are viruses that infect fungi and are known to be infectious to most fungal species. While mycoviruses have little or no influence on the fungal host in most cases, some have been shown to control the pathogenicity of the host by increasing or decreasing the virulence [1, 2]. To date, research on mycoviruses has mainly focused on plant pathogenic fungi. The hypovirulence effect of mycoviruses in fungi has led to the development of biological control agents to reduce the virulence of fungi on important crops and plant resources [2, 3]. There have been relatively few studies on the hypervirulence effect of mycoviruses, although some results have been reported recently. In 2015, the A78 virus, which infects the human pathogenic fungus, Aspergillus fumigatus, was reported to have significant mild hypervirulence effects in the moth species, Galleria mellonella , and similar virulence effects have been reported in Aspergillus fumigatus tetramycovirus-1 (AfuTmV-1) found in the same fungus . AfuTmV-1 was renamed AfuPmV-1, and its virus family name was changed from Tetramycoviridae to Polymycoviridae as more viruses showing similar characteristics to AfuTmV-1 were discovered, including BbPmV-1, which infects the insect pathogenic fungus, Beauveria bassiana  (Table 1).
Aspergillus is a fungus that belongs to the phylum Ascomycota and is related to various human diseases through infection, including aspergillosis on infection of the lungs, asthma, and allergies. Moreover, it is one of the fatal human pathogenic fungi, with a high mortality rate in patients with low immunity by opportunistic infection [7, 8]. Although studies have focused on polymycovirus-1 with regard to the increase in fungal pathogenicity in A. fumigatus, it is also important to examine other mycoviruses, such as partitivirus, chrysovirus, and victorivirus, which are infectious to other species in the genus Aspergillus, including A. fumigatus, A. foetidus, A. ochraceus, and A. niger. This study was performed to explore other mycoviruses that may increase the pathogenicity of the genus Aspergillus by comparison with polymycovirus-1, and to investigate the genetic characteristics of mycoviruses that enhance the toxicity of human pathogenic fungi.
Sequence data collection
The AfuPmV-1 genome consists of four capped double-stranded RNAs (dsRNAs), of which the largest dsRNA1 segment encodes the RNA-dependent RNA polymerase (RdRP) [5, 6]. The data were downloaded from the NCBI GenBank database (https://www.ncbi.nlm.nih.gov/genbank) and the complete nucleotide sequences of the RdRP coding region of mycoviruses that infect Aspergillus species were used. Eight sets of data published between 2006 and 2015 were studied in the analysis. Other polymycoviruses and hypervirulence-related mycoviruses data were downloaded from NCBI for comparison (Table 2).
Codon usage analysis
The programs CodonW (https://sourceforge.net/projects/codonw/) and CALcal (http://genomes.urv.es/CAIcal/) were used to conduct analyses of nucleotide composition, overall/local G + C content, relative synonymous codon usage (RSCU), effective number of codons (ENC), codon adaptation index (CAI), and correspondence analysis (COA) for each of the selected genetic data [9, 10]. Nucleotide composition analysis was conducted based on the overall frequency of nucleic acid occurrence (A%, C%, U%, and G%), total AU%/GC%, frequency of the third nucleic acid in synonymous codons (A3s, C3s, U3 s, and G3 s), and the GC3s values. The RSCU is a value represented as the ratio between predicted and observed usage rates of specific codons, assuming that all synonymous codons were used equally for an amino acid. An RSCU value of 1.0 indicates that all codons were used randomly or equally. The value would be greater (less) than 1.0 if certain codons were used more (less) frequently [11, 12]. In particular, codons with RSCU values ≥1.6 are considered over-represented codons, while those with RSCU values ≤0.6 are considered under-represented codons [13, 14]. The ENC is a simple and absolute method to measure codon usage bias in genes and genomes. The GC3s values are used in the calculation, and the ENC value ranges from 20 to 61. The ENC is inversely proportional to codon usage bias, and lower ENC values indicate stronger codon usage bias. A gene that does not show any codon usage bias would have an ENC value of 61. On the other hand, a gene that shows a high degree of codon usage bias would have an ENC value of 20, indicating that only one codon was used for an amino acid. Generally, ENC values < 35 are interpreted as showing strong codon usage bias [15, 16]. The CAI is the geometric mean of relative adaptation level and is a quantitative method to calculate the differences in codon usage bias against highly expressed known reference data. The CAI value ranges between 0 and 1, with values closer to 1 indicating a high degree of similarity between the reference data and the codon usage pattern and expression level, whereas smaller values indicate lower similarity in the codon usage pattern and expression level . The reference data of codon usage levels in humans (Homo sapiens) and fungus (Aspergillus fumigatus) were acquired from the codon usage database (http://www.kazusa.or.jp/codon/) for comparison . Furthermore, the RSCU values were examined by the COA method with the XLSTAT 2016 program for visualization of the CodonW results. The COA method is the preferred method in the area of multivariate analysis and represents the data as vectors consisting of rows and columns . Each individual data of the RdRP codon region was represented as a vector with 59 dimensions covering 59 codons except methionine (AUG) and tryptophan (UGG), which lack synonymous codons, as well as termination codons. A comparative analysis was performed including three groups: mycoviruses infecting Aspergillus (Group 1), polymycoviruses with a fungal host other than Aspergillus (Group 2), and newly reported mycoviruses that enhance host virulence (Group 3).
The phylogenetic relationships among the mycoviruses used for the analysis were inspected and phylogenetic analysis was conducted using the MEGA7 program (http://www.megasoftware.net) to infer the influence of evolutionary processes on codon usage patterns. RdRP coding sequences were aligned with the MUSCLE algorithm, and the phylogenetic tree was constructed by applying the Maximum Likelihood method and the Kimura 2-parameter substitution model. The robustness of the tree was verified with the bootstrap value set to 1000 .
Nucleotide composition features
Basic nucleotide composition analysis was conducted for the RdRP gene of mycoviruses infecting Aspergillus spp. (Table 3, Fig. 1). The mean and standard deviation (SD) of A%, C%, U%, and G% were 23.16 ± 4.00, 25.09 ± 5.13, 25.13 ± 3.54, and 26.63 ± 4.04, respectively. The mean and SD of total AU% and GC% were 48.29 ± 6.76 and 51.71 ± 6.76, respectively. The mean and SD of A3s, C3s, U3 s, G3 s, and GC3s were 22.34 ± 8.66, 34.72 ± 12.09, 34.26 ± 10.20, 32.15 ± 7.47, and 53.48 ± 13.45, respectively. AfuPmV-1 showed the lowest frequencies of nucleotides A and U and the highest frequency of nucleotide C among the Aspergillus-infecting mycoviruses included in the analysis. Similarly, the frequencies of the four types of nucleotides in the third position of synonymous codons were the lowest with A3s and U3 s and highest with C3s. As nucleotide C occurred relatively abundantly in AfuPmV-1, the GC content and GC3s indices were highest in this virus with values of 62.74 and 78.30, respectively (Fig. 1). The GC content of AfuPmV-1 was reported to be approximately 63% in a previous study of the entire or partial genome sequences of Polymycoviridae viruses . As a result of the comparison, Group 2 showed a similar pattern to that of AfuPmV-1. The GC content was between 59.97 and 61.98%, and the frequencies of nucleotides at the third position in the synonymous codons were higher in C3s and G3 s than in A3s and U3 s, respectively. Among these, BbPmV-1 is an experimentally reported virus that may be related to mild hypervirulence. The Group 3 results differed from those of Groups 1 and 2 (Table 3).
RSCU value and codon usage preference
The RSCU values of AfuPmV-1, AfuCV, and AfuPV-1 of Group1 were compared to inspect the codon usage bias according to virus species (Table 4, Fig. 2). Of codons related to the entire 18 amino acids, 18 codons were preferred in AfuPmv-1 of which 11 showed RSCU values ≥1.6. AfuCV showed 17 preferred codons, three of which had RSCU values ≥1.6. AfuPV had 20 preferred codons, three of which showed RSCU values ≥1.6. AfuPmV-1 showed similarities to AfuCV in codons CAC (His), CAG (Gln) and to AfuPV-1 in codons CUC (Leu) and AUC (Ile). The codons preferred by each virus individually were as follows. The codons solely preferred in AfuPmV-1 were GUC (Val), UAC (Tyr), GAC (Asp), UCC (Ser), CCC (Pro), ACC (Thr), GCC (Ala), UGC (Cys), CGC (Arg), and GGC (Gly) for 10 of the amino acids. Remarkably, all of these were C-ended codons (Fig. 2). Codons solely preferred in AfuCV were CUG (Leu), AUA (Ile), CCG (Pro), GCU (Ala), AGC (Ser), and AGG (Arg) for six of the amino acids, and in AfuPV-1, CUU (Leu), CAU (His), CAA (Gln), UCA (Ser), CCA (Pro), GCA (Ala), GCG (Ala), UGU (Cys), and CGU (Arg) for eight of the amino acids. The codons UUC (Phe), AAC (Asn), AAG (Lys), and GAG (Glu) showed similar preferences in all three viruses. The RSCU values and end nucleotide composition indicated that, in the RdRP coding region, C-ended codons were strongly preferred in AfuPmV-1 (15 of 18), G-ended codons were preferred in AfuCV (7 of 17), and U-ended codons were preferred in AfuPV-1 (8 of 20). Interestingly, there were no A/U-ended codons among the preferred codons of AfuPmV-1, indicating that AfuPmV-1 has a codon bias toward C- and G-ended codons.
Group 2 showed a similar codon usage pattern to that of AfuPmV-1, which indicates a preference for the C/G-ended codon. Group 3 also preferred C- or G-ended codons, but the four nucleic acids were distributed more evenly.
General codon usage pattern
The ENC value was calculated to quantitatively measure the magnitude of RdRP gene codon usage bias of the eight mycoviruses infecting Aspergillus spp. An ENC value < 35 indicates a strong codon usage bias. The results showed that the ENC value was lowest for AfuPmV-1 (40.67) and the other viruses showed ENC values > 50 (Table 3, Fig. 3). Taking into account the results of previous studies on RNA viruses in which the ENC values of Zaire Ebola virus, Chikungunya virus, Hepatitis C virus, and West Nile virus were 57.23, 55.56, 52.62, and 53.81, respectively , AfuPmV-1 appeared to have stronger codon usage bias relative to other viruses. The CAI was calculated to compare the adaptability of synonymous genetic codon usage in mutually different individuals, and the codon usage patterns were considered similar to the reference individual with CAI values closer to 1. This study referred to the CAI values of H. sapiens and A. fumigatus, which have ranges of 0.699–0.762 and 0.676–0.843, respectively. The mean and SD were 0.72 ± 0.02 and 0.74 ± 0.05 for H. sapiens and A. fumigatus, respectively. Remarkably, AfuPmV-1 showed the highest values for both H. sapiens and A. fumigatus with values of 0.762 and 0.843, respectively. These results indicated that AfuPmV-1 has the greatest similarity to the reference data in codon usage pattern and expression level, and that it has higher adaptability to human and fungal hosts compared with other mycoviruses. In Group 2, CAI and ENC values were similar to those of AfuPmV-1 (Table 3).
General trend of codon usage variation
To inspect the trends related to codon usage patterns of the Aspergillus-infecting viruses AfuPmV-1, AfuCV, and AfuPV-1, COA was performed with the RSCU values. Axis1, Axis2, Axis3, and Axis4 of the COA-RSCU explained 44.31, 18.81, 17.33, and 9.78% of the total variation, respectively. The results were based on 59 codons, excluding the three termination codons and methionine (AUG)/tryptophan (UGG) that do not have synonymous codons. Although there were exceptions in the COA results according to the third nucleotide in the codon, G- and C-ended codons formed one group and A- and U-ended codons formed another group (Fig. 4(a)). Moreover, the COA results for over-represented codons with RSCU values ≥1.6 showed that the codons strongly preferred individually by the three viruses formed separate groups (Fig. 4(b)). The observations verified that there were differences in the codon pattern preferences among the Aspergillus-infecting viruses.
Evolutionary relationship between mycoviruses
A phylogenetic tree was constructed to examine the phylogenetic relationships among the 14 mycoviruses, including AfuPmV-1. Polymycoviruses were grouped with AfuPmV-1. Polymycoviruses showed a relatively close relationship with alternaviruses (Fig. 5). To provide more information, the RdRP family for each sequence was examined from Pfam. AfuPmV-1 was assigned to RdRP_1 and all other polymycoviruses with similar codon patterns were assigned to the same RdRP_1 (pfam00680). Other mycoviruses infecting Aspergillus were classified as RdRP_1 or RdRP_4, although they were found in the same fungus (Table 5). This result suggests that the hypervirulent effects of AfuPmV-1 may be more affected by viral genome characteristics than by the effect of Aspergillus as a host.
The mechanism underlying the hypervirulence effect of polymycoviruses has yet to be determined. However, experimental studies have demonstrated the existence of mycoviruses with mild hypervirulence effects, and other mycoviruses with similar sequences are continuously being discovered. The pathogenic effects of pathogenic fungi on the hosts may be increased by infection with mycoviruses that show hypervirulence effects. Therefore, it is necessary to determine the genetic characteristics of mycoviruses with hypervirulence effects. The results of the present study showed that AfuPmV-1 has a high GC content, and all of the strongly preferred codons (RSCU value ≥1.6) ended with either a C or G nucleotide. The distinctive codon usage pattern of AfuPmV-1 compared to other mycoviruses that infect Aspergillus spp. may be related to its hypervirulence effect. These characteristics did not appear in all mycoviruses with hypervirulent effects, but were shared by polymycoviruses. Nucleotide composition and codon usage patterns of polymycoviruses may be useful in predicting hypervirulent effects of unidentified mycoviruses.
Aspergillus spp. are pathogenic fungi that cause various symptoms in humans. The hypervirulence effect of mycoviruses can increase the toxicity of Aspergillus spp. in human hosts, and thus increase the severity of symptoms. Here, AfuPmV-1 was shown to have distinct patterns in some codon usage indexes compared to other mycoviruses that infect Aspergillus spp. The distinctive codon usage pattern of AfuPmV-1 demonstrated in the present study indicated the need for monitoring of mycoviruses with similar characteristics. Research on mycoviruses has generally focused on their hypovirulence effects on fungi that infect plants. With the discovery of polymycoviruses, further research on the hypervirulence effects of mycoviruses is needed, particularly with regard to mechanisms of virulence control in mycoviruses, such as AfuPmV-1, which infects human pathogenic fungi.
Aspergillus fumigatus chrysovirus
Aspergillus fumigatus polymycovirus-1
Aspergillus fumigatus partitivirus-1
Aspergillus fumigatus tetramycovirus-1
Beauveria bassiana polymycovirus-1
Codon adaptation index
Effective number of codons
RNA-dependent RNA polymerase
Relative synonymous codon usage
Ghabrial SA, Suzuki N. Viruses of plant pathogenic fungi. Annu Rev Phytopathol. 2009;47:353–84.
Ghabrial SA, Caston JR, Jiang D, Nibert ML, Suzuki N. 50-plus years of fungal viruses. Virology. 2015;479:356–68.
Pearson MN, Beever RE, Boine B, Arthur K. Mycoviruses of filamentous fungi and their relevance to plant pathology. Mol Plant Pathol. 2009;10:115–28.
Ozkan S, Coutts Robert HA. Aspergillus fumigatus mycovirus causes mild hypervirulent effect on pathogenicity when tested on Galleria mellonella. Fungal Genet Biol. 2015;76:20–6.
Kanhayuwa L, Kotta-Loizou I, Ozkan S, Gunning AP, Goutts Robert HA. A novel mycovirus from Aspergillus fumigatus contains four unique dsRNAs as its genome and is infectious as dsRNA. Proc Natl Acad Sci. 2015;112:9100–5.
Kotta-Loizou I, Coutts Robert HA. Studies on the virome of the entomopathogenic fungus Beauveria bassiana reveal novel dsRNA elements and mild hypervirulence. PLoS Pathog. 2017. https://doi.org/10.1371/journal.ppat.1006183.
Bennett JW. An overview of the genus Aspergillus. In: Machida M, Gomi K, editors. Aspergillus: molecular biology and genomics. Poole: horizon scientific press; 2010. p. 1–17.
Kotta-Loizou I, Coutts Robert HA. Mycoviruses in Aspergilli: a comprehensive review. Front Microbiol. 2017;8:1–15.
Peden J. Analysis of codon usage, PhD [dissertation]. In: Department of Genetics: University of Nottingham; 1999.
Puigbo P, Bravo IG, Garcia-Vallve S. CAIcal: a combined set of tools to assess codon usage adaptation. Biol Direct. 2008:3–38.
Sharp PM, Li WH. Codon usage in regulatory genes in Escherichia coli does not reflect selection for ‘rare’ codons. Nucleic Acids Res. 1986;14:7737–49.
Sharp PM, Li WH. An evolutionary perspective on synonymous codon usage in unicellular organisms. J Mol Evol. 1986b;24:28–38.
Wong Emily HM, Smith DK, Rabadan R, Peiris M, Poon Leo LM. Codon usage bias and the evolution of influenza a viruses. Codon usage biases of influenza virus. BMC Evol Biol. 2010;10:253.
Butt AM, Nasrullah I, Tong Y. Genome-wide analysis of codon usage and influencing factors in chikungunya viruses. PLoS One. 2014. https://doi.org/10.1371/journal.pone.0090905.
Wright F. The ‘effective number of codons’ used in a gene. Gene. 1990;87:23–9.
Comeron JM, Aguade M. An evaluation of measures of synonymous codon usage bias. J Mol Evol. 1998;47:268–74.
Sharp PM, Li WH. The codon adaptation index-a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987;15:1281–95.
Nakamura Y, Gojobori T, Ikemura T. Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res. 2000;28:292.
Greenacre MJ. Theory and applications of correspondence analysis. London: Academic Press; 1984.
Kumar S, Stecher G, Tamra K. Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Cristina J, Fajardo A, Sonora M, Moratorio G, Musto H. A detailed comparative analysis of codon usage bias in Zika virus. Virus Res. 2016;223:147–52.
Bhatti MF, Bignell EM, Coutts Robert HA. Complete nucleotide sequences of two dsRNAs associated with a new partitivirus infecting Aspergillus fumigatus. Arch Virol. 2011;156:1677–80.
Jamal A, Bignell EM, Coutts Robert HA. Complete nucleotide sequences of four dsRNAs associated with a new chrysovirus infecting Aspergillus fumigatus. Virus Res. 2010;153:64–70.
Kozlakidis Z, Herrero N, Ozkan S, Kanhayuwa L, Jamal A, Bhatti MF, Coutts Robert HA. Sequence determination of a quadripartite dsRNA virus isolated from Aspergillus foetidus. Arch Virol. 2013;158:267–72.
Kozlakidis Z, Herrero N, Coutts Robert HA. The complete nucleotide sequence of a totivirus from Aspergillus foetidus. Arch Virol. 2013;158:263–6.
Kozlakidis Z, Herrero N, Ozkan S, Bhatti MF, Coutts Robert HA. A novel dsRNA element isolated from the Aspergillus foetidus mycovirus complex. Arch Virol. 2013;158:2625–8.
Kim JW, Choi EY, Kim YT. Intergeneric relationship between the Aspergillus ochraceous virus F and the Penicillium stoloniferum virus S. Virus Res. 2006;120:212–5.
Hammond TM, Andrewski MD, Roossinck MJ, Keller NP. Aspergillus mycoviruses are targets and suppressors of RNA silencing. Eukaryot Cell. 2008;7:350–7.
Mu F, Xie J, Cheng S, You MP, Barbetti MJ, Jia J, Wang Q, Cheng J, Fu Y, Chen T, Jiang D. Virome characterization of a collection of Sclerotinia sclerotiorum from Australia. Front Microbiol. 2018;8:2540.
Niu Y, Yuan Y, Mao J, Yang Z, Cao Q, Zhang T, Wang S, Liu D. Characterization of two novel mycoviruses from Penicillium digitatum and the related fungicide resistance analysis. Sci Rep. 2018;8:5513.
Okada R, Ichinose S, Takeshita K, Urayama S, Fukuhara T, Komatsu K, Arie T, Ishihara A, Egusa M, Kodama M, Moriyama H. Molecular characterization of a novel mycovirus in Alternaria alternata manifesting two-sided effects: Down-regulation of host growth and up-uregulation of host plant pathogenicity. Virology. 2018;519:23–32.
Lau Susanna KP, Lo George CS, Chow Franklin WN, Fan Rachel YY, Cai JJ, Yuen KY, Woo Patrick CY. Novel Partitivirus enhances virulence of and causes aberrant gene expression in Talaromyces marneffei. MBio. 2018. https://doi.org/10.1128/mBio.00947-18.
We’d like to thank those who made their invaluable data publicly available.
This work was supported by a grant from the National Research Foundation of Korea funded by the Korea government (MSIP) (No. 2016R1C1B2015511) and the Ministry of Education (No. 2017R1D1A1B03033413).
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Je, M., Kim, H. & Son, H.S. Analysis of the codon usage pattern of the RdRP gene of mycovirus infecting Aspergillus spp.. Virol J 16, 10 (2019). https://doi.org/10.1186/s12985-019-1115-y