Synergistic innovations enabled the radiation of anglerfishes in the deep open ocean.

Major ecological transitions are thought to fuel diversification, but whether they are contingent on the evolution of certain traits called key innovations1 is unclear. Key innovations are routinely invoked to explain how lineages rapidly exploit new ecological opportunities.1,2,3 However, investigations of key innovations often focus on single traits rather than considering trait combinations that collectively produce effects of interest.4 Here, we investigate the evolution of synergistic trait interactions in anglerfishes, which include one of the most species-rich vertebrate clades in the bathypelagic, or "midnight," zone of the deep sea: Ceratioidea.5 Ceratioids are the only vertebrates that possess sexual parasitism, wherein males temporarily attach or permanently fuse to females to mate.6,7 We show that the rapid transition of ancestrally benthic anglerfishes into pelagic habitats occurred during a period of major global warming 50-35 million years ago.8,9 This transition coincided with the origins of sexual parasitism, which is thought to increase the probability of successful reproduction once a mate is found in the midnight zone, Earth's largest habitat.5,6,7 Our reconstruction of the evolutionary history of anglerfishes and the loss of immune genes support that permanently fusing clades have convergently degenerated their adaptive immunity. We find that degenerate adaptive immune genes and sexual body size dimorphism, both variably present in anglerfishes outside the ceratioid radiation, likely promoted their transition into the bathypelagic zone. These results show how traits from separate physiological, morphological, and reproductive systems can interact synergistically to drive major transitions and subsequent diversification in novel environments.

Considering decoupled phenotypic diversification between ontogenetic phases in macroevolution: An example using Triggerfishes (Balistidae).

Across the Tree of Life, most studies of phenotypic disparity and diversification have been restricted to adult organisms. However, many lineages have distinct ontogenetic phases that differ from their adult forms in morphology and ecology. Focusing disproportionately on the evolution of adult forms unnecessarily hinders our understanding of the pressures shaping evolution over time. Non-adult disparity patterns are particularly important to consider for coastal ray-finned fishes, which often have juvenile phases with distinct phenotypes. These juvenile forms are often associated with sheltered nursery environments, with phenotypic shifts between adults and juvenile stages that are readily apparent in locomotor morphology. Whether this ontogenetic variation in locomotor morphology reflects a decoupling of diversification dynamics between life stages remains unknown. Here we investigate the evolutionary dynamics of locomotor morphology between adult and juvenile triggerfishes. We integrate a time-calibrated phylogenetic framework with geometric morphometric approaches and measurement data of fin aspect ratio and incidence, and reveal a mismatch between morphospace occupancy, the evolution of morphological disparity, and the tempo of trait evolution between life stages. Collectively, our results illuminate how the heterogeneity of morpho-functional adaptations can decouple the mode and tempo of morphological diversification between ontogenetic stages.

Seasonality of endemic COVID-19.

Successive waves of infection by SARS-CoV-2 have left little doubt that this virus will transition to an endemic disease. Foreknowledge of when to expect seasonal surges is crucial for healthcare and public health decision-making. However, the future seasonality of COVID-19 remains uncertain. Evaluating its seasonality is complicated due to the limited years of SARS-CoV-2 circulation, pandemic dynamics, and varied interventions. In this study, we project the expected endemic seasonality by employing a phylogenetic ancestral and descendant state approach that leverages long-term data on the incidence of circulating HCoV coronaviruses. Our projections indicate asynchronous surges of SARS-CoV-2 across different locations in the northern hemisphere, occurring between October and January in New York and between January and March in Yamagata, Japan. This knowledge of spatiotemporal surges leads to medical preparedness and enables the implementation of targeted public health interventions to mitigate COVID-19 transmission.IMPORTANCEThe seasonality of COVID-19 is important for effective healthcare and public health decision-making. Previous waves of SARS-CoV-2 infections have indicated that the virus will likely persist as an endemic pathogen with distinct surges. However, the timing and patterns of potentially seasonal surges remain uncertain, rendering effective public health policies uninformed and in danger of poorly anticipating opportunities for intervention, such as well-timed booster vaccination drives. Applying an evolutionary approach to long-term data on closely related circulating coronaviruses, our research provides projections of seasonal surges that should be expected at major temperate population centers. These projections enable local public health efforts that are tailored to expected surges at specific locales or regions. This knowledge is crucial for enhancing medical preparedness and facilitating the implementation of targeted public health interventions.

Infection with alternate frequencies of SARS-CoV-2 vaccine boosting for patients undergoing antineoplastic cancer treatments.

Patients undergoing antineoplastic therapies often exhibit reduced immune response to COVID-19 vaccination, necessitating assessment of alternate booster vaccination frequencies. However, data on reinfection risks to guide clinical decision making are limited. Here, we quantified reinfection risks for patients undergoing distinct antineoplastic therapies, given alternative frequencies of boosting with Pfizer-BioNTech BNT162b2. Integrating antibody data following vaccination with long-term antibody data from other coronaviruses in an evolutionary framework, we estimated infection probabilities based on antibody levels and calculated cumulative probabilities of breakthrough infection for alternate booster schedules over 2 years. Annual boosting reduced risks for targeted or hormonal treatments, immunotherapy, and chemotherapy-immunotherapy combinations similarly to the general population. Patients receiving no treatment or chemotherapy exhibited higher risks, suggesting that accelerated vaccination schedules should be considered. Patients treated with rituximab therapy presented the highest infection risk, suggesting that a combination of frequent boosting and additional interventions may be warranted for mitigating SARS-CoV-2 infection.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution.

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

A chromosome-level genome assembly of longnose gar, Lepisosteus osseus.

Holosteans (gars and bowfins) represent the sister lineage to teleost fishes, the latter being a clade that comprises over half of all living vertebrates and includes important models for comparative genomics and human health. A major distinction between the evolutionary history of teleosts and holosteans is that all teleosts experienced a genome duplication event in their early evolutionary history. As the teleost genome duplication occurred after teleosts diverged from holosteans, holosteans have been heralded as a means to bridge teleost models to other vertebrate genomes. However, only three species of holosteans have been genome-sequenced to date, and sequencing of more species is needed to fill sequence sampling gaps and provide a broader comparative basis for understanding holostean genome evolution. Here we report the first high quality reference genome assembly and annotation of the longnose gar (Lepisosteus osseus). Our final assembly consists of 22,709 scaffolds with a total length of 945 bp with contig N50 of 116.61 kb. Using BRAKER2, we annotated a total of 30,068 genes. Analysis of the repetitive regions of the genome reveals the genome to contain 29.12% transposable elements, and the longnose gar to be the only other known vertebrate outside of the spotted gar and bowfin to contain CR1, L2, Rex1, and Babar. These results highlight the potential utility of holostean genomes for understanding the evolution of vertebrate repetitive elements, and provide a critical reference for comparative genomic studies utilizing ray-finned fish models.

Infection by SARS-CoV-2 with alternate frequencies of mRNA vaccine boosting.

One of the most consequential unknowns of the COVID-19 pandemic is the frequency at which vaccine boosting provides sufficient protection from infection. We quantified the statistical likelihood of breakthrough infections over time following different boosting schedules with messenger RNA (mRNA)-1273 (Moderna) and BNT162b2 (Pfizer-BioNTech). We integrated anti-Spike IgG antibody optical densities with profiles of the waning of antibodies and corresponding probabilities of infection associated with coronavirus endemic transmission. Projecting antibody levels over time given boosting every 6 months, 1, 1.5, 2, or 3 years yielded respective probabilities of fending off infection over a 6-year span of >93%, 75%, 55%, 40%, and 24% (mRNA-1273) and >89%, 69%, 49%, 36%, and 23% (BNT162b2). Delaying the administration of updated boosters has bleak repercussions. It increases the probability of individual infection by SARS-CoV-2, and correspondingly, ongoing disease spread, prevalence, morbidity, hospitalization, and mortality. Instituting regular, population-wide booster vaccination updated to predominant variants has the potential to substantially forestall-and with global, widespread uptake, eliminate-COVID-19.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors.

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Placing human gene families into their evolutionary context.

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Periodic Environmental Disturbance Drives Repeated Ecomorphological Diversification in an Adaptive Radiation of Antarctic Fishes.

AbstractThe ecological theory of adaptive radiation has profoundly shaped our conceptualization of the rules that govern diversification. However, while many radiations follow classic early-burst patterns of diversification as they fill ecological space, the longer-term fates of these radiations depend on many factors, such as climatic stability. In systems with periodic disturbances, species-rich clades can contain nested adaptive radiations of subclades with their own distinct diversification histories, and how adaptive radiation theory applies in these cases is less clear. Here, we investigated patterns of ecological and phenotypic diversification within two iterative adaptive radiations of cryonotothenioid fishes in Antarctica's Southern Ocean: crocodile icefishes and notoperches. For both clades, we observe evidence of repeated diversification into disparate regions of trait space between closely related taxa and into overlapping regions of trait space between distantly related taxa. We additionally find little evidence that patterns of ecological divergence are correlated with evolution of morphological disparity, suggesting that these axes of divergence may not be tightly linked. Finally, we reveal evidence of repeated convergence in sympatry that suggests niche complementarity. These findings reflect the dynamic history of Antarctic marine habitats and may guide hypotheses of diversification dynamics in environments characterized by periodic disturbance.

Prolonged morphological expansion of spiny-rayed fishes following the end-Cretaceous.

Spiny-rayed fishes (Acanthomorpha) dominate modern marine habitats and account for more than a quarter of all living vertebrate species. Previous time-calibrated phylogenies and patterns from the fossil record explain this dominance by correlating the origin of major acanthomorph lineages with the Cretaceous-Palaeogene mass extinction. Here we infer a time-calibrated phylogeny using ultraconserved elements that samples 91.4% of all acanthomorph families and investigate patterns of body shape disparity. Our results show that acanthomorph lineages steadily accumulated throughout the Cenozoic and underwent a significant expansion of among-clade morphological disparity several million years after the end-Cretaceous. These acanthomorph lineages radiated into and diversified within distinct regions of morphospace that characterize iconic lineages, including fast-swimming open-ocean predators, laterally compressed reef fishes, bottom-dwelling flatfishes, seahorses and pufferfishes. The evolutionary success of spiny-rayed fishes is the culmination of multiple species-rich and phenotypically disparate lineages independently diversifying across the globe under a wide range of ecological conditions.

The durability of natural infection and vaccine-induced immunity against future infection by SARS-CoV-2.

The durability of vaccine-mediated immunity to SARS-CoV-2, the durations to breakthrough infection, and the optimal timings of booster vaccination are crucial knowledge for pandemic response. Here, we applied comparative evolutionary analyses to estimate the durability of immunity and the likelihood of breakthrough infections over time following vaccination by BNT162b2 (Pfizer-BioNTech), mRNA-1273 (Moderna), ChAdOx1 (Oxford-AstraZeneca), and Ad26.COV2.S (Johnson & Johnson/Janssen). We evaluated anti-Spike (S) immunoglobulin G (IgG) antibody levels elicited by each vaccine relative to natural infection. We estimated typical trajectories of waning and corresponding infection probabilities, providing the distribution of times to breakthrough infection for each vaccine under endemic conditions. Peak antibody levels elicited by messenger RNA (mRNA) vaccines mRNA-1273 and BNT1262b2 exceeded that of natural infection and are expected to typically yield more durable protection against breakthrough infections (median 29.6 mo; 5 to 95% quantiles 10.9 mo to 7.9 y) than natural infection (median 21.5 mo; 5 to 95% quantiles 3.5 mo to 7.1 y). Relative to mRNA-1273 and BNT1262b2, viral vector vaccines ChAdOx1 and Ad26.COV2.S exhibit similar peak anti-S IgG antibody responses to that from natural infection and are projected to yield lower, shorter-term protection against breakthrough infection (median 22.4 mo and 5 to 95% quantiles 4.3 mo to 7.2 y; and median 20.5 mo and 5 to 95% quantiles 2.6 mo to 7.0 y; respectively). These results leverage the tools from evolutionary biology to provide a quantitative basis for otherwise unknown parameters that are fundamental to public health policy decision-making.

On the relationship between extant innate immune receptors and the evolutionary origins of jawed vertebrate adaptive immunity.

For over half a century, deciphering the origins of the genomic loci that form the jawed vertebrate adaptive immune response has been a major topic in comparative immunogenetics. Vertebrate adaptive immunity relies on an extensive and highly diverse repertoire of tandem arrays of variable (V), diversity (D), and joining (J) gene segments that recombine to produce different immunoglobulin (Ig) and T cell receptor (TCR) genes. The current consensus is that a recombination-activating gene (RAG)-like transposon invaded an exon of an ancient innate immune VJ-bearing receptor, giving rise to the extant diversity of Ig and TCR loci across jawed vertebrates. However, a model for the evolutionary relationships between extant non-recombining innate immune receptors and the V(D)J receptors of the jawed vertebrate adaptive immune system has only recently begun to come into focus. In this review, we provide an overview of non-recombining VJ genes, including CD8ß, CD79b, natural cytotoxicity receptor 3 (NCR3/NKp30), putative remnants of an antigen receptor precursor (PRARPs), and the multigene family of signal-regulatory proteins (SIRPs), that play a wide range of roles in immune function. We then focus in detail on the VJ-containing novel immune-type receptors (NITRs) from ray-finned fishes, as recent work has indicated that these genes are at least 50 million years older than originally thought. We conclude by providing a conceptual model of the evolutionary origins and phylogenetic distribution of known VJ-containing innate immune receptors, highlighting opportunities for future comparative research that are empowered by this emerging evolutionary perspective.

Transcriptome annotation reveals minimal immunogenetic diversity among Wyoming toads, Anaxyrus baxteri.

Briefly considered extinct in the wild, the future of the Wyoming toad (Anaxyrus baxteri) continues to rely on captive breeding to supplement the wild population. Given its small natural geographic range and history of rapid population decline at least partly due to fungal disease, investigation of the diversity of key receptor families involved in the host immune response represents an important conservation need. Population decline may have reduced immunogenetic diversity sufficiently to increase the vulnerability of the species to infectious diseases. Here we use comparative transcriptomics to examine the diversity of toll-like receptors and major histocompatibility complex (MHC) sequences across three individual Wyoming toads. We find reduced diversity at MHC genes compared to bufonid species with a similar history of bottleneck events. Our data provide a foundation for future studies that seek to evaluate the genetic diversity of Wyoming toads, identify biomarkers for infectious disease outcomes, and guide breeding strategies to increase genomic variability and wild release successes.

Premetastatic shifts of endogenous and exogenous mutational processes support consolidative therapy in EGFR-driven lung adenocarcinoma.

The progression of cancer is an evolutionary process that is challenging to assess between sampling timepoints. However, investigation of cancer evolution over specific time periods is crucial to the elucidation of key events such as the acquisition of therapeutic resistance and subsequent fatal metastatic spread of therapy-resistant cell populations. Here we apply mutational signature analyses within clinically annotated cancer chronograms to detect and describe the shifting mutational processes caused by both endogenous (e.g. mutator gene mutation) and exogenous (e.g. mutagenic therapeutics) factors between tumor sampling timepoints. In one patient, we find that cisplatin therapy can introduce mutations that confer genetic resistance to subsequent targeted therapy with Erlotinib. In another patient, we trace detection of defective mismatch-repair associated mutational signature SBS3 to the emergence of known driver mutation CTNNB1 S37C. In both of these patients, metastatic lineages emerged from a single ancestral lineage that arose during therapy-a finding that argues for the consideration of local consolidative therapy over other therapeutic approaches in EGFR-positive non-small cell lung cancer. Broadly, these results demonstrate the utility of phylogenetic analysis that incorporates clinical time course and mutational signature deconvolution to inform therapeutic decision making and retrospective assessment of disease etiology.

The durability of immunity against reinfection by SARS-CoV-2: a comparative evolutionary study.

BACKGROUND: Among the most consequential unknowns of the devastating COVID-19 pandemic are the durability of immunity and time to likely reinfection. There are limited direct data on SARS-CoV-2 long-term immune responses and reinfection. The aim of this study is to use data on the durability of immunity among evolutionarily close coronavirus relatives of SARS-CoV-2 to estimate times to reinfection by a comparative evolutionary analysis of related viruses SARS-CoV, MERS-CoV, human coronavirus (HCoV)-229E, HCoV-OC43, and HCoV-NL63. METHODS: We conducted phylogenetic analyses of the S, M, and ORF1b genes to reconstruct a maximum-likelihood molecular phylogeny of human-infecting coronaviruses. This phylogeny enabled comparative analyses of peak-normalised nucleocapsid protein, spike protein, and whole-virus lysate IgG antibody optical density levels, in conjunction with reinfection data on endemic human-infecting coronaviruses. We performed ancestral and descendent states analyses to estimate the expected declines in antibody levels over time, the probabilities of reinfection based on antibody level, and the anticipated times to reinfection after recovery under conditions of endemic transmission for SARS-CoV-2, as well as the other human-infecting coronaviruses. FINDINGS: We obtained antibody optical density data for six human-infecting coronaviruses, extending from 128 days to 28 years after infection between 1984 and 2020. These data provided a means to estimate profiles of the typical antibody decline and probabilities of reinfection over time under endemic conditions. Reinfection by SARS-CoV-2 under endemic conditions would likely occur between 3 months and 5·1 years after peak antibody response, with a median of 16 months. This protection is less than half the duration revealed for the endemic coronaviruses circulating among humans (5-95% quantiles 15 months to 10 years for HCoV-OC43, 31 months to 12 years for HCoV-NL63, and 16 months to 12 years for HCoV-229E). For SARS-CoV, the 5-95% quantiles were 4 months to 6 years, whereas the 95% quantiles for MERS-CoV were inconsistent by dataset. INTERPRETATION: The timeframe for reinfection is fundamental to numerous aspects of public health decision making. As the COVID-19 pandemic continues, reinfection is likely to become increasingly common. Maintaining public health measures that curb transmission-including among individuals who were previously infected with SARS-CoV-2-coupled with persistent efforts to accelerate vaccination worldwide is critical to the prevention of COVID-19 morbidity and mortality. FUNDING: US National Science Foundation.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system.

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

The bowfin genome illuminates the developmental evolution of ray-finned fishes.

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

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.