Thuyen Truong1
Tanujaa Suriyanarayanan1
Guisheng Zeng2Thuc D. Le3Lin Liu3Jiuyong Li3Cao Tong1
Yue Wang2,4*
Chaminda J. Seneviratne1- 1Oral Sciences, Faculty of Dentistry, National University of Singapore, Singapore, Singapore
- 2Institute of Molecular and Cell Biology, Agency for Science, Technology and Research, Singapore, Singapore
- 3School of Information Technology and Mathematical Sciences, University of South Australia, Adelaide, SA, Australia
- 4Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore
Antifungal agents for the treatment of Candida albicans infections are limited. We recently discovered a novel antifungal small molecule, SM21, with promising in vivo activity. Herein, we employed the newly developed C. albicans haploid toolbox to uncover the mechanism of action of SM21. Comprehensive RNA-Seq analyses of the haploid susceptible GZY803 strain revealed significant gene expression changes related to mitochondria when exposed to SM21. Mitochondrial structure visualization and measurement of ATP generation, reactive oxygen species (ROS) levels, and the antioxidant potential of SM21-treated and untreated GZY803, mitochondrial structure defective haploid mutant (dnm1Δ), and wild-type diploid SC5314 strains confirmed defects in mitochondria. Exploiting the advantage of C. albicans haploids as a single ploidy model, we further exposed GZY803 to repetitive treatments of SM21 in order to generate resistant mutants. Three colonies designated S3, S5 and S6, which displayed resistance to SM21, were isolated. All resistant strains exhibited enhanced transcriptomic responses for peptide and protein metabolism and secreted aspartate proteases (SAPs) activity under SM21 treatment compared to the parent strain GZY803. Consistently, supplementing the resistant strains, GZY803, and SC5314 with peptone, a form of digested peptides, decreased susceptibility to SM21. The present study demonstrates the usefulness of haploid C. albicans model in antifungal drug discovery. The findings will be invaluable to develop SM21 as a novel antifungal agent, which will benefit millions of patients suffering from Candida infections.
Introduction
Candida albicans is the most prevalent fungal pathogen in humans and ranked as the fourth leading cause of nosocomial blood-stream infections (Ashman et al., 2004). Systemic candidiasis is reported to have a mortality rate as high as 40% (Wisplinghoff et al., 2004). Moreover, C. albicans also causes superficial mucosal infections such as oral and vaginal candidiasis which are highly prevalent and persistent (Webb et al., 1998; Williams and Lewis, 2011). However, currently available antifungals for clinical uses are limited. Hence, there has been an urgent need to find novel anti-Candida agents for clinical uses. Our group recently discovered a new antifungal small molecule, designated SM21, which exhibited promising activity in vitro as well as in vivo (Wong et al., 2014). However, the mechanism by which SM21 exerts its antifungal activity remained elusive.
Until recently, C. albicans was thought to be an obligate diploid organism (Jones et al., 2004; Noble and Johnson, 2007). The invariably diploid and heterozygous nature of the organism was a considerable barrier for genetic studies in identifying novel gene functions and novel drug targets (Hickman et al., 2013). In addition, diploid model of C. albicans also poses a challenge for developing resistant mutations, which is a common approach to studying antifungal resistance of the organism (Woodford and Ellington, 2007). Generating resistant mutations in diploid C. albicans is difficult, as simultaneous inactivation of both copies of a gene is often required to produce a resistant phenotype, which is an extremely rare event. For instance, it takes about 100 days on average to generate a mutation resistant to echinocandin under optimized environmental conditions (Shields et al., 2015).
Existence of haploid forms of C. albicans was discovered in 2013 by Hickman et al. (2013). Single genome of haploid C. albicans could provide numerous advantages in biological and drug discovery research. Subsequently, our group demonstrated the utility of haploid model of C. albicans in uncovering novel molecular regulators in C. albicans biofilm formation and antifungal drug resistance (Seneviratne et al., 2015; Truong et al., 2016). Some haploid C. albicans strains were found to be able to maintain the haploidy in planktonic as well as biofilm cultures (Hickman et al., 2013; Seneviratne et al., 2015). In addition, we also demonstrated that C. albicans haploids contain a highly similar proteome to that of the diploids (Truong et al., 2016). Haploid cells have been shown to be more efficient in acquiring beneficial mutations at high doses of drug treatment when compared to that of the diploids in the yeast Saccharomyces cerevisiae (Anderson et al., 2004). Hence, we hypothesized that C. albicans haploid model could be a very useful tool to uncover drug targets of novel antifungal agents and resistance mechanisms.
Taking foregoing advantages into account, in the present study, we utilized haploid model of C. albicans to study the mechanism of action of the antifungal molecule SM21 and the resistance strategies of C. albicans against it. Firstly, transcriptomic analysis of the haploid C. albicans parent strain GZY803 in response to SM21 revealed diminished gene expression in mitochondrial ATP-coupled electron transport, respiration, and enhanced cytochrome assembly factors. Downstream analysis confirmed that SM21 treatment results in mitochondrial damage and inhibition of mitochondrial activity. Secondly, we generated C. albicans SM21-resistant haploid strains and performed transcriptomic analysis on the cultures with or without SM21 treatments. RNA-Seq analysis coupled with downstream experiments suggested the association between amino acid supplement and the resistance responses of C. albicans against SM21. Hence, the present study demonstrates a new application of the C. albicans haploid model in examining the drug targets of a novel antifungal and predicting possible resistance strategy by which the fungal pathogen may develop against it. The findings will aid future development of SM21 as a novel antifungal agent against C. albicans infections.
Materials and Methods
Strains, Plasmids, and Culture Conditions
Generation of C. albicans haploid strain (GZY803) was described previously (Hickman et al., 2013). Construction of other yeast strains and plasmids used in this study are described in Tables S3, S4, respectively. Unless specified, C. albicans cells were grown at 30°C in YPD (2% yeast extract, 1% peptone, and 2% glucose), or GMM (glucose minimal medium, 6.79 g/l yeast nitrogen base without amino acids, and 2% glucose) supplemented with appropriate amino acids (uridine 80 μg/ml, arginine 40 μg/ml, and histidine 40 μg/ml) and 5-FOA (1 mg/ml) when required. Solid culture plates were obtained by adding 2% agar into the media.
Antifungal Agents
SM21 (ChemBridge) Caspofungin, Amphotericin B, Fluconazole, Ketoconazole, and Voriconazole (Sigma) were used in this study.
Antifungal Susceptibility Testing for C. albicans Cells
Antifungal susceptibility tests of C. albicans cells was performed as described previously (Truong et al., 2016) following the Clinical and Laboratory Standards Institute's (CLSI) guideline. In brief, C. albicans cell suspension of 107 cells/ml (equivalent to 0.38 McFarLand standards) was first prepared from overnight culture at 30°C. Subsequently, the culture was diluted with RPMI to yield an inoculum of ~103 cells/ml. The MIC was determined in 96-well plates and each strain was exposed to the 2-fold diluted solutions of SM21. The plates were incubated at 30–35°C for 48 h before MIC values were determined following the CLSI guideline.
Generation of SM21-Resistant Haploid Strain
To induce spontaneous resistance to SM21 in haploid C. albicans strains, an inoculum of 108 cells/ml of GZY803 was cultured in 2 μg/ml of SM21 in GMM for 3 consecutive days in 30°C. Cultures were replenished with fresh SM21 in GMM every 24 h. Subsequently, cells were pelleted and spread on GMM plates containing 4 μg/ml of SM21 and incubated at 30°C until colonies developed. All colonies developed were selected and examined for their MIC90 to SM21. Three colonies, labeled as S3, S5, and S6, which display an SM21 MIC of 4 μg/ml, were selected for further studies.
RNA Extraction and Sequencing
RNA was extracted from GZY803 and the SM21-resistant strains, S3, S5, and S6, untreated or treated with 2 μg/ml of SM21 in GMM supplemented with required amino acids for 1 or 4 h at 30°C. Experiments were performed in three biological replicates for each strain in each condition. The RNA extraction was conducted following Ribopure RNA purification protocol (Qiagen); and samples were sent to Novogene company for quality check and RNA sequencing using Illumina HiSeq 4000 platform. Preliminary analysis was performed by Novogene.
Ploidy Analysis by Flow Cytometry
Ploidy of C. albicans strains was analyzed by flow cytometry as described (Zeng et al., 2014). In brief, C. albicans cells from 30°C cultures were washed with PBS, and fixed with 70% ethanol at 4°C overnight. Cells were then treated with RNase A (5 mg/ml) at 37°C for 6–7 h or overnight, washed with PBS, and further incubated with propidium iodine (5 mg/ml) in dark at 4°C overnight. Subsequently, stained samples were sonicated once to separate clustered cells, and subjected to DNA content analysis using BD FACSCalibur (BD Biosciences). Data acquisition was performed with 10,000 cells for each sample by using the CellQuest Pro software. The acquired data were then analyzed using WinMDI (version 2.8) software. Diploid (SC5314) and haploid (GZY803) cells were used as the control.
ATP Measurement
ATP generation was measured using BacTiter-Glo™ Microbial Cell Viability Assay (Promega, Singapore). Briefly, 107 cells of mid-log phrase cultures after treatment with 2 μg/ml of SM21 for 4 h were measured for ATP generation. Untreated cells were included as a control. An aliquot of treated and untreated cells were subjected to colony forming unit (CFU) assay as described before (Truong et al., 2016) to normalize for total cells count in the cultures.
Total Antioxidant Potential Measurement
Total antioxidant potential estimation was described previously (Truong et al., 2016). C. albicans cultures at mid-log phrase (107 cells/ml) were homogenized using glass beads (0.5 mm) and a high speed vortex (Vortexgene) following 3 cycles of 3 min on and 2 min off ice. Total antioxidant capacity assay kit (Abcam) was used to estimate the antioxidant potentials of the cultures following manufacture protocol.
Florescent Microscopy Analysis
Protocol for visualization of mitochondrial structures was adapted from previous publication (Mavrianos et al., 2013). Briefly, an inoculum of 107 C. albicans cells per ml were treated with 25 nM of MitoTracker Red CMXRos or MitoSOX™ Red Mitochondrial Superoxide Indicator (Invitrogen) for 20 min and the cells were visualized under an Olympus Fluoview FV1000 TIRF confocal microscope. Images were analyzed using Olympus FV10-ASW Viewer software.
Differential Expression Analysis
For each strain, genes which showed a positive log2-fold change in both 1 h and 4 h SM21 treated conditions and with an average of larger than or equal to 1 in the two conditions as compared to the non-treated control were labeled as “higher.” Meanwhile, genes with negative log2-fold changes with an average fold change less than or equal to −1 were classified as “lower.” On the other hand, to identify genes which were commonly and similarly regulated in the resistant conditions but were differently expressed in the susceptible strain, we applied differential expression Limma method (Bioconductor) (Ritchie et al., 2015). Among the significantly altered genes with a p < 0.05, those which showed an average log2 fold change of larger than or equal to 1 were classified as “higher;” while those with an log2 average fold change of less than or equal to −1, “lower,” in the resistant strains as compared to the wild type.
Cluster Display and Gene Ontology Analysis
Cluster display of C. albicans' gene expression was generated using Cluster (version 3.0) and TreeView (version 1.1.6) software (Eisen et al., 1998). Biological processes and cellular pathways which were significantly altered in the samples were illustrated using Cytoscape (version 2.8.3) (Shannon et al., 2003) with BINGO plugin (version 2.44) (Maere et al., 2005).
Statistical Analysis
Statistical analysis was performed using Student t-test. Analysis was done using SPSS software (Version 16.0, SPSS Inc.). The level of significance was set at p < 0.05.
Results
Transcriptomics Analysis of Haploid C. albicans Exposed to SM21 Treatment
SM21 was shown to be an inhibitor of the yeast-to-hypha transition of C. albicans, having no obvious effect on bacteria. SM21 also exhibited low toxicity to human cells (Wong et al., 2014). Hence, it is likely that the mechanism of action of SM21 is through fungal-specific targets or biological pathways. In order to understand the early lethal changes in the fungal molecular pathways induced by SM21, the C. albicans haploid parent strain (GZY803) was treated with 2 μg/ml of SM21, equivalent to 2 × MIC90, and RNA samples were harvested before the exposure and at 1 and 4 h. RNA was extracted from these test cultures as well as controls without SM21 treatment. Transcriptomics analysis was performed by RNA-Seq (Wang et al., 2009). Following differential expression analysis, 91 genes were classified as “higher” and 85 as “lower” in GZY803 strain under SM21 exposure compared to the untreated control. Figure 1A illustrates the expression pattern of all detected genes in a two-dimensional hierarchical clustering. Each gene was colored according to its log fold change. Next, we mapped the differently regulated proteins to their gene ontology. The analysis revealed that the expression of genes involved in response to temperature stimuli, drug export, and cytochrome assembly was significantly higher (Figure 1B) while the expression of those associated with ATP-coupled electron transport and transmembrane transport, nucleotide metabolism, cellular respiration, and ion homeostasis was significantly lower (Figure 1C) in SM21-treated GZY803 cultures compared to the untreated controls. A full list of differently expressed RNA is summarized in Table S1.
Figure 1. Cluster display and functional annotation of differentially expressed genes in GZY803 cultures when exposed to SM21. Overnight cultures of C. albicans haploid strain GZY803 at mid-log phrase (108 cell/ml) were treated with SM21 (2 μg/m) for 1 or 4 h and RNA was extracted for sequencing analysis. (A) Clustered display of protein expression profiles of 1 or 4 h SM21-treated GZY803 strains. Biological processes, molecular functions, and cellular components that are significantly and consistently higher (B) or lower (C) in SM21-treated as compared to the untreated GZY803 strain.
Effects of SM21 on the Mitochondria of C. albicans Haploids
Transcriptomics studies above revealed that SM21 significantly alters the expression of genes associated with mitochondrial respiration, ATP synthesis, cytochrome assembly in haploid C. albicans (Figure 1, Table 1). The data suggest significant dysfunction of mitochondrial components when treated with SM21. Therefore, further downstream experiments were conducted focusing on this aspect. Firstly, mitochondrial structural changes in haploid C. albicans GZY803 and dnm1Δ strains when exposed to SM21 were examined. Dnm1 is reported to regulate mitochondrial structures in S. cerevisiae, and deletion of DMN1 results in collapsed mitochondrial network with fewer branches compared to the wild type (Bleazard et al., 1999). Hence, dnm1Δ was included as the positive control with defective mitochondrial structures. First, dnm1Δ was found to be more sensitive to SM21 with a MIC90 of 0.5 μg/ml as compared to the parent GZY803 that had a MIC90 of 1 μg/ml (Table 2). Furthermore, confocal visualization of the cultures stained with MitoTracker, a dye for mitochondria, indicated that DMN1 deletion results in compromised mitochondria with fewer branches than the wild-type GZY803. Under SM21 treatment, both GZY803 and dnm1Δ exhibited collapsed mitochondrial membrane potentials providing more evidence that SM21 targets mitochondria of C. albicans (Figure 2A). Consistently, ATP production was significantly reduced in both strains upon SM21 exposure (p < 0.05) with GZY803 showed a 1.5-fold reduction, and dnm1Δ, 3-fold, under SM21 exposure (Figure 2B). Reactive oxygen species (ROS) are the byproducts of cellular metabolism that are mostly generated at the mitochondrial level (Murphy, 2009). To examine the effect of SM21 on ROS generation and release, we stained GZY803 and dnm1Δ cells with MitoSOX, a florescent probe for ROS. Indeed, under SM21 treatment, ROS were released in both dnm1Δ and GZY803 (Figure 2C). Total antioxidant capacity measurements corroborated the observations, in which antioxidant potentials were significantly lower in dnm1Δ and GZY803 under SM21 treatment (p < 0.05) (Figure 2D). The data clearly indicate that fungal mitochondria are likely to be the target of the fungicidal effect exerted by SM21 on haploid C. albicans.
Table 1. Changes in gene expression related to complexes in electron transport chain and assembly factors in GZY803 under SM21 exposure.
Table 2. Minimum inhibition concentration of Candida albicans strains.
Figure 2. Examination of mitochondrial structures and defects in SM21-treated haploid cultures. C. albicans haploid GZY803 and dnm1Δ strains were used in the experiments. Mid-log phrase cultures (107 cell/ml) were left untreated or treated with 2 μg/ml of SM21 for 4 h and cells were collected for subsequent confocal microscope visualization and ATP or antioxidant potential measurement. MitoTracker (A) was used to visualize the mitochondrial membrane structure under no treatment or with SM21 treatment. ATP production measurement (B), ROS generation visualization via mitoSOX staining (C) and total antioxidant potential estimation (D) of treated and untreated GZY803 and dnm1Δ strains were also performed. Scale bar, 5 μm. For all graphs, the mean of at least three replicates is shown, with error bars showing SD. *p < 0.05.
Confirmation of SM21 Inhibitory Effects on Mitochondria of C. albicans Diploid Cells
In order to confirm mitochondrial defects in C. albicans diploids, we examined the changes of mitochondrial functions of a standard diploid strain, SC5314, under SM21 treatment. Consistent with the above findings, SC5314 exhibited collapsed mitochondrial membrane potentials when treated with SM21 but not with amphotericin B (Figure 3A). Under SM21 treatment, SC5314 also displayed a reduction in ATP generation (p < 0.05) (Figure 3B), an increase in ROS generation (Figure 3C), and a corresponding decrease in antioxidant potential (p < 0.05) (Figure 3D). These data, therefore, confirmed that SM21 inhibits mitochondria activity in C. albicans.
Figure 3. Changes in mitochondrial structures and functions of C. albicans diploids under SM21 exposure. C. albicans diploid strain, SC5314, at mid-log phrase (107 cell/ml) was used in the experiments. Cultures were left untreated or treated with 2 μg/ml of SM21 for 4 h and cells were collected for subsequent visualization of mitochondrial structure (A), ATP measurement (B), ROS release visualization (C), and antioxidant potential examination (D). Treatment of SC5314 cultures with amphotericin B (2 μg/ml, 4 h) was also included in confocal microscope visualization to serve as the control. Scale bar, 5 μm. For all graphs, the mean of at least three replicates is shown, with error bars showing SD. *p < 0.05.
Generation of SM21-Resistant Haploid Strains
In order to understand the drug-pathogen interaction, it is important to know the mechanisms of action of the antimicrobial as well as the resistance mechanisms expressed by the pathogen against it. Therefore, as the next step, taking advantage of the haploid C. albicans model, we generated SM21-resistant mutants from the standard haploid parent GZY803. An inoculum of 108 cells of GZY803 were grown in GMM broth supplemented with 2 μg/ml of SM21 for 3 consecutive days with the SM21-containing medium refreshed every 24 h; surviving cells were then spread on Sabouraud dextrose agar plates containing 4 μg/ml of SM21. Among the colonies developed, three colonies designated as S3, S5, and S6, displayed a SM21 MIC of 4 μg/ml and were selected for further studies. Other colonies having a MIC < 4 μg/ml were discarded. Interestingly, all three colonies resistant to SM21 exhibited similar susceptibility to other antifungals tested as the wild-type strain (Table 2). This suggested that the resistance mechanisms possessed by the resistant haploid C. albicans colonies are effective against SM21 but not to the other antifungals tested. Susceptibility of the resistant strains to SM21 as well as other antifungals were repeatedly examined in several different occasions and in parallel with any other experiments performed to ensure that resistances acquired in the strains are non-reversal. In addition, FACS ploidy examination was also performed which revealed that all of the resistant strains maintained the haploidy (Figure 4). Therefore, using the haploid model of C. albicans is a suitable strategy to uncover resistance mechanisms.
Figure 4. Ploidy examination of resistant haploid strains of C. albicans. SM21-resistant haploid cells were grown at 30°C overnight. Subsequently, mid-log phrase cultures were taken for ploidy examination by flow cytometry analysis. SC5314 and GZY803 were used as the standard for diploid and haploid genomes.
Transcriptomic Analysis on SM21-Resistant Haploid Strains
In order to obtain clues on the resistance mechanism, the resistant haploid strains were treated with SM21 together with untreated controls to obtain RNA for RNA-Seq analysis. Subsequent transcriptomic examination involved (1) analysis of genes differently expressed in each resistant haploid strain under control and SM21 treatment; and (2) determination of genes which showed a similar expression trends in all resistant strains but were differently expressed compared to that in GZY803 under the same conditions. First step of differential analyses found that genes related to mitochondrial organization, cytochrome assembly, and protein and peptide metabolism were upregulated (Figure 5), while genes involved in cellular respiration, and ATP-coupled electron and transmembrane transports, downregulated, in all resistant strains under SM21 treatment (Figure 6). This suggested that all the resistant strains were affected by SM21 inhibitory effects on mitochondrial activity. Among the strains, under SM21 treatment, S3 and S5 displayed lower expression of genes involved in ion homeostasis, and S6, aspartate and methionine biosynthesis (Figure 6). Next, we applied Limma differential gene expression analysis to identify genes commonly expressed in all the resistant strains, which, on the other hand, were differently expressed in GZY803 under the same conditions (Ritchie et al., 2015). The method identified 211 genes whose expressions were higher, and 122 genes, lower, in the resistant haploid strains as compared to GZY803 (Table S2). Figure 7A illustrated clustering of differently expressed genes in the SM21 resistant haploid strains as compared to GZY803. Functional annotation of these genes revealed that expressions of genes modulating ligase activity, aspartate endopeptidases, protein and peptide metabolism, symbiont responses, mitochondrion, and ribosome structures were higher, and expression of genes involved in biotin biosynthesis, organonitrogen catabolism, transmembrane transport, cell periphery, and plasma membrane was lower in the resistant strains than in GZY803 in un-treated or SM21 exposed conditions (Figures 7B,C).
Figure 5. Gene ontology analysis of genes with higher expression in resistant strains when exposed to SM21. SM21-resistant haploid strains at mid-log phrase (108 cell/ml) were treated with 2 μg/ml of SM21 for 1 or 4 h. RNA was extracted and sent for RNA sequencing. Differential expression analysis was first applied to determined genes which were upregulated in the resistant haploid strains under SM21 exposure. Subsequent gene ontology analysis illustrated biological processes which were significantly higher in the resistant haploid strains under SM21 treatment.
Figure 6. Gene ontology analysis of genes with lower expression in resistant strains when exposed to SM21. SM21-resistant haploid strains at mid-log phrase (108 cell/ml) were treated with 2 μg/ml of SM21 for 1 or 4 h. RNA was extracted and sent for RNA sequencing. Differential expression analysis was first applied to determine genes which were downregulated in the resistant haploid strains under SM21 exposure. Subsequent gene ontology analysis illustrated biological processes which were significantly lower in the resistant haploid strains under SM21 stimulant.
Figure 7. Cluster display and functional analysis of gene expression changes in resistant strains relative to that of GZY803. SM21-resistant haploid strains and GZY803 at mid-log phrase (108 cell/ml) were treated with 2 μg/ml of SM21 for 1 or 4 h. Subsequently, RNA were extracted and sent for sequencing. RNA from untreated cultures was also collected as a control. (A) Cluster display of genes exhibited significant changes in fold expression in the resistant strains as compared to GZY803. Biological processes which are significantly upregulated (B) or downregulated (C) in the resistant strains compared to the wild-type control are shown.
Exploring the Role of Amino Acid Supplements in SM21 Resistance
The first set of foregoing experiments found that SM21 inhibits the activities of mitochondria, which are also known to be involved in amino acid synthesis (Palmieri et al., 1997; Birsoy et al., 2015). The second set of experiments suggested that the resistant strains exhibited SM21-induced transcriptomic inhibition of mitochondrial activity but upregulated genes related to peptide metabolism under the treatment. In addition, the resistant cells also expressed higher level of genes encoding secreted aspartic protease activity as compared to GZY803. Secreted aspartic proteases (SAPs) are a family of at least 10 members (Naglik et al., 2003) which have been shown to degrade host tissues and proteins and peptides for nutrient amino acids supplement and invasion (Schaller et al., 2003). Taking the above findings into account, we hypothesized that there is a link between amino acid supplementation and SM21 resistance in C. albicans. To test the hypothesis, we supplied the wild-type and resistant C. albicans cells with peptone, a form of digested peptides from animal tissues, and examined the susceptibility changes of the strains. Indeed, supplementing the haploid GZY803 and the diploid SC5315 strains with peptone increased the MIC90 of C. albicans against SM21 from 1 to 2 μg/ml (Table 3). Adding 10% of peptone into the cultures of the resistant strains also increased their SM21 MIC90 from 4 to 8 μg/ml. On the contrary, there was no change in MIC90 when peptone supplement was provided to the amphoterin B-treated cultures of C. albicans (Table 3). In conclusion, our data have shown that improving amino acid supplement plays a role in SM21 resistance in C. albicans.
Table 3. Minimum inhibition concentration of Candida albicans strains under 10% peptone supplement.
Discussion
The recent isolation of C. albicans haploids has shown its utility in genetic and biofilm studies (Seneviratne et al., 2015; Truong et al., 2016). Herein, we provide the first demonstration of the use of haploid C. albicans model in elucidating the mechanism of action of a novel antifungal agent SM21. Initially, we obtained the transcriptomic profiles of haploid C. albicans in response to SM21. These analyses provided vital information about the possible mechanism of action of SM21, indicating fungal mitochondria as a potential target. Further downstream experiments using both haploid and diploid C. albicans provided substantial evidence on this claim. Moreover, use of haploid C. albicans helped us to efficiently generate mutations and isolate SM21-resistant strains. Hence, this proof-of-concept study has clearly demonstrated the advantages of the haploid model of C. albicans to the fungal research community. Interestingly, not only in C. albicans, uses of haploid cells have been exploited in many organisms. Due to its advantages of allowing simpler genetic manipulation, haploid cells have been used widely in recessive phenotype analysis in various species including the yeast S. cerevisae, bees, and some type of plants, in which haploid cells exist naturally as part of their life cycle (Haber, 2012; Duina et al., 2014). Even in mammals, where haploid state is only restricted to gametes, haploid embryos have been used for genetic screening (Elling et al., 2011; Yang et al., 2012). Hence, the new haploid model of C. albicans will be a useful tool for researchers working on antifungal drug discovery studies.
Using haploid C. albicans, we found that SM21 inhibits fungal mitochondria and alters mitochondrial activity. SM21 inhibited ATP synthesis in fungal cells via reducing the expression of genes in the ATP synthesis machinery (Figures 1–3). ATP is a necessary form of energy in all types of eukaryotic cells (Kaim and Dimroth, 1999; von Ballmoos et al., 2009). Therefore, inhibiting ATP generation in Candida cells would be lethal to the fungal pathogen. Limiting ATP-coupled electron transportation leads to increased electron leakage, causing accumulation of ROS inside Candida cells (Murphy, 2009; Divakaruni and Brand, 2011). Increased ROS levels causes damages to proteins, DNA, and lipids, thereby compromising C. albicans viability (Cabiscol et al., 2000). In addition, lipid peroxidation which occurs as a result of the reaction of free radicals leads to damage in cell membranes (Dantas Ada et al., 2015). In fact, this is in line with our previous observations of the activity of SM21 on C. albicans (Wong et al., 2014). High levels of ROS is also a signal for Candida cells to undergo apoptosis (Perrone et al., 2008). In addition, SM21 treatment also altered membrane potentials of C. albicans haploid and diploid cells (Figures 2, 3). Therefore, it is likely that SM21 also induces changes in mitochondrial morphology, which in turn is the signal to cell death (Ishihara et al., 2003; Kubli and Gustafsson, 2012; Song et al., 2017). However, further experiments using targeted mitochondrial proteins are warranted to obtain conclusive evidence.
Previously, it was shown that SM21 is fungal-specific and has no activity on bacteria (Wong et al., 2014). Findings in the present study show that SM21 inhibits fungal cytochrome and ATPases, explaining the foregoing observation (Figures 1–3, Table 1). Cytochromes are proteins embedded in the inner mitochondrial membrane of eukaryotes (Ow et al., 2008). Cytochromes transport protons from the mitochondrial matrix into the intermembrane space, generating an electron gradient for ATPases to synthesize ATP (Ow et al., 2008). SM21, by inhibiting cytochrome gene expression, causes reduction in proton transport and ATP synthesis in the fungal cells. Bacteria, on the other hand, do not possess any membrane-bound organelle like mitochondria, instead, rely on a different machinery to generate energy. Besides glycolysis, bacteria utilize various flavins or non-heme iron components for generating gradients across their cytoplasmic membrane for ATP synthesis (Jurtshuk, 1996). This explains why SM21 is not effective for bacteria.
Targeting fungal mitochondria was recently proposed to be a promising therapeutic approach for developing novel antifungal agents. A study by Dagley et al. (2011) demonstrated the involvement of C. albicans Sam37, a mitochondrial outer-membrane protein, in caspofungin susceptibility (Dagley et al., 2011). Similarly, mitochondrial mutants in C. glabrata were found to be associated with increased susceptibility to amphotericin B (Vandeputte et al., 2009). Likewise, inactivation of the mitochondrial protein Goa1, which is required for mitochondrial biogenesis, leads to reduced virulence in systemic candidiasis mouse model (Bambach et al., 2009; Becker et al., 2010). Furthermore, several fungal mitochondrial factors do not have close orthologs in humans (Shingu-Vazquez and Traven, 2011). Therefore, targeting mitochondrial functions could be an effective novel antifungal therapy against Candida, especially for the isolates that display resistance to existing antifungal agents. More evidence on this aspect has come from recent studies of natural compounds which inhibit the activity of C. albicans mitochondria, but not toxic to human cells (Lee et al., 2017; Tian et al., 2017). Similarly SM21 was also shown to be non-toxic to human cell lines (Wong et al., 2014), although it causes mitochondrial dysfunction in fungal cells. In addition to that, SM21 displays fungicidal properties against multidrug resistant Candida strains, such as the T 1549 C. guilliermondii isolates (Wong et al., 2014). Hence, SM21 is proven to be a substantial and effective antifungal molecule for future clinical anti-Candida drug development.
Further in this study, using the C. albicans haploid strains, we have developed SM21 haploid resistant strains and uncovered different pathways altered in the resistant cultures under SM21 treatment relative to GZY803 (Figures 5–7). From the changes observed, we examined and validated the involvement of amino acid supplements in SM21 resistance (Table 3). Mitochondria are known to be the site of many important metabolic and biosynthetic reactions, such as the tricarboxylic acid (TCA) cycle and synthesis of certain essential amino acids (Epstein et al., 2001). Indeed, failure of amino acid homeostasis results in cell death following proteasome inhibition in yeast (Suraweera et al., 2012). Hence, it is possible that amino acid supplements, by compensating for the loss at mitochondrial level, increases the survival rate under SM21 treatment. In addition, all resistant strains have been found to display higher gene expression for SAP activity than GZY803 (Figure 7). Increased SAP activity has been reported in C. albicans when exposed to antifungal agents such as fluconazole, itriconazole, and caspofungin (Copping et al., 2005). SAPs are an important virulence attribute of C. albicans and SAP-deletion strains have been shown to be less virulent and lack tissue invasion (Sanglard et al., 1997; Schaller et al., 2003). SAP is also crucial in providing nitrogen nutrient to Candida cells (Brown et al., 2014; Ramachandra et al., 2014). The information is in line with our hypothesis that higher levels of amino acid supplements benefit the resistance to SM21 in C. albicans. SAP expression is known to be induced by Efg1 and Cph1 and suppressed by Tup1 transcriptional factors (Naglik et al., 2004). The mechanism behind SAP regulation in SM21 resistance in C. albicans is an interesting aspect yet to be uncovered.
In conclusion, this study has shown, for the first time, the use of haploid C. albicans to uncover the mechanism of action of novel antifungal SM21 and C. albicans' resistance mechanism against SM21. Identification of SM21 as an inhibitor of fungal mitochondrial activity provides critical information to the development of this small molecule as a new antifungal agent in future clinical studies. Further experiments using targeted fungal mitochondrial proteins are warranted to obtain conclusive evidence on the proposed mechanism of action. The new insight given by this study on the new haploid C. albicans model will be invaluable for researchers working on antifungal drug discovery.
Availability of Data and Materials
The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE104900 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104900).
Author Contributions
CS, YW, and TT conceived, designed the experiments, and co-wrote the manuscript. TT performed the drug treatment, resistant mutation generation, RNA extraction, bioinformatics analysis, confocal imaging, and other biochemical experiments. GZ generated some fungal strains and performed the FACS ploidy examination. TS helped in the DNA extraction. JL, TL, and LL gave guidance and assisted in the bioinformatics analysis. CT gave general advices and assisted in the experiments. All authors revised the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer AC and handling Editor declared their shared affiliation.
Acknowledgments
We would like to thank Faculty of Dentistry, National University of Singapore, Institute of Molecular and Cell Biology, Agency for Science, Technology, and Research (A*STAR), and School of Computing, University of South Australia, Australia for providing facility and other administrative supports to complete the project. This study was supported by National Medical Research Council (R-221-000-086-511, and R-221-000-073-511), to CS and funding from Agency for Science, Technology, and Research (A*STAR) to YW.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2018.00164/full#supplementary-material
References
Anderson, J. B., Sirjusingh, C., and Ricker, N. (2004). Haploidy, diploidy and evolution of antifungal drug resistance in Saccharomyces cerevisiae. Genetics 168, 1915–1923. doi: 10.1534/genetics.104.033266
Ashman, R. B., Farah, C. S., Wanasaengsakul, S., Hu, Y., Pang, G., and Clancy, R. L. (2004). Innate versus adaptive immunity in Candida albicans infection. Immunol. Cell Biol. 82, 196–204. doi: 10.1046/j.0818-9641.2004.01217.x
Bambach, A., Fernandes, M. P., Ghosh, A., Kruppa, M., Alex, D., Li, D., et al. (2009). Goa1p of Candida albicans localizes to the mitochondria during stress and is required for mitochondrial function and virulence. Eukaryotic Cell 8, 1706–1720. doi: 10.1128/EC.00066-09
Becker, J. M., Kauffman, S. J., Hauser, M., Huang, L., Lin, M., Sillaots, S., et al. (2010). Pathway analysis of Candida albicans survival and virulence determinants in a murine infection model. Proc. Natl. Acad. Sci. U.S.A. 107, 22044–22049. doi: 10.1073/pnas.1009845107
Birsoy, K., Wang, T., Chen, W. W., Freinkman, E., Abu-Remaileh, M., and Sabatini, D. M. (2015). An essential role of the mitochondrial electron transport chain in cell proliferation is to enable aspartate synthesis. Cell 162, 540–551. doi: 10.1016/j.cell.2015.07.016
Bleazard, W., McCaffery, J. M., King, E. J., Bale, S., Mozdy, A., Tieu, Q., et al. (1999). The dynamin-related GTPase Dnm1 regulates mitochondrial fission in yeast. Nat. Cell Biol. 1, 298–304. doi: 10.1038/13014
Brown, A. J., Brown, G. D., Netea, M. G., and Gow, N. A. (2014). Metabolism impacts upon Candida immunogenicity and pathogenicity at multiple levels. Trends Microbiol. 22, 614–622. doi: 10.1016/j.tim.2014.07.001
Cabiscol, E., Piulats, E., Echave, P., Herrero, E., and Ros, J. (2000). Oxidative stress promotes specific protein damage in Saccharomyces cerevisiae. J. Biol. Chem. 275, 27393–27398. doi: 10.1074/jbc.M003140200
Copping, V. M., Barelle, C. J., Hube, B., Gow, N. A., Brown, A. J., and Odds, F. C. (2005). Exposure of Candida albicans to antifungal agents affects expression of SAP2 and SAP9 secreted proteinase genes. J. Antimicrob. Chemother. 55, 645–654. doi: 10.1093/jac/dki088
Dagley, M. J., Gentle, I. E., Beilharz, T. H., Pettolino, F. A., Djordjevic, J. T., Lo, T. L., et al. (2011). Cell wall integrity is linked to mitochondria and phospholipid homeostasis in Candida albicans through the activity of the post-transcriptional regulator Ccr4-Pop2. Mol. Microbiol. 79, 968–989. doi: 10.1111/j.1365-2958.2010.07503.x
Dantas Ada, S., Day, A., Ikeh, M., Kos, I., Achan, B., and Quinn, J. (2015). Oxidative stress responses in the human fungal pathogen, Candida albicans. Biomolecules 5, 142–165. doi: 10.3390/biom5010142
Divakaruni, A. S., and Brand, M. D. (2011). The regulation and physiology of mitochondrial proton leak. Physiology 26, 192–205. doi: 10.1152/physiol.00046.2010
Duina, A. A., Miller, M. E., and Keeney, J. B. (2014). Budding yeast for budding geneticists: a primer on the Saccharomyces cerevisiae model system. Genetics 197, 33–48. doi: 10.1534/genetics.114.163188
Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30, 207–210. doi: 10.1093/nar/30.1.207
Eisen, M. B., Spellman, P. T., Brown, P. O., and Botstein, D. (1998). Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. U.S.A. 95, 14863–14868. doi: 10.1073/pnas.95.25.14863
Elling, U., Taubenschmid, J., Wirnsberger, G., O'malley, R., Demers, S. P., Vanhaelen, Q., et al. (2011). Forward and reverse genetics through derivation of haploid mouse embryonic stem cells. Cell Stem Cell 9, 563–574. doi: 10.1016/j.stem.2011.10.012
Epstein, C. B., Waddle, J. A., Hale, W. T., Dave, V., Thornton, J., MacAtee, T. L., et al. (2001). Genome-wide responses to mitochondrial dysfunction. Mol. Biol. Cell 12, 297–308. doi: 10.1091/mbc.12.2.297
Haber, J. E. (2012). Mating-type genes and MAT switching in Saccharomyces cerevisiae. Genetics 191, 33–64. doi: 10.1534/genetics.111.134577
Hickman, M. A., Zeng, G., Forche, A., Hirakawa, M. P., Abbey, D., Harrison, B. D., et al. (2013). The 'obligate diploid' Candida albicans forms mating-competent haploids. Nature 494, 55–59. doi: 10.1038/nature11865
Ishihara, N., Jofuku, A., Eura, Y., and Mihara, K. (2003). Regulation of mitochondrial morphology by membrane potential, and DRP1-dependent division and FZO1-dependent fusion reaction in mammalian cells. Biochem. Biophys. Res. Commun. 301, 891–898. doi: 10.1016/S0006-291X(03)00050-0
Jones, T., Federspiel, N. A., Chibana, H., Dungan, J., Kalman, S., Magee, B. B., et al. (2004). The diploid genome sequence of Candida albicans. Proc. Natl. Acad. Sci. U.S.A. 101, 7329–7334. doi: 10.1073/pnas.0401648101
Jurtshuk, P. Jr. (1996). “Chapter 4: Bacterial metabolism,” in Medical Microbiology, 4th Edn, ed S. Baron (Galveston, TX: University of Texas Medical Branch at Galveston). Available online at: https://www.ncbi.nlm.nih.gov/books/NBK7919/
Kaim, G., and Dimroth, P. (1999). ATP synthesis by F-type ATP synthase is obligatorily dependent on the transmembrane voltage. EMBO J. 18, 4118–4127. doi: 10.1093/emboj/18.15.4118
Kubli, D. A., and Gustafsson, A. B. (2012). Mitochondria and mitophagy: the yin and yang of cell death control. Circ. Res. 111, 1208–1221. doi: 10.1161/CIRCRESAHA.112.265819
Lee, H., Hwang, J. S., and Lee, D. G. (2017). Scolopendin, an antimicrobial peptide from centipede, attenuates mitochondrial functions and triggers apoptosis in Candida albicans. Biochem. J. 474, 635–645. doi: 10.1042/BCJ20161039
Maere, S., Heymans, K., and Kuiper, M. (2005). BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21, 3448–3449. doi: 10.1093/bioinformatics/bti551
Mavrianos, J., Berkow, E. L., Desai, C., Pandey, A., Batish, M., Rabadi, M. J., et al. (2013). Mitochondrial two-component signaling systems in Candida albicans. Eukaryotic Cell 12, 913–922. doi: 10.1128/EC.00048-13
Murphy, M. P. (2009). How mitochondria produce reactive oxygen species. Biochem. J. 417, 1–13. doi: 10.1042/BJ20081386
Naglik, J., Albrecht, A., Bader, O., and Hube, B. (2004). Candida albicans proteinases and host/pathogen interactions. Cell. Microbiol. 6, 915–926. doi: 10.1111/j.1462-5822.2004.00439.x
Naglik, J. R., Challacombe, S. J., and Hube, B. (2003). Candida albicans secreted aspartyl proteinases in virulence and pathogenesis. Microbiol. Mol. Biol. Rev. 67, 400–428, table of contents. doi: 10.1128/MMBR.67.3.400-428.2003
Noble, S. M., and Johnson, A. D. (2007). Genetics of Candida albicans, a diploid human fungal pathogen. Annu. Rev. Genet. 41, 193–211. doi: 10.1146/annurev.genet.41.042007.170146
Ow, Y. P., Green, D. R., Hao, Z., and Mak, T. W. (2008). Cytochrome c: functions beyond respiration. Nat. Rev. Mol. Cell Biol. 9, 532–542. doi: 10.1038/nrm2434
Palmieri, L., De Marco, V., Iacobazzi, V., Palmieri, F., Runswick, M. J., and Walker, J. E. (1997). Identification of the yeast ARG-11 gene as a mitochondrial ornithine carrier involved in arginine biosynthesis. FEBS Lett. 410, 447–451. doi: 10.1016/S0014-5793(97)00630-3
Perrone, G. G., Tan, S. X., and Dawes, I. W. (2008). Reactive oxygen species and yeast apoptosis. Biochim. Biophys. Acta 1783, 1354–1368. doi: 10.1016/j.bbamcr.2008.01.023
Ramachandra, S., Linde, J., Brock, M., Guthke, R., Hube, B., and Brunke, S. (2014). Regulatory networks controlling nitrogen sensing and uptake in Candida albicans. PLoS ONE 9:e92734. doi: 10.1371/journal.pone.0092734
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e47. doi: 10.1093/nar/gkv007
Sanglard, D., Hube, B., Monod, M., Odds, F. C., and Gow, N. A. (1997). A triple deletion of the secreted aspartyl proteinase genes SAP4, SAP5, and SAP6 of Candida albicans causes attenuated virulence. Infect. Immun. 65, 3539–3546.
Schaller, M., Bein, M., Korting, H. C., Baur, S., Hamm, G., Monod, M., et al. (2003). The secreted aspartyl proteinases Sap1 and Sap2 cause tissue damage in an in vitro model of vaginal candidiasis based on reconstituted human vaginal epithelium. Infect. Immun. 71, 3227–3234. doi: 10.1128/IAI.71.6.3227-3234.2003
Seneviratne, C. J., Zeng, G., Truong, T., Sze, S., Wong, W., Samaranayake, L., et al. (2015). New “haploid biofilm model” unravels IRA2 as a novel regulator of Candida albicans biofilm formation. Sci. Rep. 5:12433. doi: 10.1038/srep12433
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Shields, R. K., Nguyen, M. H., and Clancy, C. J. (2015). Clinical perspectives on echinocandin resistance among Candida species. Curr. Opin. Infect. Dis. 28, 514–522. doi: 10.1097/QCO.0000000000000215
Shingu-Vazquez, M., and Traven, A. (2011). Mitochondria and fungal pathogenesis: drug tolerance, virulence, and potential for antifungal therapy. Eukaryotic Cell 10, 1376–1383. doi: 10.1128/EC.05184-11
Song, S. B., Jang, S. Y., Kang, H. T., Wei, B., Jeoun, U. W., Yoon, G. S., et al. (2017). Modulation of mitochondrial membrane potential and ROS generation by nicotinamide in a manner independent of SIRT1 and mitophagy. Mol. Cells 40, 503–514. doi: 10.14348/molcells.2017.0081
Suraweera, A., Munch, C., Hanssum, A., and Bertolotti, A. (2012). Failure of amino acid homeostasis causes cell death following proteasome inhibition. Mol. Cell 48, 242–253. doi: 10.1016/j.molcel.2012.08.003
Tian, J., Lu, Z., Wang, Y., Zhang, M., Wang, X., Tang, X., et al. (2017). Nerol triggers mitochondrial dysfunction and disruption via elevation of Ca2+ and ROS in Candida albicans. Int. J. Biochem. Cell Biol. 85, 114–122. doi: 10.1016/j.biocel.2017.02.006
Truong, T., Zeng, G., Qingsong, L., Kwang, L. T., Tong, C., Chan, F. Y., et al. (2016). Comparative ploidy proteomics of candida albicans biofilms unraveled the role of the AHP1 gene in the biofilm persistence against amphotericin, B. Mol. Cell. Proteomics 15, 3488–3500. doi: 10.1074/mcp.M116.061523
Vandeputte, P., Tronchin, G., Rocher, F., Renier, G., Berges, T., Chabasse, D., et al. (2009). Hypersusceptibility to azole antifungals in a clinical isolate of Candida glabrata with reduced aerobic growth. Antimicrob. Agents Chemother. 53, 3034–3041. doi: 10.1128/AAC.01384-08
von Ballmoos, C., Wiedenmann, A., and Dimroth, P. (2009). Essentials for ATP synthesis by F1F0 ATP synthases. Annu. Rev. Biochem. 78, 649–672. doi: 10.1146/annurev.biochem.78.081307.104803
Wang, Z., Gerstein, M., and Snyder, M. (2009). RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet. 10, 57–63. doi: 10.1038/nrg2484
Webb, B. C., Thomas, C. J., Willcox, M. D., Harty, D. W., and Knox, K. W. (1998). Candida-associated denture stomatitis. Aetiology and management: a review. Part 2. Oral diseases caused by Candida species. Aust. Dent. J. 43, 160–166. doi: 10.1111/j.1834-7819.1998.tb00157.x
Williams, D., and Lewis, M. (2011). Pathogenesis and treatment of oral candidosis. J. Oral Microbiol. 3, 5571–5582. doi: 10.3402/jom.v3i0.5771
Wisplinghoff, H., Bischoff, T., Tallent, S. M., Seifert, H., Wenzel, R. P., and Edmond, M. B. (2004). Nosocomial bloodstream infections in US hospitals: analysis of 24,179 cases from a prospective nationwide surveillance study. Clin. Infect. Dis. 39, 309–317. doi: 10.1086/421946
Wong, S. S., Kao, R. Y., Yuen, K. Y., Wang, Y., Yang, D., Samaranayake, L. P., et al. (2014). In vitro and in vivo activity of a novel antifungal small molecule against candida infections. PLoS ONE 9:e85836. doi: 10.1371/journal.pone.0085836
Woodford, N., and Ellington, M. J. (2007). The emergence of antibiotic resistance by mutation. Clin. Microbiol. Infect. 13, 5–18. doi: 10.1111/j.1469-0691.2006.01492.x
Yang, H., Shi, L., Wang, B. A., Liang, D., Zhong, C., Liu, W., et al. (2012). Generation of genetically modified mice by oocyte injection of androgenetic haploid embryonic stem cells. Cell 149, 605–617. doi: 10.1016/j.cell.2012.04.002
Keywords: Candida albicans, haploid, SM21, mechanism of action, mitochondria
Citation: Truong T, Suriyanarayanan T, Zeng G, Le TD, Liu L, Li J, Tong C, Wang Y and Seneviratne CJ (2018) Use of Haploid Model of Candida albicans to Uncover Mechanism of Action of a Novel Antifungal Agent. Front. Cell. Infect. Microbiol. 8:164. doi: 10.3389/fcimb.2018.00164
Received: 01 March 2018; Accepted: 01 May 2018;
Published: 08 June 2018.
Edited by:
Priya Uppuluri, University of California, Los Angeles and University of Texas at San Antonio, United StatesReviewed by:
Ashok K. Chaturvedi, University of Texas at San Antonio, United StatesRajendra Prasad, Jawaharlal Nehru University, India
Copyright © 2018 Truong, Suriyanarayanan, Zeng, Le, Liu, Li, Tong, Wang and Seneviratne. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yue Wang, bWNid2FuZ3lAaW1jYi5hLXN0YXIuZWR1LnNn
Chaminda J. Seneviratne, amF5YUBudXMuZWR1LnNn