1
|
Samad‐Zamini M, Schweiger W, Nussbaumer T, Mayer KF, Buerstmayr H. Time-course expression QTL-atlas of the global transcriptional response of wheat to Fusarium graminearum. PLANT BIOTECHNOLOGY JOURNAL 2017; 15:1453-1464. [PMID: 28332274 PMCID: PMC5633761 DOI: 10.1111/pbi.12729] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2016] [Revised: 01/11/2017] [Accepted: 03/16/2017] [Indexed: 05/09/2023]
Abstract
Fusarium head blight is a devastating disease of small grain cereals such as bread wheat (Triticum aestivum). The pathogen switches from a biotrophic to a nectrotrophic lifestyle in course of disease development forcing its host to adapt its defence strategies. Using a genetical genomics approach, we illustrate genome-wide reconfigurations of genetic control over transcript abundances between two decisive time points after inoculation with the causative pathogen Fusarium graminearum. Whole transcriptome measurements have been recorded for 163 lines of a wheat doubled haploid population segregating for several resistance genes yielding 15 552 at 30 h and 15 888 eQTL at 50 h after inoculation. The genetic map saturated with transcript abundance-derived markers identified of a novel QTL on chromosome 6A, besides the previously reported QTL Fhb1 and Qfhs.ifa-5A. We find a highly different distribution of eQTL between time points with about 40% of eQTL being unique for the respective assessed time points. But also for more than 20% of genes governed by eQTL at either time point, genetic control changes in time. These changes are reflected in the dynamic compositions of three major regulatory hotspots on chromosomes 2B, 4A and 5A. In particular, control of defence-related biological mechanisms concentrated in the hotspot at 4A shift to hotspot 2B as the disease progresses. Hotspots do not colocalize with phenotypic QTL, and within their intervals no higher than expected number of eQTL was detected. Thus, resistance conferred by either QTL is mediated by few or single genes.
Collapse
Affiliation(s)
- Mina Samad‐Zamini
- Institute for Biotechnology in Plant Production (IFA‐Tulln)BOKU ‐ University of Natural Resources and Life SciencesTullnAustria
| | - Wolfgang Schweiger
- Institute for Biotechnology in Plant Production (IFA‐Tulln)BOKU ‐ University of Natural Resources and Life SciencesTullnAustria
- Present address:
BIOMIN Research CenterTulln3430Austria
| | - Thomas Nussbaumer
- Plant Genome and Systems BiologyHelmholtz Zentrum MünchenNeuherbergGermany
- Present address:
Division of Computational System BiologyDepartment of Microbiology and Ecosystem ScienceUniversity of ViennaVienna1090Austria
| | - Klaus F.X. Mayer
- Plant Genome and Systems BiologyHelmholtz Zentrum MünchenNeuherbergGermany
| | - Hermann Buerstmayr
- Institute for Biotechnology in Plant Production (IFA‐Tulln)BOKU ‐ University of Natural Resources and Life SciencesTullnAustria
| |
Collapse
|
2
|
Karafiátová M, Bartoš J, Kopecký D, Ma L, Sato K, Houben A, Stein N, Doležel J. Mapping nonrecombining regions in barley using multicolor FISH. Chromosome Res 2013; 21:739-51. [DOI: 10.1007/s10577-013-9380-x] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2013] [Revised: 08/26/2013] [Accepted: 08/30/2013] [Indexed: 12/22/2022]
|
3
|
Michalak de Jimenez MK, Bassi FM, Ghavami F, Simons K, Dizon R, Seetan RI, Alnemer LM, Denton AM, Doğramacı M, Šimková H, Doležel J, Seth K, Luo MC, Dvorak J, Gu YQ, Kianian SF. A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus. Funct Integr Genomics 2013. [PMID: 23479086 DOI: 10.1007/s10142‐013‐0318‐3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."
Collapse
|
4
|
Michalak de Jimenez MK, Bassi FM, Ghavami F, Simons K, Dizon R, Seetan RI, Alnemer LM, Denton AM, Doğramacı M, Šimková H, Doležel J, Seth K, Luo MC, Dvorak J, Gu YQ, Kianian SF. A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus. Funct Integr Genomics 2013; 13:19-32. [PMID: 23479086 DOI: 10.1007/s10142-013-0318-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2012] [Revised: 02/08/2013] [Accepted: 02/12/2013] [Indexed: 10/27/2022]
Abstract
The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."
Collapse
|
5
|
Construction of whole genome radiation hybrid panels and map of chromosome 5A of wheat using asymmetric somatic hybridization. PLoS One 2012; 7:e40214. [PMID: 22815731 PMCID: PMC3398029 DOI: 10.1371/journal.pone.0040214] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2011] [Accepted: 06/06/2012] [Indexed: 11/23/2022] Open
Abstract
To explore the feasibility of constructing a whole genome radiation hybrid (WGRH) map in plant species with large genomes, asymmetric somatic hybridization between wheat (Triticum aestivum L.) and Bupleurum scorzonerifolium Willd. was performed. The protoplasts of wheat were irradiated with ultraviolet light (UV) and gamma-ray and rescued by protoplast fusion using B. scorzonerifolium as the recipient. Assessment of SSR markers showed that the radiation hybrids have the average marker retention frequency of 15.5%. Two RH panels (RHPWI and RHPWII) that contained 92 and 184 radiation hybrids, respectively, were developed and used for mapping of 68 SSR markers in chromosome 5A of wheat. A total of 1557 and 2034 breaks were detected in each panel. The RH map of chromosome 5A based on RHPWII was constructed. The distance of the comprehensive map was 2103 cR and the approximate resolution was estimated to be ∼501.6 kb/break. The RH panels evaluated in this study enabled us to order the ESTs in a single deletion bin or in the multiple bins cross the chromosome. These results demonstrated that RH mapping via protoplast fusion is feasible at the whole genome level for mapping purposes in wheat and the potential value of this mapping approach for the plant species with large genomes.
Collapse
|
6
|
Kynast RG, Riera-Lizarazu O. Development and use of oat-maize chromosome additions and radiation hybrids. Methods Mol Biol 2011; 701:259-284. [PMID: 21181536 DOI: 10.1007/978-1-61737-957-4_15] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Hybridization experiments of oat with maize require fastidious coordination of plant cultivation and flowering timing, meticulous crossing techniques, stimulation with plant growth substances, and in vitro rescue and culture of the hybrid embryos. The majority of hybrid offspring gradually lose all maize chromosomes consequently resulting in haploid oat plants. However, a minority of the offspring retain one or more maize chromosome(s) in addition to their haploid oat complements (partial hybrids). Oat haploids and partial hybrids with 1-3 maize chromosomes are partially fertile. Controlled self-fertilization of partial hybrids allows for the production of doubled haploid oat plants with an added single maize chromosome (monosomic addition) or an added pair of homologous maize chromosomes (disomic addition) among the inbred offspring. γ-Irradiation of monosomic oat-maize addition lines can be used to further dissect the maize chromosome in a given line. The lines with identified maize chromosome fragments (radiation hybrids) are the basis for establishing chromosome-specific panels. Although still in the experimental phase, the use of radiation hybrids has been useful and has widened the repertoire of maize genetics and genomics methodology.
Collapse
Affiliation(s)
- Ralf G Kynast
- Jodrell Laboratory, Royal Botanic Gardens, Kew, Richmond Surrey, UK
| | | |
Collapse
|
7
|
A Comparative BAC map for the gilthead sea bream (Sparus aurata L.). J Biomed Biotechnol 2010; 2011:329025. [PMID: 21049003 PMCID: PMC2964914 DOI: 10.1155/2011/329025] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2010] [Revised: 09/10/2010] [Accepted: 09/23/2010] [Indexed: 12/18/2022] Open
Abstract
This study presents the first comparative BAC map of the gilthead sea bream (Sparus aurata), a highly valuated marine aquaculture fish species in the Mediterranean. High-throughput end sequencing of a BAC library yielded 92,468 reads (60.6 Mbp). Comparative mapping was achieved by anchoring BAC end sequences to the three-spined stickleback (Gasterosteus aculeatus) genome. BACs that were consistently ordered along the stickleback chromosomes accounted for 14,265 clones. A fraction of 5,249 BACs constituted a minimal tiling path that covers 73.5% of the stickleback chromosomes and 70.2% of the genes that have been annotated. The N50 size of 1,485 “BACtigs” consisting of redundant BACs is 337,253 bp. The largest BACtig covers 2.15 Mbp in the stickleback genome. According to the insert size distribution of mapped BACs the sea bream genome is 1.71-fold larger than the stickleback genome. These results represent a valuable tool to researchers in the field and may support future projects to elucidate the whole sea bream genome.
Collapse
|
8
|
Yamano S, Nitta M, Tsujimoto H, Ishikawa G, Nakamura T, Endo TR, Nasuda S. Molecular mapping of the suppressor gene Igc1 to the gametocidal gene Gc3-C1 in common wheat. Genes Genet Syst 2010; 85:43-53. [PMID: 20410664 DOI: 10.1266/ggs.85.43] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Several species of the genus Aegilops, wild relatives of wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) carry gametocidal (Gc) genes. Gc genes kill the gametes without themselves by causing chromosomal breakage during post-meiotic cell divisions, and therefore are strong segregation distorters. The Gc gene Gc3-C1 derived from chromosome 3C of Ae. triuncialis (2n = 4x = 28, CCUU) induces chromosomal breakage in wheat cultivar 'Chinese Spring' (CS) but not in cultivar 'Norin 26' (N26). This cultivar-specific inhibition of Gc function is caused by a suppressor gene Igc1 located on chromosome 3B of N26. Igc1 is presumed to be a modified Gc gene without breakage function because of its homoeology to Gc3-C1. Here we report the results of linkage and physical mapping of Igc1 to help elucidate the molecular mechanisms underlying Gc action. Segregation analysis of the phenotypic data in BC(1)F(1) mapping population of the cross between (CSxN26)F(1) and CS + 3C" showed a 1:1 segregation ratio indicating that Igc1 is a dominant gene. In the linkage analysis, three molecular marker loci Xgwm285, Xgwm376, and Xcfp1886 cosegregated with the Igc1 locus. Bin mapping assigned the loci Xgwm285 and Xcfp1886 to bin C-3BS1-0.33 and Xgwm376 to bin C-3BL2-0.22. Physical mapping using Gc-induced chromosomal deletion lines of chromosome 3B of N26 revealed that the Igc1 locus resides in 52.0% or 2.1% of bins C-3BS1-0.33 and C-3BL2-0.22, respectively. Pericentromeric localization of Igc1 in chromosome 3B of N26 may have a positive effect to keep the two-component system of the Gc action. Map-based cloning approach to isolate the Igc1 may be difficult because recombination is depleted in the pericentromeric region. As is shown in this study, the combination of genetic and physical mapping offers high efficiency to identify the regions where genes are located especially in regions with low levels of recombination.
Collapse
Affiliation(s)
- Soichi Yamano
- Laboratory of Plant Genetics, Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan
| | | | | | | | | | | | | |
Collapse
|
9
|
Zhou S, Wei F, Nguyen J, Bechner M, Potamousis K, Goldstein S, Pape L, Mehan MR, Churas C, Pasternak S, Forrest DK, Wise R, Ware D, Wing RA, Waterman MS, Livny M, Schwartz DC. A single molecule scaffold for the maize genome. PLoS Genet 2009; 5:e1000711. [PMID: 19936062 PMCID: PMC2774507 DOI: 10.1371/journal.pgen.1000711] [Citation(s) in RCA: 115] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2009] [Accepted: 10/05/2009] [Indexed: 11/18/2022] Open
Abstract
About 85% of the maize genome consists of highly repetitive sequences that are interspersed by low-copy, gene-coding sequences. The maize community has dealt with this genomic complexity by the construction of an integrated genetic and physical map (iMap), but this resource alone was not sufficient for ensuring the quality of the current sequence build. For this purpose, we constructed a genome-wide, high-resolution optical map of the maize inbred line B73 genome containing >91,000 restriction sites (averaging 1 site/∼23 kb) accrued from mapping genomic DNA molecules. Our optical map comprises 66 contigs, averaging 31.88 Mb in size and spanning 91.5% (2,103.93 Mb/∼2,300 Mb) of the maize genome. A new algorithm was created that considered both optical map and unfinished BAC sequence data for placing 60/66 (2,032.42 Mb) optical map contigs onto the maize iMap. The alignment of optical maps against numerous data sources yielded comprehensive results that proved revealing and productive. For example, gaps were uncovered and characterized within the iMap, the FPC (fingerprinted contigs) map, and the chromosome-wide pseudomolecules. Such alignments also suggested amended placements of FPC contigs on the maize genetic map and proactively guided the assembly of chromosome-wide pseudomolecules, especially within complex genomic regions. Lastly, we think that the full integration of B73 optical maps with the maize iMap would greatly facilitate maize sequence finishing efforts that would make it a valuable reference for comparative studies among cereals, or other maize inbred lines and cultivars. The maize genome contains abundant repeats interspersed by low-copy, gene-coding sequences that make it a challenge to sequence; consequently, current BAC sequence assemblies average 11 contigs per clone. The iMap deals with such complexity by the judicious integration of IBM genetic and B73 physical maps, but the B73 genome structure could differ from the IBM population because of genetic recombination and subsequent rearrangements. Accordingly, we report a genome-wide, high-resolution optical map of maize B73 genome that was constructed from the direct analysis of genomic DNA molecules without using genetic markers. The integration of optical and iMap resources with comparisons to FPC maps enabled a uniquely comprehensive and scalable assessment of a given BAC's sequence assembly, its placement within a FPC contig, and the location of this FPC contig within a chromosome-wide pseudomolecule. As such, the overall utility of the maize optical map for the validation of sequence assemblies has been significant and demonstrates the inherent advantages of single molecule platforms. Construction of the maize optical map represents the first physical map of a eukaryotic genome larger than 400 Mb that was created de novo from individual genomic DNA molecules.
Collapse
Affiliation(s)
- Shiguo Zhou
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Fusheng Wei
- Department of Plant Sciences, Arizona Genomics Institute, University of Arizona, Tucson, Arizona, United States of America
| | - John Nguyen
- Departments of Mathematics, Biology, and Computer Science, University of Southern California, Los Angeles, California, United States of America
| | - Mike Bechner
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Konstantinos Potamousis
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Steve Goldstein
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Louise Pape
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Michael R. Mehan
- Departments of Mathematics, Biology, and Computer Science, University of Southern California, Los Angeles, California, United States of America
| | - Chris Churas
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Shiran Pasternak
- Cold Spring Harbor Laboratory, Cold Spring Harbor, New York, United States of America
| | - Dan K. Forrest
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
| | - Roger Wise
- Corn Insects and Crop Genetics Research, United States Department of Agriculture–Agricultural Research Service and Department of Plant Pathology, Iowa State University, Ames, Iowa, United States of America
| | - Doreen Ware
- Cold Spring Harbor Laboratory, Cold Spring Harbor, New York, United States of America
- Plant, Soil, and Nutrition Research, United States Department of Agriculture–Agricultural Research Service, Ithaca, New York, United States of America
| | - Rod A. Wing
- Department of Plant Sciences, Arizona Genomics Institute, University of Arizona, Tucson, Arizona, United States of America
| | - Michael S. Waterman
- Departments of Mathematics, Biology, and Computer Science, University of Southern California, Los Angeles, California, United States of America
| | - Miron Livny
- Computer Sciences Department, University of Wisconsin-Madison, Madison, Wisconsin, United States of America
| | - David C. Schwartz
- Laboratory for Molecular and Computational Genomics, Department of Chemistry, Laboratory of Genetics, UW Biotechnology Center, University of Wisconsin–Madison, Madison, Wisconsin, United States of America
- * E-mail:
| |
Collapse
|
10
|
Jones N, Ougham H, Thomas H, Pašakinskienė I. Markers and mapping revisited: finding your gene. THE NEW PHYTOLOGIST 2009; 183:935-966. [PMID: 19594696 DOI: 10.1111/j.1469-8137.2009.02933.x] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
This paper is an update of our earlier review (Jones et al., 1997, Markers and mapping: we are all geneticists now. New Phytologist 137: 165-177), which dealt with the genetics of mapping, in terms of recombination as the basis of the procedure, and covered some of the first generation of markers, including restriction fragment length polymorphisms (RFLPs), random amplified polymorphic DNA (RAPDs), simple sequence repeats (SSRs) and quantitative trait loci (QTLs). In the intervening decade there have been numerous developments in marker science with many new systems becoming available, which are herein described: cleavage amplification polymorphism (CAP), sequence-specific amplification polymorphism (S-SAP), inter-simple sequence repeat (ISSR), sequence tagged site (STS), sequence characterized amplification region (SCAR), selective amplification of microsatellite polymorphic loci (SAMPL), single nucleotide polymorphism (SNP), expressed sequence tag (EST), sequence-related amplified polymorphism (SRAP), target region amplification polymorphism (TRAP), microarrays, diversity arrays technology (DArT), single-strand conformation polymorphism (SSCP), denaturing gradient gel electrophoresis (DGGE), temperature gradient gel electrophoresis (TGGE) and methylation-sensitive PCR. In addition there has been an explosion of knowledge and databases in the area of genomics and bioinformatics. The number of flowering plant ESTs is c. 19 million and counting, with all the opportunity that this provides for gene-hunting, while the survey of bioinformatics and computer resources points to a rapid growth point for future activities in unravelling and applying the burst of new information on plant genomes. A case study is presented on tracking down a specific gene (stay-green (SGR), a post-transcriptional senescence regulator) using the full suite of mapping tools and comparative mapping resources. We end with a brief speculation on how genome analysis may progress into the future of this highly dynamic arena of plant science.
Collapse
Affiliation(s)
- Neil Jones
- IBERS, Aberystwyth University, Edward Llwyd Building, Penglais Campus, Aberystwyth, Ceredigion SY23 3DA, UK
| | - Helen Ougham
- IBERS, Aberystwyth University, Gogerddan Campus, Aberystwyth, Ceredigion SY23 3EB, UK
| | - Howard Thomas
- IBERS, Aberystwyth University, Edward Llwyd Building, Penglais Campus, Aberystwyth, Ceredigion SY23 3DA, UK
| | - Izolda Pašakinskienė
- Botanical Garden of Vilnius University, Kairenu 43, LT-10239 Vilnius, Lithuania
- Faculty of Natural Sciences, Department of Botany and Genetics, MK Čiurlionio g. 21, LT-03101 Vilnius, Lithuania
| |
Collapse
|