1
|
Nardi F, Boschi S, Funari R, Cucini C, Cardaioli E, Potter D, Asano SI, Toubarro D, Meier M, Paoli F, Carapelli A, Frati F. The direction, timing and demography of Popillia japonica (Coleoptera) invasion reconstructed using complete mitochondrial genomes. Sci Rep 2024; 14:7120. [PMID: 38531924 DOI: 10.1038/s41598-024-57667-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 03/20/2024] [Indexed: 03/28/2024] Open
Abstract
The Japanese beetle Popillia japonica is a pest insect that feeds on hundreds of species of wild and cultivated plants including important fruit, vegetable, and field crops. Native to Japan, the pest has invaded large areas of the USA, Canada, the Azores (Portugal), Italy, and Ticino (Switzerland), and it is considered a priority for control in the European Union. We determined the complete mitochondrial genome sequence in 86 individuals covering the entire distribution of the species. Phylogenetic analysis supports a major division between South Japan and Central/North Japan, with invasive samples coming from the latter. The origin of invasive USA samples is incompatible, in terms of the timing of the event, with a single introduction, with multiple Japanese lineages having been introduced and one accounting for most of the population expansion locally. The origin of the two invasive European populations is compatible with two different invasions followed by minimal differentiation locally. Population analyses provide the possibility to estimate the rate of sequence change from the data and to date major invasion events. Demographic analysis identifies a population expansion followed by a period of contraction prior to the invasion. The present study adds a time and demographic dimension to available reconstructions.
Collapse
Affiliation(s)
- Francesco Nardi
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy.
- NBFC, National Biodiversity Future Center, Palermo, Italy.
- BAT-Center, Interuniversity Center for Studies on Bioinspired Agro-Environmental Technology, Portici, Italy.
| | - Sara Boschi
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy
| | - Rebecca Funari
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy
| | - Claudio Cucini
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy
| | - Elena Cardaioli
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy
| | - Daniel Potter
- Department of Entomology, University of Kentucky, Lexington, USA
| | - Shin-Ichiro Asano
- Research Faculty of Agriculture, Hokkaido University, Sapporo, Japan
| | - Duarte Toubarro
- Biotechnology Centre of Azores, University of the Azores, Ponta Delgada, Portugal
| | - Michela Meier
- Servizio fitosanitario cantonale, Dipartimento delle finanze e dell'economia, Bellinzona, Switzerland
| | - Francesco Paoli
- Council for Agricultural Research and Agricultural Economy Analysis (CREA), Florence, Italy
| | - Antonio Carapelli
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy
- NBFC, National Biodiversity Future Center, Palermo, Italy
- BAT-Center, Interuniversity Center for Studies on Bioinspired Agro-Environmental Technology, Portici, Italy
| | - Francesco Frati
- Department of Life Sciences, University of Siena, via Aldo Moro 2, 53100, Siena, Italy
- NBFC, National Biodiversity Future Center, Palermo, Italy
- BAT-Center, Interuniversity Center for Studies on Bioinspired Agro-Environmental Technology, Portici, Italy
| |
Collapse
|
2
|
Bringloe TT, Parent GJ. Contrasting new and available reference genomes to highlight uncertainties in assemblies and areas for future improvement: an example with monodontid species. BMC Genomics 2023; 24:693. [PMID: 37985969 PMCID: PMC10659057 DOI: 10.1186/s12864-023-09779-3] [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/07/2023] [Accepted: 10/31/2023] [Indexed: 11/22/2023] Open
Abstract
BACKGROUND Reference genomes provide a foundational framework for evolutionary investigations, ecological analysis, and conservation science, yet uncertainties in the assembly of reference genomes are difficult to assess, and by extension rarely quantified. Reference genomes for monodontid cetaceans span a wide spectrum of data types and analytical approaches, providing the context to derive broader insights related to discrepancies and regions of uncertainty in reference genome assembly. We generated three beluga (Delphinapterus leucas) and one narwhal (Monodon monoceros) reference genomes and contrasted these with published chromosomal scale assemblies for each species to quantify discrepancies associated with genome assemblies. RESULTS The new reference genomes achieved chromosomal scale assembly using a combination of PacBio long reads, Illumina short reads, and Hi-C scaffolding data. For beluga, we identified discrepancies in the order and orientation of contigs in 2.2-3.7% of the total genome depending on the pairwise comparison of references. In addition, unsupported higher order scaffolding was identified in published reference genomes. In contrast, we estimated 8.2% of the compared narwhal genomes featured discrepancies, with inversions being notably abundant (5.3%). Discrepancies were linked to repetitive elements in both species. CONCLUSIONS We provide several new reference genomes for beluga (Delphinapterus leucas), while highlighting potential avenues for improvements. In particular, additional layers of data providing information on ultra-long genomic distances are needed to resolve persistent errors in reference genome construction. The comparative analyses of monodontid reference genomes suggested that the three new reference genomes for beluga are more accurate compared to the currently published reference genome, but that the new narwhal genome is less accurate than one published. We also present a conceptual summary for improving the accuracy of reference genomes with relevance to end-user needs and how they relate to levels of assembly quality and uncertainty.
Collapse
Affiliation(s)
- Trevor T Bringloe
- Laboratory of Genomics, Maurice Lamontagne Institute, Fisheries and Oceans Canada, Mont-Joli, QC, Canada.
| | - Geneviève J Parent
- Laboratory of Genomics, Maurice Lamontagne Institute, Fisheries and Oceans Canada, Mont-Joli, QC, Canada.
| |
Collapse
|
3
|
Morin PA, Martien KK, Lang AR, Hancock-Hanser BL, Pease VL, Robertson KM, Sattler M, Slikas E, Rosel PE, Baker CS, Taylor BL, Archer FI. Guidelines and quantitative standards for improved cetacean taxonomy using full mitochondrial genomes. J Hered 2023; 114:612-624. [PMID: 37647537 DOI: 10.1093/jhered/esad049] [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/31/2023] [Accepted: 08/28/2023] [Indexed: 09/01/2023] Open
Abstract
In many organisms, especially those of conservation concern, traditional lines of evidence for taxonomic delineation, such as morphological data, are often difficult to obtain. In these cases, genetic data are often the only source of information available for taxonomic studies. In particular, population surveys of mitochondrial genomes offer increased resolution and precision in support of taxonomic decisions relative to conventional use of the control region or other gene fragments of the mitochondrial genome. To improve quantitative guidelines for taxonomic decisions in cetaceans, we build on a previous effort targeting the control region and evaluate, for whole mitogenome sequences, a suite of divergence and diagnosability estimates for pairs of recognized cetacean populations, subspecies, and species. From this overview, we recommend new guidelines based on complete mitogenomes, combined with other types of evidence for isolation and divergence, which will improve resolution for taxonomic decisions, especially in the face of small sample sizes or low levels of genetic diversity. We further use simulated data to assist interpretations of divergence in the context of varying forms of historical demography, culture, and ecology.
Collapse
Affiliation(s)
- Phillip A Morin
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Karen K Martien
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Aimee R Lang
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Brittany L Hancock-Hanser
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Victoria L Pease
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Kelly M Robertson
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Maya Sattler
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Elizabeth Slikas
- School of Life Sciences, Center for Evolution and Medicine, Arizona State University, Tempe, AZ, United States
| | - Patricia E Rosel
- Marine Mammal and Turtle Division, Southeast Fisheries Science Center, National Marine Fisheries Service, NOAA, Lafayette, LA, United States
| | - C Scott Baker
- Marine Mammal Institute, Oregon State University, Newport, OR, United States
| | - Barbara L Taylor
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| | - Frederick I Archer
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, La Jolla, CA, United States
| |
Collapse
|
4
|
Heraghty SD, Jackson JM, Lozier JD. Whole genome analyses reveal weak signatures of population structure and environmentally associated local adaptation in an important North American pollinator, the bumble bee Bombus vosnesenskii. Mol Ecol 2023; 32:5479-5497. [PMID: 37702957 DOI: 10.1111/mec.17125] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 08/21/2023] [Accepted: 08/24/2023] [Indexed: 09/14/2023]
Abstract
Studies of species that experience environmental heterogeneity across their distributions have become an important tool for understanding mechanisms of adaptation and predicting responses to climate change. We examine population structure, demographic history and environmentally associated genomic variation in Bombus vosnesenskii, a common bumble bee in the western USA, using whole genome resequencing of populations distributed across a broad range of latitudes and elevations. We find that B. vosnesenskii exhibits minimal population structure and weak isolation by distance, confirming results from previous studies using other molecular marker types. Similarly, demographic analyses with Sequentially Markovian Coalescent models suggest that minimal population structure may have persisted since the last interglacial period, with genomes from different parts of the species range showing similar historical effective population size trajectories and relatively small fluctuations through time. Redundancy analysis revealed a small amount of genomic variation explained by bioclimatic variables. Environmental association analysis with latent factor mixed modelling (LFMM2) identified few outlier loci that were sparsely distributed throughout the genome and although a few putative signatures of selective sweeps were identified, none encompassed particularly large numbers of loci. Some outlier loci were in genes with known regulatory relationships, suggesting the possibility of weak selection, although compared with other species examined with similar approaches, evidence for extensive local adaptation signatures in the genome was relatively weak. Overall, results indicate B. vosnesenskii is an example of a generalist with a high degree of flexibility in its environmental requirements that may ultimately benefit the species under periods of climate change.
Collapse
Affiliation(s)
- Sam D Heraghty
- Department of Biological Sciences, The University of Alabama, Tuscaloosa, Alabama, USA
| | - Jason M Jackson
- Department of Biological Sciences, The University of Alabama, Tuscaloosa, Alabama, USA
| | - Jeffrey D Lozier
- Department of Biological Sciences, The University of Alabama, Tuscaloosa, Alabama, USA
| |
Collapse
|
5
|
Durante CA, Loizaga R, McCracken GR, Crespo EA, Ruzzante DE. Commerson's dolphin population structure: evidence for female phylopatry and male dispersal. Sci Rep 2022; 12:22219. [PMID: 36564404 PMCID: PMC9789068 DOI: 10.1038/s41598-022-26192-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 12/12/2022] [Indexed: 12/24/2022] Open
Abstract
A key in species conservation is understanding the amount and distribution of genetic diversity and how environmental changes that occurred in the recent past may have influenced current patterns of population structure. Commerson's dolphin, Cephalorhynchus commersonii, has two subspecies, one of which is endemic to South America (C. commersonii commersonii) and little is known about its population genetics. Our objective was to investigate the population genetics of this subspecies throughout its distribution. Using 70 skin samples and information available in GenBank, 308 mitochondrial DNA sequences and 28 species-specific microsatellites were analyzed. The species presented low genetic diversity when compared to other dolphin species, but was consistent with other species within the genus. Strong population structure based on mitochondrial DNA was exhibited throughout its entire distribution, a pattern consistent with female philopatry. However, this pattern was not detected when using microsatellites, suggesting male-mediated gene flow. Demographic tests suggested a population expansion beginning approximately 15,000 years ago, after the Last Glacial Maximum. In a climate change scenario, we recommended considering each sampling location as an independent population management unit in order to evaluate the impact of possible environmental changes on the distribution of genetic information within the species.
Collapse
Affiliation(s)
- Cristian Alberto Durante
- Laboratorio de Mamíferos Marinos, Centro Para El Estudio de Sistemas Marinos (CESIMAR) - CONICET, Bv. Brown 2915, U9120ACD Puerto Madryn, Chubut, Argentina.
- Universidad Nacional del Comahue, San Carlos de Bariloche, Rio Negro, Argentina.
| | - Rocio Loizaga
- Laboratorio de Mamíferos Marinos, Centro Para El Estudio de Sistemas Marinos (CESIMAR) - CONICET, Bv. Brown 2915, U9120ACD Puerto Madryn, Chubut, Argentina
| | | | - Enrique Alberto Crespo
- Laboratorio de Mamíferos Marinos, Centro Para El Estudio de Sistemas Marinos (CESIMAR) - CONICET, Bv. Brown 2915, U9120ACD Puerto Madryn, Chubut, Argentina
| | | |
Collapse
|
6
|
Kebke A, Samarra F, Derous D. Climate change and cetacean health: impacts and future directions. Philos Trans R Soc Lond B Biol Sci 2022; 377:20210249. [PMID: 35574848 DOI: 10.1098/rstb.2021.0249] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Climate change directly impacts the foraging opportunities of cetaceans (e.g. lower prey availability), leads to habitat loss, and forces cetaceans to move to other feeding grounds. The rise in ocean temperature, low prey availability and loss of habitat can have severe consequences for cetacean survival, particularly those species that are already threatened or those with a limited habitat range. In addition, it is predicted that the concentration of contaminants in aquatic environments will increase owing to Arctic meltwater and increased rainfall events leading to higher rates of land-based runoff in downstream coastal areas. These persistent and mobile contaminants can bioaccumulate in the ecosystem, and lead to ecotoxicity with potentially severe consequences on the reproductive organs, immune system and metabolism of marine mammals. There is a need to measure and assess the cumulative impact of multiple stressors, given that climate change, habitat alteration, low prey availability and contaminants do not act in isolation. Human-caused perturbations to cetacean foraging abilities are becoming a pervasive and prevalent threat to many cetacean species on top of climate change-associated stressors. We need to move to a greater understanding of how multiple stressors impact the metabolism of cetaceans and ultimately their population trajectory. This article is part of the theme issue 'Nurturing resilient marine ecosystems'.
Collapse
Affiliation(s)
- Anna Kebke
- School of Biological Sciences, University of Aberdeen, Aberdeen, UK
| | - Filipa Samarra
- University of Iceland's Institute of Research Centres, Vestmannaeyjar, Iceland
| | - Davina Derous
- School of Biological Sciences, University of Aberdeen, Aberdeen, UK
| |
Collapse
|
7
|
Janiak MC, Silva FE, Beck RMD, de Vries D, Kuderna LFK, Torosin NS, Melin AD, Marquès-Bonet T, Goodhead IB, Messias M, da Silva MNF, Sampaio I, Farias IP, Rossi R, de Melo FR, Valsecchi J, Hrbek T, Boubli JP. 205 newly assembled mitogenomes provide mixed evidence for rivers as drivers of speciation for Amazonian primates. Mol Ecol 2022; 31:3888-3902. [PMID: 35638312 PMCID: PMC9546496 DOI: 10.1111/mec.16554] [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: 01/25/2022] [Revised: 05/06/2022] [Accepted: 05/20/2022] [Indexed: 12/02/2022]
Abstract
Mitochondrial DNA remains a cornerstone for molecular ecology, especially for study species from which high‐quality tissue samples cannot be easily obtained. Methods using mitochondrial markers are usually reliant on reference databases, but these are often incomplete. Furthermore, available mitochondrial genomes often lack crucial metadata, such as sampling location, limiting their utility for many analyses. Here, we assembled 205 new mitochondrial genomes for platyrrhine primates, most from the Amazon and with known sampling locations. We present a dated mitogenomic phylogeny based on these samples along with additional published platyrrhine mitogenomes, and use this to assess support for the long‐standing riverine barrier hypothesis (RBH), which proposes that river formation was a major driver of speciation in Amazonian primates. Along the Amazon, Negro, and Madeira rivers, we found mixed support for the RBH. While we identified divergences that coincide with a river barrier, only some occur synchronously and also overlap with the proposed dates of river formation. The most compelling evidence is for the Amazon river potentially driving speciation within bearded saki monkeys (Chiropotes spp.) and within the smallest extant platyrrhines, the marmosets and tamarins. However, we also found that even large rivers do not appear to be barriers for some primates, including howler monkeys (Alouatta spp.), uakaris (Cacajao spp.), sakis (Pithecia spp.), and robust capuchins (Sapajus spp.). Our results support a more nuanced, clade‐specific effect of riverine barriers and suggest that other evolutionary mechanisms, besides the RBH and allopatric speciation, may have played an important role in the diversification of platyrrhines.
Collapse
Affiliation(s)
- Mareike C Janiak
- School of Science, Engineering & Environment, University of Salford, Salford, United Kingdom
| | - Felipe E Silva
- Research Group on Primate Biology and Conservation, Mamirauá Institute for Sustainable Development, Brazil.,Department of Evolutionary Biology and Ecology, Université Libre de Bruxelles, Belgium
| | - Robin M D Beck
- School of Science, Engineering & Environment, University of Salford, Salford, United Kingdom
| | - Dorien de Vries
- School of Science, Engineering & Environment, University of Salford, Salford, United Kingdom
| | - Lukas F K Kuderna
- Ilumina Inc., Hayward, CA, USA.,Institute of Evolutionary Biology (UPF-CSIC), Barcelona, Spain
| | - Nicole S Torosin
- Department of Genetics, Human Genetics Institute of New Jersey, Rutgers University, Piscataway, NJ, USA
| | - Amanda D Melin
- Department of Anthropology & Archaeology and Department of Medical Genetics, University of Calgary, Calgary, AB, Canada.,Alberta Children's Hospital Research Institute, Calgary, AB, Canada
| | - Tomàs Marquès-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), Barcelona, Spain.,Catalan Institution of Research and Advanced Studies (ICREA), Barcelona, Spain.,CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona, Spain.,Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Edifici ICTA-ICP, Cerdanyola del Vallès, Barcelona, Spain
| | - Ian B Goodhead
- School of Science, Engineering & Environment, University of Salford, Salford, United Kingdom
| | - Mariluce Messias
- Department of Biology, Universidade Federal de Rondônia, Porto Velho, Brazil
| | - Maria N F da Silva
- Coleção de Mamíferos, Instituto Nacional de Pesquisas da Amazônia, Brazil
| | | | - Izeni P Farias
- Laboratory of Evolution and Animal Genetics, Universidade Federal do Amazonas, Brazil
| | - Rogerio Rossi
- Instituto de Biociências, Universidade Federal do Mato Grosso, Brazil
| | - Fabiano R de Melo
- Department of Forestry Engineering, Universidade Federal de Viçosa, Brazil
| | - João Valsecchi
- Research Group on Primate Biology and Conservation, Mamirauá Institute for Sustainable Development, Brazil
| | - Tomas Hrbek
- Department of Biology, Trinity University, San Antonio, TX, USA
| | - Jean P Boubli
- School of Science, Engineering & Environment, University of Salford, Salford, United Kingdom.,Coleção de Mamíferos, Instituto Nacional de Pesquisas da Amazônia, Brazil
| |
Collapse
|
8
|
Ferguson SH, Yurkowski DJ, Hudson JM, Edkins T, Willing C, Watt CA. Larger body size leads to greater female beluga whale ovarian reproductive activity at the southern periphery of their range. Ecol Evol 2021; 11:17314-17322. [PMID: 34938510 PMCID: PMC8668808 DOI: 10.1002/ece3.8367] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 10/05/2021] [Accepted: 11/02/2021] [Indexed: 12/22/2022] Open
Abstract
Identification of phenotypic characteristics in reproductively successful individuals provides important insights into the evolutionary processes that cause range shifts due to environmental change. Female beluga whales (Delphinapterus leucas) from the Baffin Bay region (BB) of the Canadian Arctic in the core area of the species' geographic range have larger body size than their conspecifics at the southern range periphery in Hudson Bay (HB). We investigated the mechanism for this north and south divergence as it relates to ovarian reproductive activity (ORA = total corpora) that combines morphometric data with ovarian corpora counted from female reproductive tracts. Our study aim was to assess the relative influence of age and body size of female beluga whale on ORA in the two populations. Female beluga whale ORA increased more quickly with age (63% partial variation explained) in BB than in HB (41%). In contrast, body length in HB female beluga whales accounted for considerably more of the total variation (12% vs. 1%) in ORA compared to BB whales. We speculate that female HB beluga whale ORA was more strongly linked with body length due to higher population density resulting in food competition that favors the energetic advantages of larger body size during seasonal food limitations. Understanding the evolutionary mechanism of how ORA varies across a species' range will assist conservation efforts in anticipating and mitigating future challenges associated with a warming planet.
Collapse
Affiliation(s)
- Steven H. Ferguson
- Fisheries and Oceans CanadaFreshwater InstituteWinnipegManitobaCanada
- Biological SciencesUniversity of ManitobaWinnipegManitobaCanada
| | - David J. Yurkowski
- Fisheries and Oceans CanadaFreshwater InstituteWinnipegManitobaCanada
- Biological SciencesUniversity of ManitobaWinnipegManitobaCanada
| | - Justine M. Hudson
- Fisheries and Oceans CanadaFreshwater InstituteWinnipegManitobaCanada
| | - Tera Edkins
- Fisheries and Oceans CanadaFreshwater InstituteWinnipegManitobaCanada
| | | | - Cortney A. Watt
- Fisheries and Oceans CanadaFreshwater InstituteWinnipegManitobaCanada
- Biological SciencesUniversity of ManitobaWinnipegManitobaCanada
| |
Collapse
|