Comparison of the RNA Content of Extracellular Vesicles Derived from Paracoccidioides brasiliensis and Paracoccidioides lutzii

Paracoccidioides brasiliensis and P. lutzii cause human paracoccidioidomycosis. We have previously characterized the <200-nt RNA sub-populations contained in fungal extracellular vesicles (EVs) from P. brasiliensis Pb18 and other pathogenic fungi. We have presently used the RNA-seq strategy to compare the <200- and >200-nt RNA fractions contained in EVs isolated from culture supernatants of P. brasiliensis Pb18, Pb3, and P. lutzii Pb01. Shared mRNA sequences were related to protein modification, translation, and DNA metabolism/biogenesis, while those related to transport and oxidation-reduction were exclusive to Pb01. The presence of functional full-length mRNAs was validated by in vitro translation. Among small non-coding (nc)RNA, 15 were common to all samples; small nucleolar (sno)RNAs were enriched in P. brasiliensis EVs, whereas for P. lutzii there were similar proportions of snoRNA, rRNA, and tRNA. Putative exonic sRNAs were highly abundant in Pb18 EVs. We also found sRNA sequences bearing incomplete microRNA structures mapping to exons. RNA-seq data suggest that extracellular fractions containing Pb18 EVs can modulate the transcriptome of murine monocyte-derived dendritic cells in a transwell system. Considering that sRNA classes are involved in transcription/translation modulation, our general results may indicate that differences in virulence among fungal isolates can be related to their distinct EV-RNA content.


Introduction
Paracoccidioides brasiliensis and P. lutzii cause human paracoccidioidomycosis (PCM), which is an endemic systemic mycosis prevalent in Latin American countries, but predominantly reported in Brazil [1]. The successful establishment of this fungal infection depends on the thermo-dependent dimorphic transition from environmental conidia to the yeast phase in the lung alveoli. P. lutzii isolates are responsible for the majority of PCM cases in the Central and Northern Brazil and are represented by the highly studied Pb01 isolate [2]. P. brasiliensis is a complex of several phylogenetic groups originating from Brazilian regions and other countries [3,4]. The Pb18 and Pb3 isolates represent the phylogenetic main species S1 and a cryptic PS2 group, respectively [5]. PS2 isolates cause a milder form of murine PCM when compared to S1 isolates, apparently by stimulating a predominant Th1-type of cellular immune response [6], but the fungal features that contribute to their distinct virulence are not known. A broad analysis of genome diversity in Paracoccidioides showed that the phenotypic differences among phylogenetic groups and species are not due to differences in large protein families, but probably due to unique genes with no orthologs in other lineages [7].
Extracellular vesicles (EVs) constitute an important cellular mechanism of non-conventional secretion across all kingdoms. EV is a general term used to define spherical bi-layered-membrane structures that group, according with their origin, (a) exosomes of 30-100 nm deriving from multivesicular bodies; (b) microvesicles or ectosomes of 100-1000 nm originating from either budding or invagination of the plasma membrane; (c) apoptotic bodies that are larger than 1000 nm [8][9][10][11][12]. EVs can safely transport to the extracellular environment and distant sites a vast number of proteins, including active enzymes and toxins, lipoproteins, DNA, RNA, polysaccharides, and pigments [9]. EVs are involved in a broad array of intercellular communication and active biomolecule transfer that have implications in physiological cellular processes, infection, and disease [13,14]. EVs from microorganisms can directly interact with cells of the immune system and affect the infection process [15]. Distinct RNA specimens can be transported inside EV-like structures or associated with RNA-binding proteins and high-density lipoprotein complexes; therefore, the standardization of EV isolation and RNA purification is essential step to evaluate those populations [12,16]. After the detection of mRNA and small RNAs (sRNA) within EVs [17], the role of micro RNA (miRNA) from EVs in the genetic regulation of recipient cells has been widely addressed [18,19]. EVs also contain different classes of long non-coding RNA (lncRNAs) that are potentially able to regulate the transcription by recruiting epigenetic modifiers in the recipient cells [20].
We have previously characterized the sRNA fraction contained in EVs isolated from P. brasiliensis (isolate Pb18), Cryptococcus neoformans, Candida albicans, and also from Saccharomyces cerevisiae [21]. The sRNA fraction included sequences of various sizes shorter than 250 nt, specifically, 114 small non-coding (nc)RNA sequences of the small nucleolar (sno)RNA and transporter (t)RNA classes. We also identified 1246 conserved miRNA-like sequences, from which 20 were common to all samples and 47 had differentially represented levels. There was a small percentage of mRNA (10%) that co-fractionated in the <200 nt-enriched fraction and was also characterized. Interestingly, these sequences were unique to EVs from each species [21].
The present work aimed at comparing the RNA populations carried by EVs from two isolates (Pb18 and Pb3) representative of P. brasiliensis lineages S1 and PS2 and P. lutzii (Pb01), since differences in RNA content could be valuable to partially explain differences in virulence among isolates. We also tested the functionality of EVs mRNA and if extracellular fractions containing fungal EVs can evoke transcriptional changes in dendritic cells.

EVs Isolation
Paracoccidioides EV preparations were performed as described [21,22]. In brief, culture supernatants from two 500-mL cultures were centrifuged for 15 min at 4000× g, then for 30 min at 15,000× g, followed by filtration through a 45 µm membrane and concentration using Amicon ultrafiltration membranes (100-kDa cutoff). The cell-free, debris-free concentrated supernatant was ultracentrifuged at 100,000× g for 1 h at 4 • C and the pellets containing EVs were washed in phosphate-buffered saline (PBS). This protocol for EV isolation results in EVs with size peaks between 40 to 80 nm (suggestive of exosomes), as verified by nanoparticle-tracking analysis (NTA) for Pb18; small peaks over 150 nm are also visible. The final pellets were lyophilized for RNA extraction.

Indirect Co-Culture of Murine Monocyte-Derived CD11c+ Cells (MoDC) and Pb18
Indirect co-culture of MoDC with Pb18 was performed using a transwell system (0.4-µm membrane porosity) in 1 mL DMEM supplemented with 10% FBS in each chamber. Pb18 yeast cells (1 × 10 6 ) were seeded in the upper chamber and 1 × 10 6 CD11c+ cells (dendritic cells, MoDC) were placed in the lower chamber of a six-well plate. The plates were incubated for 48 h at 37 • C, at 5% CO 2 , the MoDCs were harvested, washed and used for RNA extraction.

RNA Isolation and Sequencing
Total RNA extraction from Paracoccidioides EVs and fractionation were carried out as described [21]. Small and large RNA were fractionated using the miRNeasy mini kit (Qiagen, Germantown, MD, USA) and the RNeasy MinElute Cleanup Kit (Qiagen), according to the manufacturer's protocol. We used an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) to test the integrity of the RNA preparations before sequencing. Total RNA was extracted from 1 × 10 6 MoDC cells using the RNeasy kit (Qiagen). The RNA was eluted in RNAse-free water and stored at −70 • C. The RNA-seq was carried out as previously described [21].

Data Analysis
The RNA sequences were analyzed by CLC Genomics Workbench© v 7.0 (Qiagen), using both the corresponding Paracoccidioides genomes from NCBI as reference (Pb18-ABKI02000000, Pb3-ABHV02000000, and Pb01-ABKH02000000) and the Saccharomyces cerevisiae non-coding database (https://www.yeastgenome.org/). All the sequencing analyses were performed in triplicate, except for the < 200nt enriched-fraction from Pb18 EV that had a duplicate evaluated. The parameters used for the alignments were: mismatch cost (2), insertion cost (3), deletion cost (3), length fraction (0.8), and similarity fraction (0.8). The abundance values are presented in reads per kilobase of exon model per million mapped reads (RPKM) [23]. For the analysis of secondary structure, the RNA-seq reads obtained for the < 200 nt fraction were directly evaluated in a PPFold plugin in the CLC Genomics Workbench v.10.1, using default parameters [24]. The sequences were up to 50-nt long.

Data Access
The RNA-seq data have been deposited at the Sequence Read Archive (SRA) database under the accession number (SRA: SRP022849).

In Vitro Translation and Protein Analysis
The in vitro translation assay was performed according to the manufacturer's instructions. Briefly, 0.6 µg total RNA from Pb18 EVs were incubated for 90 min at 30 • C with 35 µL Rabbit Reticulocyte Lysate (Promega) containing leucine/methionine (0.5 pM) and 40 U RNAseOUT (Invitrogen). The reaction product was filtered through a Millipore ® membrane (0.45 µm), diluted 5 times in PBS-NaCl 0.3 M and incubated at 4 • C with 200 µL Ni-NTA (Amersham) until the resin became red due to hemoglobin binding [25]. The suspension was centrifuged (5 min at 16,000× g at 4 • C) and the proteins were precipitated overnight (−20 • C) from the supernatant by addition of 3 volumes of cold acetone. The sample was spun down for 30 min at 16,000× g (4 • C) and washed with cold acetone. The precipitated proteins were denatured, reduced, trypsinized, and desalted as described [25]. The tryptic peptides were suspended in 0.1% formic acid (FA) and loaded to an LTQ-VelosOrbitrap (Thermo Fisher Scientific, Waltham, MA, USA) through a coupled nanoHPLC (Proxeon, Odense, Denmark). The samples were desalted and concentrated in a pre-column (10-µm C18 beads, Phenomenex, 100 µm × 2 cm, Torrance, CA, USA) and separated at 200 nL/min in a reverse-phase capillary column (5-µm beads, Phenomenex 10 cm × 75 mm). The peptides were eluted by a linear gradient from 5% to 40% of solvent A [solvent A: 5% acetonitrile (ACN)/0.1% FA; solvent B: 100% ACN/0.1% FA] over 90 min and for an additional 15 min with up to 95% of solvent B. The eluted peptides were directly injected in the mass spectrometer via a nanoelectrospray set at 2.2 kV. All analyses were performed in the positive ionization mode at the 50-2000 m/z range. The mass spectrometer was operated in the data-dependent acquisition mode to automatically switch between one orbitrap full-scan and ten ion trap tandem mass spectra. The *.raw data files were processed at MaxQuant 1.3.0.5 and the searches performed at Andromeda against a merged database of P. brasiliensis Pb18 (https: //www.ncbi.nlm.nih.gov/genome/?term=txid502780[Organism:noexp]) and Oryctolagus cuniculus of TrEMBL (downloaded at 2015/02) using a 1% false discovery rate (FDR). The search parameters included: (i) carbamidomethylation of cysteine residues as a fixed modification; (ii) oxidation of methionine residues as a variable modification; and (iii) 6 ppm and 0.5 Da for MS1 and MS2 tolerance, respectively. Proteins identified as contaminants at a reverse database and "only identified by site" were excluded from further analysis.

Results
We characterized the ncRNA and mRNA sub-populations contained in EVs isolated from culture supernatant samples of the representative Paracoccidioides isolates Pb18, Pb3, and Pb01. Our results compile the analysis of three independent biological replicates for each fungal sample. In addition, we evaluated the effect of Pb18 extracellular fractions containing EVs on the transcription pattern of a murine monocyte-derived (CD11c+) cell line co-cultivated with Pb18.

Paracoccidioides EVs Carry Functional mRNA
In a previous study, we characterized mRNA sequences in fungal EVs that were co-purified in the < 200-nt RNA fraction [21]. We presently analyzed large RNA fractions to identify the most abundant mRNA sequences contained in Paracoccidioides EVs. The reads obtained from the mRNA libraries (>200 nt) were aligned with each isolate-specific genome at NCBI (Pb18-ABKI02000000, Pb3-ABHV02000000, and Pb01-ABKH02000000). For data validation, we only considered sequences with expression values of RPKM > 100 in all biological replicates. Those sequences were individually accessed and only mRNA transcripts with reads covering at least 50% of the CDS were validated (Table 1). According to the aforementioned parameters, we validated a total of 30 mRNA sequences in EV samples from Pb18 (14), Pb3 (5), and Pb01 (18), as seen in Table 1 and Figure 1A. The sequences coding for polyubiquitin and histone H3 were found in EVs from all isolates, while those coding for calmodulin, elongation factor 1-alpha, and histone H2a were detected in two isolates ( Table 1). The remaining sequences were exclusive to either Pb18 or Pb01. For Pb3, there was only one exclusive mRNA coding for a pre-mRNA splicing factor. In general, the mRNA sequences were mainly related to protein modification, translation, and DNA metabolism/biogenesis ( Table 1). The latter function grouped most of the overlapping transcripts among the three samples. On the other hand, sequences related to transport and oxidation-reduction were only detected in Pb01. Twenty-seven out of the 30 validated sequences had over 69% of CDS coverage. Considering that we only validated mRNAs with over 50% coverage of each CDS, comparisons with our previous work were unfruitful because the validation parameters were distinct [21].    To confirm the presence of functional full-length mRNA in Paracoccidioides EVs, we carried out an in vitro translation assay. The Pb18 mRNA EV sample was translated using a rabbit reticulocyte lysate. Thereafter, the contaminating hemoglobin was removed by incubation with NI-NTA agarose and the remaining translated proteins were analyzed by LC-MS/MS. We found five proteins that were translated from EV mRNA ( Table 2). Among them, only the transcript for a heat shock protein (PADG_07715) was included in Table 1 (RPKM > 100). Together, our results show that the Paracoccidiodes EVs contain functional full-length mRNAs and that they are differentially represented in Pb18, Pb3, and Pb01. The EV RNA reads detected in the < 200-nt fraction were initially aligned with each Paracoccidioides isolate-specific genome at NCBI. We noticed that this fraction was composed mostly of short 25-nt sequences in average that aligned to a specific region of a particular mRNA exon. In our previous work we observed that over 80% of the reads found in the <200-nt RNA fraction from EVs mapped to exons not only in P. brasiliensis Pb18, but also in C. albicans and S. cerevisiae, however not in C. neoformans [21]. In fungi, exonic short interfering RNAs (ex-siRNA) have been described to match to unique exon sites in either the reverse or the forward directions. They originate directly from single-strand RNA via an RNA-dependent RNA polymerase to generate double-stranded RNA (dsRNA) and are then converted to interfering RNA by Dicer [26]. Therefore, we presently deepened our investigation of the reads that had high depth of coverage in exonic regions. Each high-coverage coding sequence was manually accessed to evaluate the alignment position of the reads. Considering the reads generated by the sequencing of both < and > 200-nt fractions, 160 mRNA sequences showed high values of depth of coverage. Among them, only 30 had reads covering more than 50% of the entire mRNA, as detailed earlier. The remaining 130 were represented by small sequences that aligned at unique sites. We therefore called them exonic sRNA. We found a total of 104 (Pb18), 19 (Pb3), and 27 (Pb01) exonic sRNA in EVs that corresponded to a specific exonic region (5 , 3 or middle) of a single mRNA (Table 3 and Figure 1A). For Pb18 EVs, 53% of those sequences were only in the forward orientation (F), while 68% and 56% were only in the reverse (R) direction in Pb3 and Pb01 EVs, respectively.
Four EV exonic sRNA sequences mapped to a common transcript in all Paracoccidioides isolates: calcium calmodulin-dependent protein kinase (PADG_07652), nucleotide binding (PADG_03535), transcription factor tfiiib complex subunit brf1 (PADG_00916), and a hypothetical protein (PADG_11439). As seen in Table 3, most of the exonic sRNAs were specific to EVs from Pb18 (89) and Pb01 (21), whereas only four were exclusive to Pb3 EVs. Ten fragments mapped to homologous transcripts in Pb18 and Pb3 EVs. Most of the target sequences have unknown functions, while the others are related to carbohydrate/protein/DNA metabolism, translation, oxidation-reduction, and the signaling process (Table 3).
In conclusion, although the number of putative exonic sRNA was 4 to 5-fold higher in Pb18 EVs, 10 target sequences seemed to be characteristic of P. brasiliensis, whereas 21 were exclusive to P. lutzii. These results reveal a previously unknown diversity in the composition of fungal EVs at the genus level. Table 3. Target mRNA of exonic sRNA sequences found in EVs from Pb18, Pb3, and Pb01 isolates. The position of sRNA in the exon is indicated as 5 , 3 , and M (middle), as well as the direction (F, forward; R, reverse). The sequences are grouped according with their functions.

Comparison of EV ncRNA Classes in Paracoccidioides EVs
The different classes of ncRNA contained in Pb18, Pb3, and Pb01 EVs were analyzed by aligning the <200-nt reads with the ncRNA database from the Saccharomyces Genome Database. The results revealed the presence of 71 different sequences of ncRNA in Paracoccidioides EVs, from which 15 were common to all isolates and 17 were shared by two of them ( Figure 1B). The most abundant class of ncRNA found in Paracoccidiodes EVs was the small nucleolar (sno)RNAs (33), followed by tRNAs (16), rRNAs (10), long ncRNAs (7), and small nuclear (sn)RNAs (4), as seen in Figure 1A. The snoRNA and tRNA were also the most abundant ncRNA populations described in our previous work not only for P. brasiliensis, but also for C. albicans, C. neoformans, and S. cerevisiae [21].
The number of total ncRNA detected in EVs varied with the isolate, being slightly more abundant in P. lutzii Pb01 (45) than in P. brasiliensis Pb3 (39) or Pb18 (35). The class profiles also differed between species. For P. brasiliensis EVs, snoRNA represented about 55% of the total, followed by rRNA (23% in Pb18 and 15% in Pb3), and 11 to 15% tRNA and snRNA ( Figure 1A). For P. lutzii EVs, we found similar proportions (20 to 31%) of snoRNA, rRNA, and tRNA, followed by 9 to 13% of snRNA and other ncRNAs ( Figure 1A).

Secondary Structure in the EVs.
In our previous publication [21], we identified putative miRNA-like (milRNA) sequences in fungal EVs that aligned with mature miRNA following a search in the miRNA database of all organisms (http://www.mirbase.org). In the present work, we searched for sRNA sequences (up to 50 nt) bearing secondary structure in order to have a more precise view of the presence of milRNA in Paracoccidioides EVs. We applied the secondary structure analysis to RNA fragments detected in EVs at high abundance (RPKM > 5000) and for further characterization we selected 42 RNA sequences bearing secondary structures with free energy values below −3.0 Kcal/mol. The sequences with the most negative values are shown in Figure 2. We performed an alignment with these 42 sequences to infer their function and identify in which strand they would align (Supplemental Table S1). Two were present in all EV samples, while 23 were found in both Pb3 and Pb01 EVs. Interestingly, 50% of the structured RNAs aligned at the complementary strand of the RNA, thus suggesting that they could act as miRNA or the transcripts (Supplemental Table S1). Most of the structured RNAs localize to transcripts that code for hypothetical proteins, however we could also identify RNAs that align to transcripts of kinesin-II, glutamine amidotransferase and lysine methyltransferase ( Figure 2). In our previous publication [21], we identified putative miRNA-like (milRNA) sequences in fungal EVs that aligned with mature miRNA following a search in the miRNA database of all organisms (http://www.mirbase.org). In the present work, we searched for sRNA sequences (up to 50 nt) bearing secondary structure in order to have a more precise view of the presence of milRNA in Paracoccidioides EVs. We applied the secondary structure analysis to RNA fragments detected in EVs at high abundance (RPKM > 5000) and for further characterization we selected 42 RNA sequences bearing secondary structures with free energy values below −3.0 Kcal/mol. The sequences with the most negative values are shown in Figure 2. We performed an alignment with these 42 sequences to infer their function and identify in which strand they would align (Supplemental Table 1). Two were present in all EV samples, while 23 were found in both Pb3 and Pb01 EVs. Interestingly, 50% of the structured RNAs aligned at the complementary strand of the RNA, thus suggesting that they could act as miRNA or the transcripts (Supplemental Table 1). Most of the structured RNAs localize to transcripts that code for hypothetical proteins, however we could also identify RNAs that align to transcripts of kinesin-II, glutamine amidotransferase and lysine methyltransferase ( Figure 2).

Paracoccidioides EVs Might Modulate the Transcriptome of Dendritic Cells
Dendritic cells make the bridge between the innate and adaptative immune responses, which ultimately define the course of the fungal infections [27]. The disease progression in paracoccidioidomycosis depends on dominant Th1 or Th2 types of immune response. The Th1-

Paracoccidioides EVs Might Modulate the Transcriptome of Dendritic Cells
Dendritic cells make the bridge between the innate and adaptative immune responses, which ultimately define the course of the fungal infections [27]. The disease progression in paracoccidioidomycosis depends on dominant Th1 or Th2 types of immune response. The Th1-driven inflammatory immune response is responsible for protection against the disease and the IL-12 expressed by dendritic cells is detrimental to stimulate this type of response [28]. In this context, we investigated if extracellular fractions containing EVs produced by P. brasiliensis Pb18 would affect gene expression in recipient dendritic cells. In order to do that, we characterized the transcriptome of murine monocyte-derived CD11c + (MoDC) cells co-cultivated for 48 h with the Pb18 yeast cells in a transwell system. We compared the sequences with those of the controls cultivated in the absence of fungal cells. In the transwell system, we had the fungal (upper compartment) and MoDC cell (bottom compartment) cultures communicating by a porous membrane that allowed the transit to the bottom well of soluble molecules and EVs up to 0.4 µm in size.
We observed that the indirect co-culture with Pb18 led to a slight alteration in the levels of mRNA expressed by MoDC cells when compared to the control ( Figure 3 and Table S2). We detected 20 upregulated and 28 downregulated mRNAs (Figure 3). We chose to follow highly stringent criteria in this analysis (similar levels of expression in both replicates, FDR ≤ 5% and log FC ≥ 2), considering we had duplicates, but not triplicates. Among the upregulated transcripts, three code for membrane proteins (Ankar, Unc13c, and Smim7) and two code for mitochondrial proteins associated to translation elongation (mt-Te and mt-Tv) (Table S2). Among the downregulated transcripts, there was 25% enrichment for transcripts related to gene expression regulation (transcription and translation). When we applied the hypergeometric test to the annotated transcripts for the gene set enrichment analysis, the p-value of the gene expression group was 2.40 × 10 5 , thus reinforcing the significance of this enrichment (Table S2). It was interesting to notice that three transcripts that code for transcription factors (Gabpb2, Pknox1, and Zfp575) are among the downregulated MoDC mRNAs in indirect co-cultures with Pb18 (Table S2), suggesting that sRNAs from the fungal EVs were potentially able to modulate the gene expression of the recipient cell. It is important to note that we have preliminary cytometry data showing that Pb18 yeast cells can actively release EVs under our transwell experimental conditions and that MoDC cells can internalize fungal EVs, as depicted in Figure S1. We observed that the indirect co-culture with Pb18 led to a slight alteration in the levels of mRNA expressed by MoDC cells when compared to the control ( Figure 3 and Table S2). We detected 20 upregulated and 28 downregulated mRNAs (Figure 3). We chose to follow highly stringent criteria in this analysis (similar levels of expression in both replicates, FDR ≤ 5% and log FC ≥ 2), considering we had duplicates, but not triplicates. Among the upregulated transcripts, three code for membrane proteins (Ankar, Unc13c, and Smim7) and two code for mitochondrial proteins associated to translation elongation (mt-Te and mt-Tv) (Table S2). Among the downregulated transcripts, there was 25% enrichment for transcripts related to gene expression regulation (transcription and translation). When we applied the hypergeometric test to the annotated transcripts for the gene set enrichment analysis, the p-value of the gene expression group was 2.40 × 10 5 , thus reinforcing the significance of this enrichment (Table S2). It was interesting to notice that three transcripts that code for transcription factors (Gabpb2, Pknox1, and Zfp575) are among the downregulated MoDC mRNAs in indirect co-cultures with Pb18 (Table S2), suggesting that sRNAs from the fungal EVs were potentially able to modulate the gene expression of the recipient cell. It is important to note that we have preliminary cytometry data showing that Pb18 yeast cells can actively release EVs under our transwell experimental conditions and that MoDC cells can internalize fungal EVs, as depicted in Figure S1.

Discussion
In the present work we have shown that the features of mRNA and sRNA sub-populations contained in EVs from P. brasiliensis Pb18, Pb3, and P. lutzii Pb01 differ considerably between species and also between representative isolates of the same species. The results suggest that Paracoccicioides

Discussion
In the present work we have shown that the features of mRNA and sRNA sub-populations contained in EVs from P. brasiliensis Pb18, Pb3, and P. lutzii Pb01 differ considerably between species and also between representative isolates of the same species. The results suggest that Paracoccicioides extracellular fractions containing EVs can modulate the transcription profile of dendritic cells. By using in vitro translation, we also reported for the first time in pathogenic fungi that the EV mRNAs are active and can be translated.
The mRNA sequences presently detected in the >200-nt RNA fraction of Paracoccidoides EVs varied in number with the isolate and only five orthologs were common to EVs from more than one isolate. We have recently compared the characteristics of mRNA sequences from EVs exported by two H. capsulatum isolates, specifically G186AR and G217B [29]. The latter isolate lacks cell wall alpha-1,3-glucan, which is a virulence factor [30]. The number of EV mRNA sequences varied from 93 in G186AR (mostly related to metabolic processes) to 31 in G217B (related to transport pathways possibly requiring vesicles, oxidation-reduction, and translation mechanisms). In the present analysis for P. brasiliensis, which is genetically related to Histoplasma [31], transcripts associated with DNA metabolism/biogenesis and protein modification were enriched in all samples, while those related with translation and transport prevailed, respectively, in Pb18 and Pb01.
While there are several pieces of evidence suggesting the role of fungal EVs on host cells, specially by interfering with the course of the immune response [32][33][34][35][36], it has recently been demonstrated that EV-like liposomal particles can cross the fungal cell wall inwards [37], thus opening the possibility that EVs can be uptaken and signal fungal cells. It has recently been demonstrated that EV-associated plant sRNA can silence virulence genes in a fungal pathogen, which agrees with the hypothesis that EVs mediate trans-kingdom regulation of gene expression [38]. On the other hand, Bielska et al. (2018) demonstrated that EVs produced by virulent Cryptococcus gattii were essential to signal quiescent strains within phagolysosomes located in distant body sites, which then became virulent; that phenomenon was apparently mediated by EV protein and RNA [39]. We have here detected active mRNA within Paracoccidioides EVs and we can envision that they might be uptaken by other neighboring or distant fungal cells and play a role in the host-fungal relationship by delivering virulence transcripts. Table 1 includes, for e.g., heat shock proteins like Hsp 70 (PAAG_08003, in Pb01 EVs) and Hsp 90-like (PADG_07715, in Pb18 EVs) that have a role in virulence of dimorphic fungi [40].
On the other hand, among the siRNA described in fungi, the exonic ex-siRNAs promote gene silencing in Mucor circinelloides and Fusarium graminearum [26,41,42]. In these species, the ex-siRNAs correspond to sense and antisense siRNAs converted into double-stranded RNA (dsRNA) via RNA-dependent RNA polymerase and Dicer. They regulate translation of the same mRNA that originated it. The presence of sense and antisense exonic sRNA-like in Paracoccidioides EVs suggests that these molecules could be involved in gene silencing via dsRNA. In the present work, we can point out the finding of exonic sRNA-like sequences, in Pb18 EVs, that map to α-amylase (PADG_04422), which is essential to the synthesis of the virulence factor α-glucan [30,43], and β-glucanase (PADG_04922), that cleaves dectin-1 ligand β-glucan present in the cell wall. We propose that they could have a role in regulating the expression of virulence genes in other fungal cells.
The comparison of the small ncRNA subtypes in Paracoccidioides EVs showed clear inter-species diversity. In P. lutzii, the percentage of snoRNA was about half that in P. brasiliensis, whereas for tRNA and other ncRNA it was, respectively, 2-and 4-fold higher. Consequently, 20 sequences were exclusive of P. lutzii EVs (mostly tRNA) versus six for Pb18 and thirteen for Pb3. In our experimental conditions, H. capsulatum EVs transport almost exclusively rRNA and tRNA [29]. The tRNA sequences prevailed in the EVs from isolate G186AR and most of them were not detected in EVs from the G217R strain. The fragments of tRNAs (the tRFs) are implicated in diverse processes in the cell, from regulation of cell viability, protein synthesis, apoptosis, to RNA metabolism, including turnover and stability [44]. The tRFs present in the EVs form Trypanosoma cruzi, the causing agent of Chagas disease, can be transferred to other parasites and promote cell communication and/or to host cells to modulate gene expression or facilitate infection [45].
We have previously detected 145 sequences in EVs from P. brasiliensis Pb18 that matched those of mature miRNA deposited in the miRNA database (mirbase). More recently, Curcio et al. [50] searched for matches of mil-RNAs already described in fungi, discarded those located in genes and looked for pre-miRNA secondary structure in the genome by considering their flanking nucleotides. They found that 11 mil-RNAs previously reported in fungal EVs matched the criteria for mil-RNA in the P. brasiliensis Pb18 genome. Besides, the authors found that in Pb18 the paralogous genes for Argonaute (1 and 2) and Dicer (1 and 2) (PADG_00716; PADG_03108; PADG_11946; PADG_07189) seem to be induced in the pathogenic yeast phase of the fungus, suggesting that mechanisms involving RNAi can be functional in this species [50]. We have presently shown partial secondary RNA structures in sRNA sequences from EVs that was differentially represented in the Paracoccidioides isolates. These sequences aligned to exons, notably in the reverse position. We found similar results for H. capsulatum [29], but the role of these molecules has to be investigated further. It was interesting to find, in Histoplasma EVs, a series of RNA-binding proteins and one of them, Snd1, is a component of the RNA-induced silencing complex (RISC) that is part of the RNAi machinery [51]. In Pb18 EVs [52] we have not found RNA-binding proteins that matched those found in Histoplasma EVs (data not shown).
Modulation of the host immune system by EV components has been described in C. neoformans, C. albicans, Malassezia sympodialis, P. brasiliensis and Sporothrix brasiliensis. The fungal EVs activate immune cells in vitro inducing the release of pro-inflammatory mediators, suggesting a role of fungal EVs in activating the immune system to respond to the fungal infection [32][33][34][35][36]. Importantly, in all studies performed so far, it was uncertain if the density of EVs used to stimulate host cells corresponded to that observed during physiologic and/or pathogenic conditions. Therefore, we opted for a transwell-based experimental system where the amount of fungal EVs corresponded to that physiologically produced by P. brasiliensis, even with the limitation that molecules not related to EVs may participate in the interaction with host cells. We observed that the Pknox1 and Gbpb2 transcription factors were highly downregulated (65-fold and 30-fold, respectively) in MoDC cells upon indirect co-coculture with P. brasiliensis in the transwell experiment. Pknox1 (Pbx/knotted 1 homeobox) belongs to the HOX family of transcription factors and is critical for the immune system homeostasis as it regulates the expression of IL-10 in macrophages and dendritic cells [53]. Gbpb2 (GA-binding protein subunit beta 2) is critical for T cell development and the expression regulation of IL-7R alpha [54]. IL-7 is associated with improved immune response against bacterial and viral infections, and also in fungal sepsis [55] and Pneumocystis infection [56]; however, a direct correlation of IL-7 and paracoccidioidomycosis has not been evaluated yet. On the other hand, IL-10 is a potent anti-inflammatory cytokine that is generally associated with increased infection in paracoccidioidomycosis [57]. Our experimental design using a transwell system did not allow distinguishing between the biological effects of EVs from other secreted molecules. However, considering our data suggesting that Paracoccidioides Pb18 produce EVs that are uptaken by MoDC in a transwell system ( Figure S1), we feel entitled to assume that EVs from P. brasiliensis Pb18 might help to modulate the MoDC response, possibly favoring the infection at the early stages of interaction with cells of the innate immune system. This modulation could be related to regulatory sRNA carried by the fungal EVs or even by surface ligands like DC-SIGN mannose ligands exposed in the EV surface [58].
Considering that sRNA classes are involved in transcription/translation modulation in a variety of systems and potentially also in Paracoccidioides, our general results may indicate that differences in virulence among fungal isolates could be related to their distinct EV RNA content. That hypothesis will hopefully be experimentally tested.
Supplementary Materials: The following files are available online, Figure S1: Uptake of EVs from Pb18 by MoDC cells in a transwell system, Table S1: List of genes from Pb18, Pb3, and Pb01 matching partially structured sRNA, as indicated. Separated lists contain those found exclusively in one isolate or shared by two or three of them, Table S2: MoDC genes differentially expressed (upregulated in green and downregulated in pink) after incubation with Pb18 cells in a transwell system, as analyzed using the RNA-seq strategy.

Conflicts of Interest:
The authors declare no conflict of interest.