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

The remarkable radiation of South American (SA) canids produced 10 extant species distributed across diverse habitats, including disparate forms such as the short-legged, hypercarnivorous bush dog and the long-legged, largely frugivorous maned wolf. Despite considerable research spanning nearly two centuries, many aspects of their evolutionary history remain unknown. Here, we analyzed 31 whole genomes encompassing all extant SA canid species to assess phylogenetic relationships, interspecific hybridization, historical demography, current genetic diversity, and the molecular bases of adaptations in the bush dog and maned wolf. We found that SA canids originated from a single ancestor that colonized South America 3.9 to 3.5 Mya, followed by diversification east of the Andes and then a single colonization event and radiation of Lycalopex species west of the Andes. We detected extensive historical gene flow between recently diverged lineages and observed distinct patterns of genomic diversity and demographic history in SA canids, likely induced by past climatic cycles compounded by human-induced population declines. Genome-wide scans of selection showed that disparate limb proportions in the bush dog and maned wolf may derive from mutations in genes regulating chondrocyte proliferation and enlargement. Further, frugivory in the maned wolf may have been enabled by variants in genes associated with energy intake from short-chain fatty acids. In contrast, unique genetic variants detected in the bush dog may underlie interdigital webbing and dental adaptations for hypercarnivory. Our analyses shed light on the evolution of a unique carnivoran radiation and how it was shaped by South American topography and climate change.

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.

Phylogeographic analyses of an endemic Neotropical fox (Lycalopex vetulus) reveal evidence of hybridization with a different canid species (L. gymnocercus).

The hoary fox (Lycalopex vetulus) is the only species of the Canidae (Mammalia: Carnivora) endemic to Brazil, and so far has been the target of few genetic studies. Using microsatellites and mtDNA markers, we investigated its present genetic diversity and population structure. We also tested the hypothesis that this species currently hybridizes with the pampas fox (L. gymnocercus), as suggested by previous mtDNA data from two individuals. We collected tissue and blood samples from animals representing most of the two species' distributions in Brazil (n = 87), including their recently discovered geographic contact zone in São Paulo state. We observed that the hoary fox exhibits high levels of genetic diversity and low levels of population structure. We identified six individuals from São Paulo state with clear evidence of hybridization based on introgressed pampas fox mtDNA and/or admixed microsatellite genotypes (three individuals bore both types of evidence). These results demonstrate the existence of admixed individuals between hoary and pampas foxes in southeastern Brazil, representing the first identified case of inter-species admixture between native South American canids. We discuss our findings in the context of the evolutionary history of these foxes and address potential conservation implications of this interspecies hybridization process.

Genomic Signatures of Divergent Ecological Strategies in a Recent Radiation of Neotropical Wild Cats.

Ecological differentiation among diverging species is an important component of the evolutionary process and can be investigated in rapid and recent radiations. Here, we use whole genome sequences of five species from the genus Leopardus, a recently diversified Neotropical lineage with species bearing distinctive morphological, ecological, and behavioral features, to investigate genome-wide diversity, comparative demographic history and signatures of positive selection. Our results show that divergent ecological strategies are reflected in genomic features, for example a generalist species shows historically larger effective population size and higher heterozygosity than habitat specialists. The demographic history of these cats seems to have been jointly driven by climate fluctuations and habitat specialization, with different ecological adaptations leading to distinct trajectories. Finally, a gene involved in vertebrate retinal neurogenesis (POU4F2) was found to be under positive selection in the margay, a cat with notoriously large eyes that are likely associated with its nocturnal and arboreal specializations.

Structure and Dynamics of Periphyton in a Neotropical Freshwater Lake, with Emphasis on Ciliates and Their Relationships with Bacterial Taxa.

Periphyton communities in freshwater systems play an essential role in biogeochemical processes, but knowledge of their structure and dynamics lags far behind other environments. We used eDNA metabarcoding of 16S and 18S rRNA markers to investigate the formation and establishment of a periphytic community, in addition to a morphology-based approach for peritrich ciliate determinations, its most abundant group. We sampled two nearby sites within a large Neotropical lake at four time points, aiming to assess whether periphyton establishment can be replicated on this local scale. Producers and denitrifiers were abundant in the community, illustrating the relevant role of biofilms in freshwater nutrient recycling. Among microeukaryotes, peritrich ciliates dominated the community, with genera Epistylis and Vorticella being the most abundant and showing a clear succession at both sites. Other ciliates were morphologically identified and, in some cases, their occurrence was strongly related to bacterial abundance. The structure of both prokaryotic and eukaryotic components of periphyton was not different, while the turnover dynamics differed between the two sites, in spite of their adjacent locations and similar abiotic properties. This indicates that the establishment of these communities can vary even on a local scale within a lake ecosystem.

Mitogenomic analysis of a late Pleistocene jaguar from North America.

The jaguar (Panthera onca) is the largest living cat species native to the Americas and one of few large American carnivorans to have survived into the Holocene. However, the extent to which jaguar diversity declined during the end-Pleistocene extinction event remains unclear. For example, Pleistocene jaguar fossils from North America are notably larger than the average extant jaguar, leading to hypotheses that jaguars from this continent represent a now-extinct subspecies (Panthera onca augusta) or species (Panthera augusta). Here, we used a hybridization capture approach to recover an ancient mitochondrial genome from a large, late Pleistocene jaguar from Kingston Saltpeter Cave, Georgia, USA, which we sequenced to 26-fold coverage. We then estimated the evolutionary relationship between the ancient jaguar mitogenome and those from other extinct and living large felids, including multiple jaguars sampled across the species' current range. The ancient mitogenome falls within the diversity of living jaguars. All sampled jaguar mitogenomes share a common mitochondrial ancestor ~400 thousand years ago, indicating that the lineage represented by the ancient specimen dispersed into North America from the south at least once during the late Pleistocene. While genomic data from additional and older specimens will continue to improve understanding of Pleistocene jaguar diversity in the Americas, our results suggest that this specimen falls within the variation of extant jaguars despite the relatively larger size and geographic location and does not represent a distinct taxon.

Genome-wide analyses reveal drivers of penguin diversification.

Penguins are the only extant family of flightless diving birds. They currently comprise at least 18 species, distributed from polar to tropical environments in the Southern Hemisphere. The history of their diversification and adaptation to these diverse environments remains controversial. We used 22 new genomes from 18 penguin species to reconstruct the order, timing, and location of their diversification, to track changes in their thermal niches through time, and to test for associated adaptation across the genome. Our results indicate that the penguin crown-group originated during the Miocene in New Zealand and Australia, not in Antarctica as previously thought, and that Aptenodytes is the sister group to all other extant penguin species. We show that lineage diversification in penguins was largely driven by changing climatic conditions and by the opening of the Drake Passage and associated intensification of the Antarctic Circumpolar Current (ACC). Penguin species have introgressed throughout much of their evolutionary history, following the direction of the ACC, which might have promoted dispersal and admixture. Changes in thermal niches were accompanied by adaptations in genes that govern thermoregulation and oxygen metabolism. Estimates of ancestral effective population sizes (Ne ) confirm that penguins are sensitive to climate shifts, as represented by three different demographic trajectories in deeper time, the most common (in 11 of 18 penguin species) being an increased Ne between 40 and 70 kya, followed by a precipitous decline during the Last Glacial Maximum. The latter effect is most likely a consequence of the overall decline in marine productivity following the last glaciation.

The evolutionary history of extinct and living lions.

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.

Genome-Wide SNPs Clarify a Complex Radiation and Support Recognition of an Additional Cat Species.

Phylogenetic reconstruction and species delimitation are often challenging in the case of recent evolutionary radiations, especially when postspeciation gene flow is present. Leopardus is a Neotropical cat genus that has a long history of recalcitrant taxonomic problems, along with both ancient and current episodes of interspecies admixture. Here, we employ genome-wide SNP data from all presently recognized Leopardus species, including several individuals from the tigrina complex (representing Leopardus guttulus and two distinct populations of Leopardus tigrinus), to investigate the evolutionary history of this genus. Our results reveal that the tigrina complex is paraphyletic, containing at least three distinct species. While one can be assigned to L. guttulus, the other two remain uncertain regarding their taxonomic assignment. Our findings indicate that the "tigrina" morphology may be plesiomorphic within this group, which has led to a longstanding taxonomic trend of lumping these poorly known felids into a single species.

The role of the environment in the spatial dynamics of an extensive hybrid zone between two neotropical cats.

Identifying factors that create and maintain a hybrid zone is of great interest to ecology, evolution and, more recently, conservation biology. Here, we investigated the role of environmental features in shaping the spatial dynamics of a hybrid zone between the southern tigrina, Leopardus guttulus, and Geoffroy's cat, L. geoffroyi, testing for exogenous selection as the main force acting on its maintenance. These Neotropical felid species are mainly allopatric, with a restricted area of sympatry in the ecotone between the Atlantic Forest and Pampa biomes. As both biomes have experienced high rates of anthropogenic habitat alteration, we also analysed the influence of habitat conversion on the hybrid zone structure. To do this, we used 13 microsatellite loci to identify potential hybrids and generated ecological niche models for them and their parental species. We compared the influence of variables on parental species and hybrid occurrence and calculated the amount of niche overlap among them. Parental species showed different habitat requirements and predicted co-occurrence was restricted to the forest-grassland mosaic of the ecotone. However, hybrids were found beyond this area, mainly in the range of L. geoffroyi. Hybrids demonstrated higher tolerance to habitat alteration than parental types, with a probability of occurrence that was positively related with mosaics of cropland areas and remnants of natural vegetation. These results indicate that exogenous selection alone does not drive the dynamics of the hybrid zone, and that habitat conversion influences its structure, potentially favouring hybrids over parental species.

Recombination-Aware Phylogenomics Reveals the Structured Genomic Landscape of Hybridizing Cat Species.

Current phylogenomic approaches implicitly assume that the predominant phylogenetic signal within a genome reflects the true evolutionary history of organisms, without assessing the confounding effects of postspeciation gene flow that can produce a mosaic of phylogenetic signals that interact with recombinational variation. Here, we tested the validity of this assumption with a phylogenomic analysis of 27 species of the cat family, assessing local effects of recombination rate on species tree inference and divergence time estimation across their genomes. We found that the prevailing phylogenetic signal within the autosomes is not always representative of the most probable speciation history, due to ancient hybridization throughout felid evolution. Instead, phylogenetic signal was concentrated within regions of low recombination, and notably enriched within large X chromosome recombination cold spots that exhibited recurrent patterns of strong genetic differentiation and selective sweeps across mammalian orders. By contrast, regions of high recombination were enriched for signatures of ancient gene flow, and these sequences inflated crown-lineage divergence times by ∼40%. We conclude that existing phylogenomic approaches to infer the Tree of Life may be highly misleading without considering the genomic architecture of phylogenetic signal relative to recombination rate and its interplay with historical hybridization.

Remarkably Complex Microbial Community Composition in Bromeliad Tank Waters Revealed by eDNA Metabarcoding.

To investigate patterns of biotic community composition at different spatial scales and biological contexts, we used environmental DNA metabarcoding to characterize eukaryotic and prokaryotic assemblages present in the phytotelmata of three bromeliad species (Aechmea gamosepala, Vriesea friburgensis, and Vriesea platynema) at a single Atlantic Forest site in southern Brazil. We sampled multiple individuals per species and multiple tanks from each individual, totalizing 30 samples. We observed very high levels of diversity in these communities, and remarkable variation across individuals and even among tanks from the same individual. The alpha diversity was higher for prokaryotes than eukaryotes, especially for A. gamosepala and V. platynema samples. Some biotic components appeared to be species-specific, while most of the biota was shared among species, but varied substantially in frequency among samples. Interestingly, V. friburgensis communities (which were sampled at nearby locations) tended to be more heterogeneous across samples, for both eukaryotes and prokaryotes. The opposite was true for V. platynema, whose samples were more broadly spaced but whose communities were more similar to each other. Our results indicate that additional attention should be devoted to within-individual heterogeneity when assessing bromeliad phytotelmata biodiversity, and highlight the complexity of the biotic assemblages gathered in these unique habitats.

Phylogenomic evidence for ancient hybridization in the genomes of living cats (Felidae).

Inter-species hybridization has been recently recognized as potentially common in wild animals, but the extent to which it shapes modern genomes is still poorly understood. Distinguishing historical hybridization events from other processes leading to phylogenetic discordance among different markers requires a well-resolved species tree that considers all modes of inheritance and overcomes systematic problems due to rapid lineage diversification by sampling large genomic character sets. Here, we assessed genome-wide phylogenetic variation across a diverse mammalian family, Felidae (cats). We combined genotypes from a genome-wide SNP array with additional autosomal, X- and Y-linked variants to sample ∼150 kb of nuclear sequence, in addition to complete mitochondrial genomes generated using light-coverage Illumina sequencing. We present the first robust felid time tree that accounts for unique maternal, paternal, and biparental evolutionary histories. Signatures of phylogenetic discordance were abundant in the genomes of modern cats, in many cases indicating hybridization as the most likely cause. Comparison of big cat whole-genome sequences revealed a substantial reduction of X-linked divergence times across several large recombination cold spots, which were highly enriched for signatures of selection-driven post-divergence hybridization between the ancestors of the snow leopard and lion lineages. These results highlight the mosaic origin of modern felid genomes and the influence of sex chromosomes and sex-biased dispersal in post-speciation gene flow. A complete resolution of the tree of life will require comprehensive genomic sampling of biparental and sex-limited genetic variation to identify and control for phylogenetic conflict caused by ancient admixture and sex-biased differences in genomic transmission.

Microbiota associated with tubes of Escarpia sp. from cold seeps in the southwestern Atlantic Ocean constitutes a community distinct from that of surrounding marine sediment and water.

As the depth increases and the light fades in oceanic cold seeps, a variety of chemosynthetic-based benthic communities arise. Previous assessments reported polychaete annelids belonging to the family Siboglinidae as part of the fauna at cold seeps, with the 'Vestimentifera' clade containing specialists that depend on microbial chemosynthetic endosymbionts for nutrition. Little information exists concerning the microbiota of the external portion of the vestimentiferan trunk wall. We employed 16S rDNA-based metabarcoding to describe the external microbiota of the chitin tubes from the vestimentiferan Escarpia collected from a chemosynthetic community in a cold seep area at the southwestern Atlantic Ocean. The most abundant operational taxonomic unit (OTU) belonged to the family Pirellulaceae (phylum Planctomycetes), and the second most abundant OTU belonged to the order Methylococcales (phylum Proteobacteria), composing an average of 21.1 and 15.4% of the total reads on tubes, respectively. These frequencies contrasted with those from the surrounding environment (sediment and water), where they represent no more than 0.1% of the total reads each. Moreover, some taxa with lower abundances were detected only in Escarpia tube walls. These data constitute on the first report of an epibiont microbial community found in close association with external surface of a cold-seep metazoan, Escarpia sp., from a chemosynthetic community in the southwestern Atlantic Ocean.

Recurrent evolution of melanism in South American felids.

Morphological variation in natural populations is a genomic test bed for studying the interface between molecular evolution and population genetics, but some of the most interesting questions involve non-model organisms that lack well annotated reference genomes. Many felid species exhibit polymorphism for melanism but the relative roles played by genetic drift, natural selection, and interspecies hybridization remain uncertain. We identify mutations of Agouti signaling protein (ASIP) or the Melanocortin 1 receptor (MC1R) as independent causes of melanism in three closely related South American species: the pampas cat (Leopardus colocolo), the kodkod (Leopardus guigna), and Geoffroy's cat (Leopardus geoffroyi). To assess population level variation in the regions surrounding the causative mutations we apply genomic resources from the domestic cat to carry out clone-based capture and targeted resequencing of 299 kb and 251 kb segments that contain ASIP and MC1R, respectively, from 54 individuals (13-21 per species), achieving enrichment of ~500-2500-fold and ~150x coverage. Our analysis points to unique evolutionary histories for each of the three species, with a strong selective sweep in the pampas cat, a distinctive but short melanism-specific haplotype in the Geoffroy's cat, and reduced nucleotide diversity for both ancestral and melanism-bearing chromosomes in the kodkod. These results reveal an important role for natural selection in a trait of longstanding interest to ecologists, geneticists, and the lay community, and provide a platform for comparative studies of morphological variation in other natural populations.

Phylogeographic analyses of the pampas cat (Leopardus colocola; Carnivora, Felidae) reveal a complex demographic history.

The pampas cat is a small felid that occurs in open habitats throughout much of South America. Previous studies have revealed intriguing patterns of morphological differentiation and genetic structure among its populations, as well as molecular evidence for hybridization with the closely related L. tigrinus. Here we report phylogeographic analyses encompassing most of its distribution (focusing particularly on Brazilian specimens, which had been poorly sampled in previous studies), using a novel dataset comprising 2,143 bp of the mitogenome, along with previously reported mtDNA sequences. Our data revealed strong population strutucture and supported a west-to-east colonization process in this species' history. We detected two population expansion events, one older (ca. 200 thousand years ago [kya]) in western South America and another more recent (ca. 60-50 kya) in eastern areas, coinciding with the expansion of savanna environments in Brazil. Analyses including L. tigrinus individuals bearing introgressed mtDNA from L. colocola showed a complete lack of shared haplotypes between species, indicating that their hybridization was ancient. Finally, we observed a close relationship between Brazilian/Uruguayan L. colocola haplotypes and those sampled in L. tigrinus, indicating that their hybridization was likely related to the demographic expansion of L. colocola into eastern South America.

Waking the undead: Implications of a soft explosive model for the timing of placental mammal diversification.

The explosive, long fuse, and short fuse models represent competing hypotheses for the timing of placental mammal diversification. Support for the explosive model, which posits both interordinal and intraordinal diversification after the KPg mass extinction, derives from morphological cladistic studies that place Cretaceous eutherians outside of crown Placentalia. By contrast, most molecular studies favor the long fuse model wherein interordinal cladogenesis occurred in the Cretaceous followed by intraordinal cladogenesis after the KPg boundary. Phillips (2016) proposed a soft explosive model that allows for the emergence of a few lineages (Xenarthra, Afrotheria, Euarchontoglires, Laurasiatheria) in the Cretaceous, but otherwise agrees with the explosive model in positing the majority of interordinal diversification after the KPg mass extinction. Phillips (2016) argues that rate transference errors associated with large body size and long lifespan have inflated previous estimates of interordinal divergence times, and further suggests that most interordinal divergences are positioned after the KPg boundary when rate transference errors are avoided through the elimination of calibrations in large-bodied and/or long lifespan clades. Here, we show that rate transference errors can also occur in the opposite direction and drag forward estimated divergence dates when calibrations in large-bodied/long lifespan clades are omitted. This dragging forward effect results in the occurrence of more than half a billion years of 'zombie lineages' on Phillips' preferred timetree. By contrast with ghost lineages, which are a logical byproduct of an incomplete fossil record, zombie lineages occur when estimated divergence dates are younger than the minimum age of the oldest crown fossils. We also present the results of new timetree analyses that address the rate transference problem highlighted by Phillips (2016) by deleting taxa that exceed thresholds for body size and lifespan. These analyses recover all interordinal divergence times in the Cretaceous and are consistent with the long fuse model of placental diversification. Finally, we outline potential problems with morphological cladistic analyses of higher-level relationships among placental mammals that may account for the perceived discrepancies between molecular and paleontological estimates of placental divergence times.

Epistylis portoalegrensis n. sp. (Ciliophora, Peritrichia): A New Freshwater Ciliate Species from Southern Brazil.

The peritrich ciliate Epistylis portoalegrensis n. sp. was found in two bodies of freshwater located in Porto Alegre, Southern Brazil. Morphological features were investigated using live and protargol-stained specimens. The zooids presented a vase to cylindrical shape narrowed at the scopula, and a mean size of 131 × 37 µm in vivo. A C-shaped macronucleus lay in the middle of the cell close to a single contractile vacuole. The oral infraciliature was typical for the genus, with all infundibular polykineties composed by three distinct rows of kinetosomes. Colonies are often nonbranched with no lateral stalk, carrying several zooids stemming from a single point. Specimens from the two sampling sites showed identical arrangement of the infraciliature, similar morphometry, identical 18S rDNA sequences, and a single nucleotide difference across the more variable ITS regions. Molecular phylogenetic analyses placed E. portoalegrensis in a well-supported clade containing other Epistylis species, within the order Vorticellida.