Comprehensive re-analysis of hairpin small RNAs in fungi reveals loci with conserved links.

RNA interference is an ancient mechanism with many regulatory roles in eukaryotic genomes, with small RNAs acting as their functional element. While there is a wide array of classes of small-RNA-producing loci, those resulting from stem-loop structures (hairpins) have received profuse attention. Such is the case of microRNAs (miRNAs), which have distinct roles in plants and animals. Fungi also produce small RNAs, and several publications have identified miRNAs and miRNA-like (mi/milRNA) hairpin RNAs in diverse fungal species using deep sequencing technologies. Despite this relevant source of information, relatively little is known about mi/milRNA features in fungi, mostly due to a lack of established criteria for their annotation. To systematically assess mi/milRNA characteristics and annotation confidence, we searched for publications describing mi/milRNA loci and re-assessed the annotations for 41 fungal species. We extracted and normalized the annotation data for 1727 reported mi/milRNA loci and determined their abundance profiles, concluding that less than half of the reported loci passed basic standards used for hairpin RNA discovery. We found that fungal mi/milRNA are generally more similar in size to animal miRNAs and were frequently associated with protein-coding genes. The compiled genomic analyses identified 25 mi/milRNA loci conserved in multiple species. Our pipeline allowed us to build a general hierarchy of locus quality, identifying more than 150 loci with high-quality annotations. We provide a centralized annotation of identified mi/milRNA hairpin RNAs in fungi which will serve as a resource for future research and advance in understanding the characteristics and functions of mi/milRNAs in fungal organisms.

Comparative transcriptome analysis reveals candidate genes related to cadmium accumulation and tolerance in two almond mushroom (Agaricus brasiliensis) strains with contrasting cadmium tolerance.

Cadmium (Cd) is a toxic metal occurring in the environment naturally. Almond mushroom (Agaricus brasiliensis) is a well-known cultivated edible and medicinal mushroom. In the past few decades, Cd accumulation in A.brasiliensis has received increasing attention. However, the molecular mechanisms of Cd-accumulation in A. brasiliensis are still unclear. In this paper, a comparative transcriptome of two A.brasiliensis strains with contrasting Cd accumulation and tolerance was performed to identify Cd-responsive genes possibly responsible for low Cd-accumulation and high Cd-tolerance. Using low Cd-accumulating and Cd-tolerant (J77) and high Cd-accumulating and Cd-sensitive (J1) A.brasiliensis strains, we investigated 0, 2 and 5 mg L-1 Cd-effects on mycelium growth, Cd-accumulation and transcriptome revealed by RNA-Seq. A total of 57,884 unigenes were obtained. Far less Cd-responsive genes were identified in J77 mycelia than those in J1 mycelia (e.g., ABC transporters, ZIP Zn transporter, Glutathione S-transferase and Cation efflux (CE) family). The higher Cd-accumulation in J1 mycelia might be due to Cd-induced upregulation of ZIP Zn transporter. Cd impaired cell wall, cell cycle, DNA replication and repair, thus decreasing J1 mycelium growth. Cd-stimulated production of sulfur-containing compounds, polysaccharides, organic acids, trehalose, ATP and NADPH, and sequestration of Cd might be adaptive responses of J1 mycelia to the increased Cd-accumulation. DNA replication and repair had better stability under 2 mg L-1 Cd, but greater positive modifications under 5 mg L-1 Cd. Better stability of DNA replication and repair, better cell wall and cell cycle stability might account for the higher Cd-tolerance of J77 mycelia. Our findings provide a comprehensive set of DEGs influenced by Cd stress; and shed light on molecular mechanism of A.brasiliensis Cd accumulation and Cd tolerance.

A method for the extraction of high quality fungal RNA suitable for RNA-seq.

Transcriptomic analysis is an OMICs technology that is becoming indispensable to understand and get a complete picture of cell functioning and adaptation to the environmental cues the cell is continuously receiving. Among the techniques available to perform transcriptomics, RNA-seq is becoming the method of choice. The quality of the RNA used for the generation of cDNA libraries and subsequent sequencing is crucial for the success of the process. Good RNA-seq performance is often limited by problems such as low RNA yield and/or integrity, RNA stability, and contamination with DNA, salts or chemicals. RNA isolation from fungi usually faces these problems and is particularly sensitive to degradation due to the high RNase activity content present in many species. Here we describe the development of a robust, highly reproducible and simple RNA purification method for filamentous fungi, which combines various strategies to get fully DNA-free RNA samples of high purity and integrity without the need to use a DNase I digestion step. The obtained RNA samples complied with all required standards to be used for RNA-seq and showed an excellent performance when subjected to Illumina-HiSeq 2500.

Characterization of Nosema ceranae Genetic Variants from Different Geographic Origins.

In recent years, large-scale colony losses of honey bees (Apis mellifera) have been reported and the infection with the microsporidia Nosema ceranae has been involved. However, the effect of N. ceranae at the colony level and its role in colony losses vary in different geographic areas. This difference may be related to the presence of multiple N. ceranae genetic variants resulting in different biological consequences. In this study, we analyzed the genetic diversity of 75 N. ceranae samples obtained from 13 countries and Hawaii through inter-sequence single repetition (ISSR) and evaluated if two of these genetic variants triggered different immune responses when infecting Apis mellifera iberiensis. The genetic diversity analysis showed that 41% of the samples had the same DNA amplification pattern, including samples from most European countries except Spain, while the remaining samples showed high variability. Infection assays were performed to analyze the infection levels and the immune response of bees infected with N. ceranae from Spain and Uruguay. The infected bees presented similar infection levels, and both isolates downregulated the expression of abaecin, confirming the ability of the microsporidia to depress the immune response. Only N. ceranae from Uruguay downregulated the expression level of imd compared to control bees. On the other hand, both genetic variants triggered different expression levels of lysozyme. As imd and lysozyme play important roles in the response to pathogens, these results could reflect differences in the biological consequences of N. ceranae variants in A. mellifera infection.

Isolation and characterization of extrachromosomal double-stranded RNA elements in Xanthophyllomyces dendrorhous.

Double-stranded RNA (dsRNA) molecules are widely found in yeasts and filamentous fungi. It has been suggested that may play important roles in the evolution of eukaryote genomes and may be a valuable tool in yeast typing. The characterization of these extrachromosomal genetic elements is usually a laborious process, especially when trying to analyze a large number of samples. In this chapter, we describe a simple method to isolate dsRNA elements from yeasts using low amounts of starting material, and their application to different Xanthophyllomyces dendrorhous strains. Furthermore, the methodologies for enzymatic and hybridization characterizations, and quantification of relative dsRNA abundance are detailed.

Primary mechanism of apoptosis induction in a leukemia cell line by fraction FA-2-b-ß prepared from the mushroom Agaricus blazei Murill

Agaricus blazei Murill is a native Brazilian mushroom which functions primarily as an anticancer substance in transplanted mouse tumors. However, the mechanism underlying this function of A. blazei Murill remains obscure. The present study was carried out to investigate the effect of fraction FA-2-b-ß, an RNA-protein complex isolated from A. blazei Murill, on human leukemia HL-60 cells in vitro. Typical apoptotic characteristics were determined by morphological methods using DNA agarose gel electrophoresis and flow cytometry. The growth suppressive effect of fraction FA-2-b-ß on HL-60 cells in vitro occurred in a dose- (5-80 mug/mL) and time-dependent (24-96 h) manner. The proliferation of HL-60 cells (1 x 10(5) cells/mL) treated with 40 mug/mL of fraction FA-2-b-ß for 24-96 h and with 5-80 mug/mL for 96 h resulted in inhibitory rates ranging from 8 to 54.5 percent, and from 4.9 to 86.3 percent, respectively. Both telomerase activity determined by TRAP-ELISA and mRNA expression of the caspase-3 gene detected by RT-PCR were increased in HL-60 cells during fraction FA-2-b-ß treatment. The rate of apoptosis correlated negatively with the decrease of telomerase activity (r = 0.926, P < 0.05), but correlated positively with caspase-3 mRNA expression (r = 0.926, P < 0.05). These data show that fraction FA-2-b-ß can induce HL-60 cell apoptosis and that the combined effect of down-regulation of telomerase activity and up-regulation of mRNA expression of the caspase-3 gene could be the primary mechanism of induction of apoptosis. These findings provide strong evidence that fraction FA-2-b-ß could be of interest for the clinical treatment of acute leukemia.

Primary mechanism of apoptosis induction in a leukemia cell line by fraction FA-2-b-ss prepared from the mushroom Agaricus blazei Murill.

Agaricus blazei Murill is a native Brazilian mushroom which functions primarily as an anticancer substance in transplanted mouse tumors. However, the mechanism underlying this function of A. blazei Murill remains obscure. The present study was carried out to investigate the effect of fraction FA-2-b-ss, an RNA-protein complex isolated from A. blazei Murill, on human leukemia HL-60 cells in vitro. Typical apoptotic characteristics were determined by morphological methods using DNA agarose gel electrophoresis and flow cytometry. The growth suppressive effect of fraction FA-2-b-ss on HL-60 cells in vitro occurred in a dose- (5-80 microg/mL) and time-dependent (24-96 h) manner. The proliferation of HL-60 cells (1 x 10(5) cells/mL) treated with 40 microg/mL of fraction FA-2-b-ss for 24-96 h and with 5-80 microg/mL for 96 h resulted in inhibitory rates ranging from 8 to 54.5%, and from 4.9 to 86.3%, respectively. Both telomerase activity determined by TRAP-ELISA and mRNA expression of the caspase-3 gene detected by RT-PCR were increased in HL-60 cells during fraction FA-2-b-ss treatment. The rate of apoptosis correlated negatively with the decrease of telomerase activity (r = 0.926, P < 0.05), but correlated positively with caspase-3 mRNA expression (r = 0.926, P < 0.05). These data show that fraction FA-2-b-ss can induce HL-60 cell apoptosis and that the combined effect of down-regulation of telomerase activity and up-regulation of mRNA expression of the caspase-3 gene could be the primary mechanism of induction of apoptosis. These findings provide strong evidence that fraction FA-2-b-ss could be of interest for the clinical treatment of acute leukemia.

Cloning, characterization and expression of a calnexin homologue from the pathogenic fungus Paracoccidioides brasiliensis.

We report the cloning of a Paracoccidioides brasiliensis cDNA, here named PbCnx, encoding the homologue of the endoplasmic reticulum molecular chaperone calnexin. Calnexin specifically recognizes monoglucosylated glycoproteins in the endoplasmic reticulum, thus being an essential component of the complex that interacts with the folded state of nascent secreted glycoproteins. The PbCnx open reading frame was found in a 1701 base pair (bp) fragment that encodes a 567 amino acid protein with an estimated mass of 62 680 Da. Northern and Southern blot hybridizations showed that PbCnx is encoded by a single, or a low number of, gene copies. PbCnx contains the hallmark KPEDWD motifs that are found in all members of the calnexin/calreticulin family proteins. A cDNA-encoding PbCnx was overexpressed as recombinant protein in Escherichia coli. The purified recombinant PbCnx was recognized by 6 out of 10 sera from PCM patients, a result that rules out its possible consideration for further use in diagnosis. Using confocal microscopy with anti-PbCnx mouse serum against yeast forms, a cytoplasmic staining pattern was observed.

A new Paracoccidioides brasiliensis 70-kDa heat shock protein reacts with sera from paracoccidioidomycosis patients.

A cDNA coding for a new member of the 70-kDa heat shock proteins (HSP70) family from the dimorphic and pathogenic fungus, Paracoccidioides brasiliensis, was cloned and characterized. The cDNA-deduced sequence coded for 655 amino acid residues and showed 95% identity to a previously described P. brasiliensis hsp70 gene. Cytoplasmic and typical nuclear localization signals, which indicate induction upon stress, were identified in the deduced peptide. The complete hsp70 cDNA coding region was cloned into a pGEX 4T-3 plasmid and expressed in Escherichia coli as a glutathione-S-transferase-tagged fusion protein. The recombinant protein reacted with a rabbit polyclonal antibody against HSP70. Western immunoblot experiments demonstrated that sera from paracoccidioidomycosis patients recognized the purified recombinant protein, suggesting an immunological role for this protein in the infectious process. The antigenicity analysis of rHSP70 detected three internal peptides that could act as activators of T-cell proliferation.

Cloning and functional analysis of the orotidine-5'-phosphate decarboxylase gene (PbrURA3) of the pathogenic fungus Paracoccidioides brasiliensis.

A genomic clone encoding the Paracoccidioides brasiliensis orotidine monophosphate decarboxylase gene (PbrURA3) was isolated by screening a subgenomic plasmid DNA library of this fungus, using a PCR amplification product of the gene as a probe. Sequence analysis revealed that the gene contains an open reading frame of 855 bp with a single intron (162 bp), and encodes a putative 285 amino acids polypeptide of estimated molecular weight 31.1 kDa and isoelectric point 6.5. The deduced amino acid sequence predicted a 73.4% identity with orotidine monophosphate decarboxylase of Aspergillus nidulans. Functionality of the gene was demonstrated by transformation into a Saccharomyces cerevisiae ura3 null mutant.

Quantitative analysis of wine yeast gene expression profiles under winemaking conditions.

Wine fermentation is a dynamic and complex process in which the yeast cell is subjected to multiple stress conditions. A successful adaptation involves changes in gene expression profiles where a large number of genes are up- or downregulated. Functional genomic approaches are commonly used to obtain global gene expression profiles, thereby providing a comprehensive view of yeast physiology. We used SAGE to quantify gene expression profiles in an industrial strain of Saccharomyces cerevisiae under winemaking conditions. The transcriptome of wine yeast was analysed at three stages during the fermentation process, mid-exponential phase, and early- and late-stationary phases. Upon correlation with the yeast genome, we found three classes of transcripts: (a) sequences that corresponded to ORFs; (b) expressed sequences from intergenic regions; and (c) messengers that did not match the published reference yeast genome. In all fermentation phases studied, the most highly expressed genes related to energy production and stress response. For many pathways, including glycolysis, different transcript levels were observed during each phase. Different isoenzymes, including hexose transporters (HXT), were differentially induced, depending on the growth phase. About 10% of transcripts matched non-annotated ORF regions within the yeast genome and could correspond to small novel genes originally omitted in the first gene annotation effort. Up to 22% of transcripts, particularly at late-stationary phase, did not match any known location within the genome. As the available reference yeast genome was obtained from a laboratory strain, these expressed sequences could represent genes only expressed by an industrial yeast strain. Further studies are necessary to identify the role of these potential genes during wine fermentation.

Cloning and expression analysis of the ornithine decarboxylase gene (PbrODC) of the pathogenic fungus Paracoccidioides brasiliensis.

We describe the isolation and sequencing of PbrODC, the gene encoding ornithine decarboxylase (ODC) in Paracoccidioides brasiliensis. The gene contains a single open reading frame made of 1413 bp with a single intron (72 bp), and encodes a 447 amino acid polypeptide with a predicted molecular weight of 50.0 kDa, an isoelectric point of 4.9 and a high similarity to other fungal ornithine decarboxylases. Functionality of the gene was demonstrated by transformation into a Saccharomyces cerevisiae odc null mutant. A phylogenetic tree generated with several fungal ODCs provided additional evidence to favour a taxonomic position for P. brasiliensis as an ascomycetous fungus, belonging to the order Onygenales. Expression of the PbrODC gene was determined by Northern analyses during growth of the mycelial and yeast forms, and through the temperature-regulated dimorphic transition between these two extreme phases. Expression of PbrODC remained constant at all stages of the fungal growth, and did not correlate with a previously observed increase in the activity of ornithine decarboxylase at the onset of the budding process in both yeast growth and mycelium-to-yeast transition. Accordingly, post-transcriptional regulation for the product of PbrODC is suggested.

The pfr1 gene from the human pathogenic fungus Paracoccidioides brasiliensis encodes a half-ABC transporter that is transcribed in response to treatment with fluconazole.

We have isolated a gene that encodes a half-ABC-transporter, designated Pfr1, from the dimorphic human pathogenic fungus Paracoccidioides brasiliensis, which has high identity with members of the ABC-superfamily involved in multidrug resistance. The pfr1 gene is predicted to encode a 827 amino acid protein that, in common with mammalian Mdr1, has a TM-NBD topology. The transcription of the pfr1 gene is induced by the triazole drug fluconazole but not by amphotericin B, suggesting a role in transport-mediated azole resistance. However, Pfr1 has greatest identity to the mitochondrial ABC transporters Mdl1 and Mdl2 from Saccharomyces cerevisiae and mammalian ABC-me, with identities of 47.2%, 40.6% and 39.5%, respectively, over the length of these proteins. Furthermore, the N-terminus of Pfr1 is rich in positively charged residues, a feature of mitochondrial targeting sequences. Considering these features, it seems likely that Pfr1 is a mitochondrial protein. Previous studies have revealed that the acquisition of azole resistance in S. cerevisiae is linked to mitochondrial loss and, conversely, that mitochondrial dysfunction can lead to the upregulation of PDR transporters mediated by the transcription factor Pdr3. Our studies suggest that a mitochondrial ABC transporter is induced as part of the cellular response to drug treatment. The promoter region of pfr1 contains a PDRE-like consensus sequence to which Pdr3 binds, which may be the element responsible for the upregulation of Pfr1 in response to fluconazole. The nucleotide binding domain of Pfr1 was expressed and purified from Escherichia coli and shown to retain ATPase activity, consistent with Pfr1 functioning as a homodimeric transport ATPase.