Fusarivirus accessory helicases present an evolutionary link for viruses infecting plants and fungi

A significant number of mycoviruses have been identified that are related to plant viruses, but their evolutionary relationships are largely unexplored. A fusarivirus, Rhizoctonia solani fusarivirus 4 (RsFV4), was identified in phytopathogenic fungus Rhizoctonia solani (R. solani) strain XY74 co-infected by an alphaendornavirus. RsFV4 had a genome of 10,833 nt (excluding the poly-A tail), and consisted of four non-overlapping open reading frames (ORFs). ORF1 encodes an 825 aa protein containing a conserved helicase domain (Hel1). ORF3 encodes 1550 aa protein with two conserved domains, namely an RNA-dependent RNA polymerase (RdRp) and another helicase (Hel2). The ORF2 and ORF4 likely encode two hypothetical proteins (520 and 542 aa) with unknown functions. The phylogenetic analysis based on Hel2 and RdRp suggest that RsFV4 was positioned within the fusarivirus group, but formed an independent branch with three previously reported fusariviruses of R. solani. Notably, the Hel1 and its relatives were phylogenetically closer to helicases of potyviruses and hypoviruses than fusariviruses, suggesting fusarivirus Hel1 formed an evolutionary link between these three virus groups. This finding provides evidence of the occurrence of a horizontal gene transfer or recombination event between mycoviruses and plant viruses or between mycoviruses. Our findings are likely to enhance the understanding of virus evolution and diversity.     

Mycoviruses mediate mycotoxin regulation in Aspergillus ochraceus

To date, no demonstration of a direct correlation between the presence of mycoviruses and the quantitative or qualitative modulation of mycotoxins has been shown. In our study, we transfected a virus-free ochratoxin A (OTA)-producing isolate of Aspergillus ochraceus with purified mycoviruses from a different A. ochraceus isolate and from Penicillium aurantiogriseum. Among the mycoviruses tested, only Aspergillus ochraceus virus (AoV), a partitivirus widespread in A. ochraceus, caused a specific interaction that led to an overproduction of OTA, which is regulated by the European Commission and is the second most important contaminant of food and feed commodities. Gene expression analysis failed to reveal a specific viral upregulation of the mRNA of genes considered to play a role in the OTA biosynthetic pathway. Furthermore, AoOTApks1, a polyketide synthase gene considered essential for OTA production, is surprisingly absent in the genome of our OTA-producing isolate. The possible biological and evolutionary implications of the mycoviral regulation of mycotoxin production are discussed.     

Evolutionary genomics of mycovirus-related dsRNA viruses reveals cross-family horizontal gene transfer and evolution of diverse viral lineages

ABSTRACT: BACKGROUND: Double-stranded (ds) RNA fungal viruses are typically isometric single-shelled particles that are classified into three families, Totiviridae, Partitiviridae and Chrysoviridae, the members of which possess monopartite, bipartite and quadripartite genomes, respectively. Recent findings revealed that mycovirus-related dsRNA viruses are more diverse than previously recognized. Although an increasing number of viral complete genomic sequences have become available, the evolution of these diverse dsRNA viruses remains to be clarified. This is particularly so since there is little evidence for horizontal gene transfer (HGT) among dsRNA viruses. RESULTS: In this study, we report the molecular properties of two novel dsRNA mycoviruses that were isolated from a field strain of Sclerotinia sclerotiorum, Sunf-M: one is a large monopartite virus representing a distinct evolutionary lineage of dsRNA viruses; the other is a new member of the family Partitiviridae. Comprehensive phylogenetic analysis and genome comparison revealed that there are at least ten monopartite, three bipartite, one tripartite and three quadripartite lineages in the known dsRNA mycoviruses and that the multipartite lineages have possibly evolved from different monopartite dsRNA viruses. Moreover, we found that homologs of the S7 Domain, characteristic of members of the genus phytoreovirus in family Reoviridae are widely distributed in diverse dsRNA viral lineages, including chrysoviruses, endornaviruses and some unclassified dsRNA mycoviruses. We further provided evidence that multiple HGT events may have occurred among these dsRNA viruses from different families. CONCLUSIONS: Our study provides an insight into the phylogeny and evolution of mycovirus-related dsRNA viruses and reveals that the occurrence of HGT between different virus species and the development of multipartite genomes during evolution are important macroevolutionary mechanisms in dsRNA viruses.     

Codivergence of mycoviruses with their hosts

Background

The associations between pathogens and their hosts are complex and can result from any combination of evolutionary events such as codivergence, switching, and duplication of the pathogen. Mycoviruses are RNA viruses which infect fungi and for which natural vectors are so far unknown. Thus, lateral transfer might be improbable and codivergence their dominant mode of evolution. Accordingly, mycoviruses are a suitable target for statistical tests of virus-host codivergence, but inference of mycovirus phylogenies might be difficult because of low sequence similarity even within families.

Methodology

We analyzed here the evolutionary dynamics of all mycovirus families by comparing virus and host phylogenies. Additionally, we assessed the sensitivity of the co-phylogenetic tests to the settings for inferring virus trees from their genome sequences and approximate, taxonomy-based host trees.

Conclusions

While sequence alignment filtering modes affected branch support, the overall results of the co-phylogenetic tests were significantly influenced only by the number of viruses sampled per family. The trees of the two largest families, Partitiviridae and Totiviridae, were significantly more similar to those of their hosts than expected by chance, and most individual host-virus links had a significant positive impact on the global fit, indicating that codivergence is the dominant mode of virus diversification. However, in this regard mycoviruses did not differ from closely related viruses sampled from non-fungus hosts. The remaining virus families were either dominated by other evolutionary modes or lacked an apparent overall pattern. As this negative result might be caused by insufficient taxon sampling, the most parsimonious hypothesis still is that host-parasite evolution is basically the same in all mycovirus families. This is the first study of mycovirus-host codivergence, and the results shed light not only on how mycovirus biology affects their co-phylogenetic relationships, but also on their presumable host range itself.     

Deep sequencing of mycovirus‐derived small RNAs from Botrytis species

RNA silencing is an ancient regulatory mechanism operating in all eukaryotic cells. In fungi, it was first discovered in Neurospora crassa, although its potential as a defence mechanism against mycoviruses was first reported in Cryphonectria parasitica and, later, in several fungal species. There is little evidence of the antiviral potential of RNA silencing in the phytopathogenic species of the fungal genus Botrytis. Moreover, little is known about the RNA silencing components in these fungi, although the analysis of public genome databases identified two Dicer‐like genes in B. cinerea, as in most of the ascomycetes sequenced to date. In this work, we used deep sequencing to study the virus‐derived small RNA (vsiRNA) populations from different mycoviruses infecting field isolates of Botrytis spp. The mycoviruses under study belong to different genera and species, and have different types of genome [double‐stranded RNA (dsRNA), (+)single‐stranded RNA (ssRNA) and (–)ssRNA]. In general, vsiRNAs derived from mycoviruses are mostly of 21, 20 and 22 nucleotides in length, possess sense or antisense orientation, either in a similar ratio or with a predominance of sense polarity depending on the virus species, have predominantly U at their 5′ end, and are unevenly distributed along the viral genome, showing conspicuous hotspots of vsiRNA accumulation. These characteristics reveal striking similarities with vsiRNAs produced by plant viruses, suggesting similar pathways of viral targeting in plants and fungi. We have shown that the fungal RNA silencing machinery acts against the mycoviruses used in this work in a similar manner independent of their viral or fungal origin.     

Double-Stranded RNA Mycovirus Infection of Aspergillus fumigatus Is Not Dependent on the Genetic Make-Up of the Host

Aspergillus fumigatus is a fungus that causes opportunistic infections in immunocompromised patients, with high morbidity and mortality. In its turn, A. fumigatus can become infected with mycoviruses. Most mycoviruses have a dsRNA genome and can cause fungal hypovirulence. For that reason, mycoviruses could theoretically be used as therapeutic tools to combat fungal infections. We determined if a certain genetic make-up of A. fumigatus was associated with the presence of mycoviruses in 86 clinical A. fumigatus isolates. Mycovirus screening was performed by isolating dsRNA from mycelial cultures using a Trizol/Chloroform method. The genetic relatedness of dsRNA infected A. fumigatus was determined by cell surface protein (CSP) typing and determination of the mating type. Sixteen (18.6%) of the 86 clinical A. fumigatus isolates contained dsRNA. The A. fumigatus collection could be divided into 11 different CSP types. DsRNA infected A. fumigatus isolates had similar CSP types as non-infected isolates. In both cases, the CSP types t01, t02, t03 and t04 were the most prevalent and the distribution comparable to the CSP types observed in other Dutch collections. Mating types MAT1-1 and MAT1-2 were evenly distributed among all A. fumigatus strains, regardless of CSP type. No difference was observed in mycovirus infections between MAT1-1 and MAT1-2 isolates. DsRNA mycovirus infections in A. fumigatus are not related to either CSP or mating type and therefore represent an interesting future therapeutic tool to combat fungal infections.     

Molecular characterization of a dsRNA mycovirus, Fusarium graminearum virus-DK21, which is phylogenetically related to hypoviruses but has a genome organization and gene expression strategy resembling those of plant potex-like viruses

Fusarium graminearum causes a serious scab disease of small grains in Korea. The nucleotide sequence of the genomic RNA of a double-stranded RNA (dsRNA) virus, Fusarium graminearum virus-DK21 (FgV-DK21), from F. graminearum strain DK21, which is associated with hypovirulence in F. graminearum, was determined and compared to the genome sequences of other mycoviruses, including Cryponectria hypoviruses. The FgV-DK21 dsRNA consists of 6,624 nucleotides, excluding the 3'-terminal poly(A) tail. The viral genome has 53- and 46-nucleotide 5' and 3' untranslated regions (UTRs), respectively, and five putative open reading frames. A phylogenetic analysis of the deduced amino acid sequence of ORF1, which encodes a putative RNA-dependent RNA polymerase, and those of other mycoviruses revealed that this organism forms a distinct virus clade with other hypoviruses, and is more distantly related to other mycoviruses (3.8 to 24.0% identity). However, pairwise sequence comparisons of the nucleotide and deduced amino acid sequences of ORFs 2 through 5 revealed no close relationships to other protein sequences currently available in GenBank. Analyses of RNA accumulation by Northern blot and primer extension indicated that these putative gene products are expressed from at least two different subgenomic RNAs (sgRNAs), in contrast to the cases in other hypoviruses. This study suggests the existence of a new, as yet unassigned, genus of mycoviruses that exhibits a potex-like genome organization and sgRNA accumulation.     

Evolutionary genomics of the Fox genes: Origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.     

Non-detection of mycoviruses in amphibian chytrid fungus (Batrachochytrium dendrobatidis) from Australia

Mycoviruses may influence the pathogenicity of disease-causing fungi. Although mycoviruses have been found in some chytrid fungi, limited testing has not detected them in Batrachochytrium dendrobatidis (Bd), the cause of the devastating amphibian disease, chytridiomycosis. Here we conducted a survey for mycovirus presence in 38 Bd isolates from Australia (n = 31), Brazil (n = 5) and South Korea (n = 2) with a combination of modern high-throughput sequencing and conventional dsRNA cellulose chromatography. Mycoviruses were not detected in any isolates. This result was unexpected, given the long evolutionary history of Bd, as well as the high prevalence of mycoviruses in related fungal species. Given our widespread sampling in Australia and the limited number of Bd introductions, we suggest that mycoviruses are uncommon or absent from Australian Bd. Testing more isolates from regions where Bd originated, as well as regions with high diversity or low fungal virulence may identify mycoviruses that could aid in disease control.     

Phylogeny and molecular evolution of the rbcL gene of St genome in Elymus sensu lato (Poaceae: Triticeae)

Analysis of the patterns and levels of diversity in duplicate gene not only traces evolutionary history of polyploids, but also provides insight into how the evolutionary process differs between lineages and between homoeologous loci within lineages. Elymus sensu lato is a group of allopolyploid species, which share a common St genome and with the different combinations of H, Y, P, and W genomes. To estimate the evolutionary process of the rbcL gene in species of Elymus s. l. and its putative dioploid relatives, 74 sequences were obtained from 21 species of Elymus s. l. together with 24 diploid taxa representing 19 basic genomes in Triticeae. Phylogeny and sequence diversity pattern analysis suggested that (1) species of Pseudoroegneria (Nevski) Á. Löve might serve as the maternal donor of the species of Elymus s. l; (2) differentiation of St genome were shown in the species of Elymus s. l. following polyploidy event; (3) divergences within the species might associate with geographic diversity and morphological variability; (4) differences in the levels and patterns of nucleotide diversity of the rbcL gene implied that the St genome lineages in the species of Elymus s. l. have differently evolutionary potentials.     

Whole genome sequencing reveals the genomic diversity,taxonomic classification,and evolutionary relationships of the genus Nocardia

Nocardia is a complex and diverse genus of aerobic actinomycetes that cause complex clinical presentations, which are difficult to diagnose due to being misunderstood. To date, the genetic diversity, evolution, and taxonomic structure of the genus Nocardia are still unclear. In this study, we investigated the pan-genome of 86 Nocardia type strains to clarify their genetic diversity. Our study revealed an open pan-genome for Nocardia containing 265,836 gene families, with about 99.7% of the pan-genome being variable. Horizontal gene transfer appears to have been an important evolutionary driver of genetic diversity shaping the Nocardia genome and may have caused historical taxonomic confusion from other taxa (primarily Rhodococcus, Skermania, Aldersonia, and Mycobacterium). Based on single-copy gene families, we established a high-accuracy phylogenomic approach for Nocardia using 229 genome sequences. Furthermore, we found 28 potentially new species and reclassified 16 strains. Finally, by comparing the topology between a phylogenomic tree and 384 phylogenetic trees (from 384 single-copy genes from the core genome), we identified a novel locus for inferring the phylogeny of this genus. The dapb1 gene, which encodes dipeptidyl aminopeptidase BI, was far superior to commonly used markers for Nocardia and yielded a topology almost identical to that of genome-based phylogeny. In conclusion, the present study provides insights into the genetic diversity, contributes a robust framework for the taxonomic classification, and elucidates the evolutionary relationships of Nocardia. This framework should facilitate the development of rapid tests for the species identification of highly variable species and has given new insight into the behavior of this genus.     

The Motilin Gene Evolved a New Function in Kangaroo Rats and Kangaroo Mice (Dipodomyinae)

The motilin receptor gene was lost in the ancestral lineage of rodents. Subsequently, the gene encoding its ligand, motilin, has experienced different evolutionary fates. Previous genomic analyses had shown that the motilin gene (MLN) became a pseudogene independently in the lineages leading to the guinea pig and the common ancestor of the mouse and rat, yet an intact, and thus potentially functional, open reading frame for the MLN was preserved in the Dipodomys ordii genome. As only a single MLN haplotype from D. ordii was available, and this sequence is from a low coverage draft genome, it is possible that the intact MLN found in the draft kangaroo rat genome is an artifact, or represents an intermediate in the process of becoming a pseudogene. In order to establish whether an intact MLN is retained in kangaroo rats despite the loss of its specific receptor, and to investigate the evolutionary mechanisms underlying the retention of this gene sequence, we isolated MLN sequences from species that represent the diversity of the Dipodomyinae [the monophyletic Dipodomyinae subfamily consists of two genera: Dipodomys (kangaroo rats) and Microdipodops (kangaroo mice)]. The results demonstrate that the MLN sequence is well conserved in Dipodomyinae, and it codes for a predicted motilin peptide sequence possessing a conserved N-terminal pharmacophore and the potential to be processed and secreted as a hormone. The observations that the MLN evolved as a functional gene during the radiation of the Dipodomyinae, species that have lost their original motilin receptor, suggest that the MLN has undergone a lineage-specific physiological adaptation to a new function.     

Aspergillus mycoviruses are targets and suppressors of RNA silencing

RNA silencing can function as a virus defense mechanism in a diverse range of eukaryotes, and many viruses are capable of suppressing the silencing machinery targeting them. However, the extent to which this occurs between fungal RNA silencing and mycoviruses is unclear. Here, three Aspergillus dsRNA mycoviruses were partially characterized, and their relationship to RNA silencing was investigated. Aspergillus virus 1816 is related to Agaricus bisporus white button mushroom virus 1 and suppresses RNA silencing through a mechanism that alters the level of small interfering RNA. Aspergillus virus 178 is related to RNA virus L1 of Gremmeniella abietina and does not appear to affect RNA silencing. The third virus investigated, Aspergillus virus 341, is distantly related to Sphaeropsis sapinea RNA virus 2. Detection of mycovirus-derived siRNA from this mycovirus demonstrates that it is targeted for degradation by the Aspergillus RNA silencing machinery. Thus, our results indicate that Aspergillus mycoviruses are both targets and suppressors of RNA silencing. In addition, they suggest that the morphological and physiological changes associated with some mycoviruses could be a result of their antagonistic relationship with RNA silencing.     

Mycoviruses of Botrytis cinerea isolates from different hosts

Botrytis cinerea (teleomorph Botryotinia fuckeliana) is a necrotrophic plant pathogenic fungus that causes grey mould and enormous economic losses worldwide in different crops. Control of B. cinerea is difficult due to the appearance of fungicide‐resistant isolates, and the diversity in virulence due to genetic variability and, perhaps, the infection with mycoviruses or fungal viruses. The discovery of mycoviruses and their possible application as biocontrol agents, as well as their use as tools to study the plant–pathogen interaction, has encouraged their study in B. cinerea. Herein, we have analysed the occurrence of mycoviruses in Spanish B. cinerea isolates to approach a better understanding of the interactions among viruses, fungi and plants in this pathosystem. Fifty‐five percent of the B. cinerea isolates analysed contained double‐stranded RNA (dsRNA) elements, and the number of dsRNA elements, their relative concentration and size were variable among isolates. Some of these dsRNAs were related to the presence of virus like rod or isometric particles, and to cellular degeneration and malformed mitochondria. We have also demonstrated that a 3 kb dsRNA present in 55% of the isolates having dsRNA elements was a mycovirus genome. Partial sequence of that mycovirus presented high identity in nucleotide and amino acid sequence with Botrytis cinerea mitovirus 1 (BcMV1). Analysis of the genetic distance within Spanish BcMV1 sequences showed the existence of different isolates of this mitovirus inside the Spanish B. cinerea population analysed. This is the first report of the variability of dsRNA elements and the partial genome sequence of a mitovirus associated with Spanish B. cinerea isolates and the genetic diversity within Spanish isolates of BcMV1.     

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications(WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin(Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species.Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression.Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.     

Loss/retention and evolution of NBS-encoding genes upon whole genome triplication of Brassica rapa

A genome triplication took place in the ancestor of Brassiceae species after the split of the Arabidopsis lineage. The postfragmentation and shuffling of the genome turned the ancestral hexaploid back to diploids and caused the radiation of Brassiceae species. The course of speciation was accompanied by the loss of duplicate genes and also influenced the evolution of retained genes. Of all the genes, those encoding NBS domains are typical R genes that confer resistance to invading pathogens. In this study, using the genome of Arabidopsis thaliana as a reference, we examined the loss/retention of orthologous NBS-encoding loci in the tripled Brassica rapa genome and discovered differential loss/retention frequencies. Further analysis indicated that loci of different retention ratios showed different evolutionary patterns. The loci of classesII and III (maintaining two and three syntenic loci, respectively, multi-loci) show sharper expansions by tandem duplications, have faster evolutionary rates and have more potential to be associated with novel gene functions. On the other hand, the loci that are retained at the minimal rate (keeping only one locus, class I, single locus) showed opposite patterns. Phylogenetic analysis indicated that recombination and translocation events were common among multi-loci in B. rapa, and differential evolutionary patterns between multi- and single-loci are likely the consequence of recombination. Investigations towards other gene families demonstrated different evolutionary characteristics between different gene families. The evolution of genes is more likely determined by the property of each gene family, and the whole genome triplication provided only a specific condition.