1
|
Ignatz EH, Rise ML, Gamperl AK. Impact of stress phenotype, elevated temperature, and bacterin exposure on male Atlantic salmon ( Salmo salar) growth, stress, and immune biomarker gene expression. Physiol Genomics 2023; 55:587-605. [PMID: 37746713 DOI: 10.1152/physiolgenomics.00055.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Revised: 08/30/2023] [Accepted: 09/21/2023] [Indexed: 09/26/2023] Open
Abstract
In this study, postsmolt male Atlantic salmon, previously identified as low responders (LRs) or high responders (HRs) based on poststress cortisol levels, had their head kidney and liver sampled at 12°C and 20°C before injection (time 0) and after injection (i.e., at 12- and 24-h postinjection, respectively) with either Forte Micro (a multivalent vaccine containing bacterin, to capture peak antibacterial responses) or an equal volume of PBS. Quantitative real-time PCR (qPCR) was then used to measure the expression of 15 biomarker genes in the head kidney and 12 genes in the liver at each temperature/sampling point. Target transcripts were chosen that were related to growth, stress, and innate antibacterial immune responses. Many temperature, phenotype, and injection effects were found for individual genes within these three broad categories, and multivariate statistical analyses (i.e., principal component analysis and permutational multivariate analysis of variance) were used to look for overall patterns in transcript expression. These analyses revealed that HR salmon at 20°C mounted a more robust response (P < 0.05) for the 10 head kidney immune-related transcripts when injected with Forte Micro than LR salmon. In contrast, the seven liver stress-related transcripts displayed a greater response (P = 0.057) in LR versus HR fish with Forte Micro at 12°C. Overall, although this research did find some differences between LR and HR fish, it does not provide strong (conclusive) evidence that the selection of a particular phenotype would have major implications for the health of salmon over the temperature range examined.NEW & NOTEWORTHY This is the first paper to describe the impact of both temperature and bacterial stimulation on head kidney and liver transcript expression in Atlantic salmon characterized as LRs versus HRs. Notably, we found that HR salmon at 20°C mounted a more robust innate antibacterial immune response than LR salmon. In addition, LR fish at 12°C may (P = 0.057) exhibit higher expression of stress-related transcripts in response to vaccine injection relative to HR fish.
Collapse
Affiliation(s)
- Eric H Ignatz
- Department of Ocean Sciences, Memorial University, St. John's, Newfoundland and Labrador, Canada
| | - Matthew L Rise
- Department of Ocean Sciences, Memorial University, St. John's, Newfoundland and Labrador, Canada
| | - A Kurt Gamperl
- Department of Ocean Sciences, Memorial University, St. John's, Newfoundland and Labrador, Canada
| |
Collapse
|
2
|
Jan K, Ahmed I, Dar NA, Farah MA, Khan FR, Shah BA, Fazio F. LC-MS/MS based characterisation and differential expression of proteins in Himalayan snow trout, Schizothorax labiatus using LFQ technique. Sci Rep 2023; 13:10134. [PMID: 37349327 PMCID: PMC10287682 DOI: 10.1038/s41598-023-35646-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2023] [Accepted: 05/22/2023] [Indexed: 06/24/2023] Open
Abstract
Molecular characterization of fish muscle proteins are nowadays considered as a key component to understand the role of specific proteins involved in various physiological and metabolic processes including their up and down regulation in the organisms. Coldwater fish specimens including snow trouts hold different types of proteins which help them to survive in highly diversified temperatures fluctuating from 0 to 20 °C. So, in current study, the liquid chromatography mass spectrometry using label free quantification technique has been used to investigate the muscle proteome profile of Schizothorax labiatus. For proteomic study, two weight groups of S. labiatus were taken from river Sindh. The proteomic analysis of group 1 revealed that a total of 235 proteins in male and 238 in female fish were recorded. However, when male and female S. labiatus were compared with each other on the basis of spectral count and abundance of peptides by ProteinLynx Global Server software, a total of 14 down-regulated and 22 up-regulated proteins were noted in this group. The highly down-regulated ones included homeodomain protein HoxA2b, retinol-binding protein 4, MHC class II beta chain and proopiomelanocortin while as the highly expressed up-regulated proteins comprised of gonadotropin I beta subunit, NADH dehydrogenase subunit 4, manganese superoxide dismutase, recombinase-activating protein 2, glycosyltransferase, chymotrypsin and cytochrome b. On the other hand, the proteomic characterisation of group 2 of S. labiatus revealed that a total of 227 proteins in male and 194 in female fish were recorded. When male and female S. labiatus were compared with each other by label free quantification, a total of 20 down-regulated and 18 up-regulated proteins were recorded. The down-regulated protein expression of group 2 comprised hepatic lipase, allograft inflammatory factor-1, NADH dehydrogenase subunit 4 and myostatin 1 while the highly expressed up-regulated proteins included glycogen synthase kinase-3 beta variant 2, glycogen synthase kinase-3 beta variant 5, cholecystokinin, glycogen synthase kinase-3 beta variant 3 and cytochrome b. Significant (P < 0.05) difference in the expression of down-regulated and up-regulated proteins was also noted between the two sexes of S. labiatus in each group. According to MS analysis, the proteins primarily concerned with the growth, skeletal muscle development and metabolism were down-regulated in river Sindh, which indicates that growth of fish during the season of collection i.e., winter was slow owing to less food availability, gonad development and low metabolic activity. While, the proteins related to immune response of fish were also noted to be down-regulated thereby signifying that the ecosystem has less pollution loads, microbial, pathogenic and anthropogenic activities. It was also found that the proteins involved in glycogen metabolism, reproductive and metabolic processes, particularly lipid metabolism were up-regulated in S. labiatus. The significant expression of these proteins may be connected to pre-spawning, gonad development and use of stored food as source of energy. The information generated in this study can be applied to future research aimed at enhancing food traceability, food safety, risk management and authenticity analysis.
Collapse
Affiliation(s)
- Kousar Jan
- Fish Nutrition Research Laboratory, Department of Zoology, University of Kashmir, Hazratbal, Srinagar, Jammu and Kashmir, 190 006, India
| | - Imtiaz Ahmed
- Fish Nutrition Research Laboratory, Department of Zoology, University of Kashmir, Hazratbal, Srinagar, Jammu and Kashmir, 190 006, India.
| | - Nazir Ahmad Dar
- Department of Biochemistry, University of Kashmir, Hazratbal, Srinagar, 190006, India
| | - Mohammad Abul Farah
- Department of Zoology, College of Science, King Saud University, Riyadh, 11451, Saudi Arabia
| | - Fatin Raza Khan
- Department of Biotechnology and Bioinformatics, University of Hyderabad, Hyderabad, India
| | - Basit Amin Shah
- Department of Biotechnology, University of Kashmir, Hazratbal, Srinagar, 190006, India
| | - Francesco Fazio
- Department of Veterinary Sciences, Polo Universitario Annunziata, University of Messina, 98168, Messina, Italy
| |
Collapse
|
3
|
Sun B, van Dissel D, Mo I, Boysen P, Haslene-Hox H, Lund H. Identification of novel biomarkers of inflammation in Atlantic salmon (Salmo salar L.) by a plasma proteomic approach. DEVELOPMENTAL AND COMPARATIVE IMMUNOLOGY 2022; 127:104268. [PMID: 34571096 DOI: 10.1016/j.dci.2021.104268] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2021] [Revised: 09/22/2021] [Accepted: 09/22/2021] [Indexed: 06/13/2023]
Abstract
Monitoring fish welfare has become a central issue for the fast-growing aquaculture industry, and finding proper biomarkers of stress, inflammation and infection is necessary for surveillance and documentation of fish health. In this study, a proteomic approach using mass spectrometry was applied to identify indicators of the acute response in Atlantic salmon blood plasma by comparing Aeromonas salmonicida subsp. salmonicida infected fish and non-infected controls. The antimicrobial proteins cathelicidin (CATH), L-plastin (Plastin-2, LCP1) and soluble toll-like receptor 5 (sTLR5) were uniquely or mainly identified in the plasma of infected fish. In addition, five immune-related proteins showed significantly increased expression in plasma of infected fish: haptoglobin, high affinity immunoglobulin Fc gamma receptor I (FcγR1, CD64), leucine-rich alpha 2 glycoprotein (LRG1), complement C4 (C4) and phospholipase A2 inhibitor 31 kDa subunit-like protein. However, various fibrinogen components, CD209 and CD44 antigen-like molecules decreased in infected fish. Selected biomarkers were further verified by Western blot analysis of plasma and real time PCR of spleen and liver, including CATH1, CATH2 and L-plastin. A significant increase of L-plastin occurred as early as 24 h after infection, and a CATH2 increase was observed from 72 h in plasma of infected fish. Real time PCR of selected genes confirmed increased transcription of CATH1 and CATH2. In addition, serum amyloid A mRNA significantly increased in liver and spleen after bacterial infection. However, transcription of L-plastin was not consistently induced in liver and spleen. The results of the present study reveal novel and promising biomarkers of the acute phase response and inflammation in Atlantic salmon.
Collapse
Affiliation(s)
- Baojian Sun
- Faculty of Veterinary Medicine, Norwegian University of Life Sciences (NMBU), Ås, Norway
| | - Dino van Dissel
- SINTEF AS, Department of Biotechnology and Nanomedicine, Trondheim, Norway
| | - Ingrid Mo
- Faculty of Veterinary Medicine, Norwegian University of Life Sciences (NMBU), Ås, Norway
| | - Preben Boysen
- Faculty of Veterinary Medicine, Norwegian University of Life Sciences (NMBU), Ås, Norway
| | - Hanne Haslene-Hox
- SINTEF AS, Department of Biotechnology and Nanomedicine, Trondheim, Norway
| | - Hege Lund
- Faculty of Veterinary Medicine, Norwegian University of Life Sciences (NMBU), Ås, Norway.
| |
Collapse
|
4
|
Multi-Approach Assessment for Stress Evaluation in Rainbow Trout Females, Oncorhynchus mykiss (Walbaum, 1792) from Three Different Farms during the Summer Season. Animals (Basel) 2021; 11:ani11061810. [PMID: 34204450 PMCID: PMC8234804 DOI: 10.3390/ani11061810] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Revised: 06/11/2021] [Accepted: 06/14/2021] [Indexed: 01/02/2023] Open
Abstract
Blood biochemistry parameters are valuable tools for monitoring fish health. Their baseline values are still undefined for a multitude of farmed fish species. In this study, changes in the blood profile of rainbow trout females (Oncorhynchus mykiss) from three farms were investigated using different biomarkers during the summer season. In the given context, the main water physicochemical parameters were investigated and twelve biochemical parameters were measured from blood samples of rainbow trout reared in the Fiad, Șoimul de Jos, and Strâmba farms. We selected these farms because the genetic background of the rainbow trout is the same, with all studied specimens coming from the Fiad farm, which has an incubation station. Forty-five samples were collected monthly (May to August) throughout summer to observe the changes in the blood profile of rainbow trout. Principal component analysis showed a clear separation both among the studied farms and months. Furthermore, significant correlations (p < 0.05) between the majority of the biochemical parameters were found, indicating that the environmental parameters can influence several blood parameters at the same time. The present study provides several useful norms for assessing the welfare of rainbow trout, indicating that the relationships among different parameters are important factors in interpreting the blood biochemical profiles.
Collapse
|
5
|
Verleih M, Borchel A, Rebl A, Brenmoehl J, Kühn C, Goldammer T. A molecular survey of programmed cell death in rainbow trout: Structural and functional specifications of apoptotic key molecules. Comp Biochem Physiol B Biochem Mol Biol 2019; 230:57-69. [DOI: 10.1016/j.cbpb.2019.01.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2018] [Revised: 01/11/2019] [Accepted: 01/18/2019] [Indexed: 12/24/2022]
|
6
|
Samaras A, Pavlidis M. Regulation of divergent cortisol responsiveness in European sea bass, Dicentrarchus labrax L. PLoS One 2018; 13:e0202195. [PMID: 30096195 PMCID: PMC6086447 DOI: 10.1371/journal.pone.0202195] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2018] [Accepted: 07/30/2018] [Indexed: 11/18/2022] Open
Abstract
Mechanisms regulating differences in cortisol responsiveness between low (LR) and high response (HR) individuals have been poorly studied. In this context, we aimed to study key regulatory processes in cortisol dynamics at the head kidneys of LR and HR European sea bass. To do so, resting plasma cortisol and ACTH concentrations were quantified in these fish. Additionally, the head kidneys of these individuals were superfused through an in vitro superfusion system and stimulated with the same amount of ACTH to assess their cortisol biosynthetic capacity. Moreover, the expression of important genes in cortisol regulation was assessed. Results showed that LR fish had lower resting cortisol concentrations than HR, although no differences existed in the circulating levels of ACTH. Additionally, the biosynthetic capacity of HR was higher than that of LR fish when in vitro stimulated with ACTH. At the molecular level, a statistically significant 3.4-fold higher expression of the ACTH receptor, mc2r, and a 2.3-fold, though not significant, higher expression of 11β-hydroxylase (cyp11b1), an enzyme involved in cortisol biosynthesis, was observed in the HR fish. Finally, a statistically significant 1.3-fold lower expression of 11β-hydroxysteroid dehydrogenase 2 (hsd11b2), an enzyme involved in cortisol inactivation, was observed in HR when compared to LR fish. Therefore, it was for the first time indicated that cortisol dynamics can also be regulated at the post-production level in the head kidney. Collectively, our results highlight the crucial role of the interrenal tissue in the regulation of differences in cortisol response between LR and HR sea bass individuals.
Collapse
Affiliation(s)
| | - Michail Pavlidis
- Department of Biology, University of Crete, Heraklion, Crete, Greece
| |
Collapse
|
7
|
Samaras A, Dimitroglou A, Sarropoulou E, Papaharisis L, Kottaras L, Pavlidis M. Repeatability of cortisol stress response in the European sea bass (Dicentrarchus labrax) and transcription differences between individuals with divergent responses. Sci Rep 2016; 6:34858. [PMID: 27703277 PMCID: PMC5050510 DOI: 10.1038/srep34858] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2016] [Accepted: 09/05/2016] [Indexed: 11/09/2022] Open
Abstract
Understanding the stress responses of organisms is of importance in the performance and welfare of farmed animals, including fish. Especially fish in aquaculture commonly face stressors, and better knowledge of their responses may assist in proper husbandry and selection of breeding stocks. European sea bass (Dicentrarchus labrax), a species with high cortisol concentrations, is of major importance in this respect. The main objectives of the present study were to assess the repeatability and consistency of cortisol stress response and to identify differences in liver transcription profiles of European sea bass individuals, showing a consistent low (LR) or high (HR) cortisol response. The progeny of six full sib families was used, and sampled for plasma cortisol after an acute stress challenge once per month, for four consecutive months. Results suggest that cortisol responsiveness was a repeatable trait with LR and HR fish showing low or high resting, free and post-stress cortisol concentrations respectively. Finally, the liver transcription profiles of LR and HR fish showed some important differences, indicating differential hepatic regulation between these divergent phenotypes. These transcription differences were related to various metabolic and immunological processes, with 169 transcripts being transcribed exclusively in LR fish and 161 exclusively in HR fish.
Collapse
Affiliation(s)
- A Samaras
- Department of Biology, University of Crete, Heraklion, Crete, Greece
| | - A Dimitroglou
- Research and Development Department, Nireus Aquaculture S.A., Greece
| | - E Sarropoulou
- Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Centre for Marine Research, Heraklion, Greece
| | - L Papaharisis
- Research and Development Department, Nireus Aquaculture S.A., Greece
| | - L Kottaras
- Research and Development Department, Nireus Aquaculture S.A., Greece
| | - M Pavlidis
- Department of Biology, University of Crete, Heraklion, Crete, Greece
| |
Collapse
|
8
|
Kocmarek AL, Ferguson MM, Danzmann RG. Co-localization of growth QTL with differentially expressed candidate genes in rainbow trout. Genome 2015; 58:393-403. [PMID: 26360524 DOI: 10.1139/gen-2015-0047] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We tested whether genes differentially expressed between large and small rainbow trout co-localized with familial QTL regions for body size. Eleven chromosomes, known from previous work to house QTL for weight and length in rainbow trout, were examined for QTL in half-sibling families produced in September (1 XY male and 1 XX neomale) and December (1 XY male). In previous studies, we identified 108 candidate genes for growth expressed in the liver and white muscle in a subset of the fish used in this study. These gene sequences were BLASTN aligned against the rainbow trout and stickleback genomes to determine their location (rainbow trout) and inferred location based on synteny with the stickleback genome. Across the progeny of all three males used in the study, 63.9% of the genes with differential expression appear to co-localize with the QTL regions on 6 of the 11 chromosomes tested in these males. Genes that co-localized with QTL in the mixed-sex offspring of the two XY males primarily showed up-regulation in the muscle of large fish and were related to muscle growth, metabolism, and the stress response.
Collapse
Affiliation(s)
- Andrea L Kocmarek
- Department of Integrative Biology, 50 Stone Rd. East, University of Guelph, Guelph, ON N1G 2W1, Canada.,Department of Integrative Biology, 50 Stone Rd. East, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - Moira M Ferguson
- Department of Integrative Biology, 50 Stone Rd. East, University of Guelph, Guelph, ON N1G 2W1, Canada.,Department of Integrative Biology, 50 Stone Rd. East, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - Roy G Danzmann
- Department of Integrative Biology, 50 Stone Rd. East, University of Guelph, Guelph, ON N1G 2W1, Canada.,Department of Integrative Biology, 50 Stone Rd. East, University of Guelph, Guelph, ON N1G 2W1, Canada
| |
Collapse
|
9
|
Rexroad CE, Vallejo RL, Liu S, Palti Y, Weber GM. Quantitative trait loci affecting response to crowding stress in an F(2) generation of rainbow trout produced through phenotypic selection. MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2013; 15:613-627. [PMID: 23709047 DOI: 10.1007/s10126-013-9512-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2013] [Accepted: 05/09/2013] [Indexed: 06/02/2023]
Abstract
Selective breeding programs for salmonids typically aim to improve traits associated with growth and disease resistance. It has been established that stressors common to production environments can adversely affect these and other traits which are important to producers and consumers. Previously, we employed phenotypic selection to create families that exhibit high or low plasma cortisol concentrations in response to crowding stress. Subsequent crosses of high × low phenotypes founded a multigenerational breeding scheme with the aim of dissecting the genetic basis for variation underlying stress response through the identification of quantitative trait loci (QTL). Multiple methods of QTL analyses differing in their assumptions of homozygosity of the causal alleles in the grandparental generation yielded similar results in the F1 generation, and the analysis of two stress response phenotype measurement indexes were highly correlated. In the current study, we conducted a genome scan with microsatellites to detect QTL in the F2 generation of two families created through phenotypic selection and having larger numbers of offspring than families screened in the previous generation. Seven suggestive and three significant QTL were detected, seven of which were not previously detected in the National Center for Cool and Cold Water Aquaculture germplasm, bringing the total number of chromosomes containing significant and suggestive stress response QTL to 4 and 15, respectively. One significant QTL which peaks at 7 cM on chromosome Omy12 spans 12 cM and explains 25 % of the phenotypic variance in family 2008052 particularly warrants further investigation. Five QTL with significant parent-of-origin effects were detected in family 2008052, including two QTL on Omy12. The 95 % confidence intervals for the remaining QTL we detected were broad, requiring validation and fine mapping with other genotyping approaches and mapping strategies. These results will facilitate identification of potential casual alleles that can be employed in strategies aimed at better understanding the genetic and physiological basis of stress responses to crowding in rainbow trout aquaculture production.
Collapse
Affiliation(s)
- Caird E Rexroad
- USDA/ARS National Center for Cool and Cold Water Aquaculture, 11861 Leetown Road, Kearneysville, WV 25430, USA.
| | | | | | | | | |
Collapse
|
10
|
Calduch-Giner JA, Sitjà-Bobadilla A, Davey GC, Cairns MT, Kaushik S, Pérez-Sánchez J. Dietary vegetable oils do not alter the intestine transcriptome of gilthead sea bream (Sparus aurata), but modulate the transcriptomic response to infection with Enteromyxum leei. BMC Genomics 2012; 13:470. [PMID: 22967181 PMCID: PMC3444936 DOI: 10.1186/1471-2164-13-470] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2012] [Accepted: 09/07/2012] [Indexed: 11/30/2022] Open
Abstract
Background Studies conducted with gilthead sea bream (Sparus aurata L.) have determined the maximum dietary replacement of fish meal and oil without compromising growth or product quality. The present study aimed to analyze the effect of the nutritional background on fish health and fish fed plant protein-based diets with fish oil (FO diet) or a blend of vegetable oils (66VO diet) were exposed for 102 days to the intestinal myxosporean parasite Enteromyxum leei, and the intestine transcriptome was analyzed with a customized oligo-microarray of 7,500 annotated genes. Results Infection prevalence was high and similar in the two diet groups, but the outcome of the disease was more pronounced in fish fed the 66VO diet. No differences were found in the transcriptome of both diet control groups, whereas the number of differentially expressed genes in infected groups was considerable. K-means clustering of these differentially expressed genes identified four expression patterns that reflected the progression of the disease with the magnitude of the fold-change being higher in infected 66VO fish. A positive correlation was found between the time of infection and the magnitude of the transcriptional change within the 66VO group, being higher in early infected animals. Within this diet group, a strong up-regulation of many components of the immune specific response was evidenced, whereas other genes related to complement response and xenobiotic metabolism were down-regulated. Conclusions The high replacement of fish oil by vegetable oils in practical fish feeds did not modify the intestine transcriptome of gilthead sea bream, but important changes were apparent when fish were exposed to the myxosporean E. leei. The detected changes were mostly a consequence rather than a cause of the different disease progression in the two diet groups. Hence, the developed microarray constitutes an excellent diagnostic tool to address changes associated with the action of intestinal pathogens, but lacks a prognostic value to predict in advance the different susceptibility of growing fish to the current pathogen.
Collapse
Affiliation(s)
- Josep A Calduch-Giner
- Nutrigenomics and Fish Growth Endocrinology Group, Department of Marine Species Biology, Culture and Pathology, Instituto de Acuicultura Torre de la Sal (IATS-CSIC), Ribera de Cabanes, Castellón, 12595, Spain
| | | | | | | | | | | |
Collapse
|