RESUMO

BACKGROUND: Zoonotic sporotrichosis is a neglected fungal disease, whereby outbreaks are primarily driven by Sporothrix brasiliensis and linked to cat-to-human transmission. To understand the emergence and spread of sporotrichosis in Brazil, the epicentre of the current epidemic in South America, we aimed to conduct whole-genome sequencing (WGS) to describe the genomic epidemiology. METHODS: In this genomic epidemiology study, we included Sporothrix spp isolates from sporotrichosis cases from Brazil, Colombia, and the USA. We conducted WGS using Illumina NovaSeq on isolates collected by three laboratories in Brazil from humans and cats with sporotrichosis between 2013 and 2022. All isolates that were confirmed to be Sporothrix genus by internal transcribed spacer or beta-tubulin PCR sequencing were included in this study. We downloaded eight Sporothrix genome sequences from the National Center for Biotechnology Information (six from Brazil, two from Colombia). Three Sporothrix spp genome sequences from the USA were generated by the US Centers for Disease Control and Prevention as part of this study. We did phylogenetic analyses and correlated geographical and temporal case distribution with genotypic features of Sporothrix spp isolates. FINDINGS: 72 Sporothrix spp isolates from 55 human and 17 animal sporotrichosis cases were included: 67 (93%) were from Brazil, two (3%) from Colombia, and three (4%) from the USA. Cases spanned from 1999 to 2022. Most (61 [85%]) isolates were S brasiliensis, and all were reported from Brazil. Ten (14%) were Sporothrix schenckii and were reported from Brazil, USA, and Colombia. For S schenckii isolates, two distinct clades were observed wherein isolates clustered by geography. For S brasiliensis isolates, five clades separated by more than 100 000 single-nucleotide polymorphisms were observed. Among the five S brasiliensis clades, clades A and C contained isolates from both human and cat cases, and clade A contained isolates from six different states in Brazil. Compared with S brasiliensis isolates, larger genetic diversity was observed among S schenckii isolates from animal and human cases within a clade. INTERPRETATION: Our results suggest that the ongoing epidemic driven by S brasiliensis in Brazil represents several, independent emergence events followed by animal-to-animal and animal-to human transmission within and between Brazilian states. These results describe how S brasiliensis can emerge and spread within a country. FUNDING: Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil; the São Paulo Research Foundation; Productivity in Research fellowships by the National Council for Scientific and Technological Development, and Ministry of Science and Technology of Brazil.

Assuntos

RESUMO

Sporothrix brasiliensis is an emerging zoonotic fungal pathogen that can be difficult to treat. Antifungal susceptibility testing was performed on the mold phase of a convenience sample of 61 Sporothrix spp. isolates from human and cat sporotrichosis cases in Brazil using the Clinical and Laboratory Standards Institute standard M38. A bimodal distribution of azole susceptibility was observed with 50% (28/56) of S. brasiliensis isolates showing elevated itraconazole minimum inhibitory concentrations ≥16 µg/mL. Phylogenetic analysis found the in vitro resistant isolates were not clonal and were distributed across three different S. brasiliensis clades. Single nucleotide polymorphism (SNP) analysis was performed to identify potential mechanisms of in vitro resistance. Two of the 28 resistant isolates (MIC ≥16 mg/L) had a polymorphism in the cytochrome P450 gene, cyp51, corresponding to the well-known G448S substitution inducing azole resistance in Aspergillus fumigatus. SNPs corresponding to other known mechanisms of azole resistance were not identified in the remaining 26 in vitro resistant isolates.

Assuntos

RESUMO

Since 2016, in Colombia, ongoing transmission of Candida auris has been reported in multiple cities. Here, we provide an updated description of C. auris genomic epidemiology and the dynamics of antifungal resistance in Colombia. We sequenced 99 isolates from C. auris cases with collection dates ranging from June 2016 to January 2021; the resulting sequences coupled with 103 previously generated sequences from C. auris cases were described in a phylogenetic analysis. All C. auris cases were clade IV. Of the 182 isolates with antifungal susceptibility data, 67 (37%) were resistant to fluconazole, and 39 (21%) were resistant to amphotericin B. Isolates predominately clustered by country except for 16 isolates from Bogotá, Colombia, which grouped with isolates from Venezuela. The largest cluster (N = 166 isolates) contained two subgroups. The first subgroup contained 26 isolates, mainly from César; of these, 85% (N = 22) were resistant to fluconazole. The second subgroup consisted of 47 isolates from the north coast; of these, 81% (N = 38) were resistant to amphotericin B. Mutations in the ERG11 and TAC1B genes were identified in fluconazole-resistant isolates. This work describes molecular mechanisms associated with C. auris antifungal resistance in Colombia. Overall, C. auris cases from different geographic locations in Colombia exhibited high genetic relatedness, suggesting continued transmission between cities since 2016. These findings also suggest a lack of or minimal introductions of different clades of C. auris into Colombia. IMPORTANCE: Candida auris is an emerging fungus that presents a serious global health threat and has caused multiple outbreaks in Colombia. This work discusses the likelihood of introductions and local transmission of C. auris and provides an updated description of the molecular mechanisms associated with antifungal resistance in Colombia. Efforts like this provide information about the evolving C. auris burden that could help guide public health strategies to control C. auris spread.

Assuntos

RESUMO

Spider venoms constitute a trove of novel peptides with biotechnological interest. Paucity of next-generation-sequencing (NGS) data generation has led to a description of less than 1% of these peptides. Increasing evidence supports the underestimation of the assembled genes a single transcriptome assembler can predict. Here, the transcriptome of the venom gland of the spider Pamphobeteus verdolaga was re-assembled, using three free access algorithms, Trinity, SOAPdenovo-Trans, and SPAdes, to obtain a more complete annotation. Assembler's performance was evaluated by contig number, N50, read representation on the assembly, and BUSCO's terms retrieval against the arthropod dataset. Out of all the assembled sequences with all software, 39.26% were common between the three assemblers, and 27.88% were uniquely assembled by Trinity, while 27.65% were uniquely assembled by SPAdes. The non-redundant merging of all three assemblies' output permitted the annotation of 9232 sequences, which was 23% more when compared to each software and 28% more when compared to the previous P. verdolaga annotation; moreover, the description of 65 novel theraphotoxins was possible. In the generation of data for non-model organisms, as well as in the search for novel peptides with biotechnological interest, it is highly recommended to employ at least two different transcriptome assemblers.

Assuntos

RESUMO

BACKGROUND: The band 9p21.3 contains an established genomic risk zone for cardiovascular disease (CVD). Since the initial 2007 Wellcome Trust Case Control Consortium study (WTCCC), the increased CVD risk associated with 9p21.3 has been confirmed by multiple studies in different continents. However, many years later there was still no confirmed report of a corresponding association of 9p21.3 with hypertension, a major CV risk factor, nor with blood pressure (BP). THEORY: In this contribution, we review the bipartite haplotype structure of the 9p21.3 risk locus: one block is devoid of protein-coding genes but contains the lead CVD risk SNPs, while the other block contains the first exon and regulatory DNA of the gene for the cell cycle inhibitor p15. We consider how findings from molecular biology offer possibilities of an involvement of p15 in hypertension etiology, with expression of the p15 gene modulated by genetic variation from within the 9p21.3 risk locus. RESULTS: We present original results from a Colombian study revealing moderate but persistent association signals for BP and hypertension within the classic 9p21.3 CVD risk locus. These SNPs are mostly confined to a 'hypertension island' that spans less than 60 kb and coincides with the p15 haplotype block. We find confirmation in data originating from much larger, recent European BP studies, albeit with opposite effect directions. CONCLUSION: Although more work will be needed to elucidate possible mechanisms, previous findings and new data prompt reconsidering the question of how variation in 9p21.3 might influence hypertension components of cardiovascular risk.

RESUMO

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.

RESUMO

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.

RESUMO

The LUFS domain (LUG/LUH, Flo8, single-strand DNA-binding protein [SSBP]) is a well-conserved and apparently ancient region found in diverse proteins and taxa. This domain, which has as its most obvious structural feature a series of three helices, has been identified in transcriptional regulator proteins of animals, plants, and fungi. Recently, in these pages (Wang et al., Protein Sci., 2019, 28:788-793), the first crystal structure of a LUFS domain was reported, for the human SSBP2, a transcriptional repressor. We briefly address how the new insights into LUFS structures might contribute to a better understanding of an important transcriptional activator of yeasts that contains the LUFS domain, Flo8, and consider how a focus on the LUFS domain and its variation could help us to understand etiologies of drug resistance in a recently emerged pathogenic fungus, Candida auris.

Assuntos

RESUMO

Background: Candida auris is a multidrug-resistant yeast associated with hospital outbreaks worldwide. During 2015-2016, multiple outbreaks were reported in Colombia. We aimed to understand the extent of contamination in healthcare settings and to characterize the molecular epidemiology of C. auris in Colombia. Methods: We sampled patients, patient contacts, healthcare workers, and the environment in 4 hospitals with recent C. auris outbreaks. Using standardized protocols, people were swabbed at different body sites. Patient and procedure rooms were sectioned into 4 zones and surfaces were swabbed. We performed whole-genome sequencing (WGS) and antifungal susceptibility testing (AFST) on all isolates. Results: Seven of the 17 (41%) people swabbed were found to be colonized. Candida auris was isolated from 37 of 322 (11%) environmental samples. These were collected from a variety of items in all 4 zones. WGS and AFST revealed that although isolates were similar throughout the country, isolates from the northern region were genetically distinct and more resistant to amphotericin B (AmB) than the isolates from central Colombia. Four novel nonsynonymous mutations were found to be significantly associated with AmB resistance. Conclusions: Our results show that extensive C. auris contamination can occur and highlight the importance of adherence to appropriate infection control practices and disinfection strategies. Observed genetic diversity supports healthcare transmission and a recent expansion of C. auris within Colombia with divergent AmB susceptibility.

Assuntos

RESUMO

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.

RESUMO

The presence of repetitive or non-unique DNA persisting over sizable regions of a eukaryotic genome can hinder the genome's successful de novo assembly from short reads: ambiguities in assigning genome locations to the non-unique subsequences can result in premature termination of contigs and thus overfragmented assemblies. Fungal mitochondrial (mtDNA) genomes are compact (typically less than 100 kb), yet often contain short non-unique sequences that can be shown to impede their successful de novo assembly in silico. Such repeats can also confuse processes in the cell in vivo. A well-studied example is ectopic (out-of-register, illegitimate) recombination associated with repeat pairs, which can lead to deletion of functionally important genes that are located between the repeats. Repeats that remain conserved over micro- or macroevolutionary timescales despite such risks may indicate functionally or structurally (e.g., for replication) important regions. This principle could form the basis of a mining strategy for accelerating discovery of function in genome sequences. We present here our screening of a sample of 11 fully sequenced fungal mitochondrial genomes by observing where exact k-mer repeats occurred several times; initial analyses motivated us to focus on 17-mers occurring more than three times. Based on the diverse repeats we observe, we propose that such screening may serve as an efficient expedient for gaining a rapid but representative first insight into the repeat landscapes of sparsely characterized mitochondrial chromosomes. Our matching of the flagged repeats to previously reported regions of interest supports the idea that systems of persisting, non-trivial repeats in genomes can often highlight features meriting further attention.

Assuntos

RESUMO

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.

Assuntos

RESUMO

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.

Assuntos

RESUMO

Selecting the values of parameters used by de novo genomic assembly programs, or choosing an optimal de novo assembly from several runs obtained with different parameters or programs, are tasks that can require complex decision-making. A key parameter that must be supplied to typical next generation sequencing (NGS) assemblers is the k-mer length, i.e., the word size that determines which de Bruijn graph the program should map out and use. The topic of assembly selection criteria was recently revisited in the Assemblathon 2 study (Bradnam et al., 2013). Although no clear message was delivered with regard to optimal k-mer lengths, it was shown with examples that it is sometimes important to decide if one is most interested in optimizing the sequences of protein-coding genes (the gene space) or in optimizing the whole genome sequence including the intergenic DNA, as what is best for one criterion may not be best for the other. In the present study, our aim was to better understand how the assembly of unicellular fungi (which are typically intermediate in size and complexity between prokaryotes and metazoan eukaryotes) can change as one varies the k-mer values over a wide range. We used two different de novo assembly programs (SOAPdenovo2 and ABySS), and simple assembly metrics that also focused on success in assembling the gene space and repetitive elements. A recent increase in Illumina read length to around 150 bp allowed us to attempt de novo assemblies with a larger range of k-mers, up to 127 bp. We applied these methods to Illumina paired-end sequencing read sets of fungal strains of Paracoccidioides brasiliensis and other species. By visualizing the results in simple plots, we were able to track the effect of changing k-mer size and assembly program, and to demonstrate how such plots can readily reveal discontinuities or other unexpected characteristics that assembly programs can present in practice, especially when they are used in a traditional molecular microbiology laboratory with a 'genomics corner'. Here we propose and apply a component of a first pass validation methodology for benchmarking and understanding fungal genome de novo assembly processes.

Assuntos

RESUMO

In recent years, readily affordable short read sequences provided by next-generation sequencing (NGS) have become longer and more accurate. This has led to a jump in interest in the utility of NGS-only approaches for exploring eukaryotic genomes. The concept of a static, 'finished' genome assembly, which still appears to be a faraway goal for many eukaryotes, is yielding to new paradigms. We here motivate an object-view concept where the raw reads are the main, fixed object, and assemblies with their annotations take a role of dynamically changing and modifiable views of that object.

Assuntos