1
|
Gorrochategui E, Le Vee M, Selmi H, Gérard A, Chaker J, Krais AM, Lindh C, Fardel O, Chevrier C, Le Cann P, Miller GW, Barouki R, Jégou B, Gicquel T, Kristensen DM, David A. High-resolution mass spectrometry identifies delayed biomarkers for improved precision in acetaminophen/paracetamol human biomonitoring. Environ Int 2023; 181:108299. [PMID: 37951015 DOI: 10.1016/j.envint.2023.108299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Revised: 10/30/2023] [Accepted: 10/30/2023] [Indexed: 11/13/2023]
Abstract
Paracetamol/acetaminophen (N-acetyl-p-aminophenol, APAP) is a top selling analgesic used in more than 600 prescription and non-prescription pharmaceuticals. To study efficiently some of the potential undesirable effects associated with increasing APAP consumption (e.g., developmental disorders, drug-induced liver injury), there is a need to improve current APAP biomonitoring methods that are limited by APAP short half-life. Here, we demonstrate using high-resolution mass spectrometry (HRMS) in several human studies that APAP thiomethyl metabolite conjugates (S-methyl-3-thioacetaminophen sulfate and S-methyl-3-thioacetaminophen sulphoxide sulfate) are stable biomarkers with delayed excretion rates compared to conventional APAP metabolites, that could provide a more reliable history of APAP ingestion in epidemiological studies. We also show that these biomarkers could serve as relevant clinical markers to diagnose APAP acute intoxication in overdosed patients, when free APAP have nearly disappeared from blood. Using in vitro liver models (HepaRG cells and primary human hepatocytes), we then confirm that these thiomethyl metabolites are directly linked to the toxic N-acetyl-p-benzoquinone imine (NAPQI) elimination, and produced via an overlooked pathway called the thiomethyl shunt pathway. Further studies will be needed to determine whether the production of the reactive hepatotoxic NAPQI metabolites is currently underestimated in human. Nevertheless, these biomarkers could already serve to improve APAP human biomonitoring, and investigate, for instance, inter-individual variability in NAPQI production to study underlying causes involved in APAP-induced hepatotoxicity. Overall, our findings demonstrate the potential of exposomics-based HRMS approach to advance towards a better precision for human biomonitoring.
Collapse
Affiliation(s)
- Eva Gorrochategui
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Marc Le Vee
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Habiba Selmi
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Anne Gérard
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Jade Chaker
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Annette M Krais
- Division of Occupational and Environmental Medicine, Lund University, Lund, Sweden
| | - Christian Lindh
- Division of Occupational and Environmental Medicine, Lund University, Lund, Sweden
| | - Olivier Fardel
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Cécile Chevrier
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Pierre Le Cann
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Gary W Miller
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY, USA
| | - Robert Barouki
- Unité UMR-S 1124 Inserm-Université Paris Descartes "Toxicologie Pharmacologie et Signalisation Cellulaire", Paris, France
| | - Bernard Jégou
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France
| | - Thomas Gicquel
- NuMeCan Institute (Nutrition, Metabolisms and Cancer), CHU Rennes, Univ Rennes, INSERM, INRAE, UMR_A 1341, UMR_S 1317, 35000 Rennes, France
| | - David M Kristensen
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France; Department of Science and Environment, Roskilde University, Roskilde, Denmark; Department of Growth and Reproduction, Copenhagen University Hospital - Rigshospitalet, University of Copenhagen, Denmark
| | - Arthur David
- Univ Rennes, Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR_S 1085, Rennes, France.
| |
Collapse
|
2
|
Gorrochategui E, Hernandez I, Pérez-Gabucio E, Lacorte S, Tauler R. Temporal air quality (NO 2, O 3, and PM 10) changes in urban and rural stations in Catalonia during COVID-19 lockdown: an association with human mobility and satellite data. Environ Sci Pollut Res Int 2022; 29:18905-18922. [PMID: 34705210 PMCID: PMC8549430 DOI: 10.1007/s11356-021-17137-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Accepted: 10/17/2021] [Indexed: 05/09/2023]
Abstract
In this study, changes in air quality by NO2, O3, and PM10 in Barcelona metropolitan area and other parts of Catalonia during the COVID-19 lockdown with respect to pre-lockdown and to previous years (2018 and 2019) were evaluated. Selected air monitoring stations included 3 urban (Gràcia, Vall d'Hebron, and Granollers), 1 control site (Fabra Observatory), 1 semi-urban (Manlleu), and 3 rural (Begur, Bellver de Cerdanya, and Juneda). NO2 lockdown levels showed a diminution, which in relative terms was maximum in two rural stations (Bellver de Cerdanya, - 63% and Begur, - 61%), presumably due to lower emissions from the ceasing hotel and ski resort activities during eastern holidays. In absolute terms and from an epidemiologic perspective, decrease in NO2, also reinforced by the high amount of rainfall registered in April 2020, was more relevant in the urban stations around Barcelona. O3 levels increased in the transited urban stations (Gràcia, + 42%, and Granollers, + 64%) due to the lower titration effect by NOx. PM10 lockdown levels decreased, mostly in Gràcia, Vall d'Hebron, and Granollers (- 35, - 39%, and - 39%, respectively) due to traffic depletion (- 90% in Barcelona's transport). Correlation among mobility index in Barcelona (- 100% in retail and recreation) and contamination was positive for NO2 and PM10 and negative for O3 (P < 0.001). Satellite images evidenced two hotspots of NO2 in Spain (Madrid and Barcelona) in April 2018 and 2019 that disappeared in 2020. Overall, the benefits of lockdown on air quality in Catalonia were evidenced with NO2, O3 and PM10 levels below WHOAQG values in most of stations opposed to the excess registered in previous years.
Collapse
Affiliation(s)
- Eva Gorrochategui
- Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), 08034, Barcelona, Spain.
| | - Isabel Hernandez
- Direcció General de Qualitat Ambiental I Canvi Climàtic, Generalitat de Catalunya, Barcelona, Spain
| | - Eva Pérez-Gabucio
- Direcció General de Qualitat Ambiental I Canvi Climàtic, Generalitat de Catalunya, Barcelona, Spain
| | - Sílvia Lacorte
- Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), 08034, Barcelona, Spain
| | - Romà Tauler
- Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), 08034, Barcelona, Spain
| |
Collapse
|
3
|
Abstract
BACKGROUND The analysis of LC-MS metabolomic datasets appears to be a challenging task in a wide range of disciplines since it demands the highly extensive processing of a vast amount of data. Different LC-MS data analysis packages have been developed in the last few years to facilitate this analysis. However, most of these strategies involve chromatographic alignment and peak shaping and often associate each "feature" (i.e., chromatographic peak) with a unique m/z measurement. Thus, the development of an alternative data analysis strategy that is applicable to most types of MS datasets and properly addresses these issues is still a challenge in the metabolomics field. RESULTS Here, we present an alternative approach called ROIMCR to: i) filter and compress massive LC-MS datasets while transforming their original structure into a data matrix of features without losing relevant information through the search of regions of interest (ROIs) in the m/z domain and ii) resolve compressed data to identify their contributing pure components without previous alignment or peak shaping by applying a Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) analysis. In this study, the basics of the ROIMCR method are presented in detail and a detailed description of its implementation is also provided. Data were analyzed using the MATLAB (The MathWorks, Inc., www.mathworks.com ) programming and computing environment. The application of the ROIMCR methodology is described in detail, with an example of LC-MS data generated in a lipidomic study and with other examples of recent applications. CONCLUSIONS The methodology presented here combines the benefits of data filtering and compression based on the searching of ROI features, without the loss of spectral accuracy. The method has the benefits of the application of the powerful MCR-ALS data resolution method without the necessity of performing chromatographic peak alignment or modelling. The presented method is a powerful alternative to other existing data analysis approaches that do not use the MCR-ALS method to resolve LC-MS data. The ROIMCR method also represents an improved strategy compared to the direct applications of the MCR-ALS method that use less-powerful data compression strategies such as binning and windowing. Overall, the strategy presented here confirms the usefulness of the ROIMCR chemometrics method for analyzing LC-MS untargeted metabolomics data.
Collapse
Affiliation(s)
- Eva Gorrochategui
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Jorsi Girona 18-25, Barcelona, 08034, Catalonia, Spain
| | - Joaquim Jaumot
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Jorsi Girona 18-25, Barcelona, 08034, Catalonia, Spain
| | - Romà Tauler
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Jorsi Girona 18-25, Barcelona, 08034, Catalonia, Spain.
| |
Collapse
|
4
|
Gorrochategui E, Li J, Fullwood NJ, Ying GG, Tian M, Cui L, Shen H, Lacorte S, Tauler R, Martin FL. Diet-sourced carbon-based nanoparticles induce lipid alterations in tissues of zebrafish (Danio rerio) with genomic hypermethylation changes in brain. Mutagenesis 2016; 32:91-103. [PMID: 27798195 PMCID: PMC5180172 DOI: 10.1093/mutage/gew050] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
With rising environmental levels of carbon-based nanoparticles (CBNs), there is an urgent need to develop an understanding of their biological effects in order to generate appropriate risk assessment strategies. Herein, we exposed zebrafish via their diet to one of four different CBNs: C60 fullerene (C60), single-walled carbon nanotubes (SWCNT), short multi-walled carbon nanotubes (MWCNTs) or long MWCNTs. Lipid alterations in male and female zebrafish were explored post-exposure in three target tissues (brain, gonads and gastrointestinal tract) using ‘omic’ procedures based in liquid chromatography coupled with mass spectrometry (LC-MS) data files. These tissues were chosen as they are often target tissues following environmental exposure. Marked alterations in lipid species are noted in all three tissues. To further explore CBN-induced brain alterations, Raman microspectroscopy analysis of lipid extracts was conducted. Marked lipid alterations are observed with males responding differently to females; in addition, there also appears to be consistent elevations in global genomic methylation. This latter observation is most profound in female zebrafish brain tissues post-exposure to short MWCNTs or SWCNTs (P < 0.05). This study demonstrates that even at low levels, CBNs are capable of inducing significant cellular and genomic modifications in a range of tissues. Such alterations could result in modified susceptibility to other influences such as environmental exposures, pathology and, in the case of brain, developmental alterations.
Collapse
Affiliation(s)
- Eva Gorrochategui
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Barcelona, 08034, Catalonia, Spain
| | | | - Nigel J Fullwood
- Biomedical and Life Sciences Division, Lancaster Environment Centre, Lancaster University, Lancaster LA1 4YQ, UK
| | - Guang-Guo Ying
- State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China
| | - Meiping Tian
- Key Lab of Urban Environment and Health, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China and
| | - Li Cui
- Key Lab of Urban Environment and Health, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China and
| | - Heqing Shen
- Key Lab of Urban Environment and Health, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China and
| | - Sílvia Lacorte
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Barcelona, 08034, Catalonia, Spain
| | - Romà Tauler
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Barcelona, 08034, Catalonia, Spain
| | - Francis L Martin
- Centre for Biophotonics and .,Biosciences, School of Pharmacy and Biomedical Sciences, Maudland Building, University of Central Lancashire, Preston PR1 2HE, UK
| |
Collapse
|
5
|
Gorrochategui E, Jaumot J, Lacorte S, Tauler R. Data analysis strategies for targeted and untargeted LC-MS metabolomic studies: Overview and workflow. Trends Analyt Chem 2016. [DOI: 10.1016/j.trac.2016.07.004] [Citation(s) in RCA: 187] [Impact Index Per Article: 23.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
6
|
Gorrochategui E, Lacorte S, Tauler R, Martin FL. Perfluoroalkylated Substance Effects in Xenopus laevis A6 Kidney Epithelial Cells Determined by ATR-FTIR Spectroscopy and Chemometric Analysis. Chem Res Toxicol 2016; 29:924-32. [PMID: 27078751 PMCID: PMC4870675 DOI: 10.1021/acs.chemrestox.6b00076] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
![]()
The effects of four perfluoroalkylated
substances (PFASs), namely,
perfluorobutanesulfonate (PFBS), perfluorooctanoic acid (PFOA), perfluorooctanesulfonate
(PFOS), and perfluorononanoic acid (PFNA) were assessed in Xenopus laevis A6 kidney epithelial cells by attenuated
total reflection Fourier-transform infrared (ATR-FTIR) spectroscopy
and chemometric analysis. Principal component analysis–linear
discriminant analysis (PCA-LDA) was used to visualize wavenumber-related
alterations and ANOVA-simultaneous component analysis (ASCA) allowed
data processing considering the underlying experimental design. Both
analyses evidenced a higher impact of low-dose PFAS-treatments (10–9 M) on A6 cells forming monolayers, while there was
a larger influence of high-dose PFAS-treatments (10–5 M) on A6 cells differentiated into dome structures. The observed
dose–response PFAS-induced effects were to some extent related
to their cytotoxicity: the EC50-values of most influential
PFAS-treatments increased (PFOS < PFNA < PFOA ≪ PFBS),
and higher-doses of these chemicals induced a larger impact. Major
spectral alterations were mainly attributed to DNA/RNA, secondary
protein structure, lipids, and fatty acids. Finally, PFOS and PFOA
caused a decrease in A6 cell numbers compared to controls, whereas
PFBS and PFNA did not significantly change cell population levels.
Overall, this work highlights the ability of PFASs to alter A6 cells,
whether forming monolayers or differentiated into dome structures,
and the potential of PFOS and PFOA to induce cell death.
Collapse
Affiliation(s)
- Eva Gorrochategui
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC) , Barcelona 08034, Catalonia, Spain
| | - Sílvia Lacorte
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC) , Barcelona 08034, Catalonia, Spain
| | - Romà Tauler
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC) , Barcelona 08034, Catalonia, Spain
| | - Francis L Martin
- Centre for Biophotonics, Lancaster Environment Centre, Lancaster University , Lancaster LA1 4YQ, U.K.,School of Pharmacy and Biomedical Sciences, University of Central Lancashire , Preston, U.K
| |
Collapse
|
7
|
Gorrochategui E, Casas J, Porte C, Lacorte S, Tauler R. Chemometric strategy for untargeted lipidomics: Biomarker detection and identification in stressed human placental cells. Anal Chim Acta 2015; 854:20-33. [DOI: 10.1016/j.aca.2014.11.010] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2014] [Revised: 10/20/2014] [Accepted: 11/06/2014] [Indexed: 10/24/2022]
|
8
|
Gorrochategui E, Casas J, Pérez-Albaladejo E, Jáuregui O, Porte C, Lacorte S. Characterization of complex lipid mixtures in contaminant exposed JEG-3 cells using liquid chromatography and high-resolution mass spectrometry. Environ Sci Pollut Res Int 2014; 21:11907-11916. [PMID: 24969426 DOI: 10.1007/s11356-014-3172-5] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2013] [Accepted: 06/06/2014] [Indexed: 06/03/2023]
Abstract
The aim of this study was to develop a method based on ultra high performance liquid chromatography coupled with mass spectrometry (UHPLC-MS) for lipid profiling in human placental choriocarcinoma (JEG-3) cells. Lipids were solid-liquid extracted from JEG-3 cells using a solution of chloroform/methanol (2:1, v/v) in a simple procedure requiring minimal sample alteration. Simultaneous separation of complex lipid mixtures in their major classes was achieved with a reversed-phase (C8) UHPLC column and a mobile phase containing methanol with 1 mM ammonium formate and 0.2 % formic acid (A)/water with 2 mM ammonium formate and 0.2 % formic acid (B). Lipids were characterized using time-of-flight (TOF) and Orbitrap under full scan and positive electrospray ionization mode with both analyzers. A total of 178 species of lipids, including 37 phosphatidylcholines (PC), 32 plasmalogen PC, 9 lyso PC, 4 lyso plasmalogen PC, 30 triacylglycerols, 22 diacylglycerols, 7 cholesterol esters, 25 phosphatidylethanolamines, and 12 sphingomyelins, were identified using TOF and Orbitrap. The identification of all lipid classes was based on exact mass characterization with an error < 5 ppm. The developed methodology was applied to study lipid alterations in human placental cells against the exposure to perfluorinated chemicals (PFCs) and tributyltin (TBT).
Collapse
Affiliation(s)
- Eva Gorrochategui
- Department of Environmental Chemistry, Institute of Environmental Assessment and Water Research (IDAEA), Consejo Superior de Investigaciones Científicas (CSIC), Jordi Girona 18-26, Barcelona, 08034, Catalonia, Spain
| | | | | | | | | | | |
Collapse
|
9
|
Gorrochategui E, Pérez-Albaladejo E, Casas J, Lacorte S, Porte C. Perfluorinated chemicals: differential toxicity, inhibition of aromatase activity and alteration of cellular lipids in human placental cells. Toxicol Appl Pharmacol 2014; 277:124-30. [PMID: 24680846 DOI: 10.1016/j.taap.2014.03.012] [Citation(s) in RCA: 93] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2014] [Revised: 03/12/2014] [Accepted: 03/16/2014] [Indexed: 01/18/2023]
Abstract
The cytotoxicity of eight perfluorinated chemicals (PFCs), namely, perfluorobutanoic acid (PFBA), perfluorohexanoic acid (PFHxA), perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), perfluorododecanoic acid (PFDoA), perfluorobutanesulfonate (PFBS), perfluorohexanesulfonate (PFHxS) and perfluorooctanesulfonate (PFOS) was assessed in the human placental choriocarcinoma cell line JEG-3. Only the long chain PFCs--PFOS, PFDoA, PFNA, PFOA--showed significant cytotoxicity in JEG-3 cells with EC50 values in the range of 107 to 647 μM. The observed cytotoxicity was to some extent related to a higher uptake of the longer chain PFCs by cells (PFDoA>PFOS≫PFNA>PFOA>PFHxA). Moreover, this work evidences a high potential of PFOS, PFOA and PFBS to act as aromatase inhibitors in placental cells with IC50s in the range of 57-80 μM, the inhibitory effect of PFBS being particularly important despite the rather low uptake of the compound by cells. Finally, exposure of JEG-3 cells to a mixture of the eight PFCs (0.6 μM each) led to a relative increase (up to 3.4-fold) of several lipid classes, including phosphatidylcholines (PCs), plasmalogen PC and lyso plasmalogen PC, which suggests an interference of PFCs with membrane lipids. Overall, this work highlights the ability of the PFC mixture to alter cellular lipid pattern at concentrations well below those that generate toxicity, and the potential of the short chain PFBS, often considered a safe substitute of PFOS, to significantly inhibit aromatase activity in placental cells.
Collapse
Affiliation(s)
- Eva Gorrochategui
- Department of Environmental Chemistry, IDAEA-CSIC, 08034 Barcelona, Catalonia, Spain
| | | | - Josefina Casas
- Department of Biomedicinal Chemistry, IQAC-CSIC, 08034 Barcelona, Catalonia, Spain
| | - Sílvia Lacorte
- Department of Environmental Chemistry, IDAEA-CSIC, 08034 Barcelona, Catalonia, Spain.
| | - Cinta Porte
- Department of Environmental Chemistry, IDAEA-CSIC, 08034 Barcelona, Catalonia, Spain.
| |
Collapse
|