Two +ssRNA mycoviruses cohabiting the fungal cultivar of leafcutter ants.

Leafcutter ants are dominant herbivores in the Neotropics and rely on a fungus (Leucoagaricus gongylophorus) to transform freshly gathered leaves into a source of nourishment rather than consuming the vegetation directly. Here we report two virus-like particles that were isolated from L. gongylophorus and observed using transmission electron microscopy. RNA sequencing identified two +ssRNA mycovirus strains, Leucoagaricus gongylophorus tymo-like virus 1 (LgTlV1) and Leucoagaricus gongylophorus magoulivirus 1 (LgMV1). Genome annotation of LgTlV1 (7401 nt) showed conserved domains for methyltransferase, endopeptidase, viral RNA helicase, and RNA-dependent RNA polymerase (RdRp). The smaller genome of LgMV1 (2636 nt) contains one open reading frame encoding an RdRp. While we hypothesize these mycoviruses function as symbionts in leafcutter farming systems, further study will be needed to test whether they are mutualists, commensals, or parasites.

The Virome of Cocoa Fermentation-Associated Microorganisms.

Theobroma cacao plantations are of significant economic importance worldwide, primarily for chocolate production. During the harvest and processing of cocoa beans, they are subjected to fermentation either by microorganisms present in the environment (spontaneous fermentation) or the addition of starter cultures, with different strains directly contributing distinct flavor and color characteristics to the beans. In addition to fungi and bacteria, viruses are ubiquitous and can affect the quality of the fermentation process by infecting fermenting organisms, destabilizing microbial diversity, and consequently affecting fermentation quality. Therefore, in this study, we explored publicly available metatranscriptomic libraries of cocoa bean fermentation in Limon Province, Costa Rica, looking for viruses associated with fermenting microorganisms. Libraries were derived from the same sample at different time points: 7, 20, and 68 h of fermentation, corresponding to yeast- and lactic acid bacteria-driven phases. Using a comprehensive pipeline, we identified 68 viral sequences that could be assigned to 62 new viral species and 6 known viruses distributed among at least nine families, with particular abundance of elements from the Lenarviricota phylum. Interestingly, 44 of these sequences were specifically associated with ssRNA phages (Fiersviridae) and mostly fungi-infecting viral families (Botourmiaviridae, Narnaviridae, and Mitoviridae). Of note, viruses from those families show a complex evolutionary relationship, transitioning from infecting bacteria to infecting fungi. We also identified 10 and 3 viruses classified within the Totiviridae and Nodaviridae families, respectively. The quantification of the virus-derived RNAs shows a general pattern of decline, similar to the dynamic profile of some microorganism genera during the fermentation process. Unexpectedly, we identified narnavirus-related elements that showed similarity to segmented viral species. By exploring the molecular characteristics of these viral sequences and applying Hidden Markov Models, we were capable of associating these additional segments with a specific taxon. In summary, our study elucidates the complex virome associated with the microbial consortia engaged in cocoa bean fermentation that could contribute to organism/strain selection, altering metabolite production and, consequently, affecting the sensory characteristics of cocoa beans.

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.

Suboptimal bioinformatic predictions of antimicrobial resistance from whole-genome sequences in multidrug-resistant Corynebacterium isolates.

Herein, we combined different bioinformatics tools and databases (BV-BRC, ResFinder, RAST, and KmerResistance) to perform a prediction of antimicrobial resistance (AMR) in the genomic sequences of 107 Corynebacterium striatum isolates for which trustable antimicrobial susceptibility (AST) phenotypes could be retrieved. Then, the reliabilities of the AMR predictions were evaluated by different metrics: area under the ROC curve (AUC); Major Error Rates (MERs) and Very Major Error Rates (VMERs); Matthews Correlation Coefficient (MCC); F1-Score; and Accuracy. Out of 15 genes that were reliably detected in the C. striatum isolates, only tetW yielded predictive values for tetracycline resistance that were acceptable considering Food and Drug Administration (FDA)'s criteria for quality (MER < 3.0% and VMER with a 95% C.I. ≤1.5-≤7.5); this was accompanied by a MCC score higher than 0.9 for this gene. Noteworthy, our results indicate that other commonly used metrics (AUC, F1-score, and Accuracy) may render overoptimistic evaluations of AMR-prediction reliabilities on imbalanced datasets. Accordingly, out of 10 genes tested by PCR on additional multidrug-resistant Corynebacterium spp. isolates (n = 18), the tetW gene rendered the best agreement values with AST profiles (94.11%). Overall, our results indicate that genome-based AMR prediction can still be challenging for MDR clinical isolates of emerging Corynebacterium spp.

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.

Rhabdoviral Endogenous Sequences Identified in the Leishmaniasis Vector Lutzomyia longipalpis Are Widespread in Sandflies from South America.

Sandflies are known vectors of leishmaniasis. In the Old World, sandflies are also vectors of viruses while little is known about the capacity of New World insects to transmit viruses to humans. Here, we relate the identification of RNA sequences with homology to rhabdovirus nucleocapsids (NcPs) genes, initially in the Lutzomyia longipalpis LL5 cell lineage, named NcP1.1 and NcP2. The Rhabdoviridae family never retrotranscribes its RNA genome to DNA. The sequences here described were identified in cDNA and DNA from LL-5 cells and in adult insects indicating that they are transcribed endogenous viral elements (EVEs). The presence of NcP1.1 and NcP2 in the L. longipalpis genome was confirmed in silico. In addition to showing the genomic location of NcP1.1 and NcP2, we identified another rhabdoviral insertion named NcP1.2. Analysis of small RNA molecules derived from these sequences showed that NcP1.1 and NcP1.2 present a profile consistent with elements targeted by primary piRNAs, while NcP2 was restricted to the degradation profile. The presence of NcP1.1 and NcP2 was investigated in sandfly populations from South America and the Old World. These EVEs are shared by different sandfly populations in South America while none of the Old World species studied presented the insertions.

Neofunctionalization driven by positive selection led to the retention of the loqs2 gene encoding an Aedes specific dsRNA binding protein.

BACKGROUND: Mosquito borne viruses, such as dengue, Zika, yellow fever and Chikungunya, cause millions of infections every year. These viruses are mostly transmitted by two urban-adapted mosquito species, Aedes aegypti and Aedes albopictus. Although mechanistic understanding remains largely unknown, Aedes mosquitoes may have unique adaptations that lower the impact of viral infection. Recently, we reported the identification of an Aedes specific double-stranded RNA binding protein (dsRBP), named Loqs2, that is involved in the control of infection by dengue and Zika viruses in mosquitoes. Preliminary analyses suggested that the loqs2 gene is a paralog of loquacious (loqs) and r2d2, two co-factors of the RNA interference (RNAi) pathway, a major antiviral mechanism in insects. RESULTS: Here we analyzed the origin and evolution of loqs2. Our data suggest that loqs2 originated from two independent duplications of the first double-stranded RNA binding domain of loqs that occurred before the origin of the Aedes Stegomyia subgenus, around 31 million years ago. We show that the loqs2 gene is evolving under relaxed purifying selection at a faster pace than loqs, with evidence of neofunctionalization driven by positive selection. Accordingly, we observed that Loqs2 is localized mainly in the nucleus, different from R2D2 and both isoforms of Loqs that are cytoplasmic. In contrast to r2d2 and loqs, loqs2 expression is stage- and tissue-specific, restricted mostly to reproductive tissues in adult Ae. aegypti and Ae. albopictus. Transgenic mosquitoes engineered to express loqs2 ubiquitously undergo developmental arrest at larval stages that correlates with massive dysregulation of gene expression without major effects on microRNAs or other endogenous small RNAs, classically associated with RNA interference. CONCLUSIONS: Our results uncover the peculiar origin and neofunctionalization of loqs2 driven by positive selection. This study shows an example of unique adaptations in Aedes mosquitoes that could ultimately help explain their effectiveness as virus vectors.

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.

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.

A repertoire of candidate effector proteins of the fungus Ceratocystis cacaofunesta.

The genus Ceratocystis includes many phytopathogenic fungi that affect different plant species. One of these is Ceratocystis cacaofunesta, which is pathogenic to the cocoa tree and causes Ceratocystis wilt, a lethal disease for the crop. However, little is known about how this pathogen interacts with its host. The knowledge and identification of possible genes encoding effector proteins are essential to understanding this pathosystem. The present work aimed to predict genes that code effector proteins of C. cacaofunesta from a comparative analysis of the genomes of five Ceratocystis species available in databases. We performed a new genome annotation through an in-silico analysis. We analyzed the secretome and effectorome of C. cacaofunesta using the characteristics of the peptides, such as the presence of signal peptide for secretion, absence of transmembrane domain, and richness of cysteine residues. We identified 160 candidate effector proteins in the C. cacaofunesta proteome that could be classified as cytoplasmic (102) or apoplastic (58). Of the total number of candidate effector proteins, 146 were expressed, presenting an average of 206.56 transcripts per million. Our database was created using a robust bioinformatics strategy, followed by manual curation, generating information on pathogenicity-related genes involved in plant interactions, including CAZymes, hydrolases, lyases, and oxidoreductases. Comparing proteins already characterized as effectors in Sordariomycetes species revealed five groups of protein sequences homologous to C. cacaofunesta. These data provide a valuable resource for studying the infection mechanisms of these pathogens in their hosts.

Gut microbiota from patients with COVID-19 cause alterations in mice that resemble post-COVID symptoms.

Long-term sequelae of coronavirus disease (COVID)-19 are frequent and of major concern. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection affects the host gut microbiota, which is linked to disease severity in patients with COVID-19. Here, we report that the gut microbiota of post-COVID subjects had a remarkable predominance of Enterobacteriaceae strains with an antibiotic-resistant phenotype compared to healthy controls. Additionally, short-chain fatty acid (SCFA) levels were reduced in feces. Fecal transplantation from post-COVID subjects to germ-free mice led to lung inflammation and worse outcomes during pulmonary infection by multidrug-resistant Klebsiella pneumoniae. transplanted mice also exhibited poor cognitive performance. Overall, we show prolonged impacts of SARS-CoV-2 infection on the gut microbiota that persist after subjects have cleared the virus. Together, these data demonstrate that the gut microbiota can directly contribute to post-COVID sequelae, suggesting that it may be a potential therapeutic target.

The virome of the invasive Asian bush mosquito Aedes japonicus in Europe.

The Asian bush mosquito Aedes japonicus is rapidly invading North America and Europe. Due to its potential to transmit multiple pathogenic arthropod-borne (arbo)viruses including Zika virus, West Nile virus, and chikungunya virus, it is important to understand the biology of this vector mosquito in more detail. In addition to arboviruses, mosquitoes can also carry insect-specific viruses that are receiving increasing attention due to their potential effects on host physiology and arbovirus transmission. In this study, we characterized the collection of viruses, referred to as the virome, circulating in Ae. japonicus populations in the Netherlands and France. Applying a small RNA-based metagenomic approach to Ae. japonicus, we uncovered a distinct group of viruses present in samples from both the Netherlands and France. These included one known virus, Ae. japonicus narnavirus 1 (AejapNV1), and three new virus species that we named Ae. japonicus totivirus 1 (AejapTV1), Ae. japonicus anphevirus 1 (AejapAV1) and Ae. japonicus bunyavirus 1 (AejapBV1). We also discovered sequences that were presumably derived from two additional novel viruses: Ae. japonicus bunyavirus 2 (AejapBV2) and Ae. japonicus rhabdovirus 1 (AejapRV1). All six viruses induced strong RNA interference responses, including the production of twenty-one nucleotide-sized small interfering RNAs, a signature of active replication in the host. Notably, AejapBV1 and AejapBV2 belong to different viral families; however, no RNA-dependent RNA polymerase sequence has been found for AejapBV2. Intriguingly, our small RNA-based approach identified an ∼1-kb long ambigrammatic RNA that is associated with AejapNV1 as a secondary segment but showed no similarity to any sequence in public databases. We confirmed the presence of AejapNV1 primary and secondary segments, AejapTV1, AejapAV1, and AejapBV1 by reverse transcriptase polymerase chain reaction (PCR) in wild-caught Ae. japonicus mosquitoes. AejapNV1 and AejapTV1 were found at high prevalence (87-100 per cent) in adult females, adult males, and larvae. Using a small RNA-based, sequence-independent metagenomic strategy, we uncovered a conserved and prevalent virome among Ae. japonicus mosquito populations. The high prevalence of AejapNV1 and AejapTV1 across all tested mosquito life stages suggests that these viruses are intimately associated with Ae. japonicus.

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.

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.

Genomic insights into the c-di-GMP signaling and biofilm development in the saprophytic spirochete Leptospira biflexa.

C-di-GMP is a bacterial second messenger with central role in biofilm formation. Spirochete bacteria from Leptospira genus present a wide diversity, with species of medical importance and environmental species, named as saprophytic. Leptospira form biofilms in the rat's reservoir kidneys and in the environment. Here, we performed genomic analyses to identify enzymatic and effector c-di-GMP proteins in the saprophytic biofilm-forming species Leptospira biflexa serovar Patoc. We identified 40 proteins through local alignments. Amongst them, 16 proteins are potentially functional diguanylate cyclases, phosphodiesterases, or hybrid proteins. We also identified nine effectors, including PilZ proteins. Enrichment analyses suggested that c-di-GMP interacts with cAMP signaling system, CsrA system, and flagella assembly regulation during biofilm development of L. biflexa. Finally, we identified eight proteins in the pathogen Leptospira interrogans serovar Copenhageni that share high similarity with L. biflexa c-di-GMP-related proteins. This work revealed proteins related to c-di-GMP turnover and cellular response in Leptospira and their potential roles during biofilm development.

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.

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.

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.