Phylogenomics and gene selection in Aspergillus welwitschiae: Possible implications in the pathogenicity in Agave sisalana.

Aspergillus welwitschiae causes bole rot disease in sisal (Agave sisalana and related species) which affects the production of natural fibers in Brazil, the main worldwide producer of sisal fibers. This fungus is a saprotroph with a broad host range. Previous research established A. welwitschiae as the only causative agent of bole rot in the field, but little is known about the evolution of this species and its strains. In this work, we performed a comparative genomics analysis of 40 Aspergillus strains. We show the conflicting molecular identity of this species, with one sisal-infecting strain sharing its last common ancestor with Aspergillus niger, having diverged only 833 thousand years ago. Furthermore, our analysis of positive selection reveals sites under selection in genes coding for siderophore transporters, Sodium­calcium exchangers, and Phosphatidylethanolamine-binding proteins (PEBPs). Herein, we discuss the possible impacts of these gene functions on the pathogenicity in sisal.

State of the Art of the Molecular Biology of the Interaction between Cocoa and Witches' Broom Disease: A Systematic Review.

Significant scientific advances to elucidate the Moniliophthora perniciosa pathosystem have been achieved in recent years, but the molecular biology of this pathogen-host interaction is still a field with many unanswered questions. In order to present insights at the molecular level, we present the first systematic review on the theme. All told, 1118 studies were extracted from public databases. Of these, 109 were eligible for the review, based on the inclusion and exclusion criteria. The results indicated that understanding the transition from the biotrophic-necrotrophic phase of the fungus is crucial for control of the disease. Proteins with strong biotechnological potential or that can be targets for pathosystem intervention were identified, but studies regarding possible applications are still limited. The studies identified revealed important genes in the M. perniciosa-host interaction and efficient molecular markers in the search for genetic variability and sources of resistance, with Theobroma cacao being the most common host. An arsenal of effectors already identified and not explored in the pathosystem were highlighted. This systematic review contributes to the understanding of the pathosystem at the molecular level, offering new insights and proposing different paths for the development of new strategies to control witches' broom disease.

BASIDIN as a New Protein Effector of the Phytopathogen Causing Witche's Broom Disease in Cocoa.

The fungus Moniliophthora perniciosa secretes protein effectors that manipulate the physiology of the host plant, but few effectors of this fungus have had their functions confirmed. We performed functional characterization of a promising candidate effector of M. perniciosa. The inoculation of rBASIDIN at 4 µmol L-1 in the mesophyll of leaflets of Solanum lycopersicum caused symptoms of shriveling within 6 h without the presence of necrosis. However, when sprayed on the plant at a concentration of 11 µmol L-1, it caused wilting symptoms only 2 h after application, followed by necrosis and cell death at 48 h. rBASIDIN applied to Theobroma cacao leaves at the same concentration caused milder symptoms. rBASIDIN caused hydrogen peroxide production in leaf tissue, damaging the leaf membrane and negatively affecting the photosynthetic rate of Solanum lycopersicum plants. Phylogenetic analysis indicated that BASIDIN has orthologs in other phytopathogenic basidiomycetes. Analysis of the transcripts revealed that BASIDIN and its orthologs are expressed in different fungal species, suggesting that this protein is differentially regulated in these basidiomycetes. Therefore, the results of applying BASIDIN allow the inference that it is an effector of the fungus M. perniciosa, with a strong potential to interfere in the defense system of the host plant.

A single unidirectional piRNA cluster similar to the flamenco locus is the major source of EVE-derived transcription and small RNAs in Aedes aegypti mosquitoes.

Endogenous viral elements (EVEs) are found in many eukaryotic genomes. Despite considerable knowledge about genomic elements such as transposons (TEs) and retroviruses, we still lack information about nonretroviral EVEs. Aedes aegypti mosquitoes have a highly repetitive genome that is covered with EVEs. Here, we identified 129 nonretroviral EVEs in the AaegL5 version of the A. aegypti genome. These EVEs were significantly associated with TEs and preferentially located in repeat-rich clusters within intergenic regions. Genome-wide transcriptome analysis showed that most EVEs generated transcripts although only around 1.4% were sense RNAs. The majority of EVE transcription was antisense and correlated with the generation of EVE-derived small RNAs. A single genomic cluster of EVEs located in a 143 kb repetitive region in chromosome 2 contributed with 42% of antisense transcription and 45% of small RNAs derived from viral elements. This region was enriched for TE-EVE hybrids organized in the same coding strand. These generated a single long antisense transcript that correlated with the generation of phased primary PIWI-interacting RNAs (piRNAs). The putative promoter of this region had a conserved binding site for the transcription factor Cubitus interruptus, a key regulator of the flamenco locus in Drosophila melanogaster Here, we have identified a single unidirectional piRNA cluster in the A. aegypti genome that is the major source of EVE transcription fueling the generation of antisense small RNAs in mosquitoes. We propose that this region is a flamenco-like locus in A. aegypti due to its relatedness to the major unidirectional piRNA cluster in Drosophila melanogaster.

Culex quinquefasciatus mosquitoes do not support replication of Zika virus.

The rapid spread of Zika virus (ZIKV) in the Americas raised many questions about the role of Culex quinquefasciatus mosquitoes in transmission, in addition to the key role played by the vector Aedes aegypti. Here we analysed the competence of Cx. quinquefasciatus (with or without Wolbachia endosymbionts) for a ZIKV isolate. We also examined the induction of RNA interference pathways after viral challenge and the production of small virus-derived RNAs. We did not observe any infection nor such small virus-derived RNAs, regardless of the presence or absence of Wolbachia. Thus, Cx. quinquefasciatus does not support ZIKV replication and Wolbachia is not involved in producing this phenotype. In short, these mosquitoes are very unlikely to play a role in transmission of ZIKV.

Sequence-independent characterization of viruses based on the pattern of viral small RNAs produced by the host.

Virus surveillance in vector insects is potentially of great benefit to public health. Large-scale sequencing of small and long RNAs has previously been used to detect viruses, but without any formal comparison of different strategies. Furthermore, the identification of viral sequences largely depends on similarity searches against reference databases. Here, we developed a sequence-independent strategy based on virus-derived small RNAs produced by the host response, such as the RNA interference pathway. In insects, we compared sequences of small and long RNAs, demonstrating that viral sequences are enriched in the small RNA fraction. We also noted that the small RNA size profile is a unique signature for each virus and can be used to identify novel viral sequences without known relatives in reference databases. Using this strategy, we characterized six novel viruses in the viromes of laboratory fruit flies and wild populations of two insect vectors: mosquitoes and sandflies. We also show that the small RNA profile could be used to infer viral tropism for ovaries among other aspects of virus biology. Additionally, our results suggest that virus detection utilizing small RNAs can also be applied to vertebrates, although not as efficiently as to plants and insects.

Functional specialization of the small interfering RNA pathway in response to virus infection.

In Drosophila, post-transcriptional gene silencing occurs when exogenous or endogenous double stranded RNA (dsRNA) is processed into small interfering RNAs (siRNAs) by Dicer-2 (Dcr-2) in association with a dsRNA-binding protein (dsRBP) cofactor called Loquacious (Loqs-PD). siRNAs are then loaded onto Argonaute-2 (Ago2) by the action of Dcr-2 with another dsRBP cofactor called R2D2. Loaded Ago2 executes the destruction of target RNAs that have sequence complementarity to siRNAs. Although Dcr-2, R2D2, and Ago2 are essential for innate antiviral defense, the mechanism of virus-derived siRNA (vsiRNA) biogenesis and viral target inhibition remains unclear. Here, we characterize the response mechanism mediated by siRNAs against two different RNA viruses that infect Drosophila. In both cases, we show that vsiRNAs are generated by Dcr-2 processing of dsRNA formed during viral genome replication and, to a lesser extent, viral transcription. These vsiRNAs seem to preferentially target viral polyadenylated RNA to inhibit viral replication. Loqs-PD is completely dispensable for silencing of the viruses, in contrast to its role in silencing endogenous targets. Biogenesis of vsiRNAs is independent of both Loqs-PD and R2D2. R2D2, however, is required for sorting and loading of vsiRNAs onto Ago2 and inhibition of viral RNA expression. Direct injection of viral RNA into Drosophila results in replication that is also independent of Loqs-PD. This suggests that triggering of the antiviral pathway is not related to viral mode of entry but recognition of intrinsic features of virus RNA. Our results indicate the existence of a vsiRNA pathway that is separate from the endogenous siRNA pathway and is specifically triggered by virus RNA. We speculate that this unique framework might be necessary for a prompt and efficient antiviral response.

Circulation of bee-infecting viruses in Brazil: a call for action.

Bees are fundamental for maintaining pollination-dependent plant populations, both economically and ecologically. In Brazil, they constitute 66.3% of pollinators, contributing to an annual market value estimated at R$ 43 billion for pollination services. Unfortunately, worldwide bee populations are declining due to parasites and pathogens, more specifically viruses, alongside climate change, habitat loss, and pesticides. In this scenario, extensive research concerning bee diversity, virus diversity and surveillance, is necessary to aid the conservation of native managed pollinators and potential wild alternatives besides mitigating the emergence and spread of viral pathogens. A decrease in pollination can be a point of economic vulnerability in a country like Brazil because of its main dependence on food exports. Here we conducted a study aiming to obtain an overview of circulating viruses in bees within Brazilian territory highlighting the need for further studies to have a more realistic view of bee-infecting viruses in Brazil.

Phylodynamics of avian influenza A(H5N1) viruses from outbreaks in Brazil.

Our study identified strains of the A/H5N1 virus in analyzed samples of subsistence poultry, wild birds, and mammals, belonging to clade 2.3.4.4b, genotype B3.2, with very high genetic similarity to strains from Chile, Uruguay, and Argentina. This suggests a migratory route for wild birds across the Pacific, explaining the phylogenetic relatedness. The Brazilian samples displayed similarity to strains that had already been previously detected in South America. Phylogeographic analysis suggests transmission of US viruses from Europe and Asia, co-circulating with other lineages in the American continent. As mutations can influence virulence and host specificity, genomic surveillance is essential to detect those changes, especially in critical regions, such as hot spots in the HA, NA, and PB2 sequences. Mutations in the PB2 gene (D701N and Q591K) associated with adaptation and transmission in mammals were detected suggesting a potential zoonotic risk. Nonetheless, resistance to neuraminidase inhibitors (NAIs) was not identified, however, continued surveillance is crucial to detect potential resistance. Our study also mapped the spread of the virus in the Southern hemisphere, identifying possible entry routes and highlighting the importance of surveillance to prevent outbreaks and protect both human and animal populations.

The Viromes of Six Ecosystem Service Provider Parasitoid Wasps.

Parasitoid wasps are fundamental insects for the biological control of agricultural pests. Despite the importance of wasps as natural enemies for more sustainable and healthy agriculture, the factors that could impact their species richness, abundance, and fitness, such as viral diseases, remain almost unexplored. Parasitoid wasps have been studied with regard to the endogenization of viral elements and the transmission of endogenous viral proteins that facilitate parasitism. However, circulating viruses are poorly characterized. Here, RNA viromes of six parasitoid wasp species are studied using public libraries of next-generation sequencing through an integrative bioinformatics pipeline. Our analyses led to the identification of 18 viruses classified into 10 families (Iflaviridae, Endornaviridae, Mitoviridae, Partitiviridae, Virgaviridae, Rhabdoviridae, Chuviridae, Orthomyxoviridae, Xinmoviridae, and Narnaviridae) and into the Bunyavirales order. Of these, 16 elements were described for the first time. We also found a known virus previously identified on a wasp prey which suggests viral transmission between the insects. Altogether, our results highlight the importance of virus surveillance in wasps as its service disruption can affect ecology, agriculture and pest management, impacting the economy and threatening human food security.

Exploring the Mycovirus Universe: Identification, Diversity, and Biotechnological Applications.

Viruses that infect fungi are known as mycoviruses and are characterized by the lack of an extracellular phase. In recent years, the advances on nucleic acids sequencing technologies have led to a considerable increase in the number of fungi-infecting viral species described in the literature, with a special interest in assessing potential applications as fungal biocontrol agents. In the present study, we performed a comprehensive review using Scopus, Web of Science, and PubMed databases to mine mycoviruses data to explore their molecular features and their use in biotechnology. Our results showed the existence of 267 mycovirus species, of which 189 are recognized by the International Committee on Taxonomy of Viruses (ICTV). The majority of the mycoviruses identified have a dsRNA genome (38.6%), whereas the Botourmiaviridae (ssRNA+) alone represents 14% of all mycoviruses diversity. Regarding fungal hosts, members from the Sclerotinicaeae appeared as the most common species described to be infected by mycoviruses, with 16 different viral families identified so far. It is noteworthy that such results are directly associated with the high number of studies and strategies used to investigate the presence of viruses in members of the Sclerotinicaeae family. The knowledge about replication strategy and possible impact on fungi biology is available for only a small fraction of the mycoviruses studied, which is the main limitation for considering these elements potential targets for biotechnological applications. Altogether, our investigation allowed us to summarize the general characteristics of mycoviruses and their hosts, the consequences, and the implications of this knowledge on mycovirus-fungi interactions, providing an important source of information for future studies.

Uncovering a Complex Virome Associated with the Cacao Pathogens Ceratocystis cacaofunesta and Ceratocystis fimbriata.

Theobroma cacao is one of the main crops of economic importance in the world as the source of raw material for producing chocolate and derivatives. The crop is the main source of income for thousands of small farmers, who produce more than 80% of the world's cocoa supply. However, the emergence, re-emergence and proliferation of pathogens, such as Ceratocystis spp., the causative agent of Ceratocystis wilt disease and canker disease, have been affecting the sustainability of many crops. Fungal control is laborious, often depending on fungicides that are expensive and/or toxic to humans, prompting researchers to look for new solutions to counteract the proliferation of these pathogens, including the use of biological agents such as mycoviruses. In this study, we investigated the diversity of microorganisms associated with the T. cacao pathogens Ceratocystis cacaofunesta and Ceratocystis fimbriata with a focus on the virome using RNA sequencing data available in public databases. We used a comprehensive bioinformatics pipeline containing several steps for viral sequence enrichment and took advantage of an integrated assembly step composed of different assemblers followed by sequence similarity searches using NCBI nonredundant databases. Our strategy was able to identify four putative C. cacaofunesta viruses (hypovirus, sclerotimonavirus, alphapartitivirus and narnavirus) and six C. fimbriata viruses (three alphaendornaviruses, one victorivirus and two mitoviruses). All the viral sequences identified showed similarity to viral genomes in public databases only at the amino acid level, likely representing new viral species. Of note, we present the first report of viruses associated with the cacao pathogens C. cacaofunesta and C. fimbriata and the second report of viral species infecting members of the Ceratocystidaceae family. Our findings highlight the need for further prospective studies to uncover the real diversity of fungus-infecting viruses that can contribute to the development of new management strategies.

Mining Public Data to Investigate the Virome of Neglected Pollinators and Other Floral Visitors.

This study reports the virome investigation of pollinator species and other floral visitors associated with plants from the south of Bahia: Aphis aurantii, Atrichopogon sp., Dasyhelea sp., Forcipomyia taiwana, and Trigona ventralis hoozana. Studying viruses in insects associated with economically important crops is vital to understand transmission dynamics and manage viral diseases that pose as threats for global food security. Using literature mining and public RNA next-generation sequencing data deposited in the NCBI SRA database, we identified potential vectors associated with Malvaceae plant species and characterized the microbial communities resident in these insects. Bacteria and Eukarya dominated the metagenomic analyses of all taxon groups. We also found sequences showing similarity to elements from several viral families, including Bunyavirales, Chuviridae, Iflaviridae, Narnaviridae, Orthomyxoviridae, Rhabdoviridae, Totiviridae, and Xinmoviridae. Phylogenetic analyses indicated the existence of at least 16 new viruses distributed among A. aurantii (3), Atrichopogon sp. (4), Dasyhelea sp. (3), and F. taiwana (6). No novel viruses were found for T. ventralis hoozana. For F. taiwana, the available libraries also allowed us to suggest possible vertical transmission, while for A. aurantii we followed the infection profile along the insect development. Our results highlight the importance of studying the virome of insect species associated with crop pollination, as they may play a crucial role in the transmission of viruses to economically important plants, such as those of the genus Theobroma, or they will reduce the pollination process. This information may be valuable in developing strategies to mitigate the spread of viruses and protect the global industry.

Corynebacterium guaraldiae sp. nov.: a new species of Corynebacterium from human infections.

Non-diphtheria Corynebacterium species (NDC) belonging to the human skin and mucosa microbiota are frequently neglected as contaminants. However, reports of human infections by Corynebacterium spp. have increased considerably in recent years. In this study, a group of six NDC isolates of urine (n = 5) and sebaceous cyst (n = 1) from two South American countries were identified at genus level or misidentified based on API® Coryne and genetic/molecular analyses. The 16S rRNA (99.09-99.56%) and rpoB (96.18-97.14%) gene sequence similarities of the isolates were higher when compared with Corynebacterium aurimucosum DSM 44532 T. Multilocus sequence analysis (MLSA) indicated that these six NDC isolates compose a distinctive phylogenetic clade. Genome-based taxonomic analysis with the whole-genome sequences was able to separate these six isolates from other known Corynebacterium type strains. Average nucleotide identity (ANI), average amino acid identity (AAI), and digital DNA-DNA hybridization (dDDH) values between closely related type strains and the six isolates were considerably lower than the currently recommended threshold values for species circumscription. Phylogenetic and genomic taxonomy analyses indicated these microorganisms as a novel Corynebacterium species, for which we formally propose the name Corynebacterium guaraldiae sp. nov. with isolate 13T (= CBAS 827T = CCBH 35012T) as type strain.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World.

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Characterization of the microbiota dynamics associated with Moniliophthora roreri, causal agent of cocoa frosty pod rot disease, reveals new viral species.

Introduction: Theobroma cacao, the cocoa tree, is a target for pathogens, such as fungi from the genera Phytophthora, Moniliophthora, Colletotrichum, Ceratocystis, among others. Some cacao pathogens are restricted to specific regions of the world, such as the Cacao swollen shoot virus (CSSV) in West African countries, while others are expanding geographically, such as Moniliophthora roreri in the Americas. M. roreri is one of the most threatening cacao pathogens since it directly attacks the cacao pods driving a significant reduction in production, and therefore economic losses. Despite its importance, the knowledge about the microenvironment of this pathogen and the cocoa pods is still poorly characterized. Methods: Herein we performed RNA sequencing of spores in differential stages of culture in a medium supplemented with cacao pod extract and mycelium collected of the susceptible variety ICT 7121 naturally infected by the pathogen to evaluate the diversity and transcriptional activity of microorganisms associated with the in vitro sporulation of M. roreri. Results: Our data revealed a great variety of fungi and bacteria associated with M. roreri, with an exceptional diversity of individuals from the genus Trichoderma sp. Interestingly, the dynamics of microorganisms from different kingdoms varied proportionally, suggesting they are somehow affected by M. roreri culture time. We also identified three sequences similar to viral genomes from the Narnaviridae family, posteriorly confirmed by phylogenetic analysis as members of the genus Narnavirus. Screening of M. roreri public datasets indicated the virus sequences circulating in samples from Ecuador, suggesting a wide spread of these elements. Of note, we did not identify traces of the viral sequences in the M. roreri genome or DNA sequencing, restricting the possibility of these sequences representing endogenized elements. Discussion: To the best of our knowledge, this is the first report of viruses infecting the fungus of the genus Moniliophthora and only the third description of viruses that are able to parasite elements from the Marasmiaceae family.

Expanding the environmental virome: Infection profile in a native rainforest tree species.

Agroforestry systems (AFS) for cocoa production combine traditional land-use practices with local biodiversity conservation, resulting in both ecological and agricultural benefits. The cacao-cabruca AFS model is widely implemented in regions of the Brazilian Atlantic Forest. Carpotroche brasiliensis (Raddi) A. Gray (Achariaceae) is a tree found in cabruca landscapes that is often used for reforestation and biotechnological applications. Despite its importance, we still lack information about viruses circulating in C. brasiliensis, particularly considering the possibility of spillover that could affect cocoa production. In our study, we analyzed the Carpotroche brasiliensis virome from Atlantic Forest and cacao-cabruca AFS regions using metatranscriptomics from several vegetative and reproductive organs. Our results revealed a diverse virome detecting near-complete or partial coding sequences of single- and double-stranded DNA and RNA viruses classified into at least six families (Botourmiaviridae, Bromoviridae, Caulimoviridae, Genomoviridae, Mitoviridae, and Rhabdoviridae) plus unclassified elements. We described with high confidence the near-complete and the partial genomes of two tentative novel viruses: Carpotroche-associated ilarvirus and Carpotroche-associated genomovirus, respectively. Interestingly, we also described sequences likely derived from a rhabdovirus, which could represent a novel member of the genus Gammanucleorhabdovirus. We observed higher viral diversity in cacao-cabruca AFS and reproductive organs of C. brasiliensis with preferential tropism to fruits, which could directly affect production. Altogether, our results provide data to better understand the virome in this unexplored agroecological interface, such as cacao-cabruca AFS and forest ecosystem, providing information on the aspects of virus-plant interactions.

An Integrative View of the Phyllosphere Mycobiome of Native Rubber Trees in the Brazilian Amazon.

The rubber tree, Hevea brasiliensis, is a neotropical Amazonian species. Despite its high economic value and fungi associated with native individuals, in its original area in Brazil, it has been scarcely investigated and only using culture-dependent methods. Herein, we integrated in silico approaches with novel field/experimental approaches and a case study of shotgun metagenomics and small RNA metatranscriptomics of an adult individual. Scientific literature, host fungus, and DNA databases are biased to fungal taxa, and are mainly related to rubber tree diseases and in non-native ecosystems. Metabarcoding retrieved specific phyllospheric core fungal communities of all individuals, adults, plantlets, and leaves of the same plant, unravelling hierarchical structured core mycobiomes. Basidiomycotan yeast-like fungi that display the potential to produce antifungal compounds and a complex of non-invasive ectophytic parasites (Sooty Blotch and Flyspeck fungi) co-occurred in all samples, encompassing the strictest core mycobiome. The case study of the same adult tree (previously studied using culture-dependent approach) analyzed by amplicon, shotgun metagenomics, and small RNA transcriptomics revealed a high relative abundance of insect parasite-pathogens, anaerobic fungi and a high expression of Trichoderma (a fungal genus long reported as dominant in healthy wild rubber trees), respectively. Altogether, our study unravels new and intriguing information/hypotheses of the foliar mycobiome of native H. brasiliensis, which may also occur in other native Amazonian trees.

The Sisal Virome: Uncovering the Viral Diversity of Agave Varieties Reveals New and Organ-Specific Viruses.

Sisal is a common name for different plant varieties in the genus Agave (especially Agave sisalana) used for high-quality natural leaf fiber extraction. Despite the economic value of these plants, we still lack information about the diversity of viruses (virome) in non-tequilana species from the genus Agave. In this work, by associating RNA and DNA deep sequencing we were able to identify 25 putative viral species infecting A. sisalana, A. fourcroydes, and Agave hybrid 11648, including one strain of Cowpea Mild Mottle Virus (CPMMV) and 24 elements likely representing new viruses. Phylogenetic analysis indicated they belong to at least six viral families: Alphaflexiviridae, Betaflexiviridae, Botourmiaviridae, Closteroviridae, Partitiviridae, Virgaviridae, and three distinct unclassified groups. We observed higher viral taxa richness in roots when compared to leaves and stems. Furthermore, leaves and stems are very similar diversity-wise, with a lower number of taxa and dominance of a single viral species. Finally, approximately 50% of the identified viruses were found in all Agave organs investigated, which suggests that they likely produce a systemic infection. This is the first metatranscriptomics study focused on viral identification in species from the genus Agave. Despite having analyzed symptomless individuals, we identified several viruses supposedly infecting Agave species, including organ-specific and systemic species. Surprisingly, some of these putative viruses are probably infecting microorganisms composing the plant microbiota. Altogether, our results reinforce the importance of unbiased strategies for the identification and monitoring of viruses in plant species, including those with asymptomatic phenotypes.