Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish.

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

High germline mutation rates, but not extreme population outbreaks, influence genetic diversity in a keystone coral predator.

Lewontin's paradox, the observation that levels of genetic diversity (π) do not scale linearly with census population size (Nc) variation, is an evolutionary conundrum. The most extreme mismatches between π and Nc are found for highly abundant marine invertebrates. Yet, the influences of new mutations on π relative to extrinsic processes such as Nc fluctuations are unknown. Here, we provide the first germline mutation rate (µ) estimate for a marine invertebrate in corallivorous crown-of-thorns sea stars (Acanthaster cf. solaris). We use high-coverage whole-genome sequencing of 14 parent-offspring trios alongside empirical estimates of Nc in Australia's Great Barrier Reef to jointly examine the determinants of π in populations undergoing extreme Nc fluctuations. The A. cf. solaris mean µ was 9.13 x 10-09 mutations per-site per-generation (95% CI: 6.51 x 10-09 to 1.18 x 10-08), exceeding estimates for other invertebrates and showing greater concordance with vertebrate mutation rates. Lower-than-expected Ne (~70,000-180,000) and low Ne/Nc values (0.0047-0.048) indicated weak influences of population outbreaks on long-term π. Our findings are consistent with elevated µ evolving in response to reduced Ne and generation time length, with important implications for explaining high mutational loads and the determinants of genetic diversity in marine invertebrate taxa.

Exploring coral speciation: Multiple sympatric Stylophora pistillata taxa along a divergence continuum on the Great Barrier Reef.

Understanding how biodiversity originates and is maintained are fundamental challenge in evolutionary biology. Speciation is a continuous process and progression along this continuum depends on the interplay between evolutionary forces driving divergence and forces promoting genetic homogenisation. Coral reefs are broadly connected yet highly heterogeneous ecosystems, and divergence with gene flow at small spatial scales might therefore be common. Genomic studies are increasingly revealing the existence of closely related and sympatric taxa within taxonomic coral species, but the extent to which these taxa might still be exchanging genes and sharing environmental niches is unclear. In this study, we sampled extensively across diverse habitats at multiple reefs of the Great Barrier Reef (GBR) and comprehensively examined genome-wide diversity and divergence histories within and among taxa of the Stylophora pistillata species complex. S. pistillata is one of the most abundant and well-studied coral species, yet we discovered five distinct taxa, with wide geographic ranges and extensive sympatry. Demographic modelling showed that speciation events have occurred with gene flow and that taxa are at different stages along a divergence continuum. We found significant correlations between genetic divergence and specific environmental variables, suggesting that niche partitioning may have played a role in speciation and that S. pistillata taxa might be differentially adapted to different environments. Conservation actions rely on estimates of species richness, population sizes and species ranges, which are biased if divergent taxa are lumped together. As coral reefs are rapidly degrading due to climate change, our study highlights the importance of recognising evolutionarily distinct and differentially adapted coral taxa to improve conservation and restoration efforts aiming at protecting coral genetic diversity.

Comparative landscape genomics has arrived with a splash.

Comparative methods are central for understanding the distribution of biodiversity. Assessing spatial variation of multiple species can identify biodiversity drivers across landscapes, including on the genetic level. Molecular ecology approaches have expanded in spatial precision, from phylogeography to landscape genetics, as have molecular tools that now allow genome-scale inferences for virtually any organism. These welcome increases in genomic data richness, however, have not been matched by growth in the multispecies dimension, and empirical studies continue to focus predominantly on single species. In this issue of Molecular Ecology, Zbinden et al. (Molecular Ecology, 2022, 32) present a landmark achievement for comparative landscape genomics, surveying an impressive 31 species of fishes, sampled from 75 locations in the White River Basin and genotyping each species for thousands of SNPs. Zbinden et al. (Molecular Ecology, 2022, 32) draw upon the statistical tools of landscape genetics to comprehensively interrogate four hypotheses-that populations of multiple freshwater fish species are isolated by river distance, barriers to dispersal, stream hierarchy or environment. Stream hierarchy overwhelmingly predicts population structure within co-distributed fish species pointing to commonalities among species at the large landscape scale (100s of km). The tight alignment between intraspecific genetic spatial structure and stream hierarchies, moreover, provides clear validation for conservation and fisheries management to use watershed divisions as distinct management units.

Moving beyond heritability in the search for coral adaptive potential.

Global environmental change is happening at unprecedented rates. Coral reefs are among the ecosystems most threatened by global change. For wild populations to persist, they must adapt. Knowledge shortfalls about corals' complex ecological and evolutionary dynamics, however, stymie predictions about potential adaptation to future conditions. Here, we review adaptation through the lens of quantitative genetics. We argue that coral adaptation studies can benefit greatly from "wild" quantitative genetic methods, where traits are studied in wild populations undergoing natural selection, genomic relationship matrices can replace breeding experiments, and analyses can be extended to examine genetic constraints among traits. In addition, individuals with advantageous genotypes for anticipated future conditions can be identified. Finally, genomic genotyping supports simultaneous consideration of how genetic diversity is arrayed across geographic and environmental distances, providing greater context for predictions of phenotypic evolution at a metapopulation scale.

Research priorities for the sustainability of coral-rich western Pacific seascapes.

Nearly a billion people depend on tropical seascapes. The need to ensure sustainable use of these vital areas is recognised, as one of 17 policy commitments made by world leaders, in Sustainable Development Goal (SDG) 14 ('Life below Water') of the United Nations. SDG 14 seeks to secure marine sustainability by 2030. In a time of increasing social-ecological unpredictability and risk, scientists and policymakers working towards SDG 14 in the Asia-Pacific region need to know: (1) How are seascapes changing? (2) What can global society do about these changes? and (3) How can science and society together achieve sustainable seascape futures? Through a horizon scan, we identified nine emerging research priorities that clarify potential research contributions to marine sustainability in locations with high coral reef abundance. They include research on seascape geological and biological evolution and adaptation; elucidating drivers and mechanisms of change; understanding how seascape functions and services are produced, and how people depend on them; costs, benefits, and trade-offs to people in changing seascapes; improving seascape technologies and practices; learning to govern and manage seascapes for all; sustainable use, justice, and human well-being; bridging communities and epistemologies for innovative, equitable, and scale-crossing solutions; and informing resilient seascape futures through modelling and synthesis. Researchers can contribute to the sustainability of tropical seascapes by co-developing transdisciplinary understandings of people and ecosystems, emphasising the importance of equity and justice, and improving knowledge of key cross-scale and cross-level processes, feedbacks, and thresholds.

Reproducibility in ecology and evolution: Minimum standards for data and code.

We call for journals to commit to requiring open data be archived in a format that will be simple and clear for readers to understand and use. If applied consistently, these requirements will allow contributors to be acknowledged for their work through citation of open data, and facilitate scientific progress.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef.

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Importance of timely metadata curation to the global surveillance of genetic diversity.

Genetic diversity within species represents a fundamental yet underappreciated level of biodiversity. Because genetic diversity can indicate species resilience to changing climate, its measurement is relevant to many national and global conservation policy targets. Many studies produce large amounts of genome-scale genetic diversity data for wild populations, but most (87%) do not include the associated spatial and temporal metadata necessary for them to be reused in monitoring programs or for acknowledging the sovereignty of nations or Indigenous peoples. We undertook a distributed datathon to quantify the availability of these missing metadata and to test the hypothesis that their availability decays with time. We also worked to remediate missing metadata by extracting them from associated published papers, online repositories, and direct communication with authors. Starting with 848 candidate genomic data sets (reduced representation and whole genome) from the International Nucleotide Sequence Database Collaboration, we determined that 561 contained mostly samples from wild populations. We successfully restored spatiotemporal metadata for 78% of these 561 data sets (n = 440 data sets with data on 45,105 individuals from 762 species in 17 phyla). Examining papers and online repositories was much more fruitful than contacting 351 authors, who replied to our email requests 45% of the time. Overall, 23% of our email queries to authors unearthed useful metadata. The probability of retrieving spatiotemporal metadata declined significantly as age of the data set increased. There was a 13.5% yearly decrease in metadata associated with published papers or online repositories and up to a 22% yearly decrease in metadata that were only available from authors. This rapid decay in metadata availability, mirrored in studies of other types of biological data, should motivate swift updates to data-sharing policies and researcher practices to ensure that the valuable context provided by metadata is not lost to conservation science forever.

Importancia de la curación oportuna de metadatos para la vigilancia mundial de la diversidad genética Resumen La diversidad genética intraespecífica representa un nivel fundamental, pero a la vez subvalorado de la biodiversidad. La diversidad genética puede indicar la resiliencia de una especie ante el clima cambiante, por lo que su medición es relevante para muchos objetivos de la política de conservación mundial y nacional. Muchos estudios producen una gran cantidad de datos sobre la diversidad a nivel genético de las poblaciones silvestres, aunque la mayoría (87%) no incluye los metadatos espaciales y temporales asociados para que sean reutilizados en los programas de monitoreo o para reconocer la soberanía de las naciones o los pueblos indígenas. Realizamos un "datatón" distribuido para cuantificar la disponibilidad de estos metadatos faltantes y para probar la hipótesis que supone que esta disponibilidad se deteriora con el tiempo. También trabajamos para reparar los metadatos faltantes al extraerlos de los artículos asociados publicados, los repositorios en línea y la comunicación directa con los autores. Iniciamos con 838 candidatos de conjuntos de datos genómicos (representación reducida y genoma completo) tomados de la colaboración internacional para la base de datos de secuencias de nucleótidos y determinamos que 561 incluían en su mayoría muestras tomadas de poblaciones silvestres. Restauramos con éxito los metadatos espaciotemporales en el 78% de estos 561 conjuntos de datos (n = 440 conjuntos de datos con información sobre 45,105 individuos de 762 especies en 17 filos). El análisis de los artículos y los repositorios virtuales fue mucho más productivo que contactar a los 351 autores, quienes tuvieron un 45% de respuesta a nuestros correos. En general, el 23% de nuestras consultas descubrieron metadatos útiles. La probabilidad de recuperar metadatos espaciotemporales declinó de manera significativa conforme incrementó la antigüedad del conjunto de datos. Hubo una disminución anual del 13.5% en los metadatos asociados con los artículos publicados y los repositorios virtuales y hasta una disminución anual del 22% en los metadatos que sólo estaban disponibles mediante la comunicación con los autores. Este rápido deterioro en la disponibilidad de los metadatos, duplicado en estudios de otros tipos de datos biológicos, debería motivar la pronta actualización de las políticas del intercambio de datos y las prácticas de los investigadores para asegurar que en las ciencias de la conservación no se pierda para siempre el contexto valioso proporcionado por los metadatos.

Influence of offshore oil and gas structures on seascape ecological connectivity.

Offshore platforms, subsea pipelines, wells and related fixed structures supporting the oil and gas (O&G) industry are prevalent in oceans across the globe, with many approaching the end of their operational life and requiring decommissioning. Although structures can possess high ecological diversity and productivity, information on how they interact with broader ecological processes remains unclear. Here, we review the current state of knowledge on the role of O&G infrastructure in maintaining, altering or enhancing ecological connectivity with natural marine habitats. There is a paucity of studies on the subject with only 33 papers specifically targeting connectivity and O&G structures, although other studies provide important related information. Evidence for O&G structures facilitating vertical and horizontal seascape connectivity exists for larvae and mobile adult invertebrates, fish and megafauna; including threatened and commercially important species. The degree to which these structures represent a beneficial or detrimental net impact remains unclear, is complex and ultimately needs more research to determine the extent to which natural connectivity networks are conserved, enhanced or disrupted. We discuss the potential impacts of different decommissioning approaches on seascape connectivity and identify, through expert elicitation, critical knowledge gaps that, if addressed, may further inform decision making for the life cycle of O&G infrastructure, with relevance for other industries (e.g. renewables). The most highly ranked critical knowledge gap was a need to understand how O&G structures modify and influence the movement patterns of mobile species and dispersal stages of sessile marine species. Understanding how different decommissioning options affect species survival and movement was also highly ranked, as was understanding the extent to which O&G structures contribute to extending species distributions by providing rest stops, foraging habitat, and stepping stones. These questions could be addressed with further dedicated studies of animal movement in relation to structures using telemetry, molecular techniques and movement models. Our review and these priority questions provide a roadmap for advancing research needed to support evidence-based decision making for decommissioning O&G infrastructure.

Deep connections: Divergence histories with gene flow in mesophotic Agaricia corals.

Largely understudied, mesophotic coral ecosystems lie below shallow reefs (at >30 m depth) and comprise ecologically distinct communities. Brooding reproductive modes appear to predominate among mesophotic-specialist corals and may limit genetic connectivity among populations. Using reduced representation genomic sequencing, we assessed spatial population genetic structure at 50 m depth in an ecologically important mesophotic-specialist species Agaricia grahamae, among locations in the Southern Caribbean. We also tested for hybridisation with the closely related (but depth-generalist) species Agaricia lamarcki, within their sympatric depth zone (50 m). In contrast to our expectations, no spatial genetic structure was detected between the reefs of Curaçao and Bonaire (~40 km apart) within A. grahamae. However, cryptic taxa were discovered within both taxonomic species, with those in A. lamarcki (incompletely) partitioned by depth and those in A. grahamae occurring sympatrically (at the same depth). Hybrid analyses and demographic modelling identified contemporary and historical gene flow among cryptic taxa, both within and between A. grahamae and A. lamarcki. These results (1) indicate that spatial connectivity and subsequent replenishment may be possible between islands of moderate geographic distances for A. grahamae, an ecologically important mesophotic species, (2) that cryptic taxa occur in the mesophotic zone and environmental selection along shallow to mesophotic depth gradients may drive divergence in depth-generalists such as A. lamarcki, and (3) highlight that gene flow links taxa within this relativity diverse Caribbean genus.

Poor data stewardship will hinder global genetic diversity surveillance.

Genomic data are being produced and archived at a prodigious rate, and current studies could become historical baselines for future global genetic diversity analyses and monitoring programs. However, when we evaluated the potential utility of genomic data from wild and domesticated eukaryote species in the world's largest genomic data repository, we found that most archived genomic datasets (86%) lacked the spatiotemporal metadata necessary for genetic biodiversity surveillance. Labor-intensive scouring of a subset of published papers yielded geospatial coordinates and collection years for only 33% (39% if place names were considered) of these genomic datasets. Streamlined data input processes, updated metadata deposition policies, and enhanced scientific community awareness are urgently needed to preserve these irreplaceable records of today's genetic biodiversity and to plug the growing metadata gap.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle.

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

Genetic and phenotypic variation exhibit both predictable and stochastic patterns across an intertidal fish metapopulation.

Interactions among selection, gene flow, and drift affect the trajectory of adaptive evolution. In natural populations, the direction and magnitude of these processes can be variable across different spatial, temporal, or ontogenetic scales. Consequently, variability in evolutionary processes affects the predictability or stochasticity of microevolutionary outcomes. We studied an intertidal fish, Bathygobius cocosensis (Bleeker, 1854), to understand how space, time, and life stage structure genetic and phenotypic variation in a species with potentially extensive dispersal and a complex life cycle (larval dispersal preceding benthic recruitment). We sampled juvenile and adult life stages, at three sites, over three years. Genome-wide SNPs uncovered a pattern of chaotic genetic patchiness, that is, weak-but-significant patchy spatial genetic structure that was variable through time and between life stages. Outlier locus analyses suggested that targets of spatially divergent selection were mostly temporally variable, though a significant number of spatial outlier loci were shared between life stages. Head shape, a putatively ecologically responsive (adaptive) phenotype in B. cocosensis also exhibited high temporal variability within sites. However, consistent spatial relationships between sites indicated that environmental similarities among sites may generate predictable phenotype distributions across space. Our study highlights the complex microevolutionary dynamics of marine systems, where consideration of multiple ecological dimensions can reveal both predictable and stochastic patterns in the distributions of genetic and phenotypic variation. Such considerations probably apply to species that possess short, complex life cycles, have large dispersal potential and fecundities, and that inhabit heterogeneous environments.

Skeletal deformities and meristic trait variations are common in the intertidal fish Bathygobius cocosensis (Perciformes-Gobiidae).

Extreme environmental conditions, such as temperature, can lead to meristic trait variation and skeletal deformities, which may have major impacts on individual fitness. As intertidal ecosystems experience rapid temperature and physicochemical fluctuations, intertidal fish living and reproducing in these environmental conditions may have phenotypes influenced by such variable environments. The impact of intertidal variability on fish development, however, has not been previously investigated. Skeletal deformities and meristic traits were assessed for Bathygobius cocosensis, a common intertidal fish living across the Indo-Pacific region, using a clearing and staining method on 72 individuals. Over 87% of individuals presented meristic variation and over 70% exhibited at least one type of skeletal deformity, mostly recorded in the caudal fin area. The unexpected prevalence of skeletal deformities among this intertidal fish population suggests that such deformities may be suitable markers to evaluate an individual's stress exposure during development and the subsequent fitness effects.

Pre-introduction introgression contributes to parallel differentiation and contrasting hybridization outcomes between invasive and native marine mussels.

Non-native species experience novel selection pressures in introduced environments and may interbreed with native lineages. Species introductions therefore provide opportunities to investigate repeated patterns of adaptation and introgression across replicated contact zones. Here, we investigate genetic parallelism between multiple introduced populations of the invasive marine mussel, Mytilus galloprovincialis, in the absence (South Africa and California) and presence of hybridization with a native congener (Mytilus planulatus in Batemans Bay and Sydney Harbour, Australia). Repeatability in post-introduction differentiation from native-range populations varied between genetically distinct Atlantic and Mediterranean lineages, with Atlantic-derived introductions displaying high differentiation (maxFST  > 0.4) and parallelism at outlier loci. Identification of long noncoding RNA transcripts (lncRNA) additionally allowed us to clarify that parallel responses are largely limited to protein-coding loci, with lncRNAs likely evolving under evolutionary constraints. Comparisons of independent hybrid zones revealed differential introgression most strongly in Batemans Bay, with an excess of M. galloprovincialis ancestry and resistance to introgression at loci differentiating parental lineages (M. planulatus and Atlantic M. galloprovincialis). Additionally, contigs putatively introgressed with divergent alleles from a closely related species, Mytilus edulis, showed stronger introgression asymmetries compared with genome-wide trends and also diverged in parallel in both Atlantic-derived introductions. These results suggest that divergent demographic histories experienced by introduced lineages, including pre-introduction introgression, influence contemporary admixture dynamics. Our findings build on previous investigations reporting contributions of historical introgression to intrinsic reproductive architectures shared between marine lineages and illustrate that interspecific introgression history can shape differentiation between colonizing populations and their hybridization with native congeners.

Building a global genomics observatory: Using GEOME (the Genomic Observatories Metadatabase) to expedite and improve deposition and retrieval of genetic data and metadata for biodiversity research.

Genetic data represent a relatively new frontier for our understanding of global biodiversity. Ideally, such data should include both organismal DNA-based genotypes and the ecological context where the organisms were sampled. Yet most tools and standards for data deposition focus exclusively either on genetic or ecological attributes. The Genomic Observatories Metadatabase (GEOME: geome-db.org) provides an intuitive solution for maintaining links between genetic data sets stored by the International Nucleotide Sequence Database Collaboration (INSDC) and their associated ecological metadata. GEOME facilitates the deposition of raw genetic data to INSDCs sequence read archive (SRA) while maintaining persistent links to standards-compliant ecological metadata held in the GEOME database. This approach facilitates findable, accessible, interoperable and reusable data archival practices. Moreover, GEOME enables data management solutions for large collaborative groups and expedites batch retrieval of genetic data from the SRA. The article that follows describes how GEOME can enable genuinely open data workflows for researchers in the field of molecular ecology.

Anthropogenic hybridization at sea: three evolutionary questions relevant to invasive species management.

Species introductions promote secondary contacts between taxa with long histories of allopatric divergence. Anthropogenic contact zones thus offer valuable contrasts to speciation studies in natural systems where past spatial isolations may have been brief or intermittent. Investigations of anthropogenic hybridization are rare for marine animals, which have high fecundity and high dispersal ability, characteristics that contrast to most terrestrial animals. Genomic studies indicate that gene flow can still occur after millions of years of divergence, as illustrated by invasive mussels and tunicates. In this context, we highlight three issues: (i) the effects of high propagule pressure and demographic asymmetries on introgression directionality, (ii) the role of hybridization in preventing introduced species spread, and (iii) the importance of postzygotic barriers in maintaining reproductive isolation. Anthropogenic contact zones offer evolutionary biologists unprecedented large scale hybridization experiments. In addition to breaking the highly effective reproductive isolating barrier of spatial segregation, they allow researchers to explore unusual demographic contexts with strong asymmetries. The outcomes are diverse, from introgression swamping to strong barriers to gene flow, and lead to local containment or widespread invasion. These outcomes should not be neglected in management policies of marine invasive species. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Twin introductions by independent invader mussel lineages are both associated with recent admixture with a native congener in Australia.

Introduced species can impose profound impacts on the evolution of receiving communities with which they interact. If native and introduced taxa remain reproductively semi-isolated, human-mediated secondary contact may promote genetic exchange across newly created hybrid zones, potentially impacting native genetic diversity and invasive species spread. Here, we investigate the contributions of recent divergence histories and ongoing (post-introduction) gene flow between the invasive marine mussel, Mytilus galloprovincialis, and a morphologically indistinguishable and taxonomically contentious native Australian taxon, Mytilus planulatus. Using transcriptome-wide markers, we demonstrate that two contemporary M. galloprovincialis introductions into south-eastern Australia originate from genetically divergent lineages from its native range in the Mediterranean Sea and Atlantic Europe, where both introductions have led to repeated instances of admixture between introduced and endemic populations. Through increased genome-wide resolution of species relationships, combined with demographic modelling, we validate that mussels sampled in Tasmania are representative of the endemic Australian taxon (M. planulatus), but share strong genetic affinities to M. galloprovincialis. Demographic inferences indicate late-Pleistocene divergence times and historical gene flow between the Tasmanian endemic lineage and northern M. galloprovincialis, suggesting that native and introduced taxa have experienced a period of historical isolation of at least 100,000 years. Our results demonstrate that many genomic loci and sufficient sampling of closely related lineages in both sympatric (e.g. Australian populations) and allopatric (e.g. northern hemisphere Mytilus taxa) ranges are necessary to accurately (a) interpret patterns of intraspecific differentiation and to (b) distinguish contemporary invasive introgression from signatures left by recent divergence histories in high dispersal marine species. More broadly, our study fills a significant gap in systematic knowledge of native Australian biodiversity and sheds light on the intrinsic challenges for invasive species research when native and introduced species boundaries are not well defined.