Excluding Loci With Substitution Saturation Improves Inferences From Phylogenomic Data.

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.].

Genomic Adaptations and Evolutionary History of the Extinct Scimitar-Toothed Cat, Homotherium latidens.

Homotherium was a genus of large-bodied scimitar-toothed cats, morphologically distinct from any extant felid species, that went extinct at the end of the Pleistocene [1-4]. They possessed large, saber-form serrated canine teeth, powerful forelimbs, a sloping back, and an enlarged optic bulb, all of which were key characteristics for predation on Pleistocene megafauna [5]. Previous mitochondrial DNA phylogenies suggested that it was a highly divergent sister lineage to all extant cat species [6-8]. However, mitochondrial phylogenies can be misled by hybridization [9], incomplete lineage sorting (ILS), or sex-biased dispersal patterns [10], which might be especially relevant for Homotherium since widespread mito-nuclear discrepancies have been uncovered in modern cats [10]. To examine the evolutionary history of Homotherium, we generated a ∼7x nuclear genome and a ∼38x exome from H. latidens using shotgun and target-capture sequencing approaches. Phylogenetic analyses reveal Homotherium as highly divergent (∼22.5 Ma) from living cat species, with no detectable signs of gene flow. Comparative genomic analyses found signatures of positive selection in several genes, including those involved in vision, cognitive function, and energy consumption, putatively consistent with diurnal activity, well-developed social behavior, and cursorial hunting [5]. Finally, we uncover relatively high levels of genetic diversity, suggesting that Homotherium may have been more abundant than the limited fossil record suggests [3, 4, 11-14]. Our findings complement and extend previous inferences from both the fossil record and initial molecular studies, enhancing our understanding of the evolution and ecology of this remarkable lineage.

A practical introduction to sequentially Markovian coalescent methods for estimating demographic history from genomic data.

A common goal of population genomics and molecular ecology is to reconstruct the demographic history of a species of interest. A pair of powerful tools based on the sequentially Markovian coalescent have been developed to infer past population sizes using genome sequences. These methods are most useful when sequences are available for only a limited number of genomes and when the aim is to study ancient demographic events. The results of these analyses can be difficult to interpret accurately, because doing so requires some understanding of their theoretical basis and of their sensitivity to confounding factors. In this practical review, we explain some of the key concepts underpinning the pairwise and multiple sequentially Markovian coalescent methods (PSMC and MSMC, respectively). We relate these concepts to the use and interpretation of these methods, and we explain how the choice of different parameter values by the user can affect the accuracy and precision of the inferences. Based on our survey of 100 PSMC studies and 30 MSMC studies, we describe how the two methods are used in practice. Readers of this article will become familiar with the principles, practice, and interpretation of the sequentially Markovian coalescent for inferring demographic history.

Expediting the sampling, decalcification, and forensic DNA analysis of large elephant ivory seizures to aid investigations and prosecutions.

The illegal ivory trade continues to drive elephant poaching. Large ivory seizures in Africa and Asia are still commonplace. Wildlife forensics is recognised as a key enforcement tool to combat this trade. However, the time and resources required to effectively test large ivory seizures is often prohibitive. This limits or delays testing, which may impede investigations and/or prosecutions. Typically, DNA analysis of an ivory seizure involves pairing and sorting the tusks, sampling the tusks, powdering the sample, decalcification, then DNA extraction. Here, we optimize the most time-consuming components of this process: sampling and decalcification. Firstly, using simulations, we demonstrate that tusks do not need to be paired to ensure an adequate number of unique elephants are sampled in a large seizure. Secondly, we determined that directly powdering the ivory using a Dremel drill with a high-speed cutter bit, instead of cutting the ivory with a circular saw and subsequently powdering the sample in liquid nitrogen with a freezer mill, produces comparable results. Finally, we optimized a rapid 2 -h decalcification protocol that produces comparable results to a standard 3-day protocol. We tested/optimised the protocols on 33 raw and worked ivory samples, and demonstrated their utility on a case study, successfully identifying 94% of samples taken from 123 tusks. Using these new rapid protocols, the entire sampling and DNA extraction process takes less than one day and requires less-expensive equipment. We expect that the implementation of these rapid protocols will promote more consistent and timely testing of ivory seizures suitable for enforcement action.

Direct transmission by injection affects competition among RNA viruses in honeybees.

The arrival of the ectoparasitic mite Varroa destructor on the western honeybee Apis mellifera saw a change in the diversity and prevalence of honeybee RNA viruses. One virus in particular, deformed wing virus (DWV) has become closely associated with V. destructor, leading many to conclude that V. destructor has affected viral virulence by changing the mode of transmission. While DWV is normally transmitted via feeding and faeces, V. destructor transmits viruses by direct injection. This change could have resulted in higher viral prevalence causing increased damage to the bees. Here we test the effect of a change in the mode of transmission on the composition and levels of honeybee RNA viruses in the absence of V. destructor. We find a rapid increase in levels of two viruses, sacbrood virus (SBV) and black queen cell virus (BQCV) after direct injection of viral extracts into honeybee pupae. In pupae injected with high levels of DWV extracted from symptomatic adult bees, DWV levels rapidly decline in the presence of SBV and BQCV. Further, we observe high mortality in honeybee pupae when injected with SBV and BQCV, whereas injecting pupae with high levels of DWV results in near 100% survival. Our results suggest a different explanation for the observed association between V. destructor and DWV. Instead of V. destructor causing an increase in DWV virulence, we hypothesize that direct virus inoculation, such as that mediated by a vector, quickly eliminates the most virulent honeybee viruses resulting in an association with less virulent viruses such as DWV.