Mitochondrial Genome Sequences of the Emerging Fungal Pathogen Candida auris.

Candida auris is an emerging fungal pathogen capable of causing invasive infections in humans. Since its first appearance around 1996, it has been isolated in countries spanning five continents. C. auris is a yeast that has the potential to cause outbreaks in hospitals, can survive in adverse conditions, including dry surfaces and high temperatures, and has been frequently misidentified by traditional methods. Furthermore, strains have been identified that are resistant to two and even all three of the main classes of antifungals currently in use. Several nuclear genome assemblies of C. auris have been published representing different clades and continents, yet until recently, the mitochondrial genomes (mtDNA chromosomes) of this species and the closely related species of C. haemulonii, C. duobushaemulonii, and C. pseudohaemulonii had not been analyzed in depth. We used reads from PacBio and Illumina sequencing to obtain a de novo reference assembly of the mitochondrial genome of the C. auris clade I isolate B8441 from Pakistan. This assembly has a total size of 28.2 kb and contains 13 core protein-coding genes, 25 tRNAs and the 12S and 16S ribosomal subunits. We then performed a comparative analysis by aligning Illumina reads of 129 other isolates from South Asia, Japan, South Africa, and South America with the B8441 reference. The clades of the phylogenetic tree we obtained from the aligned mtDNA sequences were consistent with those derived from the nuclear genome. The mitochondrial genome revealed a generally low genetic variation within clades, although the South Asian clade displayed two sub-branches including strains from both Pakistan and India. In particular, the 86 isolates from Colombia and Venezuela had mtDNA sequences that were all identical at the base level, i.e., a single conserved haplotype or mitochondrial background that exhibited characteristic differences from the Pakistan reference isolate B8441, such as a unique 25-nt insert that may affect function.

Updates and Comparative Analysis of the Mitochondrial Genomes of Paracoccidioides spp. Using Oxford Nanopore MinION Sequencing.

The mitochondrial genome of the Paracoccidioides brasiliensis reference isolate Pb18 was first sequenced and described by Cardoso et al. (2007), as a circular genome with a size of 71.3 kb and containing 14 protein coding genes, 25 tRNAs, and the large and small subunits of ribosomal RNA. Later in 2011, Desjardins et al. (2011) obtained partial assemblies of mitochondrial genomes of P. lutzii (Pb01), P. americana (Pb03), and P. brasiliensis sensu stricto (Pb18), although with a size of only 43.1 kb for Pb18. Sequencing errors or other limitations resulting from earlier technologies, and the advantages of NGS (short and long reads), prompted us to improve and update the mtDNA sequences and annotations of two Paracoccidioides species. Using Oxford Nanopore and Illumina read sequencing, we generated high-quality complete de novo mitochondrial genome assemblies and annotations for P. brasiliensis (Pb18) and P. americana (Pb03). Both assemblies were characterized by an unusually long spacer or intron region (>50 kb) between exons 2 and 3 of the nad5 gene, which was moderately conserved between Pb03 and Pb18 but not similar to other reported sequences, except for an unassigned contig in the 2011 assembly of Pb03. The reliability of the insert missing from previous mtDNA genome assemblies was confirmed by inspection of the individual Nanopore read sequences containing nad5 coding DNA, and experimentally by PCR for Pb18. We propose that the insert may aid replication initiation and may be excised to produce a smaller structural variant. The updated mtDNA genomes should enable more accurate SNP and other comparative or evolutionary analyses and primer/probe designs. A comparative analysis of the mtDNA from 32 isolates of Paracoccidioides spp., using the SNPs of the aligned mitochondrial genomes, showed groupings within the brasiliensis species complex that were largely consistent with previous findings from only five mitochondrial loci.

Genomic diversity of the human pathogen Paracoccidioides across the South American continent.

Paracoccidioidomycosis (PCM) is a life-threatening systemic mycosis widely reported in the Gran Chaco ecosystem. The disease is caused by different species from the genus Paracoccidioides, which are all endemic to South and Central America. Here, we sequenced and analyzed 31 isolates of Paracoccidioides across South America, with particular focus on isolates from Argentina and Paraguay. The de novo sequenced isolates were compared with publicly available genomes. Phylogenetics and population genomics revealed that PCM in Argentina and Paraguay is caused by three distinct Paracoccidioides genotypes, P. brasiliensis (S1a and S1b) and P. restrepiensis (PS3). P. brasiliensis S1a isolates from Argentina are frequently associated with chronic forms of the disease. Our results suggest the existence of extensive molecular polymorphism among Paracoccidioides species, and provide a framework to begin to dissect the connection between genotypic differences in the pathogen and the clinical outcomes of the disease.

Association of KIR2DL2 gene with anti-cyclic citrullinated protein antibodies for serodiagnosis in rheumatoid arthritis.

Rheumatoid arthritis is a clinical autoimmune syndrome that causes joint damage. The positive or negative anti-cyclic citrullinated protein (CCP) antibodies serodiagnosis differentiates two subsets of the disease, each with different genetic background. Previous studies have identified associations between KIR genes and rheumatoid arthritis but not with anti-CCP serodiagnosis. Therefore, we investigated the proportion of patients seropositive and seronegative to anti-CCP and its possible association with KIR (killer cell immunoglobulin-like receptor) genes. We included 100 patients with rheumatoid arthritis from western Mexico, who were determined for anti-CCP serodiagnosis by ELISA, and 16 KIR genes were genotyped by PCR-SSP. The proportion of seropositive anti-CCP patients was 83%, and they presented a higher frequency of KIR2DL2 genes than the seronegative group (73.6% vs. 46.2%, p = 0.044) which, in turn, presented a higher KIR2DL2-/KIR2DL3+ genotype frequency than the first ones (46.2% vs. 17.2%, p = 0.043). These results suggest different KIR genetic backgrounds for each subset of the disease according to anti-CCP serodiagnosis.

La artritis reumatoide es un síndrome clínico autoinmune que causa daño en las articulaciones. El serodiagnóstico positivo o negativo para anticuerpos proteicos anticíclicos citrulinados (CCP) diferencia dos subconjuntos de la enfermedad, cada uno con diferente fondo genético. Estudios previos han identificado asociaciones entre los genes killer cell immunoglobulin- like receptor (KIR) y la artritis reumatoide, pero no con el serodiagnóstico de anti-CCP. Por lo tanto, investigamos la proporción de seropositividad y seronegatividad anti-CCP y su posible asociación con genes KIR. Se incluyeron 100 pacientes con artritis reumatoide del occidente de México, a quienes se les determinó su serodiagnóstico anti-CCP por ELISA y también se les realizó genotipificación de 16 genes KIR por PCR-SSP. La proporción de pacientes seropositivos anti-CCP fue del 83% y presentaron una mayor frecuencia génica KIR2DL2 que el grupo seronegativo (73.6% vs. 46.2%, p = 0.044), estos últimos presentaron mayor frecuencia genotípica KIR2DL2-/KIR2DL3+ que los primeros (46.2% vs. 17.2%, p = 0.043). Los resultados sugieren diferente fondo genético KIR para cada subconjunto de la enfermedad, de acuerdo con el serodiagnóstico anti-CCP.

Association of kir2dl2 gene with anti-cyclic citrullinated protein antibodies for serodiagnosis in rheumatoid arthritis

Rheumatoid arthritis is a clinical autoimmune syndrome that causes joint damage. The positive or negative anti-cyclic citrullinated protein (CCP) antibodies serodiagnosis differentiates two subsets of the disease, each with different genetic background. Previous studies have identified associations between KIR genes and rheumatoid arthritis but not with anti-CCP serodiagnosis. Therefore, we investigated the proportion of patients seropositive and seronegative to anti-CCP and its possible association with KIR (killer cell immunoglobulin-like receptor) genes. We included 100 patients with rheumatoid arthritis from western Mexico, who were determined for anti-CCP serodiagnosis by ELISA, and 16 KIR genes were genotyped by PCR-SSP. The proportion of seropositive anti-CCP patients was 83%, and they presented a higher frequency of KIR2DL2 genes than the seronegative group (73.6% vs. 46.2%, p = 0.044) which, in turn, presented a higher KIR2DL2-/ KIR2DL3+ genotype frequency than the first ones (46.2% vs. 17.2%, p = 0.043). These results suggest different KIR genetic backgrounds for each subset of the disease according to anti-CCP serodiagnosis.

La artritis reumatoide es un síndrome clínico autoinmune que causa daño en las articulaciones. El serodiagnóstico positivo o negativo para anticuerpos proteicos anti-cíclicos citrulinados (CCP) diferencia dos subconjuntos de la enfermedad, cada uno con diferente fondo genético. Estudios previos han identificado asociaciones entre los genes killer cell immunoglobulin- like receptor (KIR) y la artritis reumatoide, pero no con el serodiagnóstico de anti-CCP. Por lo tanto, investigamos la proporción de seropositividad y seronegatividad anti-CCP y su posible asociación con genes KIR. Se incluyeron 100 pacientes con artritis reumatoide del occidente de México, a quienes se les determinó su serodiagnóstico anti-CCP por ELISA y también se les realizó genotipificación de 16 genes KIR por PCR-SSP. La proporción de pacientes seropositivos anti-CCP fue del 83% y presentaron una mayor frecuencia génica KIR2DL2 que el grupo seronegativo (73.6% vs. 46.2%, p = 0.044), estos últimos presentaron mayor frecuencia genotípica KIR2DL2-/KIR2DL3+ que los primeros (46.2% vs. 17.2%, p = 0.043). Los resultados sugieren diferente fondo genético KIR para cada subconjunto de la enfermedad, de acuerdo con el serodiagnóstico anti-CCP.

Draft Genome Sequences of Two Sporothrix schenckii Clinical Isolates Associated with Human Sporotrichosis in Colombia.

Sporothrix schenckii is a thermodimorphic fungal pathogen with a high genetic diversity. In this work, we present the assembly and similarity analysis of the whole-genome sequences of two clinical isolates from Colombia of S. schenckiisensu stricto.

Paracoccidioides spp. catalases and their role in antioxidant defense against host defense responses.

Dimorphic human pathogenic fungi interact with host effector cells resisting their microbicidal mechanisms. Yeast cells are able of surviving within the tough environment of the phagolysosome by expressing an antioxidant defense system that provides protection against host-derived reactive oxygen species (ROS). This includes the production of catalases (CATs). Here we identified and analyzed the role of CAT isoforms in Paracoccidioides, the etiological agent of paracoccidioidomycosis. Firstly, we found that one of these isoforms was absent in the closely related dimorphic pathogen Coccidioides and dermatophytes, but all of them were conserved in Paracoccidioides, Histoplasma and Blastomyces species. We probed the contribution of CATs in Paracoccidioides by determining the gene expression levels of each isoform through quantitative RT-qPCR, in both the yeast and mycelia phases, and during the morphological switch (transition and germination), as well as in response to oxidative agents and during interaction with neutrophils. PbCATP was preferentially expressed in the pathogenic yeast phase, and was associated to the response against exogenous H2O2. Therefore, we created and analyzed the virulence defects of a knockdown strain for this isoform, and found that CATP protects yeast cells from H2O2 generated in vitro and is relevant during lung infection. On the other hand, CATA and CATB seem to contribute to ROS homeostasis in Paracoccidioides cells, during endogenous oxidative stress. CAT isoforms in Paracoccidioides might be coordinately regulated during development and dimorphism, and differentially expressed in response to different stresses to control ROS homeostasis during the infectious process, contributing to the virulence of Paracoccidioides.

Identification and Analysis of the Role of Superoxide Dismutases Isoforms in the Pathogenesis of Paracoccidioides spp.

The ability of Paracoccidioides to defend itself against reactive oxygen species (ROS) produced by host effector cells is a prerequisite to survive. To counteract these radicals, Paracoccidioides expresses, among different antioxidant enzymes, superoxide dismutases (SODs). In this study, we identified six SODs isoforms encoded by the Paracoccidioides genome. We determined gene expression levels of representative isolates of the phylogenetic lineages of Paracoccidioides spp. (S1, PS2, PS3 and Pb01-like) using quantitative RT-PCR. Assays were carried out to analyze SOD gene expression of yeast cells, mycelia cells, the mycelia-to-yeast transition and the yeast-to-mycelia germination, as well as under treatment with oxidative agents and during interaction with phagocytic cells. We observed an increased expression of PbSOD1 and PbSOD3 during the transition process, exposure to oxidative agents and interaction with phagocytic cells, suggesting that these proteins could assist in combating the superoxide radicals generated during the host-pathogen interaction. Using PbSOD1 and PbSOD3 knockdown strains we showed these genes are involved in the response of the fungus against host effector cells, particularly the oxidative stress response, and in a mouse model of infection. Protein sequence analysis together with functional analysis of knockdown strains seem to suggest that PbSOD3 expression is linked with a pronounced extracellular activity while PbSOD1 seems more related to intracellular requirements of the fungus. Altogether, our data suggests that P. brasiliensis actively responds to the radicals generated endogenously during metabolism and counteracts the oxidative burst of immune cells by inducing the expression of SOD isoforms.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia.

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Genome update of the dimorphic human pathogenic fungi causing paracoccidioidomycosis.

Paracoccidiodomycosis (PCM) is a clinically important fungal disease that can acquire serious systemic forms and is caused by the thermodimorphic fungal Paracoccidioides spp. PCM is a tropical disease that is endemic in Latin America, where up to ten million people are infected; 80% of reported cases occur in Brazil, followed by Colombia and Venezuela. To enable genomic studies and to better characterize the pathogenesis of this dimorphic fungus, two reference strains of P. brasiliensis (Pb03, Pb18) and one strain of P. lutzii (Pb01) were sequenced [1]. While the initial draft assemblies were accurate in large scale structure and had high overall base quality, the sequences had frequent small scale defects such as poor quality stretches, unknown bases (N's), and artifactual deletions or nucleotide duplications, all of which caused larger scale errors in predicted gene structures. Since assembly consensus errors can now be addressed using next generation sequencing (NGS) in combination with recent methods allowing systematic assembly improvement, we re-sequenced the three reference strains of Paracoccidioides spp. using Illumina technology. We utilized the high sequencing depth to re-evaluate and improve the original assemblies generated from Sanger sequence reads, and obtained more complete and accurate reference assemblies. The new assemblies led to improved transcript predictions for the vast majority of genes of these reference strains, and often substantially corrected gene structures. These include several genes that are central to virulence or expressed during the pathogenic yeast stage in Paracoccidioides and other fungi, such as HSP90, RYP1-3, BAD1, catalase B, alpha-1,3-glucan synthase and the beta glucan synthase target gene FKS1. The improvement and validation of these reference sequences will now allow more accurate genome-based analyses. To our knowledge, this is one of the first reports of a fully automated and quality-assessed upgrade of a genome assembly and annotation for a non-model fungus.

Paracoccidioides brasiliensis PbP27 gene: knockdown procedures and functional characterization.

Paracoccidioides brasiliensis PbP27 gene encodes a protein localized in both the fungal cytoplasm and cell wall. The parasitic infectious form produces this protein preferentially with the gene's expression varying between the fungus phylogenetic species. The biological function of the native p27 has yet to be determined during either growth of the yeast or host infection. Therefore, in this study, through the use of antisense RNA technology and Agrobacterium tumefaciens-mediated transformation, we generated mitotically stable PbP27 mutants (PbP27 aRNA) with the goal to evaluate the role of p27 in the biology and virulence of this fungus. PbP27 expression was reduced 60-75% in mutants, as determined by real-time PCR in correlation with a decrease in p27 expression. No alterations in the growth curve or in the ability to shift from mycelia to yeast or from yeast to mycelia were observed in PbP27 aRNA strains; however, we did observe a reduction in cell vitality. Moreover, a decrease in cell viability of PbP27 aRNA yeast cells after interaction with IFN-γ-stimulated macrophages was detected. Based on these results, we propose that p27 plays a role in yeast cell architecture and represents one of the mechanisms employed by this fungus for its interaction with the monocyte/macrophage system.

Inhibition of PbGP43 expression may suggest that gp43 is a virulence factor in Paracoccidioides brasiliensis.

Glycoprotein gp43 is an immunodominant diagnostic antigen for paracoccidioidomycosis caused by Paracoccidioides brasiliensis. It is abundantly secreted in isolates such as Pb339. It is structurally related to beta-1,3-exoglucanases, however inactive. Its function in fungal biology is unknown, but it elicits humoral, innate and protective cellular immune responses; it binds to extracellular matrix-associated proteins. In this study we applied an antisense RNA (aRNA) technology and Agrobacterium tumefaciens-mediated transformation to generate mitotically stable PbGP43 mutants (PbGP43 aRNA) derived from wild type Pb339 to study its role in P. brasiliensis biology and during infection. Control PbEV was transformed with empty vector. Growth curve, cell vitality and morphology of PbGP43 aRNA mutants were indistinguishable from those of controls. PbGP43 expression was reduced 80-85% in mutants 1 and 2, as determined by real time PCR, correlating with a massive decrease in gp43 expression. This was shown by immunoblotting of culture supernatants revealed with anti-gp43 mouse monoclonal and rabbit polyclonal antibodies, and also by affinity-ligand assays of extracellular molecules with laminin and fibronectin. In vitro, there was significantly increased TNF-α production and reduced yeast recovery when PbGP43 aRNA1 was exposed to IFN-γ-stimulated macrophages, suggesting reduced binding/uptake and/or increased killing. In vivo, fungal burden in lungs of BALB/c mice infected with silenced mutant was negligible and associated with decreased lung ΙΛ-10 and IL-6. Therefore, our results correlated low gp43 expression with lower pathogenicity in mice, but that will be definitely proven when PbGP43 knockouts become available. This is the first study of gp43 using genetically modified P. brasiliensis.

Simazine transport in undisturbed soils from a vineyard at the Casablanca valley, Chile.

Simazine is a soil-active herbicide that has been applied worldwide in agricultural soils, being the second most commonly detected herbicide in groundwater and surface waters. Although its use has been restricted in many countries of Europe, it is still applied in many locations around the world in orchards, vineyards and forestry. Therefore, it is important to study its fate and transport in the environment. This paper investigates simazine transport in undisturbed bare soils from a vineyard at the Casablanca valley, Chile. In the study site, shallow groundwater tables (<1.0 m depth) and high simazine levels (>15 µg L(-1)) in the groundwater were observed and thus, there is potential for simazine to be transported further away through the saturated zone. The soils from the study site were characterized and the hydrodynamic transport parameters were determined. Column leaching experiments showed that the two-site chemical non-equilibrium model correctly represented simazine transport. It was found that 36.3% of the adsorption sites achieve instantaneous equilibrium and that the first-order kinetic rate of the non-equilibrium sites was 6.2 × 10(-3) h(-1). Hydrus 2D was used to predict the transport of simazine in the study site under natural field conditions. Simulation results showed that simazine concentrations at depths shallower than 2.1 m are above the maximum contaminant level of 4 µg L(-1) (defined by the U.S. Environmental Protection Agency). The timing of herbicide application was found to be important on simazine leaching and the main processes involved in simazine transport were degradation and adsorption, which accounted for 95.78 and 4.19% of the simulated mass of pesticide, respectively. A qualitative agreement in the timing and magnitude of simazine concentration was obtained between the simulations and the field data. Therefore, the model utilized in this investigation can be used to predict simazine transport and is a valuable tool to assess agricultural practices to minimize environmental impacts of simazine.

Involvement of the 90 kDa heat shock protein during adaptation of Paracoccidioides brasiliensis to different environmental conditions.

HSP90 is a molecular chaperone that participates in folding, stabilization, activation, and assembly of several proteins, all of which are key regulators in cell signaling. In dimorphic pathogenic fungi such as Paracoccidioides brasiliensis, the adaptation to a higher temperature, acid pH and oxidative stress, is an essential event for fungal survival and also for the establishing of the infectious process. To further understand the role of this protein, we used antisense RNA technology to generate a P. brasiliensis isolate with reduced PbHSP90 gene expression (PbHSP90-aRNA). Reduced expression of HSP90 decreased yeast cell viability during batch culture growth and increased susceptibility to acid pH environments and imposed oxidative stress. Also, PbHSP90-aRNA yeast cells presented reduced viability upon interaction with macrophages. The findings presented here suggest a protective role for HSP90 during adaptation to hostile environments, one that promotes survival of the fungus during host-pathogen interactions.

Transport of simazine in unsaturated sandy soil and predictions of its leaching under hypothetical field conditions.

The potential contamination of groundwater by herbicides is often controlled by processes in the vadose zone, through which herbicides travel before entering groundwater. In the vadose zone, both physical and chemical processes affect the fate and transport of herbicides, therefore it is important to represent these processes by mathematical models to predict contaminant movement. To simulate the movement of simazine, a herbicide commonly used in Chilean vineyards, batch and miscible displacement column experiments were performed on a disturbed sandy soil to quantify the primary parameters and processes of simazine transport. Chloride (Cl(-)) was used as a non-reactive tracer, and simazine as the reactive tracer. The Hydrus-1D model was used to estimate the parameters by inversion from the breakthrough curves of the columns and to evaluate the potential groundwater contamination in a sandy soil from the Casablanca Valley, Chile. The two-site, chemical non-equilibrium model was observed to best represent the experimental results of the miscible displacement experiments in laboratory soil columns. Predictions of transport under hypothetical field conditions using the same soil from the column experiments were made for 40 years by applying herbicide during the first 20 years, and then halting the application and considering different rates of groundwater recharge. For recharge rates smaller than 84 mm year(-1), the predicted concentration of simazine at a depth of 1 m is below the U.S. EPA's maximum contaminant levels (4 microg L(-1)). After eight years of application at a groundwater recharge rate of 180 mm year(-1) (approximately 50% of the annual rainfall), simazine was found to reach the groundwater (located at 1 m depth) at a higher concentration (more than 40 microg L(-1)) than the existing guidelines in the USA and Europe.