New Insights in Saccharomyces cerevisiae Response to the Cyanotoxin Microcystin-LR, Revealed by Proteomics and Gene Expression

Microcystins (MCs) are hepatotoxins produced by some cyanobacteria. They are cyclic peptides that inhibit the serine/threonine protein phosphatases (PPs) PP1 and PP2A, especially PP2A. The inhibition of PP2A triggers a series of molecular events, which are responsible for most MC cytotoxic and genotoxic effects on animal cells. It is also known that MCs induce oxidative stress in cells due to the production of reactive oxygen species (ROS). However, a complete characterization of the toxic effects of MCs is still not accomplished. This study aimed to clarify additional molecular mechanisms involved in MC-LR toxicity, using Saccharomyces cerevisiae as eukaryotic model organism. First, a shotgun proteomic analysis of S. cerevisiae VL3 cells response to 1 nM, 10 nM, 100 nM, and 1 μM MC-LR was undertaken and compared to the control (cells not exposed to MC-LR). This analysis revealed a high number of proteins differentially expressed related with gene translation and DNA replication stress; oxidative stress; cell cycle regulation and carbohydrate metabolism. Inference of genotoxic effects of S. cerevisiae VL3 cells exposed to different concentrations of MC-LR were evaluated by analyzing the expression of genes Apn1, Apn2, Rad27, Ntg1, and Ntg2 (from the Base Excision Repair (BER) DNA repair system) using the Real-Time RT-qPCR technique. These genes displayed alterations after exposure to MC-LR, particularly the Apn1/Apn2/Rad27, pointing out effects of MC-LR in the Base Excision Repair system (BER). Overall, this study supports the role of oxidative stress and DNA replication stress as important molecular mechanisms of MC-LR toxicity. Moreover, this study showed that even at low-concentration, MC-LR can induce significant changes in the yeast proteome and in gene expression.

be subsequently repaired by the DNA BER system". We would like to address this hypothesis with this study.
Formerly, when we performed preliminary analysis of the relative expression of yeast's genes PP2A and also of Ntg1 and Ntg2 (from BER DNA repair system), we observed more alterations of these genes at 1 nM MC-LR, opposing to the ones observed with higher MC-LR concentrations (100 nM and 1 uM) [17]. This dependence on the strength of the stimulus, i.e., the dose and/or the exposure time, was also previously observed by Humpage and Falconer (1999) [23]. The authors reported that cytokinesis was stimulated and the rate of apoptosis reduced in the presence of picomolar concentrations of MC-LR. At higher nanomolar concentrations, cytokinesis was inhibited and cell death was induced [23]. In Vero-E6 cells, MC-LR also showed a dual effect in the stimulation of cell proliferation: induction of autophagy/apoptosis or ultimately necrosis highly dependent on the dose and/or the exposure time [16,[24][25][26].
In this study, we aimed to gain more insights regarding the molecular mechanisms involved in MC-LR toxicity, using S. cerevisiae as eukaryotic model organism. Cytotoxicity was evaluated using MTT (methyl-thiazolyl-tetrazolium) assays, one of the most commonly used tests to evaluate cells viability. To identify the proteins altered by exposure to MC-LR, a shotgun proteomic analysis of S. cerevisiae VL3 cells response to 1 nM, 10 nM, 100 nM, and 1 µM MC-LR was undertaken and compared to the control (cells grow in Yeast Extract-Peptone-Dextrose (YPD) medium, without MC-LR). The genotoxic effects of S. cerevisiae VL3 cells were inferred by analyzing the gene expression profiles of genes Apn1, Apn2, Rad27, Ntg1 and Ntg2 (from the BER DNA repair system) using the Quantitative Real-Time RT-PCR (Real-Time RT-qPCR) technique.

Saccharomyces cerevisiae Cells Viability When Exposed to Different Concentrations of MC-LR
The MTT assays were conducted with S. cerevisiae commercial strain VL3. Two positive controls were used (Hydrogen peroxide, H 2 O 2 and Sodium Dodecyl-Sulphate, SDS) which showed significant differences (p ≤ 0.05) from the control (cells grown in YPD). We obtained a 59.4% decrease in relative cell viability in cells treated with 50 mM H 2 O 2 and a decrease of 89.2% when 5% SDS was added to the medium (Table 1). These results confirm that the MTT method is suitable to report cell viability in S. cerevisiae cells. However, when the MTT viability assays were used to asses cytotoxic effects of the different MC-LR concentrations (exposed during 4 h), they did not reveal significant cytotoxic effects in S. cerevisiae (p > 0.05) ( Table 1). Table 1. Relative cell viability (%) from the control of S. cerevisiae VL3 strain exposed to two positive controls (H 2 O 2 or SDS) and to different MC-LR concentrations. The values are the average of three different biological assays (n = 3). Asterisks (*) represent significant differences from the control (p < 0.05; t-Test).

Protein Expression
Proteomics analysis (shotgun approach) led to the identification of 738 proteins from a total of 15 samples analyzed (five experimental groups). The names and expression values of these proteins are reported in Supplementary Table S1. Statistics (ANOVA test) revealed changes in protein expression with exposure to MC-LR toxin. Overall, 69 proteins were significantly up-or down-regulated after exposure to four concentrations of this cyanotoxin. The expression of these proteins is reported in Figure 1, using a color gradient to represent expression levels (heatmap). All information regarding the relative expression of these proteins and statistical analysis can be found in the Supplementary  Table S2. We verified that more proteins were up-regulated (43 proteins) and less down-regulated (26 proteins). Cluster analysis of the various experimental groups (MC-LR concentrations, to which the cells where exposed to) (Figure 1, top cluster analysis) clearly separates the control group from the group exposed to the highest concentration of the toxin (1 µM). This indicates that there are more differences in the protein expression profiles between these two sample groups, than between control and 10 nM group or between 1 nM, 10 nM, and 100 nM groups. This result indicates that 1 nM, 10 nM, and 100 nM groups display less differences in protein expression ( Figure 1).

Protein Functional Analysis
Functional analysis of all differentially expressed proteins between all sample groups was performed using the Search Tool for Retrieval of Interacting Genes (STRING) [27]. This program performs database analysis of protein interactions. From the network of protein interactions predicted by the program, it is possible to infer the molecular processes and pathways in which proteins are involved. Molecular processes can be revealed as subsets of proteins (protein clusters) displaying dense interaction networks in comparison to the background proteome network.
This analysis revealed a complex protein interaction network, suggesting multiple functional relationships between the proteins affected by MC-LR toxin ( Figure 2).
The analysis revealed also two other clusters (proteins in blue and purple) constituted by proteins involved in protein folding and chaperoning activity (HSP10, MGE1, CCT6, MIA40), cell redox homeostasis and oxidative stress (TRR1), proteins involved in DNA repair in response to DNA damage (FUM) and ATP synthesis and membrane transport mediated by ATP (ATP7, ATP1, VPH1, TIM11).  Other proteins showed to be poorly related with these clusters and established few interactions with proteins from the clusters identified. These include PNC1, APE2, MVD1, CPR3, GRX1, that are proteins related with cofactor biosynthesis, proteolysis, isoprenoid biosynthesis, protein folding, cell redox and cellular response to oxidative stress, respectively.
The link between MYO2, SAR1 and SEC21, revealed in the network, reflect the participation of these proteins in common pathways, such as in endoplasmic reticulum (ER) and golgi transport and cell cycle-regulated transport of organelles and proteins.
In Table 2, we report the functions of key protein markers of the processes affected by MC-LR, as well as the toxin concentrations that induced their expression alteration. The table shows that most of these proteins were affected by more than one MC-LR concentration.

Relative Gene Expression Alterations of BER Genes in Saccharomyces cerevisiae Exposed to Different MC-LR Concentrations
The expression results of Ntg1 and Ntg2 genes of the BER system are displayed in Figure 3. In the case of Ntg1 the presence of MC-LR did not revealed significant differences from the control, although MC-LR at 1 nM and 10 nM led to a slight, but not significant, over-expression of this gene. The Ntg2 gene displayed similar alterations but showed significant statistical differences from the control in the presence of 1 and 10 nM MC-LR (p < 0.05).  Results are expressed as the mean ± standard deviation of five independent biological experiments. Asterisks (*) represent statistically significant differences between the cells exposed to MC-LR and the control cells (p < 0.05; t-Test).
Apn1, Apn2, and Rad27 genes of the BER system were significantly altered in the presence of MC-LR, however displayed different expression profiles (Figure 4). There was an up-regulation of Apn1 and Rad27 genes in response to the lower doses of MC-LR tested (1 nM and 10 nM). On the contrary, for the higher doses (100 nM and 1 µM), there is a tendency for the repression of these genes.
The Apn2 gene showed an induction in the presence of 1 nM, 10 nM, and 100 nM of MC-LR, only being repressed with the highest dose tested (1 µM). Results are expressed as the mean ± standard deviation of five independent biological experiments. Asterisks (*) represent statistically significant differences between the cells exposed to MC-LR and the control cells (p < 0.05; t-Test).

Discussion
The MTT assay is a method to determine cell viability through the reduction of 3-(4,5-dimethyl-2-thizaolyl)-2,5-diphenyl-2H-tetrazolium bromide into a colored product, formazan, by mitochondria dehydrogenases to which the cell membrane is impermeable [28]. The MTT viability assays were used in this study to assess cytotoxic effects MC-LR on yeast cells, revealing no significant cytotoxic effects in S. cerevisiae VL3 (p > 0.05).
To verify if the absence of cytotoxic effects was related to some insensitivity problem of the VL3 strain, we tested a second strain: S. cerevisiae VR5 (Fermicru). The results obtained with this second strain were similar (data not shown), thus confirming that the MTT assay does not reveal any cytotoxic effect when the S. cerevisiae cells are exposed to MC-LR. There are not many studies that use the MTT method with S. cerevisiae, nevertheless, the MTT assays are able to display decreases of S. cerevisiae viability as verified with the positive controls used (H 2 O 2 and SDS).
Previous studies, using MTT assays, revealed that there is a negative correlation between cell viability and MC-LR concentration, where there is a significant dose/time dependent cytotoxic effect of MC-LR on cell viability [28,29]. However, these studies were applied to mammal cell lines and used considerably higher concentrations of toxins than the ones tested in this study [28,29].
The proteomic analyses performed revealed a complex pattern of responses in yeast cells exposed to MC-LR, involving several molecular processes. With regard to our previous proteomics work [17], we denote the increased analytical capacity of the shotgun methodology over Two-dimensional gel Electrophoresis (2DE) based proteomics to characterize changes in yeast proteome. In our previous work, we were able to detect expression differences of only 14 proteins and functions assigned to only six of them. With the current approach, we were able to characterize in more detail the changes in the proteome of yeast cells exposed to MC-LR.
From the results obtained, it may be inferred that the molecular process most affected by the presence of MC-LR was gene translation. Indeed, the expression of many ribosomal proteins and proteins related with RNA processing, ribosome biogenesis and regulation of translation were affected ( Figure 2). While some of these proteins were down-regulated (RPS6B, RPS20, RPL32, RPL24A, HMT1, NSR1, NOP2), some others were up-regulated (RPS2, RPS23A, RPL36B), making difficult to assess the overall effects of the toxin, and to infer if gene translation was inhibited or increased during yeast exposure to the toxin. However, the decrease in the expression of the translation initiation factors RLI1 and RPG1 (Table 2 and Supplementary Table S2), may anticipate a negative effect of this toxin in gene translation. RLI1 in particular, has shown to be a primary cellular target of ROS, and indeed impairment of the functions of this protein were previously related with increase of ROS toxicity [30]. Both proteins are essential to gene translation. RLI1 is involved in ribosome biogenesis and maturation, translation initiation and ribosome recycling [30] and RPG1 has a key role in the recruitment of pre-initiation complex (PIC) to mRNA [31]. Moreover, alterations in histones and RAD23 suggest the occurrence of alterations in chromatin and DNA caused by MC-LR. RAD23 was up-regulated at high MC-LR concentrations (100 nM and 1 µM) ( Table 2). This protein is involved in proteasome-mediated degradation of misfolded proteins in the endoplasmic reticulum, but also plays a primary function in DNA repair, by supposedly mediating RAD4 stabilization and recognition of DNA damaged sites, in NER mechanism [32]. The repression of this protein has been related with increased sensitivity of cells to DNA-damaging agents [32]. Previously, we detected the increase of the intracellular ROS levels when S. cerevisiae was exposed, in similar conditions, to MC-LR [15]. Increase of ROS was verified even for the lowest MC-LR concentration (1 nM) [15]. Our previous results thus support our current hypothesis that MC-LR caused an increase in intracellular ROS. This ROS increase would explain, for instance, the down-regulation of RLI1 protein. On the same note, the increase in ROS would also explain the increase in the expression of thioredoxin reductase 1 (TRR1). This protein was up-regulated in yeast cells exposed to 100 nM and 1 µM MC-LR (Table 2). Regarding protein functions, TRR1 play a key role in the removal of superoxide radicals and thereby in maintenance of cell redox homeostasis [33].
Furthermore, ROS generation is known to be intrinsically related with mitochondria metabolism [34]. Two mitochondrial proteins, HSP10 and MGE1, were down-regulated with the exposure to the toxin (Table 2), which may indicate that mitochondria functions were affected. HSP10 is essential in protein biogenesis in mitochondria. MGE1 has also chaperonin activity and plays an important role in protein import into mitochondrial matrix. Our results are still insufficient to ascertain the involvement of mitochondria in the response of yeast to MC-LR. Notwithstanding, the role of this organelle in yeast toxicity should not be disregarded, and indeed several effects observed previously, such as the increase of ROS and the nuclear alterations like chromatin condensation and fragmentation, suggestive of apoptotic mechanisms [15], can be explained in light of the mitochondrial functions.
The effects of MC-LR in PP1 and PP2A, should be also examined and discussed, since these proteins exert a profound influence in the proteome and in the whole metabolism of cells and are specific targets of this toxin [35]. Our proteomic results do not evidence any direct alterations in PP1 and PP2A. However, our results revealed alterations in several proteins that are known to interact with PP1 and PP2A, thus giving indirect evidences that PPs can be involved in the alterations observed in yeast proteome. For instance, CCT4 was up-regulated by MC-LR (Table 2). This protein is a subunit of T-complex protein 1 (TCP1). TCP1 play chaperonin functions, i.e., participates in the folding of proteins in a process dependent of ATP. It was also demonstrated that this complex interacts with PP2A [35], being this phosphatase an important regulator of TCP1 functions. According to the recent study describing the human protein phosphatase interactions and dynamics, several other proteins identified in this work can be considered putative PP1 and PP2A interacting targets and thus be regulated by these phosphatases. Among these are the ribosomal proteins (RPL and RPS), RLI1, PGM1 and also the heat shock proteins (HSPs).
Finally, exposure to MC-LR seems also to affect carbohydrate metabolism. Several glycolytic and TCA cycle proteins were altered (PCK1, TDH2, PFK2, IDP1, GPM1, HXK2, ENO1, and SHB17). In particular, the up-regulation of TDH2, IDP1, ENO1 and SHB17 (Figure 1, Supplementary Table S2), suggest that these processes were activated in yeast cells. The increase of TCA and glycolysis may be required to increase ATP pool in yeast cells, and energy to cope with the toxic effects of MC-LR.
In light of the results observed from the MTT test and the proteome, that in general do not evidence acute adverse effects, we hypothesize that the MC-LR concentrations tested are sub-toxic or sub-lethal to yeast cells, as cells did not shown any decrease in viability in the range of concentrations tested, as verified when using the MTT test.
Regarding gene expression obtained by Real-Time RT-qPCR, we denote that Ntg1 and Ntg2 expression results differ in some extent from the results observed in our previous study [17]. No variations were observed in Ntg1 contrasting with the previous over-expression of this gene at 1 nM and the decrease in expression at 100 nM and 1 µM. However, a consistent increase in Ntg2 expression was verified in both studies for MC-LR concentrations between 1 nM and 100 nM. The differences in the expression of Ntg1 between studies are difficult to explain, and apparently are not related with the experimental conditions since these were similar in both studies.
The gene expression levels analyzed by Real-Time RT-qPCR showed that the lower MC-LR concentrations (1 nM and 10 nM) are the ones that resulted in an up-regulation of all the BER system genes. This is in accordance with the explanation that induction of genes expression may contribute to overcome the effects of the toxin when its concentration is low. However, at high toxin concentrations some genes (Apn1 and Rad27) were repressed, which may suggest that the DNA had extended deleterious effects that could not be fully overcome by this repair mechanism.
A previous study compares yeast life span with BER activity [36]. It shows that losses of enzymes Ntg1p and Ntg2p increase both spontaneous and hydrogen peroxide induced mutation frequencies, though it generally does not cause cell death. A single deletion of these enzymes showed little or no effect in life span, however, combined deletions resulted in a decreased cell survival [36]. Moreover, in this study, we did not observe cell death (using MTT test) despite the alterations in Ntg1 and Ntg2 genes.
Apn1p is considered the most important apurinic/apyrimidinic (AP) endonuclease. Without this enzyme the spontaneous mutation rate would increase, becoming more sensitive to hydrogen peroxide. The Apn2p seems to confer resistance to ROS induced damages in an alternative pathway to the Apn1p. The combined loss of both endonucleases exhibited extremely poor cell survival [36]. Our results suggest that MC-LR is able to affect the BER DNA repair system mechanism, particularly the Apn1/Apn2/Rad27 genes. These results are consistent with the hypothesis that the induction of ROS is an important molecular mechanism of MC-LR genotoxicity, and also gives more insights that H 2 O 2 is a relevant ROS induced.
It should also be noted that MC-LR concentration of 1 nM, proposed by World Health Organization (WHO) as guideline value for drinking water, induced a significant change in a high number of proteins and also in the gene expression levels of the BER system (real-time results). This agrees with the finding of He and co-workers [20] where most significant alterations in the proteome were observed in low-concentrations of MC-LR. Unfortunately, we were not able to corroborate the gene expression results at the protein level. We verify that the corresponding gene products were not present in the sub-set of the yeast proteome identified by the shotgun method (Supplementary Table S1). The reason behind the absence of these proteins from our proteomic dataset can be due to their expression in the cell that may be below the detection limits of the proteomics methodology or were underrepresented (e.g., low extraction efficiency of these proteins). RAD23 was the only protein identified whose functions can be linked to DNA repair and to the NER system. The molecular functions of this protein, however, is different from RAD27, which may explain the differences in expression verified at the gene (RAD27) and protein (RAD23) levels. However, taken together, the results attained in Apn1, Apn2, Rad27, Ntg2, and RAD23 strongly suggest that the toxin affects DNA integrity and leads to responses involving DNA repair systems.

Conclusions
In this work, we present new data concerning Saccharomyces cerevisiae VL3 response to the cyanotoxin MC-LR, revealed by proteomics and gene expression. This simple eukaryotic organism has been widely used as a eukaryotic cell model and has been particularly helpful to explore the molecular mechanisms and affected pathways of several compounds in toxicology.
The proteomics study revealed significant changes in the yeast proteome induced by this cyanotoxin. The main proteins alterations here observed were associated to gene translation, DNA replication stress; oxidative stress; cell cycle regulation and carbohydrate metabolism. Alterations in specific proteins, such as RLI1, TRR1, and RAD23, for instance, suggest an increase of ROS in yeast cells and effects them at the DNA level.
The assessment of genotoxic effects exploited using Real-Time RT-qPCR assays showed that MC-LR affects the genes of BER DNA repair system mechanism, supporting the existence of genotoxic effects induced by MC-LR, caused by intracellular ROS levels increase.
The results of this study further support that current guidelines for MCs in drinking water may need a revision since chronic exposure to low doses of MC-LR and associated risks to human health are currently still not fully elucidated.

Saccharomyces cerevisiae Cultures
S. cerevisiae VL3 (Zymaflore, Bordeaux, France) pre-culture was left growing overnight in YPD medium (1% yeast extract, 2% peptone, 2% glucose) at 20 • C, 300 rpm. Then, 20 mL culture was inoculated in five 150 mL flasks with a cell density adjusted to an initial optical density DO 600nm = 0.05, corresponding to approximately 5.5 × 10 5 cells/mL [37,38]. Each flask was treated with different pure MC-LR concentrations: without toxin (control), 1 nM, 10 nM, 100 nM, and 1 µM. The cultures were left growing for 4 hours at a climatic chamber at 20 ± 2 • C and 300 rpm (shaker IKA Labortechnik, IKA, Staufen, Germany). This exposure time was selected by our results previously obtained [15], where we observed that after 4h the cells started to recover, reaching approximately the same OD as the control.
Cells were harvested in 15 mL plastic tubes (Sarstedt) by centrifugation at 834× g on a high capacity centrifuge (Sorvall RT 6000D, Dupont, Wilmington, DE, USA) for 3 min and washed with 1 mL of RNase free sterile water (Gibco, Carlsbad, CA, USA). Cells were transferred into 1.5 mL tubes and centrifuged (Eppendorf 5415C Centrifuge, Eppendorf, Hamburg, Germany) for 3 min at 8000× g. The supernatant was discarded and the pellet was immediately frozen in liquid nitrogen and kept at −80 • C until use.

Analysis of Saccharomyces cerevisiae Viability When Exposed to Different Concentrations of MC-LR
To evaluate cells viability when exposed to different chemical compounds and MC-LR, MTT (methyl-thiazolyl-tetrazolium) assays were performed.
Two different compounds were used as positive controls: 5% sodium dodecyl sulfate (SDS) (Invitrogen, Carlsbad, CA, USA) and 50 mM hydrogen peroxide (H 2 O 2 ) (Sigma-Aldrich, St. Louis, MO, USA). SDS is a detergent commonly used to lyse cells and is thus expected to reduce cell viability. To perform the MTT assay S. cerevisiae VL3 cultures were grown in YPD with different pure MC-LR concentrations: without toxin (control), 1 nM, 10 nM, 100 nM and 1 µM, during 4 h in 6 wells microplates (Nest Biotechnology, Jiangsu, China) (for easier manipulation).
After the exposure, 1 mL of the suspension from each well was harvested to 2 mL tubes and centrifuged for 5 min at 8000× g (Eppendorf 5415C Centrifuge). The supernatant was discarded and the cells were resuspended in 100 µL of PBS (phosphate-buffered saline) buffer. Then, 100 µL of freshly prepared MTT solution was added, and the tubes incubated in the dark in the climatic chamber with agitation for 2 h. Afterwards, the suspension was centrifuged for 5 min at 8000× g and the supernatant was discarded. The pellet was resuspended in 300 µL dimethyl sulfoxide (DMSO) (Thermo Fisher, Carlsbad, CA, USA) and 100 µL of each was inoculated in a 96 well plate (Sarstedt, Nümbrecht, Germany), performing 3 replicates in total for each treatment. The absorbance was read by a microplate spectrophotometer (Thermo Fisher LabSystems) at 570 nm with a reference wavelength of 690 nm.
The relative cell viability (%) was calculated using the following formula: relative cell viability (%) = average absorbance of experimental groups average absorbance of the control group × 100 (1)

LC-MS/MS and Protein Identification
The LC-MS/MS was carried out as described previously in a nano-LC coupled to a hybrid Ion trap mass spectrometer (Linear Trap Quadropole -LTQ Orbitrap Velos Pro-ETD, Thermo Scientific, Waltham, MA, USA) [40]. Peptides were separated by reverse phase chromatography on a 20 mm × 100 µm C18 precolumn followed by a 100 mm × 75 µm C18 column (NanoSeparations, Nieuwkoop, The Netherlands) in a linear gradient of acetonitrile (2% to 95% v/v) in FA (0.1% v/v), at a flow rate of 0.3 µL/min. Full scans were performed at 30,000 resolution at a range of 380-2000 m/z. The top 20 most intense ions were isolated and fragmented with collision induced fragmentation (CID) applying normalized collision energy of 30%, isolation width of 2.0, and activation time of 10 ms and a Q-value of 0.25. In total, 15 independent LC-MS/MS runs were performed from 15 biological samples (including replicates).
Proteins were identified by searching LTQ raw data against Saccharomyces sequences in UNIPROT database using SEQUEST algorithm (Proteome Discoverer software, version 1.4, Thermo Scientific, Waltham, MA, USA) and the X! Tandem algorithm in Scaffold (version Scaffold 4.3.4, Proteome Software, Portland, OR, USA). Peptides were accepted at a probability greater than 95.0% and proteins greater than 99.9%. MS and MS/MS tolerances were set to 10 ppm and 0.6 Da. Static and dynamic modifications were carbamidomethylation and oxidation, respectively. Trypsin was selected for protein cleavage allowing for one missed cleavage.
Functional analysis of proteins was performed with the bioinformatics tool STRING [27]. This analysis was performed by selecting the organism Saccharomyces cerevisiae; the sources of evidence text mining, experiments, databases, co-expression, neighborhood, gene fusion, co-occurrence; the interaction score medium confidence (0.400) and max number of interactors (1st shell) no more than 5.

Saccharomyces cerevisiae Studied Genes and Primers Design
In this study, Real-Time RT-qPCR was used to provide quantitative measurements of gene expression. It was intended to determine how the expression of a particular gene changes in response to MC-LR exposure. Knowing that MC-LR causes oxidative DNA damages that are probably mainly repaired by the S. cerevisiae DNA BER system, the following genes were selected: Apn1, Apn2, Rad27, Ntg1, and Ntg2 to evaluate their expression levels in cells exposed to different concentrations of MC-LR. The primers for Ntg1, Ntg2 genes had been previously designed by [17], and primers for genes Apn1, Apn2, and Rad27 were designed in this study (Table 3). We selected the reference genes Alg9 to be used as internal control [42] for posterior expression normalization.

Saccharomyces cerevisiae RNA Extraction and Purification
To evaluate the alterations of the expression levels of the S. cerevisiae previously selected genes by Real-Time RT-qPCR, S. cerevisiae RNA was extracted using TRIzol method [17]. Afterwards RNA samples were purified using the High Pure RNA Isolation Kit (Roche), according to the manufacturer instructions. To confirm that purified RNA was free from DNA, the samples were subjected to a conventional PCR prior to RT-qPCR analyses.

Evaluation of RT-qPCR Parameters
Before starting the Real-Time RT-qPCR assays to determine the expression levels, the optimal conditions for RT-qPCR performance had to be verified. The linear regression coefficient (R 2 ) must be above 0.95% and efficiency must be between 80-115% [43]. The optimization method was based on the construction of standard calibration curves for Real-Time RT-qPCR using several dilutions of the DNA.
The RT-qPCR conditions were optimized for all the target genes: Apn1, Apn2, Rad27, Ntg1, Ntg2, and the reference gene Alg9, by confirming that reactions had a linear regression coefficient (R 2 ) above 0.95 and the efficiency value was between 88% and 100%, as summarized in Table 4.

Statistics
Proteomic data was analyzed by ANOVA and Student's t-test, using the program Multiple Experiment Viewer [44] for processing OMICs data. Significant differences were assumed for p < 0.05.
A normalization using the mathematical method developed by Pfaffl [45] was conducted to determine the relative expression levels of BER genes (Apn1, Apn2, Rad27, Ntg1, and Ntg2). Five independent assays were made to obtain more consistent results. The expression results where normalized using the ALG9 gene. Funding: This research was funded by the European Regional Development Fund (ERDF) through the COMPETE e Operational Competitiveness Program and national funds through FCT (Foundation for Science and Technology) grant to E. Valerio (SFRH/BPD/75922/2011). It was also funded by the Strategic Funding UIDB/04423/2020 and UIDP/04423/2020 through national funds provided by FCT-Foundation for Science and Technology and European Regional Development Fund (ERDF), in the framework of the program PT2020. This work has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 823860 (TOXICROP).

Conflicts of Interest:
There are no conflict of interest to declare.