1
|
Molecular and Physiological Responses of Citrus sinensis Leaves to Long-Term Low pH Revealed by RNA-Seq Integrated with Targeted Metabolomics. Int J Mol Sci 2022; 23:ijms23105844. [PMID: 35628662 PMCID: PMC9142915 DOI: 10.3390/ijms23105844] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Revised: 05/18/2022] [Accepted: 05/21/2022] [Indexed: 12/30/2022] Open
Abstract
Low pH-induced alterations in gene expression profiles and organic acids (OA) and free amino acid (FAA) abundances were investigated in sweet orange [Citrus sinensis (L.) Osbeck cv. Xuegan] leaves. We identified 503 downregulated and 349 upregulated genes in low pH-treated leaves. Further analysis indicated that low pH impaired light reaction and carbon fixation in photosynthetic organisms, thereby lowering photosynthesis in leaves. Low pH reduced carbon and carbohydrate metabolisms, OA biosynthesis and ATP production in leaves. Low pH downregulated the biosynthesis of nitrogen compounds, proteins, and FAAs in leaves, which might be conducive to maintaining energy homeostasis during ATP deprivation. Low pH-treated leaves displayed some adaptive responses to phosphate starvation, including phosphate recycling, lipid remodeling, and phosphate transport, thus enhancing leaf acid-tolerance. Low pH upregulated the expression of some reactive oxygen species (ROS) and aldehyde detoxifying enzyme (peroxidase and superoxidase) genes and the concentrations of some antioxidants (L-tryptophan, L-proline, nicotinic acid, pantothenic acid, and pyroglutamic acid), but it impaired the pentose phosphate pathway and VE and secondary metabolite biosynthesis and downregulated the expression of some ROS and aldehyde detoxifying enzyme (ascorbate peroxidase, aldo-keto reductase, and 2-alkenal reductase) genes and the concentrations of some antioxidants (pyridoxine and γ-aminobutyric acid), thus disturbing the balance between production and detoxification of ROS and aldehydes and causing oxidative damage to leaves.
Collapse
|
2
|
Tissue-specific enhancement of OsRNS1 with root-preferred expression is required for the increase of crop yield. J Adv Res 2022; 42:69-81. [PMID: 35609869 PMCID: PMC9788951 DOI: 10.1016/j.jare.2022.05.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Revised: 05/03/2022] [Accepted: 05/17/2022] [Indexed: 12/27/2022] Open
Abstract
INTRODUCTION Root development is a fundamental process that supports plant survival and crop productivity. One of the essential factors to consider when developing biotechnology crops is the selection of a promoter that can optimize the spatial-temporal expression of introduced genes. However, there are insufficient cases of suitable promoters in crop plants, including rice. OBJECTIVES This study aimed to verify the usefulness of a new rice root-preferred promoter to optimize the function of a target gene with root-preferred expression in rice. METHODS osrns1 mutant had defects in root development based on T-DNA insertional mutant screening and CRISPR technology. To optimize the function of OsRNS1, we generated OsRNS1-overexpression plants under two different promoters: a whole-plant expression promoter and a novel root-preferred expression promoter. Root growth, yield-related agronomic traits, RNA-seq, and reactive oxygen species (ROS) accumulation were analyzed for comparison. RESULTS OsRNS1 was found to be involved in root development through T-DNA insertional mutant analysis and gene editing mutant analysis. To understand the gain of function of OsRNS1, pUbi1::OsRNS1 was generated for the whole-plant expression, and both root growth defects and overall growth defects were found. To overcome this problem, a root-preferential overexpression line using Os1-CysPrxB promoter (Per) was generated and showed an increase in root length, plant height, and grain yield compared to wild-type (WT). RNA-seq analysis revealed that the response to oxidative stress-related genes was significantly up-regulated in both overexpression lines but was more obvious in pPer::OsRNS1. Furthermore, ROS levels in the roots were drastically decreased in pPer::OsRNS1 but were increased in the osrns1 mutants compared to WT. CONCLUSION The results demonstrated that the use of a root-preferred promoter effectively optimizes the function of OsRNS1 and is a useful strategy for improving root-related agronomic traits as well as ROS regulation.
Collapse
|
3
|
Genome-wide identification and expression pattern analysis of the ribonuclease T2 family in Eucommia ulmoides. Sci Rep 2021; 11:6900. [PMID: 33767357 PMCID: PMC7994793 DOI: 10.1038/s41598-021-86337-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Accepted: 03/11/2021] [Indexed: 11/13/2022] Open
Abstract
The 2′,3′-cycling ribonuclease (RNase) genes are catalysts of RNA cleavage and include the RNase T2 gene family. RNase T2 genes perform important roles in plants and have been conserved in the genome of eukaryotic organisms. In this study we identified 21 EURNS genes in Eucommia ulmoides Oliver (E. ulmoides) and analyzed their structure, chromosomal location, phylogenetic tree, gene duplication, stress-related cis-elements, and expression patterns in different tissues. The length of 21 predicted EURNS proteins ranged from 143 to 374 amino acids (aa), their molecular weight (MW) ranged from 16.21 to 42.38 kDa, and their isoelectric point (PI) value ranged from 5.08 to 9.09. Two classifications (class I and class III) were obtained from the conserved domains analysis and phylogenetic tree. EURNS proteins contained a total of 15 motifs. Motif 1, motif 2, motif 3, and motif 7 were distributed in multiple sequences and were similar to the conserved domain of RNase T2. EURNS genes with similar structure and the predicted EURNS proteins with conserved motif compositions are in the same group in the phylogenetic tree. The results of RT-PCR and transcription data showed that EURNS genes have tissue-specific expression and exhibited obvious trends in different developmental stages. Gene duplication analysis results indicated that segment duplication may be the dominant duplication mode in this gene family. This study provides a theoretical basis for research on the RNase T2 gene family and lays a foundation for the further study of EURNS genes.
Collapse
|
4
|
Gho YS, Choi H, Moon S, Song MY, Park HE, Kim DH, Ha SH, Jung KH. Phosphate-Starvation-Inducible S-Like RNase Genes in Rice Are Involved in Phosphate Source Recycling by RNA Decay. FRONTIERS IN PLANT SCIENCE 2020; 11:585561. [PMID: 33424882 PMCID: PMC7793952 DOI: 10.3389/fpls.2020.585561] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2020] [Accepted: 11/02/2020] [Indexed: 05/16/2023]
Abstract
The fine-tuning of inorganic phosphate (Pi) for enhanced use efficiency has long been a challenging subject in agriculture, particularly in regard to rice as a major crop plant. Among ribonucleases (RNases), the RNase T2 family is broadly distributed across kingdoms, but little has been known on its substrate specificity compared to RNase A and RNase T1 families. Class I and class II of the RNase T2 family are defined as the S-like RNase (RNS) family and have showed the connection to Pi recycling in Arabidopsis. In this study, we first carried out a phylogenetic analysis of eight rice and five Arabidopsis RNS genes and identified mono-specific class I and dicot-specific class I RNS genes, suggesting the possibility of functional diversity between class I RNS family members in monocot and dicot species through evolution. We then compared the in silico expression patterns of all RNS genes in rice and Arabidopsis under normal and Pi-deficient conditions and further confirmed the expression patterns of rice RNS genes via qRT-PCR analysis. Subsequently, we found that most of the OsRNS genes were differentially regulated under Pi-deficient treatment. Association of Pi recycling by RNase activity in rice was confirmed by measuring total RNA concentration and ribonuclease activity of shoot and root samples under Pi-sufficient or Pi-deficient treatment during 21 days. The total RNA concentrations were decreased by < 60% in shoots and < 80% in roots under Pi starvation, respectively, while ribonuclease activity increased correspondingly. We further elucidate the signaling pathway of Pi starvation through upregulation of the OsRNS genes. The 2-kb promoter region of all OsRNS genes with inducible expression patterns under Pi deficiency contains a high frequency of P1BS cis-acting regulatory element (CRE) known as the OsPHR2 binding site, suggesting that the OsRNS family is likely to be controlled by OsPHR2. Finally, the dynamic transcriptional regulation of OsRNS genes by overexpression of OsPHR2, ospho2 mutant, and overexpression of OsPT1 lines involved in Pi signaling pathway suggests the molecular basis of OsRNS family in Pi recycling via RNA decay under Pi starvation.
Collapse
Affiliation(s)
- Yun-Shil Gho
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| | - Heebak Choi
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| | - Sunok Moon
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| | - Min Yeong Song
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| | - Ha Eun Park
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| | - Doh-Hoon Kim
- Department of Life Science, College of Life Science and Natural Resources, Dong-A University, Busan, South Korea
| | - Sun-Hwa Ha
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| | - Ki-Hong Jung
- Graduate School of Biotechnology, Kyung Hee University, Yongin, South Korea
| |
Collapse
|
5
|
Diaz-Baena M, Galvez-Valdivieso G, Delgado-Garcia E, Pineda M, Piedras P. Nuclease and ribonuclease activities in response to salt stress: Identification of PvRNS3, a T2/S-like ribonuclease induced in common bean radicles by salt stress. PLANT PHYSIOLOGY AND BIOCHEMISTRY : PPB 2020; 147:235-241. [PMID: 31881432 DOI: 10.1016/j.plaphy.2019.12.016] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Revised: 12/12/2019] [Accepted: 12/13/2019] [Indexed: 05/28/2023]
Abstract
The increase in soil salinization due to global climate change could cause large losses in crop productivity affecting, among other biological processes, to germination and seedling development. We have studied how salt stress affects nucleic acid degrading activities in radicles of common bean during seedling development. In radicles of common bean, a main nuclease of 37 kDa and two ribonucleases of 17 and 19 kDa were detected. Saline stress did not alter these three activities but induced a new ribonuclease of 16 kDa. All three ribonucleases are acidic enzymes that were inhibited by Zn. The 16 and 17 kDa ribonucleases are inhibited by guanilates. In the genome of common bean, we have identified 13 genes belonging to the T2 ribonuclease family and that are grouped in the 3 classes of T2 ribonucleases. The analysis of the expression of the 3 genes belonging to Class I (PvRNS1 to 3) and the unique gene from Class II (PvRNS4) in radicles showed that PvRNS3 is highly induced under salt stress.
Collapse
Affiliation(s)
- Mercedes Diaz-Baena
- Departamento de Botánica, Ecología y Fisiología Vegetal, Plants Molecular Physiology and Biotechnology Group, Campus de Rabanales, Edif. Severo Ochoa, Universidad de Córdoba, Córdoba, Spain
| | - Gregorio Galvez-Valdivieso
- Departamento de Botánica, Ecología y Fisiología Vegetal, Plants Molecular Physiology and Biotechnology Group, Campus de Rabanales, Edif. Severo Ochoa, Universidad de Córdoba, Córdoba, Spain
| | - Elena Delgado-Garcia
- Departamento de Botánica, Ecología y Fisiología Vegetal, Plants Molecular Physiology and Biotechnology Group, Campus de Rabanales, Edif. Severo Ochoa, Universidad de Córdoba, Córdoba, Spain
| | - Manuel Pineda
- Departamento de Botánica, Ecología y Fisiología Vegetal, Plants Molecular Physiology and Biotechnology Group, Campus de Rabanales, Edif. Severo Ochoa, Universidad de Córdoba, Córdoba, Spain
| | - Pedro Piedras
- Departamento de Botánica, Ecología y Fisiología Vegetal, Plants Molecular Physiology and Biotechnology Group, Campus de Rabanales, Edif. Severo Ochoa, Universidad de Córdoba, Córdoba, Spain.
| |
Collapse
|
6
|
Nguyen-Kim H, San Clemente H, Balliau T, Zivy M, Dunand C, Albenne C, Jamet E. Arabidopsis thaliana
root cell wall proteomics: Increasing the proteome coverage using a combinatorial peptide ligand library and description of unexpected Hyp in peroxidase amino acid sequences. Proteomics 2016; 16:491-503. [DOI: 10.1002/pmic.201500129] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2015] [Revised: 10/07/2015] [Accepted: 11/10/2015] [Indexed: 01/12/2023]
Affiliation(s)
- Huan Nguyen-Kim
- Laboratoire de Recherche en Sciences Végétales, UMR 5546, UPS, Université de Toulouse; BP 42617 Castanet-Tolosan France
- UMR 5546; CNRS; BP 42617 Castanet-Tolosan France
| | - Hélène San Clemente
- Laboratoire de Recherche en Sciences Végétales, UMR 5546, UPS, Université de Toulouse; BP 42617 Castanet-Tolosan France
- UMR 5546; CNRS; BP 42617 Castanet-Tolosan France
| | - Thierry Balliau
- CNRS; PAPPSO; UMR 0320/UMR 8120 Génétique Végétale Quantitative et Evolution; Le Moulon Gif sur Yvette France
- INRA; PAPPSO; UMR 0320/UMR 8120 Génétique Végétale Quantitative et Evolution; Le Moulon Gif sur Yvette France
| | - Michel Zivy
- CNRS; PAPPSO; UMR 0320/UMR 8120 Génétique Végétale Quantitative et Evolution; Le Moulon Gif sur Yvette France
- INRA; PAPPSO; UMR 0320/UMR 8120 Génétique Végétale Quantitative et Evolution; Le Moulon Gif sur Yvette France
| | - Christophe Dunand
- Laboratoire de Recherche en Sciences Végétales, UMR 5546, UPS, Université de Toulouse; BP 42617 Castanet-Tolosan France
- UMR 5546; CNRS; BP 42617 Castanet-Tolosan France
| | - Cécile Albenne
- Laboratoire de Recherche en Sciences Végétales, UMR 5546, UPS, Université de Toulouse; BP 42617 Castanet-Tolosan France
- UMR 5546; CNRS; BP 42617 Castanet-Tolosan France
| | - Elisabeth Jamet
- Laboratoire de Recherche en Sciences Végétales, UMR 5546, UPS, Université de Toulouse; BP 42617 Castanet-Tolosan France
- UMR 5546; CNRS; BP 42617 Castanet-Tolosan France
| |
Collapse
|
7
|
Rojas H, Floyd B, Morriss SC, Bassham D, MacIntosh GC, Goldraij A. NnSR1, a class III non-S-RNase specifically induced in Nicotiana alata under phosphate deficiency, is localized in endoplasmic reticulum compartments. PLANT SCIENCE : AN INTERNATIONAL JOURNAL OF EXPERIMENTAL PLANT BIOLOGY 2015; 236:250-9. [PMID: 26025538 DOI: 10.1016/j.plantsci.2015.04.012] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2015] [Revised: 04/11/2015] [Accepted: 04/18/2015] [Indexed: 05/28/2023]
Abstract
A combined strategy of phosphate (Pi) remobilization from internal and external RNA sources seems to be conserved in plants exposed to Pi starvation. Thus far, the only ribonucleases (RNases) reported to be induced in Nicotiana alata undergoing Pi deprivation are extracellular S-like RNase NE and NnSR1. NnSR1 is a class III non S-RNase of unknown subcellular location. Here, we examine the hypothesis that NnSR1 is an intracellular RNase derived from the self-incompatibility system with specific expression in self-incompatible Nicotiana alata. NnSR1 was not induced in self-compatible Nicotiana species exposed to Pi deprivation. NnSR1 conjugated with a fluorescent protein and transiently expressed in Arabidopsis protoplasts and Nicotiana leaves showed that the fusion protein co-localized with an endoplasmic reticulum (ER) marker. Subcellular fractionation by ultracentrifugation of roots exposed to Pi deprivation revealed that the native NnSR1 migrated in parallel with the BiP protein, a typical ER marker. To our knowledge, NnSR1 is the first class III RNase reported to be localized in ER compartments. The induction of NnSR1 was detected earlier than the extracellular RNase NE, suggesting that intracellular RNA may be the first source of Pi used by the cell under Pi stress.
Collapse
Affiliation(s)
- Hernán Rojas
- Dpto Química Biológica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, 5000 Córdoba, Argentina
| | - Brice Floyd
- Dept of Genetics, Development, and Cell Biology, Iowa State University, Ames, IA 50011, USA
| | - Stephanie C Morriss
- Dept of Biochemistry, Biophysics & Molecular Biology, Iowa State University, Ames, IA 50011, USA
| | - Diane Bassham
- Dept of Genetics, Development, and Cell Biology, Iowa State University, Ames, IA 50011, USA
| | - Gustavo C MacIntosh
- Dept of Biochemistry, Biophysics & Molecular Biology, Iowa State University, Ames, IA 50011, USA.
| | - Ariel Goldraij
- Dpto Química Biológica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, 5000 Córdoba, Argentina.
| |
Collapse
|
8
|
Aguiar B, Vieira J, Cunha AE, Fonseca NA, Iezzoni A, van Nocker S, Vieira CP. Convergent evolution at the gametophytic self-incompatibility system in Malus and Prunus. PLoS One 2015; 10:e0126138. [PMID: 25993016 PMCID: PMC4438004 DOI: 10.1371/journal.pone.0126138] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2014] [Accepted: 03/30/2015] [Indexed: 12/24/2022] Open
Abstract
S-RNase-based gametophytic self-incompatibility (GSI) has evolved once before the split of the Asteridae and Rosidae. This conclusion is based on the phylogenetic history of the S-RNase that determines pistil specificity. In Rosaceae, molecular characterizations of Prunus species, and species from the tribe Pyreae (i.e., Malus, Pyrus, Sorbus) revealed different numbers of genes determining S-pollen specificity. In Prunus only one pistil and pollen gene determine GSI, while in Pyreae there is one pistil but multiple pollen genes, implying different specificity recognition mechanisms. It is thus conceivable that within Rosaceae the genes involved in GSI in the two lineages are not orthologous but possibly paralogous. To address this hypothesis we characterised the S-RNase lineage and S-pollen lineage genes present in the genomes of five Rosaceae species from three genera: M. × domestica (apple, self-incompatible (SI); tribe Pyreae), P. persica (peach, self-compatible (SC); Amygdaleae), P. mume (mei, SI; Amygdaleae), Fragaria vesca (strawberry, SC; Potentilleae), and F. nipponica (mori-ichigo, SI; Potentilleae). Phylogenetic analyses revealed that the Malus and Prunus S-RNase and S-pollen genes belong to distinct gene lineages, and that only Prunus S-RNase and SFB-lineage genes are present in Fragaria. Thus, S-RNase based GSI system of Malus evolved independently from the ancestral system of Rosaceae. Using expression patterns based on RNA-seq data, the ancestral S-RNase lineage gene is inferred to be expressed in pistils only, while the ancestral S-pollen lineage gene is inferred to be expressed in tissues other than pollen.
Collapse
Affiliation(s)
- Bruno Aguiar
- Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
- Instituto de Biologia Molecular e Celular (IBMC), Universidade do Porto, Porto, Portugal
| | - Jorge Vieira
- Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
- Instituto de Biologia Molecular e Celular (IBMC), Universidade do Porto, Porto, Portugal
| | - Ana E. Cunha
- Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
- Instituto de Biologia Molecular e Celular (IBMC), Universidade do Porto, Porto, Portugal
| | - Nuno A. Fonseca
- CRACS-INESC Porto, Rua do Campo Alegre 1021/1055, 4169–007, Porto, Portugal
- European Bioinformatics Institute (EMBL-EBI), Welcome Trust Genome Campus, CB10 1SD, Cambridge, United Kingdom
| | - Amy Iezzoni
- Michigan State University, East Lansing, Michigan, United States of America
| | - Steve van Nocker
- Michigan State University, East Lansing, Michigan, United States of America
| | - Cristina P. Vieira
- Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
- Instituto de Biologia Molecular e Celular (IBMC), Universidade do Porto, Porto, Portugal
- * E-mail:
| |
Collapse
|
9
|
Nishimura E, Jumyo S, Arai N, Kanna K, Kume M, Nishikawa JI, Tanase JI, Ohyama T. Structural and functional characteristics of S-like ribonucleases from carnivorous plants. PLANTA 2014; 240:147-59. [PMID: 24771022 DOI: 10.1007/s00425-014-2072-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2014] [Accepted: 03/24/2014] [Indexed: 05/09/2023]
Abstract
Although the S-like ribonucleases (RNases) share sequence homology with the S-RNases involved in the self-incompatibility mechanism in plants, they are not associated with this mechanism. They usually function in stress responses in non-carnivorous plants and in carnivory in carnivorous plants. In this study, we clarified the structures of the S-like RNases of Aldrovanda vesiculosa, Nepenthes bicalcarata and Sarracenia leucophylla, and compared them with those of other plants. At ten positions, amino acid residues are conserved or almost conserved only for carnivorous plants (six in total). In contrast, two positions are specific to non-carnivorous plants. A phylogenetic analysis revealed that the S-like RNases of the carnivorous plants form a group beyond the phylogenetic relationships of the plants. We also prepared and characterized recombinant S-like RNases of Dionaea muscipula, Cephalotus follicularis, A. vesiculosa, N. bicalcarata and S. leucophylla, and RNS1 of Arabidopsis thaliana. The recombinant carnivorous plant enzymes showed optimum activities at about pH 4.0. Generally, poly(C) was digested less efficiently than poly(A), poly(I) and poly(U). The kinetic parameters of the recombinant D. muscipula enzyme (DM-I) and A. thaliana enzyme RNS1 were similar. The k cat/K m of recombinant RNS1 was the highest among the enzymes, followed closely by that of recombinant DM-I. On the other hand, the k cat/K m of the recombinant S. leucophylla enzyme was the lowest, and was ~1/30 of that for recombinant RNS1. The magnitudes of the k cat/K m values or k cat values for carnivorous plant S-like RNases seem to correlate negatively with the dependency on symbionts for prey digestion.
Collapse
Affiliation(s)
- Emi Nishimura
- Major in Integrative Bioscience and Biomedical Engineering, Graduate School of Science and Engineering, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo, 162-8480, Japan
| | | | | | | | | | | | | | | |
Collapse
|