skip to main content


Title: Ancient divergence of Indian and Tibetan wolves revealed by recombination‐aware phylogenomics
Abstract

The grey wolf (Canis lupus) expanded its range across Holarctic regions during the late Pleistocene. Consequently, most grey wolves share recent (<100,000 years ago) maternal origins corresponding to a widespread Holarctic clade. However, two deeply divergent (200,000–700,000 years ago) mitochondrial clades are restricted, respectively, to the Indian subcontinent and the Tibetan Plateau, where remaining wolves are endangered. No genome‐wide analysis had previously included wolves corresponding to the mitochondrial Indian clade or attempted to parse gene flow and phylogeny. We sequenced four Indian and two Tibetan wolves and included 31 additional canid genomes to resolve the phylogenomic history of grey wolves. Genomic analyses revealed Indian and Tibetan wolves to be distinct from each other and from broadly distributed wolf populations corresponding to the mitochondrial Holarctic clade. Despite gene flow, which was reflected disproportionately in high‐recombination regions of the genome, analyses revealed Indian and Tibetan wolves to be basal to Holarctic grey wolves, in agreement with the mitochondrial phylogeny. In contrast to mitochondrial DNA, however, genomic findings suggest the possibility that the Indian wolf could be basal to the Tibetan wolf, a discordance potentially reflecting selection on the mitochondrial genome. Together, these findings imply that southern regions of Asia have been important centers for grey wolf evolution and that Indian and Tibetan wolves represent evolutionary significant units (ESUs). Further study is needed to assess whether these ESUs warrant recognition as distinct species. This question is especially urgent regarding the Indian wolf, which represents one of the world's most endangered wolf populations.

 
more » « less
NSF-PAR ID:
10377308
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
30
Issue:
24
ISSN:
0962-1083
Page Range / eLocation ID:
p. 6687-6700
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Among the three main divergent lineages of gray wolf (Canis lupus), the Holarctic lineage is the most widespread and best studied, particularly in North America and Europe. Less is known about Tibetan (also called Himalayan) and Indian wolf lineages in southern Asia, especially in areas surrounding Pakistan where all three lineages are thought to meet. Given the endangered status of the Indian wolf in neighboring India and unclear southwestern boundary of the Tibetan wolf range, we conducted mitochondrial and genome-wide sequencing of wolves from Pakistan and Kyrgyzstan. Sequences of the mitochondrial D-loop region of 81 wolves from Pakistan indicated contact zones between Holarctic and Indian lineages across the northern and western mountains of Pakistan. Reduced-representation genome sequencing of eight wolves indicated an east-to-west cline of Indian to Holarctic ancestry, consistent with a contact zone between these two lineages in Pakistan. The western boundary of the Tibetan lineage corresponded to the Ladakh region of India’s Himalayas with a narrow zone of admixture spanning this boundary from the Karakoram Mountains of northern Pakistan into Ladakh, India. Our results highlight the conservation significance of Pakistan’s wolf populations, especially the remaining populations in Sindh and Southern Punjab that represent the highly endangered Indian lineage.

     
    more » « less
  2. Abstract

    The phylogeny ofCyclops(~30 spp.), a predominantly Palearctic cold‐adapted genus, was reconstructed based on morphological and molecular characters. The morphological analysis used extensive taxon sampling from the entire Holarctic range of the genus and included 53 morphological characters. Polymorphic traits were coded by the “unordered,” “unscaled” and “scaled” methods; maximum parsimony criterion was applied in tree building. Molecular phylogenetic reconstructions utilized partial nuclear 18S and 28S ribosomal genes, mitochondrial cytochrome oxidase I and complete internal transcribed spacer regions I and II, albeit with limited taxon sampling. Bayesian inference and maximum likelihood were used in these tree reconstructions. The molecular characters were used both in combination with morphology and as an independent test of the basal relationships inferred from morphology. Monophyly of the genus received strong support in both the morphological and molecular phylogenies; the basal relationships remain unresolved. The morphology‐based phylogenies, along with the geographic distribution patterns and ecological traits, supported monophyly of theankyraeladakanusclade,scutifer‐clade (C. scutifer,C. jashnovi,C. columbianus),kolensis‐clade (C. kolensis,C. kikuchii,C. vicinus,C. furcifer,C. insignis,C. alaskaensis),abyssorum‐clade (C. abyssorums. str.,C. abyssorum larianus,C. ricae,C. sevani) anddivergens‐clade (South Carpathian “Cyclopssp. Y,”C. mauritaniae,C. divergens,C. bohater,C. lacustris). Relationships among European and North American populations ofC. scutiferandC. columbianusbased on partial sequences of the 12S mitochondrial gene showC. scutiferto be paraphyletic, suggesting two independent invasions into North America via the Bering Land Bridge from Siberia to Alaska.

     
    more » « less
  3. Abstract

    Hybrid zones typically contain novel gene combinations that can be tested by natural selection in a unique genetic context. Parental haplotypes that increase fitness can introgress beyond the hybrid zone, into the range of parental species. We used the Affymetrix canineSNPgenotyping array to identify genomic regions tagged by multiple ancestry informative markers that are more frequent in an admixed population than expected. We surveyed a hybrid zone formed in the last 100 years as coyotes expanded their range into eastern North America. Concomitant with expansion, coyotes hybridized with wolves and some populations became more wolflike, such that coyotes in the northeast have the largest body size of any coyote population. Using a set of 3102 ancestry informative markers, we identified 60 differentially introgressed regions in 44 canines across this admixture zone. These regions are characterized by an excess of exogenous ancestry and, in northeastern coyotes, are enriched for genes affecting body size and skeletal proportions. Further, introgressed wolf‐derived alleles have penetrated into SouthernUScoyote populations. Because no wolves currently exist in this area, these alleles are unlikely to have originated from recent hybridization. Instead, they probably originated from intraspecific gene flow or ancient admixture. We show that grey wolf and coyote admixture has far‐reaching effects and, in addition to phenotypically transforming admixed populations, allows for the differential movement of alleles from different parental species to be tested in new genomic backgrounds.

     
    more » « less
  4. Abstract Context

    Processes that shape genomic and ecological divergence can reveal important evolutionary dynamics to inform the conservation of threatened species.Fontaineais a genus of rainforest shrubs and small trees including critically endangered and threatened species restricted to narrow, but complex geographic and ecological regions. Several species ofFontaineaare subject to spatially explicit conditions and experience limited intra-specific gene flow, likely generating genetic differentiation and local adaptation.

    Objectives

    Here, we explored the genetic and ecological mechanisms underlying patterns of diversification in two, closely related threatenedFontaineaspecies. Our aim was to compare spatial patterns of genetic variation between the vulnerableFontainea australis(Southern Fontainea) and critically endangeredF. oraria(Coastal Fontainea), endemic to the heterogeneous subtropical region of central, eastern Australia, where large-scale clearing has severely reduced rainforest habitat to a fraction (< 1%) of its pre-European settlement extent.

    Methods

    We used a set of 10,000 reduced-representation markers to infer genetic relationships and the drivers of spatial genetic variation across the two species. In addition, we employed a combination of univariate and multivariate genome-environment association analysis using a set of topo-climatic variables to explore potential patterns of local adaptation as a factor impacting genomic divergence.

    Results

    Our study revealed that Coastal Fontainea have a close genetic relationship with Southern Fontainea. We showed that isolation by distance has played a key role in their genetic variation, indicating that vicariance can explain the spatial genetic distribution of the two species. Genotype-environment analyses showed a strong association with temperature and topographic features, suggesting adaptation to localised thermal environments. We used a multivariate redundancy analysis to identify a range of putatively adapted loci associated with local environmental conditions.

    Conclusions

    Divergent selection at the local-habitat scale as a result of dispersal limitations and environmental heterogeneity (including physical barriers) are likely contributors to adaptive divergence between the twoFontaineaspecies. Our findings have presented evidence to indicate that Southern and Coastal Fontainea were comprised of distinct genetic groups and ecotypes, that together may form a single species continuum, with further phenotype research suggested to confirm the current species boundaries. Proactive conservation actions, including assisted migration to enhance the resilience of populations lacking stress-tolerant single nucleotide polymorphisms (SNPs) may be required to secure the long-term future of both taxa. This is especially vital for the critically endangered Coastal Fontainea given projections of habitat decline for the species under future climate scenarios.

     
    more » « less
  5. Abstract

    Intraspecific polymorphism in birds, especially plumage colour polymorphism, and the mechanisms that control it are an area of active research in evolutionary biology. The black‐headed bulbul (Brachypodius atriceps) is a polymorphic species with two distinct morphs, yellow and grey. This species inhabits the mainland and virtually all continental islands of Southeast Asia where yellow morphs predominate, but on two islands in the Sunda region, Bawean and Maratua, grey morphs are common or exclusive. Here, we generated a high‐quality reference genome of a yellow individual and resequenced genomes of multiple individuals of both yellow and grey morphs to study the genetic basis of coloration and population history of the species. Using PCA and STRUCTURE analysis, we found the Maratua Island population (which is exclusively grey) to be distinct from all otherB.atricepspopulations, having been isolated c. 1.9 million years ago (Ma). In contrast, Bawean grey individuals (a subset of yellow and grey individuals on that island) are embedded within an almost panmictic Sundaic clade of yellow birds. UsingFSTanddxyto compare variable genomic segments between Maratua and yellow individuals, we located peaks of divergence and identified candidate loci involved in the colour polymorphism. Tests of selection among coding‐proteins in highFSTregions, however, did not indicate selection on the candidate genes. Overall, we report on some loci that are potentially responsible for the grey/yellow polymorphism in a species that otherwise shows little genetic diversification across most of its range.

     
    more » « less