1
|
du Plessis SJ, Blaxter M, Koepfli KP, Chadwick EA, Hailer F. Genomics Reveals Complex Population History and Unexpected Diversity of Eurasian Otters (Lutra lutra) in Britain Relative to Genetic Methods. Mol Biol Evol 2023; 40:msad207. [PMID: 37713621 PMCID: PMC10630326 DOI: 10.1093/molbev/msad207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Revised: 08/04/2023] [Accepted: 09/06/2023] [Indexed: 09/17/2023] Open
Abstract
Conservation genetic analyses of many endangered species have been based on genotyping of microsatellite loci and sequencing of short fragments of mtDNA. The increase in power and resolution afforded by whole genome approaches may challenge conclusions made on limited numbers of loci and maternally inherited haploid markers. Here, we provide a matched comparison of whole genome sequencing versus microsatellite and control region (CR) genotyping for Eurasian otters (Lutra lutra). Previous work identified four genetically differentiated "stronghold" populations of otter in Britain, derived from regional populations that survived the population crash of the 1950s-1980s. Using whole genome resequencing data from 45 samples from across the British stronghold populations, we confirmed some aspects of population structure derived from previous marker-driven studies. Importantly, we showed that genomic signals of the population crash bottlenecks matched evidence from otter population surveys. Unexpectedly, two strongly divergent mitochondrial lineages were identified that were undetectable using CR fragments, and otters in the east of England were genetically distinct and surprisingly variable. We hypothesize that this previously unsuspected variability may derive from past releases of Eurasian otters from other, non-British source populations in England around the time of the population bottleneck. Our work highlights that even reasonably well-studied species may harbor genetic surprises, if studied using modern high-throughput sequencing methods.
Collapse
Affiliation(s)
| | - Mark Blaxter
- Tree of Life, Wellcome Sanger Institute, Cambridge, UK
| | - Klaus-Peter Koepfli
- Smithsonian-Mason School of Conservation, George Mason University, Front Royal, VA, USA
- Centre for Species Survival, Smithsonian's National Zoo and Conservation Biology Institute, Washington, DC, USA
| | | | - Frank Hailer
- School of Biosciences, Cardiff University, Cardiff, UK
| |
Collapse
|
2
|
Lajmi A, Glinka F, Privman E. Optimizing ddRAD sequencing for population genomic studies with ddgRADer. Mol Ecol Resour 2023. [PMID: 37732396 DOI: 10.1111/1755-0998.13870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Revised: 08/04/2023] [Accepted: 08/31/2023] [Indexed: 09/22/2023]
Abstract
Double-digest Restriction-site Associated DNA sequencing (ddRADseq) is widely used to generate genomic data for non-model organisms in evolutionary and ecological studies. Along with affordable paired-end sequencing, this method makes population genomic analyses more accessible. However, multiple factors should be considered when designing a ddRADseq experiment, which can be challenging for new users. The generated data often suffer from substantial read overlaps and adaptor contamination, severely reducing sequencing efficiency and affecting data quality. Here, we analyse diverse datasets from the literature and carry out controlled experiments to understand the effects of enzyme choice and size selection on sequencing efficiency. The empirical data reveal that size selection is imprecise and has limited efficacy. In certain scenarios, a substantial proportion of short fragments pass below the lower size-selection cut-off resulting in low sequencing efficiency. However, enzyme choice can considerably mitigate inadvertent inclusion of these shorter fragments. A simple model based on these experiments is implemented to predict the number of genomic fragments generated after digestion and size selection, number of SNPs genotyped, number of samples that can be multiplexed and the expected sequencing efficiency. We developed ddgRADer - http://ddgrader.haifa.ac.il/ - a user-friendly webtool and incorporated these calculations to aid in ddRADseq experimental design while optimizing sequencing efficiency. This tool can also be used for single enzyme protocols such as Genotyping-by-Sequencing. Given user-defined study goals, ddgRADer recommends enzyme pairs and allows users to compare and choose enzymes and size-selection criteria. ddgRADer improves the accessibility and ease of designing ddRADseq experiments and increases the probability of success of the first population genomic study conducted in labs with no prior experience in genomics.
Collapse
Affiliation(s)
- Aparna Lajmi
- Department of Evolutionary and Environmental Biology, Institute of Evolution, University of Haifa, Haifa, Israel
| | - Felix Glinka
- Department of Evolutionary and Environmental Biology, Institute of Evolution, University of Haifa, Haifa, Israel
| | - Eyal Privman
- Department of Evolutionary and Environmental Biology, Institute of Evolution, University of Haifa, Haifa, Israel
| |
Collapse
|
3
|
Deka JR, Ali SZ, Ahamad M, Borah P, Gopi GV, Badola R, Sharma R, Hussain SA. Can Bengal Tiger ( Panthera tigris tigris) endure the future climate and land use change scenario in the East Himalayan Region? Perspective from a multiple model framework. Ecol Evol 2023; 13:e10340. [PMID: 37554398 PMCID: PMC10404654 DOI: 10.1002/ece3.10340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Revised: 07/07/2023] [Accepted: 07/10/2023] [Indexed: 08/10/2023] Open
Abstract
Large mammals are susceptible to land use and climate change, unless they are safeguarded within large, protected areas. It is crucial to comprehend the effects of these changes on mammals to develop a conservation plan. We identified ecological hotspots that can sustain an ecosystem for the endangered Bengal tiger (Panthera tigris tigris), an umbrella species. We developed three distinct ensemble species distribution models (SDMs) for the Bengal tiger in the Indian East Himalayan Region (IEHR). The first model served as the baseline and considered habitat type, climate, land cover, and anthropogenic threats. The second model focused on climate, land use, and anthropogenic threats, the third model focused on climate variables. We projected the second and third models onto two future climate scenarios: RCP 4.5 and RCP 8.5. We evaluated the threats possess to protected areas within eco-sensitive zone based on the potential tiger habitat. Finally, we compared the potential habitat with the IUCN tiger range. Our study revealed that the Brahmaputra valley will serve as the primary habitat for tigers in the future. However, considering the projected severe climate scenarios, it is anticipated that tigers will undergo a range shift towards the north and east, especially in high-altitude regions. Very high conservation priority areas, which make up 3.4% of the total area, are predominantly located in the riverine corridor of Assam. High conservation priority areas, which make up 5.5% of total area are located in Assam and Arunachal Pradesh. It is important to note that conservation priority areas outside of protected areas pose a greater threat to tigers. We recommend reassessing the IUCN Red List's assigned range map for tigers in the IEHR, as it is over-predicted. Our study has led us to conclude both land use and climate change possess threats to the future habitat of tigers. The outcomes of our study will provide crucial information on identifying habitat hotspots and facilitate appropriate conservation planning efforts.
Collapse
Affiliation(s)
| | | | | | | | | | - Ruchi Badola
- Wildlife Institute of IndiaDehradunUttarakhandIndia
| | | | | |
Collapse
|
4
|
Nittu G, Shameer TT, Nishanthini NK, Sanil R. The tide of tiger poaching in India is rising! An investigation of the intertwined facts with a focus on conservation. GEOJOURNAL 2023; 88:753-766. [PMID: 35431409 PMCID: PMC9005341 DOI: 10.1007/s10708-022-10633-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Accepted: 03/11/2022] [Indexed: 05/09/2023]
Abstract
Poaching and illegal trafficking are major threats to biodiversity, especially when endangered felids are concerned. Tigers are iconic animals, and there is huge demand for their body parts both in the national and international illegal markets. India forms the largest tiger conservation unit in the world and poaching is at its peak even though there are stringent laws and strict enforcement. In the present study, we analytically estimated the tiger seizure cases in India from 2001-2021 using newspaper archives as the main source of data. The data was geo-referenced to understand the details of seizure, demand, and locality. We statistically correlated the seizure rate with the density of tigers, tiger reserves, and various other socio-economic factors. Our result shows that skin, claws, bones, and teeth have more demand, with nails and teeth being the most preferred in local markets. The bones, flesh, and other parts were mostly seized in the border states of the north and eastern states. The intensity of seizures is very high in the states of Maharashtra, Karnataka, Tamil Nadu, and Assam. From our analysis, we predict four trade routes for the export of the seized parts: the Nepal-Bhutan border, Assam border, the Brahmaputra, and the Mumbai port. This corresponds to the five tiger conservation blocks in India, and we observed the seizure rate is high near the Western Ghats region, which has not yet been noticed. Apart from the seizure, we are unconcerned with the seizure's origin or the local trading routes. The study demonstrates the importance of identifying the source population using DNA methods and carefully enforcing the rules in area of poaching. We assert that current approaches are incapable of resolving the issue and that a more precise and effective forensic procedure capable of resolving the issue at the minute local level is critical for precisely tracing trade channels. Supplementary Information The online version contains supplementary material available at 10.1007/s10708-022-10633-4.
Collapse
Affiliation(s)
- George Nittu
- Molecular Biodiversity Lab, Department of Zoology & Wildlife Biology, Government Arts College, Udhagamandalam-643002, The Nilgiris, Tamil Nadu India
| | - Thekke Thumbath Shameer
- Molecular Biodiversity Lab, Department of Zoology & Wildlife Biology, Government Arts College, Udhagamandalam-643002, The Nilgiris, Tamil Nadu India
| | - Nanjanad Kannan Nishanthini
- Molecular Biodiversity Lab, Department of Zoology & Wildlife Biology, Government Arts College, Udhagamandalam-643002, The Nilgiris, Tamil Nadu India
| | - Raveendranathanpillai Sanil
- Molecular Biodiversity Lab, Department of Zoology & Wildlife Biology, Government Arts College, Udhagamandalam-643002, The Nilgiris, Tamil Nadu India
| |
Collapse
|
5
|
Aylward M, Sagar V, Natesh M, Ramakrishnan U. How methodological changes have influenced our understanding of population structure in threatened species: insights from tiger populations across India. Philos Trans R Soc Lond B Biol Sci 2022; 377:20200418. [PMID: 35430878 PMCID: PMC9014192 DOI: 10.1098/rstb.2020.0418] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Accepted: 03/14/2022] [Indexed: 12/26/2022] Open
Abstract
Unprecedented advances in sequencing technology in the past decade allow a better understanding of genetic variation and its partitioning in natural populations. Such inference is critical to conservation: to understand species biology and identify isolated populations. We review empirical population genetics studies of Endangered Bengal tigers within India, where 60-70% of wild tigers live. We assess how changes in marker type and sampling strategy have impacted inferences by reviewing past studies, and presenting three novel analyses including a single-nucleotide polymorphism (SNP) panel, genome-wide SNP markers, and a whole-mitochondrial genome network. At a broad spatial scale, less than 100 SNPs revealed the same patterns of population clustering as whole genomes (with the exception of one additional population sampled only in the SNP panel). Mitochondrial DNA indicates a strong structure between the northeast and other regions. Two studies with more populations sampled revealed further substructure within Central India. Overall, the comparison of studies with varied marker types and sample sets allows more rigorous inference of population structure. Yet sampling of some populations is limited across all studies, and these should be the focus of future sampling efforts. We discuss challenges in our understanding of population structure, and how to further address relevant questions in conservation genetics. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.
Collapse
Affiliation(s)
- Megan Aylward
- National Centre for Biological Sciences, TIFR, Bangalore, India, 560065
| | - Vinay Sagar
- National Centre for Biological Sciences, TIFR, Bangalore, India, 560065
| | - Meghana Natesh
- Indian Institute of Science Education and Research, Tirupati, India, 517507
| | - Uma Ramakrishnan
- National Centre for Biological Sciences, TIFR, Bangalore, India, 560065
- Senior Fellow, DBT Wellcome Trust India Alliance, Hyderabad, Telangana, India, 500034
| |
Collapse
|
6
|
Recapitulating whole genome based population genetic structure for Indian wild tigers through an ancestry informative marker panel. Heredity (Edinb) 2022; 128:88-96. [PMID: 34857925 PMCID: PMC8813985 DOI: 10.1038/s41437-021-00477-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 09/30/2021] [Accepted: 10/01/2021] [Indexed: 02/03/2023] Open
Abstract
Identification of genetic structure within wildlife populations have implications in their conservation and management. Accurately inferring population genetic structure requires whole-genome data across the geographical range of the species, which can be resource-intensive. A cheaper strategy is to employ a subset of markers that can efficiently recapitulate the population genetic structure inferred by the whole genome data. Such ancestry informative markers (AIMs), have rarely been developed for endangered species such as tigers utilizing single nucleotide polymorphisms (SNPs). Here, we first identify the population structure of the Indian tiger using whole-genome sequences and then develop an AIMs panel with a minimum number of SNPs that can recapitulate this structure. We identified four population clusters of Indian tigers with North-East, North-West, and South Indian tigers forming three separate groups, and Terai and Central Indian tigers forming a single cluster. To evaluate the robustness of our AIMs, we applied it to a separate dataset of tigers from across India. Out of 92 SNPs present in our AIMs panel, 49 were present in the new dataset. These 49 SNPs were sufficient to recapitulate the population genetic structure obtained from the whole genome data. To the best of our knowledge, this is the first-ever SNP-based AIMs panel for big cats, which can be used as a cost-effective alternative to whole-genome sequencing for detecting the biogeographical origin of Indian tigers. Our study can be used as a guideline for developing an AIMs panel for the management of other endangered species where obtaining whole genome sequences are difficult.
Collapse
|
7
|
Kannan P, Salaria S, Khan S, Mark T, Baberwal N, Bhatnagar A, Shethia Y, Thatte P, Chanchani P. Assessing Carnivore Occurrence and Community Attitudes Towards Wildlife in a Multi-Use Arid Landscape Corridor. FRONTIERS IN CONSERVATION SCIENCE 2022. [DOI: 10.3389/fcosc.2021.787431] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Small population sizes, low densities, and large area requirements make large carnivores particularly sensitive to habitat degradation and land-use change. In fragmented landscapes, many protected areas cannot accommodate viable wildlife populations in themselves, which brings the surrounding human-dominated matrix that may extend wildlife habitats or serve as corridors into focus. Such areas are typically excluded from the conservation portfolio and are subject to rapid land -use change in many areas. This study investigates the occurrence of tigers, sloth bears, leopards and striped hyenas and assesses community use of natural resources and attitudes towards wildlife in a 3,384 km2 portion of semi-arid multiple-use landscape in Western India that also serves as an important wildlife corridor. This area abuts Ranthambore Tiger Reserve, a preeminent protected area in Western India. Sign surveys spanning 1,039.22 km of trails were conducted in 94, 36 km2 grids spanning agricultural land, forests and other land use types to collate information on wildlife occurrence and associated environmental and human factors. Analysis using occupancy models revealed that tiger and sloth bear occurrence probabilities (0.093 ± 0.05), and (0.13 ± 0.02) were considerably lower than those for leopards (0.72 ± 0.22) and striped hyenas (0.91 ± 0.08). Lack of sufficient cover and limited food availability renders these multiple-use habitats poorly suited for tigers and sloth bears, while leopards and hyenas are able to adapt better to multi-use areas. Concurrently, 66 villages were surveyed across the study landscape, where data on broad socio-economic attributes of communities and their attitudes towards wildlife were assessed through questionnaire surveys. More respondents expressed negative attitudes than positive attitudes which vary as a function of education levels, occupation and land holding sizes. Ongoing landscape transformation through mining, agricultural expansion, infrastructure development, and negative attitudes towards wildlife conservation among people living in the agricultural matrix threatens the long-term functionality of these corridors. Therefore, immediate measures are needed to develop and implement corridor conservation strategies and plans, with a focus on land use planning and human-wildlife conflict mitigation. In the absence of decisive and timely action, wildlife populations may increasingly get relegated to fragmented patches, jeopardising their persistence.
Collapse
|
8
|
Khan A, Patel K, Shukla H, Viswanathan A, van der Valk T, Borthakur U, Nigam P, Zachariah A, Jhala YV, Kardos M, Ramakrishnan U. Genomic evidence for inbreeding depression and purging of deleterious genetic variation in Indian tigers. Proc Natl Acad Sci U S A 2021; 118:e2023018118. [PMID: 34848534 PMCID: PMC8670471 DOI: 10.1073/pnas.2023018118] [Citation(s) in RCA: 80] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/11/2021] [Indexed: 01/03/2023] Open
Abstract
Increasing habitat fragmentation leads to wild populations becoming small, isolated, and threatened by inbreeding depression. However, small populations may be able to purge recessive deleterious alleles as they become expressed in homozygotes, thus reducing inbreeding depression and increasing population viability. We used whole-genome sequences from 57 tigers to estimate individual inbreeding and mutation load in a small-isolated and two large-connected populations in India. As expected, the small-isolated population had substantially higher average genomic inbreeding (FROH = 0.57) than the large-connected (FROH = 0.35 and FROH = 0.46) populations. The small-isolated population had the lowest loss-of-function mutation load, likely due to purging of highly deleterious recessive mutations. The large populations had lower missense mutation loads than the small-isolated population, but were not identical, possibly due to different demographic histories. While the number of the loss-of-function alleles in the small-isolated population was lower, these alleles were at higher frequencies and homozygosity than in the large populations. Together, our data and analyses provide evidence of 1) high mutation load, 2) purging, and 3) the highest predicted inbreeding depression, despite purging, in the small-isolated population. Frequency distributions of damaging and neutral alleles uncover genomic evidence that purifying selection has removed part of the mutation load across Indian tiger populations. These results provide genomic evidence for purifying selection in both small and large populations, but also suggest that the remaining deleterious alleles may have inbreeding-associated fitness costs. We suggest that genetic rescue from sources selected based on genome-wide differentiation could offset any possible impacts of inbreeding depression.
Collapse
Affiliation(s)
- Anubhab Khan
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India;
| | - Kaushalkumar Patel
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India
| | - Harsh Shukla
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India
| | - Ashwin Viswanathan
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India
- Nature Conservation Foundation, Mysore 570017, India
| | | | | | - Parag Nigam
- Wildlife Institute of India, Dehradun 248001, India
| | | | | | - Marty Kardos
- Northwest Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration, Seattle, WA 98112;
| | - Uma Ramakrishnan
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India;
- Department of Biotechnology-Wellcome Trust India Alliance, Hyderabad 500034, India
| |
Collapse
|
9
|
Gallego-García N, Caballero S, Shaffer HB. Are genomic updates of well-studied species worth the investment for conservation? A case study of the Critically Endangered Magdalena River turtle. J Hered 2021; 112:575-589. [PMID: 34628509 DOI: 10.1093/jhered/esab063] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Accepted: 10/08/2021] [Indexed: 11/13/2022] Open
Abstract
As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used RADseq to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment.
Collapse
Affiliation(s)
- Natalia Gallego-García
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA, USA.,Departamento de Ciencias Biológicas, Laboratorio de Ecología Molecular de Vertebrados Acuáticos LEMVA, Universidad de los Andes, Bogotá, Colombia
| | - Susana Caballero
- Departamento de Ciencias Biológicas, Laboratorio de Ecología Molecular de Vertebrados Acuáticos LEMVA, Universidad de los Andes, Bogotá, Colombia
| | - H Bradley Shaffer
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA, USA.,La Kretz Center for California Conservation Science, and Institute of the Environment and Sustainability, University of California, Los Angeles, CA, USA
| |
Collapse
|
10
|
Sagar V, Kaelin CB, Natesh M, Reddy PA, Mohapatra RK, Chhattani H, Thatte P, Vaidyanathan S, Biswas S, Bhatt S, Paul S, Jhala YV, Verma MM, Pandav B, Mondol S, Barsh GS, Swain D, Ramakrishnan U. High frequency of an otherwise rare phenotype in a small and isolated tiger population. Proc Natl Acad Sci U S A 2021; 118:e2025273118. [PMID: 34518374 PMCID: PMC8488692 DOI: 10.1073/pnas.2025273118] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/09/2021] [Indexed: 11/18/2022] Open
Abstract
Most endangered species exist today in small populations, many of which are isolated. Evolution in such populations is largely governed by genetic drift. Empirical evidence for drift affecting striking phenotypes based on substantial genetic data are rare. Approximately 37% of tigers (Panthera tigris) in the Similipal Tiger Reserve (in eastern India) are pseudomelanistic, characterized by wide, merged stripes. Camera trap data across the tiger range revealed the presence of pseudomelanistic tigers only in Similipal. We investigated the genetic basis for pseudomelanism and examined the role of drift in driving this phenotype's frequency. Whole-genome data and pedigree-based association analyses from captive tigers revealed that pseudomelanism cosegregates with a conserved and functionally important coding alteration in Transmembrane Aminopeptidase Q (Taqpep), a gene responsible for similar traits in other felid species. Noninvasive sampling of tigers revealed a high frequency of the Taqpep p.H454Y mutation in Similipal (12 individuals, allele frequency = 0.58) and absence from all other tiger populations (395 individuals). Population genetic analyses confirmed few (minimal number) tigers in Similipal, and its genetic isolation, with poor geneflow. Pairwise FST (0.33) at the mutation site was high but not an outlier. Similipal tigers had low diversity at 81 single nucleotide polymorphisms (mean heterozygosity = 0.28, SD = 0.27). Simulations were consistent with founding events and drift as possible drivers for the observed stark difference of allele frequency. Our results highlight the role of stochastic processes in the evolution of rare phenotypes. We highlight an unusual evolutionary trajectory in a small and isolated population of an endangered species.
Collapse
Affiliation(s)
- Vinay Sagar
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India;
| | - Christopher B Kaelin
- Department of Genetics, Stanford University, Palo Alto, CA 94309
- HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806
| | - Meghana Natesh
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India
- Biology Department, Indian Institute of Science Education and Research, Tirupati 411008, India
| | - P Anuradha Reddy
- Laboratory for Conservation of Endangered Species, Center for Cellular & Molecular Biology, Hyderabad 500048, India
| | | | - Himanshu Chhattani
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India
| | - Prachi Thatte
- World Wide Fund for Nature - India, New Delhi 110003 India
| | - Srinivas Vaidyanathan
- Foundation for Ecological Research, Advocacy and Learning, Auroville Post, Tamil Nadu 605101 India
| | | | | | - Shashi Paul
- Odisha Forest Department, Bhubaneswar 751023, India
| | - Yadavendradev V Jhala
- Wildlife Institute of India, Dehradun 248001, India
- National Tiger Conservation Authority, Wildlife Institute of India Tiger Cell, Wildlife Institute of India, Dehradun 248001, India
| | | | | | | | - Gregory S Barsh
- Department of Genetics, Stanford University, Palo Alto, CA 94309
- HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806
| | - Debabrata Swain
- Former Member Secretary, National Tiger Conservation Authority, New Delhi 110003, India
- Former Principal Chief Conservator of Forest and Head of Forest Force, Indian Forest Service, Bhubaneswar 751023, India
| | - Uma Ramakrishnan
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore 560065, India;
- DBT - Wellcome Trust India Alliance, Hyderabad 500034, India
| |
Collapse
|
11
|
Armstrong EE, Khan A, Taylor RW, Gouy A, Greenbaum G, Thiéry A, Kang JT, Redondo SA, Prost S, Barsh G, Kaelin C, Phalke S, Chugani A, Gilbert M, Miquelle D, Zachariah A, Borthakur U, Reddy A, Louis E, Ryder OA, Jhala YV, Petrov D, Excoffier L, Hadly E, Ramakrishnan U. Recent Evolutionary History of Tigers Highlights Contrasting Roles of Genetic Drift and Selection. Mol Biol Evol 2021; 38:2366-2379. [PMID: 33592092 PMCID: PMC8136513 DOI: 10.1093/molbev/msab032] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Species conservation can be improved by knowledge of evolutionary and genetic history. Tigers are among the most charismatic of endangered species and garner significant conservation attention. However, their evolutionary history and genomic variation remain poorly known, especially for Indian tigers. With 70% of the world’s wild tigers living in India, such knowledge is critical. We re-sequenced 65 individual tiger genomes representing most extant subspecies with a specific focus on tigers from India. As suggested by earlier studies, we found strong genetic differentiation between the putative tiger subspecies. Despite high total genomic diversity in India, individual tigers host longer runs of homozygosity, potentially suggesting recent inbreeding or founding events, possibly due to small and fragmented protected areas. We suggest the impacts of ongoing connectivity loss on inbreeding and persistence of Indian tigers be closely monitored. Surprisingly, demographic models suggest recent divergence (within the last 20,000 years) between subspecies and strong population bottlenecks. Amur tiger genomes revealed the strongest signals of selection related to metabolic adaptation to cold, whereas Sumatran tigers show evidence of weak selection for genes involved in body size regulation. We recommend detailed investigation of local adaptation in Amur and Sumatran tigers prior to initiating genetic rescue.
Collapse
Affiliation(s)
| | - Anubhab Khan
- National Centre for Biological Sciences, TIFR, Bangalore, India
| | - Ryan W Taylor
- Department of Biology, Stanford University, Stanford, CA, USA.,End2End Genomics, LLC, Davis, CA, USA
| | - Alexandre Gouy
- Institute of Ecology and Evolution, University of Bern, Bern, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Gili Greenbaum
- Department of Biology, Stanford University, Stanford, CA, USA.,Department of Ecology, Evolution & Behavior, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Alexandre Thiéry
- Institute of Ecology and Evolution, University of Bern, Bern, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Jonathan T Kang
- Department of Biology, Stanford University, Stanford, CA, USA.,Genome Institute of Singapore, A*STAR, Singapore
| | | | - Stefan Prost
- Department of Biology, Stanford University, Stanford, CA, USA
| | - Gregory Barsh
- Hudsonalpha Institute, Hunstville, AL, USA.,Department of Genetics, Stanford University, Stanford, CA, USA
| | | | | | | | - Martin Gilbert
- Wildlife Conservation Society, Russia Program, New York, NY, USA.,College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Dale Miquelle
- Wildlife Conservation Society, Russia Program, New York, NY, USA
| | | | | | - Anuradha Reddy
- Laboratory for Conservation of Endangered Species, CCMB, Hyderabad, India
| | - Edward Louis
- Department of Genetics, Omaha Zoo, Omaha, NE, USA
| | - Oliver A Ryder
- San Diego Zoo, Institute for Conservation Research, Escondido, CA, USA
| | | | - Dmitri Petrov
- Department of Biology, Stanford University, Stanford, CA, USA
| | - Laurent Excoffier
- Institute of Ecology and Evolution, University of Bern, Bern, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Elizabeth Hadly
- Wildlife Conservation Society, Russia Program, New York, NY, USA
| | | |
Collapse
|
12
|
Samaha G, Wade CM, Mazrier H, Grueber CE, Haase B. Exploiting genomic synteny in Felidae: cross-species genome alignments and SNV discovery can aid conservation management. BMC Genomics 2021; 22:601. [PMID: 34362297 PMCID: PMC8348863 DOI: 10.1186/s12864-021-07899-2] [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: 09/30/2020] [Accepted: 07/14/2021] [Indexed: 11/10/2022] Open
Abstract
Background While recent advances in genomics has enabled vast improvements in the quantification of genome-wide diversity and the identification of adaptive and deleterious alleles in model species, wildlife and non-model species have largely not reaped the same benefits. This has been attributed to the resources and infrastructure required to develop essential genomic datasets such as reference genomes. In the absence of a high-quality reference genome, cross-species alignments can provide reliable, cost-effective methods for single nucleotide variant (SNV) discovery. Here, we demonstrated the utility of cross-species genome alignment methods in gaining insights into population structure and functional genomic features in cheetah (Acinonyx jubatas), snow leopard (Panthera uncia) and Sumatran tiger (Panthera tigris sumatrae), relative to the domestic cat (Felis catus). Results Alignment of big cats to the domestic cat reference assembly yielded nearly complete sequence coverage of the reference genome. From this, 38,839,061 variants in cheetah, 15,504,143 in snow leopard and 13,414,953 in Sumatran tiger were discovered and annotated. This method was able to delineate population structure but limited in its ability to adequately detect rare variants. Enrichment analysis of fixed and species-specific SNVs revealed insights into adaptive traits, evolutionary history and the pathogenesis of heritable diseases. Conclusions The high degree of synteny among felid genomes enabled the successful application of the domestic cat reference in high-quality SNV detection. The datasets presented here provide a useful resource for future studies into population dynamics, evolutionary history and genetic and disease management of big cats. This cross-species method of variant discovery provides genomic context for identifying annotated gene regions essential to understanding adaptive and deleterious variants that can improve conservation outcomes. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07899-2.
Collapse
Affiliation(s)
- Georgina Samaha
- Sydney School of Veterinary Science, Faculty of Science, The University of Sydney, Sydney, NSW, Australia.
| | - Claire M Wade
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| | - Hamutal Mazrier
- Sydney School of Veterinary Science, Faculty of Science, The University of Sydney, Sydney, NSW, Australia
| | - Catherine E Grueber
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| | - Bianca Haase
- Sydney School of Veterinary Science, Faculty of Science, The University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
13
|
Considerations for Initiating a Wildlife Genomics Research Project in South and South-East Asia. J Indian Inst Sci 2021. [DOI: 10.1007/s41745-021-00243-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
14
|
Trends in Wildlife Connectivity Science from the Biodiverse and Human-Dominated South Asia. J Indian Inst Sci 2021. [DOI: 10.1007/s41745-021-00240-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
15
|
Snow Leopard (Panthera uncia) Genetics: The Knowledge Gaps, Needs, and Implications for Conservation. J Indian Inst Sci 2021. [DOI: 10.1007/s41745-021-00236-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
16
|
Gaitán-Espitia JD, Hobday AJ. Evolutionary principles and genetic considerations for guiding conservation interventions under climate change. GLOBAL CHANGE BIOLOGY 2021; 27:475-488. [PMID: 32979891 DOI: 10.1111/gcb.15359] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Accepted: 08/26/2020] [Indexed: 05/25/2023]
Abstract
Impacts of climate change are apparent in natural systems around the world. Many species are and will continue to struggle to persist in their current location as their preferred environment changes. Traditional conservation efforts aiming to prevent local extinctions have focused on two aspects that theoretically enhance genetic diversity-population connectivity and population size-through 'passive interventions' (such as protected areas and connectivity corridors). However, the exceptionally rapid loss of biodiversity that we are experiencing as result of anthropogenic climate change has shifted conservation approaches to more 'active interventions' (such as rewilding, assisted gene flow, assisted evolution, artificial selection, genetic engineering). We integrate genetic/genomic approaches into an evolutionary biology framework in order to discuss with scientists, conservation managers and decision makers about the opportunities and risks of interventions that need careful consideration in order to avoid unwanted evolutionary outcomes.
Collapse
Affiliation(s)
- Juan D Gaitán-Espitia
- CSIRO Oceans and Atmosphere, Hobart, Tas., Australia
- The Swire Institute of Marine Science and School of Biological Sciences, The University of Hong Kong, Hong Kong SAR, China
| | | |
Collapse
|
17
|
Bit A, Thakur M, Singh SK, Joshi BD, Singh VK, Sharma LK, Tripathy B, Chandra K. Assembling mitogenome of Himalayan Black Bear (U. t. laniger) from low depth reads and its application in drawing phylogenetic inferences. Sci Rep 2021; 11:730. [PMID: 33436634 PMCID: PMC7803731 DOI: 10.1038/s41598-020-76872-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Accepted: 11/04/2020] [Indexed: 12/01/2022] Open
Abstract
The complete mitogenome of Himalayan black bear (Ursus thibetanus laniger) from Indian Himalayan region was assembled following the modified approach of mitochondrial baiting and mapping using the next-generation sequencing reads. The complete mitogenome was of 16,556 bp long, consisted of 37 genes that contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 control region. The complete base composition was 31.33% A, 15.24% G, 25.45%C, and 27.98%T and gene arrangement was similar to the other sub-species of Asiatic black bear. The relative synonymous codon usage analysis revealed the maximum abundance of Isoleucine, Tyrosine, Leucine and Threonine. The assembled mitogenome of U. t. laniger exhibited 99% similarity with the mitogenomes of Himalayan black bear available from Nepal and Tibetan Plateau-Himalaya region. The findings of the present study has proven low depth sequencing data, adequate and highly efficient in rapid recovering the mitochondrial genome by overcoming the conventional strategies of obtaining long-range PCR and subsequently drawing phylogenetic inferences.
Collapse
Affiliation(s)
- Amrita Bit
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| | - Mukesh Thakur
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India.
| | - Sujeet Kumar Singh
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| | - Bheem Dutt Joshi
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| | - Vinay Kumar Singh
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| | - Lalit Kumar Sharma
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| | - Basudev Tripathy
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| | - Kailash Chandra
- Zoological Survey of India, New Alipore, Kolkata, West Bengal, 700053, India
| |
Collapse
|
18
|
Khan A, Patel K, Bhattacharjee S, Sharma S, Chugani AN, Sivaraman K, Hosawad V, Sahu YK, Reddy GV, Ramakrishnan U. Are shed hair genomes the most effective noninvasive resource for estimating relationships in the wild? Ecol Evol 2020; 10:4583-4594. [PMID: 32551045 PMCID: PMC7297754 DOI: 10.1002/ece3.6157] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2019] [Revised: 02/10/2020] [Accepted: 02/17/2020] [Indexed: 01/01/2023] Open
Abstract
Knowledge of relationships in wild populations is critical for better understanding mating systems and inbreeding scenarios to inform conservation strategies for endangered species. To delineate pedigrees in wild populations, study genetic connectivity, study genotype-phenotype associations, trace individuals, or track wildlife trade, many identified individuals need to be genotyped at thousands of loci, mostly from noninvasive samples. This requires us to (a) identify the most common noninvasive sample available from identified individuals, (b) assess the ability to acquire genome-wide data from such samples, and (c) evaluate the quality of such genome-wide data, and its ability to reconstruct relationships between animals within a population.We followed identified individuals from a wild endangered tiger population and found that shed hair samples were the most common compared to scat samples, opportunistically found carcasses, and opportunistic invasive samples. We extracted DNA from these samples, prepared whole genome sequencing libraries, and sequenced genomes from these.Whole genome sequencing methods resulted in between 25%-98% of the genome sequenced for five such samples. Exploratory population genetic analyses revealed that these data were free of holistic biases and could recover expected population structure and relatedness. Mitochondrial genomes recovered matrilineages in accordance with long-term monitoring data. Even with just five samples, we were able to uncover the matrilineage for three individuals with unknown ancestry.In summary, we demonstrated that noninvasive shed hair samples yield adequate quality and quantity of DNA in conjunction with sensitive library preparation methods, and provide reliable data from hundreds of thousands of SNPs across the genome. This makes shed hair an ideal noninvasive resource for studying individual-based genetics of elusive endangered species in the wild.
Collapse
Affiliation(s)
- Anubhab Khan
- National Centre for Biological SciencesTIFRBangaloreIndia
| | | | - Subhadeep Bhattacharjee
- Rajasthan Forest DepartmentJaipurIndia
- Department of GEMESUniversity of JohannesburgJohannesburgSouth Africa
| | | | | | | | | | | | | | | |
Collapse
|
19
|
Hebbar P, Ravikanth G, Aravind NA. A review on the conservation genetic studies of Indian amphibians and their implications on developing strategies for conservation†. J Genet 2019. [DOI: 10.1007/s12041-019-1159-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
20
|
Hebbar P, Ravikanth G, Aravind NA. A review on the conservation genetic studies of Indian amphibians and their implications on developing strategies for conservation. J Genet 2019; 98:114. [PMID: 31819027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Amphibians show a very high level of diversity and endemism and are facing global declines from the past few decades. Studies have shown that the molecular tools can be helpful in their conservation efforts. In India, more than 80% of amphibians are endemic and most show a narrow range of distribution. Most of the Indian amphibians lack information on their genetic diversity. In this study, were view the overall trend on amphibian studies in India with the specific focus on conservation genetics. Overall, of the 173 studies, only 14 dealt with the conservation of amphibians through genetic tools and five studies estimated the genetic diversity or gene structure. Here, we discuss the gaps and provide future directions on how genetic studies can be helpful in Indian amphibian conservation.
Collapse
Affiliation(s)
- Priti Hebbar
- Suri Sehgal Centre for Biodiversity and Conservation, Ashoka Trust for Research in Ecology and the Environment (ATREE), Bengaluru 560 064, India.
| | | | | |
Collapse
|
21
|
Alcala N, Goldberg A, Ramakrishnan U, Rosenberg NA. Coalescent Theory of Migration Network Motifs. Mol Biol Evol 2019; 36:2358-2374. [PMID: 31165149 PMCID: PMC6759081 DOI: 10.1093/molbev/msz136] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Natural populations display a variety of spatial arrangements, each potentially with a distinctive impact on genetic diversity and genetic differentiation among subpopulations. Although the spatial arrangement of populations can lead to intricate migration networks, theoretical developments have focused mainly on a small subset of such networks, emphasizing the island-migration and stepping-stone models. In this study, we investigate all small network motifs: the set of all possible migration networks among populations subdivided into at most four subpopulations. For each motif, we use coalescent theory to derive expectations for three quantities that describe genetic variation: nucleotide diversity, FST, and half-time to equilibrium diversity. We describe the impact of network properties on these quantities, finding that motifs with a high mean node degree have the largest nucleotide diversity and the longest time to equilibrium, whereas motifs with low density have the largest FST. In addition, we show that the motifs whose pattern of variation is most strongly influenced by loss of a connection or a subpopulation are those that can be split easily into disconnected components. We illustrate our results using two example data sets—sky island birds of genus Sholicola and Indian tigers—identifying disturbance scenarios that produce the greatest reduction in genetic diversity; for tigers, we also compare the benefits of two assisted gene flow scenarios. Our results have consequences for understanding the effect of geography on genetic diversity, and they can assist in designing strategies to alter population migration networks toward maximizing genetic variation in the context of conservation of endangered species.
Collapse
Affiliation(s)
- Nicolas Alcala
- Department of Biology, Stanford University, Stanford, CA
| | - Amy Goldberg
- Department of Biology, Stanford University, Stanford, CA.,Department of Evolutionary Anthropology, Duke University, Durham, NC
| | - Uma Ramakrishnan
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India
| | | |
Collapse
|
22
|
Natesh M, Taylor RW, Truelove NK, Hadly EA, Palumbi SR, Petrov DA, Ramakrishnan U. Empowering conservation practice with efficient and economical genotyping from poor quality samples. Methods Ecol Evol 2019; 10:853-859. [PMID: 31511786 PMCID: PMC6738957 DOI: 10.1111/2041-210x.13173] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2018] [Accepted: 02/21/2019] [Indexed: 01/09/2023]
Abstract
Moderate- to high-density genotyping (100 + SNPs) is widely used to determine and measure individual identity, relatedness, fitness, population structure and migration in wild populations.However, these important tools are difficult to apply when high-quality genetic material is unavailable. Most genomic tools are developed for high-quality DNA sources from laboratory or medical settings. As a result, most genetic data from market or field settings is limited to easily amplified mitochondrial DNA or a few microsatellites.To enable genotyping in conservation contexts, we used next-generation sequencing of multiplex PCR products from very low-quality DNA extracted from faeces, hair and cooked samples. We demonstrated utility and wide-ranging potential application in endangered wild tigers and tracking commercial trade in Caribbean queen conch.We genotyped 100 SNPs from degraded tiger samples to identify individuals, discern close relatives and detect population differentiation. Co-occurring carnivores do not amplify (e.g. Indian wild dog/dhole) or are monomorphic (e.g. leopard). Sixty-two SNPs from conch fritters and field-collected samples were used to test relatedness and detect population structure.We provide proof of concept for a rapid, simple, cost-effective and scalable method (for both samples and number of loci), a framework that can be applied to other conservation scenarios previously limited by low-quality DNA samples. These approaches provide a critical advance for wildlife monitoring and forensics, open the door to field-ready testing, and will strengthen the use of science in policy decisions and wildlife trade.
Collapse
Affiliation(s)
- Meghana Natesh
- National Centre for Biological Sciences, TIFR, Bangalore, India
- Sastra University, Tirumalaisamudram, Thanjavur, India
| | - Ryan W. Taylor
- Department of Biology, Stanford University, Stanford, California
- End2End Genomics LLC, Davis, California
| | | | | | - Stephen R. Palumbi
- Department of Biology, Stanford University, Stanford, California
- Hopkins Marine Station, Stanford University, Pacific Grove, California
| | - Dmitri A. Petrov
- Department of Biology, Stanford University, Stanford, California
| | | |
Collapse
|