Comparative Performance of Popular Methods for Hybrid Detection using Genomic Data.

Interspecific hybridization is an important evolutionary phenomenon that generates genetic variability in a population and fosters species diversity in nature. The availability of large genome scale data sets has revolutionized hybridization studies to shift from the observation of the presence or absence of hybrids to the investigation of the genomic constitution of hybrids and their genome-specific evolutionary dynamics. Although a handful of methods have been proposed in an attempt to identify hybrids, accurate detection of hybridization from genomic data remains a challenging task. In addition to methods that infer phylogenetic networks or that utilize pairwise divergence, site pattern frequency based and population genetic clustering approaches are popularly used in practice, though the performance of these methods under different hybridization scenarios has not been extensively examined. Here, we use simulated data to comparatively evaluate the performance of four tools that are commonly used to infer hybridization events: the site pattern frequency based methods HyDe and the $D$-statistic (i.e., the ABBA-BABA test) and the population clustering approaches structure and ADMIXTURE. We consider single hybridization scenarios that vary in the time of hybridization and the amount of incomplete lineage sorting (ILS) for different proportions of parental contributions ($\gamma$); introgressive hybridization; multiple hybridization scenarios; and a mixture of ancestral and recent hybridization scenarios. We focus on the statistical power to detect hybridization and the false discovery rate (FDR) for comparisons of the $D$-statistic and HyDe, and the accuracy of the estimates of $\gamma$ as measured by the mean squared error for HyDe, structure, and ADMIXTURE. Both HyDe and the $D$-statistic are powerful for detecting hybridization in all scenarios except those with high ILS, although the $D$-statistic often has an unacceptably high FDR. The estimates of $\gamma$ in HyDe are impressively robust and accurate whereas structure and ADMIXTURE sometimes fail to identify hybrids, particularly when the proportional parental contributions are asymmetric (i.e., when $\gamma$ is close to 0). Moreover, the posterior distribution estimated using structure exhibits multimodality in many scenarios, making interpretation difficult. Our results provide guidance in selecting appropriate methods for identifying hybrid populations from genomic data. [ABBA-BABA test; ADMIXTURE; hybridization; HyDe; introgression; Patterson's $D$-statistic; Structure.].

Classes of explicit phylogenetic networks and their biological and mathematical significance.

The evolutionary relationships among organisms have traditionally been represented using rooted phylogenetic trees. However, due to reticulate processes such as hybridization or lateral gene transfer, evolution cannot always be adequately represented by a phylogenetic tree, and rooted phylogenetic networks that describe such complex processes have been introduced as a generalization of rooted phylogenetic trees. In fact, estimating rooted phylogenetic networks from genomic sequence data and analyzing their structural properties is one of the most important tasks in contemporary phylogenetics. Over the last two decades, several subclasses of rooted phylogenetic networks (characterized by certain structural constraints) have been introduced in the literature, either to model specific biological phenomena or to enable tractable mathematical and computational analyses. In the present manuscript, we provide a thorough review of these network classes, as well as provide a biological interpretation of the structural constraints underlying these networks where possible. In addition, we discuss how imposing structural constraints on the network topology can be used to address the scalability and identifiability challenges faced in the estimation of phylogenetic networks from empirical data.

On the use of median-joining networks in evolutionary biology.

Median-joining (MJ) was proposed as a method for phylogeographical analysis and is enjoying increasing popularity. Herein, we evaluate the efficacy of the approach as originally intended. We show that median-joining networks (MJNs) are theoretically untenable for evolutionary inference, and that confusion has afflicted their use for over 15 years. The approach has two obvious shortcomings: its reliance on distance-based phenetics (overall similarity instead of character transformations) and the lack of rooting (no direction or history). Given that evolution involves both change and time, and the absence of rooting removes time (ancestor-descendant relationships) from the equation, the approach cannot yield defensible evolutionary interpretations. We also examine the impact of MJ analyses on evolutionary biology via an analysis of citations and conclude that the spread of MJNs through the literature is difficult to explain, especially given the availability of character-based analyses.

Identification of major histocompatibility complex genotypes associated with resistance to an amphibian emerging infectious disease.

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90ß1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIß1 sequences. The G90ß1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

There and Back Again: The Unexpected Journeys of Metridium de Blainville, 1824 between the Old Oceans and throughout the Modern World.

AbstractMembers of the sea anemone genus Metridium are abundant in temperate rocky habitats and fouling communities. Their biogeographic history is expected to reflect changes in currents and habitats that have influenced benthic communities, such as the climate-influenced changes that occurred during the Last Glacial Maximum. More recently, however, anthropogenic influences such as shipping transportation and the creation of artificial habitat have altered and affected the composition of modern-day marine communities. Here we use sequence-capture data to examine the genetic structure of Metridium across its shallow-water distribution to (1) evaluate species boundaries within Metridium, (2) elucidate the dispersal history of Metridium between and among oceans, and (3) assess the influence of anthropogenic movement on modern-day populations. We find strong evidence for two species within Metridium: M. farcimen and M. senile. Dispersal from the Pacific to the Atlantic included a subsequent isolation of a small population in or above the Bering Sea, which has presumably moved southward. Within the native range of M. senile, admixture is prevalent even between oceans as a result of anthropogenic activities. The nonnative populations in Chile and the Falkland Islands came from at least two distinct introduction events originating from both coasts of the United States in the North Pacific and North Atlantic Oceans. Hybridization between M. senile and M. farcimen is documented as occurring in anthropogenically influenced habitats. The heavy influence from anthropogenic activities will continue to impact our understanding of marine organisms, particularly within the native range and for those that are easily transported across long distances.

A time-calibrated mitogenomic phylogeny suggests that Korean Hyalessa fuscata is a bridge between Chinese and Japanese H.maculaticollis.

The cicada species, Hyalessa fuscata and H. maculaticollis(Hemiptera: Cicadidae), share numerous morphological characters, and their status as distinct species remains controversial. We reconstructed a phylogeny based on two new mitogenomes of H. fuscata from Korea and H. maculaticollis from Japan, in combination with GenBank sequences of H. maculaticollis from China and Japan, and other closely related cicada species. Maximum likelihood and Bayesian inference phylogenies showed that H. fuscata from Korea is more closely related to H. maculaticollis from China than either is to H. maculaticollis from Japan. The time-calibrated Bayesian evolutionary analysis by sampling trees (BEAST) phylogeny indicated that the mainland and insular forms diverged approximately 1.7-2.6 million years ago. This coincides with the formation of the East China Sea land bridge between East Asia and the Japanese archipelago, which would provide a dispersal corridor for Hyalessa from the mainland via the Korean peninsula southeastward to Japan. East Asian H. fuscata is a geographic variant that may be considered synonymous with H. maculaticollis.

Digest: Frequent hybridization in Darevskia rarely leads to the evolution of asexuality.

Asexuality in vertebrates is often generated via hybridization, but is it a rare product of pervasive hybridization or a common product of rare hybridization? Freitas et al. show that hybridization is frequent among the sexual species of Darevskia, although the crossings between parents of the asexual hybrids are undetected. This study illustrates that hybridization is not extraordinary in nature, and thus scalable phylogenetic network inference methods, rather than phylogenetic trees, are needed to accurately represent the true evolutionary history.