Microbiome diversity and zoonotic bacterial pathogen prevalence in Peromyscus mice from agricultural landscapes and synanthropic habitat.

Rodents are key reservoirs of zoonotic pathogens and play an important role in disease transmission to humans. Importantly, anthropogenic land-use change has been found to increase the abundance of rodents that thrive in human-built environments (synanthropic rodents), particularly rodent reservoirs of zoonotic disease. Anthropogenic environments also affect the microbiome of synanthropic wildlife, influencing wildlife health and potentially introducing novel pathogens. Our objective was to examine the effect of agricultural development and synanthropic habitat on microbiome diversity and the prevalence of zoonotic bacterial pathogens in wild Peromyscus mice to better understand the role of these rodents in pathogen maintenance and transmission. We conducted 16S amplicon sequencing on faecal samples using long-read nanopore sequencing technology to characterize the rodent microbiome. We compared microbiome diversity and composition between forest and synanthropic habitats in agricultural and undeveloped landscapes and screened for putative pathogenic bacteria. Microbiome richness, diversity, and evenness were higher in the agricultural landscape and synanthropic habitat compared to undeveloped-forest habitat. Microbiome composition also differed significantly between agricultural and undeveloped landscapes and forest and synanthropic habitats. We detected overall low diversity and abundance of putative pathogenic bacteria, though putative pathogens were more likely to be found in mice from the agricultural landscape. Our findings show that landscape- and habitat-level anthropogenic factors affect Peromyscus microbiomes and suggest that landscape-level agricultural development may be important to predict zoonotic pathogen prevalence. Ultimately, understanding how anthropogenic land-use change and synanthropy affect rodent microbiomes and pathogen prevalence is important to managing transmission of rodent-borne zoonotic diseases to humans.

The establishment of a collaborative surveillance program with indigenous hunters to characterize primate health in Southern Guyana.

The consumption of primates is integral to the traditional subsistence strategies of many Indigenous communities throughout Amazonia. Understanding the overall health of primates harvested for food in the region is critical to Indigenous food security and thus, these communities are highly invested in long-term primate population health. Here, we describe the establishment of a surveillance comanagement program among the Waiwai, an Indigenous community in the Konashen Amerindian Protected Area (KAPA). To assess primate health in the KAPA, hunters performed field necropsies on primates harvested for food and tissues collected from these individuals were analyzed using histopathology. From 2015 to 2019, hunters conducted 127 necropsies across seven species of primates. Of this sample, 82 primates (between 2015 and 2017) were submitted for histopathological screening. Our histopathology data revealed that KAPA primates had little evidence of underlying disease. Of the tissue abnormalities observed, the majority were either due to diet (e.g., hepatocellular pigment), degenerative changes resulting from aging (e.g., interstitial nephritis, myocyte lipofusion), or nonspecific responses to antigenic stimulation (renal and splenic lymphoid hyperplasia). In our sample, 7.32% of individuals had abnormalities that were consistent with a viral etiology, including myocarditis and hepatitis. Internal parasites were observed in 53.66% of individuals and is consistent with what would be expected from a free-ranging primate population. This study represents the importance of baseline data for long-term monitoring of primate populations hunted for food. More broadly, this research begins to close a critical gap in zoonotic disease risk related to primate harvesting in Amazonia, while also demonstrating the benefits of partnering with Indigenous hunters and leveraging hunting practices in disease surveillance and primate population health assessment.

Nanoparticle-Enhanced RT-QuIC (Nano-QuIC) Diagnostic Assay for Misfolded Proteins.

Misfolded proteins associated with various neurodegenerative diseases often accumulate in tissues or circulate in biological fluids years before the clinical onset, thus representing ideal diagnostic targets. Real-time quaking-induced conversion (RT-QuIC), a protein-based seeded-amplification assay, holds great potential for early disease detection, yet challenges remain for routine diagnostic application. Chronic Wasting Disease (CWD), associated with misfolded prion proteins of cervids, serves as an ideal model for evaluating new RT-QuIC methodologies. In this study, we investigate the previously untested hypothesis that incorporating nanoparticles into RT-QuIC assays can enhance their speed and sensitivity when applied to biological samples. We show that adding 50 nm silica nanoparticles to RT-QuIC experiments (termed Nano-QuIC) for CWD diagnostics greatly improves the performance by reducing detection times 2.5-fold and increasing sensitivity 10-fold by overcoming the effect of inhibitors in complex tissue samples. Crucially, no false positives were observed with these 50 nm silica nanoparticles, demonstrating the enhanced reliability and potential for diagnostic application of Nano-QuIC in detecting misfolded proteins.

A De Novo Whole Genome Assembly and Annotation of Parelaphostrongylus tenuis.

Parelaphostrongylus tenuis causes ungulate morbidity and mortality in eastern and central North America, but no reference genome sequence exists to facilitate research. Here, we present a P. tenuis genome assembly and annotation, generated with PacBio and Illumina technologies. The assembly is 491 Mbp, with 7285 scaffolds and 185 kb N50.

Single-cell sequencing of entorhinal cortex reveals widespread disruption of neuropeptide networks in Alzheimer's disease.

INTRODUCTION: Abnormalities of neuropeptides (NPs) that play important roles in modulating neuronal activities are commonly observed in Alzheimer's disease (AD). We hypothesize that NP network disruption is widespread in AD brains. METHODS: Single-cell transcriptomic data from the entorhinal cortex (EC) were used to investigate the NP network disruption in AD. Bulk RNA-sequencing data generated from the temporal cortex by independent groups and machine learning were employed to identify key NPs involved in AD. The relationship between aging and AD-associated NP (ADNP) expression was studied using GTEx data. RESULTS: The proportion of cells expressing NPs but not their receptors decreased significantly in AD. Neurons expressing higher level and greater diversity of NPs were disproportionately absent in AD. Increased age coincides with decreased ADNP expression in the hippocampus. DISCUSSION: NP network disruption is widespread in AD EC. Neurons expressing more NPs may be selectively vulnerable to AD. Decreased expression of NPs participates in early AD pathogenesis. We predict that the NP network can be harnessed for treatment and/or early diagnosis of AD.

Alu retrotransposons and COVID-19 susceptibility and morbidity.

SARS-CoV-2 has spread rapidly across the world and is negatively impacting the global human population. COVID-19 patients display a wide variety of symptoms and clinical outcomes, including those attributed to genetic ancestry. Alu retrotransposons have played an important role in human evolution, and their variants influence host response to viral infection. Intronic Alus regulate gene expression through several mechanisms, including both genetic and epigenetic pathways. With respect to SARS-CoV-2, an intronic Alu within the ACE gene is hypothesized to be associated with COVID-19 susceptibility and morbidity. Here, we review specific Alu polymorphisms that are of particular interest when considering host response to SARS-CoV-2 infection, especially polymorphic Alu insertions in genes associated with immune response and coagulation/fibrinolysis cascade. We posit that additional research focused on Alu-related pathways could yield novel biomarkers capable of predicting clinical outcomes as well as patient-specific treatment strategies for COVID-19 and related infectious diseases.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent.

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

The mitochondrial genome and Epigenome of the Golden lion Tamarin from fecal DNA using Nanopore adaptive sequencing.

BACKGROUND: The golden lion tamarin (Leontopithecus rosalia) is an endangered Platyrrhine primate endemic to the Atlantic coastal forests of Brazil. Despite ongoing conservation efforts, genetic data on this species remains scarce. Complicating factors include limitations on sample collection and a lack of high-quality reference sequences. Here, we used nanopore adaptive sampling to resequence the L. rosalia mitogenome from feces, a sample which can be collected non-invasively. RESULTS: Adaptive sampling doubled the fraction of both host-derived and mitochondrial sequences compared to sequencing without enrichment. 258x coverage of the L. rosalia mitogenome was achieved in a single flow cell by targeting the unfinished genome of the distantly related emperor tamarin (Saguinus imperator) and the mitogenome of the closely related black lion tamarin (Leontopithecus chrysopygus). The L. rosalia mitogenome has a length of 16,597 bp, sharing 99.68% sequence identity with the L. chrysopygus mitogenome. A total of 38 SNPs between them were identified, with the majority being found in the non-coding D-loop region. DNA methylation and hydroxymethylation were directly detected using a neural network model applied to the raw signal from the MinION sequencer. In contrast to prior reports, DNA methylation was negligible in mitochondria in both CpG and non-CpG contexts. Surprisingly, a quarter of the 642 CpG sites exhibited DNA hydroxymethylation greater than 1% and 44 sites were above 5%, with concentration in the 3' side of several coding regions. CONCLUSIONS: Overall, we report a robust new mitogenome assembly for L. rosalia and direct detection of cytosine base modifications in all contexts.

Comparative analyses of two primate species diverged by more than 60 million years show different rates but similar distribution of genome-wide UV repair events.

BACKGROUND: Nucleotide excision repair is the primary DNA repair mechanism that removes bulky DNA adducts such as UV-induced pyrimidine dimers. Correspondingly, genome-wide mapping of nucleotide excision repair with eXcision Repair sequencing (XR-seq), provides comprehensive profiling of DNA damage repair. A number of XR-seq experiments at a variety of conditions for different damage types revealed heterogenous repair in the human genome. Although human repair profiles were extensively studied, how repair maps vary between primates is yet to be investigated. Here, we characterized the genome-wide UV-induced damage repair in gray mouse lemur, Microcebus murinus, in comparison to human. RESULTS: We derived fibroblast cell lines from mouse lemur, exposed them to UV irradiation, and analyzed the repair events genome-wide using the XR-seq protocol. Mouse lemur repair profiles were analyzed in comparison to the equivalent human fibroblast datasets. We found that overall UV sensitivity, repair efficiency, and transcription-coupled repair levels differ between the two primates. Despite this, comparative analysis of human and mouse lemur fibroblasts revealed that genome-wide repair profiles of the homologous regions are highly correlated, and this correlation is stronger for highly expressed genes. With the inclusion of an additional XR-seq sample derived from another human cell line in the analysis, we found that fibroblasts of the two primates repair UV-induced DNA lesions in a more similar pattern than two distinct human cell lines do. CONCLUSION: Our results suggest that mouse lemurs and humans, and possibly primates in general, share a homologous repair mechanism as well as genomic variance distribution, albeit with their variable repair efficiency. This result also emphasizes the deep homologies of individual tissue types across the eukaryotic phylogeny.

Pedigree-based and phylogenetic methods support surprising patterns of mutation rate and spectrum in the gray mouse lemur.

Mutations are the raw material on which evolution acts, and knowledge of their frequency and genomic distribution is crucial for understanding how evolution operates at both long and short timescales. At present, the rate and spectrum of de novo mutations have been directly characterized in relatively few lineages. Our study provides the first direct mutation-rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprises nearly half of the primate clade. Using high-coverage linked-read sequencing for a focal quartet of gray mouse lemurs (Microcebus murinus), we estimated the mutation rate to be among the highest calculated for a mammal at 1.52 × 10-8 (95% credible interval: 1.28 × 10-8-1.78 × 10-8) mutations/site/generation. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG sites. Despite the surprising nature of these results, we found both the rate and spectrum to be robust to the manipulation of a wide range of computational filtering criteria. We also sequenced a technical replicate to estimate a false-negative and false-positive rate for our data and show that any point estimate of a de novo mutation rate should be considered with a large degree of uncertainty. For validation, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species for which de novo mutation rates have also been estimated. These comparisons revealed general consistency of the mutation spectrum between the pedigree-based and the substitution-rate analyses for all species compared.

Warning SINEs: Alu elements, evolution of the human brain, and the spectrum of neurological disease.

Alu elements are a highly successful family of primate-specific retrotransposons that have fundamentally shaped primate evolution, including the evolution of our own species. Alus play critical roles in the formation of neurological networks and the epigenetic regulation of biochemical processes throughout the central nervous system (CNS), and thus are hypothesized to have contributed to the origin of human cognition. Despite the benefits that Alus provide, deleterious Alu activity is associated with a number of neurological and neurodegenerative disorders. In particular, neurological networks are potentially vulnerable to the epigenetic dysregulation of Alu elements operating across the suite of nuclear-encoded mitochondrial genes that are critical for both mitochondrial and CNS function. Here, we highlight the beneficial neurological aspects of Alu elements as well as their potential to cause disease by disrupting key cellular processes across the CNS. We identify at least 37 neurological and neurodegenerative disorders wherein deleterious Alu activity has been implicated as a contributing factor for the manifestation of disease, and for many of these disorders, this activity is operating on genes that are essential for proper mitochondrial function. We conclude that the epigenetic dysregulation of Alu elements can ultimately disrupt mitochondrial homeostasis within the CNS. This mechanism is a plausible source for the incipient neuronal stress that is consistently observed across a spectrum of sporadic neurological and neurodegenerative disorders.

Geogenetic patterns in mouse lemurs (genus Microcebus) reveal the ghosts of Madagascar's forests past.

Phylogeographic analysis can be described as the study of the geological and climatological processes that have produced contemporary geographic distributions of populations and species. Here, we attempt to understand how the dynamic process of landscape change on Madagascar has shaped the distribution of a targeted clade of mouse lemurs (genus Microcebus) and, conversely, how phylogenetic and population genetic patterns in these small primates can reciprocally advance our understanding of Madagascar's prehuman environment. The degree to which human activity has impacted the natural plant communities of Madagascar is of critical and enduring interest. Today, the eastern rainforests are separated from the dry deciduous forests of the west by a large expanse of presumed anthropogenic grassland savanna, dominated by the Family Poaceae, that blankets most of the Central Highlands. Although there is firm consensus that anthropogenic activities have transformed the original vegetation through agricultural and pastoral practices, the degree to which closed-canopy forest extended from the east to the west remains debated. Phylogenetic and population genetic patterns in a five-species clade of mouse lemurs suggest that longitudinal dispersal across the island was readily achieved throughout the Pleistocene, apparently ending at ∼55 ka. By examining patterns of both inter- and intraspecific genetic diversity in mouse lemur species found in the eastern, western, and Central Highland zones, we conclude that the natural environment of the Central Highlands would have been mosaic, consisting of a matrix of wooded savanna that formed a transitional zone between the extremes of humid eastern and dry western forest types.

Hybrid de novo genome assembly and centromere characterization of the gray mouse lemur (Microcebus murinus).

BACKGROUND: The de novo assembly of repeat-rich mammalian genomes using only high-throughput short read sequencing data typically results in highly fragmented genome assemblies that limit downstream applications. Here, we present an iterative approach to hybrid de novo genome assembly that incorporates datasets stemming from multiple genomic technologies and methods. We used this approach to improve the gray mouse lemur (Microcebus murinus) genome from early draft status to a near chromosome-scale assembly. METHODS: We used a combination of advanced genomic technologies to iteratively resolve conflicts and super-scaffold the M. murinus genome. RESULTS: We improved the M. murinus genome assembly to a scaffold N50 of 93.32 Mb. Whole genome alignments between our primary super-scaffolds and 23 human chromosomes revealed patterns that are congruent with historical comparative cytogenetic data, thus demonstrating the accuracy of our de novo scaffolding approach and allowing assignment of scaffolds to M. murinus chromosomes. Moreover, we utilized our independent datasets to discover and characterize sequences associated with centromeres across the mouse lemur genome. Quality assessment of the final assembly found 96% of mouse lemur canonical transcripts nearly complete, comparable to other published high-quality reference genome assemblies. CONCLUSIONS: We describe a new assembly of the gray mouse lemur (Microcebus murinus) genome with chromosome-scale scaffolds produced using a hybrid bioinformatic and sequencing approach. The approach is cost effective and produces superior results based on metrics of contiguity and completeness. Our results show that emerging genomic technologies can be used in combination to characterize centromeres of non-model species and to produce accurate de novo chromosome-scale genome assemblies of complex mammalian genomes.

The Alu neurodegeneration hypothesis: A primate-specific mechanism for neuronal transcription noise, mitochondrial dysfunction, and manifestation of neurodegenerative disease.

It is hypothesized that retrotransposons have played a fundamental role in primate evolution and that enhanced neurologic retrotransposon activity in humans may underlie the origin of higher cognitive function. As a potential consequence of this enhanced activity, it is likely that neurons are susceptible to deleterious retrotransposon pathways that can disrupt mitochondrial function. An example is observed in the TOMM40 gene, encoding a ß-barrel protein critical for mitochondrial preprotein transport. Primate-specific Alu retrotransposons have repeatedly inserted into TOMM40 introns, and at least one variant associated with late-onset Alzheimer's disease originated from an Alu insertion event. We provide evidence of enriched Alu content in mitochondrial genes and postulate that Alus can disrupt mitochondrial populations in neurons, thereby setting the stage for progressive neurologic dysfunction. This Alu neurodegeneration hypothesis is compatible with decades of research and offers a plausible mechanism for the disruption of neuronal mitochondrial homeostasis, ultimately cascading into neurodegenerative disease.

Blood transcriptomes reveal novel parasitic zoonoses circulating in Madagascar's lemurs.

Zoonotic diseases are a looming threat to global populations, and nearly 75% of emerging infectious diseases can spread among wildlife, domestic animals and humans. A 'One World, One Health' perspective offers us an ideal framework for understanding and potentially mitigating the spread of zoonoses, and the island of Madagascar serves as a natural laboratory for conducting these studies. Rapid habitat degradation and climate change on the island are contributing to more frequent contact among humans, livestock and wildlife, increasing the potential for pathogen spillover events. Given Madagascar's long geographical isolation, coupled with recent and repeated introduction of agricultural and invasive species, it is likely that a number of circulating pathogens remain uncharacterized in lemur populations. Thus, it is imperative that new approaches be implemented for de novo pathogen discovery. To this end, we used non-targeted deep sequencing of blood transcriptomes from two species of critically endangered wild lemurs (Indri indri and Propithecus diadema) to characterize blood-borne pathogens. Our results show several undescribed vector-borne parasites circulating within lemurs, some of which may cause disease in wildlife, livestock and humans. We anticipate that advanced methods for de novo identification of unknown pathogens will have broad utility for characterizing other complex disease transmission systems.

Evaluating whole transcriptome amplification for gene profiling experiments using RNA-Seq.

BACKGROUND: RNA-Seq has enabled high-throughput gene expression profiling to provide insight into the functional link between genotype and phenotype. Low quantities of starting RNA can be a severe hindrance for studies that aim to utilize RNA-Seq. To mitigate this bottleneck, whole transcriptome amplification (WTA) technologies have been developed to generate sufficient sequencing targets from minute amounts of RNA. Successful WTA requires accurate replication of transcript abundance without the loss or distortion of specific mRNAs. Here, we test the efficacy of NuGEN's Ovation RNA-Seq V2 system, which uses linear isothermal amplification with a unique chimeric primer for amplification, using white adipose tissue from standard laboratory rats (Rattus norvegicus). Our goal was to investigate potential biological artifacts introduced through WTA approaches by establishing comparisons between matched raw and amplified RNA libraries derived from biological replicates. RESULTS: We found that 93% of expressed genes were identical between all unamplified versus matched amplified comparisons, also finding that gene density is similar across all comparisons. Our sequencing experiment and downstream bioinformatic analyses using the Tuxedo analysis pipeline resulted in the assembly of 25,543 high-quality transcripts. Libraries constructed from raw RNA and WTA samples averaged 15,298 and 15,253 expressed genes, respectively. Although significant differentially expressed genes (P < 0.05) were identified in all matched samples, each of these represents less than 0.15% of all shared genes for each comparison. CONCLUSIONS: Transcriptome amplification is efficient at maintaining relative transcript frequencies with no significant bias when using this NuGEN linear isothermal amplification kit under ideal laboratory conditions as presented in this study. This methodology has broad applications, from clinical and diagnostic, to field-based studies when sample acquisition, or sample preservation, methods prove challenging.

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms.

BACKGROUND: Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods. RESULTS: We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology. CONCLUSIONS: Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Prion forensics: a multidisciplinary approach to investigate CWD at an illegal deer carcass disposal site.

Infectious prions are resistant to degradation and remain infectious in the environment for several years. Chronic wasting disease (CWD) has been detected in cervids inhabiting North America, the Nordic countries, and South Korea. CWD-prion spread is partially attributed to carcass transport and disposal. We employed a forensic approach to investigate an illegal carcass dump site connected with a CWD-positive herd. We integrated anatomic, genetic, and prion amplification methods to discover CWD-positive remains from six white-tailed deer (Odocoileus virginianus) and, using microsatellite markers, confirmed a portion originated from the CWD-infected herd. This approach provides a foundation for future studies of carcass prion transmission risk.

Natural hybridization generates mammalian lineage with species characteristics.

Most diploid species arise from single-species ancestors. Hybrid origins of new species are uncommon (except among polyploids) and are documented infrequently in animals. Examples of natural hybridization leading to speciation in mammals are exceedingly rare. Here, we show a Caribbean species of bat (Artibeus schwartzi) has a nuclear genome derived from two nonsister but congeneric species (A. jamaicensis and A. planirostris) and a mitochondrial genome that is from a third extinct or uncharacterized congener. Artibeus schwartzi is self-sustaining, morphologically distinct, and exists in near geographic isolation of its known parent species. Island effects (i.e., area, reduced habitat variability, and geographic isolation) likely have restricted gene flow from parental species into the Caribbean populations of this hybrid lineage, thus contributing to local adaptation and isolation of this newly produced taxon. We hypothesize differential rates of the development of reproductive isolation within the genus and estimate that 2.5 million years was an insufficient amount of time for the development of postzygotic isolation among the three species that hybridized to produce A. schwartzi. Reticulated evolution thus has resulted in a genomic combination from three evolutionary lineages and a transgressive phenotype that is distinct from all other known species of Artibeus. The data herein further demonstrate the phenomenon of speciation by hybridization in mammals is possible in nature.