Are rabid raccoons (Procyon lotor) ready for the rapture? Determining the geographic origin of rabies virus-infected raccoons using RADcapture and microhaplotypes.

North America is recognized for the exceptional richness of rabies virus (RV) wildlife reservoir species. Management of RV is accomplished through vaccination targeting mesocarnivore reservoir populations, such as the raccoon (Procyon lotor) in Eastern North America. Raccoons are a common generalist species, and populations may reach high densities in developed areas, which can result in contact with humans and pets with potential exposures to the raccoon variant of RV throughout the eastern United States. Understanding the spatial movement of RV by raccoon populations is important for monitoring and refining strategies supporting the landscape-level control and local elimination of this lethal zoonosis. We developed a high-throughput genotyping panel for raccoons based on hundreds of microhaplotypes to identify population structure and genetic diversity relevant to rabies management programs. Throughout the eastern United States, we identified hierarchical population genetic structure with clusters that were connected through isolation-by-distance. We also illustrate that this genotyping approach can be used to support real-time management priorities by identifying the geographic origin of a rabid raccoon that was collected in an area of the United States that had been raccoon RV-free for 8 years. The results from this study and the utility of the microhaplotype panel and genotyping method will provide managers with information on raccoon ecology that can be incorporated into future management decisions.

Old specimens for old branches: Assessing effects of sample age in resolving a rapid Neotropical radiation of squirrels.

Ultraconserved Elements (UCEs) have been useful to resolve challenging phylogenies of non-model clades, unpuzzling long-conflicted relationships in key branches of the Tree of Life at both deep and shallow levels. UCEs are often reliably recovered from historical samples, unlocking a vast number of preserved natural history specimens for analysis. However, the extent to which sample age and preservation method impact UCE recovery as well as downstream inferences remains unclear. Furthermore, there is an ongoing debate on how to curate, filter, and properly analyze UCE data when locus recovery is uneven across sample age and quality. In the present study we address these questions with an empirical dataset composed of over 3800 UCE loci from 219 historical and modern samples of Sciuridae, a globally distributed and ecologically important family of rodents. We provide a genome-scale phylogeny of two squirrel subfamilies (Sciurillinae and Sciurinae: Sciurini) and investigate their placement within Sciuridae. For historical specimens, recovery of UCE loci and mean length per locus were inversely related to sample age; deeper sequencing improved the number of UCE loci recovered but not locus length. Most of our phylogenetic inferences-performed on six datasets with alternative data-filtering strategies, and using three distinct optimality criteria-resulted in distinct topologies. Datasets containing more loci (40% and 50% taxa representativeness matrices) yielded more concordant topologies and higher support values than strictly filtered datasets (60% matrices) particularly with IQ-Tree and SVDquartets, while filtering based on information content provided better topological resolution for inferences with the coalescent gene-tree based approach in ASTRAL-III. We resolved deep relationships in Sciuridae (including among the five currently recognized subfamilies) and relationships among the deepest branches of Sciurini, but conflicting relationships remain at both genus- and species-levels for the rapid Neotropical tree squirrel radiation. Our results suggest that phylogenomic consensus can be difficult and heavily influenced by the age of available samples and the filtering steps used to optimize dataset properties.

Comparative Analysis of Annotation Pipelines Using the First Japanese White-Eye (Zosterops japonicus) Genome.

Introduced into Hawaii in the early 1900s, the Japanese white-eye or warbling white-eye (Zosterops japonicus) is now the most abundant land bird in the archipelago. Here, we present the first Z. japonicus genome, sequenced from an individual in its invasive range. This genome provides an important resource for future studies in invasion genomics. We annotated the genome using two workflows-standalone AUGUSTUS and BRAKER2. We found that AUGUSTUS was more conservative with gene predictions when compared with BRAKER2. The final number of annotated gene models was similar between the two workflows, but standalone AUGUSTUS had over 70% of gene predictions with Blast2GO annotations versus under 30% using BRAKER2. Additionally, we tested whether using RNA-seq data from 47 samples had a significant impact on annotation quality when compared with data from a single sample, as generating RNA-seq data for genome annotation can be expensive and requires well preserved tissue. We found that more data did not significantly change the number of annotated genes using AUGUSTUS but using BRAKER2 the number increased substantially. The results presented here will aid researchers in annotating draft genomes of nonmodel species as well as those studying invasion success.

Whole-Genome Sequencing of Procyonids Reveals Distinct Demographic Histories in Kinkajou (Potos flavus) and Northern Raccoon (Procyon lotor).

Here, we present the initial comparison of the nuclear genomes of the North American raccoon (Procyon lotor) and the kinkajou (Potos flavus) based on draft assemblies. These two species encompass almost 21 Myr of evolutionary history within Procyonidae. Because assemblies greatly impact downstream results, such as gene prediction and annotation, we tested three de novo assembly strategies (implemented in ALLPATHS-LG, MaSuRCA, and Platanus), some of which are optimized for highly heterozygous genomes. We discovered significant variation in contig and scaffold N50 and L50 statistics and genome completeness depending on the de novo assembler used. We compared the performance of these three assembly algorithms in hopes that this study will aid others looking to improve the quality of existing draft genome assemblies even without additional sequence data. We also estimate the demographic histories of raccoons and kinkajous using the Pairwise Sequentially Markovian Coalescent and discuss the variation in population sizes with respect to climatic change during the Pleistocene, as well as aspects of their ecology and taxonomy. Our goal is to achieve a better understanding of the evolutionary history of procyonids and to create robust genomic resources for future studies regarding adaptive divergence and selection.

Museomics of tree squirrels: a dense taxon sampling of mitogenomes reveals hidden diversity, phenotypic convergence, and the need of a taxonomic overhaul.

BACKGROUND: Tree squirrels (Sciuridae, Sciurini), in particular the highly diverse Neotropical lineages, are amongst the most rapidly diversifying branches of the mammal tree of life but also some of the least known. Negligence of this group by systematists is likely a product of the difficulties in assessing morphological informative traits and of the scarcity or unavailability of fresh tissue samples for DNA sequencing. The highly discrepant taxonomic arrangements are a consequence of the lack of phylogenies and the exclusive phenotypic-based classifications, which can be misleading in a group with conservative morphology. Here we used high-throughput sequencing and an unprecedented sampling of museum specimens to provide the first comprehensive phylogeny of tree squirrels, with a special emphasis on Neotropical taxa. RESULTS: We obtained complete or partial mitochondrial genomes from 232 historical and modern samples, representing 40 of the 43 currently recognized species of Sciurini. Our phylogenetic analyses-performed with datasets differing on levels of missing data and taxa under distinct analytical methods-strongly support the monophyly of Sciurini and consistently recovered 12 major clades within the tribe. We found evidence that the diversity of Neotropical tree squirrels is underestimated, with at least six lineages that represent taxa to be named or revalidated. Ancestral state reconstructions of number of upper premolars and number of mammae indicated that alternative conditions of both characters must have evolved multiple times throughout the evolutionary history of tree squirrels. CONCLUSIONS: Complete mitogenomes were obtained from museum specimens as old as 120 years, reinforcing the potential of historical samples for phylogenetic inferences of elusive lineages of the tree of life. None of the taxonomic arrangements ever proposed for tree squirrels fully corresponded to our phylogenetic reconstruction, with only a few of the currently recognized genera recovered as monophyletic. By investigating the evolution of two morphological traits widely employed in the taxonomy of the group, we revealed that their homoplastic nature can help explain the incongruence between phylogenetic results and the classification schemes presented so far. Based on our phylogenetic results we suggest a tentative supraspecific taxonomic arrangement for Sciurini, employing 13 generic names used in previous taxonomic classifications.

First Genome Sequence of the Gunnison's Prairie Dog (Cynomys gunnisoni), a Keystone Species and Player in the Transmission of Sylvatic Plague.

Prairie dogs (genus Cynomys) are a charismatic symbol of the American West. Their large social aggregations and complex vocalizations have been the subject of scientific and popular interest for decades. A large body of literature has documented their role as keystone species of western North America's grasslands: They generate habitat for other vertebrates, increase nutrient availability for plants, and act as a food source for mammalian, squamate, and avian predators. An additional keystone role lies in their extreme susceptibility to sylvatic plague (caused by Yersinia pestis), which results in periodic population extinctions, thereby generating spatiotemporal heterogeneity in both biotic communities and ecological processes. Here, we report the first Cynomys genome for a Gunnison's prairie dog (C. gunnisoni gunnisoni) from Telluride, Colorado (USA). The genome was constructed using a hybrid assembly of PacBio and Illumina reads and assembled with MaSuRCA and PBJelly, which resulted in a scaffold N50 of 824 kb. Total genome size was 2.67 Gb, with 32.46% of the bases occurring in repeat regions. We recovered 94.9% (91% complete) of the single copy orthologs using the mammalian Benchmarking Universal Single-Copy Orthologs database and detected 49,377 gene models (332,141 coding regions). Pairwise Sequentially Markovian Coalescent showed support for long-term stable population size followed by a steady decline beginning near the end of the Pleistocene, as well as a recent population reduction. The genome will aid in studies of mammalian evolution, disease resistance, and the genomic basis of life history traits in ground squirrels.

Phylogeny, biogeography and systematic revision of plain long-nosed squirrels (genus Dremomys, Nannosciurinae).

The plain long-nosed squirrels, genus Dremomys, are high elevation species in East and Southeast Asia. Here we present a complete molecular phylogeny for the genus based on nuclear and mitochondrial DNA sequences. Concatenated mitochondrial and nuclear gene trees were constructed to determine the tree topology, and date the tree. All speciation events within the plain-long nosed squirrels (genus Dremomys) were ancient (dated to the Pliocene or Miocene), and averaged older than many speciation events in the related Sunda squirrels, genus Sundasciurus. Within the plain long-nosed squirrels, the most recent interspecific split occurred 2.9 million years ago, older than some splits within Sunda squirrels, which dated to the Pleistocene. Our results demonstrate that the plain long-nosed squirrels are not monophyletic. The single species with a distinct distribution, the Bornean mountain ground squirrel (Dremomys everetti), which is endemic to the high mountains of Borneo, is nested within the Sunda squirrels with high support. This species diverged from its sister taxa in the Sunda squirrels 6.62 million years ago, and other plain long-nosed squirrels over 11 million years ago. Our analyses of morphological traits in these related genera support the re-classification of the Bornean mountain ground squirrel, Dremomys everetti, to the genus Sundasciurus, which changes its name to Sundasciurus everetti. Past inclusion in the plain long-nosed squirrels (Dremomys) reflects convergent evolution between these high elevation species.

In-solution hybridization for mammalian mitogenome enrichment: pros, cons and challenges associated with multiplexing degraded DNA.

Here, we present a set of RNA-based probes for whole mitochondrial genome in-solution enrichment, targeting a diversity of mammalian mitogenomes. This probes set was designed from seven mammalian orders and tested to determine the utility for enriching degraded DNA. We generated 63 mitogenomes representing five orders and 22 genera of mammals that yielded varying coverage ranging from 0 to >5400X. Based on a threshold of 70% mitogenome recovery and at least 10× average coverage, 32 individuals or 51% of samples were considered successful. The estimated sequence divergence of samples from the probe sequences used to construct the array ranged up to nearly 20%. Sample type was more predictive of mitogenome recovery than sample age. The proportion of reads from each individual in multiplexed enrichments was highly skewed, with each pool having one sample that yielded a majority of the reads. Recovery across each mitochondrial gene varied with most samples exhibiting regions with gaps or ambiguous sites. We estimated the ability of the probes to capture mitogenomes from a diversity of mammalian taxa not included here by performing a clustering analysis of published sequences for 100 taxa representing most mammalian orders. Our study demonstrates that a general array can be cost and time effective when there is a need to screen a modest number of individuals from a variety of taxa. We also address the practical concerns for using such a tool, with regard to pooling samples, generating high quality mitogenomes and detail a pipeline to remove chimeric molecules.

Taxonomic revision of the olingos (Bassaricyon), with description of a new species, the Olinguito.

We present the first comprehensive taxonomic revision and review the biology of the olingos, the endemic Neotropical procyonid genus Bassaricyon, based on most specimens available in museums, and with data derived from anatomy, morphometrics, mitochondrial and nuclear DNA, field observations, and geographic range modeling. Species of Bassaricyon are primarily forest-living, arboreal, nocturnal, frugivorous, and solitary, and have one young at a time. We demonstrate that four olingo species can be recognized, including a Central American species (Bassaricyon gabbii), lowland species with eastern, cis-Andean (Bassaricyon alleni) and western, trans-Andean (Bassaricyon medius) distributions, and a species endemic to cloud forests in the Andes. The oldest evolutionary divergence in the genus is between this last species, endemic to the Andes of Colombia and Ecuador, and all other species, which occur in lower elevation habitats. Surprisingly, this Andean endemic species, which we call the Olinguito, has never been previously described; it represents a new species in the order Carnivora and is the smallest living member of the family Procyonidae. We report on the biology of this new species based on information from museum specimens, niche modeling, and fieldwork in western Ecuador, and describe four Olinguito subspecies based on morphological distinctions across different regions of the Northern Andes.