Search

Abstract Bats are becoming recognised as new model species to study naturally evolved mammalian extended healthspan and disease tolerance. However, this research is limited by the lack of bat specific cellular resources. Here we describe an optimised protocol to develop both primary and immortalised fibroblast cell-lines from wing biopsy punches from the Egyptian fruit bat, Rousettus aegyptiacus. We show that the immortalised cell lines and primary cells show similar characteristics in their proliferative capacity and response to oxidative stress. They also exhibited a similar response in their NF-κB immune response to TLR agonists including SARS-CoV2. As wing punches can be acquired non-lethally, these methods can be used to develop primary and immortalised cells, from potentially any bat species, including those of conservation concern that cannot be sacrificed. This can expand the scope of bat species that can be studied in the future, and the development of key cellular resources required to functionally validate the regulators of bats’ unique longevity.

We present a genome assembly from an individual male Myotis mystacinus (whiskered bat; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence has a total length of 2,081.20 megabases. Most of the assembly (97.52%) is scaffolded into 23 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 16.93 kilobases in length.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Abstract Phytoplankton blooms provoke bacterioplankton blooms, from which bacterial biomass (necromass) is released via increased zooplankton grazing and viral lysis. While bacterial consumption of algal biomass during blooms is well-studied, little is known about the concurrent recycling of these substantial amounts of bacterial necromass. We demonstrate that bacterial biomass, such as bacterial alpha-glucan storage polysaccharides, generated from the consumption of algal organic matter, is reused and thus itself a major bacterial carbon source in vitro and during a diatom-dominated bloom. We highlight conserved enzymes and binding proteins of dominant bloom-responder clades that are presumably involved in the recycling of bacterial alpha-glucan by members of the bacterial community. We furthermore demonstrate that the corresponding protein machineries can be specifically induced by extracted alpha-glucan-rich bacterial polysaccharide extracts. This recycling of bacterial necromass likely constitutes a large-scale intra-population energy conservation mechanism that keeps substantial amounts of carbon in a dedicated part of the microbial loop.

We present a genome assembly from an individual female Molossus alvarezi (Chordata; Mammalia; Chiroptera; Molossidae). The genome sequence is 2.490 Gb in span. The majority of the assembly is scaffolded into 24 chromosomal pseudomolecules, with the X sex chromosomes assembled.

We present a reference genome assembly from an individual male Rhynchonycteris naso (Chordata; Mammalia; Chiroptera; Emballonuridae). The genome sequence is 2.46 Gb in span. The majority of the assembly is scaffolded into 22 chromosomal pseudomolecules, with the Y sex chromosome assembled.

We present a genome assembly from an individual male Vespertilio murinus (the particolored bat; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence is 1,925.6 megabases in span. Most of the assembly is scaffolded into 20 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 16.96 kilobases in length.

Abstract Background Marine microalgae (phytoplankton) mediate almost half of the worldwide photosynthetic carbon dioxide fixation and therefore play a pivotal role in global carbon cycling, most prominently during massive phytoplankton blooms. Phytoplankton biomass consists of considerable proportions of polysaccharides, substantial parts of which are rapidly remineralized by heterotrophic bacteria. We analyzed the diversity, activity, and functional potential of such polysaccharide-degrading bacteria in different size fractions during a diverse spring phytoplankton bloom at Helgoland Roads (southern North Sea) at high temporal resolution using microscopic, physicochemical, biodiversity, metagenome, and metaproteome analyses. Results Prominent active 0.2–3 µm free-living clades comprised Aurantivirga, “Formosa”, Cd. Prosiliicoccus, NS4, NS5, Amylibacter, Planktomarina, SAR11 Ia, SAR92, and SAR86, whereas BD1-7, Stappiaceae, Nitrincolaceae, Methylophagaceae, Sulfitobacter, NS9, Polaribacter, Lentimonas, CL500-3, Algibacter, and Glaciecola dominated 3–10 µm and > 10 µm particles. Particle-attached bacteria were more diverse and exhibited more dynamic adaptive shifts over time in terms of taxonomic composition and repertoires of encoded polysaccharide-targeting enzymes. In total, 305 species-level metagenome-assembled genomes were obtained, including 152 particle-attached bacteria, 100 of which were novel for the sampling site with 76 representing new species. Compared to free-living bacteria, they featured on average larger metagenome-assembled genomes with higher proportions of polysaccharide utilization loci. The latter were predicted to target a broader spectrum of polysaccharide substrates, ranging from readily soluble, simple structured storage polysaccharides (e.g., laminarin, α-glucans) to less soluble, complex structural, or secreted polysaccharides (e.g., xylans, cellulose, pectins). In particular, the potential to target poorly soluble or complex polysaccharides was more widespread among abundant and active particle-attached bacteria. Conclusions Particle-attached bacteria represented only 1% of all bloom-associated bacteria, yet our data suggest that many abundant active clades played a pivotal gatekeeping role in the solubilization and subsequent degradation of numerous important classes of algal glycans. The high diversity of polysaccharide niches among the most active particle-attached clades therefore is a determining factor for the proportion of algal polysaccharides that can be rapidly remineralized during generally short-lived phytoplankton bloom events. Video Abstract

We present a genome assembly from a female Plecotus auritus (Brown Long-eared bat; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence is 2163.2 megabases in span. Most of the assembly is scaffolded into 16 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 16.91 kilobases in length.

Small genes (40,000 small-gene families in ∼2.3 million phage genome contigs. We find that small genes in phage genomes are approximately 3-fold more prevalent than in host prokaryotic genomes. Our approach enriches for small genes that are translated in microbiomes, suggesting the small genes identified are coding. More than 9,000 families encode potentially secreted or transmembrane proteins, more than 5,000 families encode predicted anti-CRISPR proteins, and more than 500 families encode predicted antimicrobial proteins. By combining homology and genomic-neighborhood analyses, we reveal substantial novelty and diversity within phage biology, including small phage genes found in multiple host phyla, small genes encoding proteins that play essential roles in host infection, and small genes that share genomic neighborhoods and whose encoded proteins may share related functions.

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

We present a genome assembly from an individual male Myotis daubentonii (Daubenton's bat; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence is 2,127.8 megabases in span. Most of the assembly is scaffolded into 23 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.34 kilobases in length.

Chronic Rhinosinusitis with Nasal Polyps (CRSwNP) is a chronic inflammatory disease of the nose and paranasal sinus cavities that significantly affects well-being and social function, particularly in young adults and middle-aged populations. CRSwNP is a common health condition in the Western world, with an estimated prevalence of 3%. Despite worldwide evidence-based treatment guidelines such as the European Position Paper on Rhinosinusitis and Nasal Polyps (EPOS) 2020 and the European Forum for Research and Education in Allergy and Airway Diseases (EUFOREA) chronic rhinosinusitis (CRS) pocket guide, a significant number of patients remain undiagnosed and/or uncontrolled with repeated oral corticosteroids (OCS) treatments and/or (multiple) endoscopic sinus surgeries (ESS).

We present a genome assembly from an individual male Tadarida brasiliensis (The Brazilian free-tailed bat; Chordata; Mammalia; Chiroptera; Molossidae). The genome sequence is 2.28 Gb in span. The majority of the assembly is scaffolded into 25 chromosomal pseudomolecules, with the X and Y sex chromosomes assembled.

Abstract Background Macroalgal epiphytic microbial communities constitute a rich resource for novel enzymes and compounds, but studies so far largely focused on tag-based microbial diversity analyses or limited metagenome sequencing of single macroalgal species. Results We sampled epiphytic bacteria from specimens of Ulva sp. (green algae), Saccharina sp. (brown algae), Grateloupia sp. and Gelidium sp. (both red algae) together with seawater and sediment controls from a coastal reef in Weihai, China, during all seasons. Using 16S rRNA amplicon sequencing, we identified 14 core genera (consistently present on all macroalgae), and 14 dominant genera (consistently present on three of the macroalgae). Core genera represented ~ 0.7% of all genera, yet accounted for on average 51.1% of the bacterial abundances. Plate cultivation from all samples yielded 5,527 strains (macroalgae: 4,426) representing 1,235 species (685 potentially novel). Sequencing of selected strains yielded 820 non-redundant draft genomes (506 potentially novel), and sequencing of 23 sampled metagenomes yielded 1,619 metagenome-assembled genomes (MAGs), representing further 1,183 non-redundant genomes. 230 isolates and 153 genomes were obtained from the 28 core/dominant genera. We analyzed the genomic potential of phycosphere bacteria to degrade algal polysaccharides and to produce bioactive secondary metabolites. We predicted 4,451 polysaccharide utilization loci (PULs) and 8,810 biosynthetic gene clusters (BGCs). These were particularly prevalent in core/dominant genera. Conclusions Our metabolic annotations and analyses of MAGs and genomes provide new insights into novel species of phycosphere bacteria and their ecological niches for an improved understanding of the macroalgal phycosphere microbiome. Video Abstract

Abstract Background Blooms of marine microalgae play a pivotal role in global carbon cycling. Such blooms entail successive blooms of specialized clades of planktonic bacteria that collectively remineralize gigatons of algal biomass on a global scale. This biomass is largely composed of distinct polysaccharides, and the microbial decomposition of these polysaccharides is therefore a process of prime importance. Results In 2020, we sampled a complete biphasic spring bloom in the German Bight over a 90-day period. Bacterioplankton metagenomes from 30 time points allowed reconstruction of 251 metagenome-assembled genomes (MAGs). Corresponding metatranscriptomes highlighted 50 particularly active MAGs of the most abundant clades, including many polysaccharide degraders. Saccharide measurements together with bacterial polysaccharide utilization loci (PUL) expression data identified β-glucans (diatom laminarin) and α-glucans as the most prominent and actively metabolized dissolved polysaccharide substrates. Both substrates were consumed throughout the bloom, with α-glucan PUL expression peaking at the beginning of the second bloom phase shortly after a peak in flagellate and the nadir in bacterial total cell counts. Conclusions We show that the amounts and composition of dissolved polysaccharides, in particular abundant storage polysaccharides, have a pronounced influence on the composition of abundant bacterioplankton members during phytoplankton blooms, some of which compete for similar polysaccharide niches. We hypothesize that besides the release of algal glycans, also recycling of bacterial glycans as a result of increased bacterial cell mortality can have a significant influence on bacterioplankton composition during phytoplankton blooms. Video Abstract

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19. The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of ACE2 sequences from 410 vertebrate species, including 252 mammals, to study the conservation of ACE2 and its potential to be used as a receptor by SARS-CoV-2. We designed a five-category binding score based on the conservation properties of 25 amino acids important for the binding between ACE2 and the SARS-CoV-2 spike protein. Only mammals fell into the medium to very high categories and only catarrhine primates into the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 spike protein binding and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (frequency

Abstract There are no disease-modifying treatments available for geographic atrophy (GA), the advanced form of dry age-related macular degeneration. Current murine models fail to fully recapitulate the features of GA and thus hinder drug discovery. Here we describe a novel mouse model of retinal degeneration with hallmark features of GA. We used an 810 nm laser to create a retinal lesion with central sparing (RLCS), simulating parafoveal atrophy observed in patients with progressive GA. Laser-induced RLCS resulted in progressive GA-like pathology with the development of a confluent atrophic lesion. We demonstrate significant changes to the retinal structure and thickness in the central unaffected retina over a 24-week post-laser period, confirmed by longitudinal optical coherence tomography scans. We further show characteristic features of progressive GA, including a gradual reduction in the thickness of the central, unaffected retina and of total retinal thickness. Histological changes observed in the RLCS correspond to GA pathology, which includes the collapse of the outer nuclear layer, increased numbers of GFAP + , CD11b + and FcγRI + cells, and damage to cone and rod photoreceptors. We demonstrate a laser-induced mouse model of parafoveal GA progression, starting at 2 weeks post-laser and reaching confluence at 24 weeks post-laser. This 24-week time-frame in which GA pathology develops, provides an extended window of opportunity for proof-of-concept evaluation of drugs targeting GA. This time period is an added advantage compared to several existing models of geographic atrophy.

We present a genome assembly from an individual Eptesicus nilssonii (the northern bat; Chordata; Mammalia; Chiroptera; Vespertilionidae), derived from the placental tissue of a pregnancy that resulted a male pup. The genome sequence is 2,064.1 megabases in span. Most of the assembly is scaffolded into 26 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.04 kilobases in length.

We present a genome assembly from an individual male Pipistrellus pygmaeus (the Soprano Pipistrelle; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence is 1,895.1 megabases in span. Most of the assembly is scaffolded into 23 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.18 kilobases in length.

Summary: Bats host a range of disease-causing viruses without displaying clinical symptoms. The mechanisms behind this are a continuous source of interest. Here, we studied the antiviral response in the Egyptian fruit bat and Kuhl’s pipistrelle, representing two subordinal clades. We profiled the antiviral response in fibroblasts using RNA sequencing and compared bat with primate and rodent responses. Both bats upregulate similar genes; however, a subset of these genes is transcriptionally divergent between them. These divergent genes also evolve rapidly in sequence, have specific promoter architectures, and are associated with programs underlying tolerance and resistance. Finally, we characterized antiviral genes that expanded in bats, with duplicates diverging in sequence and expression.Our study reveals a largely conserved antiviral program across bats and points to a set of genes that rapidly evolve through multiple mechanisms. These can contribute to bat adaptation to viral infection and provide directions to understanding the mechanisms behind it.

In March 2023, the European Forum for Research and Education in Allergy and Airways diseases (EUFOREA) organized its bi-annual Summit in Brussels with expert panel members of EUFOREA, representatives of the EUFOREA patient advisory board, and the EUFOREA board and management teams. Its aim was to define the research, educational and advocacy initiatives to be developed by EUFOREA over the next 2 years until the 10th anniversary in 2025. EUFOREA is an international non-for-profit organization forming an alliance of all stakeholders dedicated to reducing the prevalence and burden of chronic allergic and respiratory diseases via research, education, and advocacy. Based on its medical scientific core competency, EUFOREA offers an evidence-supported platform to introduce innovation and education in healthcare leading to optimal patient care, bridging the gap between latest scientific evidence and daily practice. Aligned with the mission of improving health care, the expert panels of asthma, allergic rhinitis (AR), chronic rhinosinusitis (CRS) & European Position Paper on Rhinosinusitis and Nasal Polyps (EPOS), allergen immunotherapy (AIT) and paediatrics have proposed and elaborated a variety of activities that correspond to major unmet needs in the allergy and respiratory field. The current report provides a concise overview of the achievements, ambitions, and action plan of EUFOREA for the future, allowing all stakeholders in the allergy and respiratory field to be up-dated and inspired to join forces in Europe and beyond.

We present a genome assembly from an individual male Molossus nigricans (Chordata; Mammalia; Chiroptera; Molossidae). The genome sequence is 2.41 gigabases in span. The majority of the assembly is scaffolded into 24 chromosomal pseudomolecules, with the X sex chromosome assembled.

Computational thinking" (CT) is still a relatively new term in the lexicon of learning objectives and science standards. The term was popularized in an essay by Wing, who said, "To reading, writing and arithmetic, we should add computational thinking to every child's analytical ability." Agreeing with this premise, in 2013 the authors of the Next Generation Science Standards (NGSS) included "mathematical and computational thinking" as one of eight essential science and engineering practices that K-12 teachers should strive to develop in their students. There is not yet widespread agreement on the precise definition or implementation of CT, and efforts to assess CT are still maturing, even as more states adopt K-12 computer science standards. In this article we will try to summarize what CT means for a typical introductory (i.e., high school or early college) physics class. This will include a discussion of the ways that instructors may already be incorporating elements of CT in their classes without knowing it.

ABSTRACT Some marine thermophilic methanogens are able to perform energy-consuming nitrogen fixation despite deriving only little energy from hydrogenotrophic methanogenesis. We studied this process in Methanothermococcus thermolithotrophicus DSM 2095, a methanogenic archaeon of the order Methanococcales that contributes to the nitrogen pool in some marine environments. We successfully grew this archaeon under diazotrophic conditions in both batch and fermenter cultures, reaching the highest cell density reported so far. Diazotrophic growth depended strictly on molybdenum and, in contrast to other diazotrophs, was not inhibited by tungstate or vanadium. This suggests an elaborate control of metal uptake and a specific metal recognition system for the insertion into the nitrogenase cofactor. Differential transcriptomics of M. thermolithotrophicus grown under diazotrophic conditions with ammonium-fed cultures as controls revealed upregulation of the nitrogenase machinery, including chaperones, regulators, and molybdate importers, as well as simultaneous upregulation of an ammonium transporter and a putative pathway for nitrate and nitrite utilization. The organism thus employs multiple synergistic strategies for uptake of nitrogen nutrients during the early exponential growth phase without altering transcription levels for genes involved in methanogenesis. As a counterpart, genes coding for transcription and translation processes were downregulated, highlighting the maintenance of an intricate metabolic balance to deal with energy constraints and nutrient limitations imposed by diazotrophy. This switch in the metabolic balance included unexpected processes, such as upregulation of the CRISPR-Cas system, probably caused by drastic changes in transcription levels of putative mobile and virus-like elements. IMPORTANCE The thermophilic anaerobic archaeon M. thermolithotrophicus is a particularly suitable model organism to study the coupling of methanogenesis to diazotrophy. Likewise, its capability of simultaneously reducing N2 and CO2 into NH3 and CH4 with H2 makes it a viable target for biofuel production. We optimized M. thermolithotrophicus cultivation, resulting in considerably higher cell yields and enabling the successful establishment of N2-fixing bioreactors. Improved understanding of the N2 fixation process would provide novel insights into metabolic adaptations that allow this energy-limited extremophile to thrive under diazotrophy, for instance, by investigating its physiology and uncharacterized nitrogenase. We demonstrated that diazotrophic growth of M. thermolithotrophicus is exclusively dependent on molybdenum, and complementary transcriptomics corroborated the expression of the molybdenum nitrogenase system. Further analyses of differentially expressed genes during diazotrophy across three cultivation time points revealed insights into the response to nitrogen limitation and the coordination of core metabolic processes.

Abstract Background The planktonic bacterial community associated with spring phytoplankton blooms in the North Sea is responsible for a large amount of carbon turnover in an environment characterised by high primary productivity. Individual clades belonging to the Gammaproteobacteria have shown similar population dynamics to Bacteroidetes species, and are thus assumed to fill competing ecological niches. Previous studies have generated large numbers of metagenome assembled genomes and metaproteomes from these environments, which can be readily mined to identify populations performing potentially important ecosystem functions. In this study we attempt to catalogue these spring bloom-associated Gammaproteobacteria, which have thus far attracted less attention than sympatric Alphaproteobacteria and Bacteroidetes. Methods We annotated 120 non-redundant species-representative gammaproteobacterial metagenome assembled genomes from spring bloom sampling campaigns covering the four years 2010–2012 and 2016 using a combination of Prokka and PfamScan, with further confirmation via BLAST against NCBI-NR. We also matched these gene annotations to 20 previously published metaproteomes covering those sampling periods plus the spring of 2009. Results Metagenome assembled genomes with clear capacity for polysaccharide degradation via dedicated clusters of carbohydrate active enzymes were among the most abundant during blooms. Many genomes lacked gene clusters with clearly identifiable predicted polysaccharide substrates, although abundantly expressed loci for the uptake of large molecules were identified in metaproteomes. While the larger biopolymers, which are the most abundant sources of reduced carbon following algal blooms, are likely the main energy source, some gammaproteobacterial clades were clearly specialised for smaller organic compounds. Their substrates range from amino acids, monosaccharides, and DMSP, to the less expected, such as terpenoids, and aromatics and biphenyls, as well as many ‘unknowns’. In particular we uncover a much greater breadth of apparent methylotrophic capability than heretofore identified, present in several order level clades without cultivated representatives. Conclusions Large numbers of metagenome assembled genomes are today publicly available, containing a wealth of readily accessible information. Here we identified a variety of predicted metabolisms of interest, which include diverse potential heterotrophic niches of spring bloom-associated Gammaproteobacteria. Features such as those identified here could well be fertile ground for future experimental studies.

The great cats of the genus Panthera comprise a recent radiation whose evolutionary history is poorly understood. Their rapid diversification poses challenges to resolving their phylogeny while offering opportunities to investigate the historical dynamics of adaptive divergence. We report the sequence, de novo assembly, and annotation of the jaguar (Panthera onca) genome, a novel genome sequence for the leopard (Panthera pardus), and comparative analyses encompassing all living Panthera species. Demographic reconstructions indicated that all of these species have experienced variable episodes of population decline during the Pleistocene, ultimately leading to small effective sizes in present-day genomes. We observed pervasive genealogical discordance across Panthera genomes, caused by both incomplete lineage sorting and complex patterns of historical interspecific hybridization. We identified multiple signatures of species-specific positive selection, affecting genes involved in craniofacial and limb development, protein metabolism, hypoxia, reproduction, pigmentation, and sensory perception. There was remarkable concordance in pathways enriched in genomic segments implicated in interspecies introgression and in positive selection, suggesting that these processes were connected. We tested this hypothesis by developing exome capture probes targeting ~19,000 Panthera genes and applying them to 30 wild-caught jaguars. We found at least two genes (DOCK3 and COL4A5, both related to optic nerve development) bearing significant signatures of interspecies introgression and within-species positive selection. These findings indicate that post-speciation admixture has contributed genetic material that facilitated the adaptive evolution of big cat lineages.