Comparative genomics uncovers the evolutionary history, demography, and molecular adaptations of South American canids

Computationally efficient demographic history inference from allele frequencies with supervised machine learning

Inferring past demographic history of natural populations from genomic data is of central concern in many studies across research fields. Previously, our group had developed dadi, a widely used demographic history inference method based on the allele frequency spectrum (AFS) and maximum composite likelihood optimization. However, dadi's optimization procedure can be computationally expensive. Here, we developed donni (demography optimization via neural network inference), a new inference method based on dadi that is more efficient while maintaining comparable inference accuracy. For each dadi-supported demographic model, donni simulates the expected AFS for a range of model parameters then trains a set of Mean Variance Estimation neural networks using the simulated AFS. Trained networks can then be used to instantaneously infer the model parameters from future input data AFS. We demonstrated that for many demographic models, donni can infer some parameters, such as population size changes, very well and other parameters, such as migration rates and times of demographic events, fairly well. Importantly, donni provides both parameter and confidence interval estimates from input AFS with accuracy comparable to parameters inferred by dadi's likelihood optimization while bypassing its long and computationally intensive evaluation process. donni's performance demonstrates that supervised machine learning algorithms may be a promising avenue for developing more sustainable and computationally efficient demographic history inference methods.

Complex Evolutionary History With Extensive Ancestral Gene Flow in an African Primate Radiation

Understanding the drivers of speciation is fundamental in evolutionary biology, and recent studies highlight hybridization as an important evolutionary force. Using whole-genome sequencing data from 22 species of guenons (tribe Cercopithecini), one of the world's largest primate radiations, we show that rampant gene flow characterizes their evolutionary history and identify ancient hybridization across deeply divergent lineages that differ in ecology, morphology, and karyotypes. Some hybridization events resulted in mitochondrial introgression between distant lineages, likely facilitated by cointrogression of coadapted nuclear variants. Although the genomic landscapes of introgression were largely lineage specific, we found that genes with immune functions were overrepresented in introgressing regions, in line with adaptive introgression, whereas genes involved in pigmentation and morphology may contribute to reproductive isolation. In line with reports from other systems that hybridization might facilitate diversification, we find that some of the most species-rich guenon clades are of admixed origin. This study provides important insights into the prevalence, role, and outcomes of ancestral hybridization in a large mammalian radiation.

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

Genomic Basis of Adaptive Divergence in Leg Length between Ground- and Tree-Dwelling Species within a Bird Family

Hind limbs of tetrapods vary greatly in length and the variability can be associated with locomotor adaptation. Although the phenotypic evolution has been well documented, the underlying genetic basis remains poorly understood. We address this issue by integrating comparative genomics and functional prediction with a study system consisting of ground-dwelling, long-legged and tree-dwelling, short-legged species within the avian family Paridae. Genome-wide divergence and phenotypic correlation analyses jointly identified five highly divergent genomic regions that are significantly related with the difference in leg length between these two groups. Gene annotation for these regions detected three genes involved in skeletal development, that is, PTPA, BRINP1, and MIGA2, with the first one being under the strongest selection. Furthermore, four single nucleotide polymorphisms (SNPs) in the coding region of PTPA can well distinguish the two groups with distinct leg length. Among the four SNPs, one is non-synonymous mutation, and according to the prediction for protein structure and function, it can affect the 3D structure of the encoded protein by altering the corresponding amino acid's position. The alleles of PTPA were found in all sequenced species of the orders Palaeognathae and Psittaciformes, which typically take a ground locomotion style. A whole-genome scanning across bird species uncovered that the four SNPs are more likely to be present in resident passerines with increased leg length/wing length ratios (a proxy of leg-dependent locomotion efficiency). Our findings provide insight into the molecular evolution of locomotion performance based on leg morphology in birds.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Different, but the same: Inferring the hunting behaviour of the hypercarnivorous bush dog (Speothos venaticus) through finite element analysis

Cerdocyonina is a clade composed by the South-American canids in which the bush dog (Speothos venaticus) is one of the most elusive species. Known for its unique morphology within the group, this small, bear-like faced canid is the only member of the clade adapted to hypercarnivory, an almost exclusively meat-based diet currently present only in usually large, pack-hunting canids such as the grey wolf (Canis lupus). However, much of the biology of the bush dog is poorly understood, and inferences about its ecology, hunting strategies and diet are usually based on observation of captive individuals and anecdotal records, with reduced quantitative data to offer support. Here, we investigated the craniomandibular functional morphology of the bush dog through finite element analysis (FEA). FEA was employed to model the biting behaviour and to create extrinsic and intrinsic functional scenarios with different loads, corresponding to different bites used to subdue and process the prey. For comparison, the same modelling was applied to the skull of a grey wolf and a grey fox (Urocyon cinereoargenteus). Our analysis showed that the bush dog's responses to loading are more similar to the wolf's than to the fox's in most scenarios, suggesting a convergent craniomandibular functional morphology between these two hypercarnivorous species, despite their distinct phylogenetic positions and body sizes. Differences between the three taxa are noteworthy and suggested to be related to the size of the usual prey. The modelled bite force for the bush dog is relatively strong, about half of that estimated for the wolf and about 40% stronger than the fox's bite. The results strengthen with quantitative data the inferences of the bush dog as a pack-hunting predator with prey size similar to its own, such as large rodents and armadillos, being specialised in subduing and killing its prey using multiple bites. Its similarity to the wolf also confirms anecdotal accounts of predation on mammals that are much larger than itself, such as peccaries and tapirs. These data highlight the ecological specialisation of this small canid in a continent where large, pack-hunting canids are absent.