Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates

Does selection favour the maintenance of porous species boundaries?

The endpoint of speciation has been viewed as complete isolation and the absence of gene flow between species. If the influx of genes from another species is maladaptive because species have different adaptations and genetic backgrounds, selection should favour the closing of species boundaries and zero gene flow, a process known as reinforcement. Recently, numerous cases of gene flow between species have been identified, many of which involved adaptive introgression of beneficial alleles. These cases could reflect transient states on the way to closed species boundaries or the result of declining strength or efficacy of selection for reinforcement as the level of gene flow approaches zero. An alternative hypothesis, however, is that selection favours porous species boundaries that allow beneficial alleles to cross, especially in changing environments. This perspective evaluates the conditions that would be needed for selection to favour porous species boundaries and the evidence for them. A contrast is made between hybridization in sexual eukaryotes and gene transfer via homologous recombination in bacteria. Current evidence is inconclusive on whether non-zero gene flow is favoured by selection. Studies are needed that quantify selection gradients on rates of gene flow and test for evolution towards intermediate values, especially experiments that manipulate conditions and track evolution for multiple generations. Estimation of gene flow networks for more clades and regional assemblages using emerging genome data will also allow the evolutionary determinants of interspecific gene flow to be better understood.

Contrasting population genetic structure of three semi-terrestrial brachyuran crabs on the coast of the Japanese archipelago

Anthropogenic activities have reduced ecotones between the ocean and land, which is likely to threaten the population of brackish-water brachyuran crabs. To assess the current status of these crabs, we examine the population genetic structures of three semi-terrestrial brachyuran crabs widely distributed along the coast of the Japan and to clarify factors determining their genetic structures. We collected 184 Orisarma dehaani, 252 Chiromantes haematocheir, and 151 Helice tridens crabs from 36 localities of the Japanese archipelago. Genome-wide SNP data from these crabs were analyzed using MIG-seq. Bayesian clustering of STRUCTURE and DAPC analysis were used to identify genetically disturbed populations and to visualize genetic differentiation between local populations. Genetic population structure showed clear differentiation between populations on the Pacific coast of the Tohoku region and on other Japanese coasts in O. dehaani, but not in C. haematocheir or H. tridens. The inbreeding coefficient of O. dehaani was significantly higher on the Pacific coast of the Tohoku region compared to other Japanese coasts. C. haematocheir and H. tridens had homogeneous genetic structures along the Japanese coast, but showed genetic differentiation of a local population at their range limits. Thus, O. dehaani showed little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to ocean currents. Although such a regional differentiation was not found in C. haematocheir and H. tridens, one population of C. haematocheir was genetically isolated at the edge of its distribution range and likely vulnerable to environmental changes.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the closure of all U.S. black abalone fisheries since 1993. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS remains unknown. To address these uncertainties, we sequenced and analysed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Outside the inversion, genetic differentiation between sites is minimal and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Demographic inference does indicate a severe population bottleneck beginning just 15 generations in the past, but this decline is short lived, with present-day size far exceeding the pre-bottleneck status quo. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of population genetic structure, uniform diversity and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the species' designation as critically endangered. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS is unknown. To address these uncertainties, we sequenced and analyzed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Genetic divergence between sites is minimal, and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Despite this, however, demographic inference confirms a severe population bottleneck beginning around the time of WS onset, highlighting the temporal offset that may occur between a population collapse and its potential impact on genetic diversity. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of genetic structure, uniform diversity, and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

Parthenogenomics: Insights on mutation rates and nucleotide diversity in parthenogenetic Panagrolaimus nematodes

Asexual reproduction is assumed to lead to the accumulation of deleterious mutations, and reduced heterozygosity due to the absence of recombination. Panagrolaimid nematode species display different modes of reproduction. Sexual reproduction with distinct males and females, asexual reproduction through parthenogenesis in the genus Panagrolaimus, and hermaphroditism in Propanagrolaimus. Here, we compared genomic features of free-living nematodes in populations and species isolated from geographically distant regions to study diversity, and genome-wide differentiation under different modes of reproduction. We firstly estimated genome-wide spontaneous mutation rates in a triploid parthenogenetic Panagrolaimus, and a diploid hermaphroditic Propanagrolaimus via long-term mutation accumulation lines. Secondly, we calculated population genetic parameters including nucleotide diversity, and fixation index (F ST) between populations of asexually and sexually reproducing nematodes. Thirdly, we used phylogenetic network methods on sexually and asexually reproducing Panagrolaimus populations to understand evolutionary relationships between them. The estimated mutation rate was slightly lower for the asexual population, as expected for taxa with this reproductive mode. Natural polyploid asexual populations revealed higher nucleotide diversity. Despite their common ancestor, a gene network revealed a high level of genetic differentiation among asexual populations. The elevated heterozygosity found in the triploid parthenogens could be explained by the third genome copy. Given their tendentially lower mutation rates it can be hypothesized that this is part of the mechanism to evade Muller's ratchet. Our findings in parthenogenetic triploid nematode populations seem to challenge common expectations of evolution under asexuality.

Extensive introgression among strongylocentrotid sea urchins revealed by phylogenomics

Gametic isolation is thought to play an important role in the evolution of reproductive isolation in broadcast-spawning marine invertebrates. However, it is unclear whether gametic isolation commonly evolves early in the speciation process or only accumulates after other reproductive barriers are already in place. It is also unknown whether gametic isolation is an effective barrier to introgression following speciation. Here, we used whole-genome sequencing data and multiple complementary phylogenomic approaches to test whether the well-documented gametic incompatibilities among the strongylocentrotid sea urchins have limited introgression. We quantified phylogenetic discordance, inferred reticulate phylogenetic networks, and applied the Δ statistic using gene tree topologies reconstructed from multiple sequence alignments of protein-coding single-copy orthologs. In addition, we conducted ABBA-BABA tests on genome-wide single nucleotide variants and reconstructed a phylogeny of mitochondrial genomes. Our results revealed strong mito-nuclear discordance and considerable nonrandom gene tree discordance that cannot be explained by incomplete lineage sorting alone. Eight of the nine species examined demonstrated a history of introgression with at least one other species or ancestral lineage, indicating that introgression was common during the diversification of the strongylocentrotid urchins. There was strong support for introgression between four extant species pairs (Strongylocentrotus pallidus ⇔ S. droebachiensis, S. intermedius ⇔ S. pallidus, S. purpuratus ⇔ S. fragilis, and Mesocentrotus franciscanus ⇔ Pseudocentrotus depressus) and additional evidence for introgression on internal branches of the phylogeny. Our results suggest that the existing gametic incompatibilities among the strongylocentrotid urchin species have not been a complete barrier to hybridization and introgression following speciation. Their continued divergence in the face of widespread introgression indicates that other reproductive isolating barriers likely exist and may have been more critical in establishing reproductive isolation early in speciation.

Accumulation of gene copy number variations during the early phase of free-spawning abalone speciation

The genetic basis of speciation in free-spawning marine invertebrates is poorly understood. Although gene copy number variations (GCNVs) and nucleotide variations possibly trigger the speciation of these organisms, empirical evidence for such a hypothesis is limited. In this study, we searched for genomic signatures of GCNVs that may contribute to the speciation of Western Pacific abalone species. Whole-genome sequencing data suggested the existence of significant amounts of GCNVs in closely related abalones, Haliotis discus and H. madaka, in the early phase of speciation. In addition, the degree of interspecies genetic differentiation in the genes where GCNVs were estimated was higher than that in other genes, suggesting that nucleotide divergence also accumulates in the genes with GCNVs. GCNVs in some genes were also detected in other related abalone species, suggesting that these GCNVs are derived from both ancestral and de novo mutations. Our findings suggest that GCNVs have been accumulated in the early phase of free-spawning abalone speciation.

Allopatric origin, secondary contact and subsequent isolation of sympatric rockfishes (Sebastidae: Sebastes) in the north-western Pacific

Understanding how speciation occurs is central to biology. Gene flow between diverging taxa is correlated with geography and other aspects of speciation; therefore, the examination of gene flow during divergence is a potent approach to understanding the nature of speciation. Here, we inferred the speciation process of the sympatric rockfishes Sebastes steindachneri and Sebastes wakiyai in the north-western Pacific and its marginal seas based on genome-wide single nucleotide polymorphism and mitochondrial DNA data. Model-based demographic inference showed that gene flow between the two species was absent in the initial and late stages of divergence and present only in the middle stage. Population expansion occurred before or during the period of gene flow. The estimated timings of the initial divergence and population expansion fell within the Pleistocene, during which the seas currently inhabited by the two species were repeatedly isolated and reconnected. Contemporary isolation was supported by the absence of hybrids and the shared mitochondrial DNA haplotypes. Our results suggest that the two species initially diverged in allopatry, followed by secondary contact and introgression and by the completion of reproductive isolation. Given that complete isolation following secondary contact has rarely been tested or documented in marine organisms, we highlight the importance of careful consideration of alternative divergence scenarios to be tested, which should take into account the geological and environmental settings.

Speciation genomics and the role of depth in the divergence of rockfishes ( Sebastes ) revealed through Pool‐seq analysis of enriched sequences

Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth‐segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth‐segregated and have different degrees of divergence: S. chlorostictus–S. rosenblatti and S. crocotulus–S. miniatus. We were able to reliably identify “islands of divergence” in the species pair with more recent divergence (S. chlorostictus–S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.

Application of palaeogenetic techniques to historic mollusc shells reveals phylogeographic structure in a New Zealand abalone

Natural history collections worldwide contain a plethora of mollusc shells. Recent studies have detailed the sequencing of DNA extracted from shells up to thousands of years old and from various taphonomic and preservational contexts. However, previous approaches have largely addressed methodological rather than evolutionary research questions. Here we report the generation of DNA sequence data from mollusc shells using such techniques, applied to Haliotis virginea Gmelin, 1791, a New Zealand abalone, in which morphological variation has led to the recognition of several forms and subspecies. We successfully recovered near-complete mitogenomes from 22 specimens including 12 dry-preserved shells up to 60 years old. We used a combination of palaeogenetic techniques that have not previously been applied to shell, including DNA extraction optimized for ultra-short fragments and hybridization-capture of single-stranded DNA libraries. Phylogenetic analyses revealed three major, well-supported clades comprising samples from: 1) the Three Kings Islands; 2) the Auckland, Chatham and Antipodes Islands; and 3) mainland New Zealand and Campbell Island. This phylogeographic structure does not correspond to the currently recognized forms. Critically, our non-reliance on freshly collected or ethanol-preserved samples enabled inclusion of topotypes of all recognized subspecies as well as additional difficult-to-sample populations. Broader application of these comparatively cost-effective and reliable methods to modern, historical, archaeological and palaeontological shell samples has the potential to revolutionize invertebrate genetic research.

Ancient rivers shaped the current genetic diversity of the wood mouse (Apodemus speciosus) on the islands of the Seto Inland Sea, Japan

The current distributions of organisms have been shaped by both current and past geographical barriers. However, it remains unclear how past geographical factors-currently cryptic on the sea floor-affected the current distributions of terrestrial animals. Here, we examined the effects of currently cryptic ancient rivers on current genetic differentiation of the large Japanese wood mouse, Apodemus speciosus, which inhabits islands in the Seto Inland Sea, Japan. Genome-wide polymorphisms were identified by GRAS-Di (Genotyping by Random Amplicon Sequencing, Direct) analysis of 92 A. speciosus individuals. Maximum-likelihood analysis was performed with 94,142 single nucleotide polymorphisms (SNPs) identified by GRAS-Di analyses. Ancient rivers were visualized by Geographic Information System (GIS) channel analysis. Maximum-likelihood analysis showed strong support for the monophyly of each population in the islands in the Seto Inland Sea; it also showed close relationships between Innoshima-Ikuchijima, Ohmishima-Hakatajima-Oshima, Ohmishima-Hakatajima, Ohsakikamijima-Ohsakishimojima, Kamikamagarijima-Shimokamagarijima, and Kurahashijima-Etajima islands. The principal component analyses of the SNPs also supported these relationships. Furthermore, individuals from islands located on the east and west sides of the main stream of the ancient river were clustered on each side with strong support. These phylogenetic relationships were completely congruent with the paleogeographic relationships inferred from ancient rivers. In conclusion, the findings demonstrated that the current distribution of genetically distinct island lineages was shaped by ancient rivers that are currently submerged beneath the Seto Inland Sea, Japan.

Species persistence with hybridization in toad-headed lizards driven by divergent selection and low recombination

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

Susceptibility of Four Abalone Species, Haliotis gigantea, Haliotis discus discus, Haliotis discus hannai and Haliotis diversicolor, to Abalone asfa-like Virus

Abalone amyotrophia is a viral disease that causes mass mortality of juvenile Haliotis discus and H. madaka. Although the cause of this disease has yet to be identified, we had previously postulated a novel virus with partial genome sequence similarity to that of African swine fever virus is the causative agent and proposed abalone asfa-like virus (AbALV) as a provisional name. In this study, three species of juvenile abalone (H. gigantea, H. discus discus, and H. diversicolor) and four species of adult abalone (the above three species plus H. discus hannai) were experimentally infected, and their susceptibility to AbALV was investigated by recording mortality, quantitatively determining viral load by PCR, and conducting immunohistological studies. In the infection test using 7-month-old animals, H. gigantea, which was previously reported to be insusceptible to the disease, showed multiplication of the virus to the same extent as in H. discus discus, resulting in mass mortality. H. discus discus at 7 months old showed abnormal cell masses, notches in the edge of the shell and brown pigmentation inside of the shell, which are histopathological and external features of this disease, while H. gigantea did not show any of these characteristics despite suffering high mortality. Adult abalones had low mortality and viral replication in all species; however, all three species, except H. diversicolor, became carriers of the virus. In immunohistological observations, cells positive for viral antigens were detected predominantly in the gills of juvenile H. discus discus and H. gigantea, and mass mortality was observed in these species. In H. diversicolor, neither juvenile nor adult mortality from infection occurred, and the AbALV genome was not increased by experimental infection through cohabitation or injection. Our results suggest that H. gigantea, H. discus discus and H. discus hannai are susceptible to AbALV, while H. diversicolor is not. These results confirmed that AbALV is the etiological agent of abalone amyotrophia.