1
|
Rangel-Chávez CP, Galán-Vásquez E, Pescador-Tapia A, Delaye L, Martínez-Antonio A. RNA polymerases in strict endosymbiont bacteria with extreme genome reduction show distinct erosions that might result in limited and differential promoter recognition. PLoS One 2021; 16:e0239350. [PMID: 34324516 PMCID: PMC8321222 DOI: 10.1371/journal.pone.0239350] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2020] [Accepted: 06/22/2021] [Indexed: 11/26/2022] Open
Abstract
Strict endosymbiont bacteria present high degree genome reduction, retain smaller proteins, and in some instances, lack complete functional domains compared to free-living counterparts. Until now, the mechanisms underlying these genetic reductions are not well understood. In this study, the conservation of RNA polymerases, the essential machinery for gene expression, is analyzed in endosymbiont bacteria with extreme genome reductions. We analyzed the RNA polymerase subunits to identify and define domains, subdomains, and specific amino acids involved in precise biological functions known in Escherichia coli. We also perform phylogenetic analysis and three-dimensional models over four lineages of endosymbiotic proteobacteria with the smallest genomes known to date: Candidatus Hodgkinia cicadicola, Candidatus Tremblaya phenacola, Candidatus Tremblaya Princeps, Candidatus Nasuia deltocephalinicola, and Candidatus Carsonella ruddii. We found that some Hodgkinia strains do not encode for the RNA polymerase α subunit. The rest encode genes for α, β, β', and σ subunits to form the RNA polymerase. However, 16% shorter, on average, respect their orthologous in E. coli. In the α subunit, the amino-terminal domain is the most conserved. Regarding the β and β' subunits, both the catalytic core and the assembly domains are the most conserved. However, they showed compensatory amino acid substitutions to adapt to changes in the σ subunit. Precisely, the most erosive diversity occurs within the σ subunit. We identified broad amino acid substitution even in those recognizing and binding to the -10-box promoter element. In an overall conceptual image, the RNA polymerase from Candidatus Nasuia conserved the highest similarity with Escherichia coli RNA polymerase and their σ70. It might be recognizing the two main promoter elements (-10 and -35) and the two promoter accessory elements (-10 extended and UP-element). In Candidatus Carsonella, the RNA polymerase could recognize all the promoter elements except the -10-box extended. In Candidatus Tremblaya and Hodgkinia, due to the α carboxyl-terminal domain absence, they might not recognize the UP-promoter element. We also identified the lack of the β flap-tip helix domain in most Hodgkinia's that suggests the inability to bind the -35-box promoter element.
Collapse
Affiliation(s)
- Cynthia Paola Rangel-Chávez
- Biological Engineering Laboratory, Genetic Engineering Department, Center for Research and Advanced Studies of the National Polytechnic Institute, Irapuato Gto, México
| | - Edgardo Galán-Vásquez
- Departamento de Ingeniería de Sistemas Computacionales y Automatización, Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Coyoacán, Ciudad de México, CDMX, México
| | - Azucena Pescador-Tapia
- Biological Engineering Laboratory, Genetic Engineering Department, Center for Research and Advanced Studies of the National Polytechnic Institute, Irapuato Gto, México
| | - Luis Delaye
- Evolutionary Genomics Laboratory, Genetic Engineering Department, Center for Research and Advanced Studies of the National Polytechnic Institute, Irapuato Gto, México
| | - Agustino Martínez-Antonio
- Biological Engineering Laboratory, Genetic Engineering Department, Center for Research and Advanced Studies of the National Polytechnic Institute, Irapuato Gto, México
| |
Collapse
|
2
|
Hadariová L, Vesteg M, Birčák E, Schwartzbach SD, Krajčovič J. An intact plastid genome is essential for the survival of colorless Euglena longa but not Euglena gracilis. Curr Genet 2016; 63:331-341. [PMID: 27553633 DOI: 10.1007/s00294-016-0641-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2016] [Revised: 08/02/2016] [Accepted: 08/13/2016] [Indexed: 12/18/2022]
Abstract
Euglena gracilis growth with antibacterial agents leads to bleaching, permanent plastid gene loss. Colorless Euglena (Astasia) longa resembles a bleached E. gracilis. To evaluate the role of bleaching in E. longa evolution, the effect of streptomycin, a plastid protein synthesis inhibitor, and ofloxacin, a plastid DNA gyrase inhibitor, on E. gracilis and E. longa growth and plastid DNA content were compared. E. gracilis growth was unaffected by streptomycin and ofloxacin. Quantitative PCR analyses revealed a time dependent loss of plastid genes in E. gracilis demonstrating that bleaching agents produce plastid gene deletions without affecting cell growth. Streptomycin and ofloxacin inhibited E. longa growth indicating that it requires plastid genes to survive. This suggests that evolutionary divergence of E. longa from E. gracilis was triggered by the loss of a cytoplasmic metabolic activity also occurring in the plastid. Plastid metabolism has become obligatory for E. longa cell growth. A process termed "intermittent bleaching", short term exposure to subsaturating concentrations of reversible bleaching agents followed by growth in the absence of a bleaching agent, is proposed as the molecular mechanism for E. longa plastid genome reduction. Various non-photosynthetic lineages could have independently arisen from their photosynthetic ancestors via a similar process.
Collapse
Affiliation(s)
- Lucia Hadariová
- Department of Genetics, Faculty of Natural Sciences, Comenius University, Mlynská dolina G-1, 842 15, Bratislava, Slovak Republic
| | - Matej Vesteg
- Department of Biology and Ecology, Faculty of Natural Sciences, Matej Bel University, 974 01, Banská Bystrica, Slovakia
| | - Erik Birčák
- Department of Genetics, Faculty of Natural Sciences, Comenius University, Mlynská dolina G-1, 842 15, Bratislava, Slovak Republic
| | | | - Juraj Krajčovič
- Department of Genetics, Faculty of Natural Sciences, Comenius University, Mlynská dolina G-1, 842 15, Bratislava, Slovak Republic. .,Department of Biology, Faculty of Natural Sciences, University of ss. Cyril and Methodius, 917 01, Trnava, Slovakia.
| |
Collapse
|
3
|
Kim JI, Linton EW, Shin W. Taxon-rich multigene phylogeny of the photosynthetic euglenoids (Euglenophyceae). Front Ecol Evol 2015. [DOI: 10.3389/fevo.2015.00098] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
|
4
|
Lee WM, An YJ. Evidence of three-level trophic transfer of quantum dots in an aquatic food chain by using bioimaging. Nanotoxicology 2014; 9:407-12. [PMID: 25119416 DOI: 10.3109/17435390.2014.948517] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
In this study, we demonstrated the three-level trophic transfer of quantum dots (QDs) within the aquatic food chain. Using bioimaging, we observed QD transfer from protozoa (Astasia longa) to zooplankton (Moina macrocopa) to fish (Danio rerio). Bioimaging is an effective tool that can improve our understanding of the delivery of nanomaterials in vivo. Measurement with an intravital multiphoton laser scanning microscope visually proved the transfer of QDs from the first to the second and the second to the third levels. As QDs may be passed from lower organisms to humans via the food chain, our findings have implications for the safety of their use.
Collapse
Affiliation(s)
- Woo-Mi Lee
- Department of Environmental Health Science, Konkuk University , 1 Hwayang-dong, Gwangjin-gu, Seoul , Korea
| | | |
Collapse
|
5
|
Wiegert KE, Bennett MS, Triemer RE. Tracing patterns of chloroplast evolution in euglenoids: contributions from Colacium vesiculosum and Strombomonas acuminata (Euglenophyta). J Eukaryot Microbiol 2013; 60:214-21. [PMID: 23351081 DOI: 10.1111/jeu.12025] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2012] [Revised: 10/16/2012] [Accepted: 12/10/2012] [Indexed: 11/29/2022]
Abstract
The chloroplast genomes of two photosynthetic euglenoids, Colacium vesiculosum Ehrenberg (128,889 bp), and Strombomonas acuminata (Schmarda) Deflandre (144,167 bp) have been sequenced. These chloroplast genomes in combination with those of Euglena gracilis, Eutreptia viridis, and Eutreptiella gymnastica provide a snapshot of euglenoid chloroplast evolution allowing comparisons of gene content, arrangement, and expansion. The gene content of the five chloroplast genomes is very similar varying only in the presence or absence of, rrn5, roaA, psaI, psaM, rpoA, and two tRNAs. Large gene rearrangements have occurred within the C. vesiculosum and S. acuminata chloroplast genomes. Most of these rearrangements represent repositioning of entire operons rather than single genes. When compared with previously sequenced genomes, C. vesiculosum and S. acuminata chloroplast genomes more closely resemble the E. gracilis chloroplast genome in size of the genome, number of introns, and gene order than they do those of the Eutreptiales. Overall, the chloroplast genomes of these five species show an evolutionary trend toward increased intron number, a decrease in gene density, and substantial rearrangement of gene clusters.
Collapse
Affiliation(s)
- Krystle E Wiegert
- Department of Plant Biology, Michigan State University, 612 Wilson Road, East Lansing, MI 48824, USA
| | | | | |
Collapse
|
6
|
Evidence for transitional stages in the evolution of euglenid group II introns and twintrons in the Monomorphina aenigmatica plastid genome. PLoS One 2012; 7:e53433. [PMID: 23300929 PMCID: PMC3534033 DOI: 10.1371/journal.pone.0053433] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2012] [Accepted: 11/28/2012] [Indexed: 11/19/2022] Open
Abstract
Background Photosynthetic euglenids acquired their plastid by secondary endosymbiosis of a prasinophyte-like green alga. But unlike its prasinophyte counterparts, the plastid genome of the euglenid Euglena gracilis is riddled with introns that interrupt almost every protein-encoding gene. The atypical group II introns and twintrons (introns-within-introns) found in the E. gracilis plastid have been hypothesized to have been acquired late in the evolution of euglenids, implying that massive numbers of introns may be lacking in other taxa. This late emergence was recently corroborated by the plastid genome sequences of the two basal euglenids, Eutreptiella gymnastica and Eutreptia viridis, which were found to contain fewer introns. Methodology/Principal Findings To gain further insights into the proliferation of introns in euglenid plastids, we have characterized the complete plastid genome sequence of Monomorphina aenigmatica, a freshwater species occupying an intermediate phylogenetic position between early and late branching euglenids. The M. aenigmatica UTEX 1284 plastid genome (74,746 bp, 70.6% A+T, 87 genes) contains 53 intron insertion sites, of which 41 were found to be shared with other euglenids including 12 of the 15 twintron insertion sites reported in E. gracilis. Conclusions The pattern of insertion sites suggests an ongoing but uneven process of intron gain in the lineage, with perhaps a minimum of two bursts of rapid intron proliferation. We also identified several sites that represent intermediates in the process of twintron evolution, where the external intron is in place, but not the internal one, offering a glimpse into how these convoluted molecular contraptions originate.
Collapse
|
7
|
Transcriptome analysis of the Euglena gracilis plastid chromosome. Curr Genet 2009; 55:425-38. [DOI: 10.1007/s00294-009-0256-8] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2009] [Revised: 05/12/2009] [Accepted: 05/15/2009] [Indexed: 11/26/2022]
|
8
|
Kim E, Archibald JM. Diversity and Evolution of Plastids and Their Genomes. PLANT CELL MONOGRAPHS 2008. [DOI: 10.1007/978-3-540-68696-5_1] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
9
|
Krause K. From chloroplasts to “cryptic” plastids: evolution of plastid genomes in parasitic plants. Curr Genet 2008; 54:111-21. [DOI: 10.1007/s00294-008-0208-8] [Citation(s) in RCA: 74] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2008] [Revised: 07/21/2008] [Accepted: 07/22/2008] [Indexed: 10/21/2022]
|
10
|
Kim JI, Shin W. PHYLOGENY OF THE EUGLENALES INFERRED FROM PLASTID LSU rDNA SEQUENCES(1). JOURNAL OF PHYCOLOGY 2008; 44:994-1000. [PMID: 27041618 DOI: 10.1111/j.1529-8817.2008.00536.x] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
To gain insights into the phylogeny of the Euglenales, we analyzed the plastid LSU rDNA sequences from 101 strains of the photosynthetic euglenoids belonging to nine ingroup genera (Euglena, Trachelomonas, Strombomonas, Monomorphina, Cryptoglena, Colacium, Discoplastis, Phacus, and Lepocinclis) and two outgroup genera (Eutreptia and Eutreptiella). Bayesian and maximum-likelihood (ML) analyses resulted in trees of similar topologies and four major clades: a Phacus and Lepocinclis clade; a Colacium clade; a Trachelomonas, Strombomonas, Monomorphina, and Cryptoglena clade; and a Euglena clade. The Phacus and Lepocinclis clade was the sister group of all other euglenalian genera, followed by Discoplastis spathirhyncha (Skuja) Triemer and the Colacium clade, respectively, which was inconsistent with their placement based on nuclear rDNA genes. The Trachelomonas, Strombomonas, Monomorphina, and Cryptoglena clade was sister to the Euglena clade. The loricate genera, Trachelomonas and Strombomonas, were closely related to each other, while Monomorphina and Cryptoglena also grouped together. The Euglena clade formed a monophyletic lineage comprising most species from taxa formerly allocated to the subgenera Calliglena and Euglena. However, within this genus, none of the subgenera was monophyletic.
Collapse
Affiliation(s)
- Jong Im Kim
- Department of Biology, Chungnam National University, Daejeon 305-764, South Korea
| | - Woongghi Shin
- Department of Biology, Chungnam National University, Daejeon 305-764, South Korea
| |
Collapse
|
11
|
|
12
|
Transcription and transcriptional regulation in plastids. CELL AND MOLECULAR BIOLOGY OF PLASTIDS 2007. [DOI: 10.1007/4735_2007_0232] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
13
|
Chumley TW, Palmer JD, Mower JP, Fourcade HM, Calie PJ, Boore JL, Jansen RK. The Complete Chloroplast Genome Sequence of Pelargonium × hortorum: Organization and Evolution of the Largest and Most Highly Rearranged Chloroplast Genome of Land Plants. Mol Biol Evol 2006; 23:2175-90. [PMID: 16916942 DOI: 10.1093/molbev/msl089] [Citation(s) in RCA: 308] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The chloroplast genome of Pelargonium x hortorum has been completely sequenced. It maps as a circular molecule of 217,942 bp and is both the largest and most rearranged land plant chloroplast genome yet sequenced. It features 2 copies of a greatly expanded inverted repeat (IR) of 75,741 bp each and, consequently, diminished single-copy regions of 59,710 and 6,750 bp. Despite the increase in size and complexity of the genome, the gene content is similar to that of other angiosperms, with the exceptions of a large number of pseudogenes, the recognition of 2 open reading frames (ORF56 and ORF42) in the trnA intron with similarities to previously identified mitochondrial products (ACRS and pvs-trnA), the losses of accD and trnT-ggu and, in particular, the presence of a highly divergent set of rpoA-like ORFs rather than a single, easily recognized gene for rpoA. The 3-fold expansion of the IR (relative to most angiosperms) accounts for most of the size increase of the genome, but an additional 10% of the size increase is related to the large number of repeats found. The Pelargonium genome contains 35 times as many 31 bp or larger repeats than the unrearranged genome of Spinacia. Most of these repeats occur near the rearrangement hotspots, and 2 different associations of repeats are localized in these regions. These associations are characterized by full or partial duplications of several genes, most of which appear to be nonfunctional copies or pseudogenes. These duplications may also be linked to the disruption of at least 1 but possibly 2 or 3 operons. We propose simple models that account for the major rearrangements with a minimum of 8 IR boundary changes and 12 inversions in addition to several insertions of duplicated sequence.
Collapse
|
14
|
Abhiman S, Sonnhammer ELL. Large-scale prediction of function shift in protein families with a focus on enzymatic function. Proteins 2005; 60:758-68. [PMID: 16001403 DOI: 10.1002/prot.20550] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Protein function shift can be predicted from sequence comparisons, either using positive selection signals or evolutionary rate estimation. None of the methods have been validated on large datasets, however. Here we investigate existing and novel methods for protein function shift prediction, and benchmark the accuracy against a large dataset of proteins with known enzymatic functions. Function change was predicted between subfamilies by identifying two kinds of sites in a multiple sequence alignment: Conservation-Shifting Sites (CSS), which are conserved in two subfamilies using two different amino acid types, and Rate-Shifting Sites (RSS), which have different evolutionary rates in two subfamilies. CSS were predicted by a new entropy-based method, and RSS using the Rate-Shift program. In principle, the more CSS and RSS between two subfamilies, the more likely a function shift between them. A test dataset was built by extracting subfamilies from Pfam with different EC numbers that belong to the same domain family. Subfamilies were generated automatically using a phylogenetic tree-based program, BETE. The dataset comprised 997 subfamily pairs with four or more members per subfamily. We observed a significant increase in CSS and RSS for subfamily comparisons with different EC numbers compared to cases with same EC numbers. The discrimination was better using RSS than CSS, and was more pronounced for larger families. Combining RSS and CSS by discriminant analysis improved classification accuracy to 71%. The method was applied to the Pfam database and the results are available at http://FunShift.cgb.ki.se. A closer examination of some superfamily comparisons showed that single EC numbers sometimes embody distinct functional classes. Hence, the measured accuracy of function shift is underestimated.
Collapse
Affiliation(s)
- Saraswathi Abhiman
- Center for Genomics and Bioinformatics, Karolinska Institutet, Stockholm, Sweden
| | | |
Collapse
|
15
|
Cavalier-Smith T. Genomic reduction and evolution of novel genetic membranes and protein-targeting machinery in eukaryote-eukaryote chimaeras (meta-algae). Philos Trans R Soc Lond B Biol Sci 2003; 358:109-33; discussion 133-4. [PMID: 12594921 PMCID: PMC1693104 DOI: 10.1098/rstb.2002.1194] [Citation(s) in RCA: 130] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Chloroplasts originated just once, from cyanobacteria enslaved by a biciliate protozoan to form the plant kingdom (green plants, red and glaucophyte algae), but subsequently, were laterally transferred to other lineages to form eukaryote-eukaryote chimaeras or meta-algae. This process of secondary symbiogenesis (permanent merger of two phylogenetically distinct eukaryote cells) has left remarkable traces of its evolutionary role in the more complex topology of the membranes surrounding all non-plant (meta-algal) chloroplasts. It took place twice, soon after green and red algae diverged over 550 Myr ago to form two independent major branches of the eukaryotic tree (chromalveolates and cabozoa), comprising both meta-algae and numerous secondarily non-photosynthetic lineages. In both cases, enslavement probably began by evolving a novel targeting of endomembrane vesicles to the perialgal vacuole to implant host porter proteins for extracting photosynthate. Chromalveolates arose by such enslavement of a unicellular red alga and evolution of chlorophyll c to form the kingdom Chromista and protozoan infrakingdom Alveolata, which diverged from the ancestral chromalveolate chimaera. Cabozoa arose when the common ancestor of euglenoids and cercozoan chlorarachnean algae enslaved a tetraphyte green alga with chlorophyll a and b. I suggest that in cabozoa the endomembrane vesicles originally budded from the Golgi, whereas in chromalveolates they budded from the endoplasmic reticulum (ER) independently of Golgi-targeted vesicles, presenting a potentially novel target for drugs against alveolate Sporozoa such as malaria parasites and Toxoplasma. These hypothetical ER-derived vesicles mediated fusion of the perialgal vacuole and rough ER (RER) in the ancestral chromist, placing the former red alga within the RER lumen. Subsequently, this chimaera diverged to form cryptomonads, which retained the red algal nucleus as a nucleomorph (NM) with approximately 464 protein-coding genes (30 encoding plastid proteins) and a red or blue phycobiliprotein antenna pigment, and the chromobiotes (heterokonts and haptophytes), which lost phycobilins and evolved the brown carotenoid fucoxanthin that colours brown seaweeds, diatoms and haptophytes. Chromobiotes transferred the 30 genes to the nucleus and lost the NM genome and nuclear-pore complexes, but retained its membrane as the periplastid reticulum (PPR), putatively the phospholipid factory of the periplastid space (former algal cytoplasm), as did the ancestral alveolate independently. The chlorarachnean NM has three minute chromosomes bearing approximately 300 genes riddled with pygmy introns. I propose that the periplastid membrane (PPM, the former algal plasma membrane) of chromalveolates, and possibly chlorarachneans, grows by fusion of vesicles emanating from the NM envelope or PPR. Dinoflagellates and euglenoids independently lost the PPM and PPR (after diverging from Sporozoa and chlorarachneans, respectively) and evolved triple chloroplast envelopes comprising the original plant double envelope and an extra outermost membrane, the EM, derived from the perialgal vacuole. In all metaalgae most chloroplast proteins are coded by nuclear genes and enter the chloroplast by using bipartite targeting sequences--an upstream signal sequence for entering the ER and a downstream chloroplast transit sequence. I present a new theory for the four-fold diversification of the chloroplast OM protein translocon following its insertion into the PPM to facilitate protein translocation across it (of both periplastid and plastid proteins). I discuss evidence from genome sequencing and other sources on the contrasting modes of protein targeting, cellular integration, and evolution of these two major lineages of eukaryote "cells within cells". They also provide powerful evidence for natural selection's effectiveness in eliminating most functionless DNA and therefore of a universally useful non-genic function for nuclear non-coding DNA, i.e. most DNA in the biosphere, and dramatic examples of genomic reduction. I briefly argue that chloroplast replacement in dinoflagellates, which happened at least twice, may have been evolutionarily easier than secondary symbiogenesis because parts of the chromalveolate protein-targeting machinery could have helped enslave the foreign plastids.
Collapse
Affiliation(s)
- T Cavalier-Smith
- Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, UK.
| |
Collapse
|