snpAIMeR: R package for evaluating ancestry informative marker contributions in non-model population diagnostics.

MOTIVATION: Single nucleotide polymorphism (SNP) markers are increasingly popular for population genomics and inferring ancestry for individuals of unknown origin. Because large SNP datasets are impractical for rapid and routine analysis, diagnostics rely on panels of highly informative markers. Strategies exist for selecting these markers, however, resources for efficiently evaluating their performance are limited for non-model systems. RESULTS: snpAIMeR is a user-friendly R package that evaluates the efficacy of genomic markers for the cluster assignment of unknown individuals. It is intended to help minimize panel size and genotyping effort by determining the informativeness of candidate diagnostic markers. Provided genotype data from individuals of known origin, it uses leave-one-out cross-validation to determine population assignment rates for individual markers and marker combinations. AVAILABILITY AND IMPLEMENTATION: snpAIMeR is available on CRAN (https://CRAN.R-project.org/package=snpAIMeR).

Handling Logical Character Dependency in Phylogenetic Inference: Extensive Performance Testing of Assumptions and Solutions Using Simulated and Empirical Data.

Logical character dependency is a major conceptual and methodological problem in phylogenetic inference of morphological data sets, as it violates the assumption of character independence that is common to all phylogenetic methods. It is more frequently observed in higher-level phylogenies or in data sets characterizing major evolutionary transitions, as these represent parts of the tree of life where (primary) anatomical characters either originate or disappear entirely. As a result, secondary traits related to these primary characters become "inapplicable" across all sampled taxa in which that character is absent. Various solutions have been explored over the last three decades to handle character dependency, such as alternative character coding schemes and, more recently, new algorithmic implementations. However, the accuracy of the proposed solutions, or the impact of character dependency across distinct optimality criteria, has never been directly tested using standard performance measures. Here, we utilize simple and complex simulated morphological data sets analyzed under different maximum parsimony optimization procedures and Bayesian inference to test the accuracy of various coding and algorithmic solutions to character dependency. This is complemented by empirical analyses using a recoded data set on palaeognathid birds. We find that in small, simulated data sets, absent coding performs better than other popular coding strategies available (contingent and multistate), whereas in more complex simulations (larger data sets controlled for different tree structure and character distribution models) contingent coding is favored more frequently. Under contingent coding, a recently proposed weighting algorithm produces the most accurate results for maximum parsimony. However, Bayesian inference outperforms all parsimony-based solutions to handle character dependency due to fundamental differences in their optimization procedures-a simple alternative that has been long overlooked. Yet, we show that the more primary characters bearing secondary (dependent) traits there are in a data set, the harder it is to estimate the true phylogenetic tree, regardless of the optimality criterion, owing to a considerable expansion of the tree parameter space. [Bayesian inference, character dependency, character coding, distance metrics, morphological phylogenetics, maximum parsimony, performance, phylogenetic accuracy.].

Toward transparent taxonomy: an interactive web-tool for evaluating competing taxonomic arrangements.

Informative and consistent taxonomy above the species level is essential to communication about evolution, biodiversity and conservation, and yet the practice of taxonomy is considered opaque and subjective by non-taxonomist scientists and the public alike. While various proposals have tried to make the basis for the ranking and inclusiveness of taxa more transparent and objective, widespread adoption of these ideas has lagged. Here, we present TaxonomR, an interactive online decision-support tool to evaluate alternative taxonomic classifications. This tool implements an approach that quantifies the criteria commonly used in taxonomic treatments and allows the user to interactively manipulate weightings for different criteria to compare scores for taxonomic groupings under those weights. We use the butterfly taxon Argynnis to demonstrate how different weightings applied to common taxonomic criteria result in fundamentally different genus-level classifications that are predominantly used in different continents and geographic regions. These differences are objectively compared and quantified using TaxonomR to evaluate the kinds of criteria that have been emphasized in earlier classifications, and the nature of the support for current alternative taxonomic arrangements. The main role of TaxonomR is to make taxonomic decisions transparent via an explicit prioritization scheme. TaxonomR is not a prescriptive application. Rather, it aims to be a tool for facilitating our understanding of alternative taxonomic classifications that can, in turn, potentially support global harmony in biodiversity assessments through evidence-based discussion and community-wide resolution of historically entrenched taxonomic tensions.

The origin of squamates revealed by a Middle Triassic lizard from the Italian Alps.

Modern squamates (lizards, snakes and amphisbaenians) are the world's most diverse group of tetrapods along with birds 1 and have a long evolutionary history, with the oldest known fossils dating from the Middle Jurassic period-168 million years ago2-4. The evolutionary origin of squamates is contentious because of several issues: (1) a fossil gap of approximately 70 million years exists between the oldest known fossils and their estimated origin5-7; (2) limited sampling of squamates in reptile phylogenies; and (3) conflicts between morphological and molecular hypotheses regarding the origin of crown squamates6,8,9. Here we shed light on these problems by using high-resolution microfocus X-ray computed tomography data from the articulated fossil reptile Megachirella wachtleri (Middle Triassic period, Italian Alps 10 ). We also present a phylogenetic dataset, combining fossils and extant taxa, and morphological and molecular data. We analysed this dataset under different optimality criteria to assess diapsid reptile relationships and the origins of squamates. Our results re-shape the diapsid phylogeny and present evidence that M. wachtleri is the oldest known stem squamate. Megachirella is 75 million years older than the previously known oldest squamate fossils, partially filling the fossil gap in the origin of lizards, and indicates a more gradual acquisition of squamatan features in diapsid evolution than previously thought. For the first time, to our knowledge, morphological and molecular data are in agreement regarding early squamate evolution, with geckoes-and not iguanians-as the earliest crown clade squamates. Divergence time estimates using relaxed combined morphological and molecular clocks show that lepidosaurs and most other diapsids originated before the Permian/Triassic extinction event, indicating that the Triassic was a period of radiation, not origin, for several diapsid lineages.

Gauging ages of tiger swallowtail butterflies using alternate SNP analyses.

Divergence times underpin diverse evolutionary hypotheses, but conflicting age estimates across studies diminish the validity of such hypotheses. These conflicts have continued to grow as large genomics datasets become commonplace and analytical approaches proliferate. To provide more stable temporal intervals, age estimations should be interpreted in the context of both the type of data and analysis being used. Here, we use multispecies coalescent (MSC), concatenation-based, and categorical data transformation approaches on genome-wide SNP data to infer divergence ages within the Papilio glaucus group of tiger swallowtail butterflies in North America. While the SNP data supported previously recognized relationships within the group (P. multicaudata, ((P. eurymedon, P. rutulus), (P. appalachiensis, P. canadensis, P. glaucus))), estimated ages of divergence between the major lineages varied substantially among analyses. MSC produced wide credibility intervals particularly for deeper nodes, reflecting uncertainty in the coalescence times as a possible result of conflicting signal across gene trees. Concatenation, in contrast, gave narrower and more well-defined posterior distributions for the node ages; however, the higher precision of these time estimates is a likely artefact due to more simplistic underlying assumptions of this approach that do not account for conflict among gene trees. Transformed categorical data analysis gave the least precise and the most variable results, with its simple substitution model coupled with a relaxed clock tending to produce spurious results from large genome-wide datasets. While median node ages differed considerably between analyses (∼2 Mya between MSC and concatenation-based results), their corresponding credibility intervals nonetheless highlight common temporal patterns for deeper divergences in the group as well as finer-scale phylogeography. Age distributions across analyses support an origin of the group during the warm period of the early to mid-Pliocene. Late Pliocene climate aridification and cooling drove divergence between eastern and western groups that further diversified during the period of repeated Pleistocene glaciations. Our results provide a structured comparative assessment of divergence time estimates and evolutionary relationships in a well-studied group of butterflies, and support better understanding of analytical biases in divergence time estimation.

Evaluating the Performance of Probabilistic Algorithms for Phylogenetic Analysis of Big Morphological Datasets: A Simulation Study.

Reconstructing the tree of life is an essential task in evolutionary biology. It demands accurate phylogenetic inference for both extant and extinct organisms, the latter being almost entirely dependent on morphological data. While parsimony methods have traditionally dominated the field of morphological phylogenetics, a rapidly growing number of studies are now employing probabilistic methods (maximum likelihood and Bayesian inference). The present-day toolkit of probabilistic methods offers varied software with distinct algorithms and assumptions for reaching global optimality. However, benchmark performance assessments of different software packages for the analyses of morphological data, particularly in the era of big data, are still lacking. Here, we test the performance of four major probabilistic software under variable taxonomic sampling and missing data conditions: the Bayesian inference-based programs MrBayes and RevBayes, and the maximum likelihood-based IQ-TREE and RAxML. We evaluated software performance by calculating the distance between inferred and true trees using a variety of metrics, including Robinson-Foulds (RF), Matching Splits (MS), and Kuhner-Felsenstein (KF) distances. Our results show that increased taxonomic sampling improves accuracy, precision, and resolution of reconstructed topologies across all tested probabilistic software applications and all levels of missing data. Under the RF metric, Bayesian inference applications were the most consistent, accurate, and robust to variation in taxonomic sampling in all tested conditions, especially at high levels of missing data, with little difference in performance between the two tested programs. The MS metric favored more resolved topologies that were generally produced by IQ-TREE. Adding more taxa dramatically reduced performance disparities between programs. Importantly, our results suggest that the RF metric penalizes incorrectly resolved nodes (false positives) more severely than the MS metric, which instead tends to penalize polytomies. If false positives are to be avoided in systematics, Bayesian inference should be preferred over maximum likelihood for the analysis of morphological data.

New insights into patterns and rates of tooth replacement in serrasalmid and characid fishes, with implications for the subsistence fishery of Peru's remote ribereños villages.

We explored patterns, rates and unexpected socio-ecological consequences of tooth replacement in serrasalmids and characids of the Peruvian Amazon using microcomputed tomography. Of 24 specimens collected in February 2019, representing a mix of red-bellied piranha Pygocentrus nattereri, redeye piranha Serrasalmus rhombeus, silver dollar fish Ctenobrycon hauxwellianus and mojara Astyanax abramis, six individuals possessed edentulous jaw quadrants. On average, 22.9% of fish collected per day from these species featured incomplete dentition, a value three to five times higher than anticipated based on replacement rates estimated from captive fish, differences that may be driven by ontogeny, seasonality or environmental quality.

Delimitation of Alosa species (Teleostei: Clupeiformes) from the Sea of Azov: integrating morphological and molecular approaches.

Shads of the genus Alosa are essential to commercial fisheries across North America and Europe, but in some areas their species boundaries remain controversial. Traditional morphology-based taxonomy of Alosa spp. has relied heavily on the number of gill rakers and body proportions, but these can be highly variable. We use mitochondrial (mt)DNA (coI and cytb) and genome-wide single nucleotide polymorphisms (SNP) along with morphological characters to assess differentiation among endemic Ponto-Caspian shads in the Sea of Azov. Morphological species assignments based on gill-raker number were not congruent with genetic lineages shown by mtDNA and SNPs. Iterative analysis revealed that genetic lineages were associated with sampling location and several other morphometric traits (caudal peduncle depth, pre-anal length and head length). Phylogenetic analysis of the genus placed Ponto-Caspian Alosa spp. in the same evolutionary lineage as endangered Alosa spp. endemic to Greece, highlighting the importance of these findings to conservation management. We conclude that gill-raker number is not reliable for delimiting species of Alosa. This taxonomic uncertainty should be addressed by examining type material to provide a robust integrative classification for these commercially important fishes.

Sequence signatures within the genome of SARS-CoV-2 can be used to predict host source.

We conducted an in silico analysis to better understand the potential factors impacting host adaptation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in white-tailed deer, humans, and mink due to the strong evidence of sustained transmission within these hosts. Classification models trained on single nucleotide and amino acid differences between samples effectively identified white-tailed deer-, human-, and mink-derived SARS-CoV-2. For example, the balanced accuracy score of Extremely Randomized Trees classifiers was 0.984 ± 0.006. Eighty-eight commonly identified predictive mutations are found at sites under strong positive and negative selective pressure. A large fraction of sites under selection (86.9%) or identified by machine learning (87.1%) are found in genes other than the spike. Some locations encoded by these gene regions are predicted to be B- and T-cell epitopes or are implicated in modulating the immune response suggesting that host adaptation may involve the evasion of the host immune system, modulation of the class-I major-histocompatibility complex, and the diminished recognition of immune epitopes by CD8+ T cells. Our selection and machine learning analysis also identified that silent mutations, such as C7303T and C9430T, play an important role in discriminating deer-derived samples across multiple clades. Finally, our investigation into the origin of the B.1.641 lineage from white-tailed deer in Canada discovered an additional human sequence from Michigan related to the B.1.641 lineage sampled near the emergence of this lineage. These findings demonstrate that machine-learning approaches can be used in combination with evolutionary genomics to identify factors possibly involved in the cross-species transmission of viruses and the emergence of novel viral lineages.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible virus capable of infecting and establishing itself in human and wildlife populations, such as white-tailed deer. This fact highlights the importance of developing novel ways to identify genetic factors that contribute to its spread and adaptation to new host species. This is especially important since these populations can serve as reservoirs that potentially facilitate the re-introduction of new variants into human populations. In this study, we apply machine learning and phylogenetic methods to uncover biomarkers of SARS-CoV-2 adaptation in mink and white-tailed deer. We find evidence demonstrating that both non-synonymous and silent mutations can be used to differentiate animal-derived sequences from human-derived ones and each other. This evidence also suggests that host adaptation involves the evasion of the immune system and the suppression of antigen presentation. Finally, the methods developed here are general and can be used to investigate host adaptation in viruses other than SARS-CoV-2.

Characterization of an African Swine Fever Virus Field Isolate from Vietnam with Deletions in the Left Variable Multigene Family Region.

In this paper, we report the characterization of a genetically modified live-attenuated African swine fever virus (ASFV) field strain isolated from Vietnam. The isolate, ASFV-GUS-Vietnam, belongs to p72 genotype II, has six multi-gene family (MGF) genes deleted, and an Escherichia coli GusA gene (GUS) inserted. When six 6-8-week-old pigs were inoculated with ASFV-GUS-Vietnam oro-nasally (2 × 105 TCID50/pig), they developed viremia, mild fever, lethargy, and inappetence, and shed the virus in their oral and nasal secretions and feces. One of the pigs developed severe clinical signs and was euthanized 12 days post-infection, while the remaining five pigs recovered. When ASFV-GUS-Vietnam was inoculated intramuscularly (2 × 103 TCID50/pig) into four 6-8 weeks old pigs, they also developed viremia, mild fever, lethargy, inappetence, and shed the virus in their oral and nasal secretions and feces. Two contact pigs housed together with the four intramuscularly inoculated pigs, started to develop fever, viremia, loss of appetite, and lethargy 12 days post-contact, confirming horizontal transmission of ASFV-GUS-Vietnam. One of the contact pigs died of ASF on day 23 post-contact, while the other one recovered. The pigs that survived the exposure to ASFV-GUS-Vietnam via the mucosal or parenteral route were fully protected against the highly virulent ASFV Georgia 2007/1 challenge. This study showed that ASFV-GUS-Vietnam field isolate is able to induce complete protection in the majority of the pigs against highly virulent homologous ASFV challenge, but has the potential for horizontal transmission, and can be fatal in some animals. This study highlights the need for proper monitoring and surveillance when ASFV live-attenuated virus-based vaccines are used in the field for ASF control in endemic countries.

HiMAP2: Identifying phylogenetically informative genetic markers from diverse genomic resources.

Multiplexed amplicon sequencing offers a cost-effective and rapid solution for phylogenomic studies that include a large number of individuals. Selecting informative genetic markers is a critical initial step in designing such multiplexed amplicon panels, but screening various genomic data and selecting markers that are informative for the question at hand can be laborious. Here, we present a flexible and user-friendly tool, HiMAP2, for identifying, visualizing and filtering phylogenetically informative loci from diverse genomic and transcriptomic resources. This bioinformatics pipeline includes orthology prediction, exon extraction and filtering of aligned exon sequences according to user-defined specifications. Additionally, HiMAP2 facilitates exploration of the final filtered exons by incorporating phylogenetic inference of individual exon trees with raxml-ng as well as the estimation of a species tree using astral. Finally, results of the marker selection can be visualized and refined with an interactive Bokeh application that can be used to generate publication-quality figures. Source code and user instructions for HiMAP2 are available at https://github.com/popphylotools/HiMAP_v2.

Genomic and transcriptomic characterization of delta SARS-CoV-2 infection in free-ranging white-tailed deer (Odocoileus virginianus).

White-tailed deer (WTD) are susceptible to SARS-CoV-2 and represent an important species for surveillance. Samples from WTD (n = 258) collected in November 2021 from Québec, Canada were analyzed for SARS-CoV-2 RNA. We employed viral genomics and host transcriptomics to further characterize infection and investigate host response. We detected Delta SARS-CoV-2 (B.1.617.2) in WTD from the Estrie region; sequences clustered with human sequences from October 2021 from Vermont, USA, which borders this region. Mutations in the S-gene and a deletion in ORF8 were detected. Host expression patterns in SARS-CoV-2 infected WTD were associated with the innate immune response, including signaling pathways related to anti-viral, pro- and anti-inflammatory signaling, and host damage. We found limited correlation between genes associated with innate immune response from human and WTD nasal samples, suggesting differences in responses to SARS-CoV-2 infection. Our findings provide preliminary insights into host response to SARS-CoV-2 infection in naturally infected WTD.

The Spruce Budworm Genome: Reconstructing the Evolutionary History of Antifreeze Proteins.

Insects have developed various adaptations to survive harsh winter conditions. Among freeze-intolerant species, some produce "antifreeze proteins" (AFPs) that bind to nascent ice crystals and inhibit further ice growth. Such is the case of the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae), a destructive North American conifer pest that can withstand temperatures below -30°C. Despite the potential importance of AFPs in the adaptive diversification of Choristoneura, genomic tools to explore their origins have until now been limited. Here we present a chromosome-scale genome assembly for C. fumiferana, which we used to conduct comparative genomic analyses aimed at reconstructing the evolutionary history of tortricid AFPs. The budworm genome features 16 genes homologous to previously reported C. fumiferana AFPs (CfAFPs), 15 of which map to a single region on chromosome 18. Fourteen of these were also detected in five congeneric species, indicating Choristoneura AFP diversification occurred before the speciation event that led to C. fumiferana. Although budworm AFPs were previously considered unique to the genus Choristoneura, a search for homologs targeting recently sequenced tortricid genomes identified seven CfAFP-like genes in the distantly related Notocelia uddmanniana. High structural similarity between Notocelia and Choristoneura AFPs suggests a common origin, despite the absence of homologs in three related tortricids. Interestingly, one Notocelia AFP formed the C-terminus of a "zonadhesin-like" protein, possibly representing the ancestral condition from which tortricid AFPs evolved. Future work should clarify the evolutionary path of AFPs between Notocelia and Choristoneura and assess the role of the "zonadhesin-like" protein as precursor of tortricid AFPs.

Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission.

Wildlife reservoirs of broad-host-range viruses have the potential to enable evolution of viral variants that can emerge to infect humans. In North America, there is phylogenomic evidence of continual transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from humans to white-tailed deer (Odocoileus virginianus) through unknown means, but no evidence of transmission from deer to humans. We carried out an observational surveillance study in Ontario, Canada during November and December 2021 (n = 300 deer) and identified a highly divergent lineage of SARS-CoV-2 in white-tailed deer (B.1.641). This lineage is one of the most divergent SARS-CoV-2 lineages identified so far, with 76 mutations (including 37 previously associated with non-human mammalian hosts). From a set of five complete and two partial deer-derived viral genomes we applied phylogenomic, recombination, selection and mutation spectrum analyses, which provided evidence for evolution and transmission in deer and a shared ancestry with mink-derived virus. Our analysis also revealed an epidemiologically linked human infection. Taken together, our findings provide evidence for sustained evolution of SARS-CoV-2 in white-tailed deer and of deer-to-human transmission.

Megaevolutionary dynamics and the timing of evolutionary innovation in reptiles.

The origin of phenotypic diversity among higher clades is one of the most fundamental topics in evolutionary biology. However, due to methodological challenges, few studies have assessed rates of evolution and phenotypic disparity across broad scales of time to understand the evolutionary dynamics behind the origin and early evolution of new clades. Here, we provide a total-evidence dating approach to this problem in diapsid reptiles. We find major chronological gaps between periods of high evolutionary rates (phenotypic and molecular) and expansion in phenotypic disparity in reptile evolution. Importantly, many instances of accelerated phenotypic evolution are detected at the origin of major clades and body plans, but not concurrent with previously proposed periods of adaptive radiation. Furthermore, strongly heterogenic rates of evolution mark the acquisition of similarly adapted functional types, and the origin of snakes is marked by the highest rates of phenotypic evolution in diapsid history.

X-ray computed microtomography of Megachirella wachtleri.

Understanding the origin and early evolution of squamates has been a considerable challenge given the extremely scarce fossil record of early squamates and their poor degree of preservation. In order to overcome those limitations, we conducted high-resolution X-ray computed tomography (CT) studies on the fossil reptile Megachirella wachtleri (Middle Triassic, northern Italy), which revealed an important set of features indicating this is the oldest known fossil squamate in the world, predating the previous oldest record by ca. 75 million years. We also compiled a new phylogenetic data set comprising a large sample of diapsid reptiles (including morphological and molecular data) to investigate the phylogenetic relationships of early squamates and other reptile groups along with the divergence time of those lineages. The re-description of Megachirella and a new phylogenetic hypothesis of diapsid relationships are presented in a separate study. Here we present the data descriptors for the tomographic scans of Megachirella, which holds fundamental information to our understanding on the early evolution of one of the largest vertebrate groups on Earth today.

Mosasauroid phylogeny under multiple phylogenetic methods provides new insights on the evolution of aquatic adaptations in the group.

Mosasauroids were a successful lineage of squamate reptiles (lizards and snakes) that radiated during the Late Cretaceous (95-66 million years ago). They can be considered one of the few lineages in the evolutionary history of tetrapods to have acquired a fully aquatic lifestyle, similarly to whales, ichthyosaurs and plesiosaurs. Despite a long history of research on this group, their phylogenetic relationships have only been tested so far using traditional (unweighted) maximum parsimony. However, hypotheses of mosasauroid relationships and the recently proposed multiple origins of aquatically adapted pelvic and pedal features in this group can be more thoroughly tested by methods that take into account variation in branch lengths and evolutionary rates. In this study, we present the first mosasauroid phylogenetic analysis performed under different analytical methods, including maximum likelihood, Bayesian inference, and implied weighting maximum parsimony. The results indicate a lack of congruence in the topological position of halisaurines and Dallasaurus. Additionally, the genus Prognathodon is paraphyletic under all hypotheses. Interestingly, a number of traditional mosasauroid clades become weakly supported, or unresolved, under Bayesian analyses. The reduced resolutions in some consensus trees create ambiguities concerning the evolution of fully aquatic pelvic/pedal conditions under many analyses. However, when enough resolution was obtained, reversals of the pelvic/pedal conditions were favoured by parsimony and likelihood ancestral state reconstructions instead of independent origins of aquatic features in mosasauroids. It is concluded that most of the observed discrepancies among the results can be associated with different analytical procedures, but also due to limited postcranial data on halisaurines, yaguarasaurines and Dallasaurus.