Evolution of the 12-Spanner Drug:H+ Antiporter DHA1 Family in Hemiascomycetous Yeasts

Transcriptome Analysis of mfs2-Defective Penicillium digitatum Mutant to Reveal Importance of Pdmfs2 in Developing Fungal Prochloraz Resistance

Demethylation inhibitors (DMIs), including prochloraz, are popular fungicides to control citrus postharvest pathogens such as Penicillium digitatum (green mold). However, many P. digitatum strains have developed prochloraz resistance, which decreases drug efficacy. Specific major facilitator superfamily (MFS) transporter gene mfs2, encoding drug-efflux pump protein MFS2, has been identified in P. digitatum strain F6 (PdF6) to confer fungal strain prochloraz resistance. However, except for the drug-efflux pump function of MFS2, other mechanisms relating to the Pdmfs2 are not fully clear. The present study reported a transcriptome investigation on the mfs2-defective P. digitatum strain. Comparing to the wild-type strain, the mfs2-defective strain showed 717 differentially expressed genes (DEGs) without prochloraz induction, and 1221 DEGs with prochloraz induction. The obtained DEGs included multiple isoforms of MFS transporter-encoding genes, ATP-binding cassette (ABC) transporter-encoding genes, and multidrug and toxic compound extrusion (MATE) family protein-encoding genes. Many of these putative drug-efflux pump protein-encoding genes had significantly lower transcript abundances in the mfs2-defective P. digitatum strain at prochloraz induction, as compared to the wild-type strain, including twenty-two MFS transporter-encoding genes (MFS1 to MFS22), two ABC transporter-encoding genes (ABC1 and ABC2), and three MATE protein-encoding genes (MATE1 to MATE3). The prochloraz induction on special drug-efflux pump protein genes in the wild-type strain was not observed in the mfs2-defective strain, including MFS21, MFS22, ABC2, MATE1, MATE2, and MATE3. On the other hand, the up-regulation of other drug-efflux pump protein genes in the mfs2-defective strain cannot recover the fungal prochloraz resistance, including MFS23, MFS26, MFS27, MFS31, MFS33, and ABC3 to ABC8. The functional enrichment of DEGs based on Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and euKaryotic Orthologous Groups (KOG) database resources suggested some essential contributors to the mfs2-relating prochloraz resistance, including ribosome biosynthesis-related genes, oxidative phosphorylation genes, steroid biosynthesis-related genes, fatty acid and lipid metabolism-related genes, and carbon- and nitrogen-metabolism-related genes. The results indicated that the MFS2 transporter might be involved in the regulation of multiple drug-efflux pump protein gene expressions and multiple metabolism-related gene expressions, thus playing an important role in developing P. digitatum prochloraz resistance.

RNA sequencing reveals metabolic and regulatory changes leading to more robust fermentation performance during short-term adaptation of Saccharomyces cerevisiae to lignocellulosic inhibitors

BACKGROUND

The limited tolerance of Saccharomyces cerevisiae to inhibitors is a major challenge in second-generation bioethanol production, and our understanding of the molecular mechanisms providing tolerance to inhibitor-rich lignocellulosic hydrolysates is incomplete. Short-term adaptation of the yeast in the presence of dilute hydrolysate can improve its robustness and productivity during subsequent fermentation.

RESULTS

We utilized RNA sequencing to investigate differential gene expression in the industrial yeast strain CR01 during short-term adaptation, mimicking industrial conditions for cell propagation. In this first transcriptomic study of short-term adaption of S. cerevisiae to lignocellulosic hydrolysate, we found that cultures respond by fine-tuned up- and down-regulation of a subset of general stress response genes. Furthermore, time-resolved RNA sequencing allowed for identification of genes that were differentially expressed at 2 or more sampling points, revealing the importance of oxidative stress response, thiamin and biotin biosynthesis. furan-aldehyde reductases and specific drug:H⁺ antiporters, as well as the down-regulation of certain transporter genes.

CONCLUSIONS

These findings provide a better understanding of the molecular mechanisms governing short-term adaptation of S. cerevisiae to lignocellulosic hydrolysate, and suggest new genetic targets for improving fermentation robustness.

Molecular insights into the differential efflux mechanism of Rv1634 protein, a multidrug transporter of major facilitator superfamily in Mycobacterium tuberculosis

Currently, multidrug-resistant tuberculosis (MDR-TB) is a public health crisis and a major health security threat globally. In Mycobacterium tuberculosis (Mtb), major facilitator superfamily (MFS) is the largest group of secondary active transporters. Along with the transport of their natural substrates, MFS proteins were involved in a drug efflux mechanism that ultimately lead to resistance against available anti-TB drugs in Mtb. In the present study, the three-dimensional structure model of an MFS protein, Rv1634, a probable multidrug transporter from Mtb, was generated using homology modeling. The protein structure model was found in inward-open conformation having 14 transmembrane helices. In addition, a central transport channel was deduced across the protein, and a single binding pocket was identified halfway through the central cavity by structural alignment with the homologous protein structures. Further, Rv1634 protein was studied based on the differential structural behavior of apo and ligand-bound forms. All the protein systems were inserted into a phospholipid bilayer to characterize the conformational dynamics of the protein using molecular dynamics (MD) simulations. Detailed analysis of the MD trajectories showed the diverse substrate specificity of the binding pocket for the antibiotics that caused differential movement in the ciprofloxacin and norfloxacin, to which Mtb strains have now become resistant. The expulsion of the drugs outside the bacterial cell occurs through the alternating-access mechanism of N and C-terminal domains, which is intriguing and essential to the understanding the drug resistance mechanism in pathogenic bacteria.

Physiological Genomics of Multistress Resistance in the Yeast Cell Model and Factory: Focus on MDR/MXR Transporters

The contemporary approach of physiological genomics is vital in providing the indispensable holistic understanding of the complexity of the molecular targets, signalling pathways and molecular mechanisms underlying the responses and tolerance to stress, a topic of paramount importance in biology and biotechnology. This chapter focuses on the toxicity and tolerance to relevant stresses in the cell factory and eukaryotic model yeast Saccharomyces cerevisiae. Emphasis is given to the function and regulation of multidrug/multixenobiotic resistance (MDR/MXR) transporters. Although these transporters have been considered drug/xenobiotic efflux pumps, the exact mechanism of their involvement in multistress resistance is still open to debate, as highlighted in this chapter. Given the conservation of transport mechanisms from S. cerevisiae to less accessible eukaryotes such as plants, this chapter also provides a proof of concept that validates the relevance of the exploitation of the experimental yeast model to uncover the function of novel MDR/MXR transporters in the plant model Arabidopsis thaliana. This knowledge can be explored for guiding the rational design of more robust yeast strains with improved performance for industrial biotechnology, for overcoming and controlling the deleterious activities of spoiling yeasts in the food industry, for developing efficient strategies to improve crop productivity in agricultural biotechnology.

The Paralogous Genes PDR18 and SNQ2, Encoding Multidrug Resistance ABC Transporters, Derive From a Recent Duplication Event, PDR18 Being Specific to the Saccharomyces Genus

Pleiotropic drug resistance (PDR) family of ATP-binding cassette (ABC) transporters play a key role in the simultaneous acquisition of resistance to a wide range of structurally and functionally unrelated cytotoxic compounds in yeasts. Saccharomyces cerevisiae Pdr18 was proposed to transport ergosterol at the plasma membrane, contributing to the maintenance of adequate ergosterol content and decreased levels of stress-induced membrane disorganization and permeabilization under multistress challenge leading to resistance to ethanol, acetic acid and the herbicide 2,4-D, among other compounds. PDR18 is a paralog of SNQ2, first described as a determinant of resistance to the chemical mutagen 4-NQO. The phylogenetic and neighborhood analysis performed in this work to reconstruct the evolutionary history of ScPDR18 gene in Saccharomycetaceae yeasts was focused on the 214 Pdr18/Snq2 homologs from the genomes of 117 strains belonging to 29 yeast species across that family. Results support the idea that a single duplication event occurring in the common ancestor of the Saccharomyces genus yeasts was at the origin of PDR18 and SNQ2, and that by chromosome translocation PDR18 gained a subtelomeric region location in chromosome XIV. The multidrug/multixenobiotic phenotypic profiles of S. cerevisiae pdr18Δ and snq2Δ deletion mutants were compared, as well as the susceptibility profile for Candida glabrata snq2Δ deletion mutant, given that this yeast species has diverged previously to the duplication event on the origin of PDR18 and SNQ2 genes and encode only one Pdr18/Snq2 homolog. Results show a significant overlap between ScSnq2 and CgSnq2 roles in multidrug/multixenobiotic resistance (MDR/MXR) as well as some overlap in azole resistance between ScPdr18 and CgSnq2. The fact that ScSnq2 and ScPdr18 confer resistance to different sets of chemical compounds with little overlapping is consistent with the subfunctionalization and neofunctionalization of these gene copies. The elucidation of the real biological role of ScSNQ2 will enlighten this issue. Remarkably, PDR18 is only found in Saccharomyces genus genomes and is present in almost all the recently available 1,000 deep coverage genomes of natural S. cerevisiae isolates, consistent with the relevant encoded physiological function.

Genome reprogramming in Saccharomyces cerevisiae upon nonylphenol exposure

Bioaccumulative environmental estrogen, nonylphenol (NP; 4-nonylphenol), is widely used as a nonionic surfactant and can affect human health. Since genomes of Saccharomyces cerevisiae and higher eukaryotes share many structural and functional similarities, we investigated subcellular effects of NP on S. cerevisiae BY4742 cells by analyzing genome-wide transcriptional profiles. We examined effects of low (1 mg/l; <15% cell number reduction) and high (5 mg/l; >65% cell number reduction) inhibitory concentration exposures for 120 or 180 min. After 120 and 180 min of 1 mg/l NP exposure, 187 (63 downregulated, 124 upregulated) and 103 genes (56 downregulated, 47 upregulated), respectively, were differentially expressed. Similarly, 678 (168 repressed, 510 induced) and 688 genes (215 repressed, 473 induced) were differentially expressed in cells exposed to 5 mg/l NP for 120 and 180 min, respectively. Only 15 downregulated and 63 upregulated genes were common between low and high NP inhibitory concentration exposure for 120 min, whereas 16 downregulated and 31 upregulated genes were common after the 180-min exposure. Several processes/pathways were prominently affected by either low or high inhibitory concentration exposure, while certain processes were affected by both inhibitory concentrations, including ion transport, response to chemicals, transmembrane transport, cellular amino acids, and carbohydrate metabolism. While minimal expression changes were observed with low inhibitory concentration exposure, 5 mg/l NP treatment induced substantial expression changes in genes involved in oxidative phosphorylation, cell wall biogenesis, ribosomal biogenesis, and RNA processing, and encoding heat shock proteins and ubiquitin-conjugating enzymes. Collectively, these results provide considerable information on effects of NP at the molecular level.

The Zygosaccharomyces bailii transcription factor Haa1 is required for acetic acid and copper stress responses suggesting subfunctionalization of the ancestral bifunctional protein Haa1/Cup2

Background

The food spoilage yeast species Zygosaccharomyces bailii exhibits an extraordinary capacity to tolerate weak acids, in particular acetic acid. In Saccharomyces cerevisiae, the transcription factor Haa1 (ScHaa1) is considered the main player in genomic expression reprogramming in response to acetic acid stress, but the role of its homologue in Z. bailii (ZbHaa1) is unknown.

Results

In this study it is demonstrated that ZbHaa1 is a ScHaa1 functional homologue by rescuing the acetic acid susceptibility phenotype of S. cerevisiae haa1Δ. The disruption of ZbHAA1 in Z. bailii IST302 and the expression of an extra ZbHAA1 copy confirmed ZbHAA1 as a determinant of acetic acid tolerance. ZbHaa1 was found to be required for acetic acid stress-induced transcriptional activation of Z. bailii genes homologous to ScHaa1-target genes. An evolutionary analysis of the Haa1 homologues identified in 28 Saccharomycetaceae species genome sequences, including Z bailii, was carried out using phylogenetic and gene neighbourhood approaches. Consistent with previous studies, this analysis revealed a group containing pre-whole genome duplication species Haa1/Cup2 single orthologues, including ZbHaa1, and two groups containing either Haa1 or Cup2 orthologues from post-whole genome duplication species. S. cerevisiae Cup2 (alias Ace1) is a transcription factor involved in response and tolerance to copper stress. Taken together, these observations led us to hypothesize and demonstrate that ZbHaa1 is also involved in copper-induced transcriptional regulation and copper tolerance.

Conclusions

The transcription factor ZbHaa1 is required for adaptive response and tolerance to both acetic acid and copper stresses. The subfunctionalization of the single ancestral Haa1/Cup2 orthologue that originated Haa1 and Cup2 paralogues after whole genome duplication is proposed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3443-2) contains supplementary material, which is available to authorized users.

Transcriptomic responses of the basidiomycete yeast Sporobolomyces sp. to the mycotoxin patulin

BACKGROUND

Patulin is a mycotoxin produced by Penicillium expansum, the causal agent of blue mold of stored pome fruits, and several other species of filamentous fungi. This mycotoxin has genotoxic, teratogenic and immunotoxic effects in mammals, and its presence in pome fruits and derived products represents a serious health hazard. Biocontrol agents in the Pucciniomycotina, such as the yeasts Sporobolomyces sp. strain IAM 13481 and Rhodosporidium kratochvilovae strain LS11, are able to resist patulin and degrade it into the less toxic compounds desoxypatulinic acid and ascladiol.

RESULTS

In this investigation we applied a transcriptomic approach based on RNAseq to annotate the genome of Sporobolomyces sp. IAM 13481 and then study the changes of gene expression in Sporobolomyces sp. exposed to patulin. Patulin treatment leads to ROS production and oxidative stress that result in the activation of stress response mechanisms controlled by transcription factors. Upregulated Sporobolomyces genes were those involved in oxidation-reduction and transport processes, suggesting the activation of defense mechanisms to resist patulin toxicity and expel the mycotoxin out of the cells. Other upregulated genes encoded proteins involved in metabolic processes such as those of the glutathione and thioredoxin systems, which are essential to restore the cellular redox homeostasis. Conversely, patulin treatment decreased the expression of genes involved in the processes of protein synthesis and modification, such as transcription, RNA processing, translation, protein phosphorylation and biosynthesis of amino acids. Also, genes encoding proteins involved in transport of ions, cell division and cell cycle were downregulated. This indicates a reduction of metabolic activity, probably due to the high energy requirement by the cells or metabolic arrest while recovering from the insult caused by patulin toxicity.

CONCLUSIONS

Complex mechanisms are activated in a biocontrol yeast in response to patulin. The genes identified in this study can pave the way to develop i) a biodetoxification process of patulin in juices and ii) a biosensor for the rapid and cost-effective detection of this mycotoxin.

Phylogenetic and syntenic analyses of the 12-spanner drug:H(+) antiporter family 1 (DHA1) in pathogenic Candida species: evolution of MDR1 and FLU1 genes

Candida albicans and other pathogenic Candida species can develop resistance to clinical fungicides through active drug export mediated by multidrug efflux pumps, in particular by members of the drug:H(+) antiporter family 1 (DHA1). The DHA1 proteins encoded in the genomes of 31 hemiascomycetous strains from 25 species were identified and homology relationships between these proteins and the functionally characterised DHA1 in the model yeast Saccharomyces cerevisiae were established. Gene neighbourhood analysis allowed the reconstruction of sixteen DHA1 lineages conserved during the CTG complex species evolution. The evolutionary history of C. albicans MDR1 and FLU1 genes and Candida dubliniensis, Candida tropicalis and Candida parapsilosis MDR1 genes was detailed. Candida genomes show an abundant number of MDR1 and FLU1 homologues but the chromosome environment where MDR1 homologues reside was poorly conserved during evolution. Gene duplication and loss are major mechanisms underlying the evolution of the DHA1 genes in Candida species.

MFS multidrug transporters in pathogenic fungi: do they have real clinical impact?

Infections caused by opportunistic fungal pathogens have reached concerning numbers due to the increase of the immunocrompromised human population and to the development of antifungal resistance. This resistance is often attributed to the action of multidrug efflux pumps, belonging to the ATP-binding cassette (ABC) superfamily and the major facilitator superfamily (MFS). Although many studies have focused on the role of ABC multidrug efflux transporters, little is still known on the part played by the Drug:H⁺ Antiporter (DHA) family of the MFS in this context. This review summarizes current knowledge on the role in antifungal drug resistance, mode of action and phylogenetic relations of DHA transporters, from the model yeast S. cerevisiae to pathogenic yeasts and filamentous fungi. Through the compilation of the predicted DHA transporters in the medically relevant Candida albicans, C. glabrata, C. parapsilosis, C. lusitaniae, C. tropicalis, C. guilliermondii, Cryptococcus neoformans, and Aspergillus fumigatus species, the fact that only 5% of the DHA transporters from these organisms have been characterized so far is evidenced. The role of these transporters in antifungal drug resistance and in pathogen-host interaction is described and their clinical relevance discussed. Given the knowledge gathered for these few DHA transporters, the need to carry out a systematic characterization of the DHA multidrug efflux pumps in fungal pathogens, with emphasis on their clinical relevance, is highlighted.

MFS transporters required for multidrug/multixenobiotic (MD/MX) resistance in the model yeast: understanding their physiological function through post-genomic approaches

Multidrug/Multixenobiotic resistance (MDR/MXR) is a widespread phenomenon with clinical, agricultural and biotechnological implications, where MDR/MXR transporters that are presumably able to catalyze the efflux of multiple cytotoxic compounds play a key role in the acquisition of resistance. However, although these proteins have been traditionally considered drug exporters, the physiological function of MDR/MXR transporters and the exact mechanism of their involvement in resistance to cytotoxic compounds are still open to debate. In fact, the wide range of structurally and functionally unrelated substrates that these transporters are presumably able to export has puzzled researchers for years. The discussion has now shifted toward the possibility of at least some MDR/MXR transporters exerting their effect as the result of a natural physiological role in the cell, rather than through the direct export of cytotoxic compounds, while the hypothesis that MDR/MXR transporters may have evolved in nature for other purposes than conferring chemoprotection has been gaining momentum in recent years. This review focuses on the drug transporters of the Major Facilitator Superfamily (MFS; drug:H⁺ antiporters) in the model yeast Saccharomyces cerevisiae. New insights into the natural roles of these transporters are described and discussed, focusing on the knowledge obtained or suggested by post-genomic research. The new information reviewed here provides clues into the unexpectedly complex roles of these transporters, including a proposed indirect regulation of the stress response machinery and control of membrane potential and/or internal pH, with a special emphasis on a genome-wide view of the regulation and evolution of MDR/MXR-MFS transporters.

The drug:H⁺ antiporters of family 2 (DHA2), siderophore transporters (ARN) and glutathione:H⁺ antiporters (GEX) have a common evolutionary origin in hemiascomycete yeasts

BACKGROUND

The Saccharomyces cerevisiae 14-spanner Drug:H+ Antiporter family 2 (DHA2) are transporters of the Major Facilitator Superfamily (MFS) involved in multidrug resistance (MDR). Although poorly characterized, DHA2 family members were found to participate in the export of structurally and functionally unrelated compounds or in the uptake of amino acids into the vacuole or the cell. In S. cerevisiae, the four ARN/SIT family members encode siderophore transporters and the two GEX family members encode glutathione extrusion pumps. The evolutionary history of DHA2, ARN and GEX genes, encoding 14-spanner MFS transporters, is reconstructed in this study.

RESULTS

The translated ORFs of 31 strains from 25 hemiascomycetous species, including 10 pathogenic Candida species, were compared using a local sequence similarity algorithm. The constraining and traversing of a network representing the pairwise similarity data gathered 355 full size proteins and retrieved ARN and GEX family members together with DHA2 transporters, suggesting the existence of a close phylogenetic relationship among these 14-spanner major facilitators. Gene neighbourhood analysis was combined with tree construction methodologies to reconstruct their evolutionary history and 7 DHA2 gene lineages, 5 ARN gene lineages, and 1 GEX gene lineage, were identified. The S. cerevisiae DHA2 proteins Sge1, Azr1, Vba3 and Vba5 co-clustered in a large phylogenetic branch, the ATR1 and YMR279C genes were proposed to be paralogs formed during the Whole Genome Duplication (WGD) whereas the closely related ORF YOR378W resides in its own lineage. Homologs of S. cerevisiae DHA2 vacuolar proteins Vba1, Vba2 and Vba4 occur widespread in the Hemiascomycetes. Arn1/Arn2 homologs were only found in species belonging to the Saccharomyces complex and are more abundant in the pre-WGD species. Arn4 homologs were only found in sub-telomeric regions of species belonging to the Sacharomyces sensu strictu group (SSSG). Arn3 type siderophore transporters are abundant in the Hemiascomycetes and form an ancient gene lineage extending to the filamentous fungi.

CONCLUSIONS

The evolutionary history of DHA2, ARN and GEX genes was reconstructed and a common evolutionary root shared by the encoded proteins is hypothesized. A new protein family, denominated DAG, is proposed to span these three phylogenetic subfamilies of 14-spanner MFS transporters.

Structural advances for the major facilitator superfamily (MFS) transporters

The major facilitator superfamily (MFS) is one of the largest groups of secondary active transporters conserved from bacteria to humans. MFS proteins selectively transport a wide spectrum of substrates across biomembranes and play a pivotal role in multiple physiological processes. Despite intense investigation, only seven MFS proteins from six subfamilies have been structurally elucidated. These structures were captured in distinct states during a transport cycle involving alternating access to binding sites from either side of the membrane. This review discusses recent progress in MFS structure analysis and focuses on the molecular basis for substrate binding, co-transport coupling, and alternating access.

Candida albicans flu1-mediated efflux of salivary histatin 5 reduces its cytosolic concentration and fungicidal activity

Histatin 5 (Hst 5) is a salivary human antimicrobial peptide that is toxic to the opportunistic yeast Candida albicans. Fungicidal activity of Hst 5 requires intracellular translocation and accumulation to a threshold concentration for it to disrupt cellular processes. Previously, we observed that total cytosolic levels of Hst 5 were gradually reduced from intact cells, suggesting that C. albicans possesses a transport mechanism for efflux of Hst 5. Since we identified C. albicans polyamine transporters responsible for Hst 5 uptake, we hypothesized that one or more polyamine efflux transporters may be involved in the efflux of Hst 5. C. albicans FLU1 and TPO2 were found to be the closest homologs of Saccharomyces cerevisiae TPO1, which encodes a major spermidine efflux transporter, indicating that the products of these two genes may be involved in efflux of Hst 5. We found that flu1Δ/Δ cells, but not tpo2Δ/Δ cells, had significant reductions in their rates of Hst 5 efflux and had significantly higher cytoplasmic Hst 5 and Hst 5 susceptibilities than did the wild type. We also found that flu1Δ/Δ cells had reduced biofilm formation compared to wild-type cells in the presence of Hst 5. Transcriptional levels of FLU1 were not altered over the course of treatment with Hst 5; therefore, Hst 5 is not likely to induce FLU1 gene overexpression as a potential mechanism of resistance. Thus, Flu1, but not Tpo2, mediates efflux of Hst 5 and is responsible for reduction of its toxicity in C. albicans.

IONS: Identification of Orthologs by Neighborhood and Similarity-an Automated Method to Identify Orthologs in Chromosomal Regions of Common Evolutionary Ancestry and its Application to Hemiascomycetous Yeasts

Comparative sequence analysis is widely used to infer gene function and study genome evolution and requires proper ortholog identification across different genomes. We have developed a program for the Identification of Orthologs in one-to-one relationship by Neighborhood and Similarity (IONS) between closely related species. The algorithm combines two levels of evidence to determine co-ancestrality at the genome scale: sequence similarity and shared neighborhood. The method was initially designed to provide anchor points for syntenic blocks within the Génolevures project concerning nine hemiascomycetous yeasts (about 50,000 genes) and is applicable to different input databases. Comparison based on use of a Rand index shows that the results are highly consistent with the pillars of the Yeast Gene Order Browser, a manually curated database. Compared with SYNERGY, another algorithm reporting homology relationships, our method’s main advantages are its automation and the absence of dataset-dependent parameters, facilitating consistent integration of newly released genomes.