Historical museum samples enable the examination of divergent and parallel evolution during invasion.

During the Anthropocene, Earth has experienced unprecedented habitat loss, native species decline and global climate change. Concurrently, greater globalization is facilitating species movement, increasing the likelihood of alien species establishment and propagation. There is a great need to understand what influences a species' ability to persist or perish within a new or changing environment. Examining genes that may be associated with a species' invasion success or persistence informs invasive species management, assists with native species preservation and sheds light on important evolutionary mechanisms that occur in novel environments. This approach can be aided by coupling spatial and temporal investigations of evolutionary processes. Here we use the common starling, Sturnus vulgaris, to identify parallel and divergent evolutionary change between contemporary native and invasive range samples and their common ancestral population. To do this, we use reduced-representation sequencing of native samples collected recently in northwestern Europe and invasive samples from Australia, together with museum specimens sampled in the UK during the mid-19th century. We found evidence of parallel selection on both continents, possibly resulting from common global selective forces such as exposure to pollutants. We also identified divergent selection in these populations, which might be related to adaptive changes in response to the novel environment encountered in the introduced Australian range. Interestingly, signatures of selection are equally as common within both invasive and native range contemporary samples. Our results demonstrate the value of including historical samples in genetic studies of invasion and highlight the ongoing and occasionally parallel role of adaptation in both native and invasive ranges.

Emergence of an evolutionary innovation: Gene expression differences associated with the transition between oviparity and viviparity.

Our understanding of the evolution of complex biological traits is greatly advanced by examining taxa with intermediate phenotypes. The transition from oviparity (egg-laying) to viviparity (live-bearing) has occurred independently in many animal lineages, but there are few phenotypic intermediates. The lizard Saiphos equalis exhibits bimodal reproduction, with some viviparous populations, and other oviparous populations with long egg-retention, a rare trait where most of embryonic development occurs inside the mother prior to late ovipositioning. We posit that oviparous S. equalis represent an intermediate form between "true" oviparity and viviparity. We used transcriptomics to compare uterine gene expression in these two phenotypes, and provide a molecular model for the genetic control and evolution of reproductive mode. Many genes are differentially expressed throughout the reproductive cycle of both phenotypes, which have clearly different gene expression profiles overall. The differentially expressed genes within oviparous and viviparous individuals have broadly similar biological functions putatively important for sustaining embryos, including uterine remodelling, respiratory gas and water exchange, and immune regulation. These functional similarities indicate either that long egg-retention is an exaptation for viviparity, or might reflect parallel evolution of similar gravidity-related changes in gene expression in long egg-retention oviparity. In contrast, gene expression changes across the reproductive cycle of long egg-retaining oviparous S. equalis are dramatically different from those of "true" oviparous skinks (such as Lampropholis guichenoti), supporting our assertion that oviparous S. equalis exhibit an intermediate phenotype between "true" oviparity and viviparity.

The Regulation of Uterine Proinflammatory Gene Expression during Pregnancy in the Live-Bearing Lizard, Pseudemoia entrecasteauxii.

The evolutionary transition from egg-laying to live-bearing in amniote vertebrates (reptiles and mammals) requires the development of a closer association between the maternal and embryonic tissue to facilitate gas and nutrient exchange with the embryo. Because the embryo is an allograft to the father and mother, it could be considered foreign by the maternal immune system and thus be immunologically rejected during pregnancy. In eutherian ("placental") mammals, the proinflammatory genes interleukin 1B (IL1B), tumor necrosis factor (TNF) and tumor necrosis factor receptor superfamily 1A (TNFRSF1A) are tightly regulated in the pregnant uterus to prevent embryonic rejection. We tested whether inflammation is similarly regulated in pregnant viviparous reptiles by comparing the expression of IL1B, TNF, and TNFRSF1A in the pregnant and nonpregnant uterus of the viviparous lizard, Pseudemoia entrecasteauxii. We found statistically significant support for the downregulation of pregnant uterine TNF mRNA expression in P. entrecasteauxii, but no statistically significant changes in mRNA expression of TNFRSF1A or IL1B between pregnant and nonpregnant uteri. Although these genes are apparently not regulated at the transcriptional level, our immunofluorescence microscopy analyses nonetheless demonstrate that the IL1B proteins are stored intracellularly during pregnancy, possibly resulting in inhibition of inflammatory response. We therefore conclude that processes of both transcriptional (TNF) and posttranslational (IL1B) gene regulation may reduce inflammation in the pregnant uterus of this viviparous reptile. Our study is important because it demonstrates that regulating the maternal immune system to prevent embryonic rejection may be important in reptilian pregnancy as it is in mammalian pregnancy.

Comparative genomics of hormonal signaling in the chorioallantoic membrane of oviparous and viviparous amniotes.

In oviparous amniotes (reptiles, birds, and mammals) the chorioallantoic membrane (CAM) lines the inside of the egg and acts as the living point of contact between the embryo and the outside world. In livebearing (viviparous) amniotes, communication during embryonic development occurs across placental tissues, which form between the uterine tissue of the mother and the CAM of the embryo. In both oviparous and viviparous taxa, the CAM is at the interface of the embryo and the external environment and can transfer signals from there to the embryo proper. To understand the evolution of placental hormone production in amniotes, we examined the expression of genes involved in hormone synthesis, metabolism, and hormone receptivity in the CAM of species across the amniote phylogeny. We collected transcriptome data for the chorioallantoic membranes of the chicken (oviparous), the lizards Lerista bougainvillii (both oviparous and viviparous populations) and Pseudemoia entrecasteauxii (viviparous), and the horse Equus caballus (viviparous). The viviparous taxa differ in their mechanisms of nutrient provisioning: L. bougainvillii is lecithotrophic (embryonic nourishment is provided via the yolk only), but P. entrecasteauxii and the horse are placentotrophic (embryos are nourished via placental transport). Of the 423 hormone-related genes that we examined, 91 genes are expressed in all studied species, suggesting that the chorioallantoic membrane ancestrally had an endocrine function. Therefore, the chorioallantoic membrane appears to be a highly hormonally active organ in all amniotes. No genes are expressed only in viviparous species, suggesting that the evolution of viviparity has not required the recruitment of any specific hormone-related genes. Our data suggest that the endocrine function of the CAM as a placental tissue evolved in part through co-option of ancestral gene expression patterns.

Allelic expression of mammalian imprinted genes in a matrotrophic lizard, Pseudemoia entrecasteauxii.

Genomic imprinting is a process that results in the differential expression of genes depending on their parent of origin. It occurs in both plants and live-bearing mammals, with imprinted genes typically regulating the ability of an embryo to manipulate the maternal provision of nutrients. Genomic imprinting increases the potential for selection to act separately on paternally and maternally expressed genes, which increases the number of opportunities that selection can facilitate embryonic control over maternal nutrient provision. By looking for imprinting in an independent matrotrophic lineage, the viviparous lizard Pseudemoia entrecasteauxii (Scincidae), we test the hypothesis that genomic imprinting facilitates the evolution of substantial placental nutrient transport to embryos (matrotrophy). We sequenced transcriptomes from the embryonic component of lizard placentae to determine whether there are parent-of-origin differences in expression of genes that are imprinted in mammals. Of these genes, 19 had sufficiently high expression in the lizard to identify polymorphisms in transcribed sequences. We identified bi-allelic expression in 17 genes (including insulin-like growth factor 2), indicating that neither allele was imprinted. These data suggest that either genomic imprinting has not evolved in this matrotrophic skink or, if it has, it has evolved in different genes to mammals. We outline how these hypotheses can be tested. This study highlights important differences between mammalian and reptile pregnancy and the absence of any shared imprinting genes reflects fundamental differences in the way that pregnancy has evolved in these two lineages.

Evaluating the performance of anchored hybrid enrichment at the tips of the tree of life: a phylogenetic analysis of Australian Eugongylus group scincid lizards.

BACKGROUND: High-throughput sequencing using targeted enrichment and transcriptomic methods enables rapid construction of phylogenomic data sets incorporating hundreds to thousands of loci. These advances have enabled access to an unprecedented amount of nucleotide sequence data, but they also pose new questions. Given that the loci targeted for enrichment are often highly conserved, how informative are they at different taxonomic scales, especially at the intraspecific/phylogeographic scale? We investigate this question using Australian scincid lizards in the Eugongylus group (Squamata: Scincidae). We sequenced 415 anchored hybrid enriched (AHE) loci for 43 individuals and mined 1650 exons (1648 loci) from transcriptomes (transcriptome mining) from 11 individuals, including multiple phylogeographic lineages within several species of Carlia, Lampropholis, and Saproscincus skinks. We assessed the phylogenetic information content of these loci at the intergeneric, interspecific, and phylogeographic scales. As a further test of the utility at the phylogeographic scale, we used the anchor hybrid enriched loci to infer lineage divergence parameters using coalescent models of isolation with migration. RESULTS: Phylogenetic analyses of both data sets inferred very strongly supported trees at all taxonomic levels. Further, AHE loci yielded estimates of divergence times between closely related lineages that were broadly consistent with previous population-level analyses. CONCLUSIONS: Anchored-enriched loci are useful at the deep phylogeny and phylogeographic scales. Although overall phylogenetic support was high throughout the Australian Eugongylus group phylogeny, there were nonetheless some conflicting or unresolved relationships, especially regarding the placement of Pseudemoia, Cryptoblepharus, and the relationships amongst closely-related species of Tasmanian Niveoscincus skinks.

Which came first: The lizard or the egg? Robustness in phylogenetic reconstruction of ancestral states.

Changes in parity mode between egg-laying (oviparity) and live-bearing (viviparity) have occurred repeatedly throughout vertebrate evolution. Oviparity is the ancestral amniote state, and viviparity has evolved many times independently within amniotes (especially in lizards and snakes), with possibly a few reversions to oviparity. In amniotes, the shelled egg is considered a complex structure that is unlikely to re-evolve if lost (i.e., it is an example of Dollo's Principle). However, a recent ancestral state reconstruction analysis concluded that viviparity was the ancestral state of squamate reptiles (lizards and snakes), and that oviparity re-evolved from viviparity many times throughout the evolutionary history of squamates. Here, we re-evaluate support for this provocative conclusion by testing the sensitivity of the analysis to model assumptions and estimates of squamate phylogeny. We found that the models and methods used for parity mode reconstruction are highly sensitive to the specific estimate of phylogeny used, and that the point estimate of phylogeny used to suggest that viviparity is the root state of the squamate tree is far from an optimal phylogenetic solution. The ancestral state reconstructions are also highly sensitive to model choice and specific values of model parameters. A method that is designed to account for biases in taxon sampling actually accentuates, rather than lessens, those biases with respect to ancestral state reconstructions. In contrast to recent conclusions from the same data set, we find that ancestral state reconstruction analyses provide highly equivocal support for the number and direction of transitions between oviparity and viviparity in squamates. Moreover, the reconstructions of ancestral parity state are highly dependent on the assumptions of each model. We conclude that the common ancestor of squamates was oviparous, and subsequent evolutionary transitions to viviparity were common, but reversals to oviparity were rare. The three putative reversals to oviparity with the strongest phylogenetic support occurred in the snakes Eryx jayakari and Lachesis, and the lizard, Liolaemus calchaqui. Our results emphasize that because the conclusions of ancestral state reconstruction studies are often highly sensitive to the methods and assumptions of analysis, researchers should carefully consider this sensitivity when evaluating alternative hypotheses of character-state evolution.

Ancestral state reconstructions require biological evidence to test evolutionary hypotheses: A case study examining the evolution of reproductive mode in squamate reptiles.

To understand evolutionary transformations it is necessary to identify the character states of extinct ancestors. Ancestral character state reconstruction is inherently difficult because it requires an accurate phylogeny, character state data, and a statistical model of transition rates and is fundamentally constrained by missing data such as extinct taxa. We argue that model based ancestral character state reconstruction should be used to generate hypotheses but should not be considered an analytical endpoint. Using the evolution of viviparity and reversals to oviparity in squamates as a case study, we show how anatomical, physiological, and ecological data can be used to evaluate hypotheses about evolutionary transitions. The evolution of squamate viviparity requires changes to the timing of reproductive events and the successive loss of features responsible for building an eggshell. A reversal to oviparity requires that those lost traits re-evolve. We argue that the re-evolution of oviparity is inherently more difficult than the reverse. We outline how the inviability of intermediate phenotypes might present physiological barriers to reversals from viviparity to oviparity. Finally, we show that ecological data supports an oviparous ancestral state for squamates and multiple transitions to viviparity. In summary, we conclude that the first squamates were oviparous, that frequent transitions to viviparity have occurred, and that reversals to oviparity in viviparous lineages either have not occurred or are exceedingly rare. As this evidence supports conclusions that differ from previous ancestral state reconstructions, our paper highlights the importance of incorporating biological evidence to evaluate model-generated hypotheses.

The evolution of viviparity: molecular and genomic data from squamate reptiles advance understanding of live birth in amniotes.

Squamate reptiles (lizards and snakes) are an ideal model system for testing hypotheses regarding the evolution of viviparity (live birth) in amniote vertebrates. Viviparity has evolved over 100 times in squamates, resulting in major changes in reproductive physiology. At a minimum, all viviparous squamates exhibit placentae formed by the appositions of maternal and embryonic tissues, which are homologous in origin with the tissues that form the placenta in therian mammals. These placentae facilitate adhesion of the conceptus to the uterus as well as exchange of oxygen, carbon dioxide, water, sodium, and calcium. However, most viviparous squamates continue to rely on yolk for nearly all of their organic nutrition. In contrast, some species, which rely on the placenta for at least a portion of organic nutrition, exhibit complex placental specializations associated with the transport of amino acids and fatty acids. Some viviparous squamates also exhibit reduced immunocompetence during pregnancy, which could be the result of immunosuppression to protect developing embryos. Recent molecular studies using both candidate-gene and next-generation sequencing approaches have suggested that at least some of the genes and gene families underlying these phenomena play similar roles in the uterus and placenta of viviparous mammals and squamates. Therefore, studies of the evolution of viviparity in squamates should inform hypotheses of the evolution of viviparity in all amniotes, including mammals.

Relaxed clocks and inferences of heterogeneous patterns of nucleotide substitution and divergence time estimates across whales and dolphins (Mammalia: Cetacea).

Various nucleotide substitution models have been developed to accommodate among lineage rate heterogeneity, thereby relaxing the assumptions of the strict molecular clock. Recently developed "uncorrelated relaxed clock" and "random local clock" (RLC) models allow decoupling of nucleotide substitution rates between descendant lineages and are thus predicted to perform better in the presence of lineage-specific rate heterogeneity. However, it is uncertain how these models perform in the presence of punctuated shifts in substitution rate, especially between closely related clades. Using cetaceans (whales and dolphins) as a case study, we test the performance of these two substitution models in estimating both molecular rates and divergence times in the presence of substantial lineage-specific rate heterogeneity. Our RLC analyses of whole mitochondrial genome alignments find evidence for up to ten clade-specific nucleotide substitution rate shifts in cetaceans. We provide evidence that in the uncorrelated relaxed clock framework, a punctuated shift in the rate of molecular evolution within a subclade results in posterior rate estimates that are either misled or intermediate between the disparate rate classes present in baleen and toothed whales. Using simulations, we demonstrate abrupt changes in rate isolated to one or a few lineages in the phylogeny can mislead rate and age estimation, even when the node of interest is calibrated. We further demonstrate how increasing prior age uncertainty can bias rate and age estimates, even while the 95% highest posterior density around age estimates decreases; in other words, increased precision for an inaccurate estimate. We interpret the use of external calibrations in divergence time studies in light of these results, suggesting that rate shifts at deep time scales may mislead inferences of absolute molecular rates and ages.

Developmental sequences of squamate reptiles are taxon specific.

Recent studies in comparative vertebrate embryology have focused on two related questions. One concerns the existence of a phylotypic period, or indeed any period, during development in which sequence variation among taxa is constrained. The second question concerns the degree to which developmental characters exhibit a phylogenetic signal. These questions are important because they underpin attempts to understand the evolution of developmental characters and their links to adult morphology. To address these questions, we compared the sequence of developmental events spanning the so-called phylotypic period of vertebrate development in squamate reptiles (lizards and snakes), from the formation of the primary optic placode to the first appearance of scale anlagen. We used Bayesian phylogenetic ancestral state reconstruction analyses and estimates of Bayesian posterior probabilities of the rank order of developmental events to determine the level of support for phylogenetically associated variation in development. We assessed the amount of variation in event sequences by plotting the proportions of reconstructed ranks (excluding unlikely events, PP < 0.05) associated with each event. Sequence variability was the lowest towards the middle of the phylotypic period and involved three events (allantois contacts chorion, maximum number of pharyngeal slits, and appearance of the apical epidermal ridge [AER]); these events each had only two reconstructed ranks. Squamate clades also differed in the rank order of developmental events. Of the 20 events in our analyses, 12 had strongly supported (PP ≥ 0.95) sequence ranks that differed at two or more internal nodes of the tree. For example, gekkotans are distinguished by the late appearance of the allantois bud compared to all other squamates (ranks 7 and 8 vs. rank 3, respectively) and Serpentes are distinguished by the earlier completion of torsion (rank 3) compared to acrodonts and pleurodonts (ranks 7 and 5, respectively). Clade specific sequences of developmental events mean that investigators should not extend observations on the development on particular squamate species to distantly related taxa for use in comparative studies.

Evidence for determinism in species diversification and contingency in phenotypic evolution during adaptive radiation.

Adaptive radiation (AR) theory predicts that groups sharing the same source of ecological opportunity (EO) will experience deterministic species diversification and morphological evolution. Thus, deterministic ecological and morphological evolution should be correlated with deterministic patterns in the tempo and mode of speciation for groups in similar habitats and time periods. We test this hypothesis using well-sampled phylogenies of four squamate groups that colonized the New World (NW) in the Late Oligocene. We use both standard and coalescent models to assess species diversification, as well as likelihood models to examine morphological evolution. All squamate groups show similar early pulses of speciation, as well as diversity-dependent ecological limits on clade size at a continental scale. In contrast, processes of morphological evolution are not easily predictable and do not show similar pulses of early and rapid change. Patterns of morphological and species diversification thus appear uncoupled across these groups. This indicates that the processes that drive diversification and disparification are not mechanistically linked, even among similar groups of taxa experiencing the same sources of EO. It also suggests that processes of phenotypic diversification cannot be predicted solely from the existence of an AR or knowledge of the process of diversification.

Molecular phylogenetics of squirrelfishes and soldierfishes (Teleostei: Beryciformes: Holocentridae): reconciling more than 100 years of taxonomic confusion.

Squirrelfishes and soldierfishes (Holocentridae) are among the most conspicuous species in the nocturnal reef fish community. However, there is no clear consensus regarding their evolutionary relationships, which is reflected in a complicated taxonomic history. We collected DNA sequence data from multiple single copy nuclear genes and one mitochondrial gene sampled from over fifty percent of the recognized holocentrid species and infer the first species-level phylogeny of the Holocentridae. Our results strongly support the monophyly of the clades Myripristinae (soldierfishes) and Holocentrinae (squirrelfishes). The molecular phylogenies differ with regard to previous hypotheses of relationships within the Myriprisitinae, resolving a clade of cryptic reef associated and deep water non-reef dwelling lineages (Corniger+Plectrypops+Ostichthys) that is the sister lineage to a monophyletic Myripristis. Within Holocentrinae, Neoniphon and Sargocentron are strongly supported as paraphyletic, while Holocentrus is nested within Sargocentron. Using Bayesian ancestral state reconstruction methods, we demonstrate the taxonomically diagnostic characters for Neoniphon and Sargocentron likely represent character states with a complex evolutionary history that is not reflective of shared common ancestry. We propose a new classification for Holocentrinae, recognizing four lineages that are treated as genera: Sargocentron Fowler, 1904, Holocentrus Scopoli, 1777, Flameo Jordan and Evermann, 1898, and Neoniphon Castelnau, 1875.

Morphology and development of the placentae in Eulamprus quoyii group skinks (Squamata: Scincidae).

Frequent evolutionary changes in reproductive mode have produced a wide range of placental structures in viviparous squamate reptiles. Closely related species with different placental structures and resolved phylogenetic relationships are particularly useful for reconstructing how placentae might have transformed during the evolutionary process. We used light microscopy to study placental morphology in mid- to late stage embryos of four closely related species of Eulamprus, a genus of viviparous scincid lizards that we had reason to suspect may display significant interspecific variation in placental morphology. Embryos from all four species possess a chorioallantoic placenta, an omphaloplacenta and an interomphalopleuric membrane, characteristics present in other viviparous skinks. However, unlike other viviparous skinks but characteristic of oviparous skinks, the allantois expands to surround the yolk sac in each species, supplanting the omphalopleure with a larger area of chorioallantois until a chorioallantoic placenta surrounds the entire egg in one specimen that is only a few days from birth. All four Eulamprus species share the same extraembryonic membrane morphology, but the cellular morphology of the uterine epithelium in the chorioallantoic placenta and omphaloplacenta varies between species. We determined that the interomphalopleuric membrane is a shared derived character of the Eulamprus quoyii species group. New phylogenetic information indicates that variation in the chorioallantoic placenta is a result of two independent transitions, but that variation in the omphaloplacenta can be explained using a single change within the species studied. Our results indicate that E. quoyii group skinks are a valuable model for investigating the evolution of viviparity, as extraembryonic membrane development in these species shows features characteristic of both oviparous and viviparous skinks.

Accommodating heterogenous rates of evolution in molecular divergence dating methods: an example using intercontinental dispersal of Plestiodon (Eumeces) lizards.

Identifying and dating historical biological events is a fundamental goal of evolutionary biology, and recent analytical advances permit the modeling of factors known to affect both the accuracy and the precision of molecular date estimates. As the use of multilocus data sets becomes increasingly routine, it becomes more important to evaluate the potentially confounding effects of rate heterogeneity both within (e.g., codon positions) and among loci when estimating divergence times. Here, using Plestiodon lizards as a test case, we examine the effects of accommodating rate heterogeneity among data partitions on divergence time estimation. Plestiodon inhabits both East Asia and North America, yet both the geographic origin of the genus and timing of dispersal between the continents have been debated. For each of the eight independently evolving loci and a combined data set, we conduct single model and partitioned analyses. We found that extreme saturation has obscured the underlying rate of evolution in the mitochondrial DNA (mtDNA), resulting in severe underestimation of the rate in this locus. As a result, the age of the crown Plestiodon clade was overestimated by 15-17 Myr by the unpartitioned analysis of the combined loci data. However, the application of partition-specific models to the combined data resulted in ages that were fully congruent with those inferred by the individual nuclear loci. Although partitioning improved divergence date estimates of the mtDNA-only analysis, the ages were nonetheless overestimated, thus indicating an inadequacy of our current models to capture the complex nature of mtDNA evolution in over large time scales. Finally, the statistically incongruent age distributions inferred by the partitioned and unpartitioned analyses of the combined data support mutually exclusive hypotheses of the timing of intercontinental dispersal of Plestiodon from Asia to North America. Analyses that best capture the rate of evolution in the combined data set infer that this exchange occurred via Beringia ∼18.0-30 Ma.

Transcript- and annotation-guided genome assembly of the European starling.

The European starling, Sturnus vulgaris, is an ecologically significant, globally invasive avian species that is also suffering from a major decline in its native range. Here, we present the genome assembly and long-read transcriptome of an Australian-sourced European starling (S. vulgaris vAU), and a second, North American, short-read genome assembly (S. vulgaris vNA), as complementary reference genomes for population genetic and evolutionary characterization. S. vulgaris vAU combined 10× genomics linked-reads, low-coverage Nanopore sequencing, and PacBio Iso-Seq full-length transcript scaffolding to generate a 1050 Mb assembly on 6222 scaffolds (7.6 Mb scaffold N50, 94.6% busco completeness). Further scaffolding against the high-quality zebra finch (Taeniopygia guttata) genome assigned 98.6% of the assembly to 32 putative nuclear chromosome scaffolds. Species-specific transcript mapping and gene annotation revealed good gene-level assembly and high functional completeness. Using S. vulgaris vAU, we demonstrate how the multifunctional use of PacBio Iso-Seq transcript data and complementary homology-based annotation of sequential assembly steps (assessed using a new tool, saaga) can be used to assess, inform, and validate assembly workflow decisions. We also highlight some counterintuitive behaviour in traditional busco metrics, and present buscomp, a complementary tool for assembly comparison designed to be robust to differences in assembly size and base-calling quality. This work expands our knowledge of avian genomes and the available toolkit for assessing and improving genome quality. The new genomic resources presented will facilitate further global genomic and transcriptomic analysis on this ecologically important species.

Reevaluating claims of ecological speciation in Halichoeres bivittatus.

Allopatry has traditionally been viewed as the primary driver of speciation in marine taxa, but the geography of the marine environment and the larval dispersal capabilities of many marine organisms render this view somewhat questionable. In marine fishes, one of the earliest and most highly cited empirical examples of ecological speciation with gene flow is the slippery dick wrasse, Halichoeres bivittatus. Evidence for this cryptic or incipient speciation event was primarily in the form of a deep divergence in a single mitochondrial locus between the northern and southern Gulf of Mexico, combined with a finding that these two haplotypes were associated with different habitat types ("tropical" vs. "subtropical") in the Florida Keys and Bermuda, where they overlap. Here, we examine habitat assortment in the Florida Keys using a broader sampling of populations and habitat types than were available for the original study. We find no evidence to support the claim that haplotype frequencies differ between habitat types, and little evidence to support any differences between populations in the Keys. These results undermine claims of ecological speciation with gene flow in Halichoeres bivittatus. Future claims of this type should be supported by multiple lines of evidence that illuminate potential mechanisms and allow researchers to rule out alternative explanations for spatial patterns of genetic differences.

Comparing the dietary niche overlap and ecomorphological differences between invasive Hemidactylus mabouia geckos and a native gecko competitor.

Hemidactylus mabouia is one of the most successful, widespread invasive reptile species and has become ubiquitous across tropical urban settings in the Western Hemisphere. Its ability to thrive in close proximity to humans has been linked to the rapid disappearance of native geckos. However, aspects of Hemidactylus mabouia natural history and ecomorphology, often assumed to be linked with this effect on native populations, remain understudied or untested. Here, we combine data from ∂15N and ∂13C stable isotopes, stomach contents, and morphometric analyses of traits associated with feeding and locomotion to test alternate hypotheses of displacement between H. mabouia and a native gecko, Phyllodactylus martini, on the island of Curaçao. We demonstrate substantial overlap of invertebrate prey resources between the species, with H. mabouia stomachs containing larger arthropod prey as well as vertebrate prey. We additionally show that H. mabouia possesses several morphological advantages, including larger sizes in feeding-associated traits and limb proportions that could offer a propulsive locomotor advantage on vertical surfaces. Together, these findings provide the first support for the hypotheses that invasive H. mabouia and native P. martini overlap in prey resources and that H. mabouia possess ecomorphological advantages over P. martini. This work provides critical context for follow-up studies of H. mabouia and P. martini natural history and direct behavioral experiments that may ultimately illuminate the mechanisms underlying displacement on this island and act as a potential model for other systems with Hemidactylus mabouia invasions.

The mitochondrial phylogeny of an ancient lineage of ray-finned fishes (Polypteridae) with implications for the evolution of body elongation, pelvic fin loss, and craniofacial morphology in Osteichthyes.

BACKGROUND: The family Polypteridae, commonly known as "bichirs", is a lineage that diverged early in the evolutionary history of Actinopterygii (ray-finned fish), but has been the subject of far less evolutionary study than other members of that clade. Uncovering patterns of morphological change within Polypteridae provides an important opportunity to evaluate if the mechanisms underlying morphological evolution are shared among actinoptyerygians, and in fact, perhaps the entire osteichthyan (bony fish and tetrapods) tree of life. However, the greatest impediment to elucidating these patterns is the lack of a well-resolved, highly-supported phylogenetic tree of Polypteridae. In fact, the interrelationships of polypterid species have never been subject to molecular phylogenetic analysis. Here, we infer the first molecular phylogeny of bichirs, including all 12 recognized species and multiple subspecies using Bayesian analyses of 16S and cyt-b mtDNA. We use this mitochondrial phylogeny, ancestral state reconstruction, and geometric morphometrics to test whether patterns of morphological evolution, including the evolution of body elongation, pelvic fin reduction, and craniofacial morphology, are shared throughout the osteichthyan tree of life. RESULTS: Our molecular phylogeny reveals 1) a basal divergence between Erpetoichthys and Polypterus, 2) polyphyly of P. endlicheri and P. palmas, and thus 3) the current taxonomy of Polypteridae masks its underlying genetic diversity. Ancestral state reconstructions suggest that pelvic fins were lost independently in Erpetoichthys, and unambiguously estimate multiple independent derivations of body elongation and shortening. Our mitochondrial phylogeny suggested species that have lower jaw protrusion and up-righted orbit are closely related to each other, indicating a single transformation of craniofacial morphology. CONCLUSION: The mitochondrial phylogeny of polypterid fish provides a strongly-supported phylogenetic framework for future comparative evolutionary, physiological, ecological, and genetic analyses. Indeed, ancestral reconstruction and geometric morphometric analyses revealed that the patterns of morphological evolution in Polypteridae are similar to those seen in other osteichthyans, thus implying the underlying genetic and developmental mechanisms responsible for those patterns were established early in the evolutionary history of Osteichthyes. We propose developmental and genetic mechanisms to be tested under the light of this new phylogenetic framework.

Does dispersal across an aquatic geographic barrier obscure phylogeographic structure in the diamond-backed watersnake (Nerodia rhombifer)?

The impact of barriers to dispersal and gene flow is often inferred to be the primary cause of lineage divergence and phylogeographic structure in terrestrial organisms. In particular, the Mississippi River has been implicated as a barrier to gene flow in many species, including aquatic taxa. However, if barriers are permeable to organisms, then phylogeographic structure may be difficult to detect due to gene flow between lineages. Using time-calibrated Bayesian phylogenetic analyses of mtDNA, and phylogeographic coalescent simulations, we determine if the Mississippi River operates as a barrier to gene flow in the aquatic diamond-backed watersnake (Nerodia rhombifer). The phylogenetic analyses support a basal division within N. rhombifer mtDNA lineages that coincides with populations generally east and west of the Mississippi River. These results, and that of the divergence dating analyses, therefore suggest that the river was a significant barrier to gene flow in the Pleistocene ∼ 1.4 million years ago, presumably during an interglacial period when the river was much wider. However, we also detect western haplotypes in the eastern clade, and vice versa, thereby indicating that this barrier has not been complete. Nonetheless, the coalescent simulations that account for limited migration suggest that the Mississippi River was an important feature that shaped the phylogeographic history of this aquatic snake in the USA despite limited gene flow.