Complete genome sequence and evolution analysis of Eimeria stiedai RNA virus 1, a novel member of the family Totiviridae

A virus associated with the zoonotic pathogen Plasmodium knowlesi causing human malaria is a member of a diverse and unclassified viral taxon

The Apicomplexa are a phylum of single-celled eukaryotes that can infect humans and include the mosquito-borne parasite Plasmodium, the cause of malaria. Viruses that infect non-Plasmodium spp. disease-causing protozoa affect the pathogen life cycle and disease outcomes. However, only one RNA virus (Matryoshka RNA virus 1) has been identified in Plasmodium, and none have been identified in zoonotic Plasmodium species. The rapid expansion of the known RNA virosphere via metagenomic sequencing suggests that this dearth is due to the divergent nature of RNA viruses that infect protozoa. We leveraged newly uncovered data sets to explore the virome of human-infecting Plasmodium species collected in Sabah, east (Borneo) Malaysia. From this, we identified a highly divergent RNA virus in two human-infecting P. knowlesi isolates that is related to the unclassified group 'ormycoviruses'. By characterizing 15 additional ormycoviruses identified in the transcriptomes of arthropods, we show that this group of viruses exhibits a complex ecology as noninfecting passengers at the arthropod-mammal interface. With the addition of viral diversity discovered using the artificial intelligence-based analysis of metagenomic data, we also demonstrate that the ormycoviruses are part of a diverse and unclassified viral taxon. This is the first observation of an RNA virus in a zoonotic Plasmodium species. By linking small-scale experimental data to advances in large-scale virus discovery, we characterize the diversity and confirm the putative genomic architecture of an unclassified viral taxon. This approach can be used to further explore the virome of disease-causing Apicomplexa and better understand how protozoa-infecting viruses may affect parasite fitness, pathobiology, and treatment outcomes.

A virus associated with the zoonotic pathogen Plasmodium knowlesi causing human malaria is a member of a diverse and unclassified viral taxon

Apicomplexa are single-celled eukaryotes that can infect humans and include the mosquito-borne parasite Plasmodium, the cause of malaria. Increasing rates of drug resistance in human-only Plasmodium species are reducing the efficacy of control efforts and antimalarial treatments. There are also rising cases of P. knowlesi, the only zoonotic Plasmodium species that causes severe disease and death in humans. Thus, there is a need to develop additional innovative strategies to combat malaria. Viruses that infect non-Plasmodium spp. Disease-causing protozoa have been shown to affect pathogen life cycle and disease outcomes. However, only one virus (Matryoshka RNA virus 1) has been identified in Plasmodium, and none have been identified in zoonotic Plasmodium species. The rapid expansion of the known RNA virosphere using structure- and artificial intelligence-based methods suggests that this dearth is due to the divergent nature of RNA viruses that infect protozoa. We leveraged these newly uncovered data sets to explore the virome of human-infecting Plasmodium species collected in Sabah, east (Borneo) Malaysia. We identified a highly divergent RNA virus in two human-infecting P. knowlesi isolates that is related to the unclassified group 'ormycoviruses'. By characterising fifteen additional ormycoviruses identified in the transcriptomes of arthropods we show that this group of viruses exhibits a complex ecology at the arthropod-mammal interface. Through the application of artificial intelligence methods, we then demonstrate that the ormycoviruses are part of a diverse and unclassified viral taxon. This is the first observation of an RNA virus in a zoonotic Plasmodium species. By linking small-scale experimental data to large-scale virus discovery advances, we characterise the diversity and genomic architecture of an unclassified viral taxon. This approach should be used to further explore the virome of disease-causing Apicomplexa and better understand how protozoa-infecting viruses may affect parasite fitness, pathobiology, and treatment outcomes.

Transcriptome mining of RNA viruses (family Totiviridae ) in Eimeria necatrix and Eimeria stiedai

Coccidian protozoa from the genus Eimeria are widespread parasites of vertebrates, causing serious disease (coccidiosis) and economic loss most notably in poultry. Several species of Eimeria are themselves infected with small RNA viruses assigned to the family Totiviridae . In this study, the sequences of two such viruses were newly determined, one of which represents the first complete protein-coding sequence of a virus from E. necatrix , an important pathogen of chickens, and the other of which is from E. stiedai , an important pathogen of rabbits. Sequence features of the newly identified viruses, compared with those of ones reported previously, provide several insights. Phylogenetic analyses suggest that these eimerian viruses constitute a well-demarcated clade, probably deserving of recognition as a distinct genus.

Elimination of LRVs Elicits Different Responses in Leishmania spp

Leishmaniaviruses (LRVs) have been demonstrated to enhance progression of leishmaniasis, a vector-transmitted disease with a wide range of clinical manifestations that is caused by flagellates of the genus Leishmania. Here, we used two previously proposed strategies of the LRV ablation to shed light on the relationships of two Leishmania spp. with their respective viral species (L. guyanensis, LRV1 and L. major, LRV2) and demonstrated considerable difference between two studied systems. LRV1 could be easily eliminated by the expression of exogenous capsids regardless of their origin (the same or distantly related LRV1 strains, or even LRV2), while LRV2 was only partially depleted in the case of the native capsid overexpression. The striking differences were also observed in the effects of complete viral elimination with 2'C-methyladenosine (2-CMA) on the transcriptional profiles of these two Leishmania spp. While virtually no differentially expressed genes were detected after the LRV1 removal from L. guyanensis, the response of L. major after ablation of LRV2 involved 87 genes, the analysis of which suggested a considerable stress experienced even after several passages following the treatment. This effect on L. major was also reflected in a significant decrease of the proliferation rate, not documented in L. guyanensis and naturally virus-free strain of L. major. Our findings suggest that integration of L. major with LRV2 is deeper compared with that of L. guyanensis with LRV1. We presume this determines different effects of the viral presence on the Leishmania spp. infections.

IMPORTANCELeishmania spp. represent human pathogens that cause leishmaniasis, a widespread parasitic disease with mild to fatal clinical manifestations. Some strains of leishmaniae bear leishmaniaviruses (LRVs), and this has been shown to aggravate disease course. We investigated the relationships of two distally related Leishmania spp. with their respective LRVs using different strategies of virus removal. Our results suggest the South American L. guyanensis easily loses its virus with no important consequences for the parasite in the laboratory culture. Conversely, the Old-World L. major is refractory to virus removal and experiences a prominent stress if this removal is nonetheless completed. The drastically different levels of integration between the studied Leishmania spp. and their viruses suggest distinct effects of the viral presence on infections in these species of parasites.

Analyses of Leishmania-LRV Co-Phylogenetic Patterns and Evolutionary Variability of Viral Proteins

Leishmania spp. are important pathogens causing a vector-borne disease with a broad range of clinical manifestations from self-healing ulcers to the life-threatening visceral forms. Presence of Leishmania RNA virus (LRV) confers survival advantage to these parasites by suppressing anti-leishmanial immunity in the vertebrate host. The two viral species, LRV1 and LRV2 infect species of the subgenera Viannia and Leishmania, respectively. In this work we investigated co-phylogenetic patterns of leishmaniae and their viruses on a small scale (LRV2 in L. major) and demonstrated their predominant coevolution, occasionally broken by intraspecific host switches. Our analysis of the two viral genes, encoding the capsid and RNA-dependent RNA polymerase (RDRP), revealed them to be under the pressure of purifying selection, which was considerably stronger for the former gene across the whole tree. The selective pressure also differs between the LRV clades and correlates with the frequency of interspecific host switches. In addition, using experimental (capsid) and predicted (RDRP) models we demonstrated that the evolutionary variability across the structure is strikingly different in these two viral proteins.

Expression Analysis and Serodiagnostic Potential of Microneme Proteins 1 and 3 in Eimeria stiedai

Eimeria stiedai is an apicomplexan protozoan parasite that invades the liver and bile duct epithelial cells in rabbits and causes severe hepatic coccidiosis, resulting in significant economic losses in the domestic rabbit industry. Hepatic coccidiosis lacks the typical clinical symptoms and there is a lack of effective premortem tools to timely diagnose this disease. Therefore, in the present study we cloned and expressed the two microneme proteins i.e., microneme protein 1 (EsMIC1) and microneme protein 3 (EsMIC3) from E. stiedai and used them as recombinant antigens to develop a serodiagnostic method for an effective diagnosis of hepatic coccidiosis. The cDNAs encoding EsMIC1 and EsMIC3 were cloned and the mRNA expression levels of these two genes at different developmental stages of E. stiedai were determined by quantitative real-time PCR analysis (qRT-PCR). The immunoreactivity of recombinant EsMIC1 (rEsMIC1) and EsMIC3 (rEsMIC3) proteins were detected by Western blotting, and indirect enzyme-linked immunosorbent assays (ELISAs) based on these two recombinant antigens were established to evaluate their serodiagnostic potential. Our results showed that the proteins encoded by the ORFs of EsMIC1 (711 bp) and EsMIC3 (891 bp) were approximately 25.89 and 32.39 kDa in predicted molecular weight, respectively. Both EsMIC1 and EsMIC3 showed the highest mRNA expression levels in the merozoites stage of E. stiedai. Western blotting analysis revealed that both recombinant proteins were recognized by E. stiedai positive sera, and the indirect ELISAs using rEsMIC1 and rEsMIC3 were developed based on their good immunoreactivity, with 100% (48/48) sensitivity and 97.9% (47/48) specificity for rEsMIC1 with 100% (48/48) sensitivity and 100% (48/48) specificity for rEsMIC3, respectively. Moreover, rEsMIC1- and rEsMIC3-based indirect ELISA were able to detect corresponding antibodies in sera at days 6, 8, and 10 post E. stiedai infection, with the highest positive diagnostic rate (62.5% (30/48) for rEsMIC1 and 66.7% (32/48) for rEsMIC3) observed at day 10 post infection. Therefore, both EsMIC1 and EsMIC3 can be used as potential serodiagnostic candidate antigens for hepatic coccidiosis caused by E. stiedai.

Evolutionary origin of 2A-like sequences in Totiviridae genomes

In recent years there has been a significant increase in the number of new species potentially belonging to the Totiviridae family. Most of these new viruses have not yet been covered by the Committee on Taxonomy of Viruses (ICTV) official classification. In this study, a phylogenetic analysis including new sequences of Totiviridae candidates revealed a clade including Giardiavirus and a great diversity of new totiviruses, which infect arthropods, protozoa and mollusc. This expanded Giardiavirus clade comprises two monophyletic groups, one of them including Giardia lamblia virus (GLV) grouped with viruses that infect arthropods and vertebrates (GLV-like group), and the other includes the previously proposed Artivirus group (IMNV-like group). A screening of the members of the GLV-like group in search of genomic elements already described in IMNV-like group revealed the existence of sites with a high propensity to become 2 A-like oligopeptides, mainly in a specific subgroup of arthropod viruses, suggesting that these viruses preserved ancestral characteristics. The existence of these "pseudo 2 A-sites" associated to phylogenetic reconstruction indicates that these sequences appear at a decisive stage for viral evolution. If they are changed to functional 2 A-like sequences, an irreversible route to increase the genome complexity will be initiated.

Development and evaluation of novel recombinant adenovirus-based vaccine candidates for infectious bronchitis virus and Mycoplasma gallisepticum in chickens

Avian infectious bronchitis caused by the infectious bronchitis virus (IBV), and mycoplasmosis caused by Mycoplasma gallisepticum (MG) are two major respiratory diseases in chickens that have resulted in severe economic losses in the poultry industry. We constructed a recombinant adenovirus that simultaneously expresses the S1 spike glycoprotein of IBV and the TM-1 protein of MG (pBH-S1-TM-1-EGFP). For comparison, we constructed two recombinant adenoviruses (pBH-S1-EGFP and pBH-TM-1-EGFP) that express either the S1 spike glycoprotein or the TM-1 protein alone. The protective efficacy of these three vaccine constructs against challenge with IBV and/or MG was evaluated in specific pathogen free chickens. Groups of seven-day-old specific pathogen free chicks were immunized twice, two weeks apart, via the oculonasal route with the pBH-S1-TM-1-EGFP, pBH-S1-EGFP, or pBH-TM-1-EGFP vaccine candidates or the commercial attenuated infectious bronchitis vaccine strain H52 and MG vaccine strain F-36 (positive controls), and challenged with virulent IBV or MG two weeks later. Interestingly, by days 7 and 14 after the booster immunization, pBH-S1-TM-1-EGFP-induced antibody titre was significantly higher (P < 0.01) compared to attenuated commercial IBV vaccine; however, there was no significant difference between the pBH-S1-TM-1-EGFP and attenuated commercial MG vaccine groups (P > 0.05). The clinical signs, the gross, and histopathological lesions scores of the adenovirus vaccine constructs were not significantly different from that of the attenuated commercial IBV or MG vaccines (positive controls) (P > 0.05). These results demonstrate the potential of the bivalent pBH-S1-TM-1-EGFP adenovirus construct as a combination vaccine against IB and mycoplasmosis.