# Genome-Wide Mutant Screening in Yeast Reveals that the Cell Wall is a First Shield to Discriminate Light From Heavy Lanthanides

Nicolas Grosjean, Marie Le Jean, Michel Chalot, Héctor M Mora-Montes,

Jean Armengaud, Elisabeth M Gross, Damien Blaudez

► To cite this version:

Nicolas Grosjean, Marie Le Jean, Michel Chalot, Héctor M Mora-Montes, Jean Armengaud, et al. Genome-Wide Mutant Screening in Yeast Reveals that the Cell Wall is a First Shield to Discriminate Light From Heavy Lanthanides. Frontiers in Microbiology, 2022, 13, pp.881535. 10.3389/fmicb.2022.881535. [hal-03840664]

HAL Id: hal-03840664

https://hal.univ-lorraine.fr/hal-03840664

Submitted on 5 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Genome-Wide Mutant Screening in Yeast Reveals that the Cell Wall is a First Shield to Discriminate Light From Heavy Lanthanides

Nicolas Grosjean<fnref n="1" />,²†, Marie Le Jean², Michel Chalot³,⁴, Héctor M. Mora-Montes⁵, Jean Armengaud⁶, Elisabeth M. Gross² and Damien Blaudez<fnref n="1" />*

## OPEN ACCESS

**Edited by:**
Nuno Pereira Mira,
University of Lisbon, Portugal

**Reviewed by:**
Ileana Farcasanu,
University of Bucharest, Romania
Barbara Montanini,
University of Parma, Italy

*Correspondence:
Damien Blaudez
damien.blaudez@univ-lorraine.fr

¹Present address:
Nicolas Grosjean,
Biology Department,
Brookhaven National Laboratory,
Upton, NY, United States

**Specialty section:**
This article was submitted to
Microbiotechnology,
a section of the journal
Frontiers in Microbiology

Received: 22 February 2022

Accepted: 19 April 2022

Published: 19 May 2022

**Citation:**
Grosjean N, Le Jean M, Chalot M, Mora-Montes HM, Armengaud J, Gross EM and Blaudez D (2022)
Genome-Wide Mutant Screening in Yeast Reveals that the Cell Wall is a First Shield to Discriminate Light From Heavy Lanthanides.
Front. Microbiol. 13:881535.
doi: 10.3389/fmicb.2022.881535

The rapidly expanding utilization of lanthanides (Ln) for the development of new technologies, green energies, and agriculture has raised concerns regarding their impacts on the environment and human health. The absence of characterization of the underlying cellular and molecular mechanisms regarding their toxicity is a caveat in the apprehension of their environmental impacts. We performed genomic phenotyping and molecular physiology analyses of *Saccharomyces cerevisiae* mutants exposed to La and Yb to uncover genes and pathways affecting Ln resistance and toxicity. Ln responses strongly differed from well-known transition metal and from common responses mediated by oxidative compounds. Shared response pathways to La and Yb exposure were associated to lipid metabolism, ion homeostasis, vesicular trafficking, and endocytosis, which represents a putative way of entry for Ln. Cell wall organization and related signaling pathways allowed for the discrimination of light and heavy Ln. Mutants in cell wall integrity-related proteins (e.g., Kre1p, Kre6p) or in the activation of secretory pathway and cell wall proteins (e.g., Kex2p, Kex1p) were resistant to Yb but sensitive to La. Exposure of WT yeast to the serine protease inhibitor tosyl phenylalanyl chloromethyl ketone mimicked the phenotype of *kex2Δ* under Ln, strengthening these results. Our data also suggest that the relative proportions of chitin and phosphomannan could modulate the proportion of functional groups (phosphates and carboxylates) to which La and Yb could differentially bind. Moreover, we showed that *kex2Δ*, *kex1Δ*, *kre1Δ*, and *kre6Δ* strains were all sensitive to light Ln (La to Eu), while being increasingly resistant to heavier Ln. Finally, shotgun proteomic analyses identified modulated proteins in *kex2Δ* exposed to Ln, among which several plasmalemma ion transporters that were less abundant and that could play a role in Yb uptake. By combining these different approaches, we unraveled that cell wall components not only act in Ln adsorption but are also active signal effectors allowing cells to differentiate light and heavy Ln. This work paves the way for future investigations to the better understanding of Ln toxicity in higher eukaryotes.

**Keywords:** lanthanum, ytterbium, cell wall, endocytosis, signaling, deletome

<footnote n="1">¹Université de Lorraine, CNRS, LIEC, Nancy, France, ²Université de Lorraine, CNRS, LIEC, Metz, France, ³Laboratoire Chrono-Environnement, Université de Bourgogne Franche-Comté, CNRS, Besançon, France, ⁴Université de Lorraine, Nancy, France, ⁵Departamento de Biología, División de Ciencias Naturales y Exactas, Universidad de Guanajuato, Guanajuato, Mexico, ⁶Département Médicaments et Technologies pour la Santé (DMTS), Université Paris-Saclay, CEA, INRAE, Bagnols-sur-Cèze, France</footnote>
# INTRODUCTION

Lanthanides (Ln) are 15 chemical elements belonging to the rare earth elements (REEs). Two subgroups are usually distinguished based on both the atomic mass and the ionic radius, namely, light REEs (LREEs, La to Eu) and heavy REEs (HREEs, Gd to Lu). These so-called technology-critical elements are key elements in energetic transition, high technology (Hu et al., 2018), medicine, and technological military devices and are used as fertilizers and food supplements in agriculture (He et al., 2001; Hu et al., 2004; Ciacci et al., 2015). High concentrations of lanthanides can be found in soil of mining areas but also have been reported to accumulate in kidney, bone, and brain tissues after exposition to MRI contrasting agent (Ramalho et al., 2016; Gulani et al., 2017). As emerging contaminants, their large production and utilization raise the question of their impact on the environment and human health (Noack et al., 2014; Lee and Wen, 2018).

Although Ln have been considered harmless and non-essential elements (Leeuw and Academic, 2000), results from the last few years reveal the effects of Ln on living organisms. They are cofactors for pyrroloquinoline quinone-dependent methanol/ ethanol dehydrogenases for certain bacteria (Nakagawa et al., 2012; Pol et al., 2014; Skovran and Martinez-Gomez, 2015; Good et al., 2016; Wehrmann et al., 2017; Cotruvo, 2019) but also present antibacterial, antifungal, and nematocidal activities (Wakabayashi et al., 2016) by increasing membrane permeability (Weiwei et al., 2007; Yufeng et al., 2007), generate oxidative stress (Pagano et al., 2015), induce mitotic and chromosomal aberrations (Pagano et al., 2015), and impair human health (Wei et al., 2013; Noack et al., 2014; Gulani et al., 2017). However, the few existing studies mentioning toxicity remain elusive, which reflects the lack of data on cellular and molecular mechanisms causing Ln toxicity.

A few studies have hypothesized that LREEs and HREEs may induce different responses, based on their relative reactivity and chemical differences (Gonzalez et al., 2014). In addition to their different toxicity levels (HREEs > LREEs (Tyler, 2004; Gonzalez et al., 2014; Goecke et al., 2017)), only a few studies revealed the specific modes of action of these two REE groups with a single multi-scale analysis directly comparing the differential effect of two Ln (Grosjean et al., 2021). As such, specific accumulation of LREEs has been reported in ferns (Grosjean et al., 2020), while angiosperms preferentially translocate HREEs in their aboveground tissues (Liu et al., 2018; Grosjean et al., 2019), thereby suggesting different transport systems between LREEs and HREEs. Molecular evidence also supports this hypothesis in *Saccharomyces cerevisiae*, in which disruption of the Ca channel Cch1p/Mid1p restricts La but not Gd uptake (Ene et al., 2015). Therefore, it is essential to unravel the distinct molecular and cellular effects of LREEs and HREEs, since diverse exposure scenarios may occur in the environment, with potentially vastly different consequences for the biota.

The screening of a mutant collection of *S. cerevisiae* has previously proven to be a powerful approach to pinpoint the role of non-essential proteins in toxicity modulation when exposed to different environmental stressors, such as metals

(Scherens and Goffeau, 2004; Jin et al., 2008). Based on the phylogenetic relationships between *S. cerevisiae* and other eukaryotes and the strong conservation of basic stress responses, the mechanisms revealed through this approach shed light on pathways that may also be involved in other eukaryotes (Botstein et al., 1997; Kachroo et al., 2017). Recently, a pooled population of yeast mutants was screened with low Ln concentrations revealing a limited set of mutants and pathways (Pallares et al., 2021). However, while the use of pooled population for mutant screening allows the high-throughput generation of data, it also introduces several biases. For instance, mutant cells can affect one another confounding the results of a pooled screen, generating false positives (So et al., 2019), or cell interactions and competition masking minor loss-of-function phenotype (Kim et al., 2018). Therefore, complementary screening methods are needed to grasp the full response profile of *S. cerevisiae* to an acute Ln exposure. For this purpose, we leveraged an arrayed genome-wide phenotyping screen of the entire set of gene deletion mutants to identify gene products that modulate Ln cellular toxicity at high concentrations of Ln. By employing a method vastly used for other toxicants (heavy metals, chemicals, and physical stressors), it provides an opportunity to compare Ln cellular impacts to other stressors. In this regard, we hypothesized that Ln would present distinct effects from known toxicants and that LREEs and HREEs would induce specific patterns. Therefore, given their relative environmental abundances and utilizations among lanthanides, we selected La and Yb as their respective representatives. The genome-wide mutant screen allowed us to pinpoint the cellular compartments and functions affected by Ln, such as the cell wall organization that we further investigated because of a distinct response between La and Yb. Numerous genes/functions impacted by Ln are conserved in humans and provide new promising hypotheses to study Ln-mediated toxicity in humans.

# MATERIALS AND METHODS

## Yeast Strains and Chemicals

Deletion mutants as well as the wild-type *Saccharomyces cerevisiae* strains BY4741 (MATa; *his3Δ1*; *leu2Δ0*; *met15Δ0*; *ura3Δ0*) and BY4742 (MATα; *his3Δ1*; *leu2Δ0*; *lys2Δ0*; *ura3Δ0*), both isogenic to the S288C strain, were used in this study. They were purchased from EUROSCARF. The haploid deletion library consists of 4,733 mutants for non-essential genes. All Ln were purchased from Sigma-Aldrich (MO, USA) as hydrated chloride salts (LnCl₃, xH₂O). The Kex2 gene was synthesized and cloned into the pYES2 plasmid at the NotI restriction site to obtain the expression vector pYES2-KEX2, under the activity of the GAL1 promoter.

## *Saccharomyces cerevisiae* Genomic Phenotypic Screen With Lanthanides

The 4,733 mutants (BY4741 background) were used for genomic phenotyping. Two Ln (La and Yb) were selected and used as chloride forms (LaCl₃ and YbCl₃, respectively). Ln concentrations for the screen were defined by pilot experiments as 4.5 mM La
and 3.8 mM Yb, allowing discrimination in a single step of both sensitive or resistant mutants (Supplementary Figure S1). Briefly, individual deleted mutants were grown in 96-well master plates in 200 µl of YPD (10 g yeast extract, 20 g peptone, 20 g dextrose) at 28°C until the stationary phase. A Thermo Scientific™ Nunc™ Replication System (250520) was used to replica inoculate the surface of YPD agar plates supplemented with either La, Yb, or Ln-free. Plates were digitally recorded after 5 days of growth at 28°C. Ln-sensitive and Ln-resistant mutants were identified when colony size under Ln exposure was decreased or increased, respectively, compared to WT and neighboring mutants but also to the size of colonies on control plates lacking Ln. Four replications were carried out for each condition. YPD medium is commonly used in metal toxicity screen experiments in yeast. However, its high concentration in phosphates promotes the precipitation of Ln/phosphate complexes that are no longer toxic to the cells and thus justify the use of high concentrations of Ln to achieve Ln toxicity. To verify whether the relatively high concentrations of La and Yb used in this experiment had an influence on the quality of the results, a preliminary screen was performed using 350 mutants (randomly selected) on low concentrations of La (300 µM) and Yb (120 µM) in YNB agar medium devoid of inorganic phosphates and supplemented with 1 mM β-glycerophosphate. No Ln precipitation was visible in this medium. Identical results (26 sensitive or resistant mutants) were obtained for all 350 mutants between the YPD and YNB screening experiments, refuting a putative bias with high Ln concentrations in the YPD medium. Additionally, some mutants grew poorly on YNB medium making screening results more difficult to assess. For these reasons, YPD agar medium was used for the entire primary screen, as well as for the validation screen.

## Genomic Phenotypic Validation Screen

Mutants identified as Ln-sensitive or Ln-resistant in the primary screen were individually confirmed by serial dilution spot assays. Mutants were grown as previously described, and serial tenfold increment dilutions were performed. Five microliters of six dilutions were spotted on YPD agar plates containing 4.5 mM La or 3.8 mM Yb to verify Ln-sensitive mutants and 4.72 mM La and 4 mM Yb for Ln-resistant mutants. The results were observed after 5 days of growth at 28°C. Sensitivity and resistance levels were assigned to the mutants according to the number of dilutions where cells grew. Consequently, mutants exhibiting a reduction in colony-forming ability at the first, second-third, or fourth-fifth dilutions were classified as "high" (HS), "medium" (MS), or "low" (LS) sensitive, respectively. Conversely, mutant strains exhibiting an increase in colony-forming ability at the first-second, third-fourth, or fifth-sixth dilutions were classified as "low" (LR), "medium" (MR), or "high" (HR) resistance, respectively. Twenty-five mutants randomly picked in the identified mutants were also assayed in the haploid MATα strain BY4742 to validate the robustness of the screen.

## Serine Protease Inhibition Assay

To investigate the effect of the serine protease inhibitor tosyl phenylalanyl chloromethyl ketone (TPCK) on Ln resistance,

drop tests of 3-fold serial dilutions of the wild-type strain were exposed to La (3.8 mM) and Yb (3.2 mM) with the addition of 80 µM TPCK in solid YPD and grown for 5 days at 28°C.

## Chitin Staining

Wild-type and mutant cells were grown in YNB medium as described above. Cultures were inoculated at an OD₆₀₀nm of 0.05 and grown overnight, with or without the addition of 160 µM La or 8 µM Yb, concentrations corresponding to the EC₅₀ of the wild-type strain. Chitin was stained as described by Okada et al. (2010). Briefly, yeast cells were harvested at 4000 rpm, rinsed twice with deionized water, and stained with 1 mg/ml calcofluor white (Sigma-Aldrich) for 1 min. Stained cells were washed once with deionized water and resuspended in water to be visualized by fluorescence microscopy using a DAPI filter (excitation band: 300-400 nm; emission band >420 nm).

## Compositional Analysis of the Cell Wall

Cells were mechanically prepared in a Braun homogenizer and processed as described elsewhere (Mora-Montes et al., 2007). Briefly, cell homogenates were centrifuged for 10 min at 20000 x g and 4°C, and the pellet was saved and washed five times with deionized water. Then, the cell walls were cleansed by serial incubations with hot SDS, β-mercaptoethanol, and NaCl and then hydrolyzed with 2 M trifluoroacetic acid (Sigma-Aldrich). Aliquots of 20 µl of acid-hydrolyzed samples were analyzed by high-performance anion exchange chromatography with pulsed amperometric detection in a Dionex system (Thermo Fisher Scientific) using separation conditions reported elsewhere (Estrada-Mata et al., 2016). To determine the protein concentration, cleansed walls were alkali-hydrolyzed as reported (Mora-Montes et al., 2007) before quantification using the Pierce BCA Protein Assay (Thermo Fisher Scientific).

The phosphomannan content was inferred from the cell's ability to bind the Alcian blue dye (Hobson et al., 2004). Upon growth, the cell concentration was adjusted to an OD₆₀₀nm of 0.2 in deionized water, and 1 ml aliquots were used to interact with 30 µg/ml Alcian blue (Sigma-Aldrich) and analyzed as described (Hobson et al., 2004). The cell wall mannan content was quantified essentially as previously reported (Navarro-Arias et al., 2016). For O-linked or N-linked mannan trimming, cells were incubated with 1 N NaOH and gently shaken for 18 h at room temperature or for 20 h at 37°C with 25 U endoglycosidase H (New England Biolabs). In both cases, cells were pelleted by centrifugation, and the supernatants were collected, lyophilized, and analyzed by high-performance anion exchange chromatography with pulsed amperometric detection (Mora-Montes et al., 2012).

## Determination of Lanthanide Contents in Yeast by Inductively Coupled Plasma Mass Spectrometry

Mutant and wild-type (BY4741) strains were grown in liquid YNB medium lacking inorganic phosphates and supplemented with β-glycerophosphate (1 mM). Cultures were inoculated at an OD₆₀₀nm of 0.05 (160 rpm; 28°C), and once the cultures
reached an OD<sub>600nm</sub> of 0.8, Ln were added at 50 µM and 6 µM La and Yb, respectively, concentrations corresponding to the EC<sub>10</sub> of the wild-type strain, and were exposed for 1 h. Cells were harvested at 4000 rpm for 1 min at 4°C. Pellets were washed three times with ice-cold MES (20 mM)—EDTA (10 mM) buffer (pH 6) followed by three washings with ice-cold ultrapure water and dried for 48 h at 70°C. Triplicate samples were prepared. Elemental extraction was performed on aliquots (125 mg) of each sample *via* mineralization with 1.75 ml HNO<sub>3</sub> and 0.5 ml hydrogen peroxide in closed tubes placed in a block digestion system (DigiPREP, SCP Sciences, Courtaboeuf, France). A gradual heating mode was used to achieve a final temperature of 100°C (total run of 265 min). Finally, ultrapure water was added to 12.5 ml, and filtration to 1 µm was performed. Concentrations were determined in triplicate using inductively coupled plasma spectrometry (ICP-MS, X Series II Model, Thermo Fischer Scientific, Courtaboeuf, France). The validity of the analytical method was checked by means of standard reference material (Oriental basma tobacco leaves, INCT-OBTL-5, LGC Promochem, Molsheim, France).

## Proteomics Analysis

The materials and procedures used for sample preparation for proteomics analysis, tandem mass spectrometry and peptide-to-spectrum assignment and protein identification are fully described in **Appendix A**.

## Data Analysis

Enriched biological pathways in Ln-stressed conditions vs. control conditions were identified using GSEA software v2.0 (Mootha et al., 2003; Subramanian et al., 2005). To establish a gene ranking, the sensitivity value difference was used. *p*-values for enriched gene sets were defined by computing 1,000 gene set permutations. *Saccharomyces cerevisiae* GO (biological process) derived from the MSigDB format gene set list was downloaded from the GO2Msig database and used as a template. From the GO2Msig datasets, gene set sizes were restricted to a maximum of 300 and a minimum of 10. GSEA results were visualized using Cytoscape v3.4 software (Shannon et al., 2003) and the Enrichment Map v2.1.0 (Merico et al., 2010) plug-in with default settings. In computed enrichment maps, node color corresponds to sensitivity of the mutant to Ln stress (yellow for resistant and blue for sensitive), and edge thickness represents redundancies between two gene sets.

Putative human homolog(s) of yeast gene(s) and any of their associated OMIM disease phenotypes were retrieved from the SGD YeastMine platform.<fnref n="1" /> Biological process and cellular compartment (MIPS cc) analyses were performed using FunSpec<fnref n="2" /> and evaluated for statistical significance (cutoff: value of *p*<0.01). Clustering (heatmaps) was performed using the R package pheatmap. Ln-responsive mutants found in this

study were compared with those from other studies carried out on *Saccharomyces cerevisiae* deletion mutant collections on Y (Grosjean et al., 2018), Cd (Jin et al., 2008; Ruotolo et al., 2008; Serero et al., 2008; Thorsen et al., 2009), Ni (Ruotolo et al., 2008; Arita et al., 2009; Bleackley et al., 2011), Zn (Pagani et al., 2007; Wang et al., 2007; Jin et al., 2008; Bleackley et al., 2011), Al (Kakimoto et al., 2005; Tun et al., 2014), As (Haugen et al., 2004; Thorsen et al., 2006, 2009; Jin et al., 2008; Johnson et al., 2016a), Cu (Wang et al., 2007; Jin et al., 2008; Jo et al., 2008; Bleackley et al., 2011), Fe (Wang et al., 2007; Jo et al., 2008; Bleackley et al., 2011), Mn (Wang et al., 2007; Bleackley et al., 2011; Chesi et al., 2012), Co (Bleackley et al., 2011), or Cr (Jin et al., 2008; Johnson et al., 2016b) and other stressors (Bennett et al., 2001; Chang et al., 2002; Aouida et al., 2004; Serrano et al., 2004; Thorpe et al., 2004).

## Statistical Analyses

For statistical analyses, one-way ANOVA and Tukey HSD post-hoc tests were used. Two-tailed Student's *t*-test was used where appropriate. All analyses were performed using R (version 3.4.1). All reported values are the means ± standard deviation (SD) of triplicates (n = 3).

## Data Availability

The authors declare that all data related to the findings of this study are available within the article and supplementary information or are available from the corresponding author upon reasonable request. Source data are provided with this paper. The mass spectrometry proteomics data for the *kex2Δ* mutant have been deposited to the ProteomeXchange Consortium *via* the PRIDE partner repository with the dataset identifier PXD025030 and project DOI 10.6019/PXD025030.<fnref n="3" /> The mass spectrometry and proteomics dataset for the WT are available under the dataset identifier PXD010700 and project DOI 10.6019/PXD010700.

## RESULTS AND DISCUSSION

### Yeast Genomic Phenotyping Unravels a Large Set of Lanthanide-Responsive Genes

The complete set of 4,733 haploid non-essential gene deletion mutants (BY4741, MATa) was screened under sublethal La and Yb concentrations. Mutants displaying a resistant or sensitive phenotype were further confirmed. Wild-type sublethal concentrations were used to confirm sensitive mutants, while a wild-type lethal concentrations were used to confirm resistant mutants. The primary and validation screens revealed genes involved either in resistance or sensitivity to Ln based on the modified growth of the mutants toward one or both elements. They were classified into six different categories, ranging from high sensitivity (HS) to high resistance

<footnote n="1">¹http://yeastmine.yeastgenome.org/yeastmine/</footnote>

<footnote n="2">²http://funspec.med.utoronto.ca</footnote>

<footnote n="3">³http://bioinformatics.org/go2msig/</footnote>
(HR; **Supplementary Table S1**). The validity and robustness of the screen were confirmed in the opposite mating type (BY4742, MATα) using 25 mutants randomly selected among sensitive and resistant mutants (**Supplementary Table S2**). A total of 630 genes were involved in the response of yeast cells to Ln, with 1.5 times more resistant (366) than sensitive mutants (236) to at least one Ln. Three times less mutants were identified in a recent fitness analysis of a pooled mutant population performed with the diploid BY4743 strain (Pallares et al., 2021). This lower number of mutants certainly arose from the fact that low concentrations (EC20) of La and Yb were used suggesting that more toxic concentrations is more appropriate to reveal the entire set of genes involved in the tolerance to Ln. The sublethal/lethal concentrations used here reproduce conditions usually applied for metal toxicogenomics in yeast, resulting in a similar number of sensitive (Ruotolo et al., 2008) and resistant (Ruotolo et al., 2008; Grosjean et al., 2018) mutants.

## La and Yb-Responsive Genes Are Mostly Dissimilar

We identified 259 and 211 genes whose deletion caused resistance to La and Yb, respectively (Figure 1A). Five times more mutants were sensitive to La (257) than to Yb (46), suggesting additional resistance mechanisms toward La. A relatively small proportion of mutants were sensitive (13%, 39 mutants) or resistant (16%, 73 mutants) to both Ln, which reveals the existence of a core, yet minor, set of genes involved in a general Ln response. This core set of mutants is more important than the one identified when using EC20 concentrations of La and Yb (<2%) for the library screen (Pallares et al., 2021). However, in our study, a set of mutants exhibited an opposite phenotype to the two Ln, with 26 mutants being La-sensitive and Yb-resistant, while two displayed a La-resistant and Yb-sensitive phenotype (Figure 1A). Similar proportions were found when considering only the mutants displaying the strongest resistant/sensitive phenotypes (Figure 1B; Supplementary Table S3). Altogether, our data strongly highlight the dissimilarity of responses toward La and Yb in yeast.

## Yeast Lanthanide-Responsive Genes Have Counterparts Whose Mutation Trigger Human Diseases

Human homologues and putative-related diseases were retrieved from the YeastMine database for the set of mutants identified. At least one human homolog was found for 62% (391 mutants) of the knockout genes in mutants resistant or sensitive to Ln (**Supplementary Table S1**), a third of which are involved in cancers and other human diseases. Ln toxicity concerns are growing after several studies reporting Gd accumulation in bone tissues and more recently in the brain (Ramalho et al., 2016; Gulani et al., 2017). Therefore, the modified fitness toward Ln toxicity observed in mutants for genes associated with bone mineral deficiencies (SAC6) or associated with cerebral diseases (e.g., ELO3, PHO85) emphasizes that more attention should be given to Ln toxicity. The high number of yeast Ln-responsive

mutants for genes holding human counterparts sets a promising basis for further studies on the impact of Ln toxicity on human health.

## Lanthanide Responses Strongly Differ From Well-Known Transition Metal-Mediated Responses

We performed an *in silico* cross-comparison of Ln-responsive mutants with similar studies carried out with other metals and the metalloid As (Figure 2; Supplementary Table S4). Mutants were grouped into eleven clusters (C1–C11) based on their phenotypes, distinguishing three major patterns.

The first pattern groups REE-specific responsive mutants with 32% mutants that responded specifically to La (17% sensitive (C6) or resistant (C8)), to both Yb and Y (4%, C3), or to La, Yb, and Y (5%, C4). These mutations target genes involved in pathways, such as the response to pheromone-related GO biological processes (STE gene family) and chromatin silencing mutants. Twelve mutants (2%, C5) specific to Yb were involved in cell wall biogenesis. Several mutants for actin-related functions and lipid metabolic processes were specific to La (C6–C7).

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Lanthanide response patterns obtained from genomic phenotyping of a whole mutant collection of *Saccharomyces cerevisiae*. (A) Venn diagram highlighting the number of mutants being sensitive or resistant to La and/or to Yb. (B) Venn diagram highlighting the number of mutants displaying an exacerbated phenotype (highly resistant/sensitive) toward La and/or Yb.</figcaption></figure>
The second pattern corresponds to the C2 and C11 clusters. Eleven mutants (2%, C2), which are related to the ubiquitin-dependent protein catabolic process *via* the multivesicular body sorting pathway (ESCRT complex), were specifically resistant to Ni, Y, and both Ln but sensitive to most of the other metals. Most of the C11 cluster mutants were also resistant to Ni, Y, and Yb but sensitive to La and were related to the retromer complex.

The third pattern included 139 mutants (22%, clusters C1, C9, and C10) sensitive to a range of metals, including REEs, whose functions are likely linked to the general metal stress response rather than Ln-specific, such as vesicular transport systems, as well as the vacuole and its acidification.

Last, a close relationship between Ln and Y emerged. Most of the mutants we identified (46%) were also either resistant or sensitive to Y (Figure 2). Y and Yb clustered closer than the two Ln together, which could be due to the close chemical properties of Y and Yb and therefore explain common cellular functions. The identification of many specific mutants in the Yb, La, or both screens underscores the unique cellular impact of Ln over other well-known metals.

## Ln Responses Differ From Well-Known Alkaline pH and Oxidative Stress-Mediated Responses

Oxidative stress is a widespread noxious consequence of exposure to some metals and metalloids, such as Cd and As, in yeast (Brennan and Schiestl, 1996; Litwin et al., 2013). Since similarities

between Ln and As (32%) and between Ln- and Cd-responsive mutants (27%) were found (Figure 2), we further assessed whether Ln-mediated responses were common to oxidative stress-generating stressors in yeast. We compared our data with previous screens carried out on various chemical and physical stressors (Bennett et al., 2001; Chang et al., 2002; Aouida et al., 2004; Serrano et al., 2004; Thorpe et al., 2004) (Supplementary Figure S2; Supplementary Table S5). The highest similarity was found for both bleomycin and diamide but accounted for only 17% of the mutants responding to Ln and were mainly linked to genes involved in protein transport and transcription (Supplementary Figure S2). Although previous studies mentioned oxidative stress caused by Ln (Pagano et al., 2015), La and Yb clustered together as an outgroup, with more than 360 mutants specific to Ln, suggesting that the response to Ln stress is not primarily due to indirect stresses, such as pH or oxidative stress.

## The Golgi, ER, and Vacuole Are the Major Cellular Compartments Affected by Lanthanides

To highlight cellular hotspots of Ln responses, MIPS (Munich Information Center for Protein Sequences) subcellular localization analysis was performed using the identified mutants as a query (Supplementary Table S5). The disruption of genes linked to the vacuole and to the Golgi, which includes the Golgi membrane, Golgi to vacuole or to ER vesicle transport, led to cell sensitivity

<figure><img src="image_4.png" /></figure>
to both Ln, similar to Ni and Cd (Ruotolo et al., 2008), while some endosomal proteins were involved in Ln toxicity (resistant mutants). La-specific compartments were highlighted, with mutants for genes associated with the cytoplasm. Moreover, mutants for proteins of the vacuolar membrane and the nuclear envelope were sensitive to La, while impairing the actin cytoskeleton led to La resistance (Supplementary Table S6). Last, mutants for proteins of the ER membrane and the plasma membrane were specific for Yb. The lack of representation of the mitochondria was noticeable, which could be related to the low similarity between the Ln-mediated stress and chemical-induced oxidative stress.

## La and Yb Trigger Different Cellular Functions

To deepen the understanding of Ln-mediated toxicity mechanisms, the identification of the cellular processes engaged in Ln toxicity was evaluated through the analysis of biological pathways (Figure 3; Supplementary Table S7). Transport and localization functions, which include vacuolar transport, endocytosis and vesicle-mediated transport, catabolic processes, and lipid metabolism, were highly represented among La- and Yb-responsive mutants (Figure 3). Previously, endocytosis was also found as a primary pathway affected by these metals (Pallares et al., 2021) together with lipids, amino acids, and carbohydrates metabolisms (Grosjean et al., 2021).

However, relatively contrasting enrichment maps were generated for Yb and La (Figure 3), supporting that the two Ln involve diverse cellular functions in yeast. Functional categories, such as ion homeostasis, ribosome biogenesis, and ncRNA processing, were mainly specific to La. Other functions were less represented for Yb than for La, including mitotic and meiotic processes, together with metabolic processes involving organic acids and organo-nitrogen compounds. Several functions when altered caused sensitivity to La but resistance to Yb, such as cell wall organization and biogenesis and stress response pathways (Figure 3). Only a few functions enriched were mostly Yb-specific, such as signaling pathways, stress response, and catabolic processes (Supplementary Table S7). A wide range of mutants for proteins involved in translation were strictly sensitive to La, illustrated by cells lacking subunits of the ribosomal 60S and 40S subunit proteins (RPS and RPL proteins) (Figure 3; Supplementary Table S1). Given that many more mutants are sensitive to La (Figure 1), more La-responsive proteins could be required under La stress compared to Yb. Therefore, a reduced translation efficiency would be more detrimental under La stress.

This analysis allowed us to assign mutants to major functions involved in the response to Ln. This highlighted the components of the cell responses to La and Yb stress that belong to (i) the core component responsive to general Ln stress and (ii) the specific pathways involved in either La or Yb stress modulation. To further deepen our knowledge on the dissimilar response of yeast to La and Yb, we next focused our investigation on the role played by the cell wall and related signaling pathway in response to Ln.

## La and Yb Differentially Affect the Cell Wall and Related Signaling Pathway

The cell wall is the first shield against extracellular stressors, including metals. Several mutants responding to Ln exposure were knocked out for genes encoding proteins involved in the cell wall stress response and cell wall biosynthesis (Figures 3, 4; Supplementary Table S7). In a recent transcriptomics and proteomics analysis of S. cerevisiae, cell wall biosynthesis and organization-related pathways were upregulated in response to both La and Yb, supporting the strong implication of the cell wall as the first interface (Grosjean et al., 2021). However, the mutants found in the present deletome analysis displayed distinct phenotypes under La or Yb exposure, supporting the need of complementary multi-scale approaches. Mutants for mannosyltransferases (mnn10Δ, mnn9Δ, mnn2Δ, and anp1Δ) involved in cell wall mannoprotein synthesis were specifically sensitive to La, whereas mutants for glucosyltransferases involved in glycoprotein synthesis (alg3Δ, alg6Δ, and alg8Δ) were specifically resistant to Yb (Figure 4; Supplementary Table S1).

Remarkably, most of the cell wall-related mutants were resistant to Yb but sensitive to La, such as mutants in cell wall integrity-related proteins Kre1p, Kre6p, and Gas1p (Figure 4, 5; Supplementary Table S1). Mutants *kex2Δ* and *kex1Δ* also displayed this contrasted phenotype (Figures 4, 5), and Kex2p overexpression in WT cells reversed the phenotype (Figure 5B). Kex2p and Kex1p are both Ca-dependent serine proteases involved in the activation of secretory pathway and cell wall proteins (Figure 4), such as Scw4p (Fuller et al., 1989; Lesage and Bussey, 2006; Grbavac et al., 2017). Accordingly, exposure of WT yeast to the serine protease inhibitor tosyl phenylalanyl chloromethyl ketone (TPCK) resulted in phenotypes similar to those obtained with *kex2Δ* and *kex1Δ* (Figure 5C). We also found that Mkc7p, an aspartyl protease and member of the yapsin family acting downstream of Kex2p (Grbavac et al., 2017), was resistant to Yb only (Figure 4), supporting the importance of cell wall maintenance under Ln stress.

## Cell Wall Modifications Differentially Influence Ln Accumulation

We assessed whether cell wall integrity components had a role in different Ln accumulation by measuring La and Yb concentrations in WT and mutant cells (Figures 6A,B). Interestingly, *kre1Δ*, *kre6Δ*, *kex1Δ*, and *kex2Δ* accumulated two to three times less Yb (Figure 6B) than the WT, correlating with their resistance to Yb. Conversely, *kex2Δ* and *kre6Δ* accumulated significantly more La (Figure 6A), concordant with their higher sensitivity to La (Figure 5A). Consequently, these proteins and related cell wall components could be involved in (i) Ln differential adsorption followed by an internalization step or (ii) by affecting Ln uptake through the modulation of ion membrane transporters.

## La and Yb Exposure Modify Cell Wall Structure and Composition

To investigate whether cell wall composition modulates Ln interactions, the cell wall composition of WT, *kex2Δ*, and *kre6Δ*
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Functional enrichment analysis network of functions involved in the cellular response to lanthanides. Functions whose deletion renders cells either sensitive or resistant to La (outer circle) or Yb (inner circle) are shown in blue and yellow, respectively. Interlines represent gene overlap between two related functions, and edge width is proportional to the number of shared genes. The enrichment map was built using GSEA and visualized by the Enrichment map plug-in in Cytoscape.</figcaption></figure>

was analyzed under Ln exposure (Figures 6C-E; Supplementary Figure S3). The chitin content of WT cells remained constant despite Ln exposure. Interestingly, upregulation of the chitin cell wall metabolic process and the 1,3-β-D-glucan biosynthesis pathways were observed by transcriptomics and proteomics (Grosjean et al., 2021). This suggests that upregulation of these cell wall-related pathways is essential to maintain a proper cell wall structure, composition, and organization to compensate for Ln binding. Relatively negligible differences in β-1,3-glucan and cell wall proteins were observed between the different conditions and strains (Supplementary Figures S3B,D). However, independent of the absence or presence of Ln, the

chitin content at least doubled in both *kex2Δ* and *kre6Δ* relatively to WT, which supports that these cell wall modifications are a consequence of these mutations (Figures 6C,D; Supplementary Figure S3A). Additionally, while the amount of mannans remained relatively unchanged, lower amounts of phosphomannans and O-linked mannans were observed in the mutants (Figure 6E; Supplementary Figures S3C,E,F). A previous study demonstrated the adsorption of Yb and subsequent formation of Yb-phosphate crystals onto the cell wall of *S. cerevisiae* (Jiang et al., 2012). Other studies mentioned the adsorption of Ln on extracellular biopolymers, such as the cell wall of certain bacteria (Moriwaki and Yamamoto, 2013). It was further argued that
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Schematic representation of the proteins belonging to the cell wall signaling pathway and cell wall organization in the response to Ln. The data show the dual involvement of the cell wall and signalization in the response to lanthanides.</figcaption></figure>

LREEs (La) preferentially bind to phosphate groups, while HREEs (Yb) equally bind to phosphate and carboxyl groups (Ngwenya et al., 2009). Consequently, the relative proportions of chitin and phosphomannan could modulate the proportion of functional groups to which Ln could differentially bind. The lower abundance of phosphate groups containing phosphomannans could reduce La adsorption, therefore increasing the free La accumulation potential and toxicity (Figure 6A). Conversely, reduced phosphomannans combined with higher chitin content would increase the proportion of carboxyl groups, allowing greater Yb adsorption synonymous to an increased Yb resistance. To confirm the hypothesis of a selective LREE and HREE binding potential, we evaluated the growth of *kex2Δ*, *kex1Δ*, *kre1Δ*, and *kre6Δ* with 14 different Ln (Figure 5A). We observed that *kex2Δ* and, to a lesser extent, *kex1Δ*, *kre1Δ*, and *kre6Δ* were all sensitive to LREEs (La to Eu), while being gradually resistant with heavier Ln. These observations support the differential interaction of LREEs and HREEs with the cell wall based on their chemical properties.

## Cell Wall Modifications Disturb Protein Abundances

Since *kex2Δ* displayed the most contrasting phenotype, we leveraged label-free high-throughput shotgun proteomics to identify proteome modulations under Ln exposure (Figure 7;

**Supplementary Figure S4; Supplementary Table S8).** Seventy-five proteins displayed different abundances between La and Yb in *kex2Δ*, and twelve interesting candidates, because of their biological role, are shown in Figure 7. Nqm1p and Ald3p were both repressed under Yb but more abundant under La. Nqm1p is a transaldolase upregulated by various stresses, including osmotic shock (Michel et al., 2015), while Ald3p is a stress-responsive aldehyde dehydrogenase (Navarro-Aviño et al., 1999). Upregulation of these two proteins under La could be the result of osmotic stress induced by important La absorption in *kex2Δ* (Figure 6A). Interestingly, several membrane transporters had a lower abundance in *kex2Δ* than in WT. Pdr5p, a plasma membrane ABC multidrug transporter involved in cation resistance (Miyaharan et al., 1996), was specifically repressed under Yb (Figure 7). Similarly, Pho84p, a high-affinity phosphate transporter and low-affinity Mn transporter, and Ftr1p, a high-affinity iron transporter (Figure 7), were less abundant in *kex2Δ*. We can hypothesize that Yb could be taken up through these transporters, which would explain the higher resistance of *kex2Δ* to Yb by decreasing its putative non-specific uptake.

## La and Yb Trigger the Cell Wall-Related Signaling Pathway

This genome-wide deletion mutant screen further highlighted mutants for genes related to signaling pathway (Figure 4). Again, a segregation between La and Yb was observed. The cell wall
A

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

B

<figure><img src="image_10.png" /></figure>

C

<figure><img src="image_11.png" /></figure>

**FIGURE 5** | Growth of Kex and Kre mutants and effects of Kex2p overexpression or TPCK supplementation on yeast growth under Ln exposure. (A) Yeast growth was assessed on YPD medium without lanthanides (control) or supplemented with 4.0mM La, 3.8mM Ce, 5.8mM Pr, 4.2mM Nd, 3.9mM Sm, 3.9mM Eu, 3.8mM Gd, 3.8mMTb, 3.6mM Dy, 3.9mM Ho, 3.9mM Er, 3.6mM Tm, 3.6mM Yb, or 3.6mM Lu, with 10-fold serial dilutions of cultures from left to right in each panel. A representative plate (out of 3 independent experiments) is shown. Plates were incubated for 5 days at 28°C. (B) Overexpression of Kex2p confers resistance to La but sensitivity to Yb. Yeast cells were transformed with the empty plasmid (pYES2) or with the same plasmid harboring KEX2 (pYES2-KEX2). Cells were grown on YPG in the presence of either La (3.8mM) or Yb (3.2mM) or in the absence of Ln (control), and 10-fold serial dilutions were plated on YPD media. (C) Drop test of 3-fold serial dilutions of the WT strain exposed to La (3.8mM) and Yb (3.2mM) with the addition of the serine protease inhibitor tosyl phenylalanyl chloromethyl ketone (TPCK). Representative pictures (out of 3 independent experiments) are given.

stress signaling pathway was represented notably by the plasma
membrane receptor *mid2Δ* and the GTP exchange factor *rom2Δ*,
both of which are resistant to La (**Figure 4**). Additionally, *Zeo1p*,
which acts antagonistically to *Mid2p*, was significantly more
abundant in *kex2Δ* (**Figure 7**). *slg1Δ*, a cell wall integrity sensor
(Mao *et al.*, 2011), and *fks1Δ* were also resistant to both Ln.
Furthermore, *sac7Δ*, a negative regulator of the RHO1-PKC1-
MAPK cell integrity pathway, displayed the same phenotype as
*kreΔ* and *kexΔ*, suggesting the dual role of the cell wall and its
ability to respond to different Ln (**Figure 4**). Altogether, these

results suggest that the cell wall not only acts as a Ln sequestration/
adsorption compartment but is also an active signal effector
differentiating LREEs and HREEs.

CONCLUSION

To date, studies dealing with Ln toxicity are scarce, and how
organisms can cope with the toxicity of these emerging contaminants
remains an open question. In the present study, we adopted a
<figure><img src="image_12.png" /><figcaption>**FIGURE 6** | La and Yb concentrations and cell wall composition of WT and mutant strains. (A,B) Cells were grown in modified YNB and exposed to (A) La (50 µM) or (B) Yb (6 µM) for one hour. Data are the means (±SD) of three independent cultures. Significant differences from the wild-type condition are indicated by asterisks (ANOVA, Tukey HSD). (C) Cell wall composition of WT, kex2Δ, and kre6Δ in chitin, glucans, and mannans. Significant differences from control conditions in each strain are indicated by asterisks (t-test). (D) Chitin content in the cell wall of WT, kex2Δ, and kre6Δ under Ln exposure. Significant differences from the WT strain are indicated by asterisks (t-test). (E) Phosphomannan content using the Alcian blue binding test in each strain. Significant differences from the wild-type strain are indicated by asterisks (t-test). Data are the means (±SD) of three independent cultures. For all experiments: \*\*\*<0.05, \*\*\*\*<0.01, \*\*\*\*\*<0.001.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**FIGURE 7** | Heatmap displaying protein abundance changes in kex2Δ compared to the WT. Please, refer to the Supplementary text for information related to proteomics data acquisition and analysis.</figcaption></figure>

well-defined metal toxicogenomics analysis in which yeast were subjected to a high Ln-mediated stress comparatively to a recent study (Pallares et al., 2021). Interestingly, we confirmed in our study that the vesicle transport and endocytosis pathways were strongly affected by Ln. However, we report several other pathways

that were also impacted by a high Ln stress, including lipid metabolism, ion homeostasis, cell wall organization, and related signaling pathway (Figure 8), which were not reported previously (Pallares et al., 2021), extending our knowledge on the cellular and molecular responses of yeast exposed to Ln. Additionally, by using an arrayed mutant collection screening, we demonstrated that some mutants that were previously found specific to HREEs, are also shared with La (e.g., ypl056cΔ, ykr020wΔ, ypl057cΔ, ybr036cΔ), expending the yeast core response to Ln, but also supporting the need for combining experimental approaches. In a recent paper, we combined transcriptomics and proteomics analyses on *Saccharomyces cerevisiae* exposed to La and Yb and most of the pathways impacted by Ln were found in each of these complementary approaches (Grosjean et al., 2021). This multi-scale analysis (deletome, transcriptome, and proteome) of the response of yeast to Ln provides valuable target genes, proteins, and pathways involved in Ln toxicity response and detoxification that complements and expands on previous studies.

Moreover, we highlighted that the response pathways were modulated differently depending on the nature of the different Ln. The high reactivity of Ln to functional groups of the cell
<figure><img src="image_14.png" /><figcaption>**FIGURE 8** | Global pathways in *Saccharomyces cerevisiae* exposed to lanthanide stress. Schematic representation of general functions identified through mutant collection screening. Common functions between the two lanthanides are shown in red, and functions showing a different behavior are shown in blue. La is representative of LREEs, while Yb is representative of HREEs. La-P and Yb-P represent lanthanum ions bound to phosphate groups, while Yb-C represents ytterbium ions bound to carboxyl groups.</figcaption></figure>

wall makes this compartment the first layer discriminating different Ln. The cell wall stress-related signaling pathway might induce the modification of cell wall organization as well as the lipid composition of membranes in response to Ln (Figure 8). Such a close relationship between membrane lipids and adaptation to environmental stress has been shown to impact Ca influx channels (Jiang et al., 2019), by which Ln could be taken up (Ene et al., 2015). Additionally, Ln activate intracellular signaling, accumulation of phosphatidylinositol-4-phosphate at the plasma membrane, and cytoskeleton-dependent rearrangement complexes in Arabidopsis (Lee et al., 2020). Therefore, a different mechanism to regulate Ln membrane transporters should also be considered, such as through endocytosis (Figure 8). The recurrence of Mn, Zn, and Fe transporters in our deletomic and proteomic analyses are in concordance with our previous whole cell transcriptomics and proteomics analysis on wild-type S. *cerevisiae*, for which we observed inhibition of gene expression and reduced protein abundance of these transporters (e.g., Ftr1-Fet3, Zrt1,

Smf1; Grosjean et al., 2021). The identification of these cation transporters by several multi-scale and complementary approaches also calls for further investigations on these uptake pathways as a way of Ln entry in cells and their involvement in the disruption of micronutrient homeostasis.

Given the differences observed between LREEs and HREEs, engineering yeast cells to optimize REE binding (modified cell wall) and/or accumulation deserves further attention. Yeasts could be used as a biotechnology tool for REE purification in mixed metal or Ln solutions (Sun et al., 2019). Finally, in their environment, organisms including humans may be exposed to a single REE species or to a mixture of LREEs and HREEs. Considering the dissimilarities observed between La and Yb in our study, exposure to mixed Ln could display additive or antagonistic effects (Morel et al., 2021), which would require further investigation.

These results on S. *cerevisiae* are a first step in the delineation of Ln toxicity mechanisms and will enable future investigations on more complex eukaryotes, to improve risk assessment toward
Ln. The high number of yeast Ln-responsive mutants for genes holding human counterparts sets a promising basis for further studies on the impact of Ln toxicity on human health.

## DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

## AUTHOR CONTRIBUTIONS

NG: investigation, conceptualization, formal analysis, writing—original draft, and writing—review and editing. MLJ: conceptualization, formal analysis, supervision, and writing—review and editing. MC, HM-M, and JA: investigation, formal analysis, and writing—review and editing. EG: supervision and writing—review and editing. DB: investigation, conceptualization, formal analysis, supervision, project funding, and writing—review and editing. All authors contributed to the article and approved the submitted version.

## REFERENCES

Aouida, M., Pagé, N., Leduc, A., Page, N., Peter, M., and Ramotar, D. (2004). A genome-wide screen in Saccharomyces cerevisiae reveals altered transport as a mechanism of resistance to the anticancer drug bleomycin. *Cancer Res.* 64, 1102–1109. doi: 10.1158/0008-5472.CAN-03-2729

Arita, A., Zhou, X., Ellen, T. P., Liu, X., Bai, J., Rooney, J. P., et al. (2009). A genome-wide deletion mutant screen identifies pathways affected by nickel sulfate in Saccharomyces cerevisiae. *BMC Genomics* 10, 1–14. doi: 10.1186/1471-2164-10-524

Bennett, C. B., Lewis, L. K., Karthikeyan, G., Lobachev, K. S., Jin, Y. H., Sterling, J. F., et al. (2001). Genes required for ionizing radiation resistance in yeast. *Nat. Genet.* 29, 426–434. doi: 10.1038/ng778

Bleackley, M. R., Young, B. P., Loewen, C. J. R., and Macgillivray, R. T. A. (2011). High density array screening to identify the genetic requirements for transition metal tolerance in Saccharomyces cerevisiae. *Metallomics* 3, 195-205. doi: 10.1039/c0mt00035c

Botstein, D., Chervitz, S. A., and Cherry, J. M. (1997). Yeast as a model organism. *Science* 277, 1259–1260. doi: 10.1126/science.277.5330.1259

Brennan, R. J., and Schiestl, R. H. (1996). Cadmium is an inducer of oxidative stress in yeast. *Mutat. Res. - Fundam. Mol. Mech. Mutagen.* 356, 171–178. doi: 10.1016/0027-5107(96)00051-6

Chang, M., Bellaoui, M., Boone, C., and Brown, G. W. (2002). A genome-wide screen for methyl methanesulfonate- sensitive mutants reveals genes required for S phase progression in the presence of DNA damage. *Proc. Natl. Acad. Sci. U. S. A.* 99, 16934–16939. doi: 10.1073/pnas.262669299

Chesi, A., Kilaru, A., Fang, X., Cooper, A. A., and Gitler, A. D. (2012). The role of the Parkinson's disease gene PARK9 in essential cellular pathways and the manganese homeostasis network in yeast. *PLoS One* 7:e34178. doi: 10.1371/journal.pone.0034178

Ciacci, L., Reck, B. K., Nassar, N. T., and Graedel, T. E. (2015). Lost by design. *Environ. Sci. Technol.* 49, 9443–9451. doi: 10.1021/es505515z

Cotruvo, J. A. Jr. (2019). The chemistry of lanthanides in biology: recent discoveries, emerging principles, and technological applications. *ACS Cent. Sci.* 5, 1496–1506. doi: 10.1021/acscentsci.9b00642

Ene, C. D., Ruta, L. L., Nicolau, I., Popa, C. V., Iordache, V., Neagoe, A. D., et al. (2015). Interaction between lanthanide ions and Saccharomyces cerevisiae cells. *J. Biol. Inorg. Chem.* 20, 1097–1107. doi: 10.1007/s00775-015-1291-1

## FUNDING

This work was supported by the French National Research Agency through the National Program "Investissements d'Avenir" with the reference ANR-10-LABX-21-01/LABEX RESSOURCES 21 and by the Region Grand Est. HM-M is supported by Consejo Nacional de Ciencia y Tecnología, México (ref. CF-2019-6380).

## ACKNOWLEDGMENTS

We acknowledge Nadia Morin-Crini and Caroline Amiot of the PEAT² platform for the ICP-MS analyses, Luz A. López-Ramírez for the cell wall analyses, and Adam Schikora for critical review of the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.881535/ full#supplementary-material

Estrada-Mata, E., Navarro-Arias, M. J., Pérez-García, L. A., Mellado-Mojica, E., López, M. G., Csonka, K., et al. (2016). Members of the Candida parapsilosis complex and Candida albicans are differentially recognized by human peripheral blood mononuclear cells. *Front. Microbiol.* 6:1527. doi: 10.3389/fmicb.2015.01527

Fuller, R. S., Brake, A., and Thorner, J. (1989). Yeast prohormone processing enzyme (KEX2 gene product) is a Ca2+-dependent serine protease. *Proc. Natl. Acad. Sci.* 86, 1434–1438. doi: 10.1073/pnas.86.5.1434

Goecke, F., Vítová, M., Lukavský, J., Nedbalová, L., Řezanka, T., and Zachleder, V. (2017). Effects of rare earth elements on growth rate, lipids, fatty acids and pigments in microalgae. *Phycol. Res.* 65, 226–234. doi: 10.1111/pre.12180

Gonzalez, V., Vignati, D. A. L., Leyval, C., and Giamberini, L. (2014). Environmental fate and ecotoxicity of lanthanides: are they a uniform group beyond chemistry? *Environ. Int.* 71, 148–157. doi: 10.1016/j.envint.2014.06.019

Good, N. M., Vu, H. N., Suriano, C. J., Subuyuj, G. A., Skovran, E., and Martinez-Gomez, N. C. (2016). Pyrroloquinoline quinone ethanol dehydrogenase in Methylobacterium extorquens AM1 extends lanthanide-dependent metabolism to multicarbon substrates. *J. Bacteriol.* 198, 3109–3118. doi: 10.1128/JB.00478-16

Grbavac, A., Čanak, I., Stuparević, I., Teparić, R., and Mrša, V. (2017). Proteolytic processing of the Saccharomyces cerevisiae cell wall protein Scw4 regulates its activity and influences its covalent binding to glucan. *Biochim. Biophys. Acta, Mol. Cell Res.* 1864, 507–515. doi: 10.1016/j.bbamcr.2016.12.009

Grosjean, N., Blaudez, D., Chalot, M., Gross, E. M., and Le Jean, M. (2020). Identification of new hardy ferns that preferentially accumulate light rare earth elements: a conserved trait within fern species. *Environ. Chem.* 17, 191–200. doi: 10.1071/EN19182

Grosjean, N., Gross, E. M., Le Jean, M., and Blaudez, D. (2018). Global deletome profile of Saccharomyces cerevisiae exposed to the technology-critical element yttrium. *Front. Microbiol.* 9, 1–13. doi: 10.3389/fmicb.2018.02005

Grosjean, N., Le Jean, M., Armengaud, J., Schikora, A., Chalot, M., Gross, E. M., et al. (2021). Combined omics approaches reveal distinct responses between light and heavy rare earth elements in Saccharomyces cerevisiae. *J. Hazard. Mater.* 425:127830. doi: 10.1016/J.JHAZMAT.2021.127830

Grosjean, N., Le Jean, M., Berthelot, C., Chalot, M., Gross, E. M., and Blaudez, D. (2019). Accumulation and fractionation of rare earth elements are conserved traits in the Phytolacca genus. *Sci. Rep.* 9:18458. doi: 10.1038/s41598-019-54238-3

Gulani, V., Calamante, F., Shellock, F. G., Kanal, E., and Reeder, S. B., and others (2017). Gadolinium deposition in the brain: summary of evidence
and recommendations. *Lancet Neurol.* 16, 564-570, doi:10.1016/S1474-4422(17)30158-8.

Haugen, A. C., Kelley, R., Collins, J. B., Tucker, C. J., Deng, C., Afshari, C. A., et al. (2004). Integrating phenotypic and expression profiles to map arsenic-response networks. *Genome Biol.* 5:R95. doi: 10.1186/gb-2004-5-12-r95

He, M. L., Ranz, D., and Rambeck, W. A. (2001). Study on the performance enhancing effect of rare earth elements in growing and fattening pigs. *J. Anim. Physiol. Anim. Nutr. (Berl).* 85, 263-270. doi: 10.1046/j.1439-0396.2001.00327.x

Hobson, R. P., Munro, C. A., Bates, S., MacCallum, D. M., Cutler, J. E., Heinsbroek, S. E. M., et al. (2004). Loss of cell wall mannosylphosphate in Candida albicans does not influence macrophage recognition. *J. Biol. Chem.* 279, 39628-39635. doi: 10.1074/jbc.M405003200

Hu, A., Guo, J.-J., Pan, H., and Zuo, Z. (2018). Selective functionalization of methane, ethane, and higher alkanes by cerium photocatalysis. *Science* 361, 668-672. doi: 10.1126/science.aat9750

Hu, Z., Richter, H., Sparovek, G., and Schnug, E. (2004). Physiological and biochemical effects of rare earth elements on plants and their agricultural significance: a review. *J. Plant Nutr.* 27, 183-220. doi: 10.1081/PLN-120027555

Jiang, M. Y., Ohnuki, T., Tanaka, K., Kozai, N., Kamiishi, E., and Utsunomiya, S. (2012). Post-adsorption process of Yb phosphate nano-particle: formation by Saccharomyces cerevisiae. *Geochim. Cosmochim. Acta* 93, 30-46. doi: 10.1016/j.gca.2012.06.016

Jiang, Z., Zhou, X., Tao, M., Yuan, F., Liu, L., Wu, F., et al. (2019). Plant cell-surface GIPC sphingolipids sense salt to trigger Ca2+ influx. *Nature* 572, 341-346. doi: 10.1038/s41586-019-1449-z

Jin, Y. H., Dunlap, P. E., McBride, S. J., Al-Refai, H., Bushel, P. R., and Freedman, J. H. (2008). Global transcriptome and deleteme profiles of yeast exposed to transition metals. *PLoS Genet.* 4:e1000053. doi: 10.1371/journal.pgen.1000053

Jo, W. J., Loguinov, A., Chang, M., Wintz, H., Nislow, C., Arkin, A. P., et al. (2008). Identification of genes involved in the toxic response of Saccharomyces cerevisiae against iron and copper overload by parallel analysis of deletion mutants. *Toxicol. Sci.* 101, 140-151. doi: 10.1093/toxsci/kfm226

Johnson, A. J., Veljanoski, F., O'Doherty, P. J., Zaman, M. S., Petersingham, G., Bailey, T. D., et al. (2016a). Molecular insight into arsenic toxicity via the genome-wide deletion mutant screening of Saccharomyces cerevisiae. *Metallomics* 8, 228-235. doi: 10.1039/c5mt00261c

Johnson, A. J., Veljanoski, F., O'Doherty, P. J., Zaman, M. S., Petersingham, G., Bailey, T. D., et al. (2016b). Revelation of molecular basis for chromium toxicity by phenotypes of Saccharomyces cerevisiae gene deletion mutants. *Metallomics* 8, 542-550. doi: 10.1039/C6MT00039H

Kachroo, A. H., Laurent, J. M., Yellman, C. M., Meyer, A. G., Wilke, C. O., and Marcotte, E. M. (2017). Systematic humanization of yeast genes reveals conserved functions and genetic modularity. *Science* 348, 921-925. doi: 10.1016/j.antiviral.2015.06.014.Chronic

Kakimoto, M., Kobayashi, A., Fukuda, R., Ono, Y., Ohta, A., and Yoshimura, E. (2005). Genome-wide screening of aluminum tolerance in Saccharomyces cerevisiae. *Biometals* 18, 467-474. doi: 10.1007/s10534-005-4663-0

Kim, H. S., Lee, K., Kim, S. J., Cho, S., Shin, H. J., Kim, C., et al. (2018). Arrayed CRISPR screen with image-based assay reliably uncovers host genes required for coxsackievirus infection. *Genome Res.* 28, 859-868. doi: 10.1101/gr.230250.117

Lee, E., Santana, B. V. N., Samuels, E., Benitez-Fuente, F., Corsi, E., Botella, M. A., et al. (2020). Rare earth elements induce cytoskeleton-dependent and PI4P-associated rearrangement of SYT1/SYT5 endoplasmic reticulum-plasma membrane contact site complexes in Arabidopsis. *J. Exp. Bot.* 71, 3986-3998. doi: 10.1093/jxb/eraa138

Lee, J. C. K., and Wen, Z. (2018). Pathways for greening the supply of rare earth elements in China. *Nat. Sustain.* 1, 598-605. doi: 10.1038/s41893-018-0154-5

Leeuw, D., and Academic, K. (2000). Inorganic yellow-red pigments without toxic metals. *Nature* 404, 980-982. doi: 10.1038/35010082

Lesage, G., and Bussey, H. (2006). Cell wall assembly in Saccharomyces cerevisiae. *Microbiol. Mol. Biol. Rev.* 70, 317-343. doi: 10.1128/MMBR.00038-05

Litwin, I., Bocer, T., Dziadkowiec, D., and Wysocki, R. (2013). Oxidative stress and replication-independent DNA breakage induced by arsenic in Saccharomyces cerevisiae. *PLoS Genet.* 9:e1003640. doi: 10.1371/journal.pgen.1003640

Liu, W., Guo, M., Liu, C., Yuan, M., Chen, X., Huot, H., et al. (2018). Water, sediment and agricultural soil contamination from an ion- adsorption rare earth mining area. *Chemosphere* 10, 75-83. doi:10.1016/j.chemosphere.2018.10.109

Mao, K., Wang, K., Zhao, M., Xu, T., and Klionsky, D. J. (2011). Two MAPK-signaling pathways are required for mitophagy in Saccharomyces cerevisiae. *J. Cell Biol.* 193, 755-767. doi: 10.1083/jcb.201102092

Merico, D., Isserlin, R., Stueker, O., Emili, A., and Bader, G. D. (2010). Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. *PLoS One* 5:e13984. doi: 10.1371/journal.pone.0013984

Michel, S., Keller, M. A., Wamelink, M. M. C., and Ralser, M. (2015). A haploproficient interaction of the transaldolase paralogue NQM1 with the transcription factor VHR1 affects stationary phase survival and oxidative stress resistance. *BMC Genet.* 16, 13-14. doi: 10.1186/s12863-015-0171-6

Miyaharan, K., Mizunuma, M., Hirata, D., Tsuchiya, E., and Miyakawa, T. (1996). The involvement of the Saccharomyces cerevisiae multidrug resistance transporters Pdr5p and Snq2p in cation resistance. *FEBS Lett.* 399, 317-320. doi: 10.1016/S0014-5793(96)01353-1

Mootha, V. K., Lindgren, C. M., Eriksson, K.-F., Subramanian, A., Sihag, S., Lehar, J., et al. (2003). PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. *Nat. Genet.* 34, 267-273. doi: 10.1038/ng1180

Mora-Montes, H. M., Bates, S., Netea, M. G., Diaz-Jiménez, D. F., López-Romero, E., Zinker, S., et al. (2007). Endoplasmic reticulum α-glycosidases of Candida albicans are required for N glycosylation, cell wall integrity, and normal host-fungus interaction. *Eukaryot. Cell* 6, 2184-2193. doi: 10.1128/EC.00350-07

Mora-Montes, H. M., McKenzie, C., Bain, J. M., Lewis, L. E., Erwig, L. P., and Gow, N. A. R. (2012). Interactions between macrophages and cell wall oligosaccharides of Candida albicans. *Host-Fungus Interactions* (Springer), 247-260. doi: 10.1007/978-1-61779-539-8_16

Morel, E., Cui, L., Zerges, W., and Wilkinson, K. J. (2021). Mixtures of rare earth elements show antagonistic interactions in Chlamydomonas reinhardtii. *Environ. Pollut.* 287:117594. doi: 10.1016/j.envpol.2021.117594

Moriwaki, H., and Yamamoto, H. (2013). Interactions of microorganisms with rare earth ions and their utilization for separation and environmental technology. *Appl. Microbiol. Biotechnol.* 97, 1-8. doi: 10.1007/s00253-012-4519-9

Nakagawa, T., Mitsui, R., Tani, A., Sasa, K., Tashiro, S., Iwama, T., et al. (2012). A catalytic role of XoxF1 as La3+-dependent methanol dehydrogenase in Methylobacterium extorquens strain AM1. *PLoS One* 7:e50480. doi: 10.1371/journal.pone.0050480

Navarro-Arias, M. J., Defosse, T. A., Dementhon, K., Csonka, K., Mellado-Mojica, E., Dias Valério, A., et al. (2016). Disruption of protein mannosylation affects Candida guilliermondii cell wall, immune sensing, and virulence. *Front. Microbiol.* 7:1951. doi: 10.3389/fmicb.2016.01951

Navarro-Aviño, J. P., Prasad, R., Miralles, V. J., Benito, R. M., and Serrano, R. (1999). A proposal for nomenclature of aldehyde dehydrogenases in Saccharomyces cerevisiae and characterization of the stress-inducible ALD2 and ALD3 genes. *Yeast* 15, 829-842. doi: 10.1016/j.addbeh.2008.04.015

Ngwenya, B. T., Mosselmans, J. F. W., Magennis, M., Atkinson, K. D., Tourney, J., Olive, V., et al. (2009). Macroscopic and spectroscopic analysis of lanthanide adsorption to bacterial cells. *Geochim. Cosmochim. Acta* 73, 3134-3147. doi: 10.1016/j.gca.2009.03.018

Noack, C. W., Dzombak, D. A., and Karamalidis, A. K. (2014). Rare earth element distributions and trends in natural waters with a focus on groundwater. *Environ. Sci. Technol.* 48, 4317-4326. doi: 10.1021/es4053895

Okada, H., Abe, M., Asakawa-Minemura, M., Hirata, A., Qadota, H., Morishita, K., et al. (2010). Multiple functional domains of the yeast l3-beta-glucan synthase subunit Fks1p revealed by quantitative phenotypic analysis of temperature-sensitive mutants. *Genetics* 184, 1013-1024. doi: 10.1534/genetics.109.109892

Pagani, M. A., Casamayor, A., Serrano, R., Atrian, S., and Ariño, J. (2007). Disruption of iron homeostasis in Saccharomyces cerevisiae by high zinc levels: a genome-wide study. *Mol. Microbiol.* 65, 521-537. doi: 10.1111/j.1365-2958.2007.05807.x

Pagano, G., Guida, M., Tommasi, F., and Oral, R. (2015). Health effects and toxicity mechanisms of rare earth elements—knowledge gaps and research prospects. *Ecotoxicol. Environ. Saf.* 115, 40-48. doi: 10.1016/j.ecoenv.2015.01.030

Pallares, R. M., Faulkner, D., An, D. D., Hébert, S., Loguinov, A., Proctor, M., et al. (2021). Genome-wide toxicogenomic study of the lanthanides sheds
light on the selective toxicity mechanisms associated with critical materials. *Proc. Natl. Acad. Sci. U. S. A.* 118:e2025952118. doi: 10.1073/pnas.2025952118

Pol, A., Barends, T. R. M., Dietl, A., Khadem, A. F., Eygensteyn, J., Jetten, M. S. M., et al. (2014). Rare earth metals are essential for methanotrophic life in volcanic mudpots. *Environ. Microbiol.* 16, 255–264. doi: 10.1111/1462-2920.12249

Ramalho, J., Semelka, R. C., Ramalho, M., Nunes, R. H., AlObaidy, M., and Castillo, M. (2016). Gadolinium-based contrast agent accumulation and toxicity: an update. *Am. J. Neuroradiol.* 37, 1192–1198. doi: 10.3174/ajnr.A4615

Ruotolo, R., Marchini, G., and Ottonello, S. (2008). Membrane transporters and protein traffic networks differentially affecting metal tolerance: a genomic phenotyping study in yeast. *Genome Biol.* 9:R67. doi: 10.1186/gb-2008-9-4-r67

Scherens, B., and Goffeau, A. (2004). The uses of genome-wide yeast mutant collections. *Genome Biol.* 5:229. doi: 10.1186/gb-2004-5-7-229

Serero, A., Lopes, J., Nicolas, A., and Boiteux, S. (2008). Yeast genes involved in cadmium tolerance : identification of DNA replication as a target of cadmium toxicity. *DNA Repair (Amst)* 7, 1262–1275. doi: 10.1016/j.dnarep.2008.04.005

Serrano, R., Bernal, D., Simón, E., and Ariño, J. (2004). Copper and iron are the limiting factors for growth of the yeast Saccharomyces cerevisiae in an alkaline environment. *J. Biol. Chem.* 279, 19698–19704. doi: 10.1074/jbc.M313746200

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

Skovran, E., and Martinez-Gomez, N. C. (2015). Just add lanthanides. *Science* 348, 862–863. doi: 10.1126/science.aaa9091

So, R. W. L., Chung, S. W., Lau, H. H. C., Watts, J. J., Gaudette, E., Al-Azzawi, Z. A. M., et al. (2019). Application of CRISPR genetic screens to investigate neurological diseases. *Mol. Neurodegener.* 14:41. doi: 10.1186/s13024-019-0343-3

Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. *Proc. Natl. Acad. Sci.* 102, 15545–15550. doi: 10.1073/pnas.0506580102

Sun, G. L., Reynolds, E. E., and Belcher, A. M. (2019). Designing yeast as plant-like hyperaccumulators for heavy metals. *Nat. Commun.* 10, 1–12. doi: 10.1038/s41467-019-13093-6

Thorpe, G. W., Fong, C. S., Alic, N., Higgins, V. J., and Dawes, I. W. (2004). Cells have distinct mechanisms to maintain protection against different reactive oxygen species: oxidative-stress-response genes. *Proc. Natl. Acad. Sci. U. S. A.* 101, 6564–6569. doi: 10.1073/pnas.0305888101

Thorsen, M., Di, Y., Tängemo, C., Morillas, M., Ahmadpour, D., Van der Does, C., et al. (2006). The MAPK Hog1p modulates Fps1p-dependent arsenite uptake and tolerance in yeast. *Mol. Biol. Cell* 17, 4400–4410. doi: 10.1091/mbc.e06-04-0315

Thorsen, M., Perrone, G. G., Kristiansson, E., Traini, M., Ye, T., Dawes, I. W., et al. (2009). Genetic basis of arsenite and cadmium tolerance in Saccharomyces cerevisiae. *BMC Genomics* 10, 1–15. doi: 10.1186/1471-2164-10-105

Tun, N. M., O'Doherty, P. J., Chen, Z., Wu, X., Bailey, T. D., Kersaitis, C., et al. (2014). Identification of aluminium transport-related genes via genome-wide phenotypic screening of *Saccharomyces cerevisiae* *Metallomics* 6, 1558–1564. doi: 10.1039/c4mt00116h

Tyler, G. (2004). Rare earth elements in soil and plant systems - a review. *Plant Soil* 267, 191–206. doi: 10.1007/s11104-005-4888-2

Wakabayashi, T., Ymamoto, A., Kazaana, A., Nakano, Y., Nojiri, Y., and Kashiwazaki, M. (2016). Antibacterial, antifungal and nematicidal activities of rare earth ions. *Biol. Trace Elem. Res.* 174, 464–470. doi: 10.1007/s12011-016-0727-y

Wang, J., Wang, X., Fang, Y., and Zhou, B. (2007). Genome-wide screening of yeast metal homeostasis genes involved in mitochondrial functions. *Mol. Gen. Genomics.* 277, 673–683. doi: 10.1007/s00438-007-0217-1

Wehrmann, M., Billard, P., Martin-Meriadec, A., Zegeye, A., and Klebensberger, J. (2017). Functional role of lanthanides in enzymatic activity and transcriptional regulation of pyrroloquinoline quinone-dependent alcohol dehydrogenases in Pseudomonas putida KT2440. *MBio* 8, 1–14. doi: 10.1128/mBio.00570-17

Wei, B., Li, Y., Li, H., Yu, J., Ye, B., and Liang, T. (2013). Rare earth elements in human hair from a mining area of China. *Ecotoxicol. Environ. Saf.* 96, 118–123. doi: 10.1016/j.ecoenv.2013.05.031

Weiwei, H., Shihua, S., and Peidong, T. (2007). Analysis of inhibitory effect of gadolinium on Sinorhizobium fredii. *J. Radioanal. Nucl. Chem.* 25, 106–110. doi: 10.1016/S1002-0721(07)60054-4

Yufeng, Z., Lifen, Y., Kaoshan, C., and Liang, D. (2007). Effects of neodymium on growth, pectinase activity and mycelium permeability of Fusarium oxysporum. *J. Rare Earths* 25, 100–105. doi: 10.1016/S1002-0721(07) 60053-2

**Conflict of Interest:** 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.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Grosjean, Le Jean, Chalot, Mora-Montes, Armengaud, Gross and Blaudez. 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(s) 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.