Perry BW, Card DC, McGlothlin JW, Pasquesi GIM, Adams RH, Schield DR, Hales NR, Corbin AB, Demuth JP, Hoffmann FG, et al. Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome O'Connell M. Genome Biology and Evolution [Internet]. 2018;10 (8) :2110–2129. Publisher's VersionAbstract
Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.
Campbell-Staton SC, Bare A, Losos JB, Edwards SV, Cheviron ZA. Physiological and regulatory underpinnings of geographic variation in reptilian cold tolerance across a latitudinal cline. Molecular Ecology [Internet]. 2018;27 (9) :2243–2255. Publisher's VersionAbstract
Understanding the mechanisms that produce variation in thermal performance is a key component to investigating climatic effects on evolution and adaptation. However, disentangling the effects of local adaptation and phenotypic plasticity in shaping patterns of geographic variation in natural populations can prove challenging. Additionally, the physiological mechanisms that cause organismal dysfunction at extreme temperatures are still largely under debate. Using the green anole, Anolis carolinensis, we integrate measures of cold tolerance (CTmin ), standard metabolic rate, heart size, blood lactate concentration and RNAseq data from liver tissue to investigate geographic variation in cold tolerance and its underlying mechanisms along a latitudinal cline. We found significant effects of thermal acclimation and latitude of origin on variation in cold tolerance. Increased cold tolerance correlates with decreased rates of oxygen consumption and blood lactate concentration (a proxy for oxygen limitation), suggesting elevated performance is associated with improved oxygen economy during cold exposure. Consistent with these results, co-expression modules associated with blood lactate concentration are enriched for functions associated with blood circulation, coagulation and clotting. Expression of these modules correlates with thermal acclimation and latitude of origin. Our findings support the oxygen and capacity-limited thermal tolerance hypothesis as a potential contributor to variation in reptilian cold tolerance. Moreover, differences in gene expression suggest regulation of the blood coagulation cascade may play an important role in reptilian cold tolerance and may be the target of natural selection in populations inhabiting colder environments.
Ghanem M, Wang L, Zhang Y, Edwards S, Lu A, Ley D, El-Gazzar M. Core Genome Multilocus Sequence Typing: a Standardized Approach for Molecular Typing of Mycoplasma gallisepticum Fenwick B. Journal of Clinical Microbiology [Internet]. 2017;56 (1). Publisher's VersionAbstract
\textlessp\textgreater \textlessnamed-content content-type="genus-species"\textgreaterMycoplasma gallisepticum\textless/named-content\textgreater is the most virulent and economically important \textlessnamed-content content-type="genus-species"\textgreaterMycoplasma\textless/named-content\textgreater species for poultry worldwide. Currently, \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater strain differentiation based on sequence analysis of 5 loci remains insufficient for accurate outbreak investigation. Recently, whole-genome sequences (WGS) of many human and animal pathogens have been successfully used for microbial outbreak investigations. However, the massive sequence data and the diverse properties of different genes within bacterial genomes results in a lack of standard reproducible methods for comparisons among \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater whole genomes. Here, we proposed the development of a core genome multilocus sequence typing (cgMLST) scheme for \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater strains and field isolates. For development of this scheme, a diverse collection of 37 \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater genomes was used to identify cgMLST targets. A total of 425 \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater conserved genes (49.85% of \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater genome) were selected as core genome targets. A total of 81 \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater genomes from 5 countries on 4 continents were typed using \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater cgMLST. Analyses of phylogenetic trees generated by cgMLST displayed a high degree of agreement with geographical and temporal information. Moreover, the high discriminatory power of cgMLST allowed differentiation between \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater strains of the same outbreak. \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater cgMLST represents a standardized, accurate, highly discriminatory, and reproducible method for differentiation among \textlessnamed-content content-type="genus-species"\textgreaterM. gallisepticum\textless/named-content\textgreater isolates. cgMLST provides stable and expandable nomenclature, allowing for comparison and sharing of typing results among laboratories worldwide. cgMLST offers an opportunity to harness the tremendous power of next-generation sequencing technology in applied avian mycoplasma epidemiology at both local and global levels. \textless/p\textgreater
Chattopadhyay B, Garg KM, Gwee CY, Edwards SV, Rheindt FE. Gene flow during glacial habitat shifts facilitates character displacement in a Neotropical flycatcher radiation. BMC Evolutionary Biology [Internet]. 2017;17 (1) :210. Publisher's VersionAbstract
Pleistocene climatic fluctuations are known to be an engine of biotic diversification at higher latitudes, but their impact on highly diverse tropical areas such as the Andes remains less well-documented. Specifically, while periods of global cooling may have led to fragmentation and differentiation at colder latitudes, they may – at the same time – have led to connectivity among insular patches of montane tropical habitat with unknown consequences on diversification. In the present study we utilized \~5.5 kb of DNA sequence data from eight nuclear loci and one mitochondrial gene alongside diagnostic morphological and bioacoustic markers to test the effects of Pleistocene climatic fluctuations on diversification in a complex of Andean tyrant-flycatchers of the genus Elaenia. Population genetic and phylogenetic approaches coupled with coalescent simulations demonstrated disparate levels of gene flow between the taxon chilensis and two parapatric Elaenia taxa predominantly during the last glacial period but not thereafter, possibly on account of downward shifts of montane forest habitat linking the populations of adjacent ridges. Additionally, morphological and bioacoustic analyses revealed a distinct pattern of character displacement in coloration and vocal traits between the two sympatric taxa albiceps and pallatangae, which were characterized by a lack of gene flow. Our study demonstrates that global periods of cooling are likely to have facilitated gene flow among Andean montane Elaenia flycatchers that are more isolated from one another during warm interglacial periods such as the present era. We also identify a hitherto overlooked case of plumage and vocal character displacement, underpinning the complexities of gene flow patterns caused by Pleistocene climate change across the Andes.
Liu L, Zhang J, Rheindt FE, Lei F, Qu Y, Wang Y, Zhang Y, Sullivan C, Nie W, Wang J, et al. Genomic evidence reveals a radiation of placental mammals uninterrupted by the KPg boundary. PNAS [Internet]. 2017;114 (35) :E7282-E7290. Open Access Download
Termignoni-Garcia F, Jaramillo-Correa JP, Chable-Santos J, Liu M, Shultz AJ, Edwards SV, Escalante-Pliego P. Genomic footprints of adaptation in a cooperatively breeding tropical bird across a vegetation gradient. Molecular Ecology. 2017;26 (17) :4483-4496. PDF
Tigano A, Shultz AJ, Edwards SV, Robertson GJ, Friesen VL. Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird. Ecology and Evolution [Internet]. 2017;7 (7) :2370-2381. Open Access Download
Edwards SV, Cloutier AC, Baker AJ. Conserved Nonexonic Elements : A Novel Class of Marker for Phylogenomics. Systematic Biology [Internet]. 2017;66 (October) :1028–1044. Open Access DownloadAbstract
Noncoding markers have a particular appeal as tools for phylogenomic analysis because, at least in vertebrates, they appear less subject to strong variation in GC content among lineages. Thus far, ultraconserved elements (UCEs) and introns have been the most widely used noncoding markers. Here we analyze and study the evolutionary properties of a new type of noncoding marker, conserved nonexonic elements (CNEEs), which consists of noncoding elements that are estimated to evolve slower than the neutral rate across a set of species. Although they often include UCEs, CNEEs are distinct from UCEs because they are not ultraconserved, and, most importantly, the core region alone is analyzed, rather
than both the core and its flanking regions. Using a data set of 16 birds plus an alligator outgroup, and ∼3600–∼3800 loci per marker type, we found that although CNEEs were less variable than bioinformatically derived UCEs or introns and in some cases exhibited a slower approach to branch resolution as determined by phylogenomic subsampling, the quality of CNEE alignmentswas superior to those of the other markers, with fewer gaps and missing species. Phylogenetic resolution using coalescent approaches was comparable among the three marker types, with most nodes being fully and congruently resolved. Comparison of phylogenetic results across the three marker types indicated that one branch, the sister group to the
passerine+falcon clade,was resolved differently and with moderate (>70%) bootstrap support between CNEEs and UCEs or introns. Overall, CNEEs appear to be promising as phylogenomic markers, yielding phylogenetic resolution as high as for UCEs and introns but with fewer gaps, less ambiguity in alignments and with patterns of nucleotide substitution more consistent with the assumptions of commonly used methods of phylogenetic analysis.
Campbell-Staton SC, Cheviron ZA, Rochette N, Catchen J, Losos JB, Edwards SV. Winter storms drive rapid phenotypic, regulatory, and genomic shifts in the green anole lizard. Science. 2017;357 (6350) :495-497.Abstract
Extreme environmental perturbations offer opportunities to observe the effects of natural selection in wild populations. During the winter of 2013-2014, the southeastern United States endured an extreme cold event. We used thermal performance, transcriptomics, and genome scans to measure responses of lizard populations to storm-induced selection. We found significant increases in cold tolerance at the species' southern limit. Gene expression in southern survivors shifted toward patterns characteristic of northern populations. Comparing samples before and after the extreme winter, 14 genomic regions were differentiated in the surviving southern population; four also exhibited signatures of local adaptation across the latitudinal gradient and implicate genes involved in nervous system function. Together, our results suggest that extreme winter events can rapidly produce strong selection on natural populations at multiple biological levels that recapitulate geographic patterns of local adaptation.
Grayson P, Sin SYW, Sackton T, Edwards SV. Comparative genomics as a foundation for evo-devo studies in birds. In: Methods in Molecular Biology: Avian and Reptilian Developmental Biology. New York: Humana Press ; 2017.Abstract


Developmental genomics is a rapidly growing field and high quality genomes are a useful foundation for comparative developmental studies. A genome streamlines and simplifies the development of primers used to isolate putative regulatory regions for enhancer screens, cDNA probes for in situ hybridization, micro RNAs (miRNAs) or short hairpin RNAs(shRNA) forRNA interference (RNAi) knockdowns, mRNAs for misexpression studies, and in recent years, even guide RNAs (gRNAs) for CRISPR knockouts. A high quality genome also forms an essential reference onto which the data from numerous assays and experiments, including ChIPseq, ATAC-seq, and RNA-seq, can be mapped. Finally, much can be gleaned from comparative genomics alone, including identification of highly conserved putative regulatory regions. This chapter provides an overview of laboratory and bioinformatics protocols for DNA extraction, library preparation, library quantification and genome assembly, from fresh or frozen tissue to a draft avian genome. Generating a high quality draft genome can provide a developmental research group with excellent resources for their study organism, opening the doors to many additional assays and experiments.


Campbell-Staton SC, Edwards SV, Losos JB. Climate-mediated adaptation after mainland colonization of an ancestrally subtropical island lizard, Anolis carolinensis. Journal of Evolutionary Biology. 2016;29 (11) :2168-2180.Abstract
Climate-mediated evolution plays an integral role in species migration and range expansion. Gaining a clearer understanding of how climate affects demographic history and adaptation provides fundamental insight into the generation of intra-and interspecific diversity. In this study, we used the natural colonization of the green anole (Anolis carolinensis) from the island of Cuba to mainland North America to investigate the role of evolution at the niche, phenotypic and genetic levels after long-term establishment in a novel environment. The North American green anole occupies a broader range of thermal habitats than its Cuban sister species. We documented niche expansion in the mainland green anole, mediated primarily through adaptation to winter temperatures. Common garden experiments strongly suggest a genetic component to differences in thermal performance found between populations in different temperature regimes. Analysis of geographic variation in population structure based on 53 486 single nucleotide variants from RAD loci revealed increased genetic isolation between populations in different vs. similar thermal environments. Selection scans for environment-allele correlations reveal 19 genomic loci of known function that may have played a role in the physiological adaptation of A. carolinensis to temperate environments on the mainland.
Shultz AJ, Baker AJ, Hill GE, Nolan PM, Edwards SV. SNPs across time and space: population genomic signatures of founder events and epizootics in the House Finch (Haemorhous mexicanus). Ecology and Evolution [Internet]. 2016;6 :7475–7489. Publisher's VersionAbstract
Identifying genomic signatures of natural selection can be challenging against a background of demographic changes such as bottlenecks and population expansions. Here, we disentangle the effects of demography from selection in the House Finch (Haemorhous mexicanus) using samples collected before and after a pathogen-induced selection event. Using ddRADseq, we genotyped over 18,000 SNPs across the genome in native pre-epizootic western US birds, introduced birds from Hawaii and the eastern United States, post-epizootic eastern birds, and western birds sampled across a similar time span. We found 14% and 7% reductions in nucleotide diversity, respectively, in Hawaiian and pre-epizootic eastern birds relative to pre-epizootic western birds, as well as elevated levels of linkage disequilibrium and other signatures of founder events. Despite finding numerous significant frequency shifts (outlier loci) between pre-epizootic native and introduced populations, we found no signal of reduced genetic diversity, elevated linkage disequilibrium, or outlier loci as a result of the epizootic. Simulations demonstrate that the proportion of outliers associated with founder events could be explained by genetic drift. This rare view of genetic evolution across time in an invasive species provides direct evidence that demographic shifts like founder events have genetic consequences more widespread across the genome than natural selection.
Edwards SV. Inferring species trees. In: Kliman R Encyclopedia of Evolutionary Biology. Vol. 4. New York: Elsevier Inc. ; 2016. pp. 236-244.
Edwards SV, Xi Z, Janke A, Faircloth BC, McCormack JE, Glenn TC, Zhong B, Wu S, Lemmon EM, Lemmon AR. Implementing and testing the multispecies coalescent model: A valuable paradigm for phylogenomics. Molecular Phylogenetics and Evolution. 2016;94 :447-462.Abstract
In recent articles published in
Molecular Phylogenetics and Evolution
, Mark Springer and John Gatesy (S&G)
present numerous criticisms of recent implementations and testing of the multispecies coalescent (MSC)
model in phylogenomics, popularly known as ‘‘species tree” methods. After pointing out errors in align-
ments and gene tree rooting in recent phylogenomic data sets, particularly in Song et al. (2012) on mam-
mals and Xi et al. (2014) on plants, they suggest that these errors seriously compromise the conclusions
of these studies. Additionally, S&G enumerate numerous perceived violated assumptions and deficiencies
in the application of the MSC model in phylogenomics, such as its assumption of neutrality and in par-
ticular the use of transcriptomes, which are deemed inappropriate for the MSC because the constituent
exons often subtend large regions of chromosomes within which recombination is substantial. We
acknowledge these previously reported errors in recent phylogenomic data sets, but disapprove of
S&G’s excessively combative and taunting tone. We show that these errors, as well as two nucleotide
sorting methods used in the analysis of
, have little impact on the conclusions of those papers.
Moreover, several concepts introduced by S&G and an appeal to ‘‘first principles” of phylogenetics in an
attempt to discredit MSC models are invalid and reveal numerous misunderstandings of the MSC.
Contrary to the claims of S&G we show that recent computer simulations used to test the robustness
of MSC models are not circular and do not unfairly favor MSC models over concatenation. In fact, although
both concatenation and MSC models clearly perform well in regions of tree space with long branches and
little incomplete lineage sorting (ILS), simulations reveal the erratic behavior of concatenation when sub-
jected to data subsampling and its tendency to produce spuriously confident yet conflicting results in
regions of parameter space where MSC models still perform well. S&G’s claims that MSC models explain
little or none (0–15%) of the observed gene tree heterogeneity observed in a mammal data set and that
MSC models assume ILS as the only source of gene tree variation are flawed. Overall many of their crit-
icisms of MSC models are invalidated when concatenation is appropriately viewed as a special case of the
MSC, which in turn is a special case of emerging network models in phylogenomics. We reiterate that
there is enormous promise and value in recent implementations and tests of the MSC and look forward
to its increased use and refinement in phylogenomics.
Edwards SV. Phylogenomic subsampling: a brief review. Zoologica Scripta. 2016;45 (s1) :63-74. PDF
Edwards SV, Potter S, Schmitt CJ, Bragg JG, Moritz C. Reticulation, divergence, and the phylogeography–phylogenetics continuum. Proceedings of the National Academy of Sciences. 2016;113 :8025-8032.Abstract
Phylogeography, and its extensions into comparative phylogeogra-
phy, have their roots in the layering of gene trees across geography,
a paradigm that was greatly facilitated by the nonrecombining, fast
evolution provided by animal mtDNA. As phylogeography moves
into the era of next-generation sequencing, the specter of reticula-
tion at several levels
within loci and genomes in the form of re-
combination and across populations and species in the form of
has raised its head with a prominence even greater
than glimpsed during the nuclear gene PCR era. Here we explore the
theme of reticulation in comparative phylogeography, speciation
analysis, and phylogenomics, and ask how the centrality of gene
trees has fared in the next-generation era. To frame these issues,
we first provide a snapshot of multilocus phylogeographic studies
across the Carpentarian Barrier, a prominent biogeographic barrier
dividing faunas spanning the monsoon tropics in northern Australia.
We find that divergence across this barrier is evident in most spe-
cies, but is heterogeneous in time and demographic history, often
reflecting the taxonomic distinctness of lineages spanning it. We
then discuss a variety of forces generating reticulate patterns in
phylogeography, including introgression, contact zones, and the
potential selection-driven outliers on next-generation molecular
markers. We emphasize the continued need for demographic models
incorporating reticulation at the level of genomes and populations,
and conclude that gene trees, whether explicit or implicit, should
continue to play a role in the future of phylogeography.
Miller HC, O'Meally D, Ezaz T, Amemiya C, Marshall-Graves JA, Edwards S. Major Histocompatibility Complex Genes Map to Two Chromosomes in an Evolutionarily Ancient Reptile, the Tuatara Sphenodon punctatus. G3-Genes, Genomes, Genetics. 2015;5 (7) :1439-1451.Abstract
Major histocompatibility complex (MHC) genes are a central component of the vertebrate immune system and usually exist in a single genomic region. However, considerable differences in MHC organization and size exist between different vertebrate lineages. Reptiles occupy a key evolutionary position for understanding how variation in MHC structure evolved in vertebrates, but information on the structure of the MHC region in reptiles is limited. In this study, we investigate the organization and cytogenetic location of MHC genes in the tuatara (Sphenodon punctatus), the sole extant representative of the early-diverging reptilian order Rhynchocephalia. Sequencing and mapping of 12 clones containing class I and II MHC genes from a bacterial artificial chromosome library indicated that the core MHC region is located on chromosome 13q. However, duplication and translocation of MHC genes outside of the core region was evident, because additional class IMHC genes were located on chromosome 4p. We found a total of seven class I sequences and 11 class II beta sequences, with evidence for duplication and pseudogenization of genes within the tuatara lineage. The tuatara MHC is characterized by high repeat content and low gene density compared with other species and we found no antigen processing or MHC framework genes on the MHC gene-containing clones. Our findings indicate substantial differences in MHC organization in tuatara compared with mammalian and avian MHCs and highlight the dynamic nature of the MHC. Further sequencing and annotation of tuatara and other reptile MHCs will determine if the tuatara MHC is representative of nonavian reptiles in general.
Liu L, Edwards SV. Comment on “Statistical binning enables an accurate coalescent-based estimation of the avian tree”. Science. 2015;350 :171.Abstract

Mirarab et al. (Research Article, 12 December 2014, p. 1250463) introduced statistical binning to improve the signal in phylogenetic methods using the multispecies coalescent model. We show that all forms of binning—naïve, statistical, and weighted statistical—display poor performance and are statistically inconsistent in large regions of parameter space, unlike unbinned sequence data used with species tree methods.

Liu L, Xi Z, Wu S, Davis CC, Edwards SV. Estimating phylogenetic trees from genome-scale data. Annals of the New York Academy of Sciences. 2015;1360 :36–53.Abstract

The heterogeneity of signals in the genomes of diverse organisms poses challenges for traditional phylogenetic analysis. Phylogenetic methods known as “species tree” methods have been proposed to directly address one important source of gene tree heterogeneity, namely the incomplete lineage sorting that occurs when evolving lineages radiate rapidly, resulting in a diversity of gene trees from a single underlying species tree. Here we review theory and empirical examples that help clarify conflicts between species tree and concatenation methods, and misconceptions in the literature about the performance of species tree methods. Considering concatenation as a special case of the multispecies coalescent model helps explain differences in the behavior of the two methods on phylogenomic data sets. Recent work suggests that species tree methods are more robust than concatenation approaches to some of the classic challenges of phylogenetic analysis, including rapidly evolving sites in DNA sequences and long-branch attraction. We show that approaches, such as binning, designed to augment the signal in species tree analyses can distort the distribution of gene trees and are inconsistent. Computationally efficient species tree methods incorporating biological realism are a key to phylogenetic analysis of whole-genome data.

Dierickx EG, Shultz AJ, Sato F, Hiraoka T, Edwards SV. Morphological and genomic comparisons of Hawaiian and Japanese Black-footed Albatrosses (Phoebastria nigripes) using double digest RADseq: implications for conservation. Evolutionary Applications. 2015;8 :662-678.Abstract

Evaluating the genetic and demographic independence of populations of threatened species is important for determining appropriate conservation measures, but different technologies can yield different conclusions. Despite multiple studies, the taxonomic status and extent of gene flow between the main breeding populations of Black-footed Albatross (Phoebastria nigripes), a Near-Threatened philopatric seabird, are still controversial. Here, we employ double digest RADseq to quantify the extent of genomewide divergence and gene flow in this species. Our genomewide data set of 9760 loci containing 3455 single nucleotide polymorphisms yielded estimates of genetic diversity and gene flow that were generally robust across seven different filtering and sampling protocols and suggest a low level of genomic variation ( per site=similar to 0.00002-0.00028), with estimates of effective population size (N-e=similar to 500-15881) falling far below current census size. Genetic differentiation was small but detectable between Japan and Hawaii (F-ST approximate to 0.038 0.049), with no F-ST outliers. Additionally, using museum specimens, we found that effect sizes of morphological differences by sex or population rarely exceeded 4%. These patterns suggest that the Hawaiian and Japanese populations exhibit small but significant differences and should be considered separate management units, although the evolutionary and adaptive consequences of this differentiation remain to be identified.