Multiple origins of green coloration in frogs mediated by a novel biliverdin-binding serpin.

Many vertebrates have distinctive blue-green bones and other tissues due to unusually high biliverdin concentrations-a phenomenon called chlorosis. Despite its prevalence, the biochemical basis, biology, and evolution of chlorosis are poorly understood. In this study, we show that the occurrence of high biliverdin in anurans (frogs and toads) has evolved multiple times during their evolutionary history, and relies on the same mechanism-the presence of a class of serpin family proteins that bind biliverdin. Using a diverse combination of techniques, we purified these serpins from several species of nonmodel treefrogs and developed a pipeline that allowed us to assemble their complete amino acid and nucleotide sequences. The described proteins, hereafter named biliverdin-binding serpins (BBS), have absorption spectra that mimic those of phytochromes and bacteriophytochromes. Our models showed that physiological concentration of BBSs fine-tune the color of the animals, providing the physiological basis for crypsis in green foliage even under near-infrared light. Additionally, we found that these BBSs are most similar to human glycoprotein alpha-1-antitrypsin, but with a remarkable functional diversification. Our results present molecular and functional evidence of recurrent evolution of chlorosis, describe a biliverdin-binding protein in vertebrates, and introduce a function for a member of the serpin superfamily, the largest and most ubiquitous group of protease inhibitors.

Give and take: Effects of genetic admixture on mutation load in endangered Florida panthers.

Genetic admixture is a biological event inherent to genetic rescue programs aimed at the long-term conservation of endangered wildlife. Although the success of such programs can be measured by the increase in genetic diversity and fitness of subsequent admixed individuals, predictions supporting admixture costs to fitness due to the introduction of novel deleterious alleles are necessary. Here, we analyzed nonsynonymous variation from conserved genes to quantify and compare levels of mutation load (i.e. proportion of deleterious alleles and genotypes carrying these alleles) among endangered Florida panthers and non-endangered Texas pumas. Specifically, we used canonical (i.e. non-admixed) Florida panthers, Texas pumas, and F1 (canonical Florida × Texas) panthers dating from a genetic rescue program and Everglades National Park panthers with Central American ancestry resulting from an earlier admixture event. We found neither genetic drift nor selection significantly reduced overall proportions of deleterious alleles in the severely bottlenecked canonical Florida panthers. Nevertheless, the deleterious alleles identified were distributed into a disproportionately high number of homozygous genotypes due to close inbreeding in this group. Conversely, admixed Florida panthers (either with Texas or Central American ancestry) presented reduced levels of homozygous genotypes carrying deleterious alleles but increased levels of heterozygous genotypes carrying these variants relative to canonical Florida panthers. Although admixture is likely to alleviate the load of standing deleterious variation present in homozygous genotypes, our results suggest that introduced novel deleterious alleles (temporarily present in heterozygous state) in genetically rescued populations could potentially be expressed in subsequent generations if their effective sizes remain small.

Pulse magnetization elicits differential gene expression in the central nervous system of the Caribbean spiny lobster, Panulirus argus.

Diverse animals use Earth's magnetic field to guide their movements, but the neural and molecular mechanisms underlying the magnetic sense remain enigmatic. One hypothesis is that particles of the mineral magnetite (Fe3O4) provide the basis of magnetoreception. Here we examined gene expression in the central nervous system of a magnetically sensitive invertebrate, the Caribbean spiny lobster (Panulirus argus), after applying a magnetic pulse known to alter magnetic orientation behavior. Numerous genes were differentially expressed in response to the pulse, including 647 in the brain, 1256 in the subesophageal ganglion, and 712 in the thoracic ganglia. Many such genes encode proteins linked to iron regulation, oxidative stress, and immune response, consistent with possible impacts of a magnetic pulse on magnetite-based magnetoreceptors. Additionally, however, altered expression also occurred for numerous genes with no apparent link to magnetoreception, including genes encoding proteins linked to photoreception, carbohydrate and hormone metabolism, and other physiological processes. Overall, the results are consistent with the magnetite hypothesis of magnetoreception, yet also reveal that in spiny lobsters, a strong pulse altered expression of > 10% of all expressed genes, including many seemingly unrelated to sensory processes. Thus, caution is required when interpreting the effects of magnetic pulses on animal behavior.

Genomic signatures of G-protein-coupled receptor expansions reveal functional transitions in the evolution of cephalopod signal transduction.

Coleoid cephalopods show unique morphological and neural novelties, such as arms with tactile and chemosensory suckers and a large complex nervous system. The evolution of such cephalopod novelties has been attributed at a genomic level to independent gene family expansions, yet the exact association and the evolutionary timing remain unclear. In the octopus genome, one such expansion occurred in the G-protein-coupled receptors (GPCRs) repertoire, a superfamily of proteins that mediate signal transduction. Here, we assessed the evolutionary history of this expansion and its relationship with cephalopod novelties. Using phylogenetic analyses, at least two cephalopod- and two octopus-specific GPCR expansions were identified. Signatures of positive selection were analysed within the four groups, and the locations of these sequences in the Octopus bimaculoides genome were inspected. Additionally, the expression profiles of cephalopod GPCRs across various tissues were extracted from available transcriptomic data. Our results reveal the evolutionary history of cephalopod GPCRs. Unexpanded cephalopod GPCRs shared with other bilaterians were found to be mainly nervous tissue specific. By contrast, duplications that are shared between octopus and the bobtail squid or specific to the octopus' lineage generated copies with divergent expression patterns devoted to tissues outside of the brain. The acquisition of novel expression domains was accompanied by gene order rearrangement through either translocation or duplication and gene loss. Lastly, expansions showed signs of positive selection and some were found to form tandem clusters with shared conserved expression profiles in cephalopod innovations such as the axial nerve cord. Altogether, our results contribute to the understanding of the molecular and evolutionary history of signal transduction and provide insights into the role of this expansion during the emergence of cephalopod novelties and/or adaptations.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise.

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

De novo transcriptomics reveal distinct phototransduction signaling components in the retina and skin of a color-changing vertebrate, the hogfish (Lachnolaimus maximus).

Across diverse taxa, an increasing number of photoreceptive systems are being discovered in tissues outside of the eye, such as in the skin. Dermal photoreception is believed to serve a variety of functions, including rapid color change via specialized cells called chromatophores. In vitro studies of this system among color-changing fish have suggested the use of a phototransduction signaling cascade that fundamentally differs from that of the retina. Thus, the goal of this study was to identify phototransduction genes and compare their expression in the retina and skin of a color-changing fish, the hogfish Lachnolaimus maximus. De novo transcriptomics revealed the expression of genes that may underlie distinct, yet complete phototransduction cascades in L. maximus retina and skin. In contrast to the five visual opsin genes and cGMP-dependent phototransduction components expressed in the retina of L. maximus, only a single short-wavelength sensitive opsin (SWS1) and putative cAMP-dependent phototransduction components were expressed in the skin. These data suggest a separate evolutionary history of phototransduction in the retina and skin of certain vertebrates and, for the first time, indicate an expression repertoire of genes that underlie a non-retinal phototransduction pathway in the skin of a color-changing fish.

Spectral sensitivity in ray-finned fishes: diversity, ecology and shared descent.

A major goal of sensory ecology is to identify factors that underlie sensory-trait variation. One open question centers on why fishes show the greatest diversity among vertebrates in their capacity to detect color (i.e. spectral sensitivity). Over the past several decades, λmax values (photoreceptor class peak sensitivity) and chromacy (photoreceptor class number) have been cataloged for hundreds of fish species, yet the ecological basis of this diversity and the functional significance of high chromacy levels (e.g. tetra- and pentachromacy) remain unclear. In this study, we examined phylogenetic, physiological and ecological patterns of spectral sensitivity of ray-finned fishes (Actinoptergyii) via a meta-analysis of data compiled from 213 species. Across the fishes sampled, our results indicate that trichromacy is most common, ultraviolet λmax values are not found in monochromatic or dichromatic species, and increasing chromacy, including from tetra- to pentachromacy, significantly increases spectral sensitivity range. In an ecological analysis, multivariate phylogenetic latent liability modeling was performed to analyze correlations between chromacy and five hypothesized predictors (depth, habitat, diet, body coloration, body size). In a model not accounting for phylogenetic relatedness, each predictor with the exception of habitat significantly correlated with chromacy: a positive relationship in body color and negative relationships with body size, diet and depth. However, after phylogenetic correction, the only remaining correlated predictor was depth. The findings of this study indicate that phyletic heritage and depth are important factors in fish spectral sensitivity and impart caution about excluding phylogenetic comparative methods in studies of sensory trait variation.

Near absence of differential gene expression in the retina of rainbow trout after exposure to a magnetic pulse: implications for magnetoreception.

The ability to perceive the Earth's magnetic field, or magnetoreception, exists in numerous animals. Although the mechanism underlying magnetoreception has not been clearly established in any species, in salmonid fish, it is hypothesized to occur by means of crystals of magnetite associated with nervous tissue such as the brain, olfactory organ or retina. In this study, rainbow trout (Oncorhynchus mykiss) were exposed to a brief magnetic pulse known to disrupt magnetic orientation behaviour in several animals. Changes in gene expression induced by the pulse were then examined in the retina. Analyses indicated that the pulse elicited differential expression of only a single gene, gamma-crystallin M3-like (crygm3). The near absence of an effect of the magnetic pulse on gene expression in the retina stands in sharp contrast to a recent study in which 181 genes were differentially expressed in brain tissue of O. mykiss after exposure to the same pulse. Overall, our results suggest either that magnetite-based magnetoreceptors in trout are not located in the retina, or else that they are unaffected by magnetic pulses that can disrupt magnetic orientation behaviour in animals.

Genome-Wide Analysis of SNPs Is Consistent with No Domestic Dog Ancestry in the Endangered Mexican Wolf (Canis lupus baileyi).

The Mexican gray wolf (Canis lupus baileyi) was historically distributed throughout the southwestern United States and northern Mexico. Extensive predator removal campaigns during the early 20th century, however, resulted in its eventual extirpation by the mid 1980s. At this time, the Mexican wolf existed only in 3 separate captive lineages (McBride, Ghost Ranch, and Aragón) descended from 3, 2, and 2 founders, respectively. These lineages were merged in 1995 to increase the available genetic variation, and Mexican wolves were reintroduced into Arizona and New Mexico in 1998. Despite the ongoing management of the Mexican wolf population, it has been suggested that a proportion of the Mexican wolf ancestry may be recently derived from hybridization with domestic dogs. In this study, we genotyped 87 Mexican wolves, including individuals from all 3 captive lineages and cross-lineage wolves, for more than 172000 single nucleotide polymorphisms. We identified levels of genetic variation consistent with the pedigree record and effects of genetic rescue. To identify the potential to detect hybridization with domestic dogs, we compared our Mexican wolf genotypes with those from studies of domestic dogs and other gray wolves. The proportion of Mexican wolf ancestry assigned to domestic dogs was only between 0.06% (SD 0.23%) and 7.8% (SD 1.0%) for global and local ancestry estimates, respectively; and was consistent with simulated levels of incomplete lineage sorting. Overall, our results suggested that Mexican wolves lack biologically significant ancestry with dogs and have useful implications for the conservation and management of this endangered wolf subspecies.

Bringing the analysis of animal orientation data full circle: model-based approaches with maximum likelihood.

In studies of animal orientation, data are often represented as directions that can be analyzed using circular statistical methods. Although several circular statistical tests exist to detect the presence of a mean direction, likelihood-based approaches may offer advantages in hypothesis testing - especially when data are multimodal. Unfortunately, likelihood-based inference in animal orientation remains rare. Here, we discuss some of the assumptions and limitations of common circular tests and report a new R package called CircMLE to implement the maximum likelihood analysis of circular data. We illustrate the use of this package on both simulated datasets and an empirical example dataset in Chinook salmon (Oncorhynchus tshawytscha). Our software provides a convenient interface that facilitates the use of model-based approaches in animal orientation studies.

Candidate genes mediating magnetoreception in rainbow trout (Oncorhynchus mykiss).

Diverse animals use Earth's magnetic field in orientation and navigation, but little is known about the molecular mechanisms that underlie magnetoreception. Recent studies have focused on two possibilities: (i) magnetite-based receptors; and (ii) biochemical reactions involving radical pairs. We used RNA sequencing to examine gene expression in the brain of rainbow trout (Oncorhynchus mykiss) after exposure to a magnetic pulse known to disrupt magnetic orientation behaviour. We identified 181 differentially expressed genes, including increased expression of six copies of the frim gene, which encodes a subunit of the universal iron-binding and trafficking protein ferritin. Functions linked to the oxidative effects of free iron (e.g. oxidoreductase activity, transition metal ion binding, mitochondrial oxidative phosphorylation) were also affected. These results are consistent with the hypothesis that a magnetic pulse alters or damages magnetite-based receptors and/or other iron-containing structures, which are subsequently repaired or replaced through processes involving ferritin. Additionally, some genes that function in the development and repair of photoreceptive structures (e.g. crggm3, purp, prl, gcip, crabp1 and pax6) were also differentially expressed, raising the possibility that a magnetic pulse might affect structures and processes unrelated to magnetite-based magnetoreceptors.

Evolutionary and Functional Mitogenomics Associated With the Genetic Restoration of the Florida Panther.

Florida panthers are endangered pumas that currently persist in reduced patches of habitat in South Florida, USA. We performed mitogenome reference-based assemblies for most parental lines of the admixed Florida panthers that resulted from the introduction of female Texas pumas into South Florida in 1995. With the addition of 2 puma mitogenomes, we characterized 174 single nucleotide polymorphisms (SNPs) across 12 individuals. We defined 5 haplotypes (Pco1-Pco5), one of which (Pco1) had a geographic origin exclusive to Costa Rica and Panama and was possibly introduced into the Everglades National Park, Florida, prior to 1995. Haplotype Pco2 was native to Florida. Haplotypes Pco3 and Pco4 were exclusive to Texas, whereas haplotype Pco5 had an undetermined geographic origin. Phylogenetic inference suggests that haplotypes Pco1-Pco4 diverged ~202000 (95% HPDI = 83000-345000) years ago and that haplotypes Pco2-Pco4 diverged ~61000 (95% HPDI = 9000-127000) years ago. These results are congruent with a south-to-north continental expansion and with a recent North American colonization by pumas. Furthermore, pumas may have migrated from Texas to Florida no earlier than ~44000 (95% HPDI = 2000-98000) years ago. Synonymous mutations presented a greater mean substitution rate than other mitochondrial functional regions: nonsynonymous mutations, tRNAs, rRNAs, and control region. Similarly, all protein-coding genes were under predominant negative selection constraints. We directly and indirectly assessed the presence of potential deleterious SNPs in the ND2 and ND5 genes in Florida panthers prior to and as a consequence of the introduction of Texas pumas. Screenings for such variants are recommended in extant Florida panthers.

The major histocompatibility complex in Old World camelids and low polymorphism of its class II genes.

BACKGROUND: The Major Histocompatibility Complex (MHC) is a genomic region containing genes with crucial roles in immune responses. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. To counteract the high variability of pathogens, the MHC evolved into a region of considerable heterogeneity in its organization, number and extent of polymorphism. Studies of MHCs in different model species contribute to our understanding of mechanisms of immunity, diseases and their evolution. Camels are economically important domestic animals and interesting biomodels. Three species of Old World camels have been recognized: the dromedary (Camelus dromedarius), Bactrian camel (Camelus bactrianus) and the wild camel (Camelus ferus). Despite their importance, little is known about the MHC genomic region, its organization and diversity in camels. The objectives of this study were to identify, map and characterize the MHC region of Old World camelids, with special attention to genetic variation at selected class MHC II loci. RESULTS: Physical mapping located the MHC region to the chromosome 20 in Camelus dromedarius. Cytogenetic and comparative analyses of whole genome sequences showed that the order of the three major sub-regions is "Centromere - Class II - Class III - Class I". DRA, DRB, DQA and DQB exon 2 sequences encoding the antigen binding site of the corresponding class II antigen presenting molecules showed high degree of sequence similarity and extensive allele sharing across the three species. Unexpectedly low extent of polymorphism with low numbers of alleles and haplotypes was observed in all species, despite different geographic origins of the camels analyzed. The DRA locus was found to be polymorphic, with three alleles shared by all three species. DRA and DQA sequences retrieved from ancient DNA samples of Camelus dromedarius suggested that additional polymorphism might exist. CONCLUSIONS: This study provided evidence that camels possess an MHC comparable to other mammalian species in terms of its genomic localization, organization and sequence similarity. We described ancient variation at the DRA locus, monomorphic in most species. The extent of molecular diversity of MHC class II genes seems to be substantially lower in Old World camels than in other mammalian species.

The magneto-microbiome: A dataset of the metagenomic distribution of magnetotactic bacteria.

Magnetotactic bacteria (MTB) are diverse prokaryotes characterized by their ability to generate biogenic magnetic iron crystals. MTB are ubiquitous across aquatic environments, and growing evidence has indicated they may be present in association with animal microbiomes. Unfortunately, they are difficult to culture in vitro and more studies understanding their biogeographical distribution and ecological roles are needed. To provide data regarding the patterns of diversity and distribution of MTB, we screened the entire Sequence Read Archive (SRA) from the National Center for Biotechnology Information for DNA sequencing reads matching known MTB taxa. The dataset summarizes the count of reads assigned to MTB from more than 26 million SRA accessions comprising approximately 80 petabases (7.98 × 1016) of DNA. More than 396 million DNA sequencing reads were assigned to 214 MTB taxa in 691,086 (2.65 %) SRA accessions. The final dataset can be utilized by researchers to narrow their efforts in examination of both environmental and ecological roles of specific MTB or to identify potential host organisms. These data will be instrumental to further elucidating the importance and utility of these enigmatic bacteria.

Male song structure predicts offspring recruitment to the breeding population in a migratory bird.

Bird song is a classic example of a sexually selected trait, but much of the work relating individual song components to fitness has not accounted for song typically being composed of multiple, often-correlated components, necessitating a multivariate approach. We explored the role of sexual selection in shaping the complex male song of house wrens (Troglodytes aedon) by simultaneously relating its multiple components to fitness using multivariate selection analysis, which is widely used in insect and anuran studies but not in birds. The analysis revealed significant variation in the form and strength of selection acting on song across different selection episodes, from nest-site defense to recruitment of offspring to the breeding population. Males that sang more song typically employed in close communication sired more offspring that were subsequently recruited to the breeding population than those that sang more far-communication song. However, this relationship was not consistent across earlier selection episodes, as evidenced by non-linear selection acting on these song components in other contexts. Collectively, our results present a complex picture of multivariate selection on male song structure that would not be evident using univariate approaches and suggest possible trade-offs within and among song components at different points of the breeding season.

The mitochondrial genome of the pentastome parasite Raillietiella orientalis Hett, 1915 (Raillietiellida; Raillietiellidae) with notes on its phylogenetic position.

In this study we sequenced and annotated the complete mitochondrial genome of the invasive reptile parasite Raillietiella orientalis using Illumina DNA sequencing. The length of the mitogenome was 15,320 bp and had a GC content of 33.1%. The mitogenome contained 13 protein-coding genes, two ribosomal RNA genes, and 22 tRNA genes, the order of which was diagnostic of Pancrustacean mitogenomes. A phylogenetic tree constructed from the 13 protein-coding genes of R. orientalis and 26 other Pancrustacean mitogenomes supported the placement of R. orientalis as part of the monophyletic subclass Pentastomida within the Maxillopoda and sister to the subclass Branchiura.

Novel Genome Sequences of Ophidiomyces ophidiicola, the Causative Agent of Snake Fungal Disease.

Ophidiomyces ophidiicola is a globally distributed fungal pathogen of snakes. This study reports genome assemblies for three novel isolates that were derived from hosts originating in the United States, Germany, and Canada. The assemblies have a mean length of 21.4 Mbp, with coverage of 116.7×, and will contribute to wildlife disease research.

Abiotic factors shape microbial diversity in Sonoran Desert soils.

High-throughput, culture-independent surveys of bacterial and archaeal communities in soil have illuminated the importance of both edaphic and biotic influences on microbial diversity, yet few studies compare the relative importance of these factors. Here, we employ multiplexed pyrosequencing of the 16S rRNA gene to examine soil- and cactus-associated rhizosphere microbial communities of the Sonoran Desert and the artificial desert biome of the Biosphere2 research facility. The results of our replicate sampling approach show that microbial communities are shaped primarily by soil characteristics associated with geographic locations, while rhizosphere associations are secondary factors. We found little difference between rhizosphere communities of the ecologically similar saguaro (Carnegiea gigantea) and cardón (Pachycereus pringlei) cacti. Both rhizosphere and soil communities were dominated by the disproportionately abundant Crenarchaeota class Thermoprotei, which comprised 18.7% of 183,320 total pyrosequencing reads from a comparatively small number (1,337 or 3.7%) of the 36,162 total operational taxonomic units (OTUs). OTUs common to both soil and rhizosphere samples comprised the bulk of raw sequence reads, suggesting that the shared community of soil and rhizosphere microbes constitute common and abundant taxa, particularly in the bacterial phyla Proteobacteria, Actinobacteria, Planctomycetes, Firmicutes, Bacteroidetes, Chloroflexi, and Acidobacteria. The vast majority of OTUs, however, were rare and unique to either soil or rhizosphere communities and differed among locations dozens of kilometers apart. Several soil properties, particularly soil pH and carbon content, were significantly correlated with community diversity measurements. Our results highlight the importance of culture-independent approaches in surveying microbial communities of extreme environments.

Reconstruction of the Major Maternal and Paternal Lineages in the Feral New Zealand Kaimanawa Horses.

New Zealand has the fourth largest feral horse population in the world. The Kaimanawas (KHs) are feral horses descended from various domestic horse breeds released into the Kaimanawa ranges in the 19th and 20th centuries. Over time, the population size has fluctuated dramatically due to hunting, large-scale farming and forestry. Currently, the herd is managed by an annual round-up, limiting the number to 300 individuals to protect the native ecosystem. Here, we genotyped 96 KHs for uniparental markers (mitochondrial DNA, Y-chromosome) and assessed their genetic similarity with respect to other domestic horses. We show that at least six maternal and six paternal lineages contributed unequally to the KH gene pool, and today's KH population possibly represents two sub-populations. Our results indicate that three horse breeds, namely Welsh ponies, Thoroughbreds and Arabian horses had a major influence in the genetic-makeup of the extant KH population. We show that mitochondrial genetic diversity in KHs (π = 0.00687 ± 0.00355) is closer to that of the Sable Island horses (π = 0.0034 ± 0.00301), and less than other feral horse populations around the world. Our current findings, combined with ongoing genomic research, will provide insight into the population-specific genetic variation and inbreeding among KHs. This will largely advance equine research and improve the management of future breeding programs of these treasured New Zealand horse.

OptM: estimating the optimal number of migration edges on population trees using Treemix.

The software Treemix has become extensively used to estimate the number of migration events, or edges (m), on population trees from genome-wide allele frequency data. However, the appropriate number of edges to include remains unclear. Here, I show that an optimal value of m can be inferred from the second-order rate of change in likelihood (Δm) across incremental values of m. Repurposed from its original use to estimate the number of population clusters in the software Structure (ΔK), I show using simulated populations that Δm performs equally as well as current recommendations for Treemix. A demonstration of an empirical dataset from domestic dogs indicates that this method may be preferable in large, complex population histories and can prioritize migration events for subsequent investigation. The method has been implemented in a freely available R package called "OptM" and as a web application (https://rfitak.shinyapps.io/OptM/) to interface directly with the output files of Treemix.